Light beams with fractional orbital angular momentum and their vortex structure

Light emerging from a spiral phase plate with a non-integer phase step has a complicated vortex structure and is unstable on propagation. We generate light carrying fractional orbital angular momentum (OAM) not with a phase step but by a synthesis of Laguerre-Gaussian modes. By limiting the number of different Gouy phases in the superposition we produce a light beam which is well characterised in terms of its propagation. We believe that their structural stability makes these beams ideal for quantum information processes utilising fractional OAM states. © 2008 Optical Society of America OCIS codes: (350.5030) Phase; (050.4865) Optical vortices; (260.6042) Singular optics; (090.1760) Computer Holography; (230.6120) Spatial light modulators; (270.0270) Quantum optics. References and links 1. M. V. Berry, “Much ado about nothing: optical distortion lines (phase singularities, zeros, and vortices),” in International Conference on Singular Optics, M. S. Soskin, ed., Proc. SPIE 3487, 1–5 (1998). http://spie. org/x648.xml?product id=317693. 2. M. V. Berry and M. R. Dennis, “Knotted and linked phase singularities in monochromatic waves,” Proc. Royal Society of London, Series A 457(2013), 2251–2263 (2001). http://www.journals.royalsoc.ac.uk/ content/yv6lqeufbq2phg9l. 3. J. Leach, M. R. Dennis, and M. J. Padgett, “Laser beams: Knotted threads of darkness,” Nature 432(165) (2004). http://www.nature.com/nature/journal/v432/n7014/abs/432165a.html. 4. J. Leach, M. R. Dennis, J. Courtial, and M. Padgett, “Vortex knots in light,” New J. Phys. 7, 55 (2005). http: //www.iop.org/EJ/abstract/1367-2630/7/1/055/. 5. L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. P. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian modes,” Phys. Rev. A 45, 8185–8190 (1992). http://link.aps. org/abstract/PRA/v45/p8185. 6. L. Allen and M. J. Padgett, “The Poynting vector in Laguerre-Gaussian beams and the interpretation of their angular momentum density,” Opt. Commun. 184, 67–71 (2000). 7. M. W. Beijersbergen, R. P. C. Coerwinkel, M. Kristensen, and J. P. Woerdman, “Helical-wavefrontnext term laser beams produced with a spiral phaseplate,” Opt. Commun. 112, 321–327 (1994). 8. V. Y. Bazhenov, M. V. Vastnetsov, and M. S. Soskin, “Laser-beams with screw dislocations in their wave-fronts,” JETP Lett. 52, 429–431 (1990). 9. M. V. Berry, “Optical vortices evolving from helicoidal integer and fractional phase steps,” J. Opt. A 6, 259–269 (2004). http://www.iop.org/EJ/abstract/1464-4258/6/2/018. 10. J. Leach, E. Yao, and M. J. Padgett, “Observation of the vortex structure of a non-integer vortex beam,” New J. Phys. 6, 71 (2004). http://www.iop.org/EJ/abstract/1367-2630/6/1/071/. 11. K. O’Holleran, M. R. Dennis, and M. J. Padgett, “Illustrations of optical vortices in three dimensions,” Journal of European Optical Society Rapid Publications 1, 06008 (2006). https://www.jeos.org/index.php/ jeos rp/article/view/06008. #88477 $15.00 USD Received 11 Oct 2007; revised 21 Dec 2007; accepted 23 Dec 2007; published 11 Jan 2008 (C) 2008 OSA 21 January 2008 / Vol. 16, No. 2 / OPTICS EXPRESS 993 12. J. Götte, S. Franke-Arnold, R. Zambrini, and S. M. Barnett, “Quantum formulation of fractional orbital angular momentum,” J. Mod. Opt. 54, 1723–1738 (2007). http://www.informaworld.com/smpp/ content∼content=a779773614. 13. S. S. R. Oemrawsingh, A. Aiello, E. R. Eliel, G. Nienhaus, and J. P. Woerdman, “How to observe HighDimensional Two-Photon Entanglement with Only Two Detectors,” Physical Review Letters 92, 217901 (2004). http://link.aps.org/abstract/PRL/v92/e217901. 14. G. F. Calvo, A. Picón, and A. Bramon, “Measuring two-photon orbital angular momentum entanglement,” Phys. Rev. A 75, 012319 (2007). http://link.aps.org/abstract/PRA/v75/e012319. 15. A. Aiello, S. S. R. Oemrawsingh, E. R. Eliel, and J. P. Woerdman, “Nonlocality of high-dimensional two-photon orbital angular momentum states,” Phys. Rev. A 72, 052114 (2005). http://link.aps.org/abstract/ PRA/v72/e052114. 16. S. M. Barnett and D. T. Pegg, “Quantum theory of rotation angles,” Phys. Rev. A 41, 3427–3435 (1990). http: //link.aps.org/abstract/PRA/v41/p3427. 17. J. Courtial, “Self-imaging beams and the Guoy [sic] effect,” Opt. Commun. 151, 1–4 (1998). 18. R. Zambrini, L. C. Thomson, S. M. Barnett, and M. Padgett, “Momentum paradox in a vortex core,” J. Modern Opt. 52(8), 1135–1144 (2005). http://www.informaworld.com/smpp/content∼content= a713736931. 19. C. Alonzo, P. J. Rodrigo and J. Glückstad, “Helico-conical optical beams: a product of helical and conical phase fronts,” Opt. Express. 13, 1749–1760 (2005). http://www.opticsexpress.org/abstract. cfm?URI=oe-13-5-1749.


