Mathematical and numerical studies on ground states of the extended Gross-Pitaevskii equation with the Lee-Huang-Yang correction

This paper investigates the ground states of the extended Gross-Pitaevskii equation with Lee-Huang-Yang correction by deriving reduced dimensional models, establishing theoretical existence results, and employing a normalized gradient flow method to numerically characterize distinct regimes such as soliton-like and droplet-like structures.

The Gist

Imagine a giant, invisible cloud of atoms so cold that they stop acting like individual particles and start behaving like a single, giant wave. This is a Bose-Einstein Condensate (BEC), often called a "super-atom."

For decades, physicists used a standard rulebook (the Gross-Pitaevskii equation) to predict how these clouds behave. But recently, scientists discovered a new phenomenon: Quantum Droplets. These are tiny, self-contained blobs of atoms that hold themselves together without any external container, much like a drop of water on a leaf, but made of pure quantum mechanics.

This paper is a guidebook for understanding these mysterious droplets. The authors, Huang, Liu, and Ruan, did two main things: they built a mathematical theory to prove these droplets can exist, and they built a computer program to simulate what they look like.

Here is the breakdown of their work using simple analogies:

1. The New Rulebook: The "Lee-Huang-Yang" Correction

Think of the atoms in the cloud as people at a party.

• The Old Rule: The standard model says the people either like each other (attraction) or dislike each other (repulsion). If they like each other too much, the whole crowd collapses into a tiny ball. If they dislike each other, they fly apart.
• The New Twist (LHY Correction): The authors added a new rule based on the "Lee-Huang-Yang" correction. Imagine that when the crowd gets very crowded, a new force kicks in—a "quantum jitter." This jitter acts like a safety net. It pushes back just enough to stop the crowd from collapsing completely, but not so much that they fly apart.
• The Result: This balance creates a stable, self-contained "droplet" that floats in empty space.

2. The Math: Proving the Droplet Exists

Before you can build a house, you need to know the ground is solid. The authors spent a lot of time proving that these droplets are mathematically possible.

• The "Goldilocks" Zone: They found that these droplets only exist under specific conditions.
  • If the attraction is too weak, the cloud flies apart.
  • If the attraction is too strong (and the "safety net" is too weak), the cloud collapses.
  • Just Right: There is a specific "Goldilocks" zone where the attraction and the quantum jitter balance perfectly.
• The Dimensions: They checked if this works in 1D (a line), 2D (a flat sheet), and 3D (a sphere). They found that in 3D (our real world), the droplets only form if you have enough "mass" (enough atoms). If you have too few atoms, they just can't hold themselves together.

3. The Computer Simulation: The "Gradient Flow"

How do you find the shape of these droplets on a computer? You can't just guess. The authors used a method called Normalized Gradient Flow.

• The Analogy: Imagine you are blindfolded on a hilly landscape, and your goal is to find the lowest valley (the ground state). You take a step downhill, check your balance, and take another step.
• The Twist: In this quantum world, you are also holding a bucket of water (the atoms). You are not allowed to spill a drop or add any water; the total amount must stay exactly the same.
• The Method: The authors created a computer algorithm that acts like a hiker who constantly checks their water level. If they lose a drop, they instantly add one back. They keep walking downhill until they find the deepest, most stable valley. This valley represents the shape of the quantum droplet.

4. What the Droplets Look Like

The simulations revealed two distinct "personalities" for these clouds, depending on the settings:

• The Soliton (The Smooth Hill): In some conditions, the cloud looks like a smooth, rounded hill. It's dense in the middle and fades out gently at the edges.
• The Droplet (The Flat-Topped Plateau): In other conditions (stronger attraction), the cloud changes shape. It becomes a flat-topped pancake. The density is high and uniform in the center, then drops off sharply at the edges.
  • The "Flat-Top" Approximation: The authors realized that for these flat-topped droplets, you can approximate them as a simple cylinder or sphere with a flat top. This makes calculating their energy much easier, like approximating a complex mountain with a simple box.

