Interaction of solitons in tunnel-coupled optical fibers

A pair of weakly interacting optical fibers, which are described by two nonlinear Schrödinger equations with linear coupling terms, is considered in the framework of perturbation theory. It is shown that the interaction of optical solitons in this system is inelastic and that they may form a bound state (bisoliton). The stability and internal oscillations of the bisoliton are also investigated.

Semidiscrete composite solitons in arrays of quadratically nonlinear waveguides

We demonstrate that an array of discrete waveguides on a slab substrate, both featuring chi2 nonlinearity, supports stable solitons composed of discrete and continuous components. Two classes of fundamental composite soliton are identified: ones consisting of a discrete fundamental-frequency (FF) component in the waveguide array, coupled to a continuous second-harmonic (SH) component in the slab waveguide, and solitons with an inverted FF/SH structure. Twisted bound states of the fundamental solitons are found, too. In contrast with the usual systems, the intersite-centered fundamental solitons and bound states with the twisted continuous components are stable over almost the entire domain of their existence.

Interactions of solitons with positive and negative masses: Shuttle motion and coacceleration

We consider a possibility to realize self-accelerating motion of interacting states with effective positive and negative masses in the form of pairs of solitons in two-component BEC loaded in an optical-lattice (OL) potential. A crucial role is played by the fact that gap solitons may feature a negative dynamical mass, keeping their mobility in the OL. First, the respective system of coupled Gross-Pitaevskii equations (GPE) is reduced to a system of equations for envelopes of the lattice wave functions. Two generic dynamical regimes are revealed by simulations of the reduced system, viz., shuttle oscillations of pairs of solitons with positive and negative masses, and splitting of the pair. The coaccelerating motion of the interacting solitons, which keeps constant separation between them, occurs at the boundary between the shuttle motion and splitting. The position of the coacceleration regime in the system's parameter space can be adjusted with the help of an additional gravity potential, which induces its own acceleration, that may offset the relative acceleration of the two solitons, while gravity masses of both solitons remain positive. The numerical findings are accurately reproduced by a variational approximation. Collisions between shuttling or coaccelerating soliton pairs do not alter the character of the dynamical regime. Finally, regimes of the shuttle motion, coacceleration, and splitting are corroborated by simulations of the original GPE system, with the explicitly present OL potential.

Spatiotemporal dissipative solitons and vortices in a multi-transverse-mode fiber laser

We introduce a model for spatiotemporal modelocking in multimode fiber lasers, which is based on the (3+1)-dimensional cubic-quintic complex Ginzburg-Landau equation (cGLE) with conservative and dissipative nonlinearities and a 2-dimensional transverse trapping potential. Systematic numerical analysis reveals a variety of stable nonlinear modes, including stable fundamental solitons and breathers, as well as solitary vortices with winding number n = 1, while vortices with n = 2 are unstable, splitting into persistently rotating bound states of two unitary vortices. A characteristic feature of the system is bistability between the fundamental and vortex spatiotemporal solitons.

Metastable soliton necklaces supported by fractional diffraction and competing nonlinearities

We demonstrate that the fractional cubic-quintic nonlinear Schrödinger equation, characterized by its Lévy index, maintains ring-shaped soliton clusters ("necklaces") carrying orbital angular momentum. They can be built, in the respective optical setting, as circular chains of fundamental solitons linked by a vortical phase field. We predict semi-analytically that the metastable necklace-shaped clusters persist, corresponding to a local minimum of an effective potential of interaction between adjacent solitons in the cluster. Systematic simulations corroborate that the clusters stay robust over extremely large propagation distances, even in the presence of strong random perturbations.

Discrete solitons and vortices on anisotropic lattices

We consider the effects of anisotropy on solitons of various types in two-dimensional nonlinear lattices, using the discrete nonlinear Schrödinger equation as a paradigm model. For fundamental solitons, we develop a variational approximation that predicts that broad quasicontinuum solitons are unstable, while their strongly anisotropic counterparts are stable. By means of numerical methods, it is found that, in the general case, the fundamental solitons and simplest on-site-centered vortex solitons ("vortex crosses") feature enhanced or reduced stability areas, depending on the strength of the anisotropy. More surprising is the effect of anisotropy on the so-called "super-symmetric" intersite-centered vortices ("vortex squares"), with the topological charge equal to the square's size : we predict in an analytical form by means of the Lyapunov-Schmidt theory, and confirm by numerical results, that arbitrarily weak anisotropy results in dramatic changes in the stability and dynamics in comparison with the degenerate, in this case, isotropic, limit.

Three-dimensional spinning solitons in the cubic-quintic nonlinear medium

We find one-parameter families of three-dimensional spatiotemporal bright vortex solitons (doughnuts, or spinning light bullets), in bulk dispersive cubic-quintic optically nonlinear media. The spinning solitons display a symmetry-breaking azimuthal instability, which leads to breakup of the spinning soliton into a set of fragments, each being a stable nonspinning light bullet. However, in some cases the instability is developing so slowly that the spinning light bullets may be regarded as virtually stable ones, from the standpoint of an experiment with finite-size samples.