Introduction
Optical vortices are phase singularities which are formed when light waves or beams are superposed [1]. In three dimensions these points of complete destructive interference trace out connected structures such as loops and knots [2,3,4]. The study of optical vortices is closely related to the field of orbital angular momentum (OAM) of light. Light beams with an azimuthal phase structure exp(imφ ), where m is an integer number, carry OAM of mh per photon [5], which arises directly from the azimuthal component of their Poynting vector [6]. A common example for such light beams are the Laguerre-Gaussian (LG) modes, which can be produced in the laboratory using spiral phase plates [7] or computer generated holograms [8]. Optical components which generate LG beams imprint a 2πm step in the phase of the electromagnetic field. However, it is also possible to design the step so that the phase jump is not an integer multiple of 2π. This gives rise to the phenomenon of non-integer, or fractional, OAM. Within the scope of our study irrational numbers are of no particular importance and we often use 'fractional', where 'non-integer' would be more complete. The vortex structure of light emerging from a non-integer phase step was first studied by Berry [9]. It is characterised by a chain of alternating vortices which forms in a dark line associated with the edge dislocation of the phase [10]. If one looks at the vortex structure in three dimensions one can see that every pair of alternating vortices in this chain is part of a nodal line in form of a 'hairpin', that is the two vortices converge at a common 'turning point' [11]. This chain of vortices affects the intensity and phase profiles of light emerging from a fractional phase step and makes the light unstable on propagation.
A different way to produce light with fractional OAM is as a generic superposition of light modes with different values of m. Using the correspondence between optics and quantum theory we have found a representation of light with fractional OAM as a quantum state. This state can be decomposed into a basis of integer OAM states [12]. The decomposition only determines the OAM index m, which in a superposition of LG beams leaves the index for the number of concentric rings unspecified. In this paper we make use of this flexibility to find a representation of a fractional OAM state in terms of LG beams with a minimal number of different Gouy phases to increase propagational stability. We produce these beams using a spatial light modulator (SLM) and study their propagation and vortex structure. Light beams constructed in this way are an excellent realisation of non-integer OAM states and they are more stable on propagation than light emerging from fractional phase steps. Spiral phase plates and light beams carrying non-integer OAM have also been studied in the context of two photon entanglement [13,14]. Experiments in this area could benefit greatly from an enhanced propagation distance for the light beams with fractional OAM.

