
We study interactions between bright matterwave solitons which acquire
chiral transport dynamics due to an opticallyinduced densitydependent gauge
potential. Through numerical simulations, we find that the collision dynamics
feature several nonintegrable phenomena, from inelastic collisions including
population transfer and radiation losses to shortlived bound states and
soliton fission. An effective quasiparticle 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 shortrange potential, that can give rise to both attractive and
repulsive interactions.

We demonstrate a possibility to make rogue waves (RWs) in the form of the
Peregrine soliton (PS) and KuznetsovMa breathers (KMBs) effectively stable
objects, with the help of properly defined dispersion or nonlinearity
management applied to the continuouswave (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 nonlinearitymanagement 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 nonlinearitymanagement format applied periodically along the
propagation direction, creates expanding patterns featuring multiplication of
KMBs through their cascading fission.

We propose a system of sineGordon equations, with the $\mathcal{PT}$
symmetry represented by balanced gain and loss in them. The equations are
coupled by sinefield terms and firstorder derivatives. The sinusoidal
coupling stems from local interaction between adjacent particles in the coupled
FrenkelKontorova (FK) chains and related sinelattices, while the
crossderivative coupling, which was not considered before, is induced by
\emph{threeparticle} 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, kinkkink (KK) and kinkantikink (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.

We consider a twodimensional (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
anticontinuum (AC) limit of vanishing coupling. In this limit, we can
construct a wide variety of solutions including not only singlesite
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 abovementioned 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.

The existence and stability of domain walls (DWs) and bubbledroplet (BD)
states in binary mixtures of quasionedimensional ultracold Bose gases with
inter and intraspecies repulsive interactions is considered. Previously, DWs
were studied by means of coupled systems of GrossPitaevskii equations (GPEs)
with cubic terms, which model immiscible binary BoseEinstein condensates
(BECs). We address immiscible BECs with two and threebody repulsive
interactions, as well as binary TonksGirardeau (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 intraspecies
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 antiDW 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.

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 lowfrequency (HF and LF) waves. The
resulting NLSE includes a pseudostimulatedRamanscattering (pseudoSRS) term,
i.e., a spatialdomain counterpart of the SRS term, which is a known ingredient
of the temporaldomain NLSE in optics. Also included is inhomogeneity of the
spatial secondorder dispersion (SOD) and linear losses of HF waves. It is
shown that wavenumber downshift by the pseudoSRS 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

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 lowfrequency (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
pseudostimulatedRamanscattering (PSRS) term, i.e., a spatialdomain
counterpart of the SRS term which is well known as an ingredient of the
temporaldomain NLSE in optics. Also included is inhomogeneity of the spatial
secondorder 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.

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.

We considered the modulational instability of continuouswave backgrounds,
and the related generation and evolution of deterministic rogue waves in the
recently introduced paritytime (PT)symmetric system of linearlycoupled
nonlinear Schr\"odinger equations, which describes a Kerrnonlinear optical
coupler with mutually balanced gain and loss in its cores. Besides the linear
coupling, the overlapping cores are coupled through crossphasemodulation term
too. While the rogue waves, built according to the pattern of the Peregrine
soliton, are (quite naturally) unstable, we demonstrate that the focusing
crossphasemodulation interaction results in their partial stabilization. For
PTsymmetric and antisymmetric bright solitons, the stability region is found
too, in an exact analytical form, and verified by means of direct simulations.

Travelling solitary waves in the onedimensional 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 nonoscillatory tails (\textit{embedded solitons}, ESs).
Twosoliton 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 NewtonRaphson 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.

We study the mobility of smallamplitude solitons attached to moving defects
which drag the solitons across a twodimensional (2D) discrete
nonlinearSchr\"{o}dinger (DNLS) lattice. Findings are compared to the
situation when a free smallamplitude 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 PeierlsNabarro 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.

By means of a systematic numerical analysis, we demonstrate that hexagonal
lattices of parallel linearlycoupled waveguides, with the intrinsic cubic
selffocusing nonlinearity, give rise to three species of stable semidiscrete
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.

We introduce a onedimensional model of BoseEinstein condensates (BECs),
combining the doublewell potential, which is a wellknown setting for the
onset of spontaneoussymmetrybreaking (SSB) effects, and timeperiodic
modulation of the nonlinearity, which may be implemented by means of the
Feshbachresonancemanagement (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
selftrapping 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 FRMinduced 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 FRMinduced effect
is the opposite, i.e., enforced detrapping of a state which is selftrapped 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 higherorder parametric resonance may account for this effect.

We use the cubic complex GinzburgLandau equation coupled to a dissipative
linear equation as a model of lasers with an external frequencyselective
feedback. It is known that the feedback can stabilize the onedimensional (1D)
selflocalized mode. We aim to extend the analysis to 2D stripeshaped and
vortex solitons. The radius of the vortices increases linearly with their
topological charge, $m$, therefore the flatstripe 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.

Solitons of a discrete nonlinear Schr\"{o}dinger equation which includes the
nextnearestneighbor interactions are studied by means of a variational
approximation and numerical computations. A large family of multihumped
solutions, including those with a nontrivial phase structure which are a
feature particular to the nextnearestneighbor 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.

We study the spontaneous symmetry breaking (SSB) of a superfluid BoseFermi
(BF) mixture in a doublewell potential (DWP). The mixture is described by the
GrossPitaevskii 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 selfinteraction, 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 singlewell 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 ThomasFermi (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.

We consider splitting and stabilization of secondorder solitons (2soliton
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., timeperiodic modulation of the cubic coefficient, at the
frequency close to that of shape oscillations of the 2soliton. The model
applies to the light propagation in media with cubicquintic optical
nonlinearities and periodic alternation of linear loss and gain, and to BEC,
with the selffocusing quintic term accounting for the weak deviation of the
dynamics from onedimensionality, 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 2soliton under the action of the NLM. Then, using systematic
simulations and an analytical approach, we conclude that the weak quintic
nonlinearity with the selffocusing sign stabilizes the 2soliton, while the
selfdefocusing quintic nonlinearity accelerates its splitting. It is also
shown that the quintic term with the selfdefocusing/focusing sign makes the
resonant response of the 2soliton to the NLM essentially broader, in terms of
the frequency.

We study the existence and stability of localized states in the
twodimensional (2D) nonlinear Schrodinger (NLS)/GrossPitaevskii equation with
a symmetric fourwell 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 fourmode 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 largeamplitude 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.

We study the effect of an optical lattice (OL) on the groundstate properties
of onedimensional ultracold bosons with threebody attraction and twobody
repulsion, which are described by a cubicquintic GrossPitaevskii equation
with a periodic potential. Without the OL and with a vanishing twobody
interaction term, soliton solutions of the Townes type are possible only at a
critical value of the threebody interaction strength, at which an infinite
degeneracy of the groundstate occurs; a repulsive twobody 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.

We study localized modes on the surface of a threedimensional 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 threesite
"horseshoe"shaped structure, which is always unstable in the bulk, while at
the surface it is stable near the anticontinuum 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 fivesite horseshoes and pyramids, are
also considered.

We investigate effects of the application of a kick to onedimensional
matterwave solitons in a selfattractive BoseEinstein condensate trapped in a
optical lattice. The resulting soliton's dynamics is studied within the
framework of the timedependent nonpolynomial Schrodinger equation. The
crossover from the pinning to quasifree motion crucially depends on the size
of the kick, strength of the selfattraction, and parameters of the optical
lattice.

We consider a meniscus between rotating and nonrotating species in the
BoseEinstein condensate (BEC) with repulsive interatomic interactions,
confined to a pipeshaped trap. In this setting, we derive a system of coupled
onedimensional (1D) nonpolynomial Schrodinger equations (NPSEs) for two
meanfield 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 interatomic interactions.

We study the existence, stability, and mobility of fundamental discrete
solitons in two and threedimensional nonlinear Schroedinger lattices with a
combination of cubic selffocusing and quintic selfdefocusing onsite
nonlinearities. Several species of stationary solutions are constructed, and
bifurcations linking their families are investigated using parameter
continuation starting from the anticontinuum limit, and also with the help of
a variational approximation. In particular, a species of hybrid solitons,
intermediate between the site and bondcentered 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 multidimensional discrete solitons that can be
set in motion by lending them kinetic energy exceeding the appropriately
crafted PeierlsNabarro barrier; however, they eventually come to a halt, due
to radiation loss.

We consider a onedimensional model of a twocomponent BoseEinstein
condensate in the presence of periodic external potentials of opposite signs,
acting on the two species. The interaction between the species is attractive,
while intraspecies interactions may be attractive too [the system of the
rightbright (BB) type], or of opposite signs in the two components [the
gapbright (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.

We present brightdark vector solitons in quasionedimensional spinor (F=1)
BoseEinstein condensates. Using a multiscale expansion technique, we reduce
the corresponding nonintegrable system of three coupled GrossPitaevskii
equations (GPEs) to a completely integrable YajimaOikawa system. In this way,
we obtain approximate solutions for smallamplitude vector solitons of
darkdarkbright and brightbrightdark 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 quasielastic 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
smallamplitude 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 spindependent interaction is
modified.