This paper presents a numerical method for computing the ground states of the three‑dimensional Gross-Pitaevskii equation （GPE） that includes the Lee-Huang-Yang （LHY） correction.The LHY term in the Gross-Pitaevskii energy functional is of the form |ψ|5 <?fx-imagestate width="4.40266657" height="3.89466691"?> <?fx-imagestate width="4.40266657" height="3.89466691"?> and accounts for quantum fluctuations in Bose-Einstein condensates （BECs）.First， the continuous projected gradient flow is discretized with a Lagrange multiplier by using a Fourier pseudo-spectral method in space and a backward Euler scheme in time， obtaining a linearly implicit backward Euler scheme （PGF-BESP）.The key idea is to implicitly determine the Lagrange multiplier through an orthogonal condition， thereby guaranteeing mass conservation （the normalization constraint） and providing an explicit expression for the Lagrange multiplier directly from the orthogonal condition.This preserves the two fundamental properties of the continuous normalized gradient flow with a Lagrange multiplier： Mass conservation and energy dissipation.Next， we carry out a rigorous theoretical analysis of the scheme and prove that the discrete energy is monotonically non-increasing under a suitable condition on the time step size.We then perform numerical experiments to verify the theoretical results， and the results show that mass is indeed conserved and the energy decays monotonically.Finally， we use the proposed method to compute ground state configurations for different potential terms and different interaction parameters.We find that as the interaction parameter increases， the overall size of the ground state density distribution also increases， and under the optical lattice potential， multi-layer structures are observed.In summary， both theoretical and numerical results show that the proposed PGF-BESP method provides a stable， efficient and easy to implement numerical tool for ground state computations of GPE with LHY correction.
In this paper, we present a systematic study on the ground state computation of quantum droplets in homonuclear Bose-Bose mixtures, governed by the extended Gross-Pitaevskii equations (eGPEs) with Lee-Huang-Yang (LHY) corrections. This model captures the formation of self-bound droplets stabilized by the delicate balance between the attractive mean-field interaction and the repulsive quantum fluctuations. We formulate dimensionless energy functionals for both the general two-component system and the reduced single-component density-locked model. To compute the ground states efficiently, we adapt and benchmark various gradient flow discretization schemes, identifying a backward-forward sine-pseudospectral scheme based on the gradient flow with Lagrange multiplier method (GFLM-BFSP) as the robust solver for our simulations. Utilizing this method, we report three main numerical observations: (i) the density-locked model is quantitatively validated as a reliable approximation for ground state properties; (ii) the dimension-dependent convergence rates of the Thomas-Fermi approximation are established in the strong-coupling regime; and (iii) the critical particle number for self-binding in free space is numerically determined, providing a precise correction to the analytical prediction by Petrov [Phys. Rev. Lett. 115, 155302 (2015)].
We investigate the existence and dynamics of two-dimensional solitary waves in a quantum droplet environment described by the extended Gross-Pitaevskii equation featuring logarithmic mean-field and Lee-Huang-Yang interactions. In the modulationally stable regime of the background, we employ suitable multiscale asymptotic methods to derive effective nonlinear integrable models corresponding to the Kadomtsev-Petviashvili and Davey-Stewartson equations. Based on these reduced models, we construct approximate analytical solutions describing line solitons, algebraically localized lump solitons, ring solitons, and exponentially localized dromions embedded on the droplet background. The dynamical robustness of these solutions is monitored through numerical simulations. Line, lump and ring solitons stay closest to the theoretical predictions, although progressively deviate due to the emergence of small-amplitude radiation, while dromions depart from their analytical waveform the most, although they roughly maintain their shape. Our results unveil unprecedented multidimensional soliton solutions in models featuring the competition of mean-field and quantum fluctuations and as such are amenable to current ultracold atom experiments.
Abstract The static and dynamic properties of self-bound quantum droplets in a one-dimensional Bose-Bose mixture are discussed in the spirit of the Hartree–Fock-Bogoliubov theory. This latter enables us to provide beyond the Lee-Huang-Yang (LHY) quantum corrections to the equation of state at both zero and finite temperatures. In the uniform case our results for the ground-state energy and the critical temperature are confirmed through comparison with Quantum Monte-Carlo simulation and with available theoretical results. The density profiles are supported by numerical simulations of the generalized Gross-Pitaevskii equation which selfconsistently includes higher-order terms originating from the normal and anomalous fluctuations under the local density approximation. We show that the density exhibits a dip near its center in the flat-top plateau region for large interspecies interactions. We exemplify the impact of the beyond LHY corrections on the spatiotemporal evolution of the self-bound droplet in the presence of excitation induced by periodic density modulation. It is found that higher-order corrections may lead to the formation of a train of small droplets. We then extend our study for the case of inhomogeneous droplets in quasi one-dimensional Bose mixtures.
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.