In this chapter, we review the famous Landau theory of superfluidity at finite temperatures. This theory is based on coupled hydrodynamic equations for the superfluid and normal fluid components. Landau's two-fluid description is only valid when collisions among the thermal excitations making up the normal fluid are strong enough to produce local hydrodynamic equilibrium. These two-fluid equations were originally developed for liquid. He but are thought to be generic in form, describing the collision-dominated hydrodynamic region of all Bose superfluids. In this chapter, we will consider the solutions of the two-fluid equations mainly for a uniform superfluid. We discuss the existence of second sound (involving the out-of-phase motion of the superfluid and normal fluid components) as a characteristic feature of a Bose superfluid at finite temperatures. This chapter gives background material needed for Chapters 15–19. In Chapter 15, we will show that, in the appropriate limit, the Landau two-fluid equations can be derived from the ZNG coupled equations given in Chapter 3 for a trapped dilute Bose-condensed gas. In Chapters 17–19, we extend this discussion and derive the Landau–Khalatnikov two-fluid equations, which include hydrodynamic damping associated with various transport coefficients. Useful reviews of the two-fluid equations in the context of dilute spatially uniform Bose-condensed gases are given by Pethick and Smith (2008, Chapter 10) and Pitaevskii and Stringari (2003, Chapter 6). History of two-fluid equations The original discovery of superfluidity in liquid He was dramatically announced with the publication of the famous back-to-back papers of Kapitza (1938) in Moscow and Allen and Misener (1938) in Cambridge. These and subsequent experiments in the next few years showed that superfluid He could exhibit very bizarre hydrodynamic behaviour compared to classical liquids.
In this chapter, we review the famous Landau theory of superfluidity at finite temperatures. This theory is based on coupled hydrodynamic equations for the superfluid and normal fluid components. Landau's two-fluid description is only valid when collisions among the thermal excitations making up the normal fluid are strong enough to produce local hydrodynamic equilibrium. These two-fluid equations were originally developed for liquid. He but are thought to be generic in form, describing the collision-dominated hydrodynamic region of all Bose superfluids. In this chapter, we will consider the solutions of the two-fluid equations mainly for a uniform superfluid. We discuss the existence of second sound (involving the out-of-phase motion of the superfluid and normal fluid components) as a characteristic feature of a Bose superfluid at finite temperatures. This chapter gives background material needed for Chapters 15–19. In Chapter 15, we will show that, in the appropriate limit, the Landau two-fluid equations can be derived from the ZNG coupled equations given in Chapter 3 for a trapped dilute Bose-condensed gas. In Chapters 17–19, we extend this discussion and derive the Landau–Khalatnikov two-fluid equations, which include hydrodynamic damping associated with various transport coefficients. Useful reviews of the two-fluid equations in the context of dilute spatially uniform Bose-condensed gases are given by Pethick and Smith (2008, Chapter 10) and Pitaevskii and Stringari (2003, Chapter 6). History of two-fluid equations The original discovery of superfluidity in liquid He was dramatically announced with the publication of the famous back-to-back papers of Kapitza (1938) in Moscow and Allen and Misener (1938) in Cambridge. These and subsequent experiments in the next few years showed that superfluid He could exhibit very bizarre hydrodynamic behaviour compared to classical liquids.
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 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.
In the collisional region at finite temperatures, the collective modes of superfluids are described by the Landau two-fluid hydrodynamic equations reviewed in Chapter 14. In the case of trapped Bose gases, these are coupled differential equations with position-dependent coefficients associated with the local thermodynamic functions. Building on the approach initiated by Zaremba et al. (1999) for trapped atomic Bose gases, in this chapter we develop an alternative variational formulation of two-fluid hydrodynamics. This is based on the work of Zilsel (1950), originally developed to deal with superfluid He. Assuming a simple variational ansatz for the superfluid and normal fluid velocity fields, this approach reduces the problem of finding the hydrodynamic collective mode frequencies to solving coupled algebraic equations for a few variational parameters. These equations contain constants that involve spatial integrals over various equilibrium thermodynamic derivatives. Such a variational approach is both simpler and more physical than a direct attempt to solve the Landau two-fluid equations numerically. This chapter is mainly based on Taylor and Griffin (2005), Taylor (2008) and Zilsel (1950). In it, we discuss the normal modes of the non-dissipative Landau two-fluid equations for a trapped superfluid. In Section 16.3, we illustrate this formalism by deriving expressions for the frequencies of the dipole and breathing modes of a trapped Bose superfluid. In Chapters 17 and 18, we discuss an extended version of the two-fluid equations that includes hydrodynamic damping. The hydrodynamic damping of the collective modes is calculated in Chapter 19 using a generalized version of the variational approach developed in this chapter. Zilsel's variational formulation Since two-fluid hydrodynamics only describes a system in local equilibrium, all thermodynamic quantities are functions of position and time. Even in static equilibrium, in the presence of a trapping potential, most thermodynamic quantities will be position dependent.
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.
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.
In the collisional region at finite temperatures, the collective modes of superfluids are described by the Landau two-fluid hydrodynamic equations. This region can now be probed over the entire BCS--Bose-Einstein-condensate crossover in trapped Fermi superfluids with a Feshbach resonance, including the unitarity region. Building on the approach initiated by Zaremba, Nikuni, and Griffin in 1999 for trapped atomic Bose gases, we present a variational formulation of two-fluid hydrodynamic collective modes based on the work of Zilsel in 1950 developed for superfluid helium. Assuming a simple variational Ansatz for the superfluid and normal fluid velocities, the frequencies of the hydrodynamic modes are given by solutions of coupled algebraic equations, with constants only involving spatial integrals over various equilibrium thermodynamic derivatives. This variational approach is both simpler and more physical than a direct attempt to solve the Landau two-fluid differential equations. Our two-fluid results are shown to reduce to those of Pitaevskii and Stringari for a pure superfluid at T=0.