The Lee-Huang-Yang energy for a dilute gas of hard spheres: an upper bound

This paper establishes an upper bound for the ground state energy density of a dilute Bose gas of hard spheres in the thermodynamic limit that matches the celebrated Lee-Huang-Yang formula up to lower-order terms.

The Gist

Imagine a crowded dance floor inside a giant, invisible box. The dancers are bosons—a special type of quantum particle that loves to be together, almost like a synchronized swimming team. However, these dancers are also hard spheres. This means they are like rigid billiard balls: they can't occupy the same space. If two dancers get too close (closer than a tiny radius $a$), they bounce off each other instantly. They cannot overlap.

The physicists in this paper, Giulia Basti and her team, are trying to answer a very specific question: How much energy does it take to keep this dance floor moving when the room is huge and the dancers are very far apart from each other?

In the world of quantum physics, this "energy" is the ground state energy. It's the minimum energy the system has when it's as calm as possible.

The Famous Recipe (The Lee-Huang-Yang Formula)

Back in 1957, three brilliant physicists (Lee, Huang, and Yang) came up with a recipe to predict this energy. They said the energy depends on two main things:

1. The Density: How many dancers are on the floor.
2. The "Hardness": How big the billiard balls are.

Their formula had two parts:

• The Main Dish: A big, obvious chunk of energy that everyone could agree on.
• The Secret Sauce: A tiny, subtle correction (the "Lee-Huang-Yang term") that accounts for the complex way the dancers interact when they are just barely touching.

For decades, mathematicians could prove the "Main Dish" was correct. But proving the "Secret Sauce" for hard spheres (the billiard balls) was a nightmare. Previous attempts worked for "soft" interactions (like dancers who gently push each other), but the rigid, hard-core rule of the billiard balls broke the math.

The Problem: The "Hard" Wall

The difficulty with hard spheres is that they have a strict rule: Zero probability of overlap.

• If you try to calculate the energy using standard math tricks, you hit a wall (literally). The math explodes because the "hard core" condition is so abrupt.
• Previous attempts to fix this used a "patchwork" method: they built small models of the dance floor and glued them together. But this didn't quite capture the correct "Secret Sauce" constant for the hard spheres.

The Solution: A New Dance Move

The authors of this paper (Basti, Brooks, Cenatiempo, Olgiati, and Schlein) decided to stop patching things together. Instead, they worked directly on the entire infinite dance floor (the thermodynamic limit).

Here is how they did it, using a creative analogy:

1. The Jastrow Factor (The "Personal Space" Bubble)

First, they acknowledged that every dancer needs personal space. They created a mathematical "bubble" around each dancer (called a Jastrow factor) that forces other dancers to stay away if they get too close. This handles the "hard sphere" rule perfectly.

2. The Bogoliubov Transformation (The "Group Sync")

But personal space isn't enough. When the dancers move, they don't just move individually; they move in waves. If one dancer moves, the whole crowd ripples.

• The authors added a second layer of math called a Bogoliubov transformation. Think of this as a choreographer who organizes the dancers into synchronized waves.
• This handles the long-range interactions (the "ripples" across the dance floor) that the simple "personal space" bubble missed.

3. The Grand Canonical Ensemble (The "Infinite Crowd" Trick)

Usually, calculating the energy of a fixed number of dancers is hard because you have to count every single one. The authors used a clever trick: they imagined a magic door where dancers could enter and leave the box freely, as long as the average number of dancers stayed the same.

• This is called the Grand Canonical Ensemble.
• Why is this helpful? It removes the headache of "normalization" (making sure the total probability adds up to 1). In their previous attempts, the math got messy because the "bubble" and the "waves" fought over the normalization. By letting the crowd size fluctuate slightly, the math became much cleaner, allowing them to isolate the "Secret Sauce" term perfectly.

The Result: The Formula is Proven!

By combining the "Personal Space" bubble (for the hard walls) with the "Group Sync" waves (for the long-range ripples) and using the "Magic Door" trick, they finally calculated the energy.

They found that the energy density matches the famous Lee-Huang-Yang formula exactly, right down to the "Secret Sauce" term:
$$\text{Energy} \approx \text{Main Dish} + \frac{128}{15\sqrt{\pi}} (\text{Density} \times \text{Size})^{1/2}$$

Why Does This Matter?

1. It's a Victory for Rigor: For the first time, we have a mathematically airtight proof that this famous 1957 formula works for the most extreme case: hard spheres.
2. It Connects Theory to Reality: While "hard spheres" are an idealization, they are a great model for real-world atoms that repel each other strongly. This gives us confidence that our understanding of quantum gases (like those used in super-cold experiments) is correct.
3. It Solves a 60-Year Puzzle: It fills a gap in the literature that had been open since the days of the early quantum pioneers.

In short: The authors built a perfect mathematical model of a crowded dance floor of billiard balls, figured out exactly how much energy it takes to keep them dancing, and proved that the "Secret Sauce" predicted by Lee, Huang, and Yang in 1957 was correct all along.

Technical Summary

1. Problem Statement

The paper addresses a fundamental problem in quantum many-body physics: determining the ground state energy density of a dilute Bose gas consisting of $N$ hard spheres of radius $a$ in a box $\Lambda = [0, L]^3$ with periodic boundary conditions.

The goal is to rigorously derive the Lee-Huang-Yang (LHY) formula for the ground state energy per unit volume, $e(\rho)$, in the thermodynamic limit ($N, L \to \infty$ with fixed density $\rho = N/L^3$) and the dilute limit ($\rho a^3 \ll 1$). The formula is given by:
$$e(\rho) = 4\pi a \rho^2 \left[ 1 + \frac{128}{15\sqrt{\pi}}(\rho a^3)^{1/2} + \dots \right]$$
While the leading order term ($4\pi a \rho^2$) was established decades ago (Dyson for upper bound, Lieb-Yngvason for lower bound), and the LHY correction term was recently proven as a lower bound for hard spheres (Benedikter et al., [20]), a matching upper bound with the correct constant had remained an open problem. Previous upper bounds for hard spheres either failed to capture the correct LHY constant or required integrable potentials, excluding the physically relevant hard-sphere case.