5. Why This Matters

• Understanding Nature: This helps us understand how matter behaves at the coldest temperatures in the universe.
• New Materials: These droplets are stable without being trapped in a magnetic cage. This could lead to new types of sensors or even new ways to store information.
• Better Tools: The computer code the authors wrote is a powerful tool. It can handle complex shapes (like droplets in 3D space) and adapt its grid to zoom in on the sharp edges where the action happens, ensuring the results are accurate.

Summary

In short, this paper is the architectural blueprint and the construction manual for quantum droplets.

1. Theory: They proved the blueprints are valid and defined the exact conditions needed to build them.
2. Simulation: They built a digital crane (the algorithm) to construct these droplets on a computer.
3. Discovery: They found that these droplets can look like smooth hills or flat-topped pancakes, depending on how you tune the "knobs" of attraction and repulsion.

It's a beautiful mix of pure math and practical computing that helps us visualize a strange, new state of matter.

Technical Summary

1. Problem Statement

The paper investigates the ground states of the extended Gross-Pitaevskii (eGP) equation incorporating the Lee-Huang-Yang (LHY) correction. This model describes ultracold Bose gases, specifically focusing on the formation of quantum droplets, which arise from a balance between attractive mean-field interactions and repulsive quantum fluctuations.

The governing equation (in physical units) is:
$$i\hbar\partial_t\psi = \left[ -\frac{\hbar^2}{2m}\nabla^2 + V(x) + g|\psi|^2 + g_{LHY}|\psi|^3 \right] \psi$$
where:

• $g|\psi|^2$ represents the standard cubic mean-field interaction (attractive if $g<0$).
• $g_{LHY}|\psi|^3$ is the LHY correction term (repulsive, stabilizing the droplet).
• The system is subject to a mass constraint $\|\psi\|_2^2 = N$ (conservation of particle number).

The study addresses two main challenges:

1. Theoretical: Establishing rigorous existence and non-existence conditions for ground state minimizers of the associated energy functional across different spatial dimensions ($d=1, 2, 3$) and parameter regimes (free space vs. confining potentials).
2. Numerical: Developing stable and accurate algorithms to compute these ground states, particularly given the high-order nonlinearity ($|\psi|^5$ in energy) and the potential for strong localization (droplets) which challenges standard numerical resolution.

2. Methodology

A. Dimensional Reduction and Nondimensionalization

The authors start from the 3D model and derive reduced 1D and 2D equations under strong anisotropic confinement (harmonic traps).

• Scaling: The equations are nondimensionalized to unify the form across dimensions.
• Unified Form: The resulting dimensionless equation is:
  $$i\partial_t\psi = -\frac{1}{2}\nabla^2\psi + V(x)\psi + \beta|\psi|^2\psi + \lambda|\psi|^3\psi$$
  where $\beta$ and $\lambda$ are dimensionless interaction parameters.
• Reduction: By assuming Gaussian profiles in strongly confined directions, the 3D model is reduced to effective 2D (disk-shaped) and 1D (cigar-shaped) equations with renormalized coefficients.

B. Theoretical Analysis (Existence/Non-existence)

The authors analyze the minimization of the energy functional $E(\psi)$ under the mass constraint $\|\psi\|_2 = c$.

• Free Space ($V \equiv 0$):
  • If $\beta \ge 0$, no ground state exists (energy is unbounded below or infimum is 0 but not attained).
  • If $\beta < 0$:
    • 1D: Ground states exist for all masses $c > 0$.
    • 2D/3D: A critical mass $c^*$ exists. No ground state exists for $c < c^*$ (energy infimum is 0 but not attained). Ground states exist for $c > c^*$.
• Confining Potential ($V(x) \to \infty$): Ground states exist for all $c > 0$ and $\beta \in \mathbb{R}$ due to the compactness provided by the potential.
• Proof Techniques: The proofs utilize scaling arguments, interpolation inequalities (Gagliardo-Nirenberg, Sobolev), and the concentration-compactness principle.

