Two simple and high-efficiency techniques for measuring the orbital angular momentum (OAM) of paraxial laser beams are proposed and studied numerically and experimentally. One technique relies on measuring the intensity in the Fresnel zone, followed by calculating the intensity that is numerically averaged over angle at discrete radii and deriving squared modules of the light field expansion coefficients via solving a linear set of equations. With the other technique, two intensity distributions are measured in the Fourier plane of a pair of cylindrical lenses positioned perpendicularly, before calculating the first-order moments of the measured intensities. The experimental error grows almost linearly from ~1% for small fractional OAM (up to 4) to ~10% for large fractional OAM (up to 34).
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
In recent years, optical vortices and laser beams with orbital angular momentum (OAM) have been actively studied, having found many practical uses in wireless optical communications, quantum information, probing atmosphere, and micromanipulation [1,2]. Modern reviews [3,4] and monographs  dealing with these topics have recently been published. Thus, developing simple and high-efficiency techniques for measuring the OAM of such laser beams has been a topical issue. For paraxial laser beams, OAM can be measured using multi-order diffractive optical elements [6,7], diffractive transformation optics , metasurfaces , cylindrical or astigmatic lenses [10–15], interferograms [16–18], triangular apertures  and annular diffraction gratings . While not all of the above-listed techniques are suited for measuring both integer and fractional OAM, an arbitrary OAM can be measured via measuring the intensity of light with a cylindrical lens and computing intensity moments [13,14,21].
It has been proposed [21,22] that the OAM of a paraxial laser beam should be measured using intensity moments. Having been inspired by the idea, here we propose a modernized and simplified version of the technique proposed in Ref , which enables the OAM to be identified from a single measurement of the intensity distribution in an arbitrary transverse plane of the beam of interest. The technique allows both integer and fractional OAM to be measured. The only limitation of the techniques is that it is suited only for laser beams that can be fairly accurately approximated by a finite superposition of Laguerre-Gauss (LG) modes (n, 0). Another technique is based on the approach described in Refs [13,14], which describe the use of a single cylindrical lens for measuring the OAM of radially symmetric beams. Meanwhile, for arbitrary laser beams two cylindrical lenses need to be used. It is worth noting that the technique described in Refs [13,14]. is also an approximation (see Appendix). We found the conditions when this technique is almost accurate (error is less than 1%).
In this work, we suggest for the first time measuring the fractional OAM by registering only one intensity distribution and without using any optical elements. An additional advantage of the proposed method is that besides the full OAM it determines the OAM-spectrum of a beam. In addition, we generalize a proposed in [13,14] method of fractional OAM measurement for asymmetric paraxial beams.
2. Determining OAM from single measurement of intensity distribution
Similar to Ref , we assume the complex amplitude of a laser beam in a plane perpendicular to the optical axis z to be given by23]:Eq. (1) is not suitable for approximating an arbitrary function, many vortex laser beams found in practice, e.g. in optical communication lines , are described by Eq. (1).Eq. (1), with the functions in Eq. (3) having two scale parameters, (α, w0), their use for approximating the light field is preferable.
The intensity of the field in Eq. (1) is given byEq. (4) with respect to the angular variable:
Using CCD-array-aided measurements of the light field intensity I(r, φ, z) and summing it up on 2N radii rm, Eq. (5) is reduced to a set of 2N × 2N linear algebraic equations:Eq. (2) or (3), and xn = |Cn|2 are the sought-for squared modules of the expansion coefficients in Eq. (1). Solving the system (6) by, for example, inverse matrix method, we can derive the square modules of all coefficients in Eq. (1). Note that for the system (6) to have an unambiguous solution, the radial basis functions in Eq. (1) need to meet the following condition: . This condition can be satisfied by choosing the basis functions in Eq. (1), for instance, in the form:
For the basis functions given by Eq. (2) and (3), the system (6) can be unambiguously solved at n ∈ [–N, 0] or n ∈ [0, N].
Let us now derive a relationship for determining the OAM of the paraxial light field (1). May the on-axis projection Jz of the total OAM vector and the total light field energy W are given by the well-known relations :
Substituting (1) into (7) and (8) yields the normalized total OAM of the beam (1):Equation (9) shows that for the OAM to be derived it will suffice to measure just squared modules of the coefficients in Eq. (1), without knowing their phases. The squared modules of the coefficients of Eq. (1) can be derived by solving the set of Eqs. (6): xn = |Cn|2, with the powers In of all constituent beams calculated using the known expansion basis: Eq. (2) or Eq. (3), with regard to the above considerations.
2.1 Calculating fractional OAM of superposition of LG modes
Next, using Eq. (9) we can engineer a light field with the desired OAM value. For instance, a superposition of Laguerre-Gauss modes of Eq. (2) with orders (0,2) and (0, 4) can be selected, producing the normalized OAM equal to 3 and 3.5, thus generating the field in the form:
By introducing the factors 1/(πw02)1/2 and 1/(12πw02)1/2 the energy of both modes is equated. Equation (9) suggests that at C2 = 1, C4 = 1, we have Jz/W = 3, and at C2 = 1, C4 = 31/2 we have Jz/W = 3.5.
The numerical simulation also relies on the Fresnel transform. The simulation parameters are: the wavelength, λ = 532 nm, the waist radius, w0 = 1 mm, computation domain: –R ≤ x, y ≤ R, (R = 5 mm), and the increment of the coordinates, Δx = Δy = 20 µm.
Figure 1 depicts distributions of intensity [Fig. 1(a,c)] and phase [Fig. 1(b,d)] in the initial plane z = 0, for superposition of the LG modes with the normalized OAM equal to 3 [Fig. 1(a,b)] and 3.5 [Fig. 1(c,d)].
For the beams shown in Fig. 1, OAM was also calculated using Eqs. (6) and (9), assuming that the beam contained five modes with the topological charges ranging from m = 1 to m = 5. The radii rm were chosen in the form mΔr (m = 1, ..., 5), thus providing that 80% of the incident beam energy were found within a radius of 5Δr [Fig. 2(a)], i.e. . In our case 5Δr ≈1.7w0 in Fig. 1(a,b) and 5Δr ≈1.8w0 in Fig. 1(c,d).
For the beam in Fig. 1(a, b) the superposition (1) was calculated for the coefficients C1 = 0; C2 = 1; C3 = 0; C4 = 1; C5 = 0. By solving the set of Eqs. (6), the squared modules of the coefficients were found to be |C1|2 = 0.001; |C2|2 = 0.98; |C3|2 = 0.05; |C4|2 = 0.94; C5 = 0.03. With the theoretical value of the normalized OAM being Jz/W = 3.00, the value calculated from Eq. (9) was 3.01.
For the beam in Fig. 1(c, d) the superposition (1) was calculated for the coefficients C1 = 0; C2 = 1; C3 = 0; C4 = 31/2; C5 = 0. With the radii rm in Fig. 2(b) calculated in a similar way, the solution of Eqs. (6) the squared modules of the coefficients were found to be |C1|2 = −0.002; |C2|2 = 1.03; |C3|2 = −0.09; |C4|2 = 3.12; and |C5|2 = −0.06. With the theoretical value of the beam normalized OAM being Jz/W = 3.50, Eq. (9) gives a value of 3.51. The error in determining the OAM using the numerical solution of the system in Eq. (6) is less than 1%.
2.2 Experimental determining OAM of superposition of LG modes
Figure 3 shows an experimental setup for determining the OAM of various superpositions of the LG modes. Collimated and expanded laser beam from a solid state laser illuminated the display of a spatial light modulator SLM HOLOEYE PLUTO VIS (1920 × 1080 pixels, each pixel of 8 μm). The phase mask on the modulator was obtained by encoding the complex amplitudes distribution of a given superposition of the LG modes. The reflected and modulated laser beam was spatially filtered by lenses L2 (f2 = 500 mm) and L3 (f3 = 150 mm) and a diaphragm D, which blocked the light unmodulated by the modulator. Intensity distributions of the generated superpositions of the LG beams were recorded by the camera placed in the focal plane of the lens L2 (the plane conjugate to the plane of the modulator display z = 0). Figure 4 shows examples of the generated intensity distributions for different superpositions of the LG modes. Using Eqs. (6) and (9). OAM of such beams were calculated. It was supposed that the beam consists of 5 modes with their topological charges from m = 1 to m = 5. The Gaussian beam had the waist radius of ω0 = 190 μm.
For the beam in Fig. 4(a), the coefficients in the superposition in Eq. (1) are equal to C1 = 0; C2 = 1; C3 = 0; C4 = 1; C5 = 0. Solution of the system in Eq. (6) gives the following values: C1 = −0.31i; C2 = 0.95; C3 = −0.71i; C4 = 1.00; C5 = −0.21i. Theoretical value of the normalized OAM of such a superposition of LG modes is Jz/W = 3.00. Calculation by Eq. (9) using the experimentally obtained intensity distribution gives the value 3.17.
For the beam in Fig. 4(b), the coefficients in the superposition in Eq. (1) are equal to C1 = 0; C2 = 1; C3 = 0; C4 = 1/31/2; C5 = 0. Solution of the system in Eq. (6) gives the following values: C1 = −0.31i; C2 = 1.00; C3 = −0.77i; C4 = 0.75; C5 = 0.19i. Theoretical value of the normalized OAM of such a superposition of LG modes is Jz/W = 2.50. Calculation by Eq. (9) using the experimentally obtained intensity distribution gives the value 2.81.
For the beam in Fig. 4(c), the coefficients in the superposition in Eq. (1) are equal to C1 = 0; C2 = 1; C3 = 0; C4 = 0; C5 = 1/21/2. Solution of the system in Eq. (6) gives the following values: C1 = −0.24i; C2 = 1.00; C3 = −0.53i; C4 = −0.23i; C5 = 0.78. Theoretical value of the normalized OAM of such a superposition of LG modes is Jz/W = 3.00. Calculation by Eq. (9) using the experimentally obtained intensity distribution gives the value 3.23.
For the beam in Fig. 4(d), the coefficients in the superposition in Eq. (1) are equal to C1 = 1; C2 = 0; C3 = 0; C4 = 0; C5 = 1. Solution of the system in Eq. (6) gives the following values: C1 = 1.00; C2 = 0.40; C3 = −0.48i; C4 = 0.32; C5 = 0.97. Theoretical value of the normalized OAM of such a superposition of LG modes is Jz/W = 3.00. Calculation by Eq. (9) using the experimentally obtained intensity distribution gives the value 2.91.
Table 1 shows that for a superposition of small number of LG modes the OAM measurement error is about 10%. Here we investigated a simple method for determining the OAM by using a single intensity distribution. In principle, no any optical elements are required for this method. Below we study more accurate and more universal method for determining the OAM, though using this method requires measuring two intensity distributions.
3. Determining OAM from intensity measurements in focus of two cylindrical lenses
In Refs [13,14]. it was proposed that OAM should be determined from a single measurement of the intensity in the focus of a cylindrical lens. But this method is only suited for optical vortices with a circularly symmetric radial component of the beam, whose amplitude is given byEq. (10) replaced by a real-valued μ:
Below, we describe a simple generalization of the above technique allowing OAM of an arbitrary paraxial beam to be determined from two intensity distribution measurements in the Fourier planes of two cylindrical lenses positioned perpendicularly to each other and placed in identical beam branches outgoing from a splitting cube (Fig. 5).
In the Cartesian coordinates, Eq. (7) can be rewritten as
Substituting in (13) the complex amplitude E(x, y, z) expressed as a 1D Fourier transform separately for each coordinate:Eqs. (14)-(16)) are 1D Fourier transforms in the Cartesian coordinates y and x. In Eq. (16), the coordinates x, ξ, y, η are taken along the Cartesian axes in the focal plane of the cylindrical lenses. Note that with the Fourier transform of Eqs. (14) and (15) not taken along the x- and y-coordinates, these remain the same as in the cylindrical lenses' plane. Thus, OAM can be determined using either Eq. (16) via measuring two intensity distributions in the focal planes of two cylindrical lenses, or Eqs. (6) and (9) via measuring a single intensity distribution. The advantage of the approach based on Eq. (16) over that based on Eqs. (6) and (9) is that the former is suited for arbitrary paraxial light fields and not only for light fields described by superposition of a finite number of LG modes (n, 0) of Eq. (2) or BG modes of Eq. (3). With the relation (16) being approximate, for the derivation of the exact relation see the Appendix section. The exact formula (28) does not allow the OAM to be calculated from the experimental data. Meanwhile, the approximate relation (16) offers a near accurate result (with a 1 - 4% error) for the paraxial beams.
3.1 Paraxial laser beams with fractional OAM
While optical vortices can carry both integer and fractional OAM, the fractional OAM vortices can be of different nature. For example, a detailed analysis of the structure of an optical vortex with fractional topological charge was conducted in , Eq. (12). An exact formula to describe the fractional OAM of such vortices was derived in Ref . For non-integer µ, such an optical vortex was shown to carry fractional OAM and be devoid of an isolated on-axis intensity null (singularity), having instead a chain of isolated intensity nulls with topological charges + 1 and −1 near a transverse coordinate axis. Upon free-space propagation, these isolated intensity nulls can be mutually destroyed . Similar behavior is characteristic of asymmetric Bessel beams , Bessel-Gauss beams, and Lommel modes . All these beams carry fractional OAM and have a chain of isolated intensity nulls with topological charges + 1 and −1 in the transverse plane.
On the other hand, the fractional-OAM beam can be generated as a linear combination of a finite number of conventional Laguerre-Gauss modes (Eqs. (1), (2)) or Bessel-Gauss modes (Eq. (3)). In this case, as it propagates, the beam can retain its structure (with due account of scale and rotation). Such beams have a finite number of isolated intensity nulls equal to the number of terms N in Eq. (1).
Using quantum formalism, a relationship to describe the fractional OAM of the beam (12) has been derived . Below, we derive a similar relationship for conventional, rather than quantum, vortices. A fractional-OAM optical vortex can be expanded in terms of integer optical vortices:
Substituting the right-hand side of Eq. (17) into (7) yields:Eq. (8). The series in the right-hand side of Eq. (18) can be reduced to a reference series 27]:
3.2 Gaussian optical vortex with fractional topological charge
Let us now derive a relationship to describe the complex amplitude in the Fresnel zone for a Gaussian beam having passed through a spiral phase plate (SPP) with fractional topological charge. In the initial plane, the beam amplitude is given by
Then, using a Fresnel transform, we find that at a distance z from the initial plane the complex amplitude of the field (21) takes the form:Eq. (22) has been derived using a reference integral :
From Eq. (22), an intensity null is seen to occur on the optical axis (ρ = 0), while the entire intensity pattern is not circularly symmetric. With the modified Bessel functions in Eq. (22) described by a complex argument, the argument takes real values only in the far-field diffraction zone (z >> z0). The OAM of the field of Eq. (22) is the same as that of the initial field of Eq. (21) and is equal to Eq. (20).
3.3 Bessel optical vortex carrying fractional topological charge
Below, we discuss another example of an optical vortex carrying a fractional topological charge. May a SPP be illuminated by a narrow annular beam whose complex amplitude in the initial plane isEq. (20). It is worth noting that the field (26) in general has no on-axis (ρ = 0) intensity null. The intensity null in Eq. (26) occurs only when µ is integer. Actually, for non-integer µ, in Eq. (26) the only non-zero on-axis term (ρ = 0) is the term J0(0) = 1. If, however, µ is integer this term also equals zero, because the sine function in Eq. (26) takes a zero value.
An optical vortex with fractional topological charge in Eq. (12) is described by superposition of a countable number of optical vortices with integer topological charge. And vice versa, the light field (1) is composed of a finite number of optical vortices with integer topological charge. For the field (1) to carry fractional OAM it would suffice to keep just two terms. Actually, if in Eq. (1) the power of the constituent beams in Eq. (1) is normalized to unit magnitude, , and the sum of squared modules of the coefficients is set equal to unity: , Eq. (9) for two terms with topological charges n and m takes the form (W = 1):
From Eq. (27), it is seen that if a is an irrational number smaller than unity it would suffice to superimpose two optical vortices with integer topological charge to generate a fractional-OAM light beam.
3.4 Calculating OAM by using intensity distributions in foci of two cylindrical lenses
In this section, we numerically verify proposed technique for determining OAM based on Eq. (20), using as illustration the field (22), which is generated with the aid of the phase μφ formed on the liquid-crystal display of a SLM.
The modeling is based on the numerical calculation of a Fresnel transform. The simulation parameters are: wavelength of light, λ = 532 nm, Gaussian beam waist, w0 = 1 mm, focal length of the cylindrical lens, f = 1 m, calculation domain: –R ≤ x, y ≤ R, (R = 5 mm), z = f, and discretization step on both coordinates, Δx = Δy = 20 µm.
Shown in Fig. 6 are phase patterns in the initial plane z = 0 [Fig. 6(a,d,g,k)] and intensity patterns in the focus of the lens whose axis is parallel to the x-axis [Fig. 6(b,e,h,l)] and y-axis [Fig. 6(c,f,i,m)] for a Gaussian beam having passed through a SPP with the topological charge 3.00 [Fig. 6(a-c)], 3.25 [Fig. 6(d-f)], 3.50 [Fig. 6(g-i)], and 3.75 [Fig. 6(j-l)]. The intensity patterns in the focal plane of the cylindrical lens [Fig. 6(b,c,e,f,h,i,k,l)] are shown in the domain –R/2 ≤ x, y ≤ R/2.
According to Eq. (20), for the above-indicated values of the topological charge, the normalized OAM should be equal to 3.00; 3.09, 3.50, and 3.91, respectively. Actually, by substituting in Eq. (16) intensity values corresponding to the intensity distribution in Fig. 6, we obtain the following OAM values: 2.98 (1.50 + 1.48), 3.06 (1.46 + 1.60), 3.45 (1.70 + 1.75), and 3.87 (2.00 + 1.87). The respective error is 0.7%, 1%, 1.4%, and 1%. Hereafter, OAM is presented as a two-term sum in compliance with Eq. (16). This serves to demonstrate a different, albeit similar, contribution of these terms into OAM. In Refs [13,14]. the contribution of the two terms in Eq. (16) into OAM was assumed to be the same in magnitude and opposite in sign.
In a similar way, measurement of the large-valued fractional OAM can be numerically simulated using a Fresnel integral. Considering that in this case the optical vortex has a stronger divergence the computation domain R needs to be increased, while the increments Δx and Δy need to be made smaller. The simulation parameters are: wavelength, λ = 532 nm, waist radius, w0 = 1 mm, focal length of the cylindrical lens, f = 1 m, computation domain, –R ≤ x, y ≤ R, (R = 10 mm), z = f, the increment on both coordinates, Δx = Δy = 10 µm.
Figure 7 shows phase patterns [Fig. 7(a)] in the initial plane z = 0 and intensity patterns in the focus of the cylindrical lens whose axis is parallel to the x-axis [Fig. 7(b)] and y-axis [Fig. 7(c)] for a Gaussian beam having passed through a SPP with the topological charge 30.3.
According to Eq. (20), for the topological charge 30.3, the normalized OAM should be equal to 30.15. Substituting the intensity values from Fig. 7 into Eq. (16), we find that OAM is 28.87 (14.31 + 14.56), with the error amounting to 4%.
When determining the OAM of superposition of two LG modes, the error is essentially lower because such beams contain just two angular harmonics, unlike beams generated by a SPP with fractional topological charge, containing an infinite number of angular harmonics. Plots in Fig. 8 show in which way the normalized power of individual harmonics (OAM-spectrum) depends on the topological charge of (a) a Gaussian beam having passed through the SPP (12) with fractional topological charge μ = 3.5 and (b) superposition of the LG modes (10) with indices (0,2) and (0,4), whose normalized OAM also equals 3.5.
From Fig. 8 it is seen that while having different angular spectrum, beams with integer topological charge can carry the same OAM. The fewer angular harmonics in the beam, the more accurately is the OAM determined by measuring two intensities in the foci of the cylindrical lenses. The proposed technique offers a 4% accuracy when measuring a fractional OAM under 30, with the accuracy depending on the size of the computational domain. For instance, the intensity in Fig. 7 was calculated in a 10-mm × 10-mm domain. By properly increasing the computational domain, the error of OAM determination can be further decreased. However, in practice the accuracy cannot be improved in this way as when measuring low intensities under 0.01 of the maximum, the accuracy is not improved due to decreasing signal-to-noise ratio.
Now we confirm our method for asymmetric crescent-shaped Laguerre-Gaussian beams . Complex amplitude of these beams in the initial plane reads as E(x, y, z = 0) = [(x – aw0) + (y – iaw0)]nexp(–s/w02)Lmn(2s/w02), where (x, y, z) are the Cartesian coordinates, Lmn(x) is the associated Laguerre polynomial, w0 is the Gaussian beam waist radius, a is the asymmetry parameter, n is the vortex topological charge, m defines the number of light rings (or arcs), s = (x – aw0)2 + (y – iaw0)2. Figures 9(a,b) show intensity distributions of the symmetric (a = 0) and asymmetric (a = 0.25) Laguerre-Gaussian beams. We used the following parameters for calculations: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1mm, Laguerre polynomial indices (m, n) = (3, 5), calculation area -R ≤ x, y ≤ R (R = 10 mm) (4096 × 4096 points). Figures 9(c, d) show intensity patterns in the focus of the cylindrical lens whose axis is parallel to the x-axis [Fig. 9(c)] and y-axis [Fig. 9(d)]. Using these intensity patterns, the OAM was calculated. Similarly, the OAM was computed for several other values of the asymmetry parameter. Figure 9(e) shows the theoretical dependence of the OAM on the asymmetry parameter (black curve) (Eq. (22) in ) and the values computed via Eqs. (14)-(16) by using the intensity distributions in the foci of the cylindrical lenses (red circles).
3.5 Experimental determination of OAM by measuring two intensity distributions
The experimental study was conducted using an optical setup in Fig. 5. An initial Gaussian laser beam of wavelength 532 nm was expanded and collimated with a 40-µm circular pinhole PH. The collimated light beam was normally incident on a transmission SLM (a 1024 × 768-pixel HOLOEYE LC 2012, with 36-µm pixels and a 36.9 × 27.6-mm operating region). Spherical lenses of focal lengths 350, 250, and 150 mm were utilized as lenses L1, L2, and L3. Using an optical configuration composed of lenses L2, L3 and a spatial filter F, the light field from the SLM was filtered, thus cutting off an inoperative zero diffraction order. The mutually perpendicular cylindrical lenses CL1 and CL2 had a 500-mm focal length. Note that the use of long-focus cylindrical lenses in the numerical simulation and the experiment stems from the requirement of the paraxial focusing. The said lenses were imaged onto the SLM-display, with the Gaussian beam waist being about 6.5 mm. The intensity pattern in the focus of the cylindrical lenses CL1 and CL2 was recorded with 3264 × 2448-pixel CCD-cameras CCD1 and CCD2 (with a 1.67 µm × 1.67 µm pixel).
Figure 10 depicts intensity distributions obtained in the foci of the lenses CL1 and CL2 by the SLM-aided optical vortices at different values of the topological charge (3; 3.25; 3.5; 3.75, and 30.3). In this case, respective values of the normalized OAM derived from Eq. (16) based on the intensity measurements (Fig. 10) were found to be 3.02; 3.10; 3.47; 3.96, and 27.78. These results are in good agreement with the numerical simulation results.
For comparison, Table 2 gives OAM values numerically and experimentally derived using two cylindrical lenses. For the beams outgoing from a SPP with a fractional topological charge of 3.00; 3.25; 3.50; 3.75; 30.3, Table 2 gives theoretical OAM values derived from Eq. (20), alongside values obtained from Eq. (16) using numerically and experimentally derived intensity distributions.
According to Table 2, for the topological charges less than 30, the experimental error in determining the OAM of a fractional-topological-charge laser beam by measuring two intensity distributions in the foci of two cylindrical lenses rotated 90 degrees relative to each other does not exceed 8%.
Figure 11 shows theoretical dependence of the OAM on the topological charge m as well as several experimental OAM values. According to the experiments, the measurement error increases with the topological charge. Maximal error was nearly 10% for m = 34.1, while for m < 20 the error was several percents.
In this work, we have proposed and studied numerically and experimentally two simple and efficient approaches to determining OAM of paraxial laser beams. With the first one, the intensity measured in the Fresnel zone is then numerically averaged over angle at discrete radii before deriving squared modules of the light field expansion coefficients via solving a linear set of equations. With the other technique, two intensity distributions are measured in the Fourier plane of a pair of cylindrical lenses with their axes positioned perpendicularly, before calculating the first-order moments of the measured intensities. The experimental error grows almost linearly from about 1% when determining small fractional OAM (up to 4), and to about 10% for large fractional OAM (up to 34).
5 Appendix 1
Considering that cylindrical lenses perform a Fresnel transform along one Cartesian coordinate, with free space performing a Fresnel transform along the second Cartesian coordinate, Eqs. (14) and (15) are rearranged to
Then, OAM takes the form:
Taking the integrals in (30), we obtain:
Or, in a different form:
The difference between Eqs. (32) and (16), which describe the OAM, is in the third and fourth terms that cannot be derived experimentally. However, from the numerical simulation and experiment, for paraxially propagating beams, the contribution of these two terms to the OAM is seen to be insignificant (several percent).
RF Ministry of Science and Higher Education (State assignment); Russian Science Foundation (17-19-01186); Russian Foundation for Basic Research (18-29-20003).
Section 2 “Determining the OAM from a single measurement of the intensity distribution” was supported by the Russian Science Foundation (17-19-01186), Section 3 “Determining OAM from intensity measurements in the focus of two cylindrical lenses” was supported by the Russian Foundation for Basic Research (18-29-20003), Appendix 1 was supported by RF Ministry of Science and Higher Education (State assignment).
1. J. Wang, “Advances in communications using optical vortices,” Photon. Res. 4(5), B14–B28 (2016). [CrossRef]
2. Z. Wang, Y. Yan, A. Arbabi, G. Xie, C. Liu, Z. Zhao, Y. Ren, L. Li, N. Ahmed, A. J. Willner, E. Arbabi, A. Faraon, R. Bock, S. Ashrafi, M. Tur, and A. E. Willner, “Orbital angular momentum beams generated by passive dielectric phase masks and their performance in a communication link,” Opt. Lett. 42(14), 2746–2749 (2017). [CrossRef] [PubMed]
3. X. Wang, Z. Nie, Y. Liang, J. Wang, T. Li, and B. Jia, “Recent advances on optical vortex generation,” Nanophotonics 7(9), 1533–1556 (2018). [CrossRef]
4. G. J. Gbur, Singular Optics (Boca Raton: CRC Press, 545 p, 2016).
5. V. V. Kotlyar, A. A. Kovalev, and A. P. Porfirev, Vortex laser beams (Boca Raton: CRC Press, 404 p, 2018).
6. S. N. Khonina, V. V. Kotlyar, R. V. Skidanov, V. A. Soifer, P. Laakkonen, J. Turunen, and Y. Wang, “Experimental selection of spatial Gauss-Laguerre Modes,” Opt. Mem. Neural. Networks 9(1), 73–82 (2000).
7. S. N. Khonina, V. V. Kotlyar, V. A. Soifer, P. Paakkonen, J. Simonen, and J. Turunen, “An analysis of the angular momentum of a light field in terms of angular harmonics,” J. Mod. Opt. 48(10), 1543–1557 (2001). [CrossRef]
8. G. Ruffato, M. Massari, G. Parisi, and F. Romanato, “Test of mode-division multiplexing and demultiplexing in free-space with diffractive transformation optics,” Opt. Express 25(7), 7859–7868 (2017). [CrossRef] [PubMed]
9. Y. Li, X. Li, L. W. Chen, M. B. Pu, J. J. Jin, M. H. Hong, and X. G. Luo, “Orbital angular momentum multiplexing and demultiplexing by a single metasurface,” Adv. Opt. Mater. 5(2), 1600502 (2017). [CrossRef]
10. A. Ya. Bekshaev, M. S. Soskin, and M. V. Vasnetsov, “Optical vortex symmetry breakdown and decomposition of the orbital angular momentum of light beams,” J. Opt. Soc. Am. A 20(8), 1635–1643 (2003). [CrossRef] [PubMed]
11. V. Denisenko, V. Shvedov, A. S. Desyatnikov, D. N. Neshev, W. Krolikowski, A. Volyar, M. Soskin, and Y. S. Kivshar, “Determination of topological charges of polychromatic optical vortices,” Opt. Express 17(26), 23374–23379 (2009). [CrossRef] [PubMed]
13. S. N. Alperin, R. D. Niederriter, J. T. Gopinath, and M. E. Siemens, “Quantitative measurement of the orbital angular momentum of light with a single, stationary lens,” Opt. Lett. 41(21), 5019–5022 (2016). [CrossRef] [PubMed]
15. Z. Liu, S. Gao, W. Xiao, J. Yang, X. Huang, Y. Feng, J. Li, W. Liu, and Z. Li, “Measuring high-order optical orbital angular momentum with a hyperbolic gradually changing period pure-phase grating,” Opt. Lett. 43(13), 3076–3079 (2018). [CrossRef] [PubMed]
17. A. D’Errico, R. D’Amelio, B. Piccirillo, F. Cardano, and L. Marrucci, “Measuring the complex orbital angular momentum spectrum and spatial mode decomposition of structured light beams,” Optica 4(11), 1350–1357 (2017). [CrossRef]
19. L. A. Melo, A. J. Jesus-Silva, S. Chávez-Cerda, P. H. S. Ribeiro, and W. C. Soares, “Direct measurement of the topological charge in elliptical beams using diffraction by a triangular aperture,” Sci. Rep. 8(1), 6370 (2018). [CrossRef] [PubMed]
20. H. Gao, Y. Han, Y. Li, D. Zhu, M. Sun, and S. Yu, “Topological charge measurement of concentric OAM states using the phase-shift method,” J. Opt. Soc. Am. A 35(1), A40–A44 (2018). [CrossRef] [PubMed]
21. A. V. Volyar, M. V. Bretsko, Ya. E. Akimova, and Yu. A. Egorov, “Beyond the light intensity or intensity moments and measurements of the vortex spectrum in complex light beams,” Comput. Opt. 42(5), 736–743 (2018). [CrossRef]
22. A. Volyar, M. Bretsko, Y. Akimova, and Y. Egorov, “Measurement of the vortex spectrum in a vortex-beam array without cuts and gluing of the wavefront,” Opt. Lett. 43(22), 5635–5638 (2018). [CrossRef] [PubMed]
23. A. E. Siegman, Lasers (Mill Valley, CA: University Science Books, 1986).
24. Z. Xie, S. Gao, T. Lei, S. Feng, Y. Zhang, F. Li, J. Zhang, Z. Li, and X. Yuan, “Integrated (de)multiplexer for orbital angular momentum fiber communication,” Photon. Res. 6(7), 743–749 (2018). [CrossRef]
25. F. Gori, G. Guattari, and C. Padovani, “Bessel-Gauss beams,” Opt. Commun. 64(6), 491–495 (1987). [CrossRef]
26. M. V. Berry, “Optical vortices evolving from helicoidal integer and fractional phase steps,” J. Opt. A, Pure Appl. Opt. 6(2), 259–268 (2004). [CrossRef]
27. J. B. Götte, S. Franke-Arnold, R. Zambrini, and S. M. Barnett, “Quantum formulation of fractional orbital angular momentum,” J. Mod. Opt. 54(12), 1723–1738 (2007). [CrossRef]
29. A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev, Integrals and series. Volume 2: Special functions (New York: Gordon and Breach, 1986).
30. A. A. Kovalev, V. V. Kotlyar, and A. P. Porfirev, “Asymmetric Laguerre-Gaussian beams,” Phys. Rev. A (Coll. Park) 93(6), 063858 (2016). [CrossRef]