Construction of light beams with non-integer OAM
We denote the state representing light emerging from a fractional phase step with |M(α) , where M = m + μ and m is the integer part and μ lies between 0 and 1. Unlike integer OAM states, which are fully characterised by the OAM index m, these fractional OAM states also depend on the orientation of the edge dislocation α [13,15]. These states can be decomposed into the basis of integer OAM states according to where the coefficients c m [M(α)] are given by [12] The angle θ 0 is an arbitrary starting point which defines the interval θ 0 ≤ φ < θ 0 + 2π for the azimuthal angle φ [16]. The orientation of the edge dislocation α is measured from θ 0 , so that α lies always between 0 and 2π. Although arbitrary, θ 0 defines one basis set for the angle and it is important to keep in mind that non-integer OAM states in general are different for different choices of θ 0 . In the remainder of this paper we set θ 0 to −π. We should stress that in the construction of the fractional OAM state (1) we have taken great care to avoid a multi-valued behaviour as known from complex analysis and easily associated with a fractional phase [12]. In the following we construct a physical representation of Eq. (1). The LHS is represented by the light field Ψ M(α) emerging from a fractional phase step with the phase discontinuity oriented at α. The RHS is a weighted superposition of integer OAM states. As we want to use the Gouy phase to control the vortex structure we use the complete set of LG modes u m p as implementation of integer OAM states. In cylindrical coordinates the field amplitude of the LG mode is given by: where w(z) is the Gaussian spot size k denotes the wavenumber, z R the Rayleigh range, w 0 the beam waist and L |m| p are the associated Laguerre polynomials. The normalisation constants of the LG modes are given by: Of particular interest for our purposes is the Gouy phase, exp[−i(2p + |m| + 1) tan −1 (z/z R )], which describes the phase change as the beam moves through the beam waist, situated at z = 0.
On propagation from the beam waist to the far field the term tan −1 (z/z R ) changes from 0 to π/2. The decomposition of the fractional state in Eq. (1) into integer OAM states uses odd and even values of m. As p has to be an integer it is not possible to have the same Gouy phase for all modes in the superposition, but we can limit the number of different Gouy phases to two, one for even values of m and one for odd values. This can be achieved by choosing appropriate p indices for each LG mode of the superposition. A similar technique has been used earlier to achieve a self-imaging effect [17], where by using superpositions of only odd values of m or only even values of m, the same Gouy phase could be obtained for all the constituent modes. We make use of this idea later in section 6 when we discuss light beams with fractional OAM which are structurally stable on propagation.
The decomposition in Eq. (1) includes an infinite number of integer OAM states, but the optical realisation as a superposition of laboratory beams is restricted by the apertures in the experiment. For a given circular aperture of radius R we require kR |m| for propagating modes [18]. We thus denote the finite number of contributing modes in the superposition by n modes . For every value of M our distribution of the |c m [M(α)]| 2 is peaked around the nearest integer to M (see Figure 1). For μ = 1/2 the modulus square of the coefficients either side of M is equally high. The finite superposition of LG modes is thus centred around the nearest integer to M: The smallest value of m in the superposition is determined by m min = Rnd[M − n modes /2], whereas the largest value is given by m max = m min + n modes − 1. Here Rnd denotes rounding up to the nearest integer greater than M for μ ≥ 1/2 and rounding down to the nearest integer smaller than M for μ < 1/2. It is important to note that 0 ≤ μ < 1, that is μ is not negative. A negative value of M, say M = −6.5 has an integer part of m = −7 and a fractional part of μ = +0.5. By setting the mode index p for each m in the superposition Eq. (5) to we can ensure that every LG mode has one of two different Gouy phases. Here Floor(q) gives the nearest integer smaller than, or equal to, q. A graph of the values for p and the resulting sum 2p + |m| + 1 can be found in Figure 1. The distribution of the p indices depends on the total number of modes in the superposition, and therefore it is not possible to add more modes to a given superposition without having to adjust all Gouy phases. This is because for every mode the sum 2p + |m| + 1 has to be equal to |m max | + 1 or |m max |. Although a finite set of modes in the superposition Eq. (5) does in general give a non-integer mean value for the OAM, it is not an exact representation of the state from Eq. (1). However, the distribution of the coefficients |c m | 2 in Figure 1 shows that modes with an OAM index m very different from M contribute only little to the superposition. The approximation can in principle be made arbitrarily good by adding more and more modes. However, any experimental realisation will be limited by the difficulties in reproducing the superposition. More modes lead to higher p indices and therefore to a larger intensity pattern with more phase discontinuities associated with the ring structure of the LG beams. Fitting a larger structure on the finite array of programmable pixels of the SLM requires a reduction of the beam waist which in turn leads to a reduction in resolution of the phase and intensity structures on the SLM. In the Appendix A we discuss the effect of a finite number of modes in the superposition on the quantum theoretical properties of the non-integer OAM state.

