Bogoliubov Approximation for Exciton Gas Quasiparticle Spectrum

Hello! In our last two lessons, we explored the fascinating phenomena of Bose-Einstein condensation in both 3D and 2D. We saw that in 3D, a BEC forms below a critical temperature through the macroscopic occupation of the ground state. In 2D, while true long-range order is forbidden, an interacting exciton gas can still undergo a BKT transition into a superfluid state characterized by quasi-long-range order.

A common thread in these discussions is the crucial role of inter-particle interactions. The non-interacting ideal gas model is a useful starting point, but it fails to capture essential physics, such as the collective dynamics of the condensate. Today, we will develop the fundamental tool to properly describe these dynamics.

Our learning outcome is to apply the Bogoliubov approximation to derive the quasiparticle spectrum of a weakly interacting exciton gas. This will reveal that the elementary excitations in a BEC are not single particles but collective modes called quasiparticles, which have a unique energy-momentum relationship.

1. The Problem: An Interacting Bose Gas

Let's begin with the Hamiltonian for a gas of interacting bosons (excitons), which you've seen before in the context of our many-body formulation. For simplicity, we'll assume a weak, short-range (contact) repulsive interaction, which is a good model for excitons at densities where they don't form biexcitons. In second quantization, the Hamiltonian is:
$H^= \sum_{k} \frac{ℏ ^{2} k ^{2}}{2 M} a^_{k}^{†} a^_{k} + \frac{g}{2 V} \sum_{k, k^{'}, q} a^_{k + q}^{†} a^_{k^{'} - q}^{†} a^_{k^{'}} a^_{k}$
Here, $M$ is the exciton effective mass, $V$ is the volume, and $g$ is a constant representing the interaction strength.

The four-operator interaction term makes this Hamiltonian impossible to solve exactly. However, for a condensate at low temperatures, we can make a powerful approximation.

2. The Bogoliubov Approximation: Taming the Hamiltonian

In a BEC, the zero-momentum ground state is macroscopically occupied. This means the number of particles in the $k = 0$ state, $N_{0}$ , is enormous ( $N_{0} \approx N$ , the total number of particles). The key insight, developed by Nikolai Bogoliubov, is that the creation and annihilation operators for this ground state, $a^_{0}^{†}$ and $a^_{0}$ , behave more like classical numbers than quantum operators.

Why? Their commutator, $[a^_{0}, a^_{0}^{†}] = 1$ , is negligible compared to the operators themselves, whose matrix elements scale as $N_{0}$ . The situation is analogous to treating a strong electromagnetic field classically, as a c-number, rather than a quantum field of photons.

The Bogoliubov approximation consists of replacing the ground state operators with the c-number $N_{0}$ :
$a^_{0} \to N_{0}, a^_{0}^{†} \to N_{0}$
We then systematically expand the Hamiltonian, keeping only terms up to second order in the operators for the excited states ( $a^_{k}$ and $a^_{k}^{†}$ for $k \neq = 0$ ), since we assume the number of excited particles is small.

19. Bose gases

Before diving into the algebra, watch this short clip from an MIT OpenCourseWare lecture by Wolfgang Ketterle. It gives an excellent conceptual overview of the Bogoliubov approximation and the resulting quasiparticle spectrum we are about to derive.

Watch the section from 02:14 to 04:30. Focus on how replacing the condensate operator with a c-number transforms the difficult four-operator Hamiltonian into a solvable bilinear one, and note the characteristic shape of the resulting energy dispersion.

As the video explained, this approximation transforms the problem into a solvable quadratic Hamiltonian. Now, let's work through the derivation.

3. Deriving the Quasiparticle Spectrum

By applying the Bogoliubov approximation to our interaction Hamiltonian and retaining only the terms that are constant or quadratic in the excited-state operators, we arrive at a simplified, but still non-diagonal, Hamiltonian. The key challenge lies in the terms that create or annihilate pairs of excited particles from the condensate, of the form $a^_{k}^{†} a^_{- k}^{†}$ and $a^_{k} a^_{- k}$ .

To diagonalize this Hamiltonian, we introduce a canonical transformation—the Bogoliubov transformation—to a new set of bosonic operators, $b^_{k}$ and $b^_{k}^{†}$ . These new operators will represent the true elementary excitations of the interacting system: the quasiparticles.

The derivation is a cornerstone of many-body physics. The following text provides a clear, step-by-step mathematical treatment.

Bogoliubov theory of the weakly interacting Bose gas

The following text, 'Bogoliubov theory of the weakly interacting Bose gas', will guide us through the full derivation.

Please read section 4.3, starting from Eq. (4.17). Follow the derivation carefully. Your goal is to understand how the Hamiltonian is first approximated (Eq. 4.19), and then diagonalized using the Bogoliubov transformation (Eq. 4.20) to arrive at the final quasiparticle energy spectrum, \epsilon(p), given in Eq. (4.27). Pay close attention to the constraint u_p^2 - v_{-p}^2 = 1 (Eq. 4.22), which ensures the new operators remain bosons.

4. Analysis of the Bogoliubov Spectrum

Welcome back. You have now derived one of the most important results in the theory of superfluids. The energy of the quasiparticle excitations is given by:
$E_{k} = (\frac{ℏ ^{2} k ^{2}}{2 M})^{2} + 2 n g (\frac{ℏ ^{2} k ^{2}}{2 M}) = ϵ_{k} (ϵ_{k} + 2 n g)$
where $ϵ_{k} = \frac{ℏ ^{2} k ^{2}}{2 M}$ is the kinetic energy of a free exciton and $n = N / V$ is the total exciton density.

This dispersion relation contains a wealth of physics. Let's analyze its two limits, as discussed in the text you just read (see Section 4.4.3 and Figure 4.4(a) in the resource Bogoliubov theory of the weakly interacting Bose gas).

Low-Momentum Limit ( $k \to 0$ )

For small momenta, the kinetic energy term $ϵ_{k}$ is very small. We can expand the expression under the square root:
$E_{k} \approx 2 n g ϵ_{k} = 2 n g \frac{ℏ ^{2} k ^{2}}{2 M} = \frac{n g}{M} ℏ k$
This is a linear dispersion relation: $E_{k} = c_{s} ℏ k$ .
This is profound. The low-energy excitations of an interacting BEC are not particle-like ( $\propto k^{2}$ ) but are collective density waves, i.e., sound waves (phonons). The speed of these sound waves is $c_{s} = n g / M$ . This linear "Goldstone mode" is a generic feature of systems with a spontaneously broken continuous symmetry (here, the U(1) phase symmetry of the condensate).

High-Momentum Limit ( $k \to \infty$ )

For large momenta, the kinetic energy term $ϵ_{k}$ dominates the interaction term $2 n g$ .
$E_{k} = ϵ_{k}^{2} + 2 n g ϵ_{k} = ϵ_{k} 1 + \frac{2 n g}{ϵ _{k}} \approx ϵ_{k} (1 + \frac{n g}{ϵ _{k}}) = \frac{ℏ ^{2} k ^{2}}{2 M} + n g$
Here, the dispersion relation becomes quadratic again, just like a free particle. The excitation is essentially a single exciton moving so fast that it barely feels the collective effects of the condensate. The energy is simply its kinetic energy plus a constant mean-field shift ( $n g$ ) from the interaction with all other excitons.

The transition between these two regimes occurs when the kinetic energy is comparable to the interaction energy, $ϵ_{k} \approx n g$ . This defines a characteristic length scale, the healing length $ξ = ℏ/ 2 M n g$ , which is the scale over which the condensate wavefunction "heals" back to its bulk value after a local perturbation. For $k ≪ 1/ ξ$ , excitations are phonon-like; for $k ≫ 1/ ξ$ , they are particle-like.

Test your understanding!

The Bogoliubov transformation mixes creation and annihilation operators ( $a^_{k} = u_{k} b^_{k} + v_{k} b^_{- k}^{†}$ ). What does this imply about the nature of the ground state of the interacting system (the "Bogoliubov vacuum", defined by $b^_{k} ∣ G ⟩ = 0$ for all $k$ ) compared to the ground state of the non-interacting system (the "bare vacuum", where all particles are in the $k = 0$ state)?

Show answer

The transformation implies that the interacting ground state $∣ G ⟩$ is not empty of excitations with respect to the non-interacting picture. If we calculate the number of "bare" particles with momentum $k$ in the interacting ground state, we find:
$⟨ G ∣ a^_{k}^{†} a^_{k} ∣ G ⟩ = ⟨ G ∣ (u_{k} b^_{k}^{†} + v_{k} b^_{- k}) (u_{k} b^_{k} + v_{k} b^_{- k}^{†}) ∣ G ⟩ = v_{k}^{2}$
Since $v_{k} \neq = 0$ , this means the interacting ground state contains a finite number of correlated pairs of particles with opposite momenta ( $k$ and $- k$ ) that have been scattered out of the condensate by the interaction. The Bogoliubov vacuum is a sea of these virtual pairs, a much richer structure than the simple, non-interacting ground state.

5. Experimental Verification

This theoretical framework is not just a mathematical exercise; the Bogoliubov spectrum has been precisely measured in numerous systems, from ultracold atomic gases to the exciton-polariton condensates you are interested in.

How is this done? A common technique is momentum-resolved spectroscopy, such as Bragg scattering, where two laser beams impart momentum $ℏ q$ and energy $ℏ ω$ to the condensate. A strong response is detected only when $ℏ ω$ and $ℏ q$ match the dispersion relation of an allowed excitation, $E_{q}$ . This allows experimentalists to directly map out the $E_{k}$ vs. $k$ curve.

19. Bose gases

Let's return to the MIT lecture. Prof. Ketterle now explains how these quasiparticle excitations are observed experimentally.

Watch the segment from 09:48 to 13:11. This will connect the theoretical dispersion curve to what is measured in the lab, reinforcing the physical meaning of our derived spectrum.

The experimental data, such as that shown in Figure 4.5 of the Bogoliubov theory PDF, show a remarkable agreement with the predicted dispersion, beautifully confirming the validity of Bogoliubov's theory.

Conclusion

In this lesson, we have tackled the fundamental problem of a weakly interacting Bose-Einstein condensate. By treating the condensate as a classical object and diagonalizing the remaining quantum fluctuations, we have uncovered the nature of its elementary excitations.

Key Takeaways:

The Bogoliubov approximation ( $a^_{0}, a^_{0}^{†} \to N_{0}$ ) allows us to solve the Hamiltonian of a weakly interacting Bose gas by reducing it to a quadratic form.
The Hamiltonian is diagonalized by a Bogoliubov transformation, which defines new operators for collective excitations called quasiparticles.
The derived quasiparticle energy spectrum is $E_{k} = ϵ_{k} (ϵ_{k} + 2 n g)$ .
This spectrum is linear (phonon-like, $E_{k} \propto k$ ) at low momenta, reflecting the collective, sound-wave nature of long-wavelength excitations.
The spectrum becomes quadratic (free-particle-like, $E_{k} \propto k^{2}$ ) at high momenta, where excitations behave as individual particles.

Preview of the Next Lesson

We have now established the key theoretical properties of an exciton BEC: the critical temperature for its formation, the special case of the BKT transition in 2D, and the spectrum of its elementary excitations. The next logical step is to ask: how would we definitively identify an exciton BEC in an experiment?

In the next lesson, we will focus on the learning outcome: Evaluate experimental signatures of exciton BEC, such as macroscopic coherence, narrowed emission lines, and the formation of a spectral condensate. We will see how the theoretical concepts we've developed manifest in measurable optical and transport properties.

Recommended resources