Spiraling elliptic Hermite-Gaussian solitons in nonlocal nonlinear media without anisotropy

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.

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 rotate 1 . The technologies associated with OAM, including spatial light modulators and hologram design, have found their own applications ranging from optical tweezers to microscopy 2 . Spiraling solitons 3 with the OAM are usually associated with optical vortices 4,5 and the ring-shaped beams 6 . 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 lens 7 . It was discussed that the OAM has two contributions to the dynamics of elliptic beams in nonlinear self-focusing media 8 . 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 nonlinearity 8 and nonlocal nonlinearity 9 . It was claimed that OAM can result in the effective anisotropic diffraction for the spiraling elliptic beams 9 . And the deviations from the critical OAM can make the spiraling elliptic beams breathe 10 . 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 simulations 11 .
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 solitons [8][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]  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  [16][17][18] and the lead glass 19 . It is worth mentioning that the fractional Fourier transform existing the SNN media was predicted by the SMM 20 , 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 . Table 1. Parameters of the spiraling elliptic Hermite-Gaussian beam for different orders of n.  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) 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 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 work 9 . 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 . 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 approach 23 where h is the Hamiltonian density expressed as    2 2 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)   H hdxdy, we obtain the Hamiltonian of the system, which is the function of b, c, B, Q and Θ. As did in our previous work 9 , 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  (11) yields 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)

Numerical simulation
Here we take the Gaussian function as the spatial nonlocal response function 25,26  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 instability 27 .
The method of numerical simulation used here is the split-step Fourier method 28 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 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.  Fig. 1(b)] with 20% random noises (a), and the profile at the propagation distance z = 20 (b).