Twisted toroidal vortex solitons in inhomogeneous media with repulsive nonlinearity

Surface gap solitons in the Schrödinger equation with quintic nonlinearity and a lattice potential

We demonstrate the existence of surface gap solitons, a special type of asymmetric solitons, in the one-dimensional nonlinear Schrödinger equation with quintic nonlinearity and a periodic linear potential. The nonlinearity is suddenly switched in a step-like fashion in the middle of the transverse spatial region, while the periodic linear potential is chosen in the form of a simple sin 2 lattice. The asymmetric nonlinearities in this work can be realized by the Feshbach resonance in Bose-Einstein condensates or by the photorefractive effect in optics. The major peaks in the gap soliton families are asymmetric and they are located at the position of the jump in nonlinearity (at x = 0). In addition, the major peaks of the two-peak and multi-peak solitons at the position x = 0 are higher than those after that position, at x > 0. And such phenomena are more obvious when the value of chemical potential is large, or when the difference of nonlinearity values across the jump is big. Along the way, linear stability analysis of the surface gap solitons is performed and the stability domains are identified. It is found that in this model, the solitons in the first band gap are mostly stable (excepting narrow domains of instability at the edges of the gap), while those in the second band gap are mostly unstable (excepting extremely narrow domains of stability for fundamental solitons). These findings are also corroborated by direct numerical simulations.

Vortex rings in paraxial laser beams

Interference of a fundamental vortex-free Gaussian beam with a co-propagating plane wave leads to nucleation of a series of vortex rings in the planes transverse to the optical axis; the number of rings grows with vanishing amplitude of the plane wave. In contrast, such interference with a beam carrying on-axis vortex with winding number l results in the formation of |l| rings elongated and gently twisted in propagation direction. The twist handedness of the vortex lines is determined by the interplay between dynamic and geometric phases of the Gaussian beam and the twist angle grows with vanishing amplitude of the plane wave. In the counter-propagating geometry the vortex rings nucleate and twist with half-wavelength period dominated by the interference grating in propagation direction.

Triangular bright solitons in nonlinear optics and Bose-Einstein condensates

We demonstrate what we believe to be novel triangular bright solitons that can be supported by the nonlinear Schrödinger equation with inhomogeneous Kerr-like nonlinearity and external harmonic potential, which can be realized in nonlinear optics and Bose-Einstein condensates. The profiles of these solitons are quite different from the common Gaussian or sech envelope beams, as their tops and bottoms are similar to the triangle and inverted triangle functions, respectively. The self-defocusing nonlinearity gives rise to the triangle-up solitons, while the self-focusing nonlinearity supports the triangle-down solitons. Here, we restrict our attention only to the lowest-order fundamental triangular solitons. All such solitons are stable, which is demonstrated by the linear stability analysis and also clarified by direct numerical simulations. In addition, the modulated propagation of both types of triangular solitons, with the modulated parameter being the strength of nonlinearity, is also presented. We find that such propagation is strongly affected by the form of the modulation of the nonlinearity. For example, the sudden change of the modulated parameter causes instabilities in the solitons, whereas the gradual variation generates stable solitons. Also, a periodic variation of the parameter causes the regular oscillation of solitons, with the same period. Interestingly, the triangle-up and triangle-down solitons can change into each other, when the parameter changes the sign.

Dark soliton families in quintic nonlinear lattices

We prove that the dark solitons can be stable in the purely quintic nonlinear lattices, including the fundamental, tripole and five-pole solitons. These dark soliton families are generated on the periodic nonlinear backgrounds. The propagation constant affects the forms of these solitons, while the number of poles does not lead to the variation of the backgrounds. The dark solitons are stable only when the propagation constant is moderately large.

Singular solitons

We demonstrate that the commonly known concept which treats solitons as nonsingular solutions produced by the interplay of nonlinear self-attraction and linear dispersion may be extended to include modes with a relatively weak singularity at the central point, which keeps their integral norm convergent. Such states are generated by self-repulsion, which should be strong enough, represented by septimal, quintic, and usual cubic terms in the framework of the one-, two-, and three-dimensional (1D, 2D, and 3D) nonlinear Schrödinger equations (NLSEs), respectively. Although such solutions seem counterintuitive, we demonstrate that they admit a straightforward interpretation as a result of screening of an additionally introduced attractive δ-functional potential by the defocusing nonlinearity. The strength ("bare charge") of the attractive potential is infinite in 1D, finite in 2D, and vanishingly small in 3D. Analytical asymptotics of the singular solitons at small and large distances are found, entire shapes of the solitons being produced in a numerical form. Complete stability of the singular modes is accurately predicted by the anti-Vakhitov-Kolokolov criterion (under the assumption that it applies to the model), as verified by means of numerical methods. In 2D, the NLSE with a quintic self-focusing term admits singular-soliton solutions with intrinsic vorticity too, but they are fully unstable. We also mention that dissipative singular solitons can be produced by the model with a complex coefficient in front of the nonlinear term.

Semidiscrete Quantum Droplets and Vortices

