Ultrashort vortex from a Gaussian pulse – An achromatic-interferometric approach

The more than a century old Sagnac interferometer is put to first of its kind use to generate an achromatic single-charge vortex equivalent to a Laguerre-Gaussian beam possessing orbital angular momentum (OAM). The interference of counter-propagating polychromatic Gaussian beams of beam waist ωλ with correlated linear phase (ϕ0 ≥ 0.025 λ) and lateral shear (y0 ≥ 0.05 ωλ) in orthogonal directions is shown to create a vortex phase distribution around the null interference. Using a wavelength-tunable continuous-wave laser the entire range of visible wavelengths is shown to satisfy the condition for vortex generation to achieve a highly stable white-light vortex with excellent propagation integrity. The application capablitiy of the proposed scheme is demonstrated by generating ultrashort optical vortex pulses, its nonlinear frequency conversion and transforming them to vector pulses. We believe that our scheme for generating robust achromatic vortex (implemented with only mirrors and a beam-splitter) pulses in the femtosecond regime, with no conceivable spectral-temporal range and peak-power limitations, can have significant advantages for a variety of applications.

proposed for generating achromatic vortex through dispersion compensation [37][38][39][40][41] , by introducing a radial phase modulation for spiral phase plates 42 , via spin-to-orbit conversion through Geometric phase accumulation 21,[43][44][45][46] and using conical glass devices exploiting total internal reflection [47][48][49] . However, these schemes require a dedicated optical component for the task and its fabrication process often turns out to be expensive and demanding to keep the efficiency and quality of the generated achromatic vortex at an acceptable level. Moreover, the dedicated devices that works reasonably well for visible spectrum may not be suitable for other regions of the electromagnetic spectrum.
We propose and demonstrate the generation of an achromatic single-charge vortex merely from the interference of two Gaussian beams in a Sagnac interferometer with correlated linear phase and lateral shear. The SI configuration envisioned satisfies the condition for vortex generation for the entire range of visible wavelengths, enabling realization of white-light vortex. We demonstrate the practical application of the design by generating ultrashort vortex pulses. It is worth noting that high quality ultrashort single-charge vortex pulses can now be generated elegantly using inexpensive and off-the-shelf optical components. The second harmonic generation resulting in doubling of the topological charge for the generated vortex demonstrates its suitability for nonlinear optical applications 50 . Further, we also demonstrate the application capability of the generated ultrashort vortex pulse in femtosecond regime by transforming it to achromatic vector pulses [51][52][53][54][55][56][57][58] by designing a cylindrical vector beam generator using modified Mach-Zehnder interferometer.