Experimental generation of non-integer OAM light beams
Rather than using several optical elements to generate each LG modes separately we employ a single SLM programmed with a hologram that sets the phase and intensity structure for the superposition. A blazed grating is also included in the hologram to separate angularly the first diffraction order. The formula for the resulting phase distribution of the hologram in rectilinear coordinates Φ(x, y) holo is given by [4] Φ(x, y) holo = Φ(x, y) beam + Φ(x, Λ) grating mod 2π − π sinc 2 [(1 − I(x, y) beam )π] + π. (7) Here, Φ(x, y) beam is the phase profile of the superposition at the beam waist for z = 0 and Φ(x, Λ) grating is the phase profile of the blazed grating which depends on the period of the grating Λ. The two phase distributions are added modulo 2π and, after subtraction of π, multiplied by an intensity mask. In regions of low intensity the intensity mask reduces the effect of the blazed grating, which in turn leads to reduced intensity in the first diffraction order. The mapping between the phase depth and the desired intensity is not linear but rather given by the trigonometric sinc function. In Figure 2 we illustrate the algorithm of the hologram designed to generate a superposition according to Eq. (5) for 10 modes and M = 6.5. Figure 3 depicts the experimental arrangement. The hologram is implemented using a computer controllable spatial light modulator (SLM). The light beam is produced with a HeNe laser, which has been expanded with help of the lenses f 1 and f 2 to illuminate the SLM evenly (see Figure 3). The recording of the phase profiles requires a reference beam which is coupled out from the illumination over a beam splitter and which is then reflected on mirror M 1 . This reference beam can be blocked by an electronic shutter. Depending on whether the shutter is closed or open, either the reflected light from the SLM or its interference with the reference beam is imaged by lens f 3 onto a spatial filter which only transmits the first positive diffraction order. The light is then passed through a collimation lens f 4 onto the mirrors M 3 and M 4 , which are mounted on a movable track to select a propagation distance, which is then imaged by lens f 5 onto a CCD array to record the image. The process of opening and closing the shutter, moving the mirror mount to the correct position for the desired propagation distance and recording the image is automated by use of a computer running a Labview programme. The programme also generates the desired hologram. In the setup used to record the experimental results we were able to take cross sections of the intensity and phase of the beam over a propagation distance from about z = 0 to z = 4z R . By adjusting the beam waist, longer propagation distances are feasible, but at four Rayleigh ranges, the phase change induced by the Gouy phase is already close to its maximum value.

Phase and intensity profiles
The intensity and phase profiles of light emerging from a half-integer phase plate change drastically on propagation [11,12]. In contrast, our synthesised light beams show small changes in the phase and intensity profiles on propagation from the near field to the far field. To compare the phase and intensity profiles without dilation due to the increase in the Gaussian spot size we introduce scaled coordinates x/w(z) and y/w(z). In Figure 4 we compare numerical and experimental pictures of the propagation of the beam. The series of graphs shows the intensity and phase profiles at the waist plane and after propagation of 2 and 4 Rayleigh ranges (z = 2z R and z = 4z R ) for M = 6.5 and a superposition of 10 LG modes according to Eq. (5). One can see that the intensity profile shows a rotation by π/2, which we attribute to the difference in the Gouy phase between the waist plane and the far field. This difference, however, affects only the modes in the superposition for which the sum 2p + |m| + 1 = 17. The other modes satisfy the self-imaging condition as the phase difference between the waist plane and the far field is an integer multiple of 2π. The self-imaging condition is specific for this particular superposition, while a rotation of the intensity structure on propagation is characteristic for our generic superposition if a sufficiently large number of modes is included. A rotation of the intensity structure can not always be associated with the Gouy phase [19]. It is interesting to note that the 'cap' of the open ring structure in the intensity profile for z = 0 corresponds to the orientation of the fractional phase step α (in all examples we have set θ 0 = −π and α = 0). Taking more modes in the superposition into account leads generally to more rings. This is because a higher total number of modes leads to a higher p index for those modes which contribute most to the superposition according to Eq. (6). The phase profiles depend, to a greater extent, on the number of modes included in the superposition. The initial vortex structure at the beam waist is dominated by the two modes which contribute most to the superposition (for our example of M = 6.5 the LG modes with m = 6 and m = 7). As long as there are more than two modes in the superposition, one vortex of charge ±1 will form along the orientation of α. Depending on the number of modes in the superposition the high charge vortices on axis split up into simple 2π vortices distributed around  the optical axis. The picture is a bit different if the superposition contains the mode with m = 0 as this mode adds on-axis intensity and leads always to a split of the high charge vortices. This effect is more pronounced if M is small because then the mode with m = 0 contributes more to the superposition. Figure 5 shows the intensity and phase profiles for a superposition of 6, 10 and 20 modes for M = 6.5. On propagation the vortex structure changes and aligns itself roughly along a line which is rotated by π/2 from the orientation of α. This will become clearer in the figures in the next section, where we compare theoretical and experimental pictures of the vortex structure on propagation.

