Abstract
We introduce a kind of the spiraling elliptic Hermite-Gaussian solitons in nonlocal nonlinear media without anisotropy, which carries the orbital angular momentum and can rotate in the transverse. The n–th mode of the spiraling elliptic Hermite-Gaussian solitons has n holes nested in the elliptic profile. The analytical spiraling elliptic Hermite-Gaussian solitons solutions are obtained based on the variational approach, which agree well with the numerical simulations. It is found that the critical power and the critical angular velocity for the spiraling elliptic Hermite-Gaussian solitons are the same as the counterpart of the ground mode.
Similar content being viewed by others
Introduction
The nonlinear propagation of optical beams with orbital angular momentum (OAM) has been discussed in recent years. The spiraling beams carrying the OAM can exert forces and torques on the microparticles, which make them rotate1. The technologies associated with OAM, including spatial light modulators and hologram design, have found their own applications ranging from optical tweezers to microscopy2. Spiraling solitons3 with the OAM are usually associated with optical vortices4, 5 and the ring-shaped beams6. Meanwhile, the vortex-free beams with nonzero OAM are well known in linear media, such as an elliptically shaped beam focused by a tilted cylindrical lens7. It was discussed that the OAM has two contributions to the dynamics of elliptic beams in nonlinear self-focusing media8. First, it effectively strengthens the diffraction against self-focusing and can suppress collapse in Kerr media. Second, it preserves the elliptic profile of stably rotating solitons in optical media with collapse-free nonlinearities. Spiraling elliptic solitons have been found to exist in the media with saturable nonlinearity8 and nonlocal nonlinearity9. It was claimed that OAM can result in the effective anisotropic diffraction for the spiraling elliptic beams9. And the deviations from the critical OAM can make the spiraling elliptic beams breathe10. The decrease (increase) of the OAM can make the spiraling elliptic breathers converge (diffract). Introducing linear anisotropy in the nonlinear media, the OAM will not be conserved. Depending on the linear anisotropy of the media, two kinds of evolution behaviors for the dynamic breathers, rotations and molecule-like librations were predicated analytically and confirmed in numerical simulations11.
In this paper, we discuss a kind of spiraling elliptic Hermite-Gaussian solitons in nonlocal nonlinear media, the n–th mode of which has n holes nested in the elliptic profile. In particular, the fundamental mode is the spiraling elliptic solitons8,9,10,11,12. By using the variational approach, we obtained the approximate analytical solutions, which agree with the numerical simulations well. We find the critical angular velocity of such a soliton depends on the initial parameters but does not depend on its order, which has potential applications in the controlling of the optical beams.
Model
The propagation of optical beams in nonlocal cubic nonlinear media can be modeled by the following nonlocal nonlinear Schrödinger equation (NNLSE)13,14,15
where \({\nabla }_{\perp }^{2}=\frac{{\partial }^{2}}{\partial {x}^{2}}+\frac{{\partial }^{2}}{\partial {y}^{2}}\) is the transverse Laplacian operator, ψ(x, y, z) is the complex amplitude envelope, \({\rm{\Delta }}n=\iint R(x-x^{\prime} ,y-y^{\prime} ){|\psi (x^{\prime} ,y^{\prime} )|}^{2}dx^{\prime} dy^{\prime} \) is the light-induced nonlinear refractive index, z is the longitudinal coordinate, x and y are the transverse coordinates, and R is the normalized symmetrical real spatial response function of the media such that ∫∫R(x, y)dxdy = 1. In the strongly nonlocal nonlinear (SNN) media, we only need keep the first two terms of the expansion of Δn. Then, the NNLSE is simplified to the Snyder-Mitchell mode (SMM)13
where \(\gamma =-\frac{1}{2}{\partial }_{x}^{2}R(x,y){|}_{x=0,y=0}=-\frac{1}{2}{\partial }_{y}^{2}R(x,y){|}_{x=\mathrm{0,}y=0}\), P 0 = ∫∫ |ψ(r′)|2d2 r′ is the input optical power, and r′ is the transverse coordinate vector with r′ = x′e x + y′e y . Although the SMM (2) is a phenomenological model, it can keep the main features of the SNN media. For example, the theoretical predictions by the Snyder-Mitchell model13, such as the accessible solitons and the attraction of spatial solitons, have been observed in experiments in the nematic liquid crystal16,17,18 and the lead glass19. It is worth mentioning that the fractional Fourier transform existing the SNN media was predicted by the SMM20, which was also observed in the lead glass.
Optical beam carrying the orbital angular momentum (OAM) has been investigated in the nonlocal nonlinear media modeled by Eq. (1) 9, 10, 21. Hermite soliton clusters in nonlocal nonlinear media have been introduced by Buccoliero et al.22. The spiraling elliptic Hermite-Gaussian beam carrying the OAM is introduced as follows
where H n is the n–order Hermite polynomials, b(z) and c(z) are the semi-axes of the elliptic beam, A n is a parameter in connection with the amplitude of the optical beam, and ϕ = B(z)X 2 + Θ(z)XY + Q(z)Y 2 + ϑ(z) is the phase. The n–order spiraling elliptic Hermite-Gaussian beam has the elliptic profile, and has n holes aligned along the direction of the principal axis of ellipse. Figure 1 shows the lowest order spiraling elliptic Hermite-Gaussian beam as examples. It should be noted that the expression (3) of the elliptical beam is in the rotating coordinate system XYZ, where X = x cos β(z) + y sin β(z), Y = −x sin β(z) + y cos β(z), Z = z. And in the static coordinate system xyz, the optical beam (3) will rotate carrying the OAM during propagation, the angular velocity of which can be obtained as ω = dβ/dz. There exists a significant difference between the spiraling elliptic Hermite-Gaussian beam (3) and the complex-variable-function-Gaussian beams introduced in ref. 21, which is, the phase contains an additional cross term ΘXY in the former case. The fundamental mode of (3) with n = 0 corresponds to the spiraling elliptic solitons discussed in our recent work9. The optical power and the OAM can be obtained by inserting Eq. (3) into the following two formulas P 0 = ∫∫ |ψ(x′, y′)|2 dx′dy′ and \({M}_{0}={\rm{Im}}\iint {\psi }^{\ast }(x\frac{\partial \psi }{\partial y}-y\frac{\partial \psi }{\partial x})dxdy\). For the three lowest orders, as examples, the optical power and the OAM have been obtained inTable 1. And the similar procedure can be employed for other higher-order modes. As can be seen from Table 1, for the fundamental mode of the spiraling elliptic Hermite-Gaussian beams, the OAM results from the cross term ΘXY on its phase, but for other high-order mode both the cross term ΘXY and the Hermite polynomials H n (X + iY) contribute to the OAM.
Variational solution
Based on the variational approach23, Eq. (2) can be expressed as an Euler-Lagrange equation corresponding to the variational principle
with the Lagrangian density given by ref. 24
where h is the Hamiltonian density expressed as
Inserting the trial solution of spiraling elliptic Hermite-Gaussian beam (3) into the Lagrangian \(L={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }ldxdy\), and then, following the standard procedures of the variational approach23, we can obtain that
where the primes indicate derivatives with respect to the variable z. Thus it can be found that the power and the OAM of the system are conservative. In the following analytical calculations, we take the first-order mode of spiraling elliptic Hermite-Gaussian beam (3) as an example. And the similar calculations can be applied to other higher modes. The angular velocity can be obtained as
which reveals that the angular velocity ω is closely related to the OAM. Substituting the trial solution (3) into the Hamiltonian density (6) and carrying out integration \(H={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }hdxdy\), we obtain the Hamiltonian of the system, which is the function of b, c, B, Q and Θ. As did in our previous work9, by the substitution of Eqs (8), (9) and (10), the Hamiltonian is expressed by H(b′, c′, b.c), which is the summation of the generalized kinetic energy T and the generalized potential energy V. The generalized kinetic energy T is a quadratic function of the generalized velocity b′ and c′. If we assume that T = 0, i.e. b′ = c′ = 0, we can obtain the potential energy V = H(b, c) as follows
Solitons can be found as the extrema of the potential V(b, c). By setting ∂V/∂b = 0 and ∂V/∂c = 0, we can obtain the critical power and the critical OAM
respectively, and M c = P c σ c , where ρ = b/c represents the ellipticity of the elliptic beam. Substitution of the critical power (13) and the critical OAM (14) into the expression of the angular velocity (11) yields
which shows that the spiraling elliptic Hermite-Gaussian solitons make constant-angular rotations. Then the period of rotation can be obtained as T = 2πb 2 c 2/(b 2 + c 2).
For high-order mode of the spiraling elliptic Hermite-Gaussian beams, the critical power and the critical OAM can be calculated by the same process. It is found that the critical power and the critical angular velocity are the same as the counterpart of the ground mode, i.e. Eqs (13) and (15) respectively. While, it is different for the other critical parameters, for example, when w m = 20, b = 1.2, c = 0.8, we obtain σ c = 2.76925 and Θ c = 0.434028 for the second-order mode of the spiraling elliptic Hermite-Gaussian solitons.
Numerical simulation
Here we take the Gaussian function as the spatial nonlocal response function25, 26, i.e.
then the parameter γ in the SMM (2) is obtained \(\gamma =\mathrm{1/(}\pi {w}_{m}^{4})\). Although the Gaussian response function is phenomenological, which does not exist in any physical system, it can be employed to obtain the analytic solution of the NNLSE (1). In addition, for any reasonable response function the physical properties do not depend strongly on its shape. The generic properties of different types of response functions have been studied by Wyller et al. in terms of modulational instability27.
The method of numerical simulation used here is the split-step Fourier method28 using the variational solution (3) as the input. The evolution of the first-order mode of the spiraling elliptic Hermite-Gaussian solitons is shown in Fig. 2, where the parameters are b = 1.2, c = 0.8, and w m = 20. The second-order-moment beam widths based on the variational solution are obtained as \({w}_{x}={[\mathrm{(2/3)(}{b}^{-2}{\cos }^{2}{\omega }_{c}z+{c}^{-2}{\sin }^{2}{\omega }_{c}z)]}^{-\mathrm{1/2}}\) and \({w}_{y}={[\mathrm{(2}/3)({c}^{-2}{\cos }^{2}{\omega }_{c}z+{b}^{-2}{\sin }^{2}{\omega }_{c}z)]}^{-\mathrm{1/2}}\) along the x and y directions, which agrees with the numerical simulations very well as shown in Fig. 2(f). It reveals that the variational solution, including any order spiraling elliptic Hermite-Gaussian solitons is valid for the SNN media, as shown in Figs 3 and 4.
To address the stability of the spiraling elliptic Hermite-Gaussian solitons, we performed numerical simulations of Eq. (2) by employing the initial condition as [1 + εf(x, y)]ψ(x, y), where ψ(x, y), f(x, y) are spiraling elliptic Hermite-Gaussian solitons and the random function with maximum amplitude less than 1, and ε denotes the perturbation parameter. Figure 5 presents the nonlinear propagations of the first-order mode of the spiraling elliptic Hermite-Gaussian solitons with 20% random noises, where we can find the profiles remain invariant up to z = 20. Other high-order modes of the spiraling elliptic Hermite-Gaussian solitons exhibit similar dynamics with added random noises. Of course, the solitons can propagate much farther with random noises than we did in the simulation, which in fact shows the spiraling elliptic Hermite-Gaussian solitons are stable in the nonlocal nonlinear media.
Conclusion
We have introduced a kind of spiraling elliptic Hermite-Gaussian solitons in nonlocal nonlinear media, which carries the the orbital angular momentum and has the elliptic profile with n holes aligning along the direction of the principal axis of ellipse for the n–th order mode. Based on the variational approach, we obtained the approximate analytical solutions, which agree with the numerical simulations well. It was found that the critical power and the critical angular velocity are the same as the counterpart of the ground mode, irrespective of the order n.
References
Molina Terriz, G., Torres, J. P. & Torner, L. Twisted photons. Nat. Phys 3, 305 (2007).
Franke Arnold, S., Allen, L., Padgett, M. Advances in optical angular momentum. Laser & Photon. Rev. 1–15 (2008).
Desyatnikov, A. S., Kivshar, Y. S. & Torner, L. Optical vortices and vortex solitons. Prog. Opt. 47, 291 (2005).
Soskin, M. S. & Vasnetsov, M. V. Singular optics. Prog. Opt. 42, 219 (2001).
Firth, W. J. & Skryabin, D. V. Optical Solitons Carrying Orbital Angular Momentum. Phys. Rev. Lett. 79, 2450 (1997).
Soljacic, M. & Segev, M. Angular Momentum Borne on Self-Trapped Necklace-Ring Beams. Phys. Rev. Lett. 86, 420 (2001).
Courtial, J. et al. Gaussian beams with very high orbital angular momentum. Opt. Commun. 144, 210 (1997).
Desyatnikov, A. S., Buccoliero, D., Dennis, M. R. & Kivshar, Y. S. Suppression of Collapse for Spiraling Elliptic Solitons. Phys. Rev. Lett. 104, 053902 (2010).
Liang, G. & Guo, Q. Spiraling elliptic solitons in nonlocal nonlinear media without anisotropy. Phys. Rev. A 88, 043825 (2013).
Liang, G. et al. Spiraling elliptic beam in nonlocal nonlinear media. Opt. Express 23, 24612 (2015).
Liang, G., Guo, Q., Shou, Q. & Ren, Z. M. Dynamics of elliptic breathers in saturablenonlinear media with linear anisotropy. J. Opt. 16, 085205 (2014).
Liang, G. & Li, H. G. Polarized vector spiraling elliptic solitons in nonlocal nonlinear media. Opt.Commun. 352, 39 (2015).
Snyder, A. W. & Mitchell, D. J. Accessible Solitons. Science 276, 1538 (1997).
Mitchell, D. J. & Snyder, A. W. Soliton dynamics in a nonlocal medium. J. Opt. Soc. Am. B 16, 236 (1999).
Krolikowski, W., Bang, O., Rasmussen, J. J. & Wyller, J. Modulational instability in nonlocal nonlinear Kerr media. Phys. Rev. E 64, 016612 (2001).
Conti, C., Peccianti, M. & Assanto, G. Route to nonlocality and observation of accessible solitons. Phys. Rev. Lett. 91, 073901 (2003).
Conti, C., Peccianti, M. & Assanto, G. Observation of optical spatial solitons in highly nonlocal medium. Phys. Rev. Lett. 92, 113902 (2004).
Peccianti, M., Brzdkiewicz, K. A. & Assanto, G. Nonlocal spatial soliton interactions in nematic liquid crystals. Opt. Lett. 27, 1460 (2002).
Rotschild, C., Cohen, O., Manela, O., Segev, M. & Carmon, T. Solitons in Nonlinear Media with an Infinite Range of Nonlocality: First Observation of Coherent Elliptic Solitons and of Vortex-Ring Solitons. Phys. Rev. Lett. 95, 213904 (2005).
Lu, D. Q. et al. Self induced fractional Fourier transform and revivable higher order spatial solitons in strongly nonlocal nonlinear media. Phys. Rev. A. 78, 043815 (2008).
Deng, D. M., Guo, Q. & Hu, W. Complex variable function Gaussian solitons. Opt. Lett. 34, 43 (2009).
Buccoliero, D., Desyatnikov, A. S., Krolikowski, W. & Kivshar, Y. S. Laguerre and Hermite Soliton Clusters in Nonlocal Nonlinear Media. Phys. Rev. Lett. 98, 053901 (2007).
Anderson, D. Variational approach to nonlinear pulse propagation in optical fibers. Phys. Rev. A 27, 3135 (1983).
Guo, Q., Luo, B. & Chi, S. Optical beams in sub-strongly nonlocal nonlinear media: A variational solution. Opt. Commun. 259, 336 (2006).
Dai, Z. P. & Guo, Q. Spatial soliton switching in strongly nonlocal media with longitudinally increasing optical lattices. J. Opt. Soc. Am. B 28, 134 (2011).
Dai, Z. P., Yang, Z. J., Ling, X. H., Zhang, S. M. & Pang, Z. G. Interaction trajectory of solitons in nonlinear media with an arbitrary degree of nonlocality. Ann. Phys. 366, 13 (2016).
Wyller, J., Krolikowski, W., Bang, O. & Rasmussen, J. J. Generic features of modulational instability in nonlocal Kerr media. Phys. Rev. E. 66, 066615 (2002).
Agrawal G. P. In Nonlinear Fiber Optics 3rd edn 51–55 (Academic, San Diego, 2001).
Acknowledgements
This research was supported by the National Natural Science Foundation of China, Grant No. 11604199, the Key Research Fund of Higher Education of Henan Province, China (16A140030), the Scientific Research Fund of Hunan Provincial Education Department (Grant No. 17B038), and the Hunan Provincial Key Laboratory of Intelligent Information Processing and Application, Hengyang, 421002, China.
Author information
Authors and Affiliations
Contributions
Z.D. initiated the study. G.L. conducted the work in the whole process. Z.D. read the manuscript and gave the valuable suggestions. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Liang, G., Dai, Z. Spiraling elliptic Hermite-Gaussian solitons in nonlocal nonlinear media without anisotropy. Sci Rep 7, 3234 (2017). https://doi.org/10.1038/s41598-017-03669-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-03669-x
This article is cited by
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.