Results
We begin by validating the achromatic nature of the proposed Sagnac vortex generator described in Fig. 1, by experimentally generating a single-charge vortex for optical field components spanning the visible range using wavelength-tunable continuous wave laser.
As shown in Fig. 1, in a Sagnac interferometer 1-9, 59 , a controlled amount of linear phase difference with a value φ 0 at the beam waist ω λ corresponding to a spectral component with wavelength λ, is introduced between the two counter-propagating beams along x-by rotating mirror M 2 by an amount α such that α = φ λ ω λ tan 0 0 . By tilting up the mirror M 1 and tilting down the mirror M 3 a required amount of lateral shear y 0 for the beam along y-is also introduced simultaneously.
The surface plots of real and imaginary parts of the superposition of out-of-phase Gaussian beams resulting from the application of the linear phase applied along x-, the lateral shear along y-and their simultaneous correlated presence are shown in Fig. 2(a-c), respectively. The crossing of the surface of imaginary part resulting from the application of linear phase with the plane of real part with zero magnitude results in a line singularity with a phase profile shown at the base of the Fig. 2(a). Similarly, the crossing of surface of real part resulting from the application of lateral shear with the plane of imaginary part with zero magnitude results in another line singularity with a phase profile shown at the base of the Fig. 2(b). It is worth noting that though the phase profiles in both the cases exhibits a jump of an amount π across the line singularity, they are shifted by π/2 with respect to each other. As represented in Fig. 2(c), it is the 'in quadrature' nature of the two line singularities that eventually results in formation of a point singularity having a spiral phase when the linear phase and lateral shear are introduced in orthogonal directions. The equations used in the simulations are detailed in methods section and the values used are discussed below. For the rest of the wavelengths the beam waist is chosen such that, ω λ = pω λ0 with the factor p attributed to the beam size scaling for different spectral components collimated from a point-like source such as output of an optical fiber. The experimentally obtained intensities of single-charge vortices corresponding to the different spectral lines recorded merely by scanning the wavelength using He-Ne laser followed by Argon-ion laser are shown in the second row of Fig. 3. The spiral interference fringes, indicating the presence of phase singularity, due to the interference of the generated vortex beams with a spherical wavefront reference beam are shown in the third row of Fig. 3. The recorded intensities of the generated vortex beam after a free-space propagation of 5 metres are shown in the fourth row of Fig. 3. The propagation stability and integrity of the generated vortex beam behaving like a Laguerre-Gaussian beam of first order are put to test when the generated vortex is used for second-harmonic generation and vector beam generation in the femtosecond pulse regime.
Ultrashort Gaussian to vortex pulses. A 300 femtosecond pulsed fiber laser having a repetition rate of 200KHz with mean wavelength 1030 nm (SATSUMA, Amplitude Systemes fiber laser) is frequency doubled to 515 nm using a BBO crystal. The spatial intensity profile of the femtosecond pulse is given in Fig. 4(a). The single charge vortex generated using the Sagnac achromatic vortex generator is shown in Fig. 4(b). The fork pattern in the fringes at the dark core of the vortex due to interference of the generated vortex pulse with its laterally separated and tilted copy confirms the presence of phase singularity.
Interchanging the sequence of second harmonic generation and vortex generation, we demonstrate the charge doubling through frequency doubling of a 800 nm and 100 fs single-charge vortex pulse (repetition rate of 1 KHz;  Ti-Sapphire laser system from Spectra Physics). Ultrashort vortex to vector pulses. A superposition of the optical field with opposite topological charges +l and −l having orthogonal spin σ + and σ − (circular states of polarization) results in the generation of cylindrical vector beams [51][52][53] . In many interferometric schemes for cylindrical vector beam generation with a Gaussian beam at the input, vortex generators such as spiral phase plates are introduced inside the interferometer where the counter propagating beams acquire opposite topological charge [53][54][55][56][57][58] . To generate a vector beam from a vortex beam having a topological charge +l, a 'Sagnac-like' interferometer is recently proposed 60 . The experimental geometry used generates vector beams with a quasi-monochromatic light source, but could not be applied to ultrashort pulses due to unequal path delay and dispersion incompatibility between the interfering beams.
We demonstrate the generation of ultrashort cylindrical vector pulses from a scalar single-charge vortex pulse by designing a modified Mach-Zehnder interferometer. The scalar, single-charge vortex generated using a Sagnac achromatic vortex generator with a frequency doubled, 1030 nm, 300fs laser pulse is fed to the input of the Mach-Zehnder vector beam generator shown in Fig. 6. It can be inferred that the part of the pulse split using a beam splitter (BS 2 ) that travel through arm I comprising of mirrors M 4 and M 5 has a total of 3 reflections whereas the part of the pulse that travel through arm II comprising of mirrors M 6 and M 7 has a total of 4 reflections. The odd-even number of reflections ensures that the pulses have opposite topological charge at the output of the interferometer.
Additionally, the state of polarization of the pulses are made orthogonal by introducing a half wave plate (HWP 1 ) in arm I. To guarantee that they remain orthogonally polarized during their recombination and for optimum utilization of the pulse power, we used a polarizing beam splitter (PBS) at the output of the interferometer. A quarter wave plate (QWP) transforms the orthogonal linear polarizations to orthogonal spins for the pulses. As shown in Fig. 7(a-c), we observed an azimuthally polarized vector pulse after the QWP that is converted to hybrid and radially polarized pulses using half wave plates HWP 2 with its fast axis oriented along 45° with respect to the x-axis and HWP 3 oriented along 90° with respect to the x-axis, respectively.  Fig. 1; (c) the interference of the generated vortex pulse with its laterally separated tilted copy reveals the forking of fringes at the dark core of the interfering vortices. We confirm the vector nature and type of the pulses by introducing a rotating analyser where a rotating two-lobe intensity pattern shown in Fig. 7(d-f) corresponding to Fig. 7(a-c) reveals the underlying spatial distribution of the polarization of the pulses.