Dark-bright solitons in coupled nonlinear Schrödinger equations with unequal dispersion coefficients

We study a two-component nonlinear Schrödinger system with equal, repulsive cubic interactions and different dispersion coefficients in the two components. We consider states that have a dark solitary wave in one component. Treating it as a frozen one, we explore the possibility of the formation of bright-solitonic structures in the other component. We identify bifurcation points at which such states emerge in the bright component in the linear limit and explore their continuation into the nonlinear regime. An additional analytically tractable limit is found to be that of vanishing dispersion of the bright component. We numerically identify regimes of potential stability, not only of the single-peak ground state (the dark-bright soliton), but also of excited states with one or more zero crossings in the bright component. When the states are identified as unstable, direct numerical simulations are used to investigate the outcome of the instability development. Although our principal focus is on the homogeneous setting, we also briefly touch upon the counterintuitive impact of the potential presence of a parabolic trap on the states of interest.

Dissipative damping of dark solitons in optical fibers

The adiabatic evolution of the parameters of a dark soliton in the presence of dissipative fiber losses is considered. It is shown that the rate at which a dark soliton decays and spreads is only half the corresponding rate for a bright soliton.

Metastability of Quantum Droplet Clusters

We show that metastable ring-shaped clusters can be constructed from two-dimensional quantum droplets in systems described by the Gross-Pitaevskii equations augmented with Lee-Huang-Yang quantum corrections. The clusters exhibit dynamical behavior ranging from contraction to rotation with simultaneous periodic pulsations, or expansion, depending on the initial radius of the necklace pattern and phase shift between adjacent quantum droplets. We show that, using an energy-minimization analysis, one can predict equilibrium values of the cluster radius that correspond to rotation without radial pulsations. In such a regime, the clusters evolve as metastable states, withstanding abrupt variations in the underlying scattering lengths and keeping their azimuthal symmetry in the course of evolution, even in the presence of considerable perturbations.

Fermi-Pasta-Ulam phenomena and persistent breathers in the harmonic trap

We consider the long-term weakly nonlinear evolution governed by the two-dimensional nonlinear Schrödinger (NLS) equation with an isotropic harmonic oscillator potential. The dynamics in this regime is dominated by resonant interactions between quartets of linear normal modes, accurately captured by the corresponding resonant approximation. Within this approximation, we identify Fermi-Pasta-Ulam-like recurrence phenomena, whereby the normal-mode spectrum passes in close proximity of the initial configuration, and two-mode states with time-independent mode amplitude spectra that translate into long-lived breathers of the original NLS equation. We comment on possible implications of these findings for nonlinear optics and matter-wave dynamics in Bose-Einstein condensates.

Spatiotemporally localized multidimensional solitons in self-induced transparency media

"Light bullets" are multidimensional solitons which are localized in both space and time. We show that such solitons exist in two- and three-dimensional self-induced transparency media and that they are fully stable. Our approximate analytical calculation, backed and verified by direct numerical simulations, yields the multidimensional generalization of the one-dimensional sine-Gordon soliton.

Symmetric and asymmetric solitons in linearly coupled Bragg gratings

