Laguerre-Gaussian Beams with an Increased Dark Area and Autofocusing

: We introduce and investigate a novel Laguerre-Gaussian (LG) beam, different from the conventional modal LG beams, which conserve the transverse intensity structure (up to scale) on propagation. The proposed beam does not conserve its structure on free space propagation but possesses some interesting properties. This beam is Fourier-invariant, and it has an increased dark area both in the initial (waist) plane and in the far field. Thus, without changing the topological charge of the beam, varying the radial (lower) index of the associated Laguerre polynomial allows increasing or decreasing the effective diameter of the central dark spot in the intensity pattern. In addition, the beam is autofocusing, i.e., the intensity distribution at the Rayleigh distance from the waist has a shape of the light ring (at any value of the radial index) with the minimal diameter and with the maximal on-ring intensity. Such a beam can be adopted for microparticle manipulation. Increasing the dark area in the focus of a high-aperture spherical lens allows the simultaneous trapping of several absorbing microparticles into this dark area.


Introduction
For many years, vortex Laguerre-Gaussian (LG) beams have attracted a a great deal of interest due to their wide practical applications in optical trapping [1], wireless telecommunications [2], atmosphere probing [3], quantum informatics [4], atom cooling [5], optical microscopy of single molecules [6] and quantum dots [7].Besides, the LG beams have been very thoroughly studied theoretically.For instance, the possibility of the orbital angular momentum (OAM) of light was at first shown right for the LG beams [8].There are known many modifications of the LG beams, such as elegant LG beams [9], Hermite-Laguerre-Gaussian beams [10], asymmetric LG beams [11], and products of LG beams [12].The LG beams can be generated by lasers [13][14][15][16], mode converters [17], and using a spatial light modulator [18].However, despite the long history of the LG beams, their potential is still not exhausted.
In this work, we investigate a novel LG beam with its topological charge being equal to the difference between the azimuthal (upper) and radial (lower) indices of the associated Laguerre polynomial.In contrast, the radial polar coordinate is raised to a power equal to the sum of these indices.Thus, in contrast to the conventional LG beams, the topological charge of the studied beam is not equal to the power of the radial coordinate.On the one hand, it violates the modal propagation (shape invariance) of such a beam in free space.On the other hand, the beam discussed demonstrates interesting properties: autofocusing and the ability to control the diameter of the central dark intensity spot without changing the topological charge of the beam.

Fourier-Invariant Laguerre-Gaussian Beams with an Increased Dark Area
Here, we consider a coherent paraxial monochromatic light field with the following complex amplitude in the initial plane (z = 0): ( , ) exp ( ) where (r, φ) are the polar coordinates, w is the waist radius of the Gaussian beam, ( ) is the associated Laguerre polynomial.Using standard LG modes the beam (1) can be written as a product of two of them: LG LG 2 2 The structure of its light rings does not conserve under the beam propagation in free space.However, as we will see later, the beam from Equation (1) has some interesting properties useful in practice: (i) the beam from Equation ( 1) is Fourier-invariant, (ii) it is also autofocusing at the Rayleigh distance from the waist, and (iii) it has an increased dark area near the optical axis, i.e., the diameter of the central dark intensity spot can be tuned (increased) without changing the topological charge.
Equation (1) indicates that the topological charge of the optical vortex equals m − n, and the orbital angular momentum of the beam (1), normalized to the beam power, is also equal to m − n.Initially, the beam has n + 1 light rings.However, since the beam (1) does not conserve its structure on propagation, the number of rings can change.Thus, n + 1 is the maximal number of rings of the beam (1).Differentiating the intensity of the field (1), r , with respect to r and equating the derivative to zero, we find the coordinates of intensity maxima (minima are evidently zeroes of the beam (1)): ( ) where y = (r/w) 2 .Equation ( 4) is an algebraic equation of the (n + 1)th order and can be analytically solved only for small values of n only.For instance, since at n = 1, the Laguerre polynomial is equal to 1 ( ) m L y y m = − + + , Equation (2) allows obtaining the radius of the first light ring: ~2 r w m .In the general case, the radius of the first light ring depends on the radial index n, as follows from Equation (4).Below, we consider a particular case of Equation ( 1), when the topological charge of the beam is nonnegative.Replacing m n = + , the beam (1) can be rewritten as follows, ( ) exp where 0 ≥  .Equation (5) shows that, without changing the topological charge ( const =  ), the effective diameter of the dark area near the optical axis can be changed by varying only the index n.When n increases, the dark low-intensity domain near the optical axis becomes larger, whereas decreasing n reduces this domain.
The beam (1) does not keep the intensity shape during propagation in the Fresnel diffraction zone.However, it has an invariant Fourier image.Namely, where is a 2D vector, and f is the focal length of the spherical lens that performs the Fourier transform.It is seen from Equation ( 6) that the complex amplitude of the light field in the lens focus is equal to the amplitude of the light field (1) up to scale and a constant multiplier.