Discussion
Optical dispersion plays a decisive role in the generation of vortices for optical field with spectral diversity. As mentioned earlier, different approaches are designed to compensate or overcome it. In the proposed interferometric approach of introducing lateral shear and linear phase difference in orthogonal directions between superposing out-of-phase Gaussian beams, a hinge-point thus appearing plays a crucial role 61 . The source of achromatic nature of the generated vortex can be easily understood from the well-known optical white-light interferometry. The importance of zeroth-order fringe in white-light interferometry where the phase difference of all the interfering spectral components has a zero crossing to obtain a 'white fringe' on constructive interference or the 'dark fringe' on null interference with the choice of the output port of the interferometer was used to pin point the location of zero optical path difference (OPD) in optical metrology schemes 62 . From the null-interference, a correlated lateral shear and linear phase difference applied orthogonal to each other is shown to be capable of inducing a vortex around the null point for all the spectral components. The choice of Sagnac interferometer with common-path configuration right away satisfies a stable vortex pattern around zero OPD. For the chosen output of the Sagnac interferometer, the clock wise propagating Gaussian beam gets transmitted twice through beam splitter (BS) and its counterpart gets reflected twice at the BS. This asymmetry renders the output port unbalanced and satisfies the condition for out-of-phase or null interference 63 . The evenness in introducing linear phase and lateral shear between the interfering beams ensured that the hinge-point is located at x,y = 0 to generate a high-quality single-charge vortex with excellent propagation stability. The novelty of the work presented is the first of its kind use of the century-old Sagnac interferometer without any modifications or introduction of optical elements such as lenses or prisms inside the interferometer, to generate an achromatic vortex from null-interference of Gaussian beams and is therefore not compromised by the demonstrations in earlier works 25,27 . In the proposed scheme of vortex generation using the typical Sagnac interferometer, any dispersion leading to stretching of the pulse, experienced equally by both the counter-propagating beams, does not affect the spatial profile of the generated vortex, as exemplified by the cw studies. However, the conversion of the  generated achromatic vortex pulse to a vector pulse using the Mach-Zehnder interferometer requires stability of the non-commonpath interferometer and the use of achromatic wave plates. The second harmonic generation resulting in doubling of the vortex charge and the transformation of vortex pulses into vector pulses are included in our reports to substantiate that the generated vortex pulse behaves exactly like a Laguerre-Gaussian beam. We believe that our scheme for generating achromatic vortex pulses in femtosecond regime with only mirrors and beam splitter, the simplest of optical components, has an unprecedented advantage over the existing optical vortex beam generation techniques paving way for subwavelength shaping of high-power ultrashort electromagnetic pulses across the spectrum, not limited to extreme ultraviolet and X-rays 64,65 .

