ABCD-treatment of a propagating doughnut beam generated by a spiral phase plate

We apply the Collins-Huygens integral to analytically describe propagation of a doughnut beam generated by a spiral phase plate. Measured beam profiles in free space and through an ABCD-lens system illustrate excellent agreement with theory. Applications range from the creation of optical beams with angular momentum to microscopy to trapping neutral atoms. The method extends to other beam shaping components, too. © 2011 Optical Society of America OCIS codes: (050.4865) Optical vortices; (070.2580) Paraxial wave optics; (070.2590) ABCD transforms; (080.2730) Matrix methods in paraxial optics; (080.4865) Optical vortices References and links 1. T. A. Klar, S. Hell, “Subdiffraction resolution in far-field fluorescence microscopy,” Opt. Lett. 24, 954 (1999) 2. T. Kuga et al., “Novel Optical Trap of Atoms with a Doughnut Beam,” Phys. Rev. Lett. 78, 4713 (1997) 3. A. Kaplan, N. Friedman, N. Davidson, “Optimized single-beam dark optical trap,” J. Opt. Soc. Am. B 19, 1233 (2002) 4. P. Xu, X. He, J. Wang, M. Zhan, “Trapping a single atom in a blue detuned optical bottle beam trap”, Opt. Lett. 35, 2164 (2010) 5. D. W. Zhang, X.-C. Yuan, ”Optical doughnut for optical tweezers,” Opt. Lett. 28, 740 (2003) 6. L. Allen, S. M. Barnett, M. J. Padgett, Optical Angular Momentum, Taylor and Francis, (2003) 7. M.R. Dennis, K. O’Holleran, M.J. Padgett, in “Progress in Optics,” 53, 293-363 (2009) 8. V. Y. Bazhenov, M. S. Soskin, M . V . Vasnetsov, “Screw dislocations in light wavefronts”, J. Mod. Opt. 39, 985 (1992) 9. N. R. Heckenberg, R. McDuff, C. P. Smith, A. G. White, “Generation of optical phase singularities by computergenerated holograms,” Opt. Lett. 17, 221 (1992) 10. E. Abramochkin, V. Volostnikov, “Beam transformations and nontransformed beams,” Opt. Commun. 83, 123 (1991) 11. M. Padgett, J. Arlt, N. Simpson, L. Allen, “Optical tweezers and optical spanners with Laguerre-Gaussian modes,”, Am. J. Phys. 64, 77 (1996) 12. S. N. Khonina, V. V. Kotlyar, M . V. Shinkaryev, V. A. Soifer, G . V. Uspleniev, “The phase rotor filter,” J. Mod. Opt. 39, 1147 (1992) 13. M.W. Beijersbergen, et. al., “Helical-wavefront laser beams produced with a spiral phaseplate,” Opt. Commun. 112, 321 (1994) 14. J. H. Goodman, Introduction to Fourier Optics (McGrawHill, New York, 1968) 15. D. Rozas, C. T. Law, G. A. Swartzlander, “Propagation Dynamics in Optical Vortices,” J. Opt. Soc. Am. B 14, 3054 (1997) 16. Z. S. Sacks, D. Rozas, G. A. Swartzlander, “Holographic formation of optical-vortex filaments,” J. Opt. Soc. Am. 15, (1998) 17. S. A. Collins, Jr., “Lens-System Diffraction Integral Written in Terms of Matrix Optics,” J. Opt. Soc. Am. 60, 1168 (1970); A. Siegman, Lasers, University Science Books (Mill Valley 1986), p. 781 18. I. S. Gradshteyn, I. M. Rhyzhik, Table of Integrals, Series, and Products, Academic Press (San Diego 1979) #151920 $15.00 USD Received 29 Jul 2011; revised 16 Sep 2011; accepted 18 Sep 2011; published 10 Oct 2011 (C) 2011 OSA 24 October 2011 / Vol. 19, No. 22 / OPTICS EXPRESS 21205


