## Abstract

We demonstrate a simple approach, using digital holograms, to perform a complete azimuthal decomposition of an optical field. Importantly, we use a set of basis functions that are not scale dependent so that unlike other methods, no knowledge of the initial field is required for the decomposition. We illustrate the power of the method by decomposing two examples: superpositions of Bessel beams and Hermite-Gaussian beams (off-axis vortex). From the measured decomposition we show reconstruction of the amplitude, phase and orbital angular momentum density of the field with a high degree of accuracy.

©2012 Optical Society of America

## 1. Introduction

Techniques to decompose light by use of Fourier optics have been known for a long time, and have been extensively reviewed to date [1]. Historically, these techniques have been applied to pattern recognition problems, and subsequently to the problem of studying the structure and propagation characteristics of laser beams [2–5]. Despite the appropriateness of the techniques, the experiments were nevertheless rather complex. Recently this subject has been revisited by employing computer generated holograms in a mode multiplexing scheme for the modal decomposition of laser beams from fibres [6–10], and for the real-time measurement of the beam quality factor of a laser beam [11]. While these techniques have significant merit, they require a prefabricated diffractive optical element. This implies that information on the modal basis to be used, and the scale parameters of this basis, are known. To date this has been achieved by first modelling the source under study.

In this paper we consider the problem of the azimuthal decomposition of an arbitrary laser source, without any knowledge of the mode structure, the mode phases, or the scale of the amplitude distribution. We make use of a basis comprising the angular harmonics, which are independent of spatial scale, and express the spatial distribution in terms of spatially dependent coefficients in this basis. The result, as we will show, is that the complete decomposition can be achieved without any scale information. We use this to infer directly from measurements of the intensity of the superposition field, its phase, and its orbital angular momentum (OAM) density distribution. In fact, it is clear that since the entire field is known, all physical quantities associated with the field can be inferred. We illustrate the concept by executing a full azimuthal decomposition of: (i) a superposition of two OAM carrying Bessel beams, with relative phase differences and (ii) an off-axis vortex mode. A comparison of our experimental measurements to the predicted theory shows excellent agreement.

## 2. Concept and theory

The core idea is to expand our unknown field, *u*(*r,ϕ*), into a basis that is not dependent on scale. Such a basis is the angular harmonics, exp(*ilφ*), that are orthogonal over the azimuthal plane. To illustrate that the azimuthal decomposition that we employ in this work is completely general, we consider, for example, the Laguerre-Gaussian (LG) basis as our OAM basis. An expansion of an arbitrary optical field in terms of this LG basis is given by

*c*′

*denotes the complex coefficients and*

_{p,l}*R*(

_{p,l}*r*) is the radial part of the LG mode, which only depends on

*r*. One can now combine, and sum over the part of the expression that contains a

*p*-index

*r*-dependent coefficients. The coefficients are given by

*θ*, relative to some reference, which we may take to be an external source or for convenience the first mode in the series,

_{l}*l*= 0 at a specific value of

*r*. To describe

*u*(

*r,ϕ*) completely, we are required to determine the phase shift between modes and the magnitude of the coefficients. The coefficients may be found by an inner product calculation and a suitable match filter, as shown in Eq. (4). To implement this calculation in a physically realizable manner, consider the signal at the origin of the Fourier plane after a lens (for wavenumber

*k*and focal length

*f*) after the field has been modulated by a transmission function given by

*t*(

_{l}*r,ϕ*)

*r*it should consist of an annular slit centered at

*r = R*and of negligible thickness, Δ

*R.*These conditions are satisfied if the transmission function is defined as

*l*in

*t*and

_{l}*U*refer to selecting the

_{l}*l*

^{th}mode in the expansion of Eq. (2). So that

*I*

_{l}(

*R*,0), at the origin in the Fourier plane returns exactly our desired magnitude of the coefficients, sinceSecondly, we note that if we interfere our selected

*l*mode with a reference wave and consider the inner product signal, then the unknown phase of the coefficient, Δ

*θ*, can be found. If |

_{l}*g*|exp(

*iα*) represents the complex amplitude of the reference wave of the origin then the intensity, at the origin, as a function of the phase delay is given by

*u*(

*r,ϕ*)|

^{2}, the phase of the field, arg[

*u*(

*r,ϕ*)], the Poynting vector [12],

The pertinent point here is that at no stage has any assumption been made on the scale of the basis functions. This is because we have elected to scan the field with an annular slit, thereby storing the scale information in the coefficient of a scale invariant basis. The scale of any OAM basis function is set by the radial coordinate. Since the radial coordinate becomes a part of the coefficient functions, the scale parameter is removed from the basis functions, resulting in the remaining azimuthal component becoming independent of scale. We will now show how to implement the technique with a phase-only spatial light modulator (SLM), thus allowing a full azimuthal decomposition of the field with digital holograms.

