We study the ground states of the extended Gross--Pitaevskii equation with the Lee--Huang--Yang correction from both theoretical and numerical perspectives. Starting from the three-dimensional model, we derive reduced one- and two-dimensional equations through nondimensionalization and dimensional reduction. We establish existence and nonexistence results for ground states in different spatial dimensions, both in free space and under confining external potentials. For the numerical computation of ground states, we propose a normalized gradient flow method with a Lagrange multiplier. The numerical results show how the model parameters affect the ground-state profiles, and reveal different regimes in the free-space parameter plane, including no-ground-state, soliton-like, and droplet-like regions. We also introduce a simple flat-top approximation for the droplet regime and present two- and three-dimensional computations to illustrate more general localized structures.
It was recently found that the Lee-Huang-Yang (LHY) correction to the mean-field Hamiltonian of binary atomic boson condensates suppresses the collapse and creates stable localized modes (two-component quantum droplets, QDs) in two and three dimensions (2D and 3D). In particular, the LHY effect modifies the effective Gross-Pitaevskii equation (GPE) in 2D by adding a logarithmic factor to the usual cubic term. In the framework of the accordingly modified two-component GPE system, we construct 2D self-trapped modes in the form of QDs with vorticity
S embedded into each component. Because of the effect of the logarithmic factor, the QDs feature a flat-top shape, which expands with the increase of
S and norm
N. An essential finding, produced by a systematic numerical investigation and analytical estimates, is that the vortical QDs are stable (which is a critical issue for vortex solitons in nonlinear models) up to
S=5, for
N exceeding a certain threshold value, which is predicted to scale as ${N}_{\mathrm{th}}\ensuremath{\sim}{S}^{4}$ for large
S (for three-dimensional QDs, the scaling is ${N}_{\mathrm{th}}\ensuremath{\sim}{S}^{6}$). The prediction is corroborated by numerical findings. Pivots of QDs with $S\ensuremath{\ge}2$ are subject to structural instability, as specially selected perturbations can split the single pivot in a set of
S or
S+2 pivots corresponding to unitary vortices; however, the structural instability remains virtually invisible, as it occurs in a broad central ``hole'' of the vortex soliton, where values of fields are very small, and it does not cause any dynamical instability. In the condensate of
39K atoms, in which QDs with
S=0 and a quasi-2D shape were created recently, the vortical droplets may have radial size $\ensuremath{\lesssim}30\phantom{\rule{0.16em}{0ex}}\ensuremath{\mu}\mathrm{m}$, with the number of atoms in the range of
104−−105. The role of three-body losses is considered too, demonstrating that they do not prevent the creation of the vortex droplets but may produce a noteworthy effect, leading to sudden splitting of ``light'' droplets. In addition, hidden-vorticity (HV) states in QDs, with topological charges ${S}_{+}=\ensuremath{-}{S}_{\ensuremath{-}}=1$ in their components, which are prone to strong instability in other settings, have their stability region too. Unstable HV states tend to spontaneously merge into zero-vorticity solitons. Collisions of QDs, which may lead to their merger, and dynamics of elliptically deformed QDs (which form rotating elongated patterns or ones with oscillations of the eccentricity) are briefly considered too.
We study the properties of modulational instability and discrete breathers arising in a quasi-one-dimensional discrete Gross-Pitaevskii equation with Lee-Huang-Yang corrections. Conditions for modulation instability and instability regions of nonlinear plane waves are determined in parameter space. We analytically investigate the existence of different quantum droplet solutions, including intersite, onsite, front-like, flat-top, and dark localized modes, using the Page method and variational approach. Their stability is checked using linear stability analyses and numerical simulations. The analytical predictions corroborated with the numerical simulations.
This mini-review collects results predicting the creation of matter-wave solitons by the spinor system of Gross-Pitaevskii equations (GPEs) with the self-attractive cubic nonlinearity and linear first-order-derivative terms accounting for the spin-orbit coupling (SOC). In 1D, the so-predicted bright solitons are similar to usual ones, supported by the GPE in the absence of SOC. Essentially new results were recently obtained for 2D and 3D systems: SOC suppresses the collapse instability in the multidimensional GPE, creating 2D ground-state solitons and metastable 3D ones of two types: semi-vortices (SVs), with vorticities m = 1 in one component and m = 0 in the other, and mixed modes (MMs), with m = 0 and present in both components. With the Galilean invariance broken by SOC, moving solitons exist up to a certain critical velocity. The latest result predicts stable 2D “quantum droplets” of the MM type in the presence of the Lee-Huang-Yang corrections to the GPE system, induced by quantum fluctuations, in the case when the inter-component attraction dominates over the self-repulsion in each component.
We investigate, both analytically and numerically, the nonlinear dynamics of (2+1)-dimensional [(2+1)D] matter waves excited in a disk-shaped dipolar Bose-Einstein condensate (BEC) when quantum fluctuations described by the Lee-Huang-Yang (LHY) correction are taken into consideration. By using a method of multiple scales, we derive Davey-Stewartson I equations that govern the nonlinear evolution of matter-wave envelopes. We demonstrate that the system supports (2+1)D matter-wave dromions, which are superpositions of a short-wavelength excitation and a long-wavelength mean flow. We found that the stability of the matter-wave dromions can be enhanced by the LHY correction. We also found that such dromions display interesting behaviors of collision, reflection, and transmission when they interact with each other and are scattered by obstacles. The results reported here are useful not only for improving the understanding on the physical property of the quantum fluctuations in BECs, but also for possible experimental findings of new nonlinear localized excitations in systems with long-ranged interactions.
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.