🔢 mathematics

Phonon Dispersion Relation in Weakly Coupled Bose Gas Beyond Bogoliubov Approximation

This paper investigates the quasiparticle dispersion relation of a weakly coupled homogeneous Bose gas beyond the Bogoliubov approximation, demonstrating that first-order corrections in both canonical and grand-canonical settings preserve a linear phonon behavior at small momenta with a computable velocity.

Original authors: Jan Dereziński, Lorenzo Pettinari

Published 2026-10-21

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Jan Dereziński, Lorenzo Pettinari

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the coldest corners of the universe, where temperatures drop so low that atoms lose their individual identities, a strange new state of matter emerges. This is the Bose gas, a collection of particles that, instead of bouncing off one another like billiard balls, begin to move in perfect unison, behaving as a single, giant quantum wave. When this happens, the gas does not flow like a normal fluid; it becomes a superfluid, capable of moving without any friction or resistance. To understand how this superfluid behaves, physicists look at the tiny ripples that travel through it. These ripples are called phonons, and they act like sound waves moving through a solid, but on a quantum scale. For decades, a famous theory has provided a good first guess at how fast these ripples travel, but that theory is an approximation. It ignores the subtle, complex ways the particles interact with each other, treating them as if they were mostly independent. While this simple picture works well for a rough estimate, it fails to capture the precise details needed to understand the true nature of these quantum fluids, especially when scientists want to predict exactly how fast sound travels through them.

A team of researchers has now taken a significant step beyond that old approximation. They set out to calculate the precise speed of these quantum ripples, known as phonons, by accounting for the interactions between particles that the simpler theory left out. To do this, they built a more sophisticated mathematical model of the gas, one that carefully tracks how the particles influence each other in a crowded, cold environment. They approached the problem in two different ways: one method where the total number of particles is fixed, and another where the number of particles can fluctuate, similar to how a gas might exchange particles with a surrounding reservoir. By working through these complex calculations, they were able to find the first non-trivial corrections to the speed of sound in this gas. Their work reveals that while the speed of sound is indeed linear at very low energies—meaning the ripples travel at a constant speed—there is a tiny, specific adjustment to that speed caused by the particles bumping into one another.

The researchers discovered that this adjustment depends on the density of the gas and the strength of the interaction between the particles. They found that for the slowest, most gentle ripples, the speed of sound is not just a simple number but includes a small correction term that arises from the collective behavior of the entire system. This correction is surprisingly clean and well-defined, allowing them to compute a precise value for the speed of sound that includes these subtle quantum effects. Their calculations showed that the two different methods they used—one with a fixed number of particles and one with a fluctuating number—ultimately agreed with each other, confirming that their result is robust and not an artifact of a specific mathematical trick. This agreement is crucial because it suggests that the physical reality of the sound speed is independent of how we choose to count the particles in our theoretical model.

One of the most important aspects of their work is how they handled the mathematical difficulties that usually arise when dealing with infinite systems. In the real world, a gas is contained in a box, but to understand its fundamental properties, physicists often imagine the box becoming infinitely large. In this limit, certain mathematical terms can blow up or become undefined, making it impossible to get a clear answer. The researchers developed a rigorous way to separate the useful parts of the calculation from the parts that vanish or become irrelevant as the system grows. They proved that the messy, unphysical parts of their model disappear in the limit of a large system, leaving behind a clean, finite result for the speed of sound. This careful handling of the "tail" of the calculation ensures that their prediction for the speed of sound is mathematically sound and physically meaningful.

The paper also connects their findings to a different way of looking at the problem, known as the dilute limit, which is often used to describe gases that are very spread out. In this regime, the gas is so thin that the particles rarely interact, and the behavior is governed by a single number called the scattering length, which measures how likely two particles are to bounce off each other. The researchers showed that their new, more precise formula for the speed of sound matches the famous results derived by earlier physicists in this dilute limit, but with an added layer of detail. They found that the speed of sound includes a small correction that depends on the density of the gas raised to a specific power. This correction is a direct consequence of the interactions between the particles, and it confirms that even in a very thin gas, the collective behavior of the atoms plays a vital role in determining how fast sound travels.

