• We study a two-dimensional incoherently pumped exciton-polariton condensate described by an open-dissipative Gross-Pitaevskii 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/Benney-Luke type equation and its far-field expansion in terms of Kadomtsev-Petviashvili-I (KP-I) equations for right- and left-going waves. From the KP-I model, we predict the existence of vorticity-free, weakly (algebraically) localized two-dimensional dark-lump 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 lump-like and vortex-like 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 quasi-one-dimensional (stripe) solitons embedded in a two-dimensional Bose-Einstein condensate. The results of the theory are obtained both for the one-component case of dark soliton stripes, as well as for the considerably more involved case of the two-component dark-bright (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 dark-bright soliton pairs are investigated motivated by applications in a homogeneous system of two-component repulsively interacting Bose-Einstein condensate. One of the intra-species 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 dark-bright solitons that are either antisymmetric (out-of-phase) or asymmetric (mass imbalanced) with respect to their bright soliton. Both of the above dark-bright configurations coexist at the integrable limit of equal intra- and inter-species 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 dark-bright 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 intra-species 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 Bose-Einstein condensate crossing the miscibility-immiscibility threshold and vice versa, both within and in particular beyond the mean-field approximation. Increasing the interspecies repulsion leads to the filamentation of the density of each species, involving shorter wavenumbers and longer spatial scales in the many-body approach. These filaments appear to be strongly correlated and exhibit domain-wall structures. Following the reverse quench process multiple dark-antidark solitary waves are spontaneously generated and subsequently found to decay in the many-body scenario. We simulate single-shot images to connect our findings to possible experimental realizations. Finally, the growth rate of the variance of a sample of single-shots 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 self-defocusing 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 three-site discrete modes, exploring their stability semi-analytically, 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, two-component Bose-Einstein 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 (anti-bunching) 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 single-component Bragg spectroscopy, offering a useful platform on which to test beyond-mean-field 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 quasiperiodicity-induced localization of waves in strongly precompressed granular chains. We propose three different setups, inspired by the Aubry--Andr\'e (AA) model, of quasiperiodic chains; and we use these models to compare the effects of on-site and off-site quasiperiodicity in nonlinear lattices. When there is purely on-site 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 Minkowski--Bouligand 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 quasiperiodicity-induced 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 long-range interactions. We examine breather solutions, which are spatially localized and periodic in time, in a chain with algebraically-decaying interactions. It was established two decades ago [S. Flach, Phys. Rev. E 58, R4116 (1998)] that lattices with long-range 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 PT-symmetry 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 PT-symmetric, as it might be the case in realistic applications. We utilize two parameters: the first one breaks PT-symmetry 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 PT-symmetric 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 Ablowitz-Ladik (AL) model. Our model of choice will be the so-called 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 so-called accelerated imaginary-time 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 so-called spectral renormalization scheme. Finally, we apply these techniques to the so-called 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 spin-orbit (SO) coupled Bose-Einstein condensates to support multi-component solitonic states in the form of dark-bright (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 mean-field dynamics of a bent dark soliton embedded in a two-dimensional repulsively interacting Bose-Einstein condensate is explored. We examine the case of a single bent dark soliton comparing the mean-field dynamics to a correlated approach, the Multi-Configuration Time-Dependent 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 mean-field approximation "filling" of the vortex core is observed, leading in turn to the formation of filled-core vortices, instead of the mean-field vortex-antivortex 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 next-to-leading order orbital. We further demonstrate that the aforementioned beyond mean-field 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 soliton-like structures building upon higher-lying orbitals is observed. We demonstrate that signatures of the higher-lying orbital excitations emerge in the total density, and can be clearly captured by inspecting the one-body coherence. In the latter context, the localization of one-body correlations exposes the existence of the multi-orbital vortex-antidark structure.
  • We consider the motion of individual two-dimensional vortices in general radially symmetric potentials in Bose-Einstein 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 (Bogolyubov-de 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$, so-called Thomas-Fermi, 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 CH-NLS equation. We use asymptotic multiscale expansion methods to reduce this model to a Boussinesq-like equation, which is then subsequently approximated by two Korteweg-de Vries (KdV) equations for left- and right-traveling waves. We use the soliton solution of the KdV equation to construct approximate solutions of the CH-NLS 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 continuous-wave background. We also use numerical simulations to investigate the validity of the asymptotic solutions, study their evolution, and their head-on collisions. It is shown that small-amplitude dark and antidark solitons undergo quasi-elastic 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 finite-time 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, non-vanishing 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, so-called, "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 kink-antikink collisions, such as the critical velocity and the multi-bounce windows, and how they depend on the principal parameter of the model. We find that the critical velocity exhibits a non-monotonic dependence on the parameter controlling the curvature of the potential. For the multi-bounce 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 semi-analytical 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 condensed-matter 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 dark-bright solitons beyond the mean-field approximation is investigated. We first examine the case of a single dark-bright soliton and its oscillations within a parabolic trap. Subsequently, we move to the setting of collisions, comparing the mean-field 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 mean-field dark-bright solitons, placed off-center 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 mean-field approximation. A variety of excitations including dark-bright solitons in multiple (concurrently populated) orbitals is observed. Dark-antidark states and domain-wall-bright soliton complexes can also be observed to arise spontaneously in the beyond mean-field dynamics.
  • The existence, stability, and dynamics of bound pairs of symbiotic matter waves in the form of dark-bright soliton pairs in two-component mixtures of atomic Bose-Einstein 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 dark-bright 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 symmetry-broken stationary states that are mass imbalanced (asymmetric) with respect to their bright soliton counterpart. The long-time dynamics of both the stable and the unstable balanced and imbalanced dark-bright 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 soliton-impurity 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 "anti-particle character"), we find that dark solitons accelerate as they pass the first barrier and hence cannot be trapped by a second equal-height barrier. A pair of unequal barriers may lead to reflection from the second one, however trapping in the inter-barrier region cannot occur. We also give some examples of dynamical adjusting the barriers to trap the dark soliton in the inter-barrier 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 one-dimensional 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 steady-state solitary waveforms in a one-component, 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 high-dimensional 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 multi-dimensional 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 dark-bright solitons and their interactions, motivated by recent advances in atomic two-component repulsively interacting Bose-Einstein condensates. In particular, we study analytically via a two-soliton ansatz adopted within a variational formulation the interaction between two dark-bright solitons in a homogeneous environment beyond the integrable regime, by considering general inter/intra-atomic interaction coefficients. We retrieve the possibility of a fixed point in the case where the bright solitons are out of phase. As the inter-component interaction is increased, we also identify an exponential instability of the two-soliton 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 in-phase 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.