Fresnel Transform of the Laguerre-Gaussian Beam with an Increased Dark Area
To evaluate the Fresnel transform of the initial light field (3), we apply the expansion: which may be rewritten in terms of LG modes: LG ( ) The coefficients ( , , , ) Here, we use the notation proposed in [19].Namely, if ( ) A t is any power series • ex- pressions in Equation ( 9) are integer numbers and may be written in terms of Jacobi polynomials: Substituting p n = and 0 q = into Equation ( 9) we obtain the expansion of the beam (3) on the base of LG modes: LG ( ) with for odd values of k, the expansion (12) has only even terms: According to Equation ( 13), the light field ( 1) is a superposition of a finite number of the conventional LG beams.Therefore, the complex amplitude of the light field (3) can be obtained at an arbitrary propagation distance z from the initial plane.Since LG ( , ) exp (2 1) ( ) LG ( ) ( ) ( ) LG ( , ) where ( ) is the Gaussian beam radius at the distance z, If we recast Equation ( 15) and write all s-independent multipliers out of the sum, then the remained sum is given by ( ) At z = 0, the Gouy phase in Equation ( 16) is zero and the complex amplitude (15) coincides with the initial amplitude (5).In the far field, when z >> z0, the Gouy phase equals π/2 and the exponents in Equation ( 16) become equal to ( 1) 1 s − = , since s = 2p, p = 0, 1, 2, .... Thus, the sums ( 15) are coinciding at z = 0 and at z >> z0.This again proves that the light field from Equation (1) in the far field coincides with the field in the waist.At z = z0, the Gouy phase ψ(z) equals π/4 and the exponents in Equation ( 16) are equal to ( ) , or, since only the terms with even s are nonzero in Equation ( 16), ( 1) , 0,1, 2,..., . Thus, all the terms in the sum are real-valued and of positive sign (this follows from Equation ( 13)).This explains why there is a narrow large-intensity ring near the optical axis in the autofocus (z = z0).
As seen from Equation ( 5), the radial (lower) index of the Laguerre polynomial is equal to n, and, therefore, the Laguerre polynomial has also n roots and the LG beam from Equation (5) has n + 1 light rings.However, the maximal radial index of the Laguerre polynomial in Equation ( 16) is equal to 2n.This index transformation, n to 2n, can be explained easily.The polynomial part of the beam ( 5) is a product of two factors, ( exp( )) r iϕ  and 2 2 ( ) . The first one is a vortex, while the second is a polynomial in r 2 of the order 2n with a degenerate zero at the origin and n simple zeros.During the beam propagation, the vortex keeps unchanged but the degenerate zero r 2n splits into n simple factors 2 2 k r r − with different (possibly complex-valued) rk.Since the polynomial order remains the same, 2n, the beam expansion in LG basis demands for appearance the mode with radial index 2n.
Thus, in this Section, applying the Fresnel transform, we obtained the field distribution of the LG beam (1) in an arbitrary transverse plane.To do this, we decomposed the light field into a superposition of the conventional LG beams.Strictly speaking, the conventional LG beams constitute a full basis, and an arbitrary light field can, in principle be expanded into a series of the conventional LG beams.However, our expansion contains only a finite number of terms, the expansion coefficients are explicitly derived [Equation (13)].All the constituent LG beams have the same azimuthal index (i.e., a vortex factor of the same order) so that the beam (1) contains the central optical vortex, which on propagation is not split and does not move from the center.Such a decomposition remains the topological charge of the beam (1) predictable but makes its radial distribution more flexible.

