• ### Non-integrable dynamics of matter-wave solitons in a density-dependent gauge theory(1709.08037)

Feb. 22, 2018 nlin.PS, cond-mat.quant-gas
We study interactions between bright matter-wave solitons which acquire chiral transport dynamics due to an optically-induced density-dependent gauge potential. Through numerical simulations, we find that the collision dynamics feature several non-integrable phenomena, from inelastic collisions including population transfer and radiation losses to short-lived bound states and soliton fission. An effective quasi-particle model for the interaction between the solitons is derived by means of a variational approximation, which demonstrates that the inelastic nature of the collision arises from a coupling of the gauge field to velocities of the solitons. In addition, we derive a set of interaction potentials which show that the influence of the gauge field appears as a short-range potential, that can give rise to both attractive and repulsive interactions.
• ### Stabilization of the Peregrine soliton and Kuznetsov-Ma breathers by means of nonlinearity and dispersion management(1711.08220)

Feb. 8, 2018 nlin.PS
We demonstrate a possibility to make rogue waves (RWs) in the form of the Peregrine soliton (PS) and Kuznetsov-Ma breathers (KMBs) effectively stable objects, with the help of properly defined dispersion or nonlinearity management applied to the continuous-wave (CW) background supporting the RWs. In particular, it is found that either management scheme, if applied along the longitudinal coordinate, making the underlying nonlinear Schr\"odinger equation (NLSE) selfdefocusing in the course of disappearance of the PS, indeed stabilizes the global solution with respect to the modulational instability of the background. In the process, additional excitations are generated, namely, dispersive shock waves and, in some cases, also a pair of slowly separating dark solitons. Further, the nonlinearity-management format, which makes the NLSE defocusing outside of a finite domain in the transverse direction, enables the stabilization of the KMBs, in the form of confined oscillating states. On the other hand, a nonlinearity-management format applied periodically along the propagation direction, creates expanding patterns featuring multiplication of KMBs through their cascading fission.
• ### A PT-symmetric dual-core system with the sine-Gordon nonlinearity and derivative coupling(1601.07321)

May 17, 2016 nlin.PS
We propose a system of sine-Gordon equations, with the $\mathcal{PT}$ symmetry represented by balanced gain and loss in them. The equations are coupled by sine-field terms and first-order derivatives. The sinusoidal coupling stems from local interaction between adjacent particles in the coupled Frenkel-Kontorova (FK) chains and related sine-lattices, while the cross-derivative coupling, which was not considered before, is induced by \emph{three-particle} interactions, provided that the particles in the parallel FK\ chains move in different directions. Nonlinear wave structures are then studied in this model. In particular, kink-kink (KK) and kink-antikink (KA) complexes are explored by means of analytical and numerical methods. It is predicted analytically and confirmed numerically that the complexes are unstable for one sign of the sinusoidal coupling, and stable for another. Stability regions are delineated in the underlying parameter space. Unstable complexes split into free kinks/antikinks that may propagate or become stationary, depending on whether they are subject to gain or loss, respectively.
• ### Solitons and Vortices in Two-dimensional Discrete Nonlinear Schrodinger Systems with Spatially Modulated Nonlinearity(1412.4481)

