A second order upper bound to the free energy of the two dimensional Bose gas

This paper derives an explicit second-order upper bound for the free energy density of a two-dimensional Bose gas in the dilute regime below the Berezinskii-Kosterlitz-Thouless critical temperature, utilizing Bogoliubov theory to capture quasiparticle contributions with a specific dispersion relation dependent on the density and scattering length.

The Gist

Imagine a crowded dance floor filled with thousands of tiny, invisible dancers. These aren't just any dancers; they are Bose atoms (particles of a gas) that love to move in perfect unison. In physics, this is called a Bose gas.

For a long time, scientists knew exactly how these dancers behaved when the music stopped (absolute zero temperature) or when they were dancing in a huge, empty room (three dimensions). But there was a tricky scenario that remained a mystery: What happens when these dancers are in a flat, two-dimensional room (like a sheet of paper) and the room is slightly warm?

This paper, written by Florian Haberberger and Lukas Junge, solves that mystery. They figured out a precise mathematical "upper limit" for the energy of this system. Here is the story of how they did it, explained without the heavy math.

1. The Setting: A Flat Dance Floor

Imagine a giant square dance floor. The dancers are very far apart from each other (a "dilute" gas), but they still bump into each other occasionally.

• The Problem: When you heat up this floor, the dancers start moving chaotically. In the past, scientists thought that if you heated it up enough, the dancers would lose their ability to coordinate, and the "magic" of quantum mechanics (Bose-Einstein Condensation) would disappear.
• The Surprise: The authors found that even when it's warm, the dancers still manage to coordinate in small, local groups. It's like a large crowd where everyone is dancing wildly, but if you zoom in on a small patch of the floor, you see a perfect, synchronized dance happening right there.

2. The Tool: The "Jastrow Factor" (The Social Distancer)

To calculate the energy of this chaotic dance, the authors had to deal with the fact that the dancers bump into each other. In physics, these "bumps" are hard to calculate because the force gets infinitely strong when particles get too close.

They used a clever trick called a Jastrow factor.

• The Analogy: Imagine you are trying to calculate the energy of a room full of people who are constantly bumping into walls. It's a nightmare. But, what if you gave every person a personal bubble (a force field) that gently pushes them away from the walls before they actually hit them?
• The Result: This "bubble" smooths out the rough edges of the collisions. It turns a violent, impossible-to-calculate crash into a gentle, manageable nudge. This allowed the authors to simplify the math significantly.

3. The Prediction: The "Bogoliubov" Dance Step

Once they smoothed out the collisions, they applied a famous theory called Bogoliubov theory.

• The Metaphor: Think of the dancers not as individuals, but as waves on a pond. When one dancer moves, it creates a ripple that affects everyone else. These ripples are called quasiparticles.
• The authors showed that the energy of the whole system is simply the sum of the energy of the "ground state" (the calm water) plus the energy of these ripples (the waves).
• They proved that this simple "ripple" formula works even at temperatures much higher than anyone expected—right up to the point where the dance floor starts to break down (the critical temperature).

4. The "Patchwork Quilt" Strategy

How do you calculate the energy of a giant, infinite dance floor? You can't do it all at once.

• The Strategy: The authors cut the giant floor into many small, manageable squares.
• The Trick: They calculated the energy for one small square perfectly. Then, they showed that if you stitch these small squares together like a patchwork quilt, the errors cancel out, and the total energy of the giant floor is just the sum of the small squares.
• This is a bit like estimating the total weight of a million bricks by weighing one brick and multiplying by a million, but with a very sophisticated safety margin to account for the mortar between them.

5. Why This Matters

Before this paper, we had a great map of how these gases behave in 3D (like a balloon) and at absolute zero. But the "flat, warm" scenario was a blind spot.

• The Discovery: They proved that the "ripple" formula (Bogoliubov theory) is accurate even when the gas is warm and flat.
• The Impact: This gives physicists a reliable tool to predict how these quantum gases will behave in experiments. It confirms that even in a warm, chaotic environment, the underlying quantum order is still there, just hidden in small, local patches.

In a Nutshell

The authors took a chaotic, warm, flat crowd of quantum particles, gave them "personal bubbles" to smooth out their collisions, calculated the energy of a small patch, and proved that this simple calculation works for the whole system. They showed that even when things get hot and messy, the universe still follows a beautiful, predictable rhythm.

Technical Summary

1. Problem Statement

The paper addresses the rigorous derivation of the thermodynamic properties of a two-dimensional (2D) interacting Bose gas in the dilute regime (where the gas parameter $\rho a^2$ is small, with $\rho$ being the density and $a$ the scattering length).

Specifically, the authors aim to establish an upper bound for the free energy density $f(\rho, T)$ at positive temperatures $T$. The goal is to prove that the free energy is accurately described by the Bogoliubov approximation up to the Berezinskii–Kosterlitz–Thouless (BKT) critical temperature ($T_c$), a regime where Bose-Einstein Condensation (BEC) does not occur in the thermodynamic limit due to the Mermin-Wagner theorem, yet local condensation persists on short scales.

The target formula for the free energy density is:
$$f(\rho, T) \approx e(\rho) + \frac{T}{(2\pi)^2} \int_{\mathbb{R}^2} \log\left(1 - e^{-\frac{1}{T}\sqrt{p^4 + 8\pi\rho\delta p^2}}\right) dp$$
where $e(\rho)$ is the ground state energy density and $\delta$ is a diluteness parameter related to the scattering length.

