Abstract We derive and analyze the coupled equations of quadratic approximation of the Bogoliubov model for a weakly interacting Bose gas. The first equation determines the condensate density as a variational parameter and ensures the minimum of the grand thermodynamic potential. The second one provides a relation between the total number of particles and chemical potential. Their consistent theoretical analysis is performed for a number of model interaction potentials including contact (local) and nonlocal interactions, where the latter provide nontrivial dependencies in momentum space. We demonstrate that the derived equations have no solutions for the local potential, although they formally reproduce the well-known results of the Bogoliubov approach. At the same time, it is shown that these equations have the solutions for physically relevant nonlocal potentials. We show that in the regimes close to experimental realizations with ultracold atoms, the contribution of the terms originating from the quadratic part of the truncated Hamiltonian to the chemical potential can be of the same order of magnitude as from its c -number part. Due to this fact, in particular, the spectrum of single-particle excitations in the quadratic approximation acquires a gap. The issue of the gap is also discussed.
Abstract We derive and analyze the coupled equations of quadratic approximation of the Bogoliubov model for a weakly interacting Bose gas. The first equation determines the condensate density as a variational parameter and ensures the minimum of the grand thermodynamic potential. The second one provides a relation between the total number of particles and chemical potential. Their consistent theoretical analysis is performed for a number of model interaction potentials including contact (local) and nonlocal interactions, where the latter provide nontrivial dependencies in momentum space. We demonstrate that the derived equations have no solutions for the local potential, although they formally reproduce the well-known results of the Bogoliubov approach. At the same time, it is shown that these equations have the solutions for physically relevant nonlocal potentials. We show that in the regimes close to experimental realizations with ultracold atoms, the contribution of the terms originating from the quadratic part of the truncated Hamiltonian to the chemical potential can be of the same order of magnitude as from its c -number part. Due to this fact, in particular, the spectrum of single-particle excitations in the quadratic approximation acquires a gap. The issue of the gap is also discussed.
We discuss the standard approach to the problem of the low momentum limit of the spectrum for a weakly interacting Bose gas. The Bogoliubov's spectrum is shown to be obtained as a Goldstone mode thanks to the introduction of a chemical potential μ. This procedure has, however, difficulties since the breaking of the gauge symmetry implies that the corresponding chemical potential must be taken as zero, unless it is introduced before breaking the symmetry. But if this is done, after the symmetry breaking μ loses its meaning as a chemical potential. An alternative two-mode solution is suggested having two modes, one of them being the free-particle quadratic in momentum spectrum, the second bearing a gap. This gap leads to a λ-type behavior of the specific heat near the critical temperature.
We re-examine the way in which Bogoliubov's theory of a dilute Bose gas at T=0 has been extended to describe the statistical mechanics of interacting bosons at finite temperature. We show explicitly that the field-theoretic calculation of the grand partition function in this formulation amounts to a canonical trace over the eigenfunctions of the Bogoliubov Hamiltonian at fixed total number of bosons N, and that the additional trace over N that is required in the grand-canonical formalism is never carried out. This implies that what usually passes as the grand-canonical treatment of the Bogoliubov Hamiltonian is not quite grand-canonical, and is in fact a canonical one. We also show that the discontinuity in the condensate density predicted by previous formulations of this theory as the temperature T goes past the critical transition temperature Tc is a direct consequence of an inappropriate generalization of the Bogoliubov prescription to finite temperatures, and that this discontinuity disappears when this prescription is either used as a zero temperature approximation or avoided altogether. Armed with the above findings, we reformulate the statistical mechanics of interacting bosons in the canonical ensemble and derive the thermodynamics of the system. We then show how the canonical treatment can be used to setup a truly grand-canonical description of the statistical mechanics of a weakly interacting Bose gas where the average number of bosons in the system varies with temperature, unlike existing formulations where the total number of bosons N is taken to be a constant that does not depend on T. Consequences on the physics of interacting bosons are briefly discussed.
As is well-known, in the conventional formulation of Bogoliubov's theory of an interacting Bose gas, the Hamiltonian is written as a decoupled sum of contributions from different momenta of the form . Then, each of the single-mode Hamiltonians is diagonalized separately, and the resulting ground state wavefunction of the total Hamiltonian is written as a simple product of the ground state wavefunctions of each of the single-mode Hamiltonians . While this way of diagonalizing the total Hamiltonian may seem to be valid from the perspective of the standard, number non-conserving Bogoliubov's method, where the k = 0 state is removed from the Hilbert space and hence the individual Hilbert spaces where the Hamiltonians {} are diagonalized are disjoint from one another, we argue that from a number-conserving perspective this diagonalization method may not be adequate since the true Hilbert spaces where the Hamiltonians {} should be diagonalized all have the k = 0 state in common, and hence the ground state wavefunction of the total Hamiltonian may not be written as a simple product of the ground state wavefunctions of the 's. In this paper, we give a thorough review of Bogoliubov's method, and discuss a variational and number-conserving formulation of this theory in which the k = 0 state is restored to the Hilbert space of the interacting gas, and where, instead of diagonalizing the Hamiltonians separately, we diagonalize the total Hamiltonian . When this is done, we find that the ground state energy is lowered below the Bogoliubov result, and the depletion of bosons is significantly reduced with respect to the one obtained in the number non-conserving treatment. We also find that the spectrum of the usual αk excitations of Bogoliubov's method changes from a gapless one, as predicted by the standard, number non-conserving formulation of this theory, to one which exhibits a finite gap in the k → 0 limit. We discuss the presence of a gap in the spectrum of the αk's in light of Goldstone's theorem, and show that there is no contradiction with the latter.
We analyze the equilibrium properties of a weakly interacting, trapped quasi-one-dimensional Bose gas at finite temperatures and compare different theoretical approaches. We focus in particular on two stochastic theories: a number-conserving Bogoliubov (NCB) approach and a stochastic Gross-Pitaevskii equation (SGPE) that have been extensively used in numerical simulations. Equilibrium properties like density profiles, correlation functions, and the condensate statistics are compared to predictions based upon a number of alternative theories. We find that due to thermal phase fluctuations, and the corresponding condensate depletion, the NCB approach loses its validity at relatively low temperatures. This can be attributed to the change in the Bogoliubov spectrum, as the condensate gets thermally depleted, and to large fluctuations beyond perturbation theory. Although the two stochastic theories are built on different thermodynamic ensembles (NCB, canonical; SGPE, grand-canonical), they yield the correct condensate statistics in a large Bose-Einstein condensate (BEC) (strong enough particle interactions). For smaller systems, the SGPE results are prone to anomalously large number fluctuations, well known for the grand-canonical, ideal Bose gas. Based on the comparison of the above theories to the modified Popov approach, we propose a simple procedure for approximately extracting the Penrose-Onsager condensate from first- and second-order correlation functions that is both computationally convenient and of potential use to experimentalists. This also clarifies the link between condensate and quasicondensate in the Popov theory of low-dimensional systems.
In Chapter 6, we derived a generalized Gross–Pitaevskii condensate equation which is coupled to a kinetic equation for the distribution function for the thermal atoms. However, the kinetic equation in Chapter 6 is only valid in the semiclassical limit. It involves the assumption that the thermal energy k B T is much greater than the spacing between the harmonic trap energy levels ( k B T ≫ω 0 where ω 0 is the trap frequency) and also much greater than the average interaction energy ( k B T ≫ gn ). The ZNG model, based on HF excitations, is still expected to be adequate down to quite low temperatures in trapped Bose gases, as will be shown by the results in Chapter 12. However, the ZNG model will break down at very low temperatures, where the Hartree–Fock excitations must be replaced by the Bogoliubov spectrum. To deal with this, one has to derive a kinetic equation for the Bogoliubov quasiparticle excitations. This is the goal of the present chapter. In this chapter, we use the second-order Beliaev approximation to discuss the nonequilibrium dynamics of a trapped Bose-condensed gas at finite temperatures. In doing to, we combine the second-order Beliaev self-energies with the lower-order Bogoliubov excitation spectrum, including off-diagonal single-particle propagators but still omitting the anomalous correlation functions. This last condition defines the Bogoliubov–Popov approximation. In this chapter, we consider only the damping which arises from collisions. We will not explicitly calculate corrections that are second-order in g to the quasiparticle energy spectrum or to the condensate chemical potential, both of which are associated with the real parts of the second-order Beliaev self-energies. The present chapter is a natural generalization of work presented in Chapter 6 for the simpler HF excitation spectrum assumed in the ZNG analysis.
Bose system of zero spin particles is considered in the presence of the Bose–Einstein condensate in the vicinity of the phase transition point. The system is investigated in the framework of the Bogolyubov model with the separated condensate. In this model an effective Hamiltonian of the system is introduced by replacing condensate creation and annihilation operators in system Hamiltonian by n01/2 where n0 is occupation number of the condensate state. According to Bogolyubov, the grand canonical thermodynamic potential related to the effective Hamiltonian is considered as nonequilibrium thermodynamic potential. In the present paper this potential is investigated as a function of the small variable n0. With the help of the thermodynamic perturbation theory it is shown that it is expanded in a series over integer powers of n0. This corresponds to the basic idea of the Landau theory of the phase transitions of the second kind. Coefficients at terms of the first and second orders in n0 in the expansion are calculated for Bose gas in the main approximation in small interaction. Calculation of the coefficients at terms of the third and fourth orders needs accounting contributions of the thermodynamic perturbation theory at least of the 4th order and will be done elsewhere. It is established that the results obtained for Bose gas do not fit into the Landau theory of phase transitions of the second kind. Some comments that discuss the situation are given.