
We study a twodimensional incoherently pumped excitonpolariton condensate
described by an opendissipative GrossPitaevskii equation for the polariton
dynamics coupled to a rate equation for the exciton density. Adopting a
hydrodynamic approach, we use multiscale expansion methods to derive several
models appearing in the context of shallow water waves with viscosity. In
particular, we derive a Boussinesq/BenneyLuke type equation and its farfield
expansion in terms of KadomtsevPetviashviliI (KPI) equations for right and
leftgoing waves. From the KPI model, we predict the existence of
vorticityfree, weakly (algebraically) localized twodimensional darklump
solitons. We find that, in the presence of dissipation, dark lumps exhibit a
lifetime three times larger than that of planar dark solitons. Direct numerical
simulations show that dark lumps do exist, and their dissipative dynamics is
well captured by our analytical approximation. It is also shown that lumplike
and vortexlike structures can spontaneously be formed as a result of the
transverse "snaking" instability of dark soliton stripes.

In the present work, we develop an adiabatic invariant approach for the
evolution of quasionedimensional (stripe) solitons embedded in a
twodimensional BoseEinstein condensate. The results of the theory are
obtained both for the onecomponent case of dark soliton stripes, as well as
for the considerably more involved case of the twocomponent darkbright (alias
"filled dark") soliton stripes. In both cases, analytical predictions regarding
the stability and dynamics of these structures are obtained. One of our main
findings is the determination of the instability modes of the waves as a
function of the parameters of the system (such as the trap strength and the
chemical potential). Our analytical predictions are favorably compared with
results of direct numerical simulations.

The statics, stability and dynamical properties of darkbright soliton pairs
are investigated motivated by applications in a homogeneous system of
twocomponent repulsively interacting BoseEinstein condensate. One of the
intraspecies interaction coefficients is used as the relevant parameter
controlling the deviation from the integrable Manakov limit. Two different
families of stationary states are identified consisting of darkbright solitons
that are either antisymmetric (outofphase) or asymmetric (mass imbalanced)
with respect to their bright soliton. Both of the above darkbright
configurations coexist at the integrable limit of equal intra and
interspecies repulsions and are degenerate in that limit. However, they are
found to bifurcate from it in a transcritical bifurcation. The latter
interchanges the stability properties of the bound darkbright pairs rendering
the antisymmetric states unstable and the asymmetric ones stable past the
associated critical point (and vice versa before it). Finally, on the dynamical
side, it is found that large kinetic energies and thus rapid soliton collisions
are essentially unaffected by the intraspecies variation, while cases
involving near equilibrium states or breathing dynamics are significantly
modified under such a variation.

We explore the quench dynamics of a binary BoseEinstein condensate crossing
the miscibilityimmiscibility threshold and vice versa, both within and in
particular beyond the meanfield approximation. Increasing the interspecies
repulsion leads to the filamentation of the density of each species, involving
shorter wavenumbers and longer spatial scales in the manybody approach. These
filaments appear to be strongly correlated and exhibit domainwall structures.
Following the reverse quench process multiple darkantidark solitary waves are
spontaneously generated and subsequently found to decay in the manybody
scenario. We simulate singleshot images to connect our findings to possible
experimental realizations. Finally, the growth rate of the variance of a sample
of singleshots probes the degree of entanglement inherent in the system.

In the present work we consider a model that has been proposed at the
continuum level for selfdefocusing nonlinearities in atomic BECs in order to
capture phenomenologically the loss of condensate atoms to thermal ones. We
explore the model at the discrete level, using this prototypical setting
combining dispersion, nonlinearity and gain/loss to illustrate the idea that
modes associated with negative "energy" (mathematically: negative Krein
signature) can give rise to instability of excited states when gain/loss terms
are introduced in a nonlinear dynamical lattice. We showcase this idea by
considering one, two and threesite discrete modes, exploring their stability
semianalytically, and corroborating their continuation over the gain/loss
parameter numerically, as well as manifesting through direct numerical
simulations their unstable nonlinear dynamics.