## 3. Experimental methodology

The experimental realization of the technique comprises two parts: the generation and then decomposition of the field, and is shown schematically in Fig. 1(a)
. The first SLM (denoted as SLM_{1}) was programmed to produce various fields. The two cases that we considered were (i) a superposition of two OAM carrying Bessel beams [14] and (ii) an off-axis vortex mode. A Gaussian beam from a HeNe laser [Fig. 1(b)] was expanded through a 5 × telescope and directed onto the liquid crystal display of SLM_{1} (HoloEye, PLUTO-VIS, with 1920 × 1080 pixels of pitch 8 μm and calibrated for a 2π phase shift at λ~633 nm) where the hologram used to generate the field of interest was programmed. For the case where a superposition of two OAM carrying Bessel beams were studied, two annular rings modulated in their azimuthal phase [see Fig. 1(c)] were encoded onto SLM_{1}. The amplitude transmission of 0 everywhere outside the annular ring and 1 inside the ring was programmed using complex amplitude modulation for amplitude only effects on a phase-only device. The hologram takes the form of a high frequency grating that oscillates between phase values of 0 and π, and it has been shown that this results in the required amplitude transmission of Eq. (6) [15]. We shall refer to this as the checkerboard pattern: see zoomed in section of Fig. 1(c). The result is non-diffracting petal-like modes [Fig. 1(d)] that rotate as they propagate, and have been studied in detail elsewhere [15,16].

These fields were then magnified with a 10 × objective and directed to the second SLM (SLM_{2}) for executing the azimuthal decomposition. This was accomplished by executing an inner product of the incoming field with the match filter set to exp(*ilϕ*), for various *l* values, and for particular radial (*r*) positions on the field, as given by Eq. (6) and illustrated in Fig. 1(e). The typical width of the ring in the experiments was Δ*R* = 80 μm. The width of the ring, as well as the azimuthal index (*l*) encoded within the ring, are limited by the resolution of the SLM. The field at the Fourier plane is shown in Fig. 1(f) – note that we only require the intensity at the origin of this plane.

In order to measure the phase delay between the modes, we switch off the checkerboard pattern in the centre of both SLMs. This results in a portion of the initial Gaussian beam passing through the entire optical system, without “seeing” the phase holograms. Because it follows the same path as the modes that we wish to study, we can use the central peak of the beam as our reference beam, |*g*|exp(*iα*). Figure 1(g) shows the initial generating step but with a constant phase in the centre of the hologram (black disk), the resulting mixing with the non-diffracting petals [Fig. 1(h)] and after passing through the second SLM [Fig. 1(i)] the resulting interference of the waves [Fig. 1(j)]. The constant phase of the reference could be changed across several values to provide information on the quadrant of the cosine function in Eq. (10). Note that switching between the intensity measurement [Figs. 1(c)–1(f)] and the phase measurement [Figs. 1(g)–1(j)], we require only a change in the holograms loaded onto SLMs 1 and 2. This switching between holograms can be done at 60 Hz, i.e., practically real-time acquisition of the necessary data.

## 4. Results and discussion

#### 4.1. Superposition of two OAM carrying Bessel beams

To test the accuracy of the method, we first created a superposition of two Bessel beams with opposite handedness, i.e., *m* = 3 and *n* = –3 and defined generally as

*θ*, in Eq. (13) was then varied between 0 and 2π, resulting in a rotation of the intensity pattern, which can easily be measured by observing the angular displacement of any one of the petals. This measured rotation is directly related to Δ

*θ*, thus allowing the measured phase shift to be compared to the programmed value. The programmed phase shift and resulting petal rotation are shown in Fig. 2(a) , and compared to the theoretical rotation. The results of this calibration test are shown in Fig. 2(b), where it is clear that the agreement between measured and actual is very good (slope difference of 0.1%). The error bars were determined through a range of intensity values in the finite region (few pixels) around the origin of the Fourier plane. These results confirm that the setup was working correctly without spurious effects.