Introduction
Optical beams with a dark spot at the center are called doughnut beams.They have proven useful for applications in high-precision imaging of biological samples [1], for trapping dark field seeking neutral atoms in blue detuned dipole traps [2,3,4], or with optical tweezers [5].In recent years, also light beams carrying orbital angular momentum exhibiting a natural dark spot at their center have received increased interest [6,7].Methods for generating doughnut shaped beams include the use of computer-generated holograms [8,9], cylindrical lens mode converters [10,11] and the use of spiral phase plates (SPP) [12,13].
Here we focus on doughnut beams generated by transforming a Gaussian beam with an SPP.From a practical point of view an analytic description compatible with the standard ABCDformalism remains most favorable for the treatment of beam propagation.Since the transformed beam is not an eigen solution of the paraxial wave equation, a more elaborate treatment e.g. in terms of the paraxial Fresnel integral [14] is required.It was first studied theoretically in [15], and Fresnel type beam patterns were analyzed in [16].The equivalent of the paraxial Fresnel integral in terms of the ABCD formalism is the Collins-Huygens integral [17] which we explicitly calculate.We furthermore illustrate its compatibility with the ABCD formalism by profile measurements of a Gaussian beam transformed by an SPP element.

Propagation dynamics of a doughnut beam
SPP components are commercially available for any optical wavelength.They consist of a cylindrical dielectric element with optical thickness increasing linearly with the azimuth angle φ .The SPP imprints a phase factor exp(i φ ) onto incident beams, where in applications the most frequently used winding number or topological charge is = 1.In the far field, cylindrically symmetric beams such as Gaussian beams are transformed into doughnut beams since all beam elements on opposing sides of the axis contribute with opposite sign and hence interfere destructively on the axis.For an incident centered Gaussian TEM 00 beam centered on the SPP with negligible thickness and winding number at position z (Fig. 1) the field amplitude with maximal amplitude E 0 directly behind the SPP takes for cylindrical coordinates (ρ , φ , z ) the form The wave number is k = 2π/λ and w 0 and z 0 = πw 2 0 /λ designate the 1/e-waist and the half-Rayleigh length of the incoming Gaussian beam, respectively.We have R(z , and η(z ) = tan −1 (z /z 0 ) the Gouy phase.For a conventional Gaussian beam (i.e. for = 0) the field distribution E(ρ, φ , z) can be calculated once the parameters R(z) and w(z) at position z are evaluated in terms of ABCD-matrices involving all components acting on the propagating optical beam in between z and z.For a Gaussian beam transformed by a spiral phase plate (SPP) ( ≥1) this no longer the case.Collins has shown [17], however, that the ABCD formalism remains very useful in the paraxial approximation: The Collins-Huygens integral (2) allows to calculate the field distribution E(ρ, φ , z) for any known initial distribution E(ρ , φ , z ) as a function of the global ABCD coefficients characterizing propagation from z to z. Substituting Eq. ( 1) into Eq.( 2), we find (3) For the sake of clarity we introduce abbreviations which represents a Gaussian TEM 00 -mode with yielding an effective Collins curvature radius R C , and a radius ρ C characterizing the extension of the vortex structure.Re-writing Eq. ( 3) we find Azimuthal integration over φ yields an -th order Bessel function of the first kind: The integral can be expressed in terms of modified Bessel functions I m of the first kind and m-th order using the formula [18 Identifying α = 1/R 2 C and β = 1/ρ C we find the analytic expression for the field amplitude propagated through the optical system from z to z With parameters Eq. (4-5) and the analytic expression Eq. 9 the propagation of a paraxial doughnut beam through any ABCD optical system is fully characterized.The beam profile can be calculated at any position, resulting in a convenient tool to analyze the propagation and focusing properties of doughnut beams generated by an SPP.Formula Eq. ( 9) is fully equivalent with the result derived first in [15].

Comparison with an experimental doughnut beam created by an SPP
We have experimentally investigated the propagation of a doughnut beam formed by an SPP for two important cases: propagation in free space and propagation through a lens system.A commercial SPP (RPC Photonics) produced by a lithographic technique imposes a phase factor exp(iφ ) with winding number = 1 on an incident beam with wavelength of 849.9 nm.
We have used a standard setup to measure the beam profile, see Fig. 2. We filter the mode profile by a single mode optical fiber in order to prepare a clean Gaussian TEM 00 beam.After the fiber the laser power is approximately 7.5 μW and the beam is linearly polarized.The beam is then passed through the SPP to generate the donut beam.The beam profile of the generated donut beam is measured using a Spiricon camera and analyzed using its software .A typical image of the generated donut beam detected by the beam profile camera is shown in Fig. 3. From the recorded images a radial intensity distribution is extracted from a single transverse cut showing the doughnut character of the beams.The theoretical description according to Eq. ( 9) compares well at the %-level with the measurement where small asymmetries can be traced to small misalignments.

Propagation in free space
In free space, propagation of the beam through the distance z is given by the ray transfer matrix for a distance z.The parameters (5) are Insertion into Eq.9 yields for free space propagation and = 1 The variation of the 1D radial intensity distributions of the doughnut beam with distance in free space from the SPP was measured by removing the lens system from Fig. 2. The results are shown in Fig. 4 and compared with the numerically calculated intensity distributions.They agree very well at all distances.The beam profile undergoes significant changes upon propagation since the distances cover a range from the Fresnel diffraction limit to the Fraunhofer limit.At the beginning the beam shows a high peak intensity and several radial fringes.As the beam propagates the peak intensity decreases and the fringes disappear towards large radii.

Lens system
In order to illustrate the propagation of the beam through a lens system, we have chosen the configuration shown in Fig. 2(b) which is used in another experiment in our group to generate tightly focused dipole traps for neutral atoms .The ABCD-transfer matrix of the lens system is obtained by taking the product of the transfer matrices of the individual optical elements, M tot = M free • M lens 3 • M free 3 • M tel • M free 1 As before, substituting the values of the matrix elements into Eq.( 9) and Eq. ( 5) yields the intensity distribution of the light field after passing the lens system.
The beam profile at different propagation distances behind the last focusing lens is shown in Fig. 5.In order to get a good fitting between the calculated and measured intensity distributions, the theoretical distances have been simultaneously offset by 15 mm.This adjustment is justified since we neither measured the exact position of the chip of the CCD camera nor the precise focal length of each lens.Again we find very good agreement of experimental and theoretical curves.The beam profile at the focal plane of the last lens cannot be measured since the beam diameter at this position is too small to be resolved by the CCD camera.

Conclusion
We have shown that the analytic solution Eq. 9, which is easily evaluated with a computer, makes the propagation of SPP generated doughnut beams accessible by the ABCD-method.We expect the procedure to be relevant for many applications including optical microscopy, neutral atom trapping, optical tweezers, and propagation of optical angular momentum beams.Our results indicate that the slightly forgotten Collins-Huygens integral promises useful applications beyond the present SPP elements, e.g for the half-phase plates used in [3].

Fig. 1 .
Fig.1.A Gaussian beam is converted at position z by a spiral phase plate (SPP).We are seeking the field distribution of the optical beam at position z in terms of the ABCDcoefficients of its trajectory.

Fig. 3 .
Fig. 3. (a): Beam profile at z = 250 mm for free space propagation.The 1D radial intensity distribution of the doughnut beam is taken along a straight line passing through the beam center (red line).(b): Radial intensity distribution extracted from the measurement.The theoretical curve (shaded area, shifted) shows all the details of the measured profile.

Fig. 4 .
Fig. 4. (a) Intensity distributions of the SPP generated doughnut beam propagating in free space.The distribution at 250 mm is identical with Fig. 3. (b) One-dimensional radial intensity distributions with measured (black) and calculated (red) intensity distributions.

Fig. 5 .
Fig. 5. (a) Intensity distributions of the SPP generated propagating doughnut beam transformed by the lens system at different propagation distances.(b) One-dimensional radial intensity distributions with measured (black) and calculated (red) intensity distributions.