For quantum fluids, the role of quantum fluctuations may be significant in
several regimes such as when the dimensionality is low, the density is high,
the interactions are strong, or for low particle numbers. In this paper we
propose a fundamentally different regime for enhanced quantum fluctuations
without being restricted by any of the above conditions. Instead, our scheme
relies on the engineering of an effective attractive interaction in a dilute,
twocomponent BoseEinstein condensate (BEC) consisting of thousands of atoms.
In such a regime, the quantum spin fluctuations are significantly enhanced
(atom bunching with respect to the noninteracting limit) since they act to
reduce the interaction energy  a remarkable property given that spin
fluctuations are normally suppressed (antibunching) at zero temperature. In
contrast to the case of true attractive interactions, our approach is not
vulnerable to BEC collapse. We numerically demonstrate that these quantum
fluctuations are experimentally accessible by either spin or singlecomponent
Bragg spectroscopy, offering a useful platform on which to test
beyondmeanfield theories. We also develop a variational model and use it to
analytically predict the shift of the immiscibility critical point, finding
good agreement with our numerics.

We study quasiperiodicityinduced localization of waves in strongly
precompressed granular chains. We propose three different setups, inspired by
the AubryAndr\'e (AA) model, of quasiperiodic chains; and we use these models
to compare the effects of onsite and offsite quasiperiodicity in nonlinear
lattices. When there is purely onsite quasiperiodicity, which we implement in
two different ways, we show for a chain of spherical particles that there is a
localization transition (as in the original AA model). However, we observe no
localization transition in a chain of cylindrical particles in which we
incorporate quasiperiodicity in the distribution of contact angles between
adjacent cylinders by making the angle periodicity incommensurate with that of
the chain. For each of our three models, we compute the Hofstadter spectrum and
the associated MinkowskiBouligand fractal dimension, and we demonstrate that
the fractal dimension decreases as one approaches the localization transition
(when it exists). Finally, in a suite of numerical computations, we demonstrate
both localization and that there exist regimes of ballistic, superdiffusive,
diffusive, and subdiffusive transport. Our models provide a flexible set of
systems to study quasiperiodicityinduced analogs of Anderson phenomena in
granular chains that one can tune controllably from weakly to strongly
nonlinear regimes.

We study  experimentally, theoretically, and numerically  nonlinear
excitations in lattices of magnets with longrange interactions. We examine
breather solutions, which are spatially localized and periodic in time, in a
chain with algebraicallydecaying interactions. It was established two decades
ago [S. Flach, Phys. Rev. E 58, R4116 (1998)] that lattices with longrange
interactions can have breather solutions in which the spatial decay of the
tails has a crossover from exponential to algebraic decay. In this Letter, we
revisit this problem in the setting of a chain of repelling magnets with a mass
defect and verify, both numerically and experimentally, the existence of
breathers with such a crossover.

The subject of PTsymmetry and its areas of application have been blossoming
over the past decade. Here, we consider a nonlinear Schr\"odinger model with a
complex potential that can be tuned controllably away from being PTsymmetric,
as it might be the case in realistic applications. We utilize two parameters:
the first one breaks PTsymmetry but retains a proportionality between the
imaginary and the derivative of the real part of the potential; the second one,
detunes from this latter proportionality. It is shown that the departure of the
potential from the PT symmetric form does not allow for the numerical
identification of exact stationary solutions. Nevertheless, it is of crucial
importance to consider the dynamical evolution of initial beam profiles. In
that light, we define a suitable notion of optimization and find that even for
non PTsymmetric cases, the beam dynamics, both in 1D and 2D although prone to
weak growth or decay suggests that the optimized profiles do not change
significantly under propagation for specific parameter regimes.

In the present work, we examine the potential robustness of extreme wave
events associated with large amplitude fluctuations of the Peregrine soliton
type, upon departure from the integrable analogue of the discrete nonlinear
Schr\"odinger (DNLS) equation, namely the AblowitzLadik (AL) model. Our model
of choice will be the socalled Salerno model, which interpolates between the
AL and the DNLS models. We find that rogue wave events essentially are
drastically distorted even for very slight perturbations of the homotopic
parameter connecting the two models off of the integrable limit. Our results
suggest that the Peregrine soliton structure is a rather sensitive feature of
the integrable limit, which may not persist under "generic" perturbations of
the limiting integrable case.

