---
abstract: 'We study the nonlinear optical propagation of two different classes of space-varying polarized light beams – radially symmetric vector beams and Poincaré beams with lemon and star topologies – in a rubidium vapour cell. Unlike Laguerre-Gauss and other types of beams that experience modulational instabilities, we observe that their propagation is not marked by beam breakup while still exhibiting traits such as nonlinear confinement and self-focusing. Our results suggest that by tailoring the spatial structure of the polarization, the effects of nonlinear propagation can be effectively controlled. These findings provide a novel approach to transport high-power light beams in nonlinear media with controllable distortions to their spatial structure and polarization properties.'
author:
- Frédéric Bouchard
- Hugo Larocque
- 'Alison M. Yao'
- Christopher Travis
- Israel De Leon
- Andrea Rubano
- Ebrahim Karimi
- 'Gian-Luca Oppo'
- 'Robert W. Boyd'
title: Polarization shaping for control of nonlinear propagation
---

*Introduction:* Light beams that can propagate without significant change to their spatial profile are of interest for modern optical technologies and high-power laser systems. Self-trapped light filaments, or spatial solitons, are formed when their spreading due to linear diffraction is carefully balanced by a self-focusing (Kerr) nonlinearity that causes the beam to narrow. Due to their potential to carry an increased information content, there has been significant interest in the formation of spatial solitons carrying orbital angular momentum (OAM) [@firth:97; @desyatnikov:01; @musslimani:01; @soljacic:01; @bigelow:02; @bigelow:04]. OAM-carrying beams are characterized by an azimuthal phase dependence of the form $\exp{(i\ell{\varphi})}$ [@allen:92], where the integer $\ell$ corresponds to the topological charge of the phase singularity present at the beam centre (e.g. see [@yao:11] and references therein). Such beams include the Laguerre-Gauss (LG) modes [@barnett:07], which are solutions to the paraxial wave equation.

It is well known that in (2+1) dimensions, spatial solitons are unstable in homogeneous Kerr media [@rasmussen:86]. One way to increase their stability is to use a saturable self-focusing medium to prevent the catastrophic collapse due to self-focusing. By using an intensity dependent nonlinear refractive index, $n_2$, it becomes possible to balance out the effects of self-focusing and diffraction. However, even in the case of saturable self-focusing media, it is known that optical beams carrying OAM will fragment into several solitons possessing particle-like attributes [@desyatnikov:01; @malomed:05]. In particular, a scalar beam carrying an OAM value of $\ell$ is predicted to break up into $2\ell$ daughter solitons [@firth:97; @bigelow:04].

In comparison with scalar ring solitons that carry a definite non-zero OAM, it has been suggested that stability can be increased by using two beams with opposite OAM to produce a beam with a net zero OAM [@soljacic:01; @bigelow:02]. For example, [“necklace”]{} (petal) beams, which consist of a scalar superposition of two modes with equal and opposite OAM, have been shown to exhibit quasi-stable propagation in a self-focusing medium although they expand upon propagation [@soljacic:98]. For vectorial superpositions of beams carrying OAM, however, the resulting *vector solitons* have been shown theoretically to exhibit quasi-stable propagation for much larger distances than the corresponding scalar vortex solitons and necklace beams [@desyatnikov:01; @bigelow:02].

Vector vortex beams are fully correlated solutions to the vector paraxial wave equation that have space-varying polarization distributions. Cylindrical vector (CV) beams are a subclass of vector beams with an axially symmetric polarization profile about the beam’s propagation axis [@zhan:09; @galvez:12]. Examples, include *radial*, *azimuthal* and *spiral* polarization distributions. Another class of light beams with non-uniform polarization structure is that of full Poincaré beams [@beckley:10; @galvez:12]. These are of intrinsic interest because they carry polarization singularities. Such beams typically consist of a superposition of two orthogonally polarized LG modes of different orbital angular momenta [@beckley:10; @galvez:12] and thus carry a net value of angular momenta.

In this Letter, we demonstrate, both experimentally and numerically, the stable propagation of space-varying polarized light beams in a saturable self-focusing nonlinear medium. More specifically, our study focuses on vector vortex and Poincaré beams traveling through rubidium vapour. We compare the intensity and polarization distributions of the beams at the entrance and exit of the nonlinear cell. This allows us to see how beam-break up is affected both by the net OAM of the beam and by its polarization distribution.

