Polarization dependent nanostructuring of silicon with femtosecond vortex pulse

We fabricated conical nanostructures on silicon with a tip dimension of ∼ 70 nm using a single twisted femtosecond light pulse carrying orbital angular momentum (` =±1). The height of the nano-cone, encircled by a smooth rim, increased from ∼ 350 nm to ∼ 1 μm with the pulse energy and number of pulses, whereas the apex angle remained constant. The nano-cone height was independent of the helicity of the twisted light; however, it is reduced for linear polarization compared to circular at higher pulse energies. Fluid dynamics simulations show nano-cones formation when compressive forces arising from the radial inward motion of the molten material push it perpendicular to the surface and undergo re-solidification. Simultaneously, the radial outward motion of the molten material re-solidifies after reaching the cold boundary to form a rim. Overlapping of two irradiated spots conforms to the fluid dynamics model. © 2017 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). [http://dx.doi.org/10.1063/1.4999219]

APL Photonics 2, 086104 (2017)   crater while an OV pulse formed a nano-cone encircled by a rim, irrespective of laser polarization.We simulate the action of the laser pulse on fluid dynamics of the molten silicon and show that the radial inward (outward) motion of the molten silicon forms the nano-cone (rim).We demonstrate the nano-cone height to decrease at higher pulse energies for linear compared to that of circular polarization.This behaviour was attributed to an asymmetry in the radial inward fluid motion that arises from transient light-plasma interaction.
Linearly polarized light from a Ti:Sapphire laser (operating at 1 kHz repetition rate and producing 800 nm, 45 fs pulses) was transformed into the left-or right-circularly polarized beam (s = ±1) by a quarter-wave plate (QWP).A q-plate with a topological charge of q = 1/2 produced OV beams with = ±1 and inverted polarization (s = ∓1). 26A linearly polarized OV beam (s = 0, = ±1) was produced by a second QWP and a polarizer.Electrically detuning the q-plate generated Gaussian beams.A 0.25 NA (16×) aspheric objective focused the femtosecond pulse onto a p-type silicon surface [(100) orientation] mounted on a three-axis translation stage.
For a Gaussian beam, the focal spot size was measured to be 1.8 ± 0.1 µm using the knife-edge method.For an OV beam, the spot size, defined as the diameter of the maximum intensity, was measured to be 2.8 ± 0.2 µm compared to a theoretical value of 2.3 µm. 29The ablated regions were characterized by a scanning electron microscope (SEM) with the electron beam perpendicular to the sample (zero tilt) and at an angle of 88 • and by an atomic force microscope (AFM) in non-contact mode.
Figure 1  OV pulses ( = 1, s = +1), respectively.The Gaussian pulse produced a rim with a uniform height of ∼ 50 nm around a ∼ 70 nm deep ablation crater.Besides the rim, the OV pulse also produced a nanocone in the middle of the ablation crater protruding ∼ 450 nm above the surface.Similar structures were observed for intrinsic and n-type silicon.1][32] A small asymmetry in the rim shape was due to a slight misalignment of the beam through the q-plate.
Evolution of the ablation crater features with pulse energy for a single femtosecond vortex pulse is shown in Fig. 2. Four key observations can be made.
First, the nano-cone height was independent of the handedness of the twisted light ( = ±1) with spin (s = ∓1), corresponding to a total angular momentum, j = + s = 0. Formation of a nano-cone with a single femtosecond OV pulse differs from previous reports in three aspects.(i) Only pico-and nano-second pulses produced micro-needles, 11 while femtosecond pulses produced annular ablation craters with no features in the centre.(ii) Ablation with a single twisted nanosecond OV pulse with j = 0 produced no central protuberance. 18It required overlapping of several twisted light pulses to produce a nanostructure.(iii) The tightly focused circularly polarized femtosecond vortex pulse ablated the material differently depending on the polarization handedness of light. 33econd, the nano-cone height increased with pulse energy at two different rates: ∼ 0.6 nm/nJ and 0.2 nm/nJ below and above 600 nJ, respectively, for left-and right-circular polarizations ( j = 0).
Third, the nano-cone height increased at the same rate for linear ( j = +1) and circular ( j = 0) polarizations up to 600 nJ.Beyond this energy, the nano-cone height decreased for linear polarization at a rate of ∼1.3 nm/nJ (highlighted by an ellipse).This is discussed below in terms of an asymmetry in compressive forces that causes nano-cone formation.
Fourth, rim height and crater depth increased with pulse energy at the same rate for all three polarizations.
Figure 3 shows variation of the nano-cone width, measured at the surface, and the apex angle with pulse energy.The nano-cone width varied from 1.4 µm to 1.9 µm with two different rates  -0.6 nm/nJ (0.18 nm/nJ) below (above) 700 nJ, similar to the height.However, the apex angle (measured with an AFM and SEM) remained nearly constant at 105 • , decreasing marginally at higher energies.The radius of curvature of the nano-cone was about 70 nm.The insets show SEM images of the nano-cone at a pulse energy of 370 nJ.While the top view (left panel) suggests a spiral conical structure, it could not be confirmed by the side view (right panel) even when the handedness of the twisted light was changed.In contrast, spiral needles of several microns in length were produced with pico-and nano-second pulses. 11,17,18This could be due to significantly lower height of the nano-cone and shorter duration of the OV pulse, thereby obscuring the effects of the twisted light.
The height of the nano-cone induced by the OV pulse of fixed energy increased with the number of laser pulses as shown in Fig. 4.However, the aspect ratio, defined as the ratio of the height of the nano-cone to its width, varied marginally from ∼ 0.38 for a single pulse to ∼ 0.45 for 5 pulses.This was in contrast to the narrowing of the micro-needle observed with nanosecond light pulses. 18SEM images show that with increased number of laser pulses, more molten material was pushed to the outer rim forming nano droplets.Also, the inside of the ablation crater developed ripples.
Polarization dependent rim formation on silicon irradiated by a single femtosecond Gaussian pulse has been recently observed [34][35][36] and discussed in terms of fluid dynamics 34 which can be extended to the vortex pulse.The annular intensity distribution of the OV beam is mapped onto the spatial profile of the free carriers generated by linear and two-photon absorption [Fig.5(a)].Picoseconds later, when energy of the free carriers is transferred to the lattice, the temperature distribution retains the shape of OV beam intensity distribution [Fig.5(b)].When melting occurs on nanoseconds, only a thin layer of molten silicon is formed, determined by the skin depth d s .Thermocapillary [37][38][39] and/or hydrodynamic forces 40 displace the molten silicon radially outward to the periphery and also radially inward to the centre of the crater [Fig.5(c)].The outward motion leads to rim formation similar to a Gaussian beam.The inward motion leads to a compressive force that pushes the molten material away from the surface.Rapid expansion causes resolidification into a nano-cone [Fig.5(d)].The temporal evolution of the silicon micro-needle formation was recently captured with an ultra-high speed camera. 11o gain further insight, we simulated the action of the short laser pulse on the fluid dynamics of the molten silicon.We solved the continuity and Navier-Stokes equations using the methodology of Schwarz-Selinger et al. 37 under the assumption of incompressibility, large viscosity, and steadystate flow.To simulate the effect of laser polarization with an asymmetric distribution, we have generalized the calculation to the full two-dimensional transverse coordinate (r and φ or x and y) and extended it to include Laguerre-Gaussian intensity profile.The boundary conditions are such that the fluid velocity must vanish at the bottom of the molten material and that the fluid velocity at the surface is given in relation to the surface tension, i.e., v = 0 at z = h; Here, η is the dynamic viscosity of the fluid, α is the surface tension, and r, φ, and z are the cylindrical coordinates.
Radial velocity and azimuthal velocity ( r and v φ ) were obtained by solving the Navier-Stokes equation.Substituting them in the continuity equation and integrating with respect to z gives For the case in which the melt depth is independent of the transverse coordinates, r and φ, the above expression takes the simple form of, is the transverse Laplacian operator.For simple materials, the surface tension may be given in terms of the temperature at the surface, where ∂ T α is often constant (independent of the temperature).Since the temperature distribution is proportional to the intensity profile of the incoming short pulse, i.e., T (r, φ) ∝ I(r, φ), the mass transport in the fluid can be approximated by v z (r, φ) ∝ ∇ 2 ⊥ I. To calculate the mass transport for various intensity distributions of the laser pulse, we used Laguerre-Gaussian beams given by the following intensity profiles, Here, A is a constant that depends on the power of the laser pulse, L | | p (.) are the associated Laguerre polynomials, and p are the azimuthal and radial mode indices, respectively, ρ = r/w 0 , and 0 is the beam waist at the focus.
Figure 5(e) shows the mass transport for a Gaussian beam with the temperature distribution shown at the bottom.The longitudinal velocity, v z (r, φ), is maximum at the centre of the interaction region pointing towards the thin layer of molten silicon.The resultant displacement causes the molten material to accumulate at the edges where v z (r, φ) is pointing away from the silicon surface.This fluid motion leads to formation of a rim around the ablation crater.
When an asymmetry was introduced in the temperature distribution along a given direction, it is reflected in the mass transport giving rise to an asymmetric rim height and width as shown in Fig. 5(f).Such an asymmetry was observed for linearly polarized light, whereas it did not exist for circular polarization. 34The polarization dependence was explained in terms of transient plasmonics. 41,42A dense free carrier plasma is created by the leading edge of an intense femtosecond pulse.The trailing edge of the pulse then interacts with the plasma giving rise to local field enhancement.The resultant asymmetry in the carrier distribution eventually translates into an asymmetric distribution of lattice temperature and hence thermocapillary/hydrodynamics forces.The field enhancement occurs along or perpendicular to the laser polarization depending on the plasma density. 42or an OV beam, the annular intensity profile translates into a doughnut shaped temperature distribution as shown in Fig. 5(g).Where temperature is maximum, mass transport points towards the molten silicon.However, v z (r, φ) is maximum at the centre, pointing away from the silicon surface.The resultant displacement causes molten silicon to move radially both outward and inwards forming the rim and nano-cone, respectively.
We now address the decrease in nano-cone height at higher pulse energies for the linearly polarized vortex pulse, shown in Fig. 2. Local field enhancement resulting from light-plasma interaction leads to an asymmetry in the radial inward motion of the molten material.The resultant asymmetric compressive force causes the nano-cone height to be less.The local field enhancement vis-a-vis asymmetric compressive force becomes prominent at higher pulse energies.As a result, the height of the nano-cone with linear polarization deviates from that of a circularly polarized vortex pulse.
The role of fluid dynamics in the nano-cone formation is evident from SEM images of Fig. 6.Silicon was irradiated with a single vortex pulse at two different spots.At 16 µm separation [Fig.6(a)], the ablation crater morphology was unaffected by the presence of the other.At 15 µm spacing, molten silicon from the second spot collided with the resolidified rim of the first spot deforming the rim at the point of contact [Fig.6(b)].When the two spots overlapped (10 µm separation), the nano-cones remained unaffected, but the rim of the first spot vanished in the overlap region [Fig.6(c)].When the spacing was reduced further, the molten material was propelled past the first nano-cone [Fig.6(d)] and eventually disrupted it [Fig.6(e)].Intactness of the nano-cones even when the irradiated spots overlapped enables one to fabricate arrays of silicon nano-cones similar to those shown in Fig. 6(f).
To conclude, we demonstrated fabrication of a nano-cone on silicon with a single femtosecond OV pulse.Its formation was described in terms of fluid dynamics, and its polarization dependence was attributed to transient light-plasma interaction.Our results highlight key aspects of light-matter interaction.First, the fluid dynamics model suggests that any tightly focussed annular beam will produce a nano-cone.Tighter focusing by higher numerical aperture objectives can lead to higher compressive forces and can further increase the height of the nano-cone.Second, transient plasmonics influences the laser processed structures.Intense linearly polarized femtosecond light pulses distort the free-carrier density distribution due to local field effects that is imprinted on the subsequent steps of the ablation/modification process. 34,41,43Its effects can be minimized by using circularly polarized light.
5][46][47] For example, an array of nano-antennas was used to show enhancement of plasmon assisted high harmonic generation from crystalline silicon substrate. 48Also, nanostructured ZnO and silicon were used to demonstrate customized high-harmonic wave fields and to confine them to diffraction-limited spot sizes. 49tructured light enables to not only alter the surface topography of a target material and possibly change its phase but also investigate the target by generating high harmonics.This all-optical technique can provide an opportunity to control and probe light-induced dynamical changes to material properties.

FIG. 3 .
FIG. 3. Nano-cone width (left axis) and the apex angle (right axis) as a function of pulse energy for a single femtosecond OV pulse ( = +1, s = 1).The insets show SEM images of a nano-cone at a pulse energy of 370 nJ; top view (left panel) and side view at a tilt angle of 88 • (right panel).

FIG. 5 .
FIG. 5. (a)-(d) Schematic showing rim and nano-cone formation with an OV beam at different stages of the interaction.(a) Intensity distribution of the incident pulse.k and E are propagation and polarization directions, respectively.(b) Spatial profile of the temperature distribution d s is the melt layer thickness.(c) Displacement of molten silicon indicated by arrows and melt layer velocity profile due to the exerted forces.(d) Nano-cone (rim) formation by radial inward (outward) motion of molten silicon.(e)-(g) Mass transport v z (r, φ) for several temperature distributions (shown at the bottom) given by (e) symmetric Gaussian profile, (f) an asymmetric Gaussian profile, and (g) a symmetric Laguerre-Gaussian profile with a mode index of = 1 and p = 0.The asymmetry in (f) was obtained by multiplying the Gaussian distribution by (1 + y 2 ), where we consider the first order symmetric distortion.