pith. sign in

arxiv: 2604.04460 · v3 · submitted 2026-04-06 · 🧮 math-ph · math.MP

Ground states and droplet regimes of the extended Gross-Pitaevskii equation with Lee-Huang-Yang correction

Weijie Huang , Yang Liu , Xinran Ruan This is my paper

Pith reviewed 2026-05-10 19:52 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords ground statesGross-Pitaevskii equationLee-Huang-Yang correctionquantum dropletsexistence resultsdimensional reductionnormalized gradient flow
0
0 comments X

The pith

The extended Gross-Pitaevskii equation with Lee-Huang-Yang correction admits ground states only in restricted parameter regimes that separate into no-solution, soliton-like, and droplet-like regions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives reduced one- and two-dimensional models from the three-dimensional extended Gross-Pitaevskii equation that includes the Lee-Huang-Yang correction. It then proves existence and nonexistence of ground states in free space and under confining potentials across these dimensions. Numerical solutions obtained via a normalized gradient flow method map the free-space parameter plane into distinct regimes and introduce a flat-top approximation that describes droplet profiles. These results characterize when localized states can form and how their shapes depend on model parameters.

Core claim

In free space the model parameters divide the plane into three regimes with no ground states, soliton-like states, and droplet-like states; existence theorems hold under external potentials and in lower dimensions; a flat-top approximation captures the droplet profiles; and two- and three-dimensional numerical computations illustrate the localized structures.

What carries the argument

The normalized gradient flow method with Lagrange multiplier, combined with the flat-top approximation that models the droplet regime.

If this is right

  • Ground states cease to exist outside certain parameter intervals in free space.
  • External confining potentials can create ground states where free-space solutions are absent.
  • The flat-top approximation accurately represents the shape of droplet-like ground states.
  • Ground-state profiles vary systematically with dimension and with the strength of the Lee-Huang-Yang term.
  • The normalized gradient flow algorithm reliably computes the minimizers in all regimes where they exist.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The identified regimes suggest concrete parameter windows for laboratory searches for quantum droplets in trapped gases.
  • Because the reduced models inherit the droplet behavior, one- and two-dimensional experiments can test predictions originally derived in three dimensions.
  • The same existence framework could be applied to time-dependent or multi-component versions of the equation to study droplet collisions or mixtures.
  • The flat-top approximation offers a quick way to estimate droplet size and energy without solving the full partial differential equation.

Load-bearing premise

The dimensional reduction from the three-dimensional model to one- and two-dimensional equations preserves the essential physics of the Lee-Huang-Yang correction across the studied parameter regimes.

What would settle it

A stable localized ground-state solution found numerically or observed experimentally in the free-space parameter region predicted to have no ground states would contradict the nonexistence result.

Figures

Figures reproduced from arXiv: 2604.04460 by Weijie Huang, Xinran Ruan, Yang Liu.