Next, the programmed phase shift Δ*θ* was measured using the interference technique and analyzed using Eq. (10). The phase shift in the *n* = –3 beam was scanned through 0 to 2π and the interference monitored, as shown in Fig. 2(c). By removing the offset phase of the reference beam (corresponding to the maximum in the interference signal), the measured phase could be compared to the programmed phase. The results are also shown in Fig. 2(d), and again the agreement is good (slope difference of 3%), but with some uncertainty near a programmed phase shift of 0.60π and 1.60π; these discontinuities occur at maxima and minima in the interference intensity [Fig. 2(c)], and are due to the flatness of the change in the measured intensity for a change in phase value (three values were used to uniquely determine the phase). To illustrate further that the technique works, we constructed a superposition field with *m* = + 3 and *n* = –2, to deliberately break the symmetry of the field. The annular ring on SLM_{2}, of width 80 μm, was scanned through 28 radial positions (maximum number to scan our field) and the azimuthal phase varied across *l*∈[–10,10]. The only non-zero components (>0.5% of the modal power) were those corresponding to *l* = + 3 and *l* = –2, which contained ~95% of the modal power. With all the unknown terms in Eq. (4) measured experimentally, we were able to reconstruct the intensity and phase of the initial field, and compare it to the theoretical prediction (based on our own hologram function on SLM_{1}). The results for the intensity of the field are shown in Figs. 3(a)
and 3(d) while that of the phase of the field is given in Figs. 3(b) and 3(e). There is very good agreement between the theoretically predicted field and the field reconstructed from the azimuthal decomposition.

An interesting application of this technique is the determination of the complete OAM density of the field. It has been known for some time now that photons carry orbital angular momentum ([17] and references therein), yet it remains a challenge to measure the local and global OAM of a field in a manner that is quantitative. Previously we have attempted such measurements without the phase information, and found success for particular structures of optical fields [12,13]. With a full azimuthal decomposition of the field, including the phase delays, we can infer the OAM density of *any* field directly from our measurement. The theoretical and experimental results of this are shown in Figs. 3(c) and 3(f), respectively, and are in agreement.

#### 4.2. Off-axis vortex mode

Next we apply our decomposition technique to the example of a Gaussian beam with an off-axis vortex. The results are presented in Fig. 4 .

The same measurement procedure, as that discussed above, was implemented on the field, whose experimentally recorded intensity profile is given in Fig. 4(a). The reconstructed field is shown in Fig. 4(b) illustrating that this technique is successful in the reconstruction of optical fields which possess an on-axis intensity. The problem in this case is that there is no dark core in the centre of the field under investigation where we could pass a Gaussian reference beam through, for the implementation of the phase measurements. In this case we use a section of the field under investigation as our reference (assuming no singularities, a single pixel size region will have a uniform phase).

#### 4.3. Alignment sensitivity

Although one can represent the beam in terms of expansions with respect to any beam-axis, the coefficients in the expansions are sensitive to a shift of the beam-axis. Here, our concern is the fidelity of the reproduced beam rather than the fidelity of the expansion coefficients. Naturally we do not want the beam-axis to be too far off, because that would increase the higher-order coefficients and reduce the fidelity of the reconstructed beam, but small shifts in the beam axis does not pose a serious problem and this is illustrated in Fig. 5(a)
. The hologram on SLM_{2} was displaced with respects to the propagation axis of the field under investigation for step sizes of a single pixel (i.e. 8 µm) and the fluctuation in the on-axis intensity of the inner product, from which the unknown field is reconstructed, shows very little variation (the percentage error is roughly 0.5% per pixel shift). It is fairly easy to get the alignment of our system within such distances from the actual beam axis. If the displacement is much greater (10 times greater) then the alignment does pose a problem and the on-axis intensity of the inner product varies significantly, evident in Fig. 5(b).

For the experimental measurements presented here we aligned each SLM to a common propagation axis by initially centering each SLM on a Gaussian mode that propagated through the system. More general alignment techniques, as used by other decomposition methods, could also be employed [6,7].

## 5. Conclusion

We have outlined an improved method for the azimuthal decomposition of an arbitrary field that requires no scale information of the basis functions, and allows for the phases to be easily measured. It is significant to note that using only minor changes to the digital holograms, we can switch the setup between measuring azimuthal weights to azimuthal phase delays – no optical re-alignment is necessary, and no path length adjustment is required for the phase measurements. Moreover, even though our system is not automated and a scan of all measurements would take 4 hours, in an automated system the holograms can be displayed at fast refresh rates (60 Hz), so that the measurement per radial position takes about 1 s. Using superpositions of OAM carrying Bessel beams as an example, we have shown excellent reconstruction of both the intensity and phase of the field, as well as the OAM density everywhere in the field. We have also shown accurate reconstruction of an off-axis vortex beam, illustrating the versatility of the method. This technique will surely be of relevance to those working in OAM, both at the classical and quantum level, as well as to the emerging field of mode multiplexing for future telecommunications ([18], and references therein), where mode generation and decomposition are necessary tools.