Numerical Simulation
In this Section, using the Fresnel transform of the LG beam from Equation (1), we compute the transverse intensity and phase distributions at different distances from the initial (waist) plane.We use the following parameters: wavelength λ = 532 nm, waist radius of the Gaussian beam w = 0.5 mm, azimuthal (upper) and radial (lower) indices of the associated Laguerre polynomial are m = 4 and n = 3, respectively.Figure 1 illustrates the intensity (columns 1 and 3) and phase (columns 2 and 4) distributions of the beam from Equation (1) at the following distances from the waist z/z0: 0 (Figure 1a,b), 1/4 (Figure 1c,d), 1/2 (Figure 1e,f), 3/4 (Figure 1g,h), 1 (Figure 1i,j), 1.5 (Figure 1k,l), 2 (Figure 1m,n), 3 (Figure 1o,p), 4 (Figure 1q,r), 10 (Figure 1s,t).According to Figure 1, at the Rayleigh distance z0 from the waist, a single light ring is generated with its radius nearly equal to 0.4 mm.This is equivalent to focusing a Gaussian beam with a waist radius of 0.5 mm by a spherical lens with a very weak numerical aperture, nearly NA = 1.41w/z0 = 0.0005.Increasing the radial index n allows increasing NA several times.We note that a low numerical focusing aperture leads to a large depth of field since the depth of focus (DOF) is inverse proportional to the squared numerical aperture.We also note that in the initial plane, the first light ring has a radius almost twice as large as the waist radius and equal to nearly 0.9 mm.The conventional LG beam with a waist radius of 0.5 mm would have such a radius at the topological charge of 3, whereas the LG beam from Equation (1) has the topological charge m − n = 4 − 3 = 1.Therefore, modeling confirms that the beams (1) have an enlarged dark area in the initial plane (Figure 1a), and at the Rayleigh distance, a focal spot is formed in the form of a narrow light ring (Figure 1i) almost without side lobes.Although to be precise, focusing (maximum intensity on the ring) takes place at a distance of 0.75 z0 (Figure 1g) and 2.38 times greater than the maximum intensity on the ring in the initial plane (Figure 1a).At the Rayleigh distance, the intensity is slightly less (Figure 1i) and 2.37 times greater than in the initial plane.It can also be seen from Figure 1 that the far-field intensity (Figure 1s) coincides with the initial intensity (Figure 1a).This confirms the Fourier invariance property (4) of beams (1).As to the phase evolution, it is seen in Figure 1 that in all transverse planes, the phase contains an (m − n)th-order singularity in the center since all the terms in Equation ( 15) are proportional to exp(ilφ).In addition, the vortex phase rotates on propagation, and the wavefront acquires an additional parabolic shape due to converging and diverging.
Figure 1 was obtained by the numerical Fresnel transform, implemented as a convolution using the fast Fourier transform (implemented in the SciPy library for the Python language).To confirm the expansion (13) into a superposition of the conventional LG beams, Figure 2 illustrates the same patterns as in Figure 1e,f,o,p, but computed analytically via the Laguerre polynomials and the expansion coefficients (13).It is seen that   1e,f (a,b) and as in Figure 1o,p (c,d), but obtained analytically by Equation ( 7) using the expansion coefficients (8).Now when Equation ( 13) is tested, we can use Equation (15) for computing the longitudinal intensity distributions of the beam from Equation (1) and investigate how the beam transforms when the radial and azimuthal indices vary.Figure 3 illustrates these distributions for the following parameters: wavelength λ = 532 nm, waist radius w = 0.5 mm, radial and azimuthal indices (n, m) are, respectively, (0, 1) (Figure 3a), (0, 4) (Figure 3b), (3, 4) (Figure 3c), (3, 7) (Figure 3d), (6, 7) (Figure 3e), (6, 10) (Figure 3f).For these values, the topological charge of the beams in the left column of Figure 1 is 1, whereas, for the beams from the right column, it is equal to 4. Since n = 0 in Figure 3a,b, these beams reduce to the conventional single-ringed LG beams of the topological charge m.These beams do not possess the autofocusing property.At n = 3 and n = 6, Figure 3 confirms that the beams are autofocusing, although the pronounced focus, when the intensity decays more than two times from its maximum, is only at (n, m) = (6, 7).The depth-of-field, in this case, is almost equal to the Rayleigh length: DOF ≈ 0.9z0.In addition, Figure 3 confirms that the radial index affects the size of the dark area inside the main light ring, but in the initial plane and the far field.On the contrary, in the plane of the autofocus, the beam is narrower than the conventional LG beam in its waist (the beam in Figure 3e is narrower than the beam in Figure 3c, which is narrower than the conventional LG beam from Figure 3a).To demonstrate to what extent we can increase or decrease the dark area in focus, Figure 4 depicts the beam (1) with a unit topological charge but with small and large index values: (n, m) = (1, 2) (Figure 4a-c) and (n, m) = (20, 21) (Figure 4d-f).At (n, m) = (1, 2), the focal distance decreases to 0.7z0, and the ring diameter is nearly equal to the ring diameter of the conventional LG beam with the charge of 1 in its waist, i.e., 0.7 mm (2 1/2 w).At (n, m) = (20, 21), the beam is much narrower in the focal plane: its main ring has a diameter of nearly 0.2 mm.The depth-of-field, in this case, is nearly equal to the half of Rayleigh length: DOF ≈ 0.5z0.We believe that further increasing of the values n and m will lead to further decreasing (theoretically unlimited) of the ring radius in the autofocus plane.Still, in this case, the light beam becomes nonparaxial, and we cannot use Equation ( 15) and the Fresnel transform to estimate the ring diameter in the focal plane.