Figure 5
Figure 5. Figure 5: shows the [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows that, for fixed [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows the distribution of [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 5
Figure 5. Figure 5 [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows the computed ground states for these two values [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows three-dimensional visualizations of the computed g [PITH_FULL_IMAGE:figures/full_fig_p022_5.png] view at source ↗
read the original abstract

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies ground states of the extended Gross-Pitaevskii equation with Lee-Huang-Yang correction. Starting from the 3D model, it derives reduced 1D and 2D equations via nondimensionalization and dimensional reduction. It establishes existence and nonexistence theorems for ground states in free space and with confining potentials across dimensions. Numerically, a normalized gradient flow with Lagrange multiplier is used to compute profiles, revealing regimes in the free-space parameter plane (no-ground-state, soliton-like, droplet-like) and introducing a flat-top approximation for droplets, with additional 2D/3D computations shown.

Significance. If the dimensional reductions are rigorously justified and the theorems hold with the stated assumptions, the paper provides a useful bridge between variational analysis and numerics for quantum droplet models. The regime classification and flat-top approximation offer practical insights for ultracold-atom experiments, while the gradient-flow method supplies a reproducible computational tool. These elements strengthen the contribution to the study of density-dependent nonlinear Schrödinger equations.

major comments (2)
  1. [Section 2 (dimensional reduction)] The derivation of the reduced 1D/2D models (Section 2) applies transverse integration to the LHY term, yielding a density-dependent coefficient. No quantitative error bounds or scaling validation are supplied for the remainder when crossing from soliton-like to droplet-like regimes or when the external potential is removed. This assumption is load-bearing for the nonexistence theorems and free-space regime boundaries, as these may not carry over to the parent 3D functional without control on the approximation error.
  2. [Section 3 (existence/nonexistence theorems)] The existence/nonexistence results (Section 3) rely on variational methods applied to the reduced functionals. The full proofs, including handling of the density-dependent LHY coefficient and any a priori estimates needed for the free-space cases, are not fully detailed in the provided text; verification of the central claims requires these steps, especially given the reduction step above.
minor comments (2)
  1. [Section 4 (numerical method)] The description of the normalized gradient flow could include a brief statement on the discretization scheme and stopping criteria used for the numerical results.
  2. [Figures in Section 5] Figure captions for the regime diagrams would benefit from explicit parameter values or axis labels to facilitate direct comparison with the flat-top approximation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comments point by point below, indicating where revisions will be made.

read point-by-point responses

  1. Referee: [Section 2 (dimensional reduction)] The derivation of the reduced 1D/2D models (Section 2) applies transverse integration to the LHY term, yielding a density-dependent coefficient. No quantitative error bounds or scaling validation are supplied for the remainder when crossing from soliton-like to droplet-like regimes or when the external potential is removed. This assumption is load-bearing for the nonexistence theorems and free-space regime boundaries, as these may not carry over to the parent 3D functional without control on the approximation error.

    Authors: We agree that the dimensional reduction in Section 2 is performed formally via transverse integration, resulting in an effective density-dependent coefficient for the LHY term in the reduced models. This is a standard approach in the literature for deriving lower-dimensional effective equations for quantum droplets, but we acknowledge that no quantitative error bounds or scaling validation are provided for the approximation error, particularly when crossing regimes or removing the external potential. Consequently, the nonexistence theorems and free-space regime boundaries established for the reduced functionals may not directly apply to the parent 3D model without further control on the remainder. In the revised manuscript, we will add a clarifying remark in Section 2 discussing the formal nature of the reduction, its limitations, and the implications for the applicability of our results to the full 3D functional. Deriving rigorous quantitative bounds, however, constitutes a substantial separate analysis that lies outside the scope of the present work, which centers on the reduced models and their properties. revision: partial

  2. Referee: [Section 3 (existence/nonexistence theorems)] The existence/nonexistence results (Section 3) rely on variational methods applied to the reduced functionals. The full proofs, including handling of the density-dependent LHY coefficient and any a priori estimates needed for the free-space cases, are not fully detailed in the provided text; verification of the central claims requires these steps, especially given the reduction step above.

    Authors: The existence and nonexistence theorems in Section 3 are established using variational methods on the reduced energy functionals in one and two dimensions. The arguments rely on the concentration-compactness principle together with direct estimates that account for the density-dependent LHY term. While the manuscript outlines the principal steps, we accept that the handling of the variable coefficient and the a priori estimates for the free-space cases could be presented with greater detail to facilitate verification. In the revised version we will expand these proofs by inserting the missing intermediate estimates and explicitly addressing the density-dependent coefficient, thereby rendering the arguments self-contained. revision: yes

Circularity Check

0 steps flagged

No circularity: standard derivation from known 3D model via nondimensionalization and variational methods

full rationale

The paper begins from the established 3D extended Gross-Pitaevskii equation with Lee-Huang-Yang correction, performs nondimensionalization and transverse integration to obtain reduced 1D/2D models, then applies standard variational analysis for existence/nonexistence results and a normalized gradient flow scheme for numerics. These steps rely on classical PDE techniques without self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations that reduce the central claims to unverified inputs. The dimensional reduction is presented as an approximation whose validity is assumed under tight confinement, but no evidence exists that any theorem or regime boundary is forced by construction from the paper's own fitted quantities or prior self-referential results. This matches the default expectation for a self-contained mathematical physics derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard variational existence theory for nonlinear Schrödinger-type equations and convergence properties of normalized gradient flows; the dimensional reduction step invokes domain-specific assumptions about confinement that are not detailed in the abstract.

axioms (1)
  • domain assumption Standard assumptions for dimensional reduction in trapped quantum gases
    Invoked to obtain the 1D and 2D models from the 3D equation via nondimensionalization.

pith-pipeline@v0.9.0 · 5446 in / 1388 out tokens · 61016 ms · 2026-05-10T19:52:11.928965+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Antoine, A

    X. Antoine, A. Levitt, and Q. Tang, Efficient spectral computation of the stationary states of rotating Bose–Einstein condensates by preconditioned non linear conjugate gradient methods, J. Comput. Phys. 343 (2017), 92–109

    work page 2017

  2. [2]

    Bao and Y

    W. Bao and Y. Cai, Mathematical theory and numerical methods for Bose–Einste in conden- sation, Kinet. Relat. Models 6 (2013), 1–135

    work page 2013

  3. [3]

    Bao and Q

    W. Bao and Q. Du, Computing the ground state solution of Bose–Einstein conde nsates by a normalized gradient flow , SIAM J. Sci. Comput. 25 (2004), 1674–1697

    work page 2004

  4. [4]

    Bao and F

    W. Bao and F. Y. Lim, Computing ground states of spin-1 Bose–Einstein condensat es by the normalized gradient flow , SIAM J. Sci. Comput. 30 (2008), 1925–1948

    work page 2008

  5. [5]

    W. Bao, Q. Tang, and Z. Xu, Numerical methods and comparison for computing dark and bright solitons in the nonlinear Schr¨ odinger equation, J. Comput. Phys. 235 (2013), 423–445

    work page 2013

  6. [6]

    Ben Abdallah, F

    N. Ben Abdallah, F. M´ ehats, C. Schmeiser, and R. M. Weish¨ aupl,The nonlinear Schr¨ odinger equation with a strongly anisotropic harmonic potential , SIAM J. Math. Anal. 37 (2005), 189–199

    work page 2005

  7. [7]

    Cabrera, L

    C. Cabrera, L. Tanzi, J. Sanz, B. Naylor, P. Thomas, P. Cheiney , and L. Tarruell, Quantum liquid droplets in a mixture of Bose–Einstein condensates , Science 359 (2018), 301–304

    work page 2018

  8. [8]

    Chang, W.-W

    S.-M. Chang, W.-W. Lin, and S.-F. Shieh, Gauss–Seidel-type methods for energy states of a multi-component Bose–Einstein condensate , J. Comput. Phys. 202 (2005), 367–390

    work page 2005

  9. [9]

    Danaila and B

    I. Danaila and B. Protas, Computation of ground states of the Gross–Pitaevskii funct ional via Riemannian optimization , SIAM J. Sci. Comput. 39 (2017), B1102–B1129

    work page 2017

  10. [10]

    Danaila and P

    I. Danaila and P. Kazemi, A new Sobolev gradient method for direct minimization of the Gross–Pitaevskii energy with rotation , SIAM J. Sci. Comput. 32 (2010), 2447–2467. 24

    work page 2010

  11. [11]

    B. Feng, L. Cao, and J. Liu, Existence of stable standing waves for the Lee–Huang–Yang corrected dipolar Gross–Pitaevskii equation , Appl. Math. Lett. 115 (2021), 106952

    work page 2021

  12. [12]

    Ferrier-Barbut, H

    I. Ferrier-Barbut, H. Kadau, M. Schmitt, M. Wenzel, and T. Pf au, Observation of quantum droplets in a strongly dipolar Bose gas , Phys. Rev. Lett. 116 (2016), 215301

    work page 2016

  13. [13]

    E. P. Gross, Structure of a quantized vortex in boson systems , Il Nuovo Cim. 20 (1961), 454–477

    work page 1961

  14. [14]

    Hecht, O

    F. Hecht, O. Pironneau, A. Le Hyaric, and K. Ohtsuka, FreeFEM++ manual, Laboratoire Jacques Louis Lions (2005), 1–188

    work page 2005

  15. [15]

    T. D. Lee, K. Huang, and C. N. Yang, Eigenvalues and eigenfunctions of a Bose system of hard spheres and its low-temperature properties , Phys. Rev. 106 (1957), 1135–1145

    work page 1957

  16. [16]

    Liu and Y

    W. Liu and Y. Cai, Normalized gradient flow with Lagrange multiplier for compu ting ground states of Bose–Einstein condensates , SIAM J. Sci. Comput. 43 (2021), B219–B242

    work page 2021

  17. [17]

    Luo and A

    Y. Luo and A. Stylianou, Ground states for a nonlocal mixed order cubic-quartic Gros s– Pitaevskii equation , J. Math. Anal. Appl. 496 (2021), 124802

    work page 2021

  18. [18]

    Luo and A

    Y. Luo and A. Stylianou, On 3d dipolar Bose–Einstein condensates involving quantum fluc- tuations and three-body interactions , Discret. Contin. Dyn. Syst. Ser. B 26 (2021), 3455– 3477

    work page 2021

  19. [19]

    D. S. Petrov, Quantum mechanical stabilization of a collapsing Bose–Bos e mixture , Phys. Rev. Lett. 115 (2015), 155302

    work page 2015

  20. [20]

    L. P. Pitaevskii, Vortex lines in an imperfect Bose gas , J. Exp. Theor. Phys. 13 (1961), 451–454

    work page 1961

  21. [21]

    Schmitt, M

    M. Schmitt, M. Wenzel, F. B¨ ottcher, I. Ferrier-Barbut, and T. Pfau, Self-bound droplets of a dilute magnetic quantum liquid , Nature 539 (2016), 259–262

    work page 2016

  22. [22]

    Semeghini, G

    G. Semeghini, G. Ferioli, L. Masi, C. Mazzinghi, L. Wolswijk, F. Minar di, M. Modugno, G. Modugno, M. Inguscio, and M. Fattori, Self-bound quantum droplets of atomic mixtures in free space , Phys. Rev. Lett. 120 (2018), 235301

    work page 2018

  23. [23]

    X. Wu, Z. Wen, and W. Bao, A regularized Newton method for computing ground states of Bose–Einstein condensates , J. Sci. Comput. 73 (2017), 303–329

    work page 2017

  24. [24]

    Zhuang and J

    Q. Zhuang and J. Shen, Efficient SAV approach for imaginary time gradient flows with applications to one- and multi-component Bose-Einstein co ndensates, J. Comput. Phys. 396 (2019), 72–88. 25

    work page 2019