Abstract
Control of orbital angular momentum (OAM) in optical fields has seen tremendous growth of late, with a myriad of tools existing for their creation and detection. What has been lacking is the ability to arbitrarily modify the OAM spectrum of a superposition in amplitude and phase, especially if a priori knowledge of the initial OAM spectrum is absent. Motivated by a quasi-mapping that exists between the position and OAM of Laguerre-Gaussian modes, we propose an approach for a single-step modulation of a field’s OAM spectrum. We outline the concept and implement it through the use of binary ring apertures encoded on spatial light modulators. We show that complete control of the OAM spectrum is achievable in a single step, fostering applications in classical and quantum information processing that utilise the OAM basis.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
It has been known for a long time that light can carry angular momentum in the form of spin (polarization) and orbital (OAM) components, but it wasn’t until approximately 25 years ago that OAM carrying fields could be easily created in the laboratory [1]. Since then, many approaches have been used to create structured light fields that contain scalar or vector superpositions of OAM [2,3], driven by the many applications such fields have cultivated [4,5]. In a similar vein, the detection of OAM has likewise received much attention with deterministic approaches based on conformal mapping solutions in so-called mode sorters [6,7]. What has remained somewhat lacking is the ability to arbitrarily modify some initial (perhaps unknown) OAM spectrum in amplitude and phase. So far, only limited tools exist based primarily on the principle of exploiting some property such as parity (odd/even modes) and symmetry (symmetric/antisymmetric modes) and routing the resulting OAM into paths [8–11]. But such approaches have not yet been able to shape the OAM spectrum because of the complexity and the limited versatility of the techniques. Further, once separated, each path would have to be appropriately modified and then recombined again later through an inverse system, thus requiring a multi-element and multi-step process. Often such approaches employ interferometers, further reducing the up-take due to stability issues. An approach to modify the OAM spectrum directly has been to control the range of the field’s azimuthal angle (the conjugate Fourier variable of OAM) [12], however this mostly only affects the dispersion of the OAM spectrum. Similarly, it was shown that phase dislocations, discontinuities and delays affect the OAM spectrum in distinct ways which has applications in imaging [13], but a method of tailoring the OAM spectrum in this way has not been demonstrated.
Altogether, arbitrary shaping of the OAM spectrum of light remains an open challenge. Yet, many applications require this, for example, the quantum Fourier transform as well as crucial steps in quantum computational algorithms such as that of Shor [14] and Grover [15] and general quantum information processing simulated with classical light [16–18]. In general, the task is to take an arbitrary initial OAM spectrum (specified by the coefficients $c_\ell$) and modify it according to $\sum c_\ell {\left |{\ell }\right \rangle } \rightarrow \sum \tilde {c}_\ell {\left |{\ell }\right \rangle }$. For instance, a key step in the Shor algorithm is to execute the transformation $c_\ell \rightarrow c_\ell \exp (2 \pi i a^{\ell }/N)$ for each $ {\left |{\ell }\right \rangle }$ in the initial OAM spectrum where $N$ is the number one wishes to find the factors of and $a<N$. Presently, there is no protocol for modifying an arbitrary initial OAM spectrum in such a manner and in one step. One should note the distinction between generating an engineered/tailored OAM spectrum (as in [19]) versus shaping an OAM spectrum that was previously generated.
Here, we outline and demonstrate an approach to shaping the OAM spectrum in a single step. We circumvent the difficulties of manipulating the spectrum directly in OAM space by recognizing a quasi-mapping between the OAM and position spaces for Laguerre-Gaussian (LG) and Perfect Vortex (PV) beams, where each $\ell$ maps to a radial position $r_\ell$. Such a mapping enables the manipulation of the OAM content through manipulation of the light field at specific spatial positions, facilitating easy OAM spectrum modulating in amplitude and phase using conventional devices such as spatial light modulators (SLMs). We then implement our scheme experimentally with PV beams and show arbitrary conversion of an initial OAM spectrum to some desired spectrum. In what follows, we give theoretical and experimental credibility to such a scheme.
2. Concept
We begin by phrasing the problem we wish to solve. Suppose we are given an arbitrary OAM-containing light field $U(r,\phi )$, where $(r,\phi )$ are transverse cylindrical coordinates. We can expand this field as,
Our proposed solution to the problem as stated above is a general procedure which is outlined visually in Fig. 1. Firstly, each OAM mode $ {\left |{\ell }\right \rangle }$ is associated to a unique spatial (radial) position $r_\ell$ via a quasi-mapping. In general, this mapping will be unique for the chosen set of vortex modes: $A_\ell (r)$. This allows $f(\ell )$ to be replaced by $f(r_\ell )$. Since only a spatially varying function is needed, many optical devices could implement this transformation, such as spatial light modulators. Secondly, annular ring apertures centred at $r_\ell$ and of a certain width $\rho _\ell - \rho _{\ell -1}$ (to be defined later) are assigned to each OAM mode one wishes to shape. Lastly, one specifies a transmission function $t_\ell$ for each ring aperture which will be directly related to the “amount” by which the corresponding OAM spectrum coefficient is modified ($t_\ell \sim f(\ell$)). The collection of ring apertures can then be applied to the field as a whole so as to shape the OAM spectrum in a single shot. Note that this scheme only requires the field’s inherent scale parameter as input. In the following, we will describe the specifics of each step of the process in detail.
2.1 Position-OAM mapping
As stated earlier, directly manipulating the OAM of light in an arbitrary way is difficult, particularly if the OAM spectrum is unknown. Instead, we exploit a mapping from OAM to position space $ {\left |{\ell }\right \rangle } \leftrightarrow r_\ell$, as depicted by the solid line in step 1 of Fig. 1. By way of example, consider the Laguerre-Gaussian modes $\textrm {LG}_p^{\ell }(r)$ which, at the waist plane $z=0$, are given by
Analogous quasi-mappings can be found for other vortex mode sets (numerically if not analytically). Since the vortex ring of light is typically the global maximum of the field amplitude, for many vortex modes the position-OAM mapping can be found from
2.2 Construction of annular ring apertures
The next step in the procedure is to construct the annular ring apertures centred at $r_\ell$, the boundaries of which are depicted by the horizontal lines in steps 2 and 3 of Fig. 1. It seems reasonable to assemble $ {\left |{\ell }\right \rangle }$’s annular ring about $r_\ell$ in such a way that it contains as much of the associated mode as possible without overlapping with adjacent rings. One has a certain degree of freedom here in choosing how to achieve this. One possibility is to choose the boundaries of the annular ring to be where the given mode $ {\left |{\ell }\right \rangle }$ first intersects with the adjacent modes: $ {\left |{\ell \pm 1}\right \rangle }$. In the case of LG modes, for example, the point of intersection of $ {\left |{\ell }\right \rangle }$ with $ {\left |{\ell +1}\right \rangle }$ can be calculated analytically (with all modes normalised to unit amplitude) and is given by,
2.3 Specifying the ring aperture transmission function
Once Eqs. (4) and (6) are determined for the specific vortex modes chosen, all that is left is to specify the transmission function $t_\ell$ for each ring aperture which is directly related to the desired transformation function $f(\ell )$. Physically speaking, we modify the amount and phase of light within the annular rings which, owing to the position-OAM mapping, modifies the amplitude and phase of the corresponding OAM modes. This is shown in steps 2 and 3 in Fig. 1 where the gray level within each ring aperture corresponds to $|t_\ell |$ and the colour corresponds to $\arg (t_\ell )$. For example, if one wishes to lower the contribution of $ {\left |{\ell }\right \rangle }$ to the OAM spectrum relative to the others, we encode $|t_\ell | < 1$. If one wishes to apply an inter-modal phase shift to $ {\left |{\ell }\right \rangle }$, we encode $\arg (t_\ell ) \in [-\pi ,\pi ]$.
By encoding the ring apertures on phase-only SLMs, the implementation of this step is straightforward: the amount of light within an annular ring can be modified by encoding a diffraction grating of variable height, thus displacing a controllable amount of light away from this region. Similarly, the phase of light within an annular ring can be modified by encoding a constant phase relative to some reference mode. In this implementation, light can only be removed from annular regions and so this imposes the restriction $|t_\ell | \leq 1$. In addition, due to the discrete nature of the SLM pixels, one should ensure that the pixelated approximation of an annular ring aperture is good or else performance will suffer.
Although it’s clear that the transmission function of the ring apertures and the transformation function are related, the main factor which determines their exact relationship is the OAM mode overlap. In general, OAM modes in the field will have some spatial overlap and as a consequence, this will inhibit the ability to modify each spectrum coefficient independently. The higher the spatial overlap, the more coupled the OAM modes become and the more difficult it is to determine the relationship between $t_\ell$ and $f(\ell )$ in our scheme. Conversely, if the OAM modes in the field are sufficiently spatially distinct, then $t_\ell = f(\ell )$ and one can directly encode the desired transformation (which is the ideal scenario).
3. Simulated OAM spectral shaping
We can test the proposed OAM shaping scheme by simulating the generation, modification and detection of superpositions of OAM modes. An example of such a simulation for LG modes is shown in Figs. 2(a)–2(f). Here, the following uniform ($c_\ell =c=1/\sqrt {11}$) OAM spectrum was generated
As mentioned earlier, the spatial overlap between adjacent LG-OAM modes (as can be seen in Fig. 2(a)) is the culprit of this unpredictability. In fact, in the limit that $\ell \rightarrow \infty$, adjacent LG-OAM modes overlap one another completely. Hence, modifying the amplitude and phase in one annular region affects the amplitude and phase of other OAM modes in a complex way. This means that unless some explicit relationship can be found, one may have to resort to trial and error to find appropriate transmission functions that will enable the intended shaping. For our example, we found by trial and error that if we wanted $f(\ell ) \propto 2^{\ell }$ for modes up to $\ell = 5$, we would have to encode $t_\ell \propto 2^{\ell }$ for ring apertures up to $r_{20}$. In other words, it is required to adjust the light in annular regions extending past $\ell = 5$ even though such modes were not present in the original spectrum.
We expect that the OAM shaping for other vortex modes (such as hypergeometric-Gauss and Bessel-Gauss) will exhibit similar behaviour to LGs since the OAM modes are distributed in a similar way. However, we can ask if there is a set of vortex modes that are sufficiently spatially separated to minimise modal overlap and thus allow a more predictable implementation.
In this regard, we turn to the so-called Perfect Vortex (PV) modes which are topical in the optical community owing to an annular field structure whose width is OAM independent [21]. We’ll see that if we utilise PVs as the OAM carriers, we can engineer an OAM-containing field in such a way that the OAM modes are sufficiently non-overlapping to enable $t_\ell = f(\ell )$. Further, there is the freedom to place OAM modes in any desired order and so one is not restricted to the inherent hierarchy in “regular” vortex modes of having larger OAM values further from the optical axis. As an added benefit, one can now differentiate between $ {\left |{\ell }\right \rangle }$ and $ {\left |{-\ell }\right \rangle }$ and shape these two modes independently.
PVs are not particularly exotic; in fact, they are known to be the Fourier transform of the well-studied Bessel-Gauss (BG) beams,
where $k_r$ is the radial wavenumber of light, $w_0$ is the Gaussian waist radius and $J_\ell (\cdot )$ is the Bessel function of order $\ell$. As such, the electric field of PV beams can be derived by Fourier transforming Eq. (9) and is described by [22],The basic idea for engineering non-overlapping OAM modes with PVs is to assign a single OAM mode to a single PV and to separate each PV from one another. Since PVs are generated from BGs, $k,w_0$ and $f$ will typically be global constants for all modes. This is because it is difficult to engineer a superposition of BGs where each mode has a different wavelength, width or Fourier lens. Hence, by Eqs. (11) and (12), the PV ring thickness is fixed and the ring radius is controlled by $k_r$. To force the PVs to be non-overlapping, the minimal-overlap condition can be phrased as follows: the difference in ring radii between adjacent rings should not be less than twice the ring thickness, or equivalently $\Delta R \geq 2 T$. Substituting Eqs. (11) and (12), this condition is reduced to,
Any choice of ring spacing satisfying this inequality will enable one to modify the OAM modes independently. In what follows we choose $\Delta k_r = 6/w_0$ to make the separation of rings more apparent. Deriving a closed form expression for the intersection of $ {\left |{\ell }\right \rangle }$ with $ {\left |{\ell +1}\right \rangle }$ for PVs (analogous to Eq. (5)) is elusive. However, when $R\gg T$ (or equivalently when $k_r w_0 \gg 1$), the intersection point is approximately the midpoint of the two ring radii.We now showcase (through simulation) the OAM spectrum shaping with PVs, which is shown in Figs. 2(g)–2(l). Similar to the flat LG-OAM spectrum shown earlier, we consider the following superposition,
where $R_\ell = R_{0} + \ell \Delta R$ is the OAM-position mapping with $R_0$ being the ring radius of some reference mode ($ {\left |{0}\right \rangle }$ in our case) and $\Delta R$ satisfying Eq. (13). This field amplitude is shown in Fig. 2(h) where each PV in the superposition is normalised to unit amplitude. The vertical dashed lines in Fig. 2(g) denote the boundaries between PV-OAM modes. As we see, the approximation that the intersection point is the midpoint of the two ring radii is good enough in this case. The ring apertures are encoded with the same transmission functions as to those used earlier for LG modes (only that now $ {\left |{\ell }\right \rangle }$ and $ {\left |{-\ell }\right \rangle }$ are distinguishable).We observe from the shaped coefficients that the transformation function $f(\ell )$ is precisely the encoded transmittance of the ring apertures $t_\ell$ (as desired). Furthermore, the experimenter now has the freedom to choose how to structure the OAM-position mapping as it is no longer imposed by nature. In our example, the OAM-position space mapping is linear and of the form
Note that both $R_0$ and $\Delta R$ can be manipulated as desired. It is also conceivable to manufacture a non-linear mapping.4. Experimental OAM spectral shaping
We now present the experimental setup used to verify the OAM shaping scheme with PVs. Similar to the simulations performed earlier, the experiment comprises generation, modification and detection phases. To this end, we built the experimental setup shown schematically in Fig. 3. The generation step of the setup creates a superposition of PV-OAM modes from an expanded and collimated He-Ne laser beam using a Spatial Light Modulator (Holoeye Pluto) employing complex amplitude modulation [23]. Using a 4f lens system, this field is then relayed to a second SLM which performs the modification and detection steps simultaneously (otherwise a third SLM would be required). This can be done since these steps occur at the same plane and the transmission functions of the holograms used for modification and detection are linear.
In order to measure the OAM spectrum coefficients $c_\ell$, we chose to utilise a standard modal decomposition setup which has been shown to be effective in detecting the OAM of PVs [24]. Specifically, to measure $|c_\ell |$ (we do not measure inter-modal phases here for brevity) the detection system needs to perform the following overlap integral,
4.1 Correction factors
Equation (17) suggests that the raw OAM spectrum we measure with this decomposition method requires post processing in the form of correction factors (this is true more generally when decomposing into a non-orthonormal basis such as the OAM eigenstates). Note that this is an implementation artefact of measuring the spectrum coefficients through an optical inner product [24] and is separate from the spectrum shaping process. Without this correction step, the values for $|c_\ell |$ measured via this approach will differ from the actual values. We point it out here because this artefact of optical modal decomposition when projecting onto the OAM eigenstates is not fully appreciated.
Experimentally, $|\gamma _\ell |^{2}$ controls the magnitude of the on-axis intensity that we measure. For the case of PVs, it turns out that $|\gamma _\ell |^{2}$ is proportional to the energy of each PV annulus in the transverse plane. Since each PV mode in the superposition in Eq. (14) has the same thickness but a different radius and thus a different area, the transverse energy of each OAM carrier is different. This means that we would measure different on-axis intensities for different OAM modes even though they may have the same weight in the spectrum. Hence, before shaping can be performed, we should ensure to normalise the energy in each OAM mode so as to measure the true spectrum. One way to do this is to apply the factors $1/|\gamma _\ell |$ to the raw coefficients $|c_\ell |$. These factors can be determined in one of two ways: numerically or empirically.
To numerically correct the measured on-axis signal we compute $|\gamma _\ell |$ explicitly from the radial integral in Eq. (17). These values correspond to the dashed line Fig. 4(b). This approach to correcting the spectrum will work assuming the generation and detection system is already well aligned and calibrated.
Empirically correcting the raw OAM spectrum involves sending a uniform superposition through the detection system. The correction factors are then precisely the inverse of each measured spectrum coefficient. This is also an effective way to correct for small systematic errors in the generation and/or detection system, such as a non-uniform beam from the laser, errors in determining the on-axis camera pixel coordinates or misalignment of the detection hologram. Both types of corrections are shown in Fig. 4(b). We see that the numerical and empirical corrections are in good agreement, indicating that our generation and detection systems are well aligned and calibrated. In our experiment, the source of systematic error was the fact that the laser beam at SLM$_1$ had some residual Gaussian intensity variation, resulting in modes near the optical axis having slightly more power than they should.
4.2 Results
Once the setup is calibrated and the correction factors found, more exotic OAM shaping can be pursued. We show four such examples in Fig. 5: where an initial uniform spectrum (A) is transformed to a triangular-shaped spectrum and a U-shaped spectrum and where an initial sinusoidal spectrum (B) is transformed to a flat spectrum and a U-shaped spectrum, all in a single step. Here, the initial spectrum is known and $|t_\ell |$ was determined by dividing the desired coefficients by the initial coefficients. Figure 6, on the other hand, shows three examples of shaping where the initial spectrum is not known. Here, a random OAM spectrum was generated and three different OAM filters were applied: a low-pass, high-pass and prime-pass filer.
In principle, there are no limitations to how the shaping can be done, one need only adjust the transmission function of the ring apertures accordingly, remembering that the amplitude of spectrum coefficients can only be made smaller in our implementation.
4.3 Inter-modal phase shaping
Although we do not explicitly show inter-modal phase shaping in our previous examples, the procedure is very similar to how the amplitude of $c_\ell$ is manipulated and is straightforward to implement if utilising SLMs. The difference is that $t_\ell$ becomes complex and so one also encodes a constant phase value $\arg \,(t_\ell )$ within the ring aperture which is defined relative to an arbitrary reference mode in the spectrum.
To showcase this, we present qualitative evidence (instead of doing the full modal decomposition) that our scheme has shaped the inter-modal phases. In this regard, we experimentally generated the uniform LG and PV superpositions as in Fig. 2 (having no inter-modal phases) and encoded $\arg (t_\ell ) = (-1)^{|\ell |-1} \pi /2$ and $|t_\ell | = 1$ within each ring aperture. Since the amplitude of the field is unaffected at the modification plane, we can observe the phase changes at the Fourier plane. The results are shown in Fig. 7. The insets correspond to simulated intensities at the Fourier plane having the OAM spectrum as shown alongside the figure. We argue that the high fidelity between the theoretical and experimental images is qualitative evidence for the shaping of inter-modal phases.
Similar to the shaping of $|c_\ell |$ for LG modes, the spatial overlap between LG-OAM modes inhibits the ability to modify the inter-modal phases independently. Thus, it is difficult to determine the exact relationship between $\arg (t_\ell )$ (what we encode) and $\arg (f(\ell ))$ (the transformation we want). In addition, shaping the inter-modal phases also has an effect on the amplitude of $c_\ell$ for the case of LGs; phase discontinuities at the ring aperture boundaries invariably affect the amplitude of the beam. For PVs, on the other hand, the engineered spatial independence of each OAM mode facilitates $\arg (f(\ell )) = \arg (t_\ell )$ and the amplitude of the OAM spectrum is unchanged.
5. Perspectives
Overall, our OAM shaping protocol from start to finish involves 3 steps: mapping OAM states to specific spatial positions, constructing binary ring apertures for those regions and then modifying the amount and phase of light within those regions. Step 1 is specific to the particular vortex modes chosen and step 3 depends on how one wishes to shape the OAM spectrum; these two steps are not subject to much change. Of the three, step 2 is the most malleable. If the OAM modes are sufficiently spatially separated (as is the case with the engineered PV superposition) then choosing ring apertures in such a way as to enclose the entire mode is the most fruitful approach: there is a direct equivalence between the desired shaping (defined by $f(\ell$)) and what is encoded on the ring aperture (defined by $t_\ell$). With regards to “regular” vortex modes, perhaps future work could explore different ways of constructing the ring apertures than the approach we considered here. It may be that there is a more optimal choice which will enable the shaping of these vortex modes to be more predictable.
Since it is significantly more difficult to add/inject light into specific annular regions, this meant that we only considered the case where the amplitude of the OAM coefficients is lessened. One way to increase the relative contribution of a particular OAM mode is to decrease the relative contribution of all other modes. This is equivalent but has the potential disadvantage of significantly reducing the overall power in the beam. A solution to this would be to utilise a gain medium in conjunction with a SLM. In this way, individual modes in the OAM spectrum could be both raised and lowered in amplitude.
6. Conclusion
To summarise, we proposed and demonstrated a scheme for single-step modulation of a light field’s OAM spectrum which utilises a quasi-mapping between the position and OAM spaces. This grants one the ability to shape the OAM spectrum using conventional devices (such as SLMs). Using our approach, we have shown that complete control of the amplitude and phase of the OAM spectrum is achievable in a single step. While most vortex mode sets have an inevitable spatial overlap, PV modes can be engineered with minimal spatial overlap, thus facilitating greater control over the OAM spectrum. Since OAM is being utilised more and more for applications in classical and quantum information theory, we conjecture that the ability to arbitrarily modify the OAM spectrum could drive new techniques that would enhance some of these technologies.
Funding
Council for Scientific and Industrial Research, South Africa.
References
1. L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. P. Woerdman, “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A 45(11), 8185–8189 (1992). [CrossRef]
2. A. Forbes, A. Dudley, and M. McLaren, “Creation and detection of optical modes with spatial light modulators,” Adv. Opt. Photonics 8(2), 200–227 (2016). [CrossRef]
3. C. Rosales-Guzmán and A. Forbes, How to Shape Light with Spatial Light Modulators (SPIE, 2017).
4. H. Rubinsztein-Dunlop, A. Forbes, M. V. Berry, M. R. Dennis, D. L. Andrews, M. Mansuripur, C. Denz, C. Alpmann, P. Banzer, T. Bauer, E. Karimi, L. Marrucci, M. Padgett, M. Ritsch-Marte, N. M. Litchinitser, N. P. Bigelow, C. Rosales-Guzman, A. Belmonte, J. P. Torres, T. W. Neely, M. Baker, R. Gordon, A. B. Stilgoe, J. Romero, A. G. White, R. Fickler, A. E. Willner, G. Xie, B. McMorran, and A. M. Weiner, “Roadmap on structured light,” J. Opt. 19(1), 013001 (2017). [CrossRef]
5. M. J. Padgett, “Orbital angular momentum 25 years on,” Opt. Express 25(10), 11265–11274 (2017). [CrossRef]
6. G. C. G. Berkhout, M. P. J. Lavery, M. J. Padgett, and M. W. Beijersbergen, “Measuring orbital angular momentum superpositions of light by mode transformation,” Opt. Lett. 36(10), 1863–1865 (2011). [CrossRef]
7. G. Ruffato, M. Girardi, M. Massari, E. Mafakheri, B. Sephton, P. Capaldo, A. Forbes, and F. Romanato, “A compact diffractive sorter for high-resolution demultiplexing of orbital angular momentum beams,” Sci. Rep. 8(1), 10248 (2018). [CrossRef]
8. G. C. G. Berkhout, M. P. J. Lavery, J. Courtial, M. W. Beijersbergen, and M. J. Padgett, “Efficient sorting of orbital angular momentum states of light,” Phys. Rev. Lett. 105(15), 153601 (2010). [CrossRef]
9. M. P. J. Lavery, A. Dudley, A. Forbes, J. Courtial, and M. J. Padgett, “Robust interferometer for the routing of light beams carrying orbital angular momentum,” New J. Phys. 13(9), 093014 (2011). [CrossRef]
10. Y. Zhang, F. S. Roux, T. Konrad, M. Agnew, J. Leach, and A. Forbes, “Engineering two-photon high-dimensional states through quantum interference,” Sci. Adv. 2(2), e1501165 (2016). [CrossRef]
11. J. Liu and J. Wang, “Demonstration of reconfigurable joint orbital angular momentum mode and space switching,” Sci. Rep. 6(1), 37331 (2016). [CrossRef]
12. B. Jack, M. J. Padgett, and S. Franke-Arnold, “Angular diffraction,” New J. Phys. 10(10), 103013 (2008). [CrossRef]
13. L. Torner, J. P. Torres, and S. Carrasco, “Digital spiral imaging,” Opt. Express 13(3), 873–881 (2005). [CrossRef]
14. P. W. Shor, “Algorithms for quantum computation: Discrete logarithms and factoring,” in “Proceedings 35th Annual Symposium on Foundations of Computer Science,” (IEEE, 1994), pp. 124–134.
15. B. Perez-Garcia, R. I. Hernandez-Aranda, A. Forbes, and T. Konrad, “The first iteration of grover’s algorithm using classical light with orbital angular momentum,” J. Mod. Opt. 65(16), 1942–1948 (2018). [CrossRef]
16. T. Konrad and A. Forbes, “Quantum mechanics and classical light,” Contemp. Phys. 60(1), 1–22 (2019). [CrossRef]
17. A. Forbes, A. Aiello, and B. Ndagano, “Classically entangled light,” Prog. Opt. 64, 99–153 (2019). [CrossRef]
18. E. Toninelli, B. Ndagano, A. Vallés, B. Sephton, I. Nape, A. Ambrosio, F. Capasso, M. J. Padgett, and A. Forbes, “Concepts in quantum state tomography and classical implementation with intense light: a tutorial,” Adv. Opt. Photonics 11(1), 67–134 (2019). [CrossRef]
19. F. Li, T. Xu, W. Zhang, X. Qiu, X. Lu, and L. Chen, “Optical images rotation and reflection with engineered orbital angular momentum spectrum,” Appl. Phys. Lett. 113(16), 161109 (2018). [CrossRef]
20. C. Schulze, A. Dudley, D. Flamm, M. Duparré, and A. Forbes, “Measurement of the orbital angular momentum density of light by modal decomposition,” New J. Phys. 15(7), 073025 (2013). [CrossRef]
21. A. S. Ostrovsky, C. Rickenstorff-Parrao, and V. Arrizón, “Generation of the “perfect” optical vortex using a liquid-crystal spatial light modulator,” Opt. Lett. 38(4), 534–536 (2013). [CrossRef]
22. P. Vaity and L. Rusch, “Perfect vortex beam: Fourier transformation of a bessel beam,” Opt. Lett. 40(4), 597–600 (2015). [CrossRef]
23. V. Arrizon, U. Ruiz, R. Carrada, and L. A. González, “Pixelated phase computer holograms for the accurate encoding of scalar complex fields,” J. Opt. Soc. Am. A 24(11), 3500–3507 (2007). [CrossRef]
24. J. Pinnell, V. Rodríguez-Fajardo, and A. Forbes, “Quantitative orbital angular momentum measurement of perfect vortex beams,” Opt. Lett. 44(11), 2736–2739 (2019). [CrossRef]