Abstract
Turing patterns in self-focussing nonlinear optical cavities pumped by beams carrying orbital angular momentum (OAM) $m$ are shown to rotate with an angular velocity $\omega = 2m/R^2$ on rings of radii $R$. We verify this prediction in 1D models on a ring and for 2D Laguerre-Gaussian and top-hat pumps with OAM. Full control over the angular velocity of the pattern in the range $- 2m/R^2 \le \omega \le 2m/R^2$ is obtained by using cylindrical vector beam pumps that consist of orthogonally polarized eigenmodes with equal and opposite OAM. Using Poincaré beams that consist of orthogonally polarized eigenmodes with different magnitudes of OAM, the resultant angular velocity is $\omega = (m_L + m_R)/R^2$, where $m_L, m_R$ are the OAMs of the eigenmodes, assuming good overlap between the eigenmodes. If there is no, or very little, overlap between the modes then concentric Turing pattern rings, each with angular velocity $\omega = 2m_{L,R}/R^2$ will result. This can lead to, for example, concentric, counter-rotating Turing patterns creating an optical peppermill-type structure. Full control over the speeds of multiple rings has potential applications in particle manipulation and stretching, atom trapping, and circular transport of cold atoms and BEC wavepackets.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Pattern formation is ubiquitous in nonlinear dynamical systems, the most famous examples being Turing patterns in reaction-diffusion systems [1] and Benjamin-Feir modulational instability in water waves [2]. In propagation through self-focusing nonlinear media with optical vortex beams, an azimuthal modulational instability leads to beam fragmentation as described in [3–6]. In optical cavities, instead, spontaneous spatial pattern formation results due to the Turing instability deriving from the interplay of nonlinearity, cavity detuning and a spatial coupling, such as diffraction or dispersion. Such systems are very well described by the Lugiato-Lefever equation for optical cavities with plane-plane mirrors [7]. In this paper we consider the effect of pumping optical cavities containing a self-focussing Kerr medium with beams carrying optical angular momentum (OAM) and show the formation of rotating Turing structures. We derive analytical expressions that fully describe these two-dimensional rotating Turing structures in single field (scalar) Kerr resonators and confirm our predictions numerically using pumps consisting of Laguerre-Gaussian modes or top-hat beams carrying OAM. In particular, we demonstrate that the number of intensity peaks in the patterns is determined by the power and OAM of the input pump, as well as the cavity detuning. We also demonstrate that the angular velocity $\omega$ of the patterns is fully determined by the OAM $m$ of the pump and the radius $R$ of the ring structure according to $\omega = 2 m/R^2$. Spatial structures rotating on a transverse ring including cavity solitons can be considered as slow light pulses with fully controllable speed and structure for use in optical quantum memories and delay lines. These studies complete early investigations that focused on optical parametric oscillators, semiconductor heterostructures and photorefractive materials, respectively [8–10].
Fully-structured light consisting of a vector superposition of two scalar OAM-carrying Laguerre-Gaussian eigenmodes with orthogonal circular polarizations [11–13], has attracted increasing attention for a number of applications [14–17]. The inclusion of a second field component in the light-matter interaction inside the cavity offers further degrees of control in the shape and polarization of the pump and the resultant nonlinear structures. In particular we show how the use of fully-structured light to pump the cavity allows us full control over the angular velocity of the Turing structures. Using numerical simulations we demonstrate how biasing cylindrical vector (CV) beam pumps - orthogonally polarized eigenmodes with equal and opposite OAM ($m_L = -m_R$) - allows us to produce patterns with angular velocity $-2 m/R^2 \le \omega \le 2 m/R^2$ and that Poincaré pumps - orthogonally polarized eigenmodes with different magnitudes of OAM - produce patterns with angular velocity $\omega = (m_L + m_R)/R^2$. Applications of these rotating structures to particle manipulation, optical beam shaping and photonic devices is discussed. Finally we give examples of fields with counter-rotating Turing patterns that we refer to as an optical peppermill, that may be of particular interest in trapping, manipulating and deforming biological specimens.
2. Lugiato-Lefever case
We start with the description of Kerr media in plane-plane optical cavities through the well-known Lugiato-Lefever equation (LLE) in two transverse dimensions [7]:
where $E$ is the intracavity field, $P$ is the amplitude of the input pump, $\theta$ is the detuning between the input pump and the closest cavity resonance, $\beta$ is proportional to the Kerr coefficient of the nonlinear material, and the term with the transverse Laplacian $\nabla ^2$ describes diffraction and can be written in either Cartesian or polar coordinates. The time scale has been normalised by $\tau _p$ the mean lifetime of photons in the cavity given by $2L/cT$ for a unidirectional ring cavity and by $4L/cT$ for a Fabry-Perot cavity, with $L$ being the cavity length, $T$ the (intensity) transmission coefficient of the cavity mirrors and $c$ the speed of light in vacuum. The transverse spatial scale ($x$, $y$) has been normalised byIn order to derive some analytical results, we note that LG modes with $m > 0$ can be considered as rings of fixed radius $R_m$ normalised via (2). We therefore express the transverse Laplacian in polar coordinates $(R, \varphi )$ and, as $R$ can be considered a constant, we can write the LLE Eq. (1) in one angular transverse dimension:
2.1 Homogeneous stationary states
For any value of $m$, the homogeneous stationary solutions $F_s$ are obtained from:
where $I_s$ is the intensity of the stationary solution $I_s=|F_s|^2=|E_s|^2$. Once the stationary intensity $I_s$ is selected, the amplitude and phase of the pump field are obtained from (7) implicitly. We highlight the radial dependent detuning term that comes from the OAM associated with the helical phase of the stationary solution (5) and note that for $m\neq 0$, homogeneous stationary states of Eq. (6) correspond to stationary states for the field $E$ that are not homogeneous in the phase $\varphi$.2.2 Turing instabilities on the ring: $m=0$
For $m=0$, $F=E$ and Eq. (6) is equivalent to Eq. (3). Both Eqs. (1) and (3) are well known to display a Turing instability of the homogeneous stationary state. In order to analyze the stability of the solutions we introduce a small perturbation with wavevector $k$: $E = E_s e^{\lambda (k)}\delta E$ and neglect terms nonlinear in $\delta E, \delta E^*$. Performing a linear stability analysis (LSA) we find that above the Turing instability, both Eqs. (1) and (3) with $m=0$ experience the linear growth of a perturbation with wavevector $k$ given by [7,18] :
where $\Delta = \theta + k^2$. It is clear from Eq. (8) that if the square root is imaginary, there are no instabilities since both eigenvalues have negative real part. For a real eigenvalue to be positive, the quantity in the square root has to be larger than one. The instability boundary where the square root in (8) is exactly equal to one, provides a relation between the detuning $\Delta$ (which contains the wavevector $k$) and the stationary intensity $I_s$: This shows that there is an instability threshold in the stationary intensity given by $I_s^c=1/\beta$. For a given $I_s>I_s^c=1/\beta$ the most unstable wavevector is obtained by finding the maximum of the square root in (8) when changing $\Delta$: Above threshold, $N$ peaks appear along the ring separated by a distance given by, or close to, the wavelength of the Turing structure $\Lambda _c = 2 \pi / k_c$. Note that for a ring of circumference $2 \pi R$, the number of peaks is $N = 2 \pi R/\Lambda _c = R \sqrt {2 \beta I_s - \theta }$, and exactly $N$ peaks fit inside a ring of radius $R$ to satisfy the periodic boundary conditions.2.3 Rotating solutions: $m\neq 0$
We now consider the case of pumps carrying OAM, i.e. $m \neq 0$. Above the instability threshold these are seen numerically to form patterns that move around the ring at a constant angular velocity. We start by rearranging Eq. (6) such that the first order derivatives are on the l.h.s.:
Among these travelling waves solutions we can identify Turing patterns for $m \neq 0$ arising from Eq. (15) by analogy with the $m=0$ case. Travelling wave Turing patterns are rotating solutions of Eq. (11) satisfying (8) upon the redefinition of the space-dependent detuning
and most unstable wavevector: The wavelength of the Turing structure $\Lambda _c = 2 \pi / k_c$ and the number of peaks now both depend on the OAM and the radius of the ring. Note that for detuning $\theta$ different from $2 \beta I_s$ and for $m^2/R^2$ small, e.g. for small magnitudes of OAM and large radii, the critical wavevectors from (10) and (17), and hence the number of peaks, are approximately the same. We also note that, at difference from the case of pure propagation [3–6], the number of intensity peaks is crucially dependent on the cavity detuning $\theta$ and the input pump intensity $I_s$. Historically, expressions for the angular velocity similar to (14) had been obtained and applied to rotating domain walls in optical parametric oscillators [8] and used as numerical ansatz for self-trapped necklace-ring beams in a self-focusing nonlinear Schrödinger equation [20].The present analysis is confined to Turing patterns close to the threshold of instability of the homogeneous stationary state. In Appendix B we show that it is possible to obtain equations that describe rotating Turing patterns for generic values of the pump $P_m$, detuning and OAM.
3. Numerical simulations
Although the analysis in the previous section assumed a quasi-1D geometry (rings of fixed radius), all of our numerical simulations are performed in 2D.
3.1 Laguerre-Gaussian pumps
We start by using the well-known split-step Fourier method to numerically integrate Eq. (1) using a Laguerre-Gaussian pump with radial index $p = 0$ [21]:
For our given parameters, $I_s = 1.44, \theta = 1, \beta = 2/3, w_0 = 15.0$ (normalised units), we find that $11$ peaks form on a ring of radius $R = 11.0 \pm 0.5$, as shown in the top panel of Fig. 2. This is in accordance with the closest integer value from our predicted value, using Eq. (18), of $10.5 \pm 0.5$.
From Eq. (14) we predict that the peaks should rotate counter-clockwise at a constant angular velocity $\omega = 0.0164 \pm 0.0003$. To calculate the angular velocity in the numerical results we plot the time evolution of the field at radius $R$, as shown in the bottom panel of Fig. 2, where the diagonal red lines correspond to the peaks of intensity. Angular velocity $\omega = \Delta \varphi / \Delta t = \Delta s / (R \Delta t)$ where $\Delta s$ is the distance a peak travels around the circumference of the circle in a time $\Delta t$. Numerically we find $\omega = 0.0164 \pm 0.0003$, in complete agreement with our analytical results.
We repeated our simulations for an LG pump with $m = -1$, and found that the now peaks rotated clockwise at the same speed, $\omega = 0.0164 \pm 0.0003$, again as predicted from Eq. (14). However, when we measured the angular velocities for $m = 2$, 3, 4 & 5 numerically we found that in each case $\omega = 0.0164 \pm 0.0003$, as shown by the blue line in Fig. 3 i.e. the angular velocity was constant and apparently independent of the OAM of the beam. We can explain this result by noting that the radius of maximum amplitude of the field is OAM-dependent: the red line in Fig. 3 shows our numerical results are in good agreement with the predicted angular velocities calculated using the measured values of $R$ for each OAM.
For any LG mode the radius of maximum amplitude is $r_{max} = w_0 \sqrt {\left | m \right |/2}$. Substituting this into the angular velocity we find $\omega = 2 m /r^2 = \pm 4/w_0^2$. We therefore expect the angular velocity of LG modes to be independent of $m$ but inversely proportional to the beam waist, $w_0$.
To confirm this we numerically integrated Eq. (1) using LG pumps with $m = 1$ and beam waists of $10.0, 15.0, 20.0$, which formed rings of bright spots at radii $19.23, 28.11, 37.5$, and with a beam waist of $25.0$, which formed two rings of bright spots at radii $39.0, 53.0$. Figure 4 shows that we have very good agreement between the angular velocity of the rings measured numerically (blue line) and the predicted values from Eq. (14) using the measured radii (red line).
These measurements simultaneously confirm the direct proportionality of the angular velocity to the OAM index $m$ and inverse proportionality to the square of the radius of the input ring (see Eq. (14). We note that in previous studies of the effect of LG pumps on rotating cavity solitons in semiconductor microresonators [9] and on rotating patterns in photorefractive media in single mirror feedback configurations [10], the pump radius and the OAM index where changed simultaneously leading to weak dependencies of the angular velocity on the OAM index (see Fig. 3. We believe that both these investigations provide support to the universality of Eq. (14) when one includes the changing radius of the input pump of LG modes with different OAM.
3.2 Top-hat pumps carrying OAM
We next consider the case where the input pump has a top-hat shape amplitude multiplied by an azimuthal phase:
Here $P_m$ is a spatially independent complex amplitude, $S$ and $r_t$ control the steepness of the sides and the radius of the top-hat, respectively, and $m$ is an integer corresponding to the topological charge of the optical vortex, as before.Diffraction due to the finite size of the pump induces concentric rings, whose amplitude decreases from the outer ring inwards for $m=0$, as shown in Fig. 5(a). The amplitude of the outermost ring increases with the steepness of the sides of the pump and this can allow its intensity to trigger the Turing instability, see the red line in Fig. 6 in a comparison with $I_s^c=1/\beta$, and an azimuthal pattern forms on the ring. Once the pattern has formed on the outer ring (left panel Fig. 7), we observe a sequence of azimuthal instabilities from the outer to the inner ring. Each patterned ring forms a number of peaks separated by the critical wavelength corresponding to the radius of the particular diffraction ring. The final patterns for $m = 0$ are stationary and close to the centre of the pump they have a hexagonal structure, see Fig. 5(b), typical of that found for plane-waves [7]. Patterns on the outermost diffraction rings, however, are similar to the well-known daisy or sunflower patterns as observed, for example, in VCSELs with an electronic pump with a steep oxide confinement [22–26]. When the pump carries OAM, i.e. $m \ne 0$, the phase at the centre of the pump is undefined and hence the field at the origin has to be zero, as is typical for Laguerre-Gaussian modes [27]. The physical effect of this on-axis vortex is the induction of diffractive rings close to the centre of the beam, see the dashed blue line in Fig. 6. In this case the intensity of the inner ring may exceed $I_s^c$ and undergo an azimuthal Turing instability. The distance between the concentric rings is close to the critical wavelength $\Lambda _c$ along the radial direction. Once the Turing pattern has formed on the inner ring we observe a sequence of azimuthal instabilities taking place from the inner to the outer ring. As before, azimuthal peaks are separated by the critical wavelength corresponding to the radial length of the particular diffraction ring, modified by the factor $m^2/R^2$ as shown in Eq. (17). Note: The effect of the factor $m^2/R^2$ on the wavelength is stronger towards the centre of the beam where the radius is smaller. Moreover, the size of the central vortex, and hence radius of the first diffraction ring, increases with increasing OAM, $m$. These effects alter the radial wavevector in the vicinity of the vortex and can prevent regular patterns from forming on the innermost rings.
As mentioned above, the steepness of the top-hat pump, determined by $S$ in Eq. (20), and also of the optical vortex, affects the relative intensity of the diffraction rings. By careful choice of steepness, and/or the OAM of the pump we can control if the Turing instability first occurs on the inner or the outer ring, or even on both simultaneously. This is demonstrated in Fig. 7 for a top-hat pump with OAM $m = 1$ and steepnesses of $S = 3.0, 1.8, 1.9$ from left to right. For $S = 3.0$ the sides of the pump are steeper than the vortex and so the pattern forms first on the outer ring (left image). For $S = 1.8$ the sides of the vortex are steeper than the pump and so the pattern forms first on the inner ring (middle image). For $S = 1.9$ the steepnesses are almost balanced and so the pattern can form on the inner and outer rings almost simultaneously (right image). Ring by ring azimuthal instabilities can occur from the outside to the inside, or vice versa depending on where the pattern is first formed, or may even collide.
Whatever the mechanism of ring excitation (innermost or outermost), for $m \ne 0$ Turing patterns appear on the first excited ring and start to rotate before the next ring is excited. This leads to concentric Turing pattern rings where the peaks are seen to rotate at exactly the angular frequency given by (14) with $R$ equal to the radius of the specific diffraction ring. It is interesting to compare the final phase and final far-field distributions for the $m=0$ and $m=1$ cases, as shown in Fig. 8. In both cases there are phase peaks where the local intensity also displays peaks. For $m \neq 0$, the central singularity and the line of phase discontinuity corresponding to the OAM of the input pump are clearly visible. Note that there are no further phase discontinuities in spite of the presence of rotating rings. For $m=0$ the hexagonal structure of the near-field is reflected in the peaks of the far-field distribution (see Fig. 8(b)) while the intensity peaks in the near field for $m \neq 0$ are separated by the critical wavelength $\Lambda _c$ but do not display any specific 2D geometry since they are located on rotating rings (see Fig. 8(d)). Each ring is decoupled from the rest of the structure, meaning that they behave as independent 1D azimuthal structures although embedded in a fully 2D field. These rings maintain a fixed structure and rotation speed ad infinitum, unlike the case of optical vortices propagating in Kerr media.
In Fig. 9 we plot the angular velocity of each ring that forms versus its radius for top-hat pumps with $m = 1$ (red/cyan), $m = 2$ (blue/magenta) and $m = 3$ (green/orange). The numerical results (solid lines) show excellent agreement with the analytical results (dashed line) calculated using (14) provided with the measured radii of the rotating rings. This confirms that scalar pumps carrying OAM $m$ form independent Turing patterns on concentric rings of radius $R$ each rotating with constant angular velocity $\omega = 2 m/R^2$ (see Fig. 5). For completion, we note that the dynamics leading to the asymptotic ring rotation for $m \ne 0$ is much faster than that of $m = 0$ where the formation of hexagonal Turing structures is slowed down by the circular symmetry of the input pump.
4. Fully-structured pumps
Vector, or fully structured light (FSL), beams [11–13] have attracted increasing attention for a number of applications. These beams consist of a vector superposition of two scalar orbital angular momentum (OAM) carrying Laguerre-Gaussian (LG) eigenmodes with orthogonal circular polarizations:
4.1 Cylindrical vector beam pumps
If the two modes have equal but opposite OAM the resultant beam is know as a cylindrical vector (CV) beam [11,13]
In Fig. 11 we plot the angular velocity (normalised to $2 m/R^2$) versus the bias parameter $\gamma$ for cylindrical vector (CV) beams with $|m| = 1$ (red), $|m| = 2$ (blue) and $|m| = 3$ (green). By controlling the bias ($\gamma$) between the two eigenmodes of the CV beam, we can fine tune the angular velocity of the output field from
(Recall that LG modes all have the same angular velocity for any given beam waist $w_0$.)4.2 Poincaré pumps
If the two modes have different magnitudes of OAM, the resultant beam is know as a Poincaré beam [12]:
4.3 Optical peppermill
Up until now we have only considered vector beams with some degree of spatial overlap. By considering eigenmodes with significantly different transverse radii such that there is very little interaction between them we can create, for example, counter-rotating rings of spots such as the optical peppermill shown in Fig. 13. In this case the pump consists of a left -circularly polarized mode with $m_L = -1$ and right-circularly polarized mode with $m_R = 8$. The output field has two rings of peaks, with the outer ring rotating counter-clockwise and the inner clockwise with the same angular velocity $\omega = 0.0375\pm 0.0015$. In principle we can produce even more complex superpositions of modes as shown:
5. Conclusion
We have demonstrated formation and rotation of spatio-temporal patterns in self-focussing nonlinear optical cavities pumped by beams carrying orbital angular momentum, $m$. For scalar pumps we see the formation of a ring, or concentric rings, around an optical vortex that rotate at angular velocity $\omega$. Using a 1D Lugiato-Lefever model we find that $\omega = \frac {2 m}{R^2}$, where $R$ is the radius of each ring. For a 1D azimuthal model this formula is exact but we confirm numerically that these angular velocities extend to the 2D case and demonstrate this using input pumps that are Laguerre-Gaussian modes and top-hat shaped pumps with OAM. We note that the radius of maximum intensity of an LG mode scales with the OAM s.t. LG beams with the same beam waist $w_0$ have the same angular velocity, $\omega = \pm 4/w_0^2$. This means that we can control the angular velocity of the patterns by the choice of:
Note that the numer of independent concentric rings that can form inside the top-hat depends on its diameter and the Turing pattern wavelength. Our analysis confirms earlier results on rotating domain walls in optical parametric oscillators and self-trapped necklace-ring beams in a self-focusing nonlinear Schrödinger equation.Further control over the angular velocity of the pattern can be achieved using vector pumps with orthogonally polarization eigenmodes with good spatial overlap.
- • Using cylindrical vector beams, that have eigenmodes with equal and opposite OAM $m$, controlling the relative weightings of the eigenmodes, the bias, allows the angular velocity to range from $- 2 m/R^2 \le \omega \le 2 m/R^2$.
- • Using Poincaré beams, that have eigenmodes with different magnitudes of OAM $m_L, m_R$, the resultant angular velocity is $\omega = (m_L + m_R)/R^2$.
Appendix A. Retarded time transformation
Equation (11) is the LL equation on a ring for the field $F(\varphi , t)$. By using the angular velocity (14) it is possible to introduce retarded time transformations
s.t we can writeAppendix B. Rotating Turing patterns away from threshold
Here we investigate azimuthal Turing patterns on a ring due to pump fields carrying OAM well above threshold. From the analysis close to threshold, Turing patterns are spatially modulated structures with wavelength $\Lambda _c=2\pi /k_c$ where $k_c$ is the critical wavevector given by Eq. (17). We consider spatially modulated solutions of the ring LLE (3) of the form
Funding
Leverhulme Trust (RPG-2017-048); Horizon 2020 Framework Programme (721465).
Acknowledgments
We acknowledge support from the Leverhulme Trust Research Project Grant No. RPG-2017-048, the European Training Network ColOpt, which is funded by the European Union (EU) Horizon 2020 programme under the Marie Sklodowska-Curie action, grant agreement 721465, and the Engineering and Physical Sciences Research Council DTA Grant No. EP/M506643/1.
Disclosures
The authors declare that there are no conflicts of interest related to this article.
References
1. A. M. Turing, “The Chemical Basis of Morphogenesis,” Phil. Trans. R. Soc. Lond. B 237(641), 37–72 (1952). [CrossRef]
2. T. B. Benjamin and J. E. Feir, “The disintegration of wave trains on deep water. Part 1. Theory,” J. Fluid Mech. 27(3), 417–430 (1967). [CrossRef]
3. W. J. Firth and D. V. Skryabin, “Optical Solitons Carrying Orbital Angular Momentum,” Phys. Rev. Lett. 79(13), 2450–2453 (1997). [CrossRef]
4. M. S. Bigelow, P. Zerom, and R. W. Boyd, “Breakup of Ring Beams Carrying Orbital Angular Momentum in Sodium Vapor,” Phys. Rev. Lett. 92(8), 083902 (2004). [CrossRef]
5. A. Vinçotte and L. Bergé, “Femtosecond Optical Vortices in Air,” Phys. Rev. Lett. 95(19), 193901 (2005). [CrossRef]
6. L. T. Vuong, T. D. Grow, A. Ishaaya, A. L. Gaeta, G. W. ’t Hooft, E. R. Eliel, and G. Fibich, “Collapse of Optical Vortices,” Phys. Rev. Lett. 96(13), 133901 (2006). [CrossRef]
7. L. A. Lugiato and R. Lefever, “Spatial Dissipative Structures in Passive Optical Systems,” Phys. Rev. Lett. 58(21), 2209–2211 (1987). [CrossRef]
8. G.-L. Oppo, A. J. Scroggie, and W. J. Firth, “Characterization, dynamics and stabilization of diffractive domain walls and dark ring cavity solitons in parametric oscillators,” Phys. Rev. E 63(6), 066209 (2001). [CrossRef]
9. R. Kheradmand, L. A. Lugiato, G. Tissoni, M. Brambilla, and H. Tajalli, “Rotating and Fugitive Cavity Solitons in semiconductor microresonators,” Opt. Express 11(26), 3612–3621 (2003). [CrossRef]
10. V. Caullet, N. Marsal, D. Wolfersberger, and M. Sciamanna, “Vortex Induced Rotation Dynamics of Optical Patterns,” Phys. Rev. Lett. 108(26), 263903 (2012). [CrossRef]
11. Q. Zhan, “Cylindrical vector beams: from mathematical concepts to applications,” Adv. Opt. Photonics 1(1), 1–57 (2009). [CrossRef]
12. A. M. Beckley, T. G. Brown, and M. A. Alonso, “Full Poincaré beams,” Opt. Express 18(10), 10777–10785 (2010). [CrossRef]
13. E. J. Galvez, S. Khadka, W. H. Schubert, and S. Nomoto, “Poincaré-beam patterns produced by nonseparable superpositions of Laguerre-Gauss and polarization modes of light,” Appl. Opt. 51(15), 2925–2934 (2012). [CrossRef]
14. A. V. Nesterov and V. G. Niziev, “Laser beams with axially symmetric polarization,” J. Phys. D: Appl. Phys. 33(15), 1817–1822 (2000). [CrossRef]
15. L. Novotny, M. R. Beversluis, K. S. Youngworth, and T. G. Brown, “Longitudinal Field Modes Probed by Single Molecules,” Phys. Rev. Lett. 86(23), 5251–5254 (2001). [CrossRef]
16. R. Dorn, S. Quabis, and G. Leuchs, “Sharper Focus for a Radially Polarized Light Beam,” Phys. Rev. Lett. 91(23), 233901 (2003). [CrossRef]
17. F. Bouchard, H. Larocque, A. M. Yao, C. Travis, I. De Leon, A. Rubano, E. Karimi, G.-L. Oppo, and R. W. Boyd, “Polarization Shaping for Control of Nonlinear Propagation,” Phys. Rev. Lett. 117(23), 233903 (2016). [CrossRef]
18. A. J. Scroggie, W. J. Firth, G. S. McDonald, M. Tlidi, R. Lefever, and L. A. Lugiato, “Pattern Formation in a Passive Kerr Cavity,” Chaos, Solitons Fractals 4(8-9), 1323–1354 (1994). [CrossRef]
19. W. J. Firth and A. J. Scroggie, “Optical Bullet Holes: Robust Controllable Localized States of a Nonlinear Cavity,” Phys. Rev. Lett. 76(10), 1623–1626 (1996). [CrossRef]
20. M. Soljačić and M. Segev, “Integer and Fractional Angular Momentum Borne on Self-Trapped Necklace-Ring Beams,” Phys. Rev. Lett. 86(3), 420–423 (2001). [CrossRef]
21. S. M. Barnett and R. Zambrini, “Orbital Angular Momentum of Light,” in “Quantum Imaging,” pg. 284, K. I. Kolobov, ed. (Springer, Singapore, 2007).
22. H. Li, T. L. Lucas, J. G. McInerney, and R. A. Morgan, “Transverse modes and patterns of electrically pumped vertical-cavity surface-emitting semiconductor lasers,” Chaos, Solitons Fractals 4(8-9), 1619–1636 (1994). [CrossRef]
23. Y. G. Zhao and J. McInerney, “Transverse-mode control of vertical-cavity surface-emitting lasers,” IEEE J. Quantum Electron. 32(11), 1950–1958 (1996). [CrossRef]
24. S. F. Pereira, M. B. Willemsen, M. P. van Exter, and J. P. Woerdman, “Pinning of daisy modes in optically pumped vertical-cavity surface-emitting lasers,” Appl. Phys. Lett. 73(16), 2239–2241 (1998). [CrossRef]
25. C. Degen, I. Fischer, and W. Elßer, “Transverse modes in oxide confined VCSELs: Influence of pump profile, spatial hole burning, and thermal effects,” Opt. Express 5(3), 38–47 (1999). [CrossRef]
26. F. Lemke, C. Kropla, A. Mischok, R. Brückner, H. Fröb, and K. Leo, “Nonlinearity-induced Laguerre-Gauss modes in organic vertical cavity lasers,” Appl. Phys. Lett. 111(6), 063303 (2017). [CrossRef]
27. A. M. Yao and M. J. Padgett, “Orbital angular momentum: origins, behavior and applications,” Adv. Opt. Photonics 3(2), 161–204 (2011). [CrossRef]
28. J. B. Geddes, J. V. Moloney, E. M. Wright, and W. J. Firth, “Polarisation patterns in a nonlinear cavity,” Opt. Commun. 111(5-6), 623–631 (1994). [CrossRef]
29. C. J. Gibson, P. Bevington, G.-L. Oppo, and A. M. Yao, “Control of polarization rotation in nonlinear propagation of fully structured light,” Phys. Rev. A 97(3), 033832 (2018). [CrossRef]
30. J. Guck, R. Ananthakrishnan, H. Mahmood, T. J. Moon, C. C. Cunningham, and J. Käs, “The Optical Stretcher: A Novel Laser Tool to Micromanipulate Cells,” Biophys. J. 81(2), 767–784 (2001). [CrossRef]
31. G. Baio, G. R. M. Robb, A. M. Yao, and G.-L. Oppo are preparing a manuscript to be called “Optomechanical transport of cold atoms induced by structured light.”