Vortex structure
Phase profiles at different propagation distances give only a limited picture of the spatial form of the vortices. In Figure 6 we show a three dimensional plot of the nodal lines traced by the vortices. The graph compares the numerical calculation with the experimental results. Both plots are produced by searching the phase profiles for vortex positions at different propagation distances. The vortices are traced from the waist plane at z = 0 to a propagation distance of z = 2.5z R . One can see that the experimental result is in very good agreement with the numerical simulation and that the topological aspects of the vortex structure are correctly reproduced. The graphs show quite clearly the line of vortices which forms on propagation roughly 90 degrees rotated from the orientation of the cap. Another feature is the presence of 'hairpin' curves, a pair of vortices which are joined at a turning point. Shorter hairpins appear in the dark rings and extend only for a small fraction of a Rayleigh range, but one hairpin connects the vortex in the direction of α with a vortex from the ring structure and extends to a propagation distance of over one Rayleigh range. This vortex pair is important for the formation of a new vortex on changing M through the half-integer value. Interestingly, a new hairpin is also formed at the same position, z ≈ 2z R . The good agreement between the numerical and experimental results in the vortex structure is an indication that the desired and generated beam conform to a high degree as the vortex structure is very susceptible to perturbations [11]. The numerical calculations predict for example the existence of either a vortex of charge 2 or two vortices with charge 1 in the centre. The experimental results show that these two vortices start to spiral around each other on propagation, which indicates a small perturbation of the charge 1 vortices.

Structurally stable light beams with fractional OAM
Using a superposition of two fractional OAM states with a π difference in the orientation α, it is possible to create a beam which, when decomposed into OAM states, has either only odd or only even m. By adjusting the p values of the superposition modes, a uniform Gouy phase can be obtained. The corresponding light beam, apart from dilation, propagates unchanged and is hence structurally stable. We construct such a beam by setting For an even integer part of M the state |M + (α) gives a superposition of even modes, while for an odd integer part of M the decomposition contains only odd modes. The normalisation of these states depend on μ as only for μ = 1/2 fractional states are orthogonal if the difference in orientation is π. Because of the interference between the two consituent beams, the total OAM mean value is different from the single states |M(α) [10], and is given bȳ We can prepare a stable approximation to the superposition (8) by limiting the number of contributing modes, as described in section 2, and choosing the values of p to produce a unique Gouy phase. The resulting beam should be stable on propagation. For our example of M = 6.5 and 10 modes in the representation of |M(α) and |M(α + π) only the 5 even modes for m = 2, 4, 6, 8, 10 are left. Including 10 even modes in the representation of |M(α) would lead to a much higher number of rings. As with Figure 4 we adopt a radial coordinate such that divergence of the beam is surpressed. One can see clearly in Figure 7 that the scaled profiles remains unchanged on propagation. The phase profile shows roughly an m-fold rotational symmetry. This is because the superposition is now dominated by the mode with m equal to the integer part of M. The positions of the phase singularities remain unchanged as predicted.

Conclusion
It is well-established that beams with vortices can carry OAM and in particular that the LG modes carry an OAM of mh per photon, where m is the charge of the on-axis vortex [5]. Vortices of fractional strength can be incorporated only by also introducing a discontinuity in the beam. This follows directly and simply from the requirement that the field amplitude must be everywhere single valued. It is the discontinuity that leads to the break-up of the central vortex structure and to the consequent change in the beam profile on propagation [9,10].
We have shown that stable fractional OAM beams can be generated by using the radial degree of freedom. The changes upon propagation exhibited by conventional fractional OAM beams can be attributed to interference between the contributing LG modes with different Gouy phases. These differences can be suppressed, however, by selecting LG modes with appropriate radial indices, p, for each contributing value of m [17]. Our work has demonstrated the existence of simple superpositions of LG modes resulting in fractional OAM, but which undergo only a π/2 rotation on propagation from the beam waist to the far field, or even maintain their structural stability completely.
The structural stability of our beams is attractive for the generation or measurement of fractional OAM states in both the classical and quantum regimes. One potential application is in the measurement and manipulation of fractional OAM for entangled photon pairs [13].