In Chapter 17, we derived two-fluid hydrodynamic equations that include damping related to transport coefficients. Our entire analysis was based on the coupled ZNG equations for the condensate and in the thermal cloud. These involved a generalized GP equation for the condensate and a kinetic equation for the thermal atoms. A crucial role is played by the C 12 collision term in the kinetic equation, which describes the interactions between atoms in the condensate and in the thermal cloud. Our analysis of the deviation from the diffusive local equilibrium solution of the kinetic equation was based on the Chapman–Enskog approach, extensively developed for classical gases and first applied to Bose-condensed gases by Kirkpatrick and Dorfman (1983, 1985a). This approach required a careful treatment of the novel feature relating to the C 12 collisions both in the kinetic equation describing the thermal atoms and also in the source term Γ 12 in the generalized GP equation for the condensate. Using the Chapman– Enskog approach to solve the kinetic equation for a trapped Bose gas, we obtained explicit expressions for the function ψ( p, r , t ) that describes the deviation from diffusive local equilibrium, as defined by (17.25) and (17.39). This deviation can be related to various transport coefficients, as discussed in Chapter 17. These transport coefficients are determined by the solutions of the three integral equations (17.40)–(17.42) for the three contributions to the deviation function ψ( p, r , t ) in (17.39). In Section 18.1, we will solve these integral equations and obtain explicit expressions for the thermal conductivity k , the shear viscosity η and the four second viscosity coefficients ζ i .
In Chapter 17, we derived two-fluid hydrodynamic equations that include damping related to transport coefficients. Our entire analysis was based on the coupled ZNG equations for the condensate and in the thermal cloud. These involved a generalized GP equation for the condensate and a kinetic equation for the thermal atoms. A crucial role is played by the C 12 collision term in the kinetic equation, which describes the interactions between atoms in the condensate and in the thermal cloud. Our analysis of the deviation from the diffusive local equilibrium solution of the kinetic equation was based on the Chapman–Enskog approach, extensively developed for classical gases and first applied to Bose-condensed gases by Kirkpatrick and Dorfman (1983, 1985a). This approach required a careful treatment of the novel feature relating to the C 12 collisions both in the kinetic equation describing the thermal atoms and also in the source term Γ 12 in the generalized GP equation for the condensate. Using the Chapman– Enskog approach to solve the kinetic equation for a trapped Bose gas, we obtained explicit expressions for the function ψ( p, r , t ) that describes the deviation from diffusive local equilibrium, as defined by (17.25) and (17.39). This deviation can be related to various transport coefficients, as discussed in Chapter 17. These transport coefficients are determined by the solutions of the three integral equations (17.40)–(17.42) for the three contributions to the deviation function ψ( p, r , t ) in (17.39). In Section 18.1, we will solve these integral equations and obtain explicit expressions for the thermal conductivity k , the shear viscosity η and the four second viscosity coefficients ζ i .
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 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 Chapter 17, we derived the Landau–Khalatnikov two-fluid hydrodynamics which describes the collision-dominated region of a trapped Bose condensate interacting with a thermal cloud. In this chapter, we use these equations to discuss the damping of the hydrodynamic collective modes in a trapped Bose gas at finite temperatures. We derive variational expressions based on these equations for both the frequency and the damping of collective modes. This extends the analysis in Chapter 16 in which a variational approach was developed to calculate the hydrodynamic two-fluid oscillation frequencies in the non-dissipative limit. A novel feature of our treatment is the introduction of frequency-dependent transport coefficients, which produce a natural cutoff eliminating the collisionless region in the low-density tail of the thermal cloud. Our expression for the damping in trapped superfluid Bose gases is a natural generalization of the approach used by Landau and Lifshitz (1959) for uniform classical fluids. This chapter is mainly based on Nikuni and Griffin (2004). In Chapters 15 and 17, we derived a closed set of two-fluid hydrodynamic equations for a trapped Bose-condensed gas starting from the simplified microscopic model describing the coupled dynamics of the condensate and noncondensate atoms given in Chapter 3. These hydrodynamic two-fluid equations include dissipative terms associated with the shear viscosity, the thermal conductivity and the four second viscosity coefficients. Explicit formulas for these transport coefficients were derived in Section 18.1. Our goal in this chapter is to find a general expression for the damping of the two-fluid modes in terms of these transport coefficients. We emphasize that the damping of hydrodynamic two-fluid oscillations is completely different in nature from the Landau and Beliaev damping of oscillations in the collisionless region which is treated in Chapters 12 and 13.
Starting from the quantum kinetic equation for the noncondensate atoms and the generalized Gross--Pitaevskii equation for the condensate, we derive the two-fluid hydrodynamic equations of a trapped Bose gas at finite temperatures. We follow the standard Chapman--Enskog procedure, starting from a solution of the kinetic equation corresponding to the complete local equilibrium between the condensate and the noncondensate components. Our hydrodynamic equations are shown to reduce to a form identical to the well-known Landau-Khalatnikov two-fluid equations, with hydrodynamic damping due to the deviation from local equilibrium. The deviation from local equilibrium within the thermal cloud gives rise to dissipation associated with shear viscosity and thermal conduction. In addition, we show that effects due to the deviation from the diffusive local equilibrium between the condensate and the noncondensate (recently considered by Zaremba, Nikuni, and Griffin) can be described by four frequency-dependent second viscosity transport coefficients. We also derive explicit formulas for all the transport coefficients. These results are used to introduce two new characteristic relaxation times associated with hydrodynamic damping. These relaxation times give the rate at which local equilibrium is reached and hence determine whether one is in the two-fluid hydrodynamic region.
Griffin, Wu and Stringari have derived the hydrodynamic equations of a trapped dilute Bose gas above the Bose-Einstein transition temperature. We give the extension which includes hydrodynamic damping, following the classic work of Uehling and Uhlenbeck based on the Chapman-Enskog procedure. Our final result is a closed equation for the velocity fluctuations δv which includes the hydrodynamic damping due to the shear viscosity η and the thermal conductivity κ. Following Kavoulakis, Pethick and Smith, we introduce a spatial cutoff in our linearized equations when the density is so low that the hydrodynamic description breaks down. Explicit expressions are given for η and κ, which are position-dependent through dependence on the local fugacity when one includes the effect of quantum degeneracy of the trapped gas. We also discuss a trapped Bose-condensed gas, generalizing the work of Zaremba, Griffin and Nikuni to include hydrodynamic damping due to the (non-condensate) normal fluid.
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.