*Theory:* Light beams with spatially inhomogeneous polarization distributions can be obtained by a superposition of two spatial transverse modes, $E_1$ and $E_2$, with orthogonal polarizations $$\begin{aligned}
\label{eq:superposition}
    \mathbf{E}(r,{\varphi},z)=E_1(r,{\varphi},z)\,\mathbf{e}_1+E_2(r,{\varphi},z)\,\mathbf{e}_2,\end{aligned}$$ where $\mathbf{e}_1$ and $\mathbf{e}_2$ are orthonormal polarization vectors, and $r,{\varphi},z$ are the cylindrical coordinates. Here we adopt the circular polarization basis, i.e. $\mathbf{e}_1=\mathbf{e}_\text{L}$ and $\mathbf{e}_2=\mathbf{e}_\text{R}$, and the LG basis for the spatial transverse modes [@beckley:10]. If the two beams have equal but opposite OAM, the polarization of the beam varies along the azimuthal coordinate and, if the beams are equally weighted, spans the equator of the Poincaré sphere. These radially symmetric *vector vortex* beams [@zhan:09] can have polarization distributions that are radial, azimuthal (see Fig. \[fig:beams\]-(a), (b)) or spiral. If $E_1$ and $E_2$ carry a zero and a non-zero OAM value, respectively, the resulting full Poincaré beam [@beckley:10] has a polarization that varies in both the angular and radial coordinates [@galvez:12] and covers all polarization states on the Poincaré sphere [@beckley:10]. The state of elliptic polarization varies with position [@dennis:09] and polarization singularities occur at C-points where the azimuth is not defined and the polarization is circular [@nye:83b], and along L-lines where the polarization is linear and its handedness is not defined [@nye:83a]. C-points can have three fundamental polarization topologies classified by the index $\eta$ and the number of polarization lines that terminate at the singularity: these topologies are known as “star” ($\eta=-1/2$, three lines), “lemon” ($\eta=1/2$, one line) and “monstar” ($\eta=1/2$, infinitely many, with three straight, lines) [@nye:99; @dennis:09]. Examples of lemon and star topologies are shown in Fig.\[fig:beams\] (c) and (d).

![(a) Radial and (b) azimuthal vector beams. (c) Lemon and (d) star Poincaré beams. Red (blue) ellipses correspond to left (right) circular polarization.[]{data-label="fig:beams"}](fig1.png){width="0.9\columnwidth"}

We simulate propagation through the medium using a two-dimensional nonlinear Schrödinger equation with saturable self-focusing nonlinerity, derivable from the two-level model, under the slowly varying envelope and paraxial approximations and normalized to dimensionless quantities, $\rho=r/w_0$ and $\zeta={z}/({2z_R})$, where $w_0$ is the beam waist and $z_R$ is the Rayleigh range of the beam [@firth:97; @bigelow:02; @desyatnikov:01]. As we are dealing with vector beams, our model consists of two coupled equations that interact through the cross phase modulation (XPM) term characterized by the parameter $\nu$, which takes the value of $2$ for circularly polarized beams [@Agrawal] $$\begin{aligned}
\frac{\partial E_1}{\partial \zeta} = \frac{i}{2}\nabla^2_\bot E_1 +i\mu \; \frac{|E_1|^2+\nu \, |E_2|^2}{1+\sigma(|E_1|^2+ \nu \, |E_2|^2)} E_1, 
    \label{eq:coupled} \\
 \frac{\partial E_2}{\partial \zeta} = \frac{i}{2}\nabla^2_\bot E_2 +i\mu \; \frac{|E_2|^2+\nu \, |E_1|^2}{1+\sigma(|E_2|^2+ \nu \, |E_1|^2)} E_2. \nonumber\end{aligned}$$ The parameters of importance are the nonlinear parameter, $\mu$, and the saturation parameter, $\sigma$, given by: $$\mu = \frac{k_0^2 n_2 P_0}{n_0}  \; ; \qquad\quad \sigma = \frac{2P_0}{I_{sat} w_0^2} \, ,
  \label{eq:gammasigma}$$ where $k_0$ is the free-space wavenumber, $n_0$ and $n_2$ are the linear and nonlinear refractive indices ($n_2 > 0$ for self-focusing), $I_{sat}$ is the saturation intensity, and $P_0$ is the power of the incident laser beam. In the simulations reported below, we have selected $\mu=386$, $\sigma=51.7$ that reproduce the experimental configuration of the natural rubidium (Rb) cell. We performed numerical integrations of the propagation Eqs. (\[eq:coupled\]), using the split-step method with fast Fourier transforms and parameters corresponding to the experiments performed.