C. Numerical Method

To compute ground states, the authors propose a Normalized Gradient Flow with Discrete Normalization (GFDN) combined with a Finite Element Method (FEM).

• Gradient Flow: An imaginary-time evolution equation is used:
  $$\partial_t \phi = \frac{1}{2}\nabla^2\phi - V\phi - \beta|\phi|^2\phi - \lambda|\phi|^3\phi + \mu(t)\phi$$
  where $\mu(t)$ is a Lagrange multiplier enforcing the mass constraint.
• Discretization Strategy:
  • Time: Semi-discrete scheme using an implicit-explicit (IMEX) approach. The linear part is treated implicitly, while nonlinear terms are handled semi-implicitly to ensure stability and preserve non-negativity.
  • Space: Finite Element Method (FEM) with linear elements. This is chosen over finite differences for its ability to handle adaptive mesh refinement, which is crucial for resolving the sharp transition layers of quantum droplets.
  • Symmetry: For symmetric cases (radial/spherical), the problem is reduced to 1D using finite differences to save computational cost.

3. Key Results

A. Theoretical Findings

• Existence Thresholds: The paper rigorously proves that in free space, ground states only exist if the attractive interaction ($\beta < 0$) is strong enough relative to the mass (in 2D/3D) or for any mass (in 1D).
• Non-existence: For small masses in 2D/3D free space, the system cannot form a stable self-bound droplet; the energy infimum is zero but is not attained by any function in the Sobolev space.

B. Numerical Findings

• Parameter Dependence:
  • Increasing the attractive strength ($|\beta|$) or decreasing the LHY repulsion ($\lambda$) leads to more localized, higher-density profiles.
  • Increasing the mass ($c$) broadens the profile but maintains a "flat-top" structure characteristic of droplets.
• Phase Diagram: The authors map the $(\beta, \lambda)$ parameter plane into three distinct regimes:
  1. No Ground State: Parameters where the system is unstable.
  2. Soliton-like: Localized states with smooth decay (no flat top).
  3. Droplet-like: States with a high-density core and a distinct flat-top profile.
• Flat-Top Approximation: In the droplet regime (strong attraction/weak LHY), the ground state is approximated by a constant density on a support domain $D$. The authors derive an analytical estimate for the peak density and energy, showing that the relative error decreases as $|\beta| \to \infty$ or $\lambda \to 0$.
• 2D/3D Structures:
  • Free Space: Spherically symmetric droplets with sharp boundaries.
  • Optical Lattice: The method successfully captures multi-peak structures where the lattice potential creates separated localized islands.
  • Anisotropic Traps: The method resolves deformation of droplets (e.g., cigar-shaped or disk-shaped) under anisotropic harmonic confinement.

4. Significance and Contributions

1. Rigorous Mathematical Foundation: The paper provides the first comprehensive existence/non-existence theory for the eGP equation with LHY correction across all dimensions and potential types, clarifying the stability conditions for quantum droplets.
2. Robust Numerical Framework: The proposed GFDN-FEM scheme effectively handles the high-order nonlinearity and the stiffness of the problem. The use of adaptive meshing is critical for accurately capturing the "flat-top" density profiles and sharp interfaces of quantum droplets, which are often missed by uniform grid methods.
3. Physical Insight: The construction of the phase diagram and the validation of the flat-top approximation provide a clear physical understanding of the transition between soliton-like and droplet-like states.
4. Computational Efficiency: By reducing symmetric problems to 1D and utilizing adaptive FEM for general cases, the authors demonstrate a scalable approach for studying complex 3D quantum droplet configurations.

In summary, this work bridges the gap between the theoretical analysis of quantum droplet stability and the practical numerical simulation of their complex spatial structures, offering a unified framework for future studies in ultracold atomic physics.