2. Methodology

The authors employ a sophisticated variational approach within the grand canonical ensemble on the bosonic Fock space $\mathcal{F}$. The core of their methodology involves constructing a highly refined trial state $\Psi$ that captures correlations at multiple length scales.

A. The Trial State Construction

The trial state is defined as:
$$\Psi = Z^{-1} J W(\rho_0) T \Omega$$
where:

• $\Omega$ is the vacuum vector.
• $W(\rho_0)$ is a Weyl operator generating a coherent state modeling a Bose-Einstein condensate (BEC) with density $\rho_0$.
• $J$ is a Jastrow factor operator implemented via multiplication by a function $f_\ell(x)$, which enforces the hard-core condition ($\psi=0$ if $|x_i - x_j| \le a$) and describes short-range correlations up to a scale $\ell \gg a$.
• $T$ is a Bogoliubov transformation (unitary operator) designed to capture long-range correlations up to and slightly beyond the healing length $\xi \sim (\rho a)^{-1/2}$.
• $Z$ is a normalization constant.

B. Handling Hard Spheres

The primary challenge with hard spheres is that standard Bogoliubov transformations (which are quasi-free states) do not naturally satisfy the hard-core condition. The authors solve this by:

1. Using the Jastrow factor $J$ to strictly enforce the hard-core condition on short scales ($a \ll |x| \lesssim \ell$).
2. Using the Bogoliubov transformation $T$ to model the "soft" correlations on larger scales ($\ell \lesssim |x| \lesssim \ell_0$), where $\ell_0$ is slightly larger than the healing length.
3. Carefully choosing the kernel of the Bogoliubov transformation to diagonalize the effective quadratic Hamiltonian emerging from the interaction, ensuring the correct LHY term is captured.

C. Key Technical Innovations

• Direct Thermodynamic Limit: Unlike previous works that constructed trial states by patching together small boxes, this work works directly in the thermodynamic box, avoiding issues with boundary conditions and finite-size effects.
• Operator Identities: The authors derive a crucial identity for the action of the annihilation operator $a_x$ on the trial state $\Psi$ (Lemma 2.3):
  $$a_x \Psi = \sqrt{\rho_0} J(x) \Psi + J(x) \int dy \, (\gamma^{-1} * \sigma)(x-y) a_y^* J(y) \Psi$$
  This identity allows them to compute expectations of the number operator and kinetic energy without needing to cancel huge normalization factors between the numerator and denominator (a common difficulty in Jastrow-based calculations).
• A-priori Bounds: They establish rigorous bounds on the local number of particles and excitations (Lemmas 4.1–4.6) to control error terms arising from the truncation of correlations and the approximation of the hard-core potential.

3. Key Contributions

1. First Rigorous Upper Bound for Hard Spheres: The paper provides the first rigorous upper bound for the ground state energy of a hard-sphere Bose gas that matches the Lee-Huang-Yang formula up to lower-order terms.
2. Correct Constant: The derived bound includes the exact coefficient $\frac{128}{15\sqrt{\pi}}$ for the $(\rho a^3)^{1/2}$ term, resolving a long-standing gap in the literature.
3. Hybrid Approach: The successful combination of a Jastrow factor (for hard-core enforcement) and a Bogoliubov transformation (for long-range correlations) in the thermodynamic limit. This demonstrates that the "soft" potential techniques can be adapted to hard spheres by separating length scales.
4. Generalization of Techniques: The methods developed, particularly the handling of the Jastrow-Bogoliubov hybrid state in the grand canonical ensemble, offer a robust framework for future studies of strongly interacting quantum gases.

4. Main Results

Theorem 1.1: For sufficiently small density $\rho > 0$, there exists a constant $C > 0$ and a small $\delta > 0$ such that the ground state energy density $e(\rho)$ satisfies:
$$e(\rho) \le 4\pi a \rho^2 \left[ 1 + \frac{128}{15\sqrt{\pi}}(\rho a^3)^{1/2} + C(\rho a^3)^{1/2+\delta} \right]$$

Combined with the existing lower bound from [20], this result proves the validity of the Lee-Huang-Yang formula for hard spheres in the dilute limit:
$$e(\rho) = 4\pi a \rho^2 \left[ 1 + \frac{128}{15\sqrt{\pi}}(\rho a^3)^{1/2} + o((\rho a^3)^{1/2}) \right]$$

5. Significance

• Mathematical Physics: This work closes a significant chapter in the rigorous mathematical theory of Bose gases. It confirms that the heuristic predictions made by Lee, Huang, and Yang in 1957 hold true even for the most singular interaction potential (hard spheres), which had been the most difficult case to treat rigorously.
• Physical Relevance: Hard spheres are a standard model for dilute atomic gases (like Helium-4 or alkali atoms in the limit of strong repulsion). Establishing the LHY correction is crucial for understanding the equation of state, stability, and properties of quantum fluids.
• Methodological Advancement: The paper demonstrates that the "Bogoliubov + Jastrow" strategy is viable in the thermodynamic limit, overcoming previous obstacles related to the non-integrability of the hard-core potential. This opens the door to proving higher-order corrections (e.g., the Wu term) for hard spheres.

In summary, this paper represents a major breakthrough in the rigorous analysis of quantum many-body systems, successfully bridging the gap between heuristic physical predictions and mathematical proof for one of the most fundamental models in condensed matter physics.