![Experimental set-up. A CW wavelength-tunable ($760$ nm - $790$ nm) Toptica diode laser is coupled to a polarization-maintaining single-mode optical fiber. The power of the laser beam is adjusted by means of a half-wave ($\lambda/2$) plate followed by a polarizing beam splitter (PBS) in order to reach the saturation threshold ($I_\text{sat}=5$ Wcm${}^{-2}$) for the $D_2$ resonance line of Rb. A combination of a $\lambda/2$ and a $\lambda/4$ plate is used to generate an arbitrary polarization state, which is then sent to a $q$-plate (QP), with a topical charge of $q=\pm 1/2$, to generate vector vortex and Poincaré beams. This is then focused ($w_0=60~\mu$m) into a thermally controlled Rb atomic vapour cell by means of a lens ($f=150$ mm). The transmitted beam is analyzed with a combination of a $\lambda/4$ plate, a $\lambda/2$ plate and a PBS (polarization tomography), and attenuated using neutral density filters before being recorded by a CCD camera.[]{data-label="fig:exp"}](fig2.pdf){width="0.8\columnwidth"}

![image](fig3.png){width="2\columnwidth"}

*Experiment:* We use a spatially filtered, linearly polarized, tunable CW single-mode diode laser (Toptica DL pro 780, 760 nm – 790 nm) together with a set of half- and quarter-wave plates to generate a Gaussian beam with an arbitrary polarization state, $\cos{\left(\theta/2\right)}\,\mathbf{e}_\text{L}+\exp{(i\chi)}\,\sin{\left(\theta/2\right)}\,\mathbf{e}_\text{R}$, where $\theta$ and $\chi$ are set by the orientation of the waveplates. This beam is converted into a space-varying polarized light beam using a $q$-plate – a slab of patterned liquid crystal – that couples optical spin to OAM [@marrucci:06]. The *unitary* action of a $q$-plate in the circular polarization basis is described by $$\begin{aligned}
\label{eq:qPlateUnitary}
\hat{\text{U}}_q\cdot\left[\begin{matrix} 
    \mathbf{e}_\text{L} \\
    \mathbf{e}_\text{R} \\
    \end{matrix} \right] = \cos{\left(\frac{\delta}{2} \right)}\left[\begin{matrix}
    \mathbf{e}_\text{L} \\
    \mathbf{e}_\text{R} \\
    \end{matrix} \right]+ i \sin{\left(\frac{\delta}{2}\right)} \left[\begin{matrix}
    \mathbf{e}_\text{R} e^{+2 i \left( q {\varphi}+ \alpha_0 \right)} \\
    \mathbf{e}_\text{L} e^{-2 i \left( q {\varphi}+ \alpha_0 \right)} \\
    \end{matrix} \right],\end{aligned}$$ where $q$ is a half-integer number corresponding to the topological charge of the liquid crystal pattern, $\alpha_0$ is the azimuthal orientation of the liquid crystal elements at ${\varphi}=0$ (laboratory frame) and $\delta$ is the optical retardation of the $q$-plate. The parameter $\delta$ can be experimentally adjusted by applying an electric field onto the plate in such a way that the resulting optical retardation corresponds to a half ($\delta=\pi$) or quarter ($\delta=\pi/2$) wavelength [@slussarenko:11]. The generated beam is then focused by a $150$-mm-focal-length lens into a $9$-cm-long cell containing Rb atomic vapour. A detailed depiction of this experimental apparatus is provided in Fig. \[fig:exp\]. The intensity of the incident beam and the temperature of the atomic vapour are set so that the medium exhibits saturable Kerr nonlinearities (powers in the vicinity of $7.44$ mW and a temperature of $95^\circ$C). In order to observe beam breakup of OAM-carrying beams, the laser’s output wavelength was tuned near the $D_2$ transition line of Rb. Moreover, the laser was blue-detuned by less than $0.7$ GHz from the $5S_{1/2}(F=3) \rightarrow 5 P_{3/2} (F=4)$ hyperfine transition. The beam exiting the Rb cell is then imaged using a lens with a focal length of $200$ mm (not shown in the experimental setup) and its polarization distribution is reconstructed using polarization tomography. The tomography is performed using an appropriate sequence of a quarter-wave plate, a half-wave plate, a polarizer and a spatially resolving detector (CCD camera) set at an exposure time of $0.1$ s.

In order to generate vector vortex beams (radial, azimuthal and spiral), a linearly polarized Gaussian beam is sent to a tuned ($\delta=\pi$) $q$-plate of topological charge $q=1/2$ [@cardano:12]. Apart from a global phase, the generated beam will be given by $\left(\text{LG}_{0,-1}(r,{\varphi},z)\, \mathbf{e}_\text{L} + \exp{\left(i \beta \right)}\,\text{LG}_{0,1}(r,{\varphi},z)\,\mathbf{e}_\text{R} \right)/\sqrt{2}$, where $\beta \, \equiv 4 \alpha_0$ depends on the orientation of the $q$-plate with respect to the input polarization. The generated beams correspond to radial ($\beta=0$), azimuthal ($\beta=\pi$) and spiral ($\beta=\pm \pi/2$) vector vortex beams.