By providing a rigorous derivation of these corrections, the researchers have offered a clearer picture of the quantum world. They have shown that the speed of sound in a Bose gas is not a fixed, unchangeable constant but a value that shifts slightly depending on how the particles interact. This shift is small, but it is real and measurable. Their work bridges the gap between the simple, idealized models of the past and the complex, interacting reality of quantum fluids. It demonstrates that even in a system where particles seem to move in perfect harmony, the subtle push and pull between them leaves a distinct fingerprint on the way energy travels through the medium. This level of precision is essential for understanding the fundamental properties of superfluids and for interpreting experiments that probe the quantum nature of matter.

The researchers did not stop at just finding the correction; they also verified that their results hold up under different conditions and assumptions. They checked that their findings were consistent with the known behavior of the system at very low momenta, where the ripples are long and gentle. They confirmed that the speed of sound remains linear, meaning the relationship between the energy of the ripple and its momentum stays a straight line, but with a slightly different slope than the simple theory predicted. This slope is the speed of sound, and their calculation provides the exact value of that slope, including the tiny adjustments caused by particle interactions. The work stands as a testament to the power of careful mathematical analysis in revealing the hidden details of the quantum world, showing that even the most subtle effects can be understood and predicted with the right tools.

In the end, this paper provides a definitive answer to a question that has lingered since the early days of quantum mechanics: how exactly do particles in a superfluid affect the speed of sound? The answer is that they do affect it, and the effect is calculable. The researchers have shown that the speed of sound is determined by a combination of the density of the gas and the strength of the interactions between the particles, with a specific correction term that arises from the complex dance of quantum mechanics. This correction is not just a theoretical curiosity; it is a fundamental property of the system that can be tested in experiments. By providing a rigorous and precise calculation, the researchers have opened the door to a deeper understanding of quantum fluids, bringing us one step closer to mastering the strange and wonderful behavior of matter at its coldest.

Technical Summary: Phonon Dispersion Relation in Weakly Coupled Bose Gas Beyond Bogoliubov Approximation

Problem Statement
The paper addresses the calculation of the quasiparticle dispersion relation for a homogeneous Bose gas at zero temperature, specifically seeking corrections beyond the standard Bogoliubov approximation. While the Bogoliubov theory successfully predicts a linear dispersion relation for low-momentum phonons (sound waves) and the existence of maxons and rotons at higher momenta, it is an approximation that discards specific interaction terms and relies on a parameter $\nu$ (related to density and chemical potential) whose precise fixing is non-trivial beyond the leading order.

The authors aim to compute the leading non-trivial corrections to the dispersion relation in the weak-coupling regime. A primary challenge identified is the definition of the dispersion relation itself in a finite volume, where it is ill-defined, versus the thermodynamic limit where it emerges. Furthermore, in the lowest non-trivial order of perturbation theory, the dispersion shift acquires a small imaginary part, indicating quasiparticle instability. This work focuses specifically on the real part of the energy shift to determine the corrected phonon velocity.

Methodology
The authors employ a rigorous perturbative approach in both the canonical (fixed particle number $N$ ) and grand-canonical (fixed chemical potential $\mu$ ) ensembles.

Hamiltonian Rearrangement: Instead of using the traditional "c-number substitution" (replacing the zero-mode operator $a_0$ with $\sqrt{N_0}$ ), the authors reorganize the Hamiltonian. They introduce an "extended space" formalism (inspired by Arnowitt and Girardeau) where the zero-mode occupation number is allowed to take negative values, embedding the physical $N$ -particle space into a larger Fock space. This allows for a consistent second-quantized treatment without ad-hoc substitutions.
Perturbative Expansion: The interaction potential is rescaled as $\hat{v} \to \kappa \hat{v}$ with $\kappa$ small. The Hamiltonian is decomposed into a quadratic Bogoliubov part ( $H_{bg}$ ) and interaction terms of various orders in $\kappa$ ( $\sqrt{\kappa}H_3$ , $\kappa H_2$ , $\kappa H_4$ , etc.).
Dispersion Shift Calculation: The energy shift for a one-quasiparticle state $b^*_k \Omega_{bg}$ $b_{k}^{*} Ω_{b g}$ is computed using Rayleigh–Schrödinger perturbation theory. The shift $\delta(k)$ $δ (k)$ is the difference between the excited state energy shift and the ground state energy shift. It consists of two main contributions:
- Hellmann-Feynman term ( $\delta_4$ ): The expectation value of the quartic interaction terms.
- Fermi Golden Rule term ( $\delta_3$ ): The second-order contribution arising from the cubic interaction terms.
Infrared Regularization: A critical technical aspect is the handling of infrared singularities. Both $\delta_3$ and $\delta_4$ individually exhibit infrared divergences (behaving as $O(k^{-1})$ ) in the limit of small momentum $k$ . The authors demonstrate that these singularities cancel exactly when the terms are summed, yielding a finite, linear correction to the dispersion relation.
Thermodynamic Limit: The authors rigorously justify discarding "tail terms" (remainders arising from the extended space formalism) which vanish exponentially in the thermodynamic limit ( $L \to \infty$ ) due to the localization of the Bogoliubov vacuum.

