We explore the dynamics and interactions of multiple bright droplets and bubbles, as well as the interactions of kinks with droplets and with antikinks, in the extended one-dimensional Gross–Pitaevskii model including the Lee–Huang–Yang correction. Existence regions are identified for the one-dimensional droplets and bubbles in terms of their chemical potential, verifying the stability of the droplets and exposing the instability of the bubbles. The limiting case of the droplet family is a stable kink. The interactions between droplets demonstrate in-phase (out-of-phase) attraction (repulsion), with the so-called Manton’s method explicating the observed dynamical response, and mixed behavior for intermediate values of the phase shift. Droplets bearing different chemical potentials experience mass-exchange phenomena. Individual bubbles exhibit core expansion and mutual attraction prior to their destabilization. Droplets interacting with kinks are absorbed by them, a process accompanied by the emission of dispersive shock waves and gray solitons. Kink–antikink interactions are repulsive, generating counter-propagating shock waves. Our findings reveal dynamical features of droplets and kinks that can be detected in current experiments.
We explore the dynamics and interactions of multiple bright droplets and bubbles, as well as the interactions of kinks with droplets and with antikinks, in the extended one-dimensional Gross–Pitaevskii model including the Lee–Huang–Yang correction. Existence regions are identified for the one-dimensional droplets and bubbles in terms of their chemical potential, verifying the stability of the droplets and exposing the instability of the bubbles. The limiting case of the droplet family is a stable kink. The interactions between droplets demonstrate in-phase (out-of-phase) attraction (repulsion), with the so-called Manton’s method explicating the observed dynamical response, and mixed behavior for intermediate values of the phase shift. Droplets bearing different chemical potentials experience mass-exchange phenomena. Individual bubbles exhibit core expansion and mutual attraction prior to their destabilization. Droplets interacting with kinks are absorbed by them, a process accompanied by the emission of dispersive shock waves and gray solitons. Kink–antikink interactions are repulsive, generating counter-propagating shock waves. Our findings reveal dynamical features of droplets and kinks that can be detected in current experiments.
We unravel stationary states in the form of dark soliton stripes, bubbles, and kinks embedded in a two-dimensional droplet-bearing setting emulated by an extended Gross-Pitaevskii approach. The existence of these configurations is corroborated through an effectively reduced potential picture demonstrating their concrete parametric regions of existence. The excitation spectra of such configurations are analyzed within the Bogoliubov--de Gennes framework exposing the destabilization of dark soliton stripes and bubbles, while confirming the stability of droplets, and importantly unveiling spectral stability of the kink against transverse excitations. Additionally, a variational approach is constructed providing access to the transverse stability analysis of the dark soliton stripe for arbitrary chemical potentials and widths of the structure. This is subsequently compared with the stability analysis outcome demonstrating very good agreement at small wave numbers. Dynamical destabilization of dark soliton stripes via the snake instability is showcased, while bubbles are found to feature both a splitting into a gray soliton pair and a transverse instability thereof. These results shed light on unexplored stability and instability properties of nonlinear excitations in environments featuring a competition of mean-field repulsion and beyond-mean-field attraction that can be probed by state-of-the-art experiments.
We study nonlinear holes, bubbles, kink fronts, and quantum-droplet-like states in an elongated Bose mixture governed by an extended Gross-Pitaevskii equation with attractive cubic mean-field interactions and repulsive quartic Lee-Huang-Yang contribution. We classify the constant-amplitude backgrounds: the Bogoliubov-de Gennes spectrum shows that the upper branch is modulationally stable, whereas the lower branch is unstable to long-wavelength perturbations. Grand-potential and pressure determine the vacuum-finite-density coexistence point and the unique background supporting stationary kink and antikink fronts. Symmetric holes on the stable background have a density deficit and integral width that diverge near coexistence as the hole separates into two nonlinear fronts. A multiplicative kink-antikink construction connects bell-shaped states to broad flat-top droplets and approaches the stationary quantum-droplet family as the front separation increases. An antikink-kink sum with equal signs generates bubble-like density depletions without a phase jump, whereas their difference produces dark-hole-like profiles with a π phase difference between the asymptotic backgrounds. Bubble depth increases with front separation, distinguishing gray-like and dark-like profiles. Real-time simulations show that deep bubbles persist under weak random perturbations, demonstrating finite-time robustness without establishing spectral stability. In the asymmetric two-component model, localized relative-density modes can induce spinor-driven breakup on a modulationally stable background, while robust holes persist near the kink limit. Collisions of kink-antikink droplet-like states display phase-dependent merging, asymmetric particle transfer, and effective reflection. These results connect nonlinear fronts with bubbles, phase-jump holes, and self-bound droplets in beyond-mean-field Bose mixtures.
Abstract We address the existence and stability of one-dimensional (1D) holes and kinks and two-dimensional (2D) vortex-holes nested in extended binary Bose mixtures, which emerge in the presence of Lee–Huang–Yang (LHY) quantum corrections to the mean-field energy, along with self-bound quantum droplets. We consider both the symmetric system with equal intra-species scattering lengths and atomic masses, modeled by a single (scalar) LHY-corrected Gross–Pitaevskii equation (GPE), and the general asymmetric case with different intra-species scattering lengths, described by two coupled (spinor) GPEs. We found that in the symmetric setting, 1D and 2D holes can exist in a stable form within a range of chemical potentials that overlaps with that of self-bound quantum droplets, but that extends far beyond it. In this case, holes are found to be always stable in 1D and they transform into pairs of stable out-of-phase kinks at the critical chemical potential at which localized droplets turn into flat-top states, thereby revealing the connection between localized and extended nonlinear states. In contrast, we found that the spinor nature of the asymmetric systems may lead to instability of 1D holes, which tend to break into two gray states moving in the opposite directions. Remarkably, such instability arises due to spinor nature of the system and it affects only holes nested in extended modulationally-stable backgrounds, while localized quantum droplet families remain completely stable, even in the asymmetric case, while 1D holes remain stable only close to the point where they transform into pairs of kinks. We also found that symmetric systems allow fully stable 2D vortex-carrying single-charge states at moderate amplitudes, while unconventional instabilities appear also at high amplitudes. Symmetry also strongly inhibits instabilities for double-charge vortex-holes, which thus exhibit unexpectedly robust evolutions at low amplitudes.
We derive a channel-resolved Lee-Huang-Yang correction and construct an extended GrossPitaevskii model for one-dimensional polar spin-1 quantum droplets. The fluctuation contribution separates into density and spin channels, which supports self-bound droplets even when the spinindependent mean-field interaction is repulsive. Stationary solutions exhibit a continuous crossover from soliton-like to flat-top droplets, accompanied by saturation of the chemical potential and peak density as the particle number increases. Within the parameter range examined here, linear Bogoliubov analysis together with weak-perturbation dynamics supports the stability of both droplet types. A quadratic-Zeeman quench reveals a finite-size crossover in breathing dynamics and distinct nonequilibrium roles of the density and spin fluctuation channels. Representative head-on collisions further show that the finite-size crossover modulates phase-sensitive nonlinear scattering, with in-phase impact producing coalescence-like central retention and out-of-phase impact favoring quasi-elastic separation. The analysis clarifies how density and spin fluctuations shape equilibrium structure and nonequilibrium response in low-dimensional polar spinor droplets.
We unravel the existence and stability properties of dark soliton solutions as they extend from the regime of trapped quantum droplets towards the Thomas-Fermi limit in homonuclear symmetric Bose mixtures. Leveraging a phase-plane analysis, we identify the regimes of existence of different types of quantum droplets and subsequently examine the possibility of black and gray solitons and kink-type structures in this system. Moreover, we employ the Landau dynamics approach to extract an analytical estimate of the oscillation frequency of a single dark soliton in the relevant extended Gross--Pitaevskii model. Within this framework, we also find that the single soliton immersed in a droplet is stable, while multisoliton configurations exhibit parametric windows of oscillatory instabilities. Our results pave the way for studying dynamical features of nonlinear multisoliton excitations in a droplet environment in contemporary experimental settings.
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 intraspecies 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\"odinger 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 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.