Several methods exist for finding ground (as well as excited) states of
nonlinear waves equations. In this paper we first introduce two modifications
of the socalled accelerated imaginarytime evolution method (AITEM). In our
first modification, time integration of the underlying gradient flow is done
using exponential time differencing instead of using more standard methods. In
the second modification, we present a generalization of the gradient flow
model, motivated by the work of Nesterov, as well as that of Candes and
collaborators. Additionally, we consider combinations of these methods with the
socalled spectral renormalization scheme. Finally, we apply these techniques
to the socalled Squared Operator Method, enabling convergence to excited
states. Various examples are shown to illustrate the effectiveness of these new
schemes, comparing them to standard ones established in the literature. In most
cases, we find significant reductions in the number of iterations needed to
reach convergence.

In the present work, we explore the potential of spinorbit (SO) coupled
BoseEinstein condensates to support multicomponent solitonic states in the
form of darkbright (DB) solitons. In the case where Raman linear coupling
between components is absent, we use a multiscale expansion method to reduce
the model to the integrable Mel'nikov system. The soliton solutions of the
latter allow us to reconstruct approximate traveling DB solitons for the
reduced SO coupled system. For small values of the formal perturbation
parameter, the resulting waveforms propagate undistorted, while for large
values thereof, they shed some dispersive radiation, and subsequently distill
into a robust propagating structure. After quantifying the relevant radiation
effect, we also study the dynamics of DB solitons in a parabolic trap,
exploring how their oscillation frequency varies as a function of the bright
component mass and the Raman laser wavenumber.

The beyond meanfield dynamics of a bent dark soliton embedded in a
twodimensional repulsively interacting BoseEinstein condensate is explored.
We examine the case of a single bent dark soliton comparing the meanfield
dynamics to a correlated approach, the MultiConfiguration TimeDependent
Hartree method for Bosons. Dynamical snaking of this bent structure is
observed, signaling the onset of fragmentation which becomes significant during
the vortex nucleation. In contrast to the meanfield approximation "filling" of
the vortex core is observed, leading in turn to the formation of filledcore
vortices, instead of the meanfield vortexantivortex pairs. The resulting
smearing effect in the density is a rather generic feature, occurring when
solitonic structures are exposed to quantum fluctuations. Here, we show that
this filling owes its existence to the dynamical building of an antidark
structure developed in the nexttoleading order orbital. We further
demonstrate that the aforementioned beyond meanfield dynamics can be
experimentally detected using the variance of single shot measurements.
Additionally, a variety of excitations including vortices, oblique dark
solitons, and open ring dark solitonlike structures building upon higherlying
orbitals is observed. We demonstrate that signatures of the higherlying
orbital excitations emerge in the total density, and can be clearly captured by
inspecting the onebody coherence. In the latter context, the localization of
onebody correlations exposes the existence of the multiorbital
vortexantidark structure.

We consider the motion of individual twodimensional vortices in general
radially symmetric potentials in BoseEinstein condensates. We find that
although in the special case of the parabolic trap there is a logarithmic
correction in the dependence of the precession frequency $\omega$ on the
chemical potential $\mu$, this is no longer true for a general potential $V(r)
\propto r^p$. Our calculations suggest that for $p>2$, the precession frequency
scales with $\mu$ as $\omega \sim \mu^{2/p}$. This theoretical prediction is
corroborated by numerical computations, both at the level of spectral
(Bogolyubovde Gennes) stability analysis by identifying the relevant
precession mode dependence on $\mu$, but also through direct numerical
computations of the vortex evolution in the large $\mu$, socalled
ThomasFermi, limit. Additionally, the dependence of the precession frequency
on the radius of an initially displaced from the center vortex is examined and
the corresponding predictions are tested against numerical results.

We study a deformation of the defocusing nonlinear Schr\"odinger (NLS)
equation, the defocusing Camassa Holm NLS, hereafter referred to as CHNLS
equation. We use asymptotic multiscale expansion methods to reduce this model
to a Boussinesqlike equation, which is then subsequently approximated by two
Kortewegde Vries (KdV) equations for left and righttraveling waves. We use
the soliton solution of the KdV equation to construct approximate solutions of
the CHNLS system. It is shown that these solutions may have the form of either
dark or antidark solitons, namely dips or humps on top of a stable
continuouswave background. We also use numerical simulations to investigate
the validity of the asymptotic solutions, study their evolution, and their
headon collisions. It is shown that smallamplitude dark and antidark solitons
undergo quasielastic collisions.