Method
Interferometric generation of achromatic vortex from Gaussian beam. Let us assume a Gaussian spatial profile of a source of light diverging from the ouput of a fiber as where ∼ A is the maximum amplitude, ω  0 radius of the light spot having a Gaussian profile and  x and  y are the spatial coordinates. The spectral component corresponding to the central wavelength λ 0 in the light beam collimated from such a point-like source using a lens of focal length f can be written as We can treat the profile of Eλ 0 (x,y) as well a Gaussian such that where A is the maximum amplitude, ω λ0 is the beam waist and x and y are the spatial coordinates. The field distribution corresponding to any other spectral component having wavelength λ = pλ 0 can be written as The parameters y 0 and φ 0 controlling the lateral shear and linear phase for the two beams are described in Fig. 2. As discussed, the asymmetry in number of transmissions and reflections at the beam splitter for the counter-propagating beams offers the condition for an out-of-phase superposition of the two Gaussian beams at the chosen output of the interferometer. The optical field distribution at the output of the Sagnac interferometer can be represented as Restricting the variation of the real and imaginary parts of the argument of 'sinh' term to the linear region of the hyperbolic sine function and sine function respectively, we can simplify Eq. (7) as To understand the influence of difference in the ratio of lateral shear y 0 to the beam waist ω λ of different spectral components on the symmetry of the generated single-charge vortex, a parameter plot of lateral shear y 0 (in units of ω λ ) against linearly increasing phase φ 0 at ω λ (in units of 2π) is shown in Fig. 8. The white line represents the locus of points satisfying the condition, ω π φ = λ y sinh(2 / ) 2 0 0 0 on which the generated vortices for the central wavelength of λ 0 = 545 nm are canonical with a symmetric dark core. For the same parameter settings of y 0 and φ 0 , the vortices generated for other spectral components can be calculated by finding the locus in parameter space with condition ω π φ = λ p y sinh(2 ( / )) 2 0 0 0 . The red, green and blue lines describe the locus of the single-charge vortices that can be generated for λ = 632.8 nm of He-Ne laser and λ = 514 nm and λ = 457 nm from Argon-ion laser with the same parameters that generated canonical vortex for the central wavelength of λ 0 = 545 nm The slight asymmetry for the core of vortices generated for λ = 632.8 nm and λ = 457 nm compared to the symmetric core for central wavelength of λ 0 = 545 nm can be inferred. In the realm of first-order optical modes comprising of Hermite-Gaussian and Laguerre-Gaussian modes that are represented on a modal/Padgett sphere 29,30 equivalent to the Poincaré sphere representation of the state of polarization, the stretching of the symmetric vortex core, for wavelengths neighbouring the central wavelength can be attributed to fractional orbital angular momentum (charge) in a way the fractional spin angular momentum is attributed to the elliptical state of polarization.
By scaling the lateral shear y 0 introduced for different spectral components of the field such that y 0 (λ)=y 0 p, by introducing dispersion inside the interferometer, it could be possible to maintain the single-charge vortex canonical for all the spectral components, if needed. In our experimental design using Sagnac interferometer, we preferred not to introduce dispersion for the lateral shear. For many practical applications such as the generation of high-power ultrashort vortex and vector pulses, this aspect can be safely ignored due to negligibly small asymmetry for the dark core in the vortices corresponding to spectral components spanning the visible spectrum as shown in Fig. 3.

Interferometric generation of cylindrical vector beam from a scalar vortex beam.
We briefly describe the method of generating a cylindrical vector beam from a scalar single-charge vortex beam. As shown in Fig. 6, at the output of the modified Mach-Zehnder interferometer, the quarter wave plate QWP introduced is oriented such that the states of polarization of the vortex beams with topological charge l = +1 and −1 are transformed to left (σ − ) and right (σ + ) circular states, respectively. Their on-axis, out-of-phase superposition generating an azimuthally polarized vector beam can be written in Jones matrix formalism as x azim y azim The out-of-phase superposition can be attributed to the choice of the output port of the interferometer. At the introduction of half-wave plate HWP 2 following QWP, with its fast-axis oriented at an angle θ 2 with respect to x axis, the transformation of the azimuthally-polarized optical field represented in Eq. (9) can be obtained by a rotation operation performed through R(θ) defined as  At the introduction of half-wave plate HWP 3 following HWP 2 , with its fast-axis oriented at an angle θ 3 with respect to x axis, the transformed optical field can be written as For θ 2 = π/4 and θ 3 = π/2, the optical field distribution represented in Eq. (13) describes a radially polarized optical field given by x rad y rad x azim y azim 3 2 The sinusoidal variation of x and y components of the cylindrical vector optical field as a function of azimuthal coordinate is observed through rotating analyser as described in Fig. 7 in order to confirm their generation experimentally.