A Gapless Theory of Bose-Einstein Condensation in Dilute Gases at Finite Temperature

TL;DR

This paper presents a comprehensive gapless theoretical framework for Bose-Einstein condensation in dilute gases at finite temperature, addressing both low and high energy divergences and providing a systematic approach for excitation analysis.

Contribution

It introduces a consistent gapless theory applicable to trapped and homogeneous gases, resolving issues with infra-red and ultra-violet divergences and improving upon previous approximations.

Findings

The excitation spectrum is confirmed to be gapless.

The theory remains finite beyond quadratic order.

Perturbation theory is justified away from the critical region.

Abstract

In this paper we develop a gapless theory of BEC which can be applied to both trapped and homogeneous gases at zero and finite temperature. The many-body Hamiltonian for the system is written in a form which is approximately quadratic with higher order cubic and quartic terms. The quadratic part is diagonalized exactly by transforming to a quasiparticle basis, while the non-quadratic terms are dealt with using first and second order perturbation theory. The conventional treatment of these terms, based on factorization approximations, is shown to be inconsistent. Infra-red divergences can appear in individual terms of the perturbation expansion, but we show analytically that the total contribution beyond quadratic order is finite. The resulting excitation spectrum is gapless and the energy shifts are small for a dilute gas away from the critical region, justifying the use of perturbation theory. Ultra-violet divergences can appear if a contact potential is used to describe particle interactions. We show that the use of this potential as an approximation to the two-body T-matrix leads naturally to a high-energy renormalization. The theory developed in this paper is therefore well-defined at both low and high energy and provides a systematic description of Bose-Einstein condensation in dilute gases. It can therefore be used to calculate the energies and decay rates of the excitations of the system at temperatures approaching the phase transition.