We consider the asymptotic behavior of the solutions of a nonlinear
Schr\"odinger (NLS) model incorporating linear and nonlinear gain/loss. First,
we describe analytically the dynamical regimes (depending on the gain/loss
strengths), for finitetime collapse, decay, and global existence of solutions
in the dynamics. Then, for all the above parametric regimes, we use direct
numerical simulations to study the dynamics corresponding to algebraically
decaying initial data. We identify crucial differences between the dynamics of
vanishing initial conditions, and those converging to a finite constant
background: in the former (latter) case we find strong (weak) collapse or
decay, when the gain/loss parameters are selected from the relevant regimes.
One of our main results, is that in all the above regimes, nonvanishing
initial data transition through spatiotemporal, algebraically decaying
waveforms. While the system is nonintegrable, the evolution of these waveforms
is reminiscent to the evolution of the Peregrine rogue wave of the integrable
NLS limit. The parametric range of gain and loss for which this phenomenology
persists is also touched upon.

We explore a variant of the $\phi^6$ model originally proposed in Phys. Rev.
D 12, 1606 (1975) as a prototypical, socalled, "bag" model where domain walls
play the role of quarks within hadrons. We examine the steady state of the
model, namely an apparent bound state of two kink structures. We explore its
linearization, and we find that, as a function of a parameter controlling the
curvature of the potential, an {\it effectively arbitrary} number of internal
modes may arise in the point spectrum of the linearization about the domain
wall profile. We explore some of the key characteristics of kinkantikink
collisions, such as the critical velocity and the multibounce windows, and how
they depend on the principal parameter of the model. We find that the critical
velocity exhibits a nonmonotonic dependence on the parameter controlling the
curvature of the potential. For the multibounce windows, we find that their
range and complexity decreases as the relevant parameter decreases (and as the
number of internal modes in the model increases). We use a modified collective
coordinate method [in the spirit of recent works such as Phys. Rev. D 94,
085008 (2016)] in order to capture the relevant phenomenology in a
semianalytical manner.

The study of energy transport properties in heterogeneous materials has
attracted scientific interest for more than a century, and it continues to
offer fundamental and rich questions. One of the unanswered challenges is to
extend Anderson theory for uncorrelated and fully disordered lattices in
condensedmatter systems to physical settings in which additional effects
compete with disorder. Specifically, the effect of strong nonlinearity has been
largely unexplored experimentally, partly due to the paucity of testbeds that
can combine the effect of disorder and nonlinearity in a controllable manner.
Here we present the first systematic experimental study of energy transport and
localization properties in simultaneously disordered and nonlinear granular
crystals. We demonstrate experimentally that disorder and nonlinearity 
which are known from decades of studies to individually favor energy
localization  can in some sense "cancel each other out", resulting in the
destruction of wave localization. We also report that the combined effect of
disorder and nonlinearity can enable the manipulation of energy transport speed
in granular crystals from subdiffusive to superdiffusive ranges.

The dynamics of darkbright solitons beyond the meanfield approximation is
investigated. We first examine the case of a single darkbright soliton and its
oscillations within a parabolic trap. Subsequently, we move to the setting of
collisions, comparing the meanfield approximation to that involving multiple
orbitals in both the dark and the bright component. Fragmentation is present
and significantly affects the dynamics, especially in the case of slower
solitons and in that of lower atom numbers. It is shown that the presence of
fragmentation allows for bipartite entanglement between the distinguishable
species. Most importantly the interplay between fragmentation and entanglement
leads to the splitting of each of the parent meanfield darkbright solitons,
placed offcenter within the parabolic trap, into a fast and a slow daughter
solitary wave. The latter process is in direct contrast to the predictions of
the meanfield approximation. A variety of excitations including darkbright
solitons in multiple (concurrently populated) orbitals is observed.
Darkantidark states and domainwallbright soliton complexes can also be
observed to arise spontaneously in the beyond meanfield dynamics.

The existence, stability, and dynamics of bound pairs of symbiotic matter
waves in the form of darkbright soliton pairs in twocomponent mixtures of
atomic BoseEinstein condensates is investigated. Motivated by the tunability
of the atomic interactions in recent experiments, we explore in detail the
impact that changes in the interaction strengths have on these bound pairs by
considering significant deviations from the Manakov (integrable) limit. It is
found that darkbright soliton pairs exist as stable configurations in a wide
parametric window spanning both the miscible and the immiscible regime of
interactions. Outside this parameter interval two unstable regions are
identified and are associated with a supercritical and a subcritical pitchfork
bifurcation, respectively. Dynamical manifestation of these instabilities gives
rise to a redistribution of the bright density between the dark solitons, and
also to symmetrybroken stationary states that are mass imbalanced (asymmetric)
with respect to their bright soliton counterpart. The longtime dynamics of
both the stable and the unstable balanced and imbalanced darkbright soliton
pairs is analyzed.