Key Contributions and Results

Canonical Formalism (Theorem 1.1): The authors derive the dispersion relation up to order $\kappa$ in the canonical setting. They prove that for small momenta, the corrected dispersion relation is linear:
$\omega(k) \approx |k| c_{can}(\nu, \hat{v})$
where the corrected sound velocity $c_{can}$ is given by:
$c_{can} = \sqrt{\nu} - \kappa \hat{v}(0)\sqrt{\nu} \frac{1}{16} \int \frac{dp}{(2\pi)^3} \frac{\hat{v}(p)^2}{\hat{v}(0)^2} \frac{|p|^4}{\omega_{bg}(p)^3}$
The calculation explicitly shows the cancellation of infrared divergences between the Hellmann-Feynman and Fermi Golden Rule terms.
Grand-Canonical Formalism (Theorem 1.2): A parallel derivation is performed in the grand-canonical setting. While the intermediate steps differ (specifically regarding the treatment of the chemical potential $\mu$ and the zero-mode displacement), the final result for the phonon velocity agrees with the canonical result up to order $\kappa$ :
$c_{can}(\nu, \hat{v}) = c_{gcan}(\mu, \hat{v})(1 + o(\kappa))$
This establishes the equivalence of the two ensembles for this specific observable (dispersion relation) at this order of perturbation, a non-trivial result given that dispersion relations are finer objects than extensive thermodynamic quantities.
Connection to Dilute Limit (Corollary 1.3): By relating the weak coupling parameter $\kappa$ to the scattering length $a$ via the Born series, the authors recast their result in the dilute limit ( $na^3 \to 0$ ). They recover the well-known Beliaev result for the phonon velocity:
$c \approx \sqrt{4\pi n_0 a} \left( 1 + \frac{7}{6\pi^2} \sqrt{n_0 (4\pi a)^3} \right)$
This confirms that their perturbative expansion in the weak-coupling regime is consistent with established results in the dilute limit.

Significance and Claims
The paper claims to provide a mathematically rigorous derivation of the leading-order corrections to the phonon dispersion relation in a weakly coupled Bose gas. Its significance lies in:

Rigorous Justification: Unlike many physics derivations that rely on diagrammatic heuristics or Green's functions without full control of the thermodynamic limit, this work provides a rigorous operator-theoretic proof, including estimates for the remainders and justification for discarding specific terms.
Resolution of Infrared Issues: The work explicitly demonstrates how the infrared divergences, which appear in standard perturbative calculations, cancel out in the sum of the Hellmann-Feynman and Fermi Golden Rule contributions, ensuring a well-defined phonon velocity.
Ensemble Equivalence: It provides a concrete example where the canonical and grand-canonical approaches yield identical results for a non-extensive quantity (the dispersion relation) beyond the leading order, clarifying the relationship between the parameters $\nu$ , $\mu$ , and $\mu_0$ .
Methodological Innovation: The use of the "extended space" formalism to avoid direct c-number substitution offers a cleaner mathematical framework for handling the zero-mode in the presence of interactions, which the authors argue is superior to traditional approaches for this specific problem.

The authors note that while the weak-coupling approximation is mathematically clean, it may not be the best physical approximation for systems with hard-core interactions (like liquid helium), where the dilute limit is more standard. However, their method successfully bridges the gap, recovering the dilute limit results from the weak-coupling expansion.

More like this