We consider a two-dimensional (2D) generalization of a recently proposed model [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 anti-continuum (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 what was 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 (analytically, whenever possible). Typical scenarios of instability development are exhibited through direct simulations.
• ### Domain walls and bubble-droplets in immiscible binary Bose gases(1410.1557)

The existence and stability of domain walls (DWs) and bubble-droplet (BD) states in binary mixtures of quasi-one-dimensional ultracold Bose gases with inter- and intra-species repulsive interactions is considered. Previously, DWs were studied by means of coupled systems of Gross-Pitaevskii equations (GPEs) with cubic terms, which model immiscible binary Bose-Einstein condensates (BECs). We address immiscible BECs with two- and three-body repulsive interactions, as well as binary Tonks--Girardeau (TG) gases, using systems of GPEs with cubic and quintic nonlinearities for the binary BEC, and coupled nonlinear Schr\"{o}dinger equations with quintic terms for the TG gases. Exact DW\ solutions are found for the symmetric BEC mixture, with equal intra-species scattering lengths. Stable asymmetric DWs in the BEC mixtures with dissimilar interactions in the two components, as well as of symmetric and asymmetric DWs in the binary TG gas, are found by means of numerical and approximate analytical methods. In the BEC system, DWs can be easily put in motion by phase imprinting. Combining a DW and anti-DW on a ring, we construct BD states for both the BEC and TG models. These consist of a dark soliton in one component (the "bubble"), and a bright soliton (the "droplet") in the other. In the BEC system, these composite states are mobile too.
• ### Damped solitons in an extended nonlinear Schrodinger equation with a spatial stimulated Raman scattering and decreasing dispersion(1401.4890)

Jan. 20, 2014 nlin.PS, physics.optics
Dynamics of solitons is considered in the framework of an extended nonlinear Schrodinger equation (NLSE), which is derived from a system of the Zakharov's type for the interaction between high- and low-frequency (HF and LF) waves. The resulting NLSE includes a pseudo-stimulated-Raman-scattering (pseudo-SRS) term, i.e., a spatial-domain counterpart of the SRS term, which is a known ingredient of the temporal-domain NLSE in optics. Also included is inhomogeneity of the spatial second-order dispersion (SOD) and linear losses of HF waves. It is shown that wavenumber downshift by the pseudo-SRS may be compensated by upshift provided by SOD whose local strength is an exponentially decaying function of the coordinate. An analytical soliton solution with a permanent shape is found in an approximate form, and is verified by comparison with numerical results
• ### Soliton dynamics in an extended nonlinear Schrodinger equation with a spatial counterpart of the stimulated Raman scattering(1306.4550)

June 19, 2013 nlin.PS
Dynamics of solitons is considered in the framework of the extended nonlinear Schrodinger equation (NLSE), which is derived from a system of Zakharov's type for the interaction between high- and low-frequency (HF and LF) waves, in which the LF field is subject to diffusive damping. The model may apply to the propagation of HF waves in plasmas. The resulting NLSE includes a pseudo-stimulated-Raman-scattering (PSRS) term, i.e., a spatial-domain counterpart of the SRS term which is well known as an ingredient of the temporal-domain NLSE in optics. Also included is inhomogeneity of the spatial second-order diffraction (SOD). It is shown that the wavenumber downshift of solitons, caused by the PSRS, may be compensated by an upshift provided by the SOD whose coefficient is a linear function of the coordinate. An analytical solution for solitons is obtained in an approximate form. Analytical and numerical results agree well, including the predicted balance between the PSRS and the linearly inhomogeneous SOD.
• ### Soliton interaction mediated by cascaded four wave mixing with dispersive waves(1305.5909)

May 25, 2013 nlin.PS, physics.optics
We demonstrate that trapping of dispersive waves between two optical solitons takes place when resonant scattering of the waves on the solitons leads to nearly perfect reflections. The momentum transfer from the radiation to solitons results in their mutual attraction and a subsequent collision. The spectrum of the trapped radiation can either expand or shrink in the course of the propagation, which is controlled by arranging either collision or separation of the solitons.
• ### Instabilities, solitons, and rogue waves in PT-coupled nonlinear waveguides(1304.7369)

April 27, 2013 nlin.PS, physics.optics
We considered the modulational instability of continuous-wave backgrounds, and the related generation and evolution of deterministic rogue waves in the recently introduced parity-time (PT)-symmetric system of linearly-coupled nonlinear Schr\"odinger equations, which describes a Kerr-nonlinear optical coupler with mutually balanced gain and loss in its cores. Besides the linear coupling, the overlapping cores are coupled through cross-phase-modulation term too. While the rogue waves, built according to the pattern of the Peregrine soliton, are (quite naturally) unstable, we demonstrate that the focusing cross-phase-modulation interaction results in their partial stabilization. For PT-symmetric and antisymmetric bright solitons, the stability region is found too, in an exact analytical form, and verified by means of direct simulations.
• ### Variational approximations for travelling solitons in a discrete nonlinear Schr\"{o}dinger equation(1201.1426)

Jan. 6, 2012 nlin.PS
Travelling solitary waves in the one-dimensional discrete nonlinear Schr\"{o}dinger equation (DNLSE) with saturable onsite nonlinearity are studied. A variational approximation (VA) for the solitary waves is derived in an analytical form. The stability is also studied by means of the VA, demonstrating that the solitons are stable, which is consistent with previously published results. Then, the VA is applied to predict parameters of travelling solitons with non-oscillatory tails (\textit{embedded solitons}, ESs). Two-soliton bound states are considered too. The separation distance between the solitons forming the bound state is derived by means of the VA. A numerical scheme based on the discretization of the equation in the moving coordinate frame is derived and implemented using the Newton--Raphson method. In general, a good agreement between the analytical and numerical results is obtained. In particular, we demonstrate the relevance of the analytical prediction of characteristics of the embedded solitons.
• ### Dragging two-dimensional discrete solitons by moving linear defects(1107.0637)

We study the mobility of small-amplitude solitons attached to moving defects which drag the solitons across a two-dimensional (2D) discrete nonlinear-Schr\"{o}dinger (DNLS) lattice. Findings are compared to the situation when a free small-amplitude 2D discrete soliton is kicked in the uniform lattice. In agreement with previously known results, after a period of transient motion the free soliton transforms into a localized mode pinned by the Peierls-Nabarro potential, irrespective of the initial velocity. However, the soliton attached to the moving defect can be dragged over an indefinitely long distance (including routes with abrupt turns and circular trajectories) virtually without losses, provided that the dragging velocity is smaller than a certain critical value. Collisions between solitons dragged by two defects in opposite directions are studied too. If the velocity is small enough, the collision leads to a spontaneous symmetry breaking, featuring fusion of two solitons into a single one, which remains attached to either of the two defects.
• ### Spatiotemporal vortex solitons in hexagonal arrays of waveguides(1105.5769)

May 29, 2011 nlin.PS, physics.optics
By means of a systematic numerical analysis, we demonstrate that hexagonal lattices of parallel linearly-coupled waveguides, with the intrinsic cubic self-focusing nonlinearity, give rise to three species of stable semi-discrete complexes (which are continuous in the longitudinal direction), with embedded vorticity S: triangular modes with S=1, hexagonal ones with S=2, both centered around an empty central core, and compact triangles with S=1, which do not not include the empty site. Collisions between stable triangular vortices are studied too. These waveguiding lattices can be realized in optics and BEC.
• ### Control of the symmetry breaking in double-well potentials by the resonant nonlinearity management(1102.1245)

Feb. 9, 2011 cond-mat.quant-gas
We introduce a one-dimensional model of Bose-Einstein condensates (BECs), combining the double-well potential, which is a well-known setting for the onset of spontaneous-symmetry-breaking (SSB) effects, and time-periodic modulation of the nonlinearity, which may be implemented by means of the Feshbach-resonance-management (FRM) technique. Both cases of the nonlinearity which is repulsive or attractive on the average are considered. In the former case, the main effect produced by the application of the FRM is spontaneous self-trapping of the condensate in either of the two potential wells in parameter regimes where it would remain untrapped in the absence of the management. In the weakly nonlinear regime, the frequency of intrinsic oscillations in the FRM-induced trapped state is very close to half the FRM frequency, suggesting that the effect is accounted for by a parametric resonance. In the case of the attractive nonlinearity, the FRM-induced effect is the opposite, i.e., enforced detrapping of a state which is self-trapped in its unmanaged form. In the latter case, the frequency of oscillations of the untrapped mode is close to a quarter of the driving frequency, suggesting that a higher-order parametric resonance may account for this effect.
• ### From one- to two-dimensional solitons in the Ginzburg-Landau model of lasers with frequency selective feedback(1102.1595)

Feb. 8, 2011 nlin.PS
We use the cubic complex Ginzburg-Landau equation coupled to a dissipative linear equation as a model of lasers with an external frequency-selective feedback. It is known that the feedback can stabilize the one-dimensional (1D) self-localized mode. We aim to extend the analysis to 2D stripe-shaped and vortex solitons. The radius of the vortices increases linearly with their topological charge, $m$, therefore the flat-stripe soliton may be interpreted as the vortex with $m=\infty$, while vortex solitons can be realized as stripes bent into rings. The results for the vortex solitons are applicable to a broad class of physical systems. There is a qualitative agreement between our results and those recently reported for models with saturable nonlinearity.
• ### Variational approximations in discrete nonlinear Schr\"odinger equations with next-nearest-neighbor couplings(1102.0296)

Feb. 1, 2011 nlin.PS
Solitons of a discrete nonlinear Schr\"{o}dinger equation which includes the next-nearest-neighbor interactions are studied by means of a variational approximation and numerical computations. A large family of multi-humped solutions, including those with a nontrivial phase structure which are a feature particular to the next-nearest-neighbor interaction model, are accurately predicted by the variational approximation. Bifurcations linking solutions with the trivial and nontrivial phase structures are also captured remarkably well, including a prediction of critical parameter values.
• ### Spontaneous symmetry breaking of Bose-Fermi mixtures in double-well potentials(1004.3520)

We study the spontaneous symmetry breaking (SSB) of a superfluid Bose-Fermi (BF) mixture in a double-well potential (DWP). The mixture is described by the Gross-Pitaevskii equation (GPE) for the bosons, coupled to an equation for the order parameter of the Fermi superfluid, which is derived from the respective density functional in the unitarity limit (a similar model applies to the BCS regime too). Straightforward SSB in the degenerate Fermi gas loaded into a DWP is impossible, as it requires an attractive self-interaction, while the intrinsic nonlinearity in the Fermi gas is repulsive. Nonetheless, we demonstrate that the symmetry breaking is possible in the mixture with attraction between fermions and bosons, like 40K and 87Rb. Numerical results are represented by dependencies of asymmetry parameters for both components on particle numbers of the mixture, N_F and N_B, and by phase diagrams in the (N_F,N_B) plane, which displays regions of symmetric and asymmetric ground states. The dynamical picture of the SSB, induced by a gradual transformation of the single-well potential into the DWP, is reported too. An analytical approximation is proposed for the case when GPE for the boson wave function may be treated by means of the Thomas-Fermi (TF) approximation. Under a special linear relation between N_F and N_B, the TF approximation allows us to reduce the model to a single equation for the fermionic function, which includes competing repulsive and attractive nonlinear terms. The latter one directly displays the mechanism of the generation of the effective attraction in the Fermi superfluid, mediated by the bosonic component of the mixture.
• ### Stabilization and destabilization of second-order solitons against perturbations in the nonlinear Schr\"{o}dinger equation(0910.0176)

Oct. 1, 2009 nlin.PS
We consider splitting and stabilization of second-order solitons (2-soliton breathers) in a model based on the nonlinear Schr\"{o}dinger equation (NLSE), which includes a small quintic term, and weak resonant nonlinearity management (NLM), i.e., time-periodic modulation of the cubic coefficient, at the frequency close to that of shape oscillations of the 2-soliton. The model applies to the light propagation in media with cubic-quintic optical nonlinearities and periodic alternation of linear loss and gain, and to BEC, with the self-focusing quintic term accounting for the weak deviation of the dynamics from one-dimensionality, while the NLM can be induced by means of the Feshbach resonance. We propose an explanation to the effect of the resonant splitting of the 2-soliton under the action of the NLM. Then, using systematic simulations and an analytical approach, we conclude that the weak quintic nonlinearity with the self-focusing sign stabilizes the 2-soliton, while the self-defocusing quintic nonlinearity accelerates its splitting. It is also shown that the quintic term with the self-defocusing/focusing sign makes the resonant response of the 2-soliton to the NLM essentially broader, in terms of the frequency.
• ### A two-dimensional paradigm for symmetry breaking: the nonlinear Schrodinger equation with a four-well potential(0904.0255)

April 1, 2009 nlin.PS
We study the existence and stability of localized states in the two-dimensional (2D) nonlinear Schrodinger (NLS)/Gross-Pitaevskii equation with a symmetric four-well potential. Using a fourmode approximation, we are able to trace the parametric evolution of the trapped stationary modes, starting from the corresponding linear limits, and thus derive the complete bifurcation diagram for the families of these stationary modes. The predictions based on the four-mode decomposition are found to be in good agreement with the numerical results obtained from the NLS equation. Actually, the stability properties coincide with those suggested by the corresponding discrete model in the large-amplitude limit. The dynamics of the unstable modes is explored by means of direct simulations. Finally, while we present the full analysis for the case of the focusing nonlinearity, the bifurcation diagram for the defocusing case is briefly considered too.
• ### Ultracold Bosons with 3-Body Attractive Interactions in an Optical Lattice(0902.4734)

Feb. 27, 2009 cond-mat.stat-mech
We study the effect of an optical lattice (OL) on the ground-state properties of one-dimensional ultracold bosons with three-body attraction and two-body repulsion, which are described by a cubic-quintic Gross-Pitaevskii equation with a periodic potential. Without the OL and with a vanishing two-body interaction term, soliton solutions of the Townes type are possible only at a critical value of the three-body interaction strength, at which an infinite degeneracy of the ground-state occurs; a repulsive two-body interaction makes such localized solutions unstable. We show that the OL opens a stability window around the critical point when the strength of the periodic potential is above a critical threshold. We also consider the effect of an external parabolic trap, studying how the stability of the solitons depends on matching between minima of the periodic potential and the minimum of the parabolic trap.
• ### Surface Solitons in Three Dimensions(0806.1087)

Aug. 22, 2008 nlin.PS
We study localized modes on the surface of a three-dimensional dynamical lattice. The stability of these structures on the surface is investigated and compared to that in the bulk of the lattice. Typically, the surface makes the stability region larger, an extreme example of that being the three-site "horseshoe"-shaped structure, which is always unstable in the bulk, while at the surface it is stable near the anti-continuum limit. We also examine effects of the surface on lattice vortices. For the vortex placed parallel to the surface this increased stability region feature is also observed, while the vortex cannot exist in a state normal to the surface. More sophisticated localized dynamical structures, such as five-site horseshoes and pyramids, are also considered.
• ### Dynamics of kicked matter-wave solitons in an optical lattice(0807.3432)

July 22, 2008 cond-mat.other
We investigate effects of the application of a kick to one-dimensional matter-wave solitons in a self-attractive Bose-Einstein condensate trapped in a optical lattice. The resulting soliton's dynamics is studied within the framework of the time-dependent nonpolynomial Schrodinger equation. The crossover from the pinning to quasi-free motion crucially depends on the size of the kick, strength of the self-attraction, and parameters of the optical lattice.
• ### Effects of Axial Vorticity in Elongated Mixtures of Bose-Einstein Condensates(0806.0262)

We consider a meniscus between rotating and nonrotating species in the Bose-Einstein condensate (BEC) with repulsive inter-atomic interactions, confined to a pipe-shaped trap. In this setting, we derive a system of coupled one-dimensional (1D) nonpolynomial Schrodinger equations (NPSEs) for two mean-field wave functions. Using these equations, we analyze the phase separation/mixing in the pipe with periodic axial boundary conditions, i.e. in a toroidal configuration. We find that the onset of the mixing, in the form of suction, i.e., filling the empty core in the vortical component by its nonrotating counterpart, crucially depends on the vorticity of the first component, and on the strengths of the inter-atomic interactions.
• ### Multistable Solitons in Higher-Dimensional Cubic-Quintic Nonlinear Schroedinger Lattices(0804.0497)

April 3, 2008 nlin.PS
We study the existence, stability, and mobility of fundamental discrete solitons in two- and three-dimensional nonlinear Schroedinger lattices with a combination of cubic self-focusing and quintic self-defocusing onsite nonlinearities. Several species of stationary solutions are constructed, and bifurcations linking their families are investigated using parameter continuation starting from the anti-continuum limit, and also with the help of a variational approximation. In particular, a species of hybrid solitons, intermediate between the site- and bond-centered types of the localized states (with no counterpart in the 1D model), is analyzed in 2D and 3D lattices. We also discuss the mobility of multi-dimensional discrete solitons that can be set in motion by lending them kinetic energy exceeding the appropriately crafted Peierls-Nabarro barrier; however, they eventually come to a halt, due to radiation loss.
• ### Trapping of two-component matter-wave solitons by mismatched optical lattices(0803.3297)

March 22, 2008 nlin.PS
We consider a one-dimensional model of a two-component Bose-Einstein condensate in the presence of periodic external potentials of opposite signs, acting on the two species. The interaction between the species is attractive, while intra-species interactions may be attractive too [the system of the right-bright (BB) type], or of opposite signs in the two components [the gap-bright (GB) model]. We identify the existence and stability domains for soliton complexes of the BB and GB types. The evolution of unstable solitons leads to the establishment of oscillatory states. The increase of the strength of the nonlinear attraction between the species results in symbiotic stabilization of the complexes, despite the fact that one component is centered around a local maximum of the respective periodic potential.
• ### Bright-Dark Soliton Complexes in Spinor Bose-Einstein Condensates(0705.1324)

Sept. 27, 2007 cond-mat.other
We present bright-dark vector solitons in quasi-one-dimensional spinor (F=1) Bose-Einstein condensates. Using a multiscale expansion technique, we reduce the corresponding nonintegrable system of three coupled Gross-Pitaevskii equations (GPEs) to a completely integrable Yajima-Oikawa system. In this way, we obtain approximate solutions for small-amplitude vector solitons of dark-dark-bright and bright-bright-dark types, in terms of the $m_{F}=+1,-1,0$ spinor components, respectively. By means of numerical simulations of the full GPE system, we demonstrate that these states indeed feature soliton properties, i.e., they propagate undistorted and undergo quasi-elastic collisions. It is also shown that, in the presence of a parabolic trap of strength $\omega$, the bright component(s) is (are) guided by the dark one(s), and, as a result, the small-amplitude vector soliton as a whole performs harmonic oscillations of frequency $\omega/ \sqrt{2}$ in the shallow soliton limit. We investigate numerically deviations from this prediction, as the depth of the solitons is increased, as well as when the strength of the spin-dependent interaction is modified.