Motivated by the recent theoretical study of (bright) soliton diode effects
in systems with multiple scatterers, as well as by experimental investigations
of solitonimpurity interactions, we consider some prototypical case examples
of interactions of dark solitons with a pair of scatterers. In a way
fundamentally opposite to the case of bright solitons (but consonant to their
"antiparticle character"), we find that dark solitons accelerate as they pass
the first barrier and hence cannot be trapped by a second equalheight barrier.
A pair of unequal barriers may lead to reflection from the second one, however
trapping in the interbarrier region cannot occur. We also give some examples
of dynamical adjusting the barriers to trap the dark soliton in the
interbarrier region, yet we show that this can only occur over finite time
horizons, with the dark soliton always escaping eventually, contrary again to
what is potentially the case with bright solitons.

We study the dynamics of isospectral granular chains that are highly tunable
due to the nonlinear Hertz contact law interaction between the granular
particles. The system dynamics can thus be tuned easily from being linear to
strongly nonlinear by adjusting the initial compression applied to the chain.
In particular, we introduce both discrete and continuous spectral
transformation schemes to generate a family of granular chains that are
isospectral in their linear limit. Inspired by the principle of supersymmetry
in quantum systems, we also introduce a methodology to add or remove certain
eigenfrequencies, and we demonstrate numerically that the corresponding
physical system can be constructed in the setting of onedimensional granular
crystals. In the linear regime, we highlight the commonalities and differences
in the elastic wave transmission characteristics of such isospectral systems,
and emphasize that the presented mathematical framework allows one to suitably
tailor the wave transmission through a general class of granular chains.
Moreover, we show how the dynamic response of these structures deviates from
its linear limit as we introduce Hertzian nonlinearity in the chain and how
nonlinearity breaks the notion of linear isospectrality.

In this work we employ a recently proposed bifurcation analysis technique,
the deflated continuation algorithm, to compute steadystate solitary waveforms
in a onecomponent, two dimensional nonlinear Schr\"odinger equation with a
parabolic trap and repulsive interactions. Despite the fact that this system
has been studied extensively, we discover a wide variety of previously unknown
branches of solutions. We analyze the stability of the newly discovered
branches and discuss the bifurcations that relate them to known solutions both
in the near linear (Cartesian, as well as polar) and in the highly nonlinear
regimes. While deflated continuation is not guaranteed to compute the full
bifurcation diagram, this analysis is a potent demonstration that the algorithm
can discover new nonlinear states and provide insights into the energy
landscape of complex highdimensional Hamiltonian dynamical systems.

It is known that multidimensional complex potentials obeying
$\mathcal{PT}$symmetry may possess all real spectra and continuous families of
solitons. Recently it was shown that for multidimensional systems these
features can persist when the parity symmetry condition is relaxed so that the
potential is invariant under reflection in only a single spatial direction. We
examine the existence, stability and dynamical properties of localized modes
within the cubic nonlinear Schr\"odinger equation in such a scenario of
partially $\mathcal{PT}$symmetric potential.

In this work we present a systematic theoretical analysis regarding
darkbright solitons and their interactions, motivated by recent advances in
atomic twocomponent repulsively interacting BoseEinstein condensates. In
particular, we study analytically via a twosoliton ansatz adopted within a
variational formulation the interaction between two darkbright solitons in a
homogeneous environment beyond the integrable regime, by considering general
inter/intraatomic interaction coefficients. We retrieve the possibility of a
fixed point in the case where the bright solitons are out of phase. As the
intercomponent interaction is increased, we also identify an exponential
instability of the twosoliton state, associated with a subcritical pitchfork
bifurcation. The latter gives rise to an asymmetric partition of the bright
soliton mass and dynamically leads to spontaneous splitting of the bound pair.
In the case of the inphase bright solitons, we explain via parsing the
analytical approximations and monitoring the direct dynamics why no such pair
is identified, despite its prediction by the variational analysis.