Stability of vortex solitons in thermal nonlinear media with cylindrical symmetry

We analyze the salient features of vortex-ring solitons supported by cylindrically symmetric media with nonlocal thermal nonlinearity. We discover the existence of a maximum allowed topological charge for such vortex solitons to be stable on propagation: Only vortex-ring solitons with topological charge m<=2 are found to be stable. This remarkable result holds independently of the radius of the sample.

In many optical materials, a high-power light beam modifies the refractive index with the local light intensity. Nevertheless, in some materials the response can be highly nonlocal, so that the nonlinear contribution to the refractive index depends on the intensity distribution in the entire transverse plane, on a characteristic spatial scale that depends on the type of nonlocality. Different types of nonlocalities are encountered in actual materials, such as photorefractive, liquid crystals or plasmas. Such nonlocalities can have a profound effect on the properties of nonlinear excitations [1][2][3][4]. In two transverse dimensions, nonlocality stabilizes not only fundamental solitons, but also more complex structures, such as multipolemode solitons, studied theoretically in Refs. [5][6][7] and observed in Ref. [8], rotating azimuthally modulated beams connecting multipole and vortex soliton families [6,9], and vortex-ring solitons [10][11][12][13][14][15]. Notice that, because of their particular shape, the latter tend to be highly prone to azimuthal instabilities in most local materials [16]. Vortex solitons with single topological charge have been recently observed in a selffocusing thermal medium [17]. In such materials long-range soliton interactions take place [18], while the geometry of the sample affects the beam trajectory [19] and determines the shape of the induced refractive index profile [20]. Thermal nonlinearities occur in many optical media, even though they become dominant only under specific conditions in materials with high enough light absorption coefficient, large thermo-optic effects, etc. Thermal nonlinearities are intrinsically nonlocal, with the range of nonlocality being determined by the geometry of the sample and its dimensions.
The fact that nonlocality stabilizes vortex solitons is well established [10][11][12]. However, the shape of the response function of the medium, which depends on the physical mechanism of nonlocality, and its inner scale, appears to be crucial factor to set a maximal possible charge of stable vortex solitons. For example, no limits have been found for the charge of stable nonlinear vortices in materials with a Gaussian-shaped response function [11], while in media described by a Helmholtz-type response function (such as liquid crystals) only singlecharge vortex solitons are stable [12]. In this paper we discover that in thermal nonlinear medium with a cylindrical symmetry only vortex solitons with topological charge can be stable under propagation. This finding is the outcome of a rigorous linear stability analysis and it is confirmed by accurate, direct simulations of the propagation of perturbed vortex soliton stationary solutions.

m ≤
We consider the propagation of a laser beam along the ξ axis in a focusing thermal medium of circular cross-section described by the following system of equations: A is the dimensionless light field amplitude; is proportional to the nonlinear change of the refractive index 0 n ; are the optical absorption, thermo-optic, and thermal conductivity coefficients, respectively; the radial and longitudinal coordinates are scaled to the characteristic beam radius and the , respectively; φ is the azimuthal angle. We solved the system (1) with the boundary conditions 2 0 0 k r , r R q n → = 0, where R is the radius of the sample. Such boundary conditions are consistent with typical experimental conditions where a laser beam is focused to a typical radius of tens of microns and the rod radius amounts to a few millimeters. We assume that the boundary of the medium is kept at a fixed temperature. Notice that a similar physical model was considered by Kruglov and co-workers in Ref. [10], where the very fact that azimuthal instabilities of nonlinear vortex beams can be suppressed in thermal nonlinear media was predicted.
Under such conditions light launched along the axis of cylindrical rod and experiencing slight absorption, raises the temperature of the medium, while heat diffusion causes temperature redistribution in the entire sample that reaches a steady state after a relatively long time interval, which is proportional to and inversely proportional to the thermal conductivity coefficient. Increasing the temperature in a material with positive thermo-optic coefficient leads to a refractive index growth and facilitates light trapping in the heated regions, not only for simplest bell-shaped beams but also for topologically nontrivial beams carrying phase singularities. where the response function of the thermal medium is given by 0 for and 0 for r . As mentioned above, the shape of the response function is set by the geometry of the sample and depends on its transverse extent. Note that this is in a sharp contrast to materials with Helmholtz [7,12] and Gaussian [11] response functions possessing an intrinsic nonlocality scale. In our computer simulations reported here we set , a value which closely resembles actual experimental conditions. However, it is worth stressing that all the results obtained remain qualitatively valid for other values.
We searched for stationary solutions of Eq. (1) in the form , where is a real function describing the profile of radially symmetric nonlinear wave solutions, b is the propagation constant, and m is the winding number, or topological charge. The profiles of fundamental solitons ( and vortices are found numerically with a standard relaxation algorithm. To analyze the soliton stability, we look for perturbed solutions in the form , where are the perturbation components that can grow with a complex rate δ upon propagation, and k is the azimuthal perturbation index. Substitution of the light field in such form into Eq. (1) and linearization yields the eigenvalue problem: is the refractive index perturbation corresponding to the azimuthal index k , while at r and k G r at r . We solved the eigenvalue problem (2) numerically. . At fixed b , the vortex-ring radius gradually increases with growth of m , while the peak amplitude decreases. The energy flow and the integral radius, defined as have larger values for solitons with higher topological charges m at fixed b . While soliton profiles are strongly localized (provided that b is not too small), the temperature and the refractive index distributions are much broader and cover the whole cross-section of the rod. Such refractive index distribution forms a waveguide that traps light near the axis of the cylindrical sample. Notice that the refractive index profile features an almost flat plateau around the point for vortex solitons with ( Fig. 1(b)), while for fundamental soliton it exhibits a cone-like profile instead. The width of the plateau increases with growth of the topological charge. Increasing b results in a monotonic growth of the energy flow U for any m (Fig. 2(a)) and leads to a progressive soliton localization, so that the integral soliton radius monotonically decreases with b (Fig. 2(b)).
Besides the simplest solitons shown in Fig. 1, we found for any value of topological charge a number of higher-order soliton families for which the field has one or multiple nodes along the radial direction (not shown here). Our stability analysis has proved that all such solutions are unstable. This is in clear contrast, e.g., to media with a Gaussian response function where ring-like solitons with nodes can be stable [11].  The rigorous linear stability analysis that we have conducted shows that fundamental solitons are stable in the entire existence domain, which is consistent with the so-called Vakhitov-Kolokolov stability criterion, since for any b value. The central result of this paper is that vortex solitons can be stable only if their topological charges fulfil . Our results indicate that such vortex solitons are stable in the entire existence domain, irrespectively of the radius R of the medium. Therefore, this is a rather unique constraint. For example, recall that in media with a Helmholtz-type response function only charge-one vortex solitons can be stable [12], while in the case of Gaussian response there are no restrictions on the charge of stable vortex solitons [11]. On physical grounds, the shrinking of the stability / dU db > 0 2 m ≤