To generate the lemon and star Poincaré beams, a circularly polarized Gaussian beam is sent to a perfectly detuned ($\delta=\pi/2$) $q$-plate with $q=1/2$ and $q=-1/2$, respectively [@cardano:13]. Here, $\beta=2\alpha_0$ and does not affect the polarization topology but does cause a rotation in the polarization pattern. Thus, we choose $\beta = 0$, resulting in an output beam of the form $(\text{LG}_{0,0}(r,{\varphi},z)\, \mathbf{e}_\text{L} +\text{LG}_{0,2q}(r,{\varphi},z)\, \mathbf{e}_\text{R})/\sqrt{2}$, again omitting any global phase factors. Note also that monstar topologies cannot be readily generated in the laboratory. This scheme allows us to switch between different structured light beams without altering their intensities.

*Analysis:* It is well-known that azimuthal modulational instabilities associated with the helical phase structure found in beams carrying OAM result in their filamentation as they propagate through self-focusing media [@firth:97; @desyatnikov:01; @bigelow:02]. It has been shown both analytically and numerically, however, that the dominant low-frequency perturbations that typically disrupt ring solitons are inhibited for vector solitons with no net OAM [@desyatnikov:01; @bigelow:02]. Here we confirm this experimentally and numerically by propagating a beam of the form $\cos{\gamma}\;\text{LG}_{0,-1}(r,{\varphi},z)\; \mathbf{e}_\text{L}+\sin{\gamma}\; \exp{(i\beta)}\;\text{LG}_{0,1}(r,{\varphi},z)\;\mathbf{e}_\text{R}$ through a self-focusing medium. When $\gamma = 0,\pi/2$, the beam is a scalar LG mode of $\ell=-1 , +1$ with left, right-hand circular polarization. For $\gamma = \pi/4$ we have a vector vortex beam with linear polarization and for $\gamma = \pi/8, 3 \pi/8$, a vector vortex beam of left, right-hand elliptical polarization, respectively, following the topology of the vector beam with $\gamma = \pi/4$. Note that each beam has the same total intensity.

A comparison of the experimental results with the simulations based on Eqs. (\[eq:coupled\]) is shown in Fig. \[fig:Ratio\] for a vector beam defined by $\beta = -\pi/2$. From our results we can see that the nonlinearity counterbalances diffraction up until the point at which the beams fragment. As expected, we see that the scalar OAM-carrying beams ($\gamma=0, \, \pi/2$) are starting to break-up at the exit of the nonlinear cell. The vector beam ($\gamma = \pi/4$), on the other hand, seems almost unperturbed, both in terms of its amplitude and polarization distribution. Indeed, we can see numerically that the fragmentation point occurs much later for vector beams, $\sim 16$ cm, than for scalar beams, $\sim 9$ cm. For the elliptically polarized vector beams ($\gamma = \pi/8, 3 \pi/8$) the stability length is between the cases of scalar and linear vector vortex beams. In this case the effect of the nonlinear propagation is also evident in the change of polarization distributions: the initial spiral distribution has now become almost azimuthal or radial, respectively (see Fig. \[fig:beams\]). This suggests that vector beams allow for greater control over the nonlinearity.

In spite of the approximations used in the derivation of Eqs. (\[eq:coupled\]), small differences in polarization distributions between the experiment and the numerical simulation are only really evident in the biased mode cases of $\gamma = \pi/8, 3 \pi/8$, thus demonstrating the robustness of the vector vortex beam ($\gamma = \pi/4$).

It has been proposed that the increase in stability seen in vector beams is due to the fact that they carry no net OAM [@bigelow:02]. We therefore repeated our analysis using lemon and star polarization topologies which have a net OAM.

![\[fig:Poincare\] Intensity and polarization distributions of lemon and star topologies after propagating through the Rb cell, numerical (upper row) and experimental (lower row). Red (blue) ellipses correspond to left (right) circular polarization. ](fig4.png){width="\columnwidth"}

The experimental and numerical results are shown in Fig. \[fig:Poincare\] where we have plotted the intensity and the polarization distributions for the lemon and star topologies after propagating through the Rb cell. These results again show that beam break-up has been inhibited when compared to scalar OAM-carrying beams. This result demonstrates that the increased stability is not simply due to a net OAM of zero. As in the case of linear vector vortex beams ($\gamma=\pi/4$) above, the polarization distribution remains unaltered after nonlinear propagation and even preserves the number and kind of polarization singularities.

