Divergence of an orbital-angular-momentum-carrying beam upon propagation Miles Padgett1, Filippo M. Miatto2, Martin Lavery1, Anton Zeilinger3;4, and Robert W. Boyd1;2;5 1School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom 2Dept. of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, Ontario, K1N 6N5 Canada 3Institute for Quantum Optics and Quantum Information (IQOQI), Austrian Academy of Sciences, Boltzmanngasse 3, A-1090 Vienna, Austria 4Vienna Center for Quantum Science and Technology, Faculty of Physics, University of Vienna, Boltzmanngasse 5, A-1090 Vienna, Austria and 5The Institute of Optics, University of Rochester, Rochester, New York 14627, USA (Dated: February 12, 2018) There is recent interest in the use of light beams carrying orbital angular momentum (OAM) for creating multiple channels within free-space optical communication systems. One limiting issue is that, for a given beam size at the transmitter, the beam divergence angle increases with increasing OAM, thus requiring a larger aperture at the receiving optical system if the efficiency of detection is to be maintained. Confusion exists as to whether this divergence scales linarly with, or with the square root of, the beam's OAM. We clarify how both these scaling laws are valid, depending upon whether it is the radius of the Gaussian beam waist or the rms intensity which is kept constant while varying the OAM. I. INTRODUCTION intensity of one of these modes, r(Imax), is given as [3] r ` Over the past 20 years there has been a growing in- r(Imax) = j jw(z); (2) terest in the orbital angular momentum (OAM) of light, 2 which is carried by any optical beam possessing helical p 2 2 1 2 where w(z) = w0 1 + z =zR and zR = 2 kw0 is the phase fronts [1]. For phase fronts described by an az- Rayleigh range (here k = 2π/λ). One notes that for imuthally dependent phase factor exp(i`φ), the OAM is high OAM values the size of the LG mode is much larger equivalent to `¯h per photon [2]. This OAM of a light than the beam waist. The restriction of our study to the beam can be derived in a number of ways, but for our paraxial regime also implies that w(z) λ and hence purposes it can be seen to arise from the local direction that β (see Eq. (1)) remains small for all radii at which of the Poynting vector. At all points within the beam the the beam has a non-negligible intensity. Poynting vector is perpendicular to the phase fronts and In terms of their complex amplitudes, LG beams with hence possesses an azimuthal momentum component [3]. different mode indices are orthogonal to each other and As confirmed by measurement [4], helical phase fronts this suggests their use in communication applications lead to the local Poynting vector being skewed with re- where the orthogonality of the modes creates indepen- spect to the beam axis by an angle β, which is a func- dent channels. A simple optical communication system tion of r, the radial distance from the beam axis and based on OAM was demonstrated over 10 years ago [5], k0 = 2π/λ: but recently the idea has been revisited with reports of both high data rates [6] and outdoor demonstrations [7]. ` β = j j ; (1) One issue of critical importance is the degree to which k0r the use of OAM-carrying beams require the use of larger aperture optics. Confusion exists as to whether the far- One analytically simple, and widely used, form of field divergence of a beam carrying OAM scales linearly helically-phased beams is the Laguerre-Gaussian (LG) or with the square root of its angular momentum. arXiv:1410.8722v1 [quant-ph] 31 Oct 2014 modes having a field amplitude given as the product of a Gaussian term with beam waist w0, (defined as the radius at which the electric field of the Gaussian term is II. DIVERGENCE OF A BEAM lower by a factor 1=e), a generalised Laguerre polynomial with indices ` and p, and the coordinate r raised to the The skew angle of the Pointing vector calculated with j j power ` (see Eq. (3)). respect to the beam axis might be thought to approx- j j In this work we restrict ourselves to the paraxial regime imate the angular divergence of the beam. However, and small numerical aperture where the size of any spa- an additional contribution to the divergence that has to tial features in the modes or their superpositions is sig- be considered is the normal diffractive spreading arising nificantly larger than the optical wavelength. For p = 0, from the finite beam diameter. Rather than the beam the LG modes are single-ringed annular modes with zero waist w0, or the radius of maximum intensity r(Imax), the on-axis intensity and an azimuthal phase term exp(i`φ), divergence of a light beam is governed by the standard unless also ` is zero, in which case one has the funda- deviation of its spatial distribution, rrms [8, 9]. The cal- mental Gaussian mode. The radius of maximum optical culation of these quantities is straightforward. We start 2 from the definition of a LG mode from which we find 3 that the intensity distribution for the lowest-order radial ) mode (p = 0) is given by 0 2 !2j`j 2 p2r 2r2 I`(r; φ, z) = exp ; w(z)2π ` ! w(z) −w(z)2 r rms j j r (I ) (3) 1 max R (units of w beam size which is normalized such that I`(r; φ, z) rdrdφ = 1. From this, we can calculate the radial position of the maximum of intensity, by solving the equation @I`=@r = 1 2 3 4 5 6 7 8 9 10 0, obtaining (as anticipated in eq. (2) [3]) | | r FIG. 1: A comparison between the radius in the pupil plane ` of maximum intensity and the rms radius of the spatial dis- r(Imax) = j jw(z); (4) 2 tribution of a p = 0 LG mode. For large `, the two radii get p asymptotically close, as r(Imax) is proportional to j`j, while However, as r(Imax) is zero for ` = 0, clearly one cannot r is proportional to pj`j + 1. generally use this quantity to derive the angular spread rms of the beam as it propagates. As has been reported pre- viously, a better quantity is the standard deviation of the spatial distribution of the beam, rrms, which is a func- field, the divergence angle of a p = 0 LG mode is given tion of z and it is given by the square root of the radial by variance of the intensity distribution [10]: r ` + 1 λ s Z 1 α` = j j : (11) 2 2 2πw0 rrms(z) = 2π r I`(r; φ, z) rdr (5) 0 Eq. (11) is a convenient form for expressing the ` r ` + 1 dependence of the beam divergence whilst holding w0 = j j w(z); (6) 2 constant, revealing the approximate square root scaling, see figure 2. Alternatively the divergence angle can be In the plane of the beam waist, this simplifies to written in terms of other quantities, for instance, using r Eq. (7), one can instead express the ` dependence of the ` + 1 beam divergence whilst holding r constant, as r (0) = j j w ; (7) rms rms 2 0 ` + 1 Fig. 1 shows the ` dependence of r(I ) and r (0) α` = j j ; (12) max rms k0rrms(0) in the pupil plane. One notes that Under this condition w0 itself becomes a function of ` and s − = r ` + 1 in fact scales as ( ` + 1) 1 2. In this case, corresponding rms = ; (8) j j j j to a fixed rrms, we see that the beam divergence scales r(Imax) ` j j linearly with `, see figure 2. which tends to 1 as ` increases, so these two radii be- come equal for large j`,j see Fig. 1. III. DISCUSSION Given the radius rrms(z) as a function of propagation distance, z, we can calculate the corresponding diver- gence angle of a p = 0 LG mode, α`, to be Typical of a practical system is that the aperture of the launch optics is fixed by design, with a limiting radius @rrms(z) R. Since the size of the LG beam can be much bigger α`(z) = arctan (9) @z than w0, if wishing to design a system incorporating a r ! wide range of different OAM modes with a constant beam ` + 1 w2 z 0 waist, then one must choose a waist w0 R. LG beams = arctan j j 2 : (10) 2 zR w(z) of a fixed beam waist are produced by a cylindrical-lens mode converter [15]. Using a mode converter of this type This divergence angle is a function of z, but it asymp- to produce the modes one would observe the square root totically reaches a limiting value (when the beam is suf- scaling of the far-field beam divergence (11). ficiently far from the waist). Within the paraxial regime, A more frequently used and flexible method of produc- for small α`, we can equate the tangent of the angle with ing beams carrying OAM is the illumination of a diffrac- the angle and for limz!1 z=w(z) = zR=w0, so in the far tive optical element with an expanded Gaussian beam 3 10 9 a scaling approximately proportional to ` . This linear 8 j j 7 scaling of the beam size has previously been noted within 6 ` +1 optical tweezers and microscope systems [12].