We consider a binary bosonic condensate with weak mean-field (MF) residual repulsion, loaded in an array of nearly one-dimensional traps coupled by transverse hopping. With the MF force balanced by the effectively one-dimensional attraction, induced in each trap by the Lee-Hung-Yang correction (produced by quantum fluctuations around the MF state), stable on-site- and intersite-centered semidiscrete quantum droplets (QDs) emerge in the array, as fundamental ones and self-trapped vortices, with winding numbers, at least, up to five, in both tightly bound and quasicontinuum forms. The application of a relatively strong trapping potential leads to squeezing transitions, which increase the number of sites in fundamental QDs and eventually replace vortex modes by fundamental or dipole ones. The results provide the first realization of stable semidiscrete vortex QDs, including ones with multiple vorticity.

Purely Kerr nonlinear model admitting flat-top solitons

We elaborate one- and two-dimensional (1D and 2D) models of media with self-repulsive cubic nonlinearity, whose local strength is subject to spatial modulation that admits the existence of flat-top solitons of various types, including fundamental ones, 1D multipoles, and 2D vortices. Previously, solitons of this type were only produced by models with competing nonlinearities. The present setting may be implemented in optics and Bose-Einstein condensates. The 1D version gives rise to an exact analytical solution for stable flat-top solitons, and generic families may be predicted by means of the Thomas-Fermi approximation. Stability of the obtained flat-top solitons is analyzed by means of the linear-stability analysis and direct simulations. Fundamental solitons and 1D multipoles with k=1 and 2 nodes, as well as vortices with winding number m=1, are completely stable. For multipoles with k≥3 and vortices with m≥2, alternating stripes of stability and instability are identified in their parameter spaces.

Numerical realization of the variational method for generating self-trapped beams

We introduce a numerical variational method based on the Rayleigh-Ritz optimization principle for predicting two-dimensional self-trapped beams in nonlinear media. This technique overcomes the limitation of the traditional variational approximation in performing analytical Lagrangian integration and differentiation. Approximate soliton solutions of a generalized nonlinear Schrödinger equation are obtained, demonstrating robustness of the beams of various types (fundamental, vortices, multipoles, azimuthons) in the course of their propagation. The algorithm offers possibilities to produce more sophisticated soliton profiles in general nonlinear models.

Topological Vortex and Knotted Dissipative Optical 3D Solitons Generated by 2D Vortex Solitons

We predict a new class of three-dimensional (3D) topological dissipative optical one-component solitons in homogeneous laser media with fast saturable absorption. Their skeletons formed by vortex lines where the field vanishes are tangles, i.e., N_{c} knotted or unknotted, linked or unlinked closed lines and M unclosed lines that thread all the closed lines and end at the infinitely far soliton periphery. They are generated by embedding two-dimensional laser solitons or their complexes in 3D space after their rotation around an unclosed, infinite vortex line with topological charge M_{0} (N_{c}, M, and M_{0} are integers). With such structure propagation, the "hula-hoop" solitons form; their stability is confirmed numerically. For the solitons found, all vortex lines have unit topological charge: the number of closed lines N_{c}=1 and 2 (unknots, trefoils, and Solomon knots links); unclosed vortex lines are unknotted and unlinked, their number M=1, 2, and 3.

Localized dark solitons and vortices in defocusing media with spatially inhomogeneous nonlinearity

Recent studies have demonstrated that defocusing cubic nonlinearity with local strength growing from the center to the periphery faster than r^{D}, in space of dimension D with radial coordinate r, supports a vast variety of robust bright solitons. In the framework of the same model, but with a weaker spatial-growth rate ∼r^{α} with α≤D, we test here the possibility to create stable localized continuous waves (LCWs) in one-dimensional (1D) and 2D geometries, localized dark solitons (LDSs) in one dimension, and localized dark vortices (LDVs) in two dimensions, which are all realized as loosely confined states with a divergent norm. Asymptotic tails of the solutions, which determine the divergence of the norm, are constructed in a universal analytical form by means of the Thomas-Fermi approximation (TFA). Global approximations for the LCWs, LDSs, and LDVs are constructed on the basis of interpolations between analytical approximations available far from (TFA) and close to the center. In particular, the interpolations for the 1D LDS, as well as for the 2D LDVs, are based on a deformed-tanh expression, which is suggested by the usual 1D dark-soliton solution. The analytical interpolations produce very accurate results, in comparison with numerical findings, for the 1D and 2D LCWs, 1D LDSs, and 2D LDVs with vorticity S=1. In addition to the 1D fundamental LDSs with the single notch and 2D vortices with S=1, higher-order LDSs with multiple notches are found too, as well as double LDVs, with S=2. Stability regions for the modes under consideration are identified by means of systematic simulations, the LCWs being completely stable in one and two dimensions, as they are ground states in the corresponding settings. Basic evolution scenarios are identified for those vortices that are unstable. The settings considered in this work may be implemented in nonlinear optics and in Bose-Einstein condensates.

Rotating vortex clusters in media with inhomogeneous defocusing nonlinearity

