Two-soliton collisions in a near-integrable lattice system

Discrete and Semi-Discrete Multidimensional Solitons and Vortices: Established Results and Novel Findings

This article presents a concise survey of basic discrete and semi-discrete nonlinear models, which produce two- and three-dimensional (2D and 3D) solitons, and a summary of the main theoretical and experimental results obtained for such solitons. The models are based on the discrete nonlinear Schrödinger (DNLS) equations and their generalizations, such as a system of discrete Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang corrections, the 2D Salerno model (SM), DNLS equations with long-range dipole-dipole and quadrupole-quadrupole interactions, a system of coupled discrete equations for the second-harmonic generation with the quadratic (χ(2)) nonlinearity, a 2D DNLS equation with a superlattice modulation opening mini-gaps, a discretized NLS equation with rotation, a DNLS coupler and its PT-symmetric version, a system of DNLS equations for the spin-orbit-coupled (SOC) binary Bose-Einstein condensate, and others. The article presents a review of the basic species of multidimensional discrete modes, including fundamental (zero-vorticity) and vortex solitons, their bound states, gap solitons populating mini-gaps, symmetric and asymmetric solitons in the conservative and PT-symmetric couplers, cuspons in the 2D SM, discrete SOC solitons of the semi-vortex and mixed-mode types, 3D discrete skyrmions, and some others.

How close are integrable and nonintegrable models: A parametric case study based on the Salerno model

In the present work we revisit the Salerno model as a prototypical system that interpolates between a well-known integrable system (the Ablowitz-Ladik lattice) and an experimentally tractable, nonintegrable one (the discrete nonlinear Schrödinger model). The question we ask is, for "generic" initial data, how close are the integrable to the nonintegrable models? Our more precise formulation of this question is, How well is the constancy of formerly conserved quantities preserved in the nonintegrable case? Upon examining this, we find that even slight deviations from integrability can be sensitively felt by measuring these formerly conserved quantities in the case of the Salerno model. However, given that the knowledge of these quantities requires a deep physical and mathematical analysis of the system, we seek a more "generic" diagnostic towards a manifestation of integrability breaking. We argue, based on our Salerno model computations, that the full spectrum of Lyapunov exponents could be a sensitive diagnostic to that effect.

Discrete vortices on spatially nonuniform two-dimensional electric networks

Two-dimensional arrays of nonlinear electric oscillators are considered theoretically where nearest neighbors are coupled by relatively small constant but nonequal capacitors. The dynamics is approximately reduced to a weakly dissipative defocusing discrete nonlinear Schrödinger equation with translationally noninvariant linear dispersive coefficients. Behavior of quantized discrete vortices in such systems is shown to depend strongly on the spatial profile of the internode coupling as well as on the ratio between time-increasing healing length and lattice spacings. In particular, vortex clusters can be stably trapped for some initial period of time by a circular barrier in the coupling profile, but then, due to gradual dissipative broadening of vortex cores, they lose stability and suddenly start to move.

On the effect of discreteness in the modulation instability for the Salerno model

A Salerno model with first-and second-neighbor couplings is derived for the nonlinear transmission lines. We revisit the problem of modulation instability in the Salerno model. We derive the expression for the modulation instability gain and use them to explore the role of discreteness. We show that discreteness has an impact on the mechanism by which wave trains of soliton type can be generated in the Salerno model. We also show that second-neighbor couplings have an effect on the signal voltage.

Non-standard Hubbard models in optical lattices: a review

Originally, the Hubbard model was derived for describing the behavior of strongly correlated electrons in solids. However, for over a decade now, variations of it have also routinely been implemented with ultracold atoms in optical lattices, allowing their study in a clean, essentially defect-free environment. Here, we review some of the vast literature on this subject, with a focus on more recent non-standard forms of the Hubbard model. After giving an introduction to standard (fermionic and bosonic) Hubbard models, we discuss briefly common models for mixtures, as well as the so-called extended Bose-Hubbard models, that include interactions between neighboring sites, next-neighbor sites, and so on. The main part of the review discusses the importance of additional terms appearing when refining the tight-binding approximation for the original physical Hamiltonian. Even when restricting the models to the lowest Bloch band is justified, the standard approach neglects the density-induced tunneling (which has the same origin as the usual on-site interaction). The importance of these contributions is discussed for both contact and dipolar interactions. For sufficiently strong interactions, the effects related to higher Bloch bands also become important even for deep optical lattices. Different approaches that aim at incorporating these effects, mainly via dressing the basis, Wannier functions with interactions, leading to effective, density-dependent Hubbard-type models, are reviewed. We discuss also examples of Hubbard-like models that explicitly involve higher p orbitals, as well as models that dynamically couple spin and orbital degrees of freedom. Finally, we review mean-field nonlinear Schrödinger models of the Salerno type that share with the non-standard Hubbard models nonlinear coupling between the adjacent sites. In that part, discrete solitons are the main subject of consideration. We conclude by listing some open problems, to be addressed in the future.

Phase controlling of collisions between solitons in the two-dimensional complex Ginzburg-Landau equation without viscosity

We present a systematic analysis of the outcome of soliton collisions upon variation of the relative phase φ of the solitons, in the two-dimensional cubic-quintic complex Ginzburg-Landau equation in the absence of viscosity. Three generic outcomes are identified: merger of the solitons into a single one, creation of an extra soliton, and quasielastic interaction. The velocities of the merger soliton and the extra soliton can be effectively controlled by φ. In addition, the range of the outcome of creating an extra soliton decreases to zero with the reduction of gain or the increasing of loss. The above features have potential applications in optical switching and logic gates based on interaction of optical solitons.