Conclusions
Here, we investigated a new type of Laguerre-Gaussian beam with autofocusing.At a fixed value of the topological charge l of these beams, changing the radial index n of the Laguerre polynomial allows controlling the focusing degree of these beams at the Rayleigh distance from the waist.Although the numerical aperture of such focusing is not high (inversely proportional to the Rayleigh distance z0), the depth of field, in this case, is significant (proportional to z0).Another interesting property of these beams, in contrast to the conventional Laguerre-Gaussian beams, is that at the zero topological charge (l = 0), the beams (5) generate in the focus of a spherical lens a light ring.The diameter of this ring can be adjusted by changing the radial index n and the ratio of the lens focal length and the Rayleigh distance (f/z0).This property of these beams can be used for simultaneous optical trapping of several metallic microparticles into the increased dark area of the light ring in the focus of the spherical lens.We note that the Bessel or Bessel-Gaussian [20] beams are also able to generate a light ring in focus at a zero topological charge.However, the Bessel and Bessel-Gaussian beams are not Fourier-invariant and generate in the focus of one light ring.The investigated here LG beams are Fourier-invariant and can generate in the focus several light rings with nearly the same intensity.
Figure 2a-d are quite identical to Figure 1e,f,o,p, respectively.

Figure 2 .
Figure 2. The same intensity (a,c) and phase (b,d) distributions as in Figure1e,f (a,b) and as in Figure1o,p (c,d), but obtained analytically by Equation (7) using the expansion coefficients(8).

Figure 3 .
Figure 3. Longitudinal intensity distributions of the beam from Equation (1) with the topological charge of 1 (a,c,e) and 4 (b,d,f) when the radial (lower) index n of the associated Laguerre polynomial is as 0 (a,b), 3 (c,d), 6 (e,f).Below each distribution, a longitudinal intensity cross-section is shown.Horizontal and vertical scale marks denote, respectively, z0/2 and 1 mm.Transverse dashed line denotes the autofocus plane z = z0.

Figure 4 .
Figure 4. Transverse initial (a,d) and focal (b,e) intensity distributions, as well as the longitudinal intensity distributions (c,e) of the beam from Equation (1) with the (n, m) = (1, 2) (a-c) and with the (n, m) = (20, 21) (d,f).Curves (b,c,e,f) show the intensity cross-sections.Horizontal scale marks in the transverse patterns denote 1 mm.Horizontal and vertical scale marks in longitudinal patterns denote z0/2 and 1 mm, respectively.