We show that media with inhomogeneous defocusing cubic nonlinearity growing toward the periphery can support a variety of stable vortex clusters nested in a common localized envelope. Nonrotating symmetric clusters are built from an even number of vortices with opposite topological charges, located at equal distances from the origin. Rotation makes the clusters strongly asymmetric, as the centrifugal force shifts some vortices to the periphery, while others approach the origin, depending on the topological charge. We obtain such asymmetric clusters as stationary states in the rotating coordinate frame, identify their existence domains, and show that the rotation may stabilize some of them.

Rotating and Precessing Dissipative-Optical-Topological-3D Solitons

We predict and study a new type of three-dimensional soliton: asymmetric rotating and precessing stable topological-dissipative-optical localized structures in homogeneous media with saturable amplification and absorption. The crucial factor determining their dynamics is the ratio of the diffusion coefficients characterizing the frequency dispersion and angular selectivity (dichroism) of the scheme. These vortex solitons exist and are stable for overcritical values of the selectivity coefficients and can be realized in lasers of large sizes with saturable absorption.

Dynamics of dipoles and vortices in nonlinearly coupled three-dimensional field oscillators

The dynamics of a pair of harmonic oscillators represented by three-dimensional fields coupled with a repulsive cubic nonlinearity is investigated through direct simulations of the respective field equations and with the help of the finite-mode Galerkin approximation (GA), which represents the two interacting fields by a superposition of 3+3 harmonic-oscillator p-wave eigenfunctions with orbital and magnetic quantum numbers l=1 and m=1, 0, -1. The system can be implemented in binary Bose-Einstein condensates, demonstrating the potential of the atomic condensates to emulate various complex modes predicted by classical field theories. First, the GA very accurately predicts a broadly degenerate set of the system's ground states in the p-wave manifold, in the form of complexes built of a dipole coaxial with another dipole or vortex, as well as complexes built of mutually orthogonal dipoles. Next, pairs of noncoaxial vortices and/or dipoles, including pairs of mutually perpendicular vortices, develop remarkably stable dynamical regimes, which feature periodic exchange of the angular momentum and periodic switching between dipoles and vortices. For a moderately strong nonlinearity, simulations of the coupled-field equations agree very well with results produced by the GA, demonstrating that the dynamics is accurately spanned by the set of six modes limited to l=1.

Precession and nutation dynamics of nonlinearly coupled non-coaxial three-dimensional matter wave vortices

Nonlinearity is the driving force for numerous important effects in nature typically showing transitions between different regimes, regular, chaotic or catastrophic behavior. Localized nonlinear modes have been the focus of intense research in areas such as fluid and gas dynamics, photonics, atomic and solid state physics etc. Due to the richness of the behavior of nonlinear systems and due to the severe numerical demands of accurate three-dimensional (3D) numerical simulations presently only little knowledge is available on the dynamics of complex nonlinear modes in 3D. Here, we investigate the dynamics of 3D non-coaxial matter wave vortices that are trapped in a parabolic potential and interact via a repulsive nonlinearity. Our numerical simulations demonstrate the existence of an unexpected and fascinating nonlinear regime that starts immediately when the nonlinearity is switched-on and is characterized by a smooth dynamics representing torque-free precession with nutations. The reported motion is proven to be robust regarding various effects such as the number of particles, dissipation and trap deformations and thus should be observable in suitably designed experiments. Since our theoretical approach, i.e., coupled nonlinear Schrödinger equations, is quite generic, we expect that the obtained novel dynamical behavior should also exist in other nonlinear systems.

Stable Solitons in Three Dimensional Free Space without the Ground State: Self-Trapped Bose-Einstein Condensates with Spin-Orbit Coupling

By means of variational methods and systematic numerical analysis, we demonstrate the existence of metastable solitons in three dimensional (3D) free space, in the context of binary atomic condensates combining contact self-attraction and spin-orbit coupling, which can be engineered by available experimental techniques. Depending on the relative strength of the intra- and intercomponent attraction, the stable solitons feature a semivortex or mixed-mode structure. In spite of the fact that the local cubic self-attraction gives rise to the supercritical collapse in 3D, and hence the setting produces no true ground state, the solitons are stable against small perturbations, motion, and collisions.

Solitons and vortices in two-dimensional discrete nonlinear Schrödinger systems with spatially modulated nonlinearity

We consider a two-dimensional (2D) generalization of a recently proposed model [Gligorić et al., Phys. Rev. E 88, 032905 (2013)], which gives rise to bright discrete solitons supported by the defocusing nonlinearity whose local strength grows from the center to the periphery. We explore the 2D model starting from the anticontinuum (AC) limit of vanishing coupling. In this limit, we can construct a wide variety of solutions including not only single-site excitations, but also dipole and quadrupole ones. Additionally, two separate families of solutions are explored: the usual "extended" unstaggered bright solitons, in which all sites are excited in the AC limit, with the same sign across the lattice (they represent the most robust states supported by the lattice, their 1D counterparts being those considered as 1D bright solitons in the above-mentioned work), and the vortex cross, which is specific to the 2D setting. For all the existing states, we explore their stability (also analytically, when possible). Typical scenarios of instability development are exhibited through direct simulations.