Supersolitons: solitonic excitations in atomic soliton chains

We show that, by tuning interactions in nonintegrable vector nonlinear Schrödinger equations modeling Bose-Einstein condensates and other relevant physical systems, it is possible to achieve a regime of elastic particlelike collisions between solitons. This would allow one to construct a Newton's cradle with solitons and supersolitons: localized collective excitations in solitary-wave chains.

Radiationless energy exchange in three-soliton collisions

We revisit the problem of the three-soliton collisions in the weakly perturbed sine-Gordon equation and develop an effective three-particle model allowing us to explain many interesting features observed in numerical simulations of the soliton collisions. In particular, we explain why collisions between two kinks and one antikink are observed to be practically elastic or strongly inelastic depending on relative initial positions of the kinks. The fact that the three-soliton collisions become more elastic with an increase in the collision velocity also becomes clear in the framework of the three-particle model. The three-particle model does not involve internal modes of the kinks, but it gives a qualitative description to all the effects observed in the three-soliton collisions, including the fractal scattering and the existence of short-lived three-soliton bound states. The radiationless energy exchange between the colliding solitons in weakly perturbed integrable systems takes place in the vicinity of the separatrix multi-soliton solutions of the corresponding integrable equations, where even small perturbations can result in a considerable change in the collision outcome. This conclusion is illustrated through the use of the reduced three-particle model.

Moving solitons in the discrete nonlinear Schrödinger equation

Using the method of asymptotics beyond all orders, we evaluate the amplitude of radiation from a moving small-amplitude soliton in the discrete nonlinear Schrödinger equation. When the nonlinearity is of the cubic type, this amplitude is shown to be nonzero for all velocities and therefore small-amplitude solitons moving without emitting radiation do not exist. In the case of a saturable nonlinearity, on the other hand, the radiation is found to be completely suppressed when the soliton moves at one of certain isolated "sliding velocities." We show that a discrete soliton moving at a general speed will experience radiative deceleration until it either stops and remains pinned to the lattice or--in the saturable case--locks, metastably, onto one of the sliding velocities. When the soliton's amplitude is small, however, this deceleration is extremely slow; hence, despite losing energy to radiation, the discrete soliton may spend an exponentially long time traveling with virtually unchanged amplitude and speed.

Discrete solitons and vortices in the two-dimensional Salerno model with competing nonlinearities

An anisotropic lattice model in two spatial dimensions, with on-site and intersite cubic nonlinearities (the Salerno model), is introduced, with emphasis on the case in which the intersite nonlinearity is self-defocusing, competing with on-site self-focusing. The model applies, for example, to a dipolar Bose-Einstein condensate trapped in a deep two-dimensional (2D) optical lattice. Soliton families of two kinds are found in the model: ordinary ones and cuspons, with peakons at the border between them. Stability borders for the ordinary solitons are found, while all cuspons (and peakons) are stable. The Vakhitov-Kolokolov criterion does not apply to cuspons, but for the ordinary solitons it correctly identifies the stability limits. In direct simulations, unstable solitons evolve into localized pulsons. Varying the anisotropy parameter, we trace a transition between the solitons in 1D and 2D versions of the model. In the isotropic model, we also construct discrete vortices of two types, on-site and intersite centered (vortex crosses and squares, respectively), and identify their stability regions. In simulations, unstable vortices in the noncompeting model transform into regular solitons, while in the model with the competing nonlinearities they evolve into localized vortical pulsons, which maintain their topological character. Bound states of regular solitons and vortices are constructed too, and their stability is identified.

Solitons in the Salerno model with competing nonlinearities

We consider a lattice equation (Salerno model) combining onsite self-focusing and intersite self-defocusing cubic terms, which may describe a Bose-Einstein condensate of dipolar atoms trapped in a strong periodic potential. In the continuum approximation, the model gives rise to solitons in a finite band of frequencies, with sechlike solitons near one edge, and an exact peakon solution at the other. A similar family of solitons is found in the discrete system, including a peakon; beyond the peakon, the family continues in the form of cuspons. Stability of the lattice solitons is explored through computation of eigenvalues for small perturbations, and by direct simulations. A small part of the family is unstable (in that case, the discrete solitons transform into robust pulsonic excitations); both peakons and cuspons are stable. The Vakhitov-Kolokolov criterion precisely explains the stability of regular solitons and peakons, but does not apply to cuspons. In-phase and out-of-phase bound states of solitons are also constructed. They exchange their stability at a point where the bound solitons are peakons. Mobile solitons, composed of a moving core and background, exist up to a critical value of the strength of the self-defocusing intersite nonlinearity. Colliding solitons always merge into a single pulse.

Nonintegrable Schrodinger discrete breathers

In an extensive numerical investigation of nonintegrable translational motion of discrete breathers in nonlinear Schrödinger lattices, we have used a regularized Newton algorithm to continue these solutions from the limit of the integrable Ablowitz-Ladik lattice. These solutions are shown to be a superposition of a localized moving core and an excited extended state (background) to which the localized moving pulse is spatially asymptotic. The background is a linear combination of small amplitude nonlinear resonant plane waves and it plays an essential role in the energy balance governing the translational motion of the localized core. Perturbative collective variable theory predictions are critically analyzed in the light of the numerical results.