2. Methodology

The authors employ a variational approach based on the Gibbs variational principle:
$$f(\rho, T) = \inf_{\Gamma} \left( \frac{\text{Tr}(H\Gamma) - T S(\Gamma)}{|\Omega|} \right)$$
To derive the upper bound, they construct a specific trial state $\Gamma_0$ that minimizes the free energy functional. The construction involves several sophisticated steps:

A. Localization and Grand Canonical Setup

• The thermodynamic box $\Omega$ is partitioned into smaller squares of side length $L$.
• On each small box, a grand canonical trial state is constructed with periodic boundary conditions.
• These local states are later patched together to form a global state with Dirichlet boundary conditions on the full thermodynamic box.

B. Construction of the Trial State

The trial state $\Gamma_0$ is built via a sequence of unitary transformations applied to a reference Gibbs state $\Gamma_D$:

1. Jastrow Factor ($F$): A correlation factor $F = \prod f(x_i - x_j)$ is introduced to handle the short-range singularities of the interaction potential $v$. This "softens" the potential, replacing it with a regular effective potential $\tilde{v}$ while preserving the scattering length.
2. Weyl Transformation ($W$): Introduces a condensate of $N_0(1+\eta)$ particles in the zero-momentum mode to account for the macroscopic occupation expected at the chosen length scale $L$.
3. Bogoliubov Transformation ($e^B$): A quadratic unitary transformation that diagonalizes the quadratic part of the Hamiltonian, incorporating correlations between excited modes.

The resulting state is:
$$\Gamma_0 = \frac{F W e^B \Gamma_D e^{-B} W^* F}{\text{Tr}(F W e^B \Gamma_D e^{-B} W^* F)}$$

C. Key Technical Innovations

• Entropy Control: Unlike previous works (e.g., for ground state energy), handling positive temperature requires controlling the entropy change induced by the Jastrow factor. The authors use a Fannes-type inequality (Lemma 12) to bound the difference in entropy between the trial state and the uncorrelated state, relying on the trace norm distance.
• Length Scale Selection: The side length $L$ of the localization boxes is chosen as $L = \rho^{-1/2} Y^{-\alpha}$ (where $Y \sim |\log(\rho a^2)|^{-1}$). This scale is crucial as it corresponds to the scale where local BEC is expected to occur even above the global $T_c$.
• Soft Potential Analysis: The authors rigorously analyze the effective potential generated by the Jastrow factor, showing it behaves like a delta function with a specific strength related to the scattering length.

3. Key Contributions and Results

Main Theorem (Theorem 1)

The authors prove that for a 2D Bose gas with a positive, radial, compactly supported potential, there exists a constant $C$ such that for temperatures $T \leq c T_c$ and small diluteness:
$$f(\rho, T) \leq 2\pi\rho^2\delta \left(1 + \left(\gamma + \frac{1}{4} + \frac{\log \pi}{2}\right)\delta \right) + \frac{T}{(2\pi)^2} \int_{\mathbb{R}^2} \log\left(1 - e^{-\frac{1}{T}\sqrt{p^4 + 8\pi\rho\delta p^2}}\right) dp + \text{Error Terms}$$
where:

• $\delta = \frac{2}{|\log(\rho a^2)| + \log|\log(\rho a^2)|}$.
• $\gamma$ is the Euler-Mascheroni constant.
• The error term is of order $O(\rho^2 \delta (\delta |\log \delta| + T/T_c)^2)$.

Significance of the Result

1. Validity up to $T_c$: The result is valid for temperatures comparable to the BKT critical temperature $T_c = \frac{4\pi\rho}{|\log \delta|}$. Previous rigorous results in 3D were restricted to much lower temperatures ($T \sim \rho a$).
2. Bogoliubov Spectrum: The derivation confirms that the excitation spectrum remains of the Bogoliubov form $\sqrt{p^4 + 8\pi\rho\delta p^2}$ even in the absence of global BEC, validating the physical intuition that the free energy is a local quantity.
3. Second-Order Precision: The result captures the second-order correction to the ground state energy (the Mora-Castin term) and the thermal contribution simultaneously.
4. 2D Specificity: The proof is not a simple adaptation of 3D methods. The logarithmic behavior of the 2D scattering length requires a distinct analysis, particularly regarding the choice of the localization scale and the handling of the Jastrow factor.

4. Comparison with Existing Literature

• Ground State ($T=0$): The ground state energy term matches the rigorous results by Lieb-Yngvason (leading order) and Fournais et al. (second order).
• Finite Temperature: Previous works (e.g., Mayer [24]) provided bounds for higher temperatures but with lower precision or restricted to the ideal gas regime. This paper bridges the gap by providing a rigorous upper bound that includes interaction effects up to the critical temperature.
• 3D vs. 2D: In 3D, the temperature range for such rigorous bounds is currently limited to $T \ll T_c$. This work extends the rigorous validity of the Bogoliubov approximation to the critical regime in 2D.

5. Conclusion

Haberberger and Junge provide a rigorous mathematical proof that the Bogoliubov approximation correctly describes the free energy density of a 2D dilute Bose gas up to the BKT critical temperature. By constructing a sophisticated trial state involving Jastrow correlations, Weyl shifts, and Bogoliubov rotations, and by carefully controlling the entropy and localization errors, they establish an upper bound that captures both the universal ground-state corrections and the thermal quasiparticle contributions. This work resolves a central open problem in mathematical physics regarding the thermodynamics of 2D Bose gases in the interacting regime.