*Conclusion:* We have compared the propagation of scalar OAM-carrying beams with two different classes of beams with non-uniform transverse polarization distributions in a saturable self-focusing nonlinear medium of Rb vapour. With respect to scalar vortex beams (LG modes), we found that beam break up can be inhibited while nonlinear confinement, self-focusing and polarization distributions are not altered for specific cases of non-uniform spatial polarization, both with and without net OAM. This suggests that the spatial structure of the polarization plays an important role in preventing beam fragmentation. These findings provide a novel approach to transport high-power light beams in nonlinear media with controllable distortions to their spatial structure and polarization properties.

*Acknowledgment:* F.B., H.L., E.K. and R.W.B. acknowledge the support of the Canada Excellence Research Chairs (CERC) program. A. R. acknowledges funding from the European Union (FP7-PEOPLE-2012-CIG, PCIG12-GA-2012-326499-FOXIDUET). E.K. acknowledges the support of the Canada Research Chairs (CRC) program and Canada Foundation for Innovation (CFI).

[26]{}

W. J. Firth and D. Skryabin, Phys. Rev. Lett. [**79**]{}, 2450 (1997).

A. S. Desyatnikov and Y. S. Kivshar, Phys. Rev. Lett. [**87**]{}, 033901 (2001);

Z. H. Musslimani, M. Soljačić, M. Segev, and D. N. Christodoulides Phys. Rev. E [**63**]{}, 066608 (2001).

M. Soljačić, and M. Segev, Phys. Rev. Lett. [**86**]{}, 420 (2001).

M. S. Bigelow, Q-H. Park, and R. W. Boyd, Phys. Rev. E [**66**]{}, 046631 (2002).

M. S. Bigelow, P. Zerom, and R. W. Boyd, Phys. Rev. Lett. [**92**]{}, 083902 (2004).

L. Allen, M. W. Beijersbergen, R. Spreeuw, and J. Woerdman, Phys. Rev. A [**45**]{}, 8185 (1992).

A. M. Yao and M. J. Padgett, Adv. in Optics and Photonics [**3**]{}, 161 (2011).

S. M. Barnett and R. Zambrini, *Quantum Imaging*, edited by K. I. Kolobov (Singapore: Springer) p 284 (2007)

J. J. Rasmussen and K. Rypdal, Phys. Scr. [**33**]{}, 481 (1986).

B. A. Malomed, D. Mihalache, F. Wise, and L. Torner, J. Opt. B: Quant. and Semiclassical Optics [**7**]{}, R53 (2005).

M. Soljačić, S. Sears, and M. Segev, Phys. Rev. Lett. [**81**]{}, 4851 (1998).

Q. Zhan, Adv. in Optics and Photonics [**1**]{}, 1 (2009).

E. J. Galvez, S. Khadka, W. H. Schubert, and S. Nomoto, App. Opt. [**51**]{}, 2925 (2012).

A. M. Beckley, T. G. Brown, and M. A. Alonso, Opt. Express [**18**]{}, 10777 (2010).

M. R. Dennis, K. O’Holleran, and M. J. Padgett, Prog. Opt. [**53**]{}, 293 (2009).

J. F. Nye, Proc. R Soc. Lond. A [**389**]{}, 279 (1983).

J. F. Nye, Proc. R Soc. Lond. A [**387**]{}, 105 (1983).

J. F. Nye, *Natural focusing and fine structure of light: caustics and wave dislocations* (CRC Press, 1999).

G.P. Agrawal, *Nonlinear fiber optics*, Fourth edition, Academic press (2007).

M. Born and E. Wolf, *Principles of Optics*, Pergamon (1959).

B. E. A. Saleh and M. C. Teich, *Fundamentals of Photonics*, Second edition, Wiley-Interscience (2007).

L. Marrucci, C. Manzo, and D. Paparo, Phys. Rev. Lett. [**96**]{}, 163905 (2006).

S. Slussarenko, A. Murauski, T. Du, V. Chigrinov, L. Marrucci, and E. Santamato, Opt. Express [**19**]{}, 4085 (2011).

F. Cardano, E. Karimi, S. Slussarenko, L. Marrucci, C. de Lisio, and E. Santamato, App. Optics [**51**]{}, C1 (2012).

F. Cardano, E. Karimi, L. Marrucci, C. de Lisio, and E. Santamato, Opt. Express [**21**]{}, 8815 (2013).