We demonstrate that a symmetric system of two linearly coupled waveguides, with Kerr nonlinearity and resonant grating in both of them, gives rise to a family of symmetric and antisymmetric solitons in an exact analytical form, a part of which exists outside of the bandgap in the system's spectrum, i.e., they may be regarded as embedded solitons (ES's, i.e., the ones partly overlapping with the continuous spectrum). Parameters of the family are the soliton's amplitude and velocity. Asymmetric ES's, unlike the regular (nonembedded) gap solitons (GS's), do not exist in the system. Moreover, ES's exist even in the case when the system's spectrum contains no bandgap. The main issue is the stability of the solitons. We demonstrate that some symmetric ES's are stable, while all the antisymmetric solitons are unstable; an explanation is given to the latter property, based on the consideration of the system's Hamiltonian. We produce a full stability diagram, which comprises both embedded and regular solitons, quiescent and moving. A stability region for ES's is found around the point where the constant of the linear coupling between the two cores is equal to the Bragg-reflectivity coefficient accounting for the linear conversion between the right- and left-traveling waves in each core, i.e., the ES's are the "most endemic" solitary solitons in this system. The stability region quickly shrinks with the increase of the soliton's velocity c, and completely disappears when c exceeds half the maximum velocity. Collisions between stable moving solitons of various types are also considered, with a conclusion that the collisions are always quasielastic.

Stabilization of dark solitons in the cubic ginzburg-landau equation

The existence and stability of exact continuous-wave and dark-soliton solutions to a system consisting of the cubic complex Ginzburg-Landau (CGL) equation linearly coupled with a linear dissipative equation is studied. We demonstrate the existence of vast regions in the system's parameter space associated with stable dark-soliton solutions, having the form of the Nozaki-Bekki envelope holes, in contrast to the case of the conventional CGL equation, where they are unstable. In the case when the dark soliton is unstable, two different types of instability are identified. The proposed stabilized model may be realized in terms of a dual-core nonlinear optical fiber, with one core active and one passive.

Robust soliton clusters in media with competing cubic and quintic nonlinearities

Systematic results are reported for dynamics of circular patterns ("necklaces"), composed of fundamental solitons and carrying orbital angular momentum, in the two-dimensional model, which describes the propagation of light beams in bulk media combining self-focusing cubic and self-defocusing quintic nonlinearities. Semianalytical predictions for the existence of quasistable necklace structures are obtained on the basis of an effective interaction potential. Then, direct simulations are run. In the case when the initial pattern is far from an equilibrium size predicted by the potential, it cannot maintain its shape. However, a necklace with the initial shape close to the predicted equilibrium survives very long evolution, featuring persistent oscillations. The quasistable evolution is not essentially disturbed by a large noise component added to the initial configuration. Basic conclusions concerning the necklace dynamics in this model are qualitatively the same as in a recently studied one which combines quadratic and self-defocusing cubic nonlinearities. Thus, we infer that a combination of competing self-focusing and self-defocusing nonlinearities enhances the robustness not only of vortex solitons but also of vorticity-carrying necklace patterns.

Raman-induced optical shocks in nonlinear fibers

It is demonstrated that the stimulated Raman scattering in optical fibers supports kink solitons (shocks) not only in the recently considered case of anomalous dispersion but also in the normal-dispersion regime. In both cases there is a family of the kinks depending on an arbitrary (positive) wave number of the carrying wave. Although all the kinks are unstable because of the inherent instability of the cw background behind them, in the normal-dispersion regime the instability, induced solely by the stimulated Raman scattering, is essentially weaker than that in the case of anomalous dispersion.

Ring dark solitons and vortex necklaces in Bose-Einstein condensates

We introduce the concept of ring dark solitons in Bose-Einstein condensates. We show that relatively shallow rings are not subject to the snake instability, but a deeper ring splits into a robust ringlike cluster of vortex pairs, which performs oscillations in the radial and azimuthal directions, following the dynamics of the original ring soliton.

Unbreakable PT symmetry of solitons supported by inhomogeneous defocusing nonlinearity

We consider bright solitons supported by a symmetric inhomogeneous defocusing nonlinearity growing rapidly enough toward the periphery of the medium, combined with an antisymmetric gain-loss profile. Despite the absence of any symmetric modulation of the linear refractive index, which is usually required to establish a parity-time (PT) symmetry in the form of a purely real spectrum of modes, we show that the PT symmetry is never broken in the present system, and that the system always supports stable bright solitons, i.e., fundamental and multi-pole ones. This fact is connected to the nonlinearizability of the underlying evolution equation. The increase of the gain-loss strength results, in lieu of the PT symmetry breaking, in merger of pairs of different soliton branches, such as fundamental and dipole, or tripole and quadrupole ones. The fundamental and dipole solitons remain stable at arbitrarily large values of the gain-loss coefficient.

Modulational instability and soliton generation in chiral Bose-Einstein condensates with zero-energy nonlinearity

By means of analytical and numerical methods, we address the modulational instability (MI) in chiral condensates governed by the Gross-Pitaevskii equation including the current nonlinearity. The analysis shows that this nonlinearity partly suppresses the MI driven by the cubic self-focusing, although the current nonlinearity is not represented in the system's energy (although it modifies the momentum), hence it may be considered as zero-energy nonlinearity. Direct simulations demonstrate generation of trains of stochastically interacting chiral solitons by MI. In the ring-shaped setup, the MI creates a single traveling solitary wave. The sign of the current nonlinearity determines the direction of propagation of the emerging solitons.

Localized modes in dissipative lattice media: an overview

We give an overview of recent theoretical studies of the dynamics of one- and two-dimensional spatial dissipative solitons in models based on the complex Ginzburg-Landau equations with the cubic-quintic combination of loss and gain terms, which include imaginary, real or complex spatially periodic potentials. The imaginary potential represents periodic modulation of the local loss and gain. It is shown that the effective gradient force, induced by the inhomogeneous loss distribution, gives rise to three generic propagation scenarios for one-dimensional dissipative solitons: transverse drift, persistent swing motion, and damped oscillations. When the lattice-average loss/gain value is zero, and the real potential has spatial parity opposite to that of the imaginary component, the respective complex potential is a realization of the parity-time symmetry. Under the action of lattice potentials of the latter type, one-dimensional solitons feature motion regimes in the form of the transverse drift and persistent swing. In the two-dimensional geometry, three types of axisymmetric radial lattices are considered, namely those based solely on the refractive-index modulation, or solely on the linear-loss modulation, or on a combination of both. The rotary motion of solitons in such axisymmetric potentials can be effectively controlled by varying the strength of the initial tangential kick.