## References and links

**1. **J. W. Goodman, *Introduction to Fourier Optics* (McGraw-Hill Publishing Company, 1968).

**2. **E. Tervonen, J. Turunen, and A. Friberg, “Transverse laser mode structure determination from spatial coherence measurements: experimental results,” Appl. Phys. B **49**(5), 409–414 (1989). [CrossRef]

**3. **A. Cutolo, T. Isernia, I. Izzo, R. Pierri, and L. Zeni, “Transverse mode analysis of a laser beam by near- and far-field intensity measurements,” Appl. Opt. **34**(34), 7974–7978 (1995). [CrossRef] [PubMed]

**4. **M. Santarsiero, F. Gori, R. Borghi, and G. Guattari, “Evaluation of the modal structure of light beams composed of incoherent mixtures of Hermite-Gaussian modes,” Appl. Opt. **38**(25), 5272–5281 (1999). [CrossRef] [PubMed]

**5. **X. Xue, H. Wei, and A. G. Kirk, “Intensity-based modal decomposition of optical beams in terms of Hermite-Gaussian functions,” J. Opt. Soc. Am. A **17**(6), 1086–1091 (2000). [CrossRef] [PubMed]

**6. **D. Flamm, O. A. Schmidt, C. Schulze, J. Borchardt, T. Kaiser, S. Schröter, and M. Duparré, “Measuring the spatial polarization distribution of multimode beams emerging from passive step-index large-mode-area fibers,” Opt. Lett. **35**(20), 3429–3431 (2010). [CrossRef] [PubMed]

**7. **T. Kaiser, D. Flamm, S. Schröter, and M. Duparré, “Complete modal decomposition for optical fibers using CGH-based correlation filters,” Opt. Express **17**(11), 9347–9356 (2009). [CrossRef] [PubMed]

**8. **J. W. Nicholson, A. D. Yablon, S. Ramachandran, and S. Ghalmi, “Spatially and spectrally resolved imaging of modal content in large-mode-area fibers,” Opt. Express **16**(10), 7233–7243 (2008). [CrossRef] [PubMed]

**9. **D. B. S. Soh, J. Nilsson, S. Baek, C. Codemard, Y. Jeong, and V. Philippov, “Modal power decomposition of beam intensity profiles into linearly polarized modes of multimode optical fibers,” J. Opt. Soc. Am. A **21**(7), 1241–1250 (2004). [CrossRef] [PubMed]

**10. **M. Paurisse, L. Lévèque, M. Hanna, F. Druon, and P. Georges, “Complete measurement of fiber modal content by wavefront analysis,” Opt. Express **20**(4), 4074–4084 (2012). [CrossRef] [PubMed]

**11. **O. A. Schmidt, C. Schulze, D. Flamm, R. Brüning, T. Kaiser, S. Schröter, and M. Duparré, “Real-time determination of laser beam quality by modal decomposition,” Opt. Express **19**(7), 6741–6748 (2011). [CrossRef] [PubMed]

**12. **I. A. Litvin, A. Dudley, and A. Forbes, “Poynting vector and orbital angular momentum density of superpositions of Bessel beams,” Opt. Express **19**(18), 16760–16771 (2011). [CrossRef] [PubMed]

**13. **A. Dudley, I. A. Litvin, and A. Forbes, “Quantitative measurement of the orbital angular momentum density of light,” Appl. Opt. **51**(7), 823–833 (2012). [CrossRef] [PubMed]

**14. **R. Vasilyeu, A. Dudley, N. Khilo, and A. Forbes, “Generating superpositions of higher-order Bessel beams,” Opt. Express **17**(26), 23389–23395 (2009). [CrossRef] [PubMed]

**15. **R. Rop, A. Dudley, C. Lopez-Mariscal, and A. Forbes, “Measuring the rotation rates of superpositions of higher-order Bessel beams,” J. Mod. Opt. **59**(3), 259–267 (2012). [CrossRef]

**16. **R. Rop, I. A. Litvin, and A. Forbes, “Generation and propagation dynamics of obstructed and unobstructed rotating orbital angular momentum-carrying Helicon beams,” J. Opt. **14**(3), 035702 (2012). [CrossRef]

**17. **A. M. Yao and M. J. Padgett, “Orbital angular momentum: origins, behaviour and applications,” Adv. Opt. Photon. **3**(2), 161–204 (2011). [CrossRef]

**18. **G. Li and X. Liu, “Focus Issue: Space multiplexed optical transmission,” Opt. Express **19**(17), 16574–16575 (2011). [CrossRef] [PubMed]