In Chapter 3, we introduced a simple but reasonable approximation for the nonequilibrium dynamics of a Bose-condensed gas based on a generalized GP equation coupled to a kinetic equation. In Chapters 4–7, we turn to the question of how to derive such coupled equations for the condensate and noncondensate components in a way that gives a deeper understanding of the ZNG theory. Chapters 4–7 involve an introduction to Green's functions and field theoretic techniques for nonequilibrium problems. These provide the natural language and formalism to deal with the many subtle aspects of a Bose-condensed gas at finite temperatures. These four chapters are fairly technical. This chapter is mainly based on Kadanoff and Baym (1962) and Imamović-Tomasović (2001). Readers who are not interested in these questions can go straight to Chapter 8, which begins the discussion of applications of the ZNG coupled equations given in Chapter 3. Overview of Green's function approach To derive a microscopic theory of the nonequilibrium behaviour of a dilute weakly interacting Bose-condensed gas at finite temperatures, there are several different approaches available in the literature. We will use the wellknown Kadanoff–Baym (KB) nonequilibrium Green's function method. The generalization of this formalism to a Bose-condensed system was first considered by Kane and Kadanoff (1965), whose goal was to derive the Landau two-fluid hydrodynamic equations for a system with a Bose broken symmetry. The general problem consists of how to calculate the nonequilibrium response of a system induced by an external (space- and time-dependent) perturbation. In response to such an external perturbation, many interesting physical phenomena appear, including the excitation of collective modes and various transport processes.
In Chapter 3, we introduced a simple but reasonable approximation for the nonequilibrium dynamics of a Bose-condensed gas based on a generalized GP equation coupled to a kinetic equation. In Chapters 4–7, we turn to the question of how to derive such coupled equations for the condensate and noncondensate components in a way that gives a deeper understanding of the ZNG theory. Chapters 4–7 involve an introduction to Green's functions and field theoretic techniques for nonequilibrium problems. These provide the natural language and formalism to deal with the many subtle aspects of a Bose-condensed gas at finite temperatures. These four chapters are fairly technical. This chapter is mainly based on Kadanoff and Baym (1962) and Imamović-Tomasović (2001). Readers who are not interested in these questions can go straight to Chapter 8, which begins the discussion of applications of the ZNG coupled equations given in Chapter 3. Overview of Green's function approach To derive a microscopic theory of the nonequilibrium behaviour of a dilute weakly interacting Bose-condensed gas at finite temperatures, there are several different approaches available in the literature. We will use the wellknown Kadanoff–Baym (KB) nonequilibrium Green's function method. The generalization of this formalism to a Bose-condensed system was first considered by Kane and Kadanoff (1965), whose goal was to derive the Landau two-fluid hydrodynamic equations for a system with a Bose broken symmetry. The general problem consists of how to calculate the nonequilibrium response of a system induced by an external (space- and time-dependent) perturbation. In response to such an external perturbation, many interesting physical phenomena appear, including the excitation of collective modes and various transport processes.
In Chapter 4, we introduced the Kadanoff–Baym equations of motion for the imaginary-time nonequilibrium Green's functions for a Bose gas, as given by (4.59) and (4.60). In this chapter, we will use the generalization of these equations of motion to find the equivalent equations of motion for the real-time Green's functions. These can be written in a natural way in the form of a kinetic equation. Using a simple Hartree–Fock approximation, we show how the coupled equations for the condensate and thermal cloud given in Chapter 3 emerge naturally from the Kadanoff–Baym (KB) formalism. This chapter is based on Imamović-Tomasović and Griffin (2001) and Imamović-Tomasović (2001), building on the pioneering work of Kane and Kadanoff (1965). In this chapter and Chapter 7 we review the KB formalism. However, we also encourage the reader to read the original account given by Kadanoff and Baym (1962). The goals and accomplishments of their seminal book are beautifully captured by the following quote from p. 138: Our rather elaborate Green's function arguments provide a means of describing transport phenomena in a self-contained way, starting from a dynamical approximation, i.e. an approximation for G 2 ( U ) in terms of G 1 ( U ). These calculations require no extra assumptions. The theory provides at the same time a description of the transport processes that occur and a determination of the quantities which appear in the transport equations. A closely related way of treating the nonequilibrium dynamics of a Bose-condensed gas is based on the two-particle irreducible (2PI) effective action together with the Schwinger–Keldysh closed-time path formalism. Berges (2004) gives a detailed review of this approach, which allows one to derive the nonequilibrium action on the basis of controllable approximations.
In Chapters 11–13, we gave a detailed discussion of the dynamics of a trapped Bose gas at finite temperatures in a region where the collisions described by the C 12 and C 22 terms in the kinetic equation (3.42) do not play the central role. In this “collisionless” region, the dominant interaction effects are associated with the self-consistent fields which both the condensate and noncondensate atoms feel. Thus the dynamics can be understood to a first approximation by neglecting the C 12 and C 22 collision integrals in the kinetic equation and, at the next stage, treating them as a weak perturbation on the collisionless dynamics. In the rest of this book (Chapters 15–19), we turn to the study of the coupled ZNG equations in the opposite limit, where the C 12 and C 22 collision integrals completely determine the dynamics of the thermal cloud. Specifically, the collisions lead to the thermal cloud being in local hydrodynamic equilibrium, and hence this regime is described by the equations of collisional hydrodynamics. Its characteristic feature is that the nonequilibrium behaviour of the thermal cloud atoms can be completely described in terms of a few differential equations involving coarse-grained variables that are dependent on position and time, analogous to the condensate variables n c ( r , t ) and v c ( r , t ). In the present chapter, assuming that the thermal cloud distribution function f ( p, r , t ) is given by the Bose distribution describing partial local equilibrium, (15.16), we show how the ZNG coupled equations lead precisely to Landau's two-fluid equations, reviewed in Chapter 14. This equivalence is not obvious, mainly because Landau's equations are expressed in terms of thermodynamic variables, which are not used in a more microscopic analysis such as that used in the ZNG approach.
In contrast with Chapter 2, in this chapter we include the dynamics of the thermal cloud. As noted in Chapter 1, we treat the noncondensate atoms using the simplest microscopic model approximation that captures the important physics. In particular, we consider only temperatures high enough ( T ≥ 0.4 T BEC ) that the noncondensate atoms can be described by a particle-like Hartree–Fock (HF) spectrum. To extend the analysis to very low temperatures is in principle straightforward (see Chapter 7). However, the details are more complicated since the excitations of the thermal cloud take on a collective aspect (i.e. a Bogoliubov-type quasiparticle spectrum must be used). In trapped Bose gases, the HF single-particle spectrum gives a good approximation down to much lower temperatures than in the case of uniform Bose gases, as first emphasized by Giorgini et al . (1997). In Section 3.1, we derive a generalized form of the Gross–Pitaevskii equation for the Bose order parameter Φ( r , t ) that is valid at finite temperatures. It involves terms that are coupled to the noncondensate component (the thermal cloud) and thus its solution in general requires one to know the equations of motion for the dynamics of the noncondensate atoms. In Section 3.2 we restrict ourselves to finite temperatures high enough that the noncondensate atoms can be described by a quantum kinetic equation for the single-particle distribution function f ( p, r , t ). A detailed microscopic derivation of this kind of kinetic equation is given in Chapters 6 and 7 using the Kadanoff–Baym Green's function formalism. A characteristic feature of a Bose-condensed gas is that the kinetic equation governing f ( p, r , t ) involves a collision integral C 12 [ f , Φ] describing collisions between condensate and noncondensate atoms.
Single-particle Green's functions, density response functions and other correlation functions are calculated in many different ways in the literature on Bose-condensed gases. An in-depth comparison and classification of different approaches was first given in the classic paper by Hohenberg and Martin (1965), with an emphasis on the various exact identities (conservation laws) that are satisfied. A key feature to be included in any theory is that a Bose broken symmetry leads to a hybridization of the single-particle excitations and the collective density fluctuations in such a way that the two excitation spectra become identical. This key feature is demonstrated in Chapter 5 of Griffin (1993). How to relate and assess various approximations for correlation functions in a Bose superfluid has been a topic of continual interest (and some controversy) since the late 1950s. These questions were largely resolved by the early 1960s at a conceptual level but the detailed applications of the theory were limited to dilute Bose-condensed gases at T = 0. Since it was difficult to relate the theory to the properties of superfluid. He at a quantitative level, this formalism based on a Bose broken symmetry was of little interest to experimentalists. The creation of superfluid Bose condensed gases in 1995 changed all this and has given new life to the many body theory of dilute weakly interacting Bose condensed gases. Various approximations for the Beliaev single-particle self-energies Σ αβ were derived in Sections 4.3 and 4.4. The discussion in Chapter 4 was somewhat abstract.
In Chapter 15, we showed that in the limit of short collision times the coupled equations of motion for the condensate and noncondensate atoms lead to Landau's non-dissipative two-fluid hydrodynamics. However the approach used in Chapter 15 was not based on a small expansion parameter, in contrast with the more systematic Chapman–Enskog procedure used to derive hydrodynamic damping in the kinetic theory of classical gases. In the present chapter, we generalize the procedure of Chapter 15 to trapped Bose-condensed gases, in order to derive two-fluid hydrodynamic equations that include dissipation due to transport processes. We solve the kinetic equation by expanding the nonequilibrium single-particle distribution function f ( p, r , t ) around the distribution function f (0) ( p, r , t ) that describes complete local equilibrium between the condensate and the noncondensate components. All hydrodynamic damping effects are included by taking into account deviations from the local equilibrium distribution function f (0) . Our discussion for a trapped Bose gas is a natural extension of the pioneering work of Kirkpatrick and Dorfman (1983, 1985a) for a uniform Bose-condensed gas. This chapter is mainly based on their work as well as on Nikuni and Griffin (2001a,b). We will prove that, with appropriate definitions of various thermodynamic variables, our two-fluid hydrodynamic equations including damping have precisely the structure of those first derived by Landau and Khalatnikov for superfluid He. In particular, the damping associated with the collisional exchange of atoms between the condensate and noncondensate components, which is discussed at length in Chapter 15, is now expressed in terms of frequency-dependent second viscosity coefficients. This special type of damping is a characteristic signature of a dilute Bose superfluid and exists in addition to the hydrodynamic damping associated with the shear viscosity and thermal conductivity of the normal fluid.
In this chapter we describe the numerical methods that can be used to solve the ZNG equations in the context of a dynamical simulation. These equations consist of a generalized GP equation (3.21) for the condensate and a Boltzmann equation (3.42) for the thermal component. The fact that the two equations are coupled makes their numerical solution more complex than when either is considered on its own. Indeed, the distinct quantum and classical aspects of the problem require specifically tailored numerical methods. Although most of these methods are well established and described elsewhere (Taha and Ablowitz, 1984; Sanz-Serna and Calvo, 1994), we provide in this chapter a detailed pedagogical discussion that will serve as a guide to those interested in carrying out such calculations for trapped Bose gases. This chapter is based on the papers of Jackson and Zaremba (2002a,b). There are two main parts to the numerical problem. The first is developing a method for solving the time-dependent GP equation for an arbitrary three-dimensional geometry. This we take up in Section 11.1. Second, a method is needed for solving the Boltzmann equation that accounts for the dynamics of the thermal component. Here one must deal both with the Hamiltonian dynamics of the thermal atoms, as they move in the self-consistent mean field of the condensate and thermal cloud, and with the collisions that take place between the thermal atoms themselves (the C 22 collisions) and between the thermal atoms and the condensate (the C 12 collisions). The methods used to account for these two distinct collisional processes are taken up in Section 11.3. As we shall see, collisions play an important role and cannot be neglected even when the dynamical behaviour is dominated by mean-field interactions.