arxiv:cond-mat/ v1 5 Nov 1997
|
|
- Ruby Craig
- 5 years ago
- Views:
Transcription
1 Hydrodynamic damping in trapped Bose gases T. Nikuni Department of Physics, Tokyo Institute of Technology, Oh-okayama, Tokyo, 152 Japan and Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 arxiv:cond-mat/ v1 5 Nov 1997 A. Griffin Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 (February 29, 2008) Abstract 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 κ. 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. A novel type of contribution arises which depends on the gradient of the shear viscosity. It leads to certain normal mode solutions (with δv = const) becoming unstable (i.e. exponential growth with time), which signals the breakdown of the Chapman-Enskog assumption that deviations from local thermal equilibrium are small. 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. PACS numbers: Fi, Jp, Db Typeset using REVTEX 1
2 I. INTRODUCTION Recent observation of Bose-Einstein condensation in trapped atomic gases has stimulated interest in the collective oscillations of non-uniform Bose gases [1,2]. Ignoring damping effects, the hydrodynamic equations of a trapped Bose gas have been recently derived above [3] and below [4] the Bose-Einstein transition temperature (T BEC ). These hydrodynamic equations describe low frequency phenomena (ω 1/τ c, where τ c is the mean time between collisions of atoms), when collisions are sufficiently strong to ensure local thermodynamic equilibrium. These results can be used to derive closed equations for the velocity fluctuations in a trapped gas, whose solution gives the normal modes of oscillation. In particular, above T BEC, Ref. [3] gives explicit results for the frequencies of surface, monople, and coupled monopole-quadrupole modes. Available experimental data on the frequencies and damping of collective oscillations in trapped gases [1,2] is mainly for the collisionless regime, rather than the hydrodynamic regime we consider in this paper. However sufficiently high densities which allow one to probe the hydrodynamic region at finite temperatures are expected in the near future (see also Refs. [5,6]). Thus it is worthwhile to extend the discussion of hydrodynamic modes in Refs. [3,4] to include damping. In the present paper, we do this by following the standard Chapman-Enskog procedure, as first generalized for quantum gases by Uehling and Uhlenbeck [7,8]. In Section II, we derive the hydrodynamic equations of a trapped Bose gas above T BEC which include damping due to shear viscosity η and thermal conductivity κ. We also obtain integral equations which give expressions for the transport coefficients η and κ, and in Section III we solve for these. Our approximation corresponds to a lowest-order polynomial approximation of the Chapman- Enskog results such as obtained by Uehling [8]. The transport coefficients have a position dependence when we include the effect of quantum degeneracy of the trapped Bose gas. In Section IV, we evaluate the hydrodynamic damping of the normal modes in a trapped Bose gas discussed in Ref. [3]. A novel contribution arises from the position dependence of the shear viscosity coefficient, which leads to the surface and coupled monopole-surface modes becoming unstable. This means that, for these modes, we cannot expand around the local thermal equilibrium solution assumed in the Chapman-Enskog approach. This in turn 2
3 means that the usual hydrodynamic equations must be used with care in discussing trapped Bose gases. In Section V, we extend these results to the Bose-condensed region just below T BEC. Zaremba, Griffin and Nikuni [4] have recently given a microscopic derivation of the two-fluid hydrodynamic equations of motion for a reasonable model of a trapped weakly-interacting Bose-condensed gas. We generalize the results of Ref. [4] to include hydrodynamic damping. Finally in Section VI, we give some concluding remarks. II. CHAPMAN-ENSKOG METHOD In Section II-IV, we limit ourselves to a non-condensed Bose gas above T BEC, as in Ref. [3]. The atoms are described by the semi-classical kinetic equation for the distribution function f(r,p, t) of a Bose gas [9] [ t + p ] m r U(r, t) p f(r,p, t) = f t. (1) coll Here U(r, t) = U 0 (r) + 2gn(r, t) includes an external potential U 0 (r) as well as the Hartree- Fock (HF) self-consistent mean field 2gn(r, t), where n(r, t) is the local density. In usual discussions [7,8], the HF field is omitted but it will play a crucial role when we discuss a Bose-condensed gas below T BEC (see Section V). The quantum collision integral in the right hand side of (1) is given by [9] f ( ) = 4πg 2 dp 2 dp 3 dp 4 p 2 t coll (2π) 3 (2π) 3 (2π) 3δ(p + p 2 p 3 p 4 )δ 2m + p2 2 2m p2 3 2m p2 4 2m [(1 + f)(1 + f 2 )f 3 f 4 ff 2 (1 + f 3 )(1 + f 4 )], (2) where f f(r,p, t), f i f(r,p i, t). The interaction strength g = 4π h 2 a/m is determined by the s-wave scattering length a. The conservation laws are obtained by multiplying (1) by 1, p and p 2 and integrating over p. In all three cases, the integrals of the collision term in (2) vanish and one finds the general hydrodynamic equations: n + (nv) = 0, t (3a) ( ) mn t + v v µ + P µν + n U = 0, x ν x µ ε t + (εv) + Q + D µνp µν = 0, (3b) (3c) 3
4 where n(r, t) = dp (2π) 3f(r,p, t), (4) n(r, t)v(r, t) = dp p f(r,p, t), (5) (2π) 3 m P µν (r, t) = m ( ) ( ) dp pµ (2π) 3 m v pν µ m v ν f(r,p, t), (6) ε(r, t) = dp (2π) 3 1 2m (p mv)2 f(r,p, t), (7) Q(r, t) = ( ) dp 1 p (p (2π) 3 mv)2 2m m v f(r,p, t), (8) D µν = 1 2 ( vµ + v ) ν. (9) x ν x µ The first approximation for the distribution function is the local equilibrium form f (0) (r,p, t) = {exp[β(p mv) 2 /2m + U µ] 1} 1, (10) where the thermodynamic variables β,v and µ all depend on r and t, and U(r, t) has been defined after (1). This expression ensures that the collision integral (2) vanishes. This gives the lowest-order approximation in l/l, where l is the mean free path and L is a characteristic wavelength. If one uses (10) to calculate the quantities from (4) to (8), one finds that Q = 0 and n(r, t) = 1 Λ 3g 3/2(z(r, t)), (11) P µν (r, t) = δ µν P(r, t), P(r, t) = k BT Λ g 5/2(z(r, t)) = 2 ε(r, t). (12) 3 3 Here z(r, t) e β(r,t)[µ(r,t) U(r,t)] is the local fugacity, Λ(r, t) = [2π h 2 /mk B T(r, t)] 1/2 is the local thermal de Broglie wavelength and g n (z) = l=1 z l /l n are the well-known Bose-Einstein functions. The equilibrium value of the fugacity is given by z 0 (r) = e β 0(µ U(r)), with U(r) = U 0 (r) + 2gn 0 (r). Putting Q = 0 and using (12) in Eqs. (3), one obtains 4
5 n + (nv) = 0, t (13a) [ ] v mn + (v )v = P n U, (13b) t ε t + 5 (εv) = v P. 3 (13c) The linearization of the equations in (13) around equilibrium leads to the equations (4-6) in Ref. [3] if we ignore the HF field in U(r, t). Solutions of these equations describe undamped oscillations, some example of which are discussed in Ref. [3] (see also Section IV). In order to obtain damping of the oscillations, we have to consider the deviation of the distribution function from the local equilibrium form (10). One assumes a solution of the quantum Boltzmann equation (1) of the form [7,10]: f(r,p, t) = f (0) (r,p, t) + f (0) (r,p, t)[1 + f (0) (r,p, t)]ψ(r,p, t), (14) where ψ expresses a small deviation from local equilibrium. Working to first order in ψ, the collision integral in (2) can be reduced to ( ) 4πg 2 dp 2 dp 3 p 2 dp (2π) 3 (2π) 3 4 δ(p + p 2 p 3 p 4 )δ 2m + p2 2 2m p2 3 2m p2 4 2m f (0) f (0) 2 (1 + f (0) 3 )(1 + f (0) 4 )(ψ 3 + ψ 4 ψ 2 ψ) ˆL[ψ], (15) where ψ i ψ(r,p i, t). In the left hand side of (1), we approximate f by f (0). The various derivatives of v(r, t), µ(r, t), T(r, t) and U(r, t) with respect to r and t can be rewritten using the lowest-order hydrodynamic equations given in (13). The resulting linearized equation for ψ is (for details, see Appendix) { [ u T mu 2 T 2k B T 5g ] 5/2(z) + m ( 2g 3/2 (z) k B T D µν u µ u ν 1 ) } 3 δ µνu 2 f (0) (1 + f (0) ) = ˆL[ψ], (16) where the thermal velocity u is defined by mu p mv and the strain tensor D µν is defined in (9). The linearized collision operator ˆL is defined by (15). This equation can be shown to have a unique solution for ψ if we impose the constraints dpf (0) (1 + f (0) )ψ = dpp µ f (0) (1 + f (0) )ψ = dpp 2 f (0) (1 + f (0) )ψ = 0, (17) which mean physically that n,v and ε [see (4),(5) and (7)] are determined only by the first term f (0) in (14). Since the left hand of (16) and also f (0) only depend on the relative 5
6 thermal velocity u, the function ψ will also depend only on u (i.e., not separately on p and v). It is convenient to introduce dimensionless velocity variables ( ) m 1/2 u ξ. (18) 2k B T With these dimensionless velocity variables, (16) becomes π 3 ( ) 1/2 [ 2kB T T ξ ξ 2 5g ] 5/2(z) + 2D 8a 2 mkbt 2 2 µν (ξ µ ξ ν 13 ) m T 2g 3/2 (z) δ µνξ 2 f(0) (ξ)[1 + f (0) (ξ)] with f (0) (ξ) = (z 1 e ξ2 1) 1 and ˆL [ψ] dξ 2 dξ 3 dξ 4 δ(ξ + ξ 2 ξ 3 ξ 4 )δ ( ξ 2 + ξ2 2 ξ2 3 ) ξ2 4 = ˆL [ψ], (19) f (0) f (0) 2 (1 + f (0) 3 )(1 + f (0) 4 )(ψ 3 + ψ 4 ψ 2 ψ). (20) For a more detailed discussion of the mathematical structure of (19) and (20), we refer to the treatment of the analogous equations for classical gases (see for example, Ref. [11]). The most general solution of the integral equation (19) is of the form [7] π 3 ( ) 1/2 ( 2kB T T ξ ψ = A(ξ) + 2D 8a 2 mkbt 2 2 µν ξ µ ξ ν 1 m T 3 δ µνξ )B(ξ) 2. (21) The functions A(ξ) and B(ξ) obey the following integral equations: [ ˆL [ξa(ξ)] = ξ ξ 2 5g ] 5/2(z) f (0) (1 + f (0) ), 2g 3/2 (z) [( ˆL ξ µ ξ ν 1 ] 3 δ µνξ )B(ξ) 2 = (ξ µ ξ ν 1 ) 3 δ µνξ 2 f (0) (1 + f (0) ). (22a) (22b) For (21) to satisfy the constraints given in (17), we must also have dξξ 2 A(ξ)f (0) (1 + f (0) ) = 0. (23) Using a solution of the form (21) in conjunction with (14), one can calculate the heat current density Q in (8) and the pressure tensor P µν in (6). One finds these have the form Q = κ T, (24a) [ P µν = δ µν P 2η D µν 1 ] 3 (TrD)δ µν, (24b) 6
7 where the thermal conductivity κ and the generalized shear viscosity coefficient η are given in terms of the functions A(ξ) and B(ξ): κ = k ( ) 1/2 B 2kB T dξξ 4 A(ξ)f (0) (1 + f (0) ), 48a 2 m (25a) η = m ( ) 1/2 2kB T dξξ 4 B(ξ)f (0) (1 + f (0) ). 120a 2 m (25b) Introducing (24) into the general hydrodynamic equations (3), one obtains the following hydrodynamic equations with the effect of shear viscosity and heat conduction included: n t mn + (nv) = 0, (26a) ( ) t + v v µ + P + n U = { [ 2η D µν 1 ]} x µ x µ x ν 3 (TrD)δ µν, (26b) [ ε t + (εv) + ( v)p = (κ T) + 2η D µν 1 2 µν] 3 (TrD)δ. (26c) The form of these equations can be shown to agree with those originally obtained by Uehling and Uhlenbeck [7]. We note that η and κ are position-dependent, but only through their dependence on the equilibrium value of the fugacity z = z 0 (r). One slight generalization we have made over the derivation in Ref. [7] is that we have included the Hartree-Fock mean field. III. THE TRANSPORT COEFFICIENTS A. The thermal conductivity To find the thermal conductivity as given by (25a), we can introduce a simple ansatz for the form of the function A(ξ) [8,11]: [ A(ξ) = A ξ 2 5g ] 5/2(z). (27) 2g 3/2 (z) The constant A is determined by multiplying (22a) by ξ[ξ 2 5g 5/2 (z)/2g 3/2 (z)] and integrating over ξ: A = [ dξξ 2 ξ 2 5g ] 2 5/2(z) f (0) (1 + f (0) ) 2g 3/2 (z) { dξ ( ξ 2 5g 5/2(z) 2g 3/2 (z) )ξ ˆL [( 7 ξ 2 5g 5/2(z) 2g 3/2 (z) ) ξ ]} 1
8 [ = 15π3/2 7 4I A 2 g 7/2(z) 5g2 5/2 (z) ], (28) 2g 3/2 (z) where the integral I A is defined by [ I A dξ ξ 2 5g ] [( 5/2(z) ξ 2g 3/2 (z) ˆL ξ 2 5g ) 5/2(z) 2g 3/2 (z) ] ξ = dξξ 2 ξ ˆL [ξ 2 ξ]. (29) In order to evaluate the integral in (29), it is convenient to introduce the change of variables ξ = 1 2 (ξ 0 + ξ ), ξ 2 = 1 2 (ξ 0 ξ ), ξ 3 = 1 2 (ξ 0 + ξ ), ξ 4 = 1 2 (ξ 0 ξ ). (30) Then we introduce transformations from ξ xξ yξ z to ξ θ φ and ξ xξ yξ z to ξ θ φ, where θ, θ and φ, φ are the polar and azimuthal angles with respect to the vector ξ 0. One obtains the following expression for I A, I A = 4 2π 3 I A, 1 1 I A(z) = dξ 0 dξ dy dy ξ0ξ 4 7 F(ξ 0, ξ, y, y ; z)[y 2 + y 2 2y 2 y 2 ], (31) where y = cos θ, y = cosθ. Here the function F is defined by with F f (0) f (0) 2 (1 + f (0) 3 )(1 + f (0) 4 ) = 0 z 2 e (ξ2 0 +ξ 2 ) (1 ze ξ2 )(1 ze ξ2 2)(1 ze ξ2 3)(1 ze ξ2 4), (32) ξ 2 = 1 2 (ξ ξ 0ξ y + ξ 2 ), ξ 2 2 = 1 2 (ξ2 0 2ξ 0ξ y + ξ 2 ), ξ 2 3 = 1 2 (ξ ξ 0ξ y + ξ 2 ), ξ 2 4 = 1 2 (ξ2 0 2ξ 0ξ y + ξ 2 ). (33) Inserting the expression in (27) into (25a) and carrying out the integration, we obtain the following expression for the thermal conductivity κ: κ = 75k ( ) 1/2 B mkb T π 1/2 [ 7 64a 2 m π 16I A(z) 2 g 7/2(z) 5g2 5/2 (z) ] 2, (34) 2g 3/2 (z) where the function I A(z) is defined in (31). 8
9 B. The shear viscosity In evaluating the shear viscosity in (25b), the simplest consistent approximation [8,11] is to use B(ξ) B. The constant B can be determined by multiplying (22b) by (ξ µ ξ ν δ µν ξ 2 /3) and integrating over ξ, { ( B = dξ ξ µ ξ ν 1 ) 2 { ( 3 δ µνξ 2 f (0) (1 + f )} (0) dξ ξ µ ξ ν 1 ) 3 δ µνξ 2 ˆL [ξ µ ξ ν 1 ]} 1 3 δ µνξ 2 = 5π3/2 g 5/2 (z) 2I B. The function I B is defined by I B = ( dξ ξ µ ξ ν 1 ) 3 δ µνξ 2 ˆL [ξ µ ξ ν 1 ] 3 δ µνξ 2 = 2 2π 3 I B, dξ(ξ µ ξ ν )ˆL [ξ µ ξ ν ] (36) (35) where 1 1 I B (z) = dξ 0 dξ dy dy F(ξ 0, ξ, y, y ; z)ξ0 2 ξ 7 (1 + y 2 + y 2 3y 2 y 2 ). (37) This involves the same function F as defined in (32). Using the ansatz B(ξ) B in conjunction with (35) in (25b), we can carry out the integral. Our final the expression for the viscosity η is η = a 2 ( ) 1/2 mkb T π 1/2 π 8I B(z) g2 5/2 (z). (38) C. High-temperature limit The formulas in (34) and (38) give κ and η as a function of the fugacity z. We note that the integrals I A and I B also depend on z. These four-dimensional integrals in (31) and (37) can be evaluated numerically. In Figs. 1 and 2, we plot both κ and η as a function of the fugacity z. These graphs are valid for any trapping potential U 0 (r) since the latter only enters into the equilibrium fugacity z = exp[β 0 (µ U 0 (r) 2gn 0 (r))]. In the high-temperature limit, where the fugacity z = e β(µ U(r)) is small, the local distribution function f (0) can be expanded in terms of z. To third order in z, the function F in (32) reduces to 9
10 F = z 2 e (ξ2 0 +ξ 2) + 2z 3 e 3 2 (ξ2 0 +ξ 2) [cosh(ξ 0 ξ y ) + cosh(ξ 0 ξ y )] + O(z 4 ). (39) The integrals in (31) and (37) can be evaluated analytically, I B(z) = 2I A(z) = π 1/2 z z 9 3. (40) 16 2 We thus obtain the following explicit expressions for the transport coefficients to first order in the fugacity z: κ = 1 ( ) ( ) 1/2 75 kb kb T 1 + z a 2 πm ( ) ( ) 1/2 75 kb kb T (1 0.07n 8 64 a 2 0 Λ 3 ), (41a) πm η = 1 ( ) ( ) 1/2 5 m kb T 1 + z a 2 πm ( ) ( ) 1/2 5 m kb T ( n0 Λ 3). (41b) 8 16 a 2 πm The local equilibrium density n 0 (r) in the high-temperature limit is the classical result n 0 = Λ 3 e β(µ U(r)). The terms first order in z in (41a) and (41b) give the first order corrections to the classical results due to Bose statistics. If we ignore the HF mean field 2gn 0 (r), these transport coefficients reduce to the expressions first obtained for a uniform gas by Uehling [8]. The Bose quantum corrections to the classical results in both η and κ depend on the local fugacity z and, through this, on position in a trapped gas. The density-independent terms in Eqs. (41a) and (41b) are 8 times smaller than the wellknown Chapman-Enskog expressions for classical hard spheres. This is due to the difference (see also Ref. [8]) in the quantum binary scattering cross-section for Bosons when correctly calculated using symmetrized wavefunctions (σ = 8πa 2, instead of πa 2 ). IV. HYDRODYNAMIC DAMPING OF NORMAL MODES The linearized version of hydrodynamic equations in (26) are δn t + (n 0δv) = 0, mn 0 δp t t δv = δp δn U 2gn 0 δn + η 3 ( δv) + η 2 δv + R, (42a) (42b) = 5 3 (P 0δv) δv P (κ T), (42c) 3 10
11 where the vector R is defined by R µ 2 η [ D µν 1 ] x ν 3 (TrD)δ µν. (43) We shall see that the contribution R leads to quite novel effects in trapped quantum gases, which do not seem to have been studied before. In the absence of an external field (uniform system), the thermal conductivity and the viscosity are position-independent and then the hydrodynamic equations (42a, b, c) reduce to the standard textbook results (where R 0). Taking the time derivative of (42b) and using (42a) and (42c) gives a closed equation for the velocity fluctuations m 2 δv t 2 = 5P 0 3n 0 ( δv) (δv U 0 ) 2 3 ( δv) U 0 +2g [n 0 ( δv)] 4 3 g( δv) n [ ( η n 0 3 δv ) + η 2 δv t t + R ] t 2 3n 0 2 (κ δt), (44) where we have used the relation P 0 = n 0 U valid in equilibrium. Since we assume the dissipative terms to be small, we use the lowest-order hydrodynamic equation for δt in (44) (see (A7) in Appendix), δt t = 2 3 T 0 δv. (45) Equation (44) generalizes that given by (13) in Ref. [3] to include damping due to viscosity and thermal conductivity as well as the HF mean field 2gn. In the following discussion of normal modes above T BEC, we follow Ref. [3] and ignore the HF mean field. We look for normal mode oscillations with the convention δv(r, t) = δv ω (r)e iωt. We first consider an isotropic trap U 0 (r) = 1 2 mω2 0 r2. For a monopole oscillation described by δv ω (r) r, the undamped frequency is found to be ω = 2ω 0 [3]. Using this solution in (44), we note that δv = 0 and 2 δv = 0. Moreover, we find D µν δ µν and hence the contribution R in (43) also vanishes. Thus there is no shear viscous contribution to the damping of this mode. Since δv =const, δt is independent of position (see (45)), and thus the last term in (44) involving the thermal conductivity coefficient κ does not contribute. We conclude the monopole mode is isothermal and has no damping for an isotropic harmonic trap. As noted in Ref. [3], this surprising result was first derived for a classical gas by Boltzmann in
12 The undamped solution δv = δv ω (r)e iωt for divergence-free surface modes ( δv = 0) is given by [3] δv ω (r) [r l Y lm ], (46) with the dispersion relation ω = lω 0, l = 1, 2,. The last term in (44) vanishes because the zero-divergence solutions involve isothermal oscillations [3]. The only contribution to damping comes from the R/ t term in (44). Using the solution (46) in (43), we find that R = 2n 0 C(l 1)mω0 2 δv, (47) where C = (z/n 0 k B T)( η/ z). We thus conclude that the hydrodynamic damping of the surface modes (46) is entirely related to the position dependence of the viscosity coefficient η though the dependence on the local fugacity z. Using (46) and (47) in (44) leads to mω 2 δv = lmω 2 0δv + 2iωC(l 1)mω 2 0δv. (48) For T T BEC, we can use the approximate expression (41) for η. In this case, which only includes the lowest order correction from Bose statistics, C is independent of position and is given by h 3 C = ma 2 (k B T) 2. (49) In this case, the solution of (48) is trivial, giving ω = lω 0 iγ, with ( ) 2 Γ = (l 1)Cω0 2 h hω0 = 0.116(l 1). (50) ma 2 k B T For l = 1, Γ = 0; this means that the center-of-mass mode has no hydrodynamic damping. For l > 1, Γ is always negative and the solution (46) is unstable (with an amplitude exponentially increasing with time e Γ t ). More generally, we note that from Fig. 2, η/ z is always negative and hence so is C(z). In an anisotropic trap potential described by U 0 (r) = 1 2 m(ω2 s2 + ω 2 z z2 ), s 2 x 2 + y 2, (51) the monopole mode is coupled to the quadrupole mode, as discussed in Ref. [3]. The undamped solution has the form δv ω = χ, with χ = αs 2 + βz 2. Putting this solution into (44) leads to two coupled equations for α, β: 12
13 ( 10 mω 2 α = mω 2 3 α + 2 ) 3 β + iω 2 3 Cmω2 (α β), ( 4 mω 2 β = mωz 2 3 α + 8 ) 3 β iω 4 3 Cmω2 z (α β), (52) where the function C defined below (47) is position independent for T T BEC (see (49)). The solution of these equations gives ω = ω ± iγ ±, with ω 2 ± = 1 3 ) (5ω 2 + 4ω2z ± 25ω ω4 z 32ω2 ω2 z, (53) and Γ ± = C 6 ω 2 + 2ω2 z ± 5ω4 22ω ω 2 z 2 + 8ωz 4 25ω ω4 z 32ω2 ω2 z. (54) For an isotropic trap, Γ + = 0 and Γ reduces to (50) with l = 2, corresponding to the monopole and quadrupole excitations, respectively. For ω ω z, one can show that Γ +, Γ < 0, which means that the both solutions are always unstable. As we discuss in more detail in Section VI, these unstable solutions indicate that in a trapped Bose gas, the local equilibrium distribution function (10) is not always a good starting point in solving the quantum Boltzmann equation and thus the Chapman-Enskog solution is not valid. More generally, however, we expect that there will be class of solutions of (44) which will be stable and in this case, our explicit expressions for κ and η [see (34) and (38)] will allow one to calculate the damping of these modes. A systematic study to classify which modes are stable and which are unstable would be useful. V. GENERALIZATION TO A SUPERFLUID GAS In this section, we extend the hydrodynamic equations derived in Section II to the case of a superfluid trapped gas by generalizing the results of Ref. [4] to include hydrodynamic damping. The two-fluid equations in Ref. [4] consist of equations of motion for the condensate and a set of hydrodynamic equations for the non-condensate. The latter description is only valid for in the semiclassical region k B T hω 0, gn 0. The time-dependent Hartree- Fock-Popov equation of motion for the condensate wavefunction can be rewritten in terms of 13
14 a pair of hydrodynamic equations for the condensate density n c and the superfluid velocity v s, The potential φ is defined by [4] n c t = (n cv s ), (55a) [ vs m t + 1 ] 2 v2 s = φ, (55b) φ(r, t) h2 2 [n c (r, t)] 1/2 2m[n c (r, t)] 1/2 + U 0 (r) + 2gñ(r, t) + gn c (r, t), (56) where ñ(r, t) is the non-condensate density describing the excited atoms. The first term in (56) is sometimes called the quantum pressure term and vanishes in uniform systems. It is associated with spatially varying amplitude of the Bose order parameter. The equations of motion for the non-condensate can be derived from the kinetic equation (1) for the distribution function f(r,p, t) of the excited atoms. The HF mean field in U(r, t) is now 2gn(r, t) = 2g[ñ(r, t) + n c (r, t)], where ñ(n c ) represents the local density of the non-condensate (condensate) atoms. Apart from this change, the expression for the collision integral in (2) is still applicable at finite temperatures below T BEC. Thus the hydrodynamic equations for the non-condensate can be derived using precisely the same procedure developed in Section II for T > T BEC. We obtain equations analogous to (26) in the following form: ñ t + (ñv n) = 0, ( ) mñ t + v n (57a) v nµ + P + ñ U = { [ 2η D µν 1 ]} x µ x µ x ν 3 (TrD)δ µν, (57b) ε t + ( εv n) + ( v n ) P = (κ T) + 2η [ D µν 1 3 (TrD)δ µν] 2. (57c) Here ñ, P and ε are given by the same expressions as n, P and ε in (11) and (12). Dµν is defined as in (9), where now the velocity v is the normal fluid velocity v n. The expressions for the transport coefficients given in (34) and (38) are directly applicable below T BEC for any trapping potential, apart from change in the form of the HF field in the local fugacity z. We might note that the results for κ and η by Kirkpatrick and Dorfman [10,12] for a uniform gas below T BEC have an additional factor (1 + bn c Λ 3 ) 1, where b is of order 1. This 14
15 factor arises from collisions between condensate and excited atoms, which are important at very low temperatures (T T BEC ). Our present analysis ignores this contribution. In this regard, we note that (1) and (2) are based on validity of a simple semi-classical Hartree-Fock description of excited atoms (ignoring the off-diagonal or anomalous self-energy terms). This has been shown [13] to give an excellent approximation for the thermodynamic properties of a trapped Bose gas down to much lower temperatures (T T BEC ) than in a uniform gas. If we neglect the first term in (56), (55) and (57) can be shown [4] to be equivalent to the two-fluid equations of Landau [14], except that our equations do not include any effect from the bulk viscosities. It has been shown by Kirkpatrick and Dorfman [10] for a uniform superfluid Bose gas that the bulk viscosities vanish in the finite temperature region (k B T gn 0 ) we have been considering. Our present results for a trapped Bose gas are consistent with this and moreover justify the omission of any dissipative terms in the condensate equation of motion (55b). VI. CONCLUDING REMARKS Summarizing our main results, starting from the microscopic equations of motion, we have derived a closed equation (44) for the velocity fluctuations of a trapped Bose gas, which includes damping due to hydrodynamic processes. We have also obtained explicit expressions for the shear viscosity η and the thermal conductivity κ for a trapped Bose gas, which depend on position through the local fugacity z 0. We have given a detailed derivation for a trapped gas above T BEC, and more briefly discussed the generalization of these results to a Bose-condensed gas below T BEC in Section V. For illustration, in Section IV, we used these equations to discuss some of the normal modes of oscillation for a normal Bose gas above T BEC. We showed that the monopole mode of an isotropic trap was undamped. In contrast, we found that the undamped surface mode solutions discussed in Ref. [3] are unstable solutions of (44) since Γ in (50) leads to exponential growth rather than decay. A similar result is obtained for the coupled quadrupole-monopole modes of an anisotropic trap. We note that these unstable modes are either isothermal or characterized by a spatially constant temperature fluctuations. It is the contribution R defined in (43) which makes these solutions unstable and this is entirely due 15
16 to the position-dependent quantum correction to the shear viscosity coefficient through its dependence on the local fugacity z. In the classical limit, these modes have zero damping. As soon as the first correction due to Bose statistics is included (as in (41b)), these solutions become unstable. These results show that for certain normal mode oscillations, the deviation from local equilibrium is not small, which removes the basis of the Chapman-Enskog procedure [7,11] for solving the Boltzmann equation given by (1) and (2). We recall that the entire analysis in this paper is built on the assumption that the deviation from the local distribution function f (0) (r,p, t) in (14) can be treated as a small perturbation. The precise reasons why the local equilibrium distribution is not stable for certain modes in a trap requires further study. Several authors have noted the similarity between some of the normal modes in a trapped gas above T BEC and those of a Bose condensate at T = 0 [15,16]. The fact that certain modes above T BEC are unstable when hydrodynamic damping is included emphasizes the need for further studies on the damping of the analogous condensate modes in the superfluid phase. We note that the general form of the hydrodynamic equations given in Eq. (26) can be derived more generally using conservation laws for conserved quantities in conjunction with local thermal equilibrium. This approach is developed in many textbooks [17] and can be formally extended to deal with inhomogeneous superfluids in an external potential (see, for example, Ref. [18]). This method is, of course, phenomenological in that it introduces various linear response coefficients [relating densities (thermodynamic derivatives) and currents (transport coefficients) to fields] whose evaluation requires some specific microscopic model. While very useful, this approach must be ultimately justified by a fully microscopic approach, such as used in this paper and in Refs. [3,4,7,10]. In particular, for certain kinds of oscillations, our results imply that the standard general derivation of hydrodynamic equations becomes invalid in trapped Bose gases. The advantage of our approach (see also Refs. [3,4]) is that we can investigate these problems by going back to a more microscopic description given by kinetic equations, such as (1) and (2). In conclusion, we see that hydrodynamical oscillations of a Bose gas trapped in a parabolic potential are quite unusual at a very fundamental level. Much further work is needed but we hope that our present results will encourage experimental efforts to probe 16
17 the hydrodynamic regime. ACKNOWLEDGMENTS T.N. was supported by a grant from the Japan Society for the Promotion of Science (JSPS) and A.G. by a grant from NSERC Canada. APPENDIX A: We briefly sketch the derivation of the left hand side of the kinetic equation in (16). Using (10) in the left hand side of (1), one has [ t + p ] m r r U(r, t) p f (0) (r,p, t) [ ( 1 = z t + p ) ( m r z + mv2 2k B T 2 t + p ) m r T + mu ( k B T t + p ) m r v + ] ru(r, t) u f (0) (1 + f (0) ). k B T (A1) Using the lowest-order hydrodynamic equations in (13) and the expressions for the density n in (11) and the pressure P in (12), one finds where γ 1 1 k B g T Λ 3 1/2 (z), and n t = 3n T 2T t + γk BT z z t = 3n 2T v T γk BT v z n v, z P t = 5 P T 2 T t + nk BT z z t = 5 P 2 T v T nk BT v z 5 P( v). z 3 (A5) (A2) (A3) (A4) One may combine these equations to obtain z t = v z = z ( P nk B T v 5 2 ) P T T, (A6) T t = 2 T v (v )T. 3 (A7) 17
18 The analogous equation for v/ t is given directly by (13b). Using these results in (A1), one finds that it reduces to ( t + p ) m r r U p f (0) { ( 1 mu 2 = T u T 2k B T 5P ) + m 2nk B T k B T [ ]} u (u )v u2 3 v f (0) (1 + f (0) ), (A8) where we recall u p/m v. This can be rewritten in the form shown on the left hand side of (16). 18
19 REFERENCES [1] D.S. Jin et al., Phys. Rev. Lett. 77, 420 (1996); M.-O. Mewes et al., Phys. Rev. Lett., 77, 992(1996). [2] D.S. Jin et al., Phys. Rev. Lett. 78, 764 (1997). [3] A. Griffin, W.C. Wu and S. Stringari, Phys. Rev. Lett. 78, 1838 (1997). [4] E. Zaremba, A. Griffin and T. Nikuni, Phys. Rev. Lett., submitted; cond-mat/ [5] M.R. Andrews et al., Phys. Rev. Lett. 79, 553 (1997). [6] A useful discussion of the conditions which must be satisfied to probe the hydrodynamic region is given by M. Kavoulakis, C.J. Pethick and H. Smith, cond-mat/ [7] E.A. Uehling and G.E. Uhlenbeck, Phys. Rev. 43, 552 (1933). [8] E.A. Uehling, Phys. Rev. 46, 917 (1934). [9] L.P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, N.Y., 1962), Ch. 6. [10] T.R. Kirkpatrick and J.R. Dorfman, Journ. Low. Temp. Phys. 58, 304 (1985); ibid. 399 (1985). [11] For a detailed review of the Chapman-Enskog method as applied to classical gases, see J.H. Ferziger and H.G. Kaper, Mathematical Theory of Transport Processes in Gases (North-Holland, London, 1972). [12] At T BEC (when z = 1 and n c = 0), our integral equations (22) are equivalent to Eqs. (24) in the second paper in Ref. [10], and hence our results for κ and η should coincide at T BEC with Eqs. (26) of this reference. However, our numerical factors differ slightly. [13] S. Giorgini, L. Pitaevskii and S. Stringari, Journ. Low. Temp. Phys., in press, condmat/ [14] I.M. Khalatnikov, An Introduction to the Theory of Superfluidity (Benjamin, N.Y., 1965). [15] S. Stringari, Phys. Rev. Lett. 77, 902 (1996). [16] Yu. Kagan et al., Phys. Rev. A 55, R18 (1997). [17] For example, see P.M. Chaikin and T.C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, N.Y., 1995), Ch. 8 and expecially p.460ff. [18] See, for example, T.L. Ho and V.B. Shenoy, cond-mat/
20 FIGURE CAPTIONS Fig.1: The thermal conductivity κ as a function of the fugacity z. The values are normalized to the classical gas results at z = 0. (see (41a)) Fig.2: The viscosity coefficient η as a function of the fugacity z. The values are normalized to the classical gas results at z = 0. (see (41b)) 20
21 Figure 1 κ(z)/κ(z=0) z
22 Figure 2 η(z)/η(z=0) z
arxiv:cond-mat/ v3 4 Jan 1998
Hydrodynamic damping in trapped Bose gases T. Nikuni Department of Physics, Tokyo Institute of Technology, Oh-okayama, Tokyo, 152 Japan and Department of Physics, University of Toronto, Toronto, Ontario,
More informationDynamics of a Trapped Bose-Condensed Gas at Finite Temperatures
Dynamics of a Trapped Bose-Condensed Gas at Finite Temperatures Tetsuro Nikuni, Tokyo University of Science Collaborators: Allan Griffin, University of Toronto (Canada) Eugene Zaremba, Queen s University
More informationBOSE-CONDENSED GASES AT FINITE TEMPERATURES
BOSE-CONDENSED GASES AT FINITE TEMPERATURES The discovery of Bose Einstein condensation (BEC) in trapped ultracold atomic gases in 1995 has led to an explosion of theoretical and experimental research
More informationStrongly correlated systems in atomic and condensed matter physics. Lecture notes for Physics 284 by Eugene Demler Harvard University
Strongly correlated systems in atomic and condensed matter physics Lecture notes for Physics 284 by Eugene Demler Harvard University January 25, 2011 2 Chapter 12 Collective modes in interacting Fermi
More informationThe Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs
The Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs i ( ) t Φ (r, t) = 2 2 2m + V ext(r) + g Φ (r, t) 2 Φ (r, t) (Mewes et al., 1996) 26/11/2009 Stefano Carignano 1 Contents 1 Introduction
More informationLandau damping of transverse quadrupole oscillations of an elongated Bose-Einstein condensate
PHYSICAL REVIEW A 67, 053607 2003 Landau damping of transverse quadrupole oscillations of an elongated Bose-Einstein condensate M. Guilleumas 1 and L. P. Pitaevskii 2,3 1 Departament d Estructura i Constituents
More informationThe Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs
The Gross-Pitaevskii Equation and the Hydrodynamic Expansion of BECs RHI seminar Pascal Büscher i ( t Φ (r, t) = 2 2 ) 2m + V ext(r) + g Φ (r, t) 2 Φ (r, t) 27 Nov 2008 RHI seminar Pascal Büscher 1 (Stamper-Kurn
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 24 Jul 2001
arxiv:cond-mat/010751v1 [cond-mat.stat-mech] 4 Jul 001 Beyond the Thomas-Fermi Approximation for Nonlinear Dynamics of Trapped Bose-Condensed Gases Alexander L. Zubarev and Yeong E. Kim Department of Physics,
More informationLandau Theory of Fermi Liquids : Equilibrium Properties
Quantum Liquids LECTURE I-II Landau Theory of Fermi Liquids : Phenomenology and Microscopic Foundations LECTURE III Superfluidity. Bogoliubov theory. Bose-Einstein condensation. LECTURE IV Luttinger Liquids.
More informationarxiv:cond-mat/ v1 28 Feb 1996
Bose-Einstein Condensation of Atoms in a Trap T.T. Chou, Chen Ning Yang 2 and L.H. Yu 3 arxiv:cond-mat/960253v 28 Feb 996 Department of Physics, University of Georgia, Athens, Georgia 30602 2 Institute
More informationRoton Mode in Dipolar Bose-Einstein Condensates
Roton Mode in Dipolar Bose-Einstein Condensates Sandeep Indian Institute of Science Department of Physics, Bangalore March 14, 2013 BECs vs Dipolar Bose-Einstein Condensates Although quantum gases are
More informationarxiv:cond-mat/ v1 2 Apr 1997
Thermodynamics of a trapped Bose-condensed gas S. Giorgini 1,2, L. P. Pitaevskii 2,3,4 and S. Stringari 2 1 ECT, Villa Tambosi, Strada delle Tabarelle 286, I-38050 Villazzano, Italy 2 Dipartimento di Fisica,
More informationShock waves in the unitary Fermi gas
Shock waves in the unitary Fermi gas Luca Salasnich Dipartimento di Fisica e Astronomia Galileo Galilei, Università di Padova Banff, May 205 Collaboration with: Francesco Ancilotto and Flavio Toigo Summary.
More information5. Gross-Pitaevskii theory
5. Gross-Pitaevskii theory Outline N noninteracting bosons N interacting bosons, many-body Hamiltonien Mean-field approximation, order parameter Gross-Pitaevskii equation Collapse for attractive interaction
More informationLow-dimensional Bose gases Part 1: BEC and interactions
Low-dimensional Bose gases Part 1: BEC and interactions Hélène Perrin Laboratoire de physique des lasers, CNRS-Université Paris Nord Photonic, Atomic and Solid State Quantum Systems Vienna, 2009 Introduction
More informationAnisotropic Spin Diffusion in Trapped Boltzmann Gases
Anisotropic Spin Diffusion in Trapped Boltzmann Gases W. J. Mullin 1 and R. J. Ragan 2 1 Physics Department, Hasbrouck Laboratory University of Massachusetts, Amherst, MA 01003 USA 2 Physics Department,
More informationSuperfluidity in bosonic systems
Superfluidity in bosonic systems Rico Pires PI Uni Heidelberg Outline Strongly coupled quantum fluids 2.1 Dilute Bose gases 2.2 Liquid Helium Wieman/Cornell A. Leitner, from wikimedia When are quantum
More informationChapter 14. Ideal Bose gas Equation of state
Chapter 14 Ideal Bose gas In this chapter, we shall study the thermodynamic properties of a gas of non-interacting bosons. We will show that the symmetrization of the wavefunction due to the indistinguishability
More informationINTERACTING BOSE GAS AND QUANTUM DEPLETION
922 INTERACTING BOSE GAS AND QUANTUM DEPLETION Chelagat, I., *Tanui, P.K., Khanna, K.M.,Tonui, J.K., Murunga G.S.W., Chelimo L.S.,Sirma K. K., Cheruiyot W.K. &Masinde F. W. Department of Physics, University
More informationMany-body physics 2: Homework 8
Last update: 215.1.31 Many-body physics 2: Homework 8 1. (1 pts) Ideal quantum gases (a)foranidealquantumgas,showthatthegrandpartitionfunctionz G = Tre β(ĥ µ ˆN) is given by { [ ] 1 Z G = i=1 for bosons,
More informationMonte Carlo Simulation of Bose Einstein Condensation in Traps
Monte Carlo Simulation of Bose Einstein Condensation in Traps J. L. DuBois, H. R. Glyde Department of Physics and Astronomy, University of Delaware Newark, Delaware 19716, USA 1. INTRODUCTION In this paper
More informationPublication II American Physical Society. Reprinted with permission.
II Publication II Tomoya Isoshima, Jukka Huhtamäki, and Martti M. Salomaa, Precessional motion of a vortex in a finite-temperature Bose-Einstein condensate, Phys. Rev. A69, 063601 (2004). 2004 American
More informationarxiv:cond-mat/ v1 13 Mar 1998
Stability of Solution of the Nonlinear Schrödinger Equation for the Bose-Einstein Condensation arxiv:cond-mat/9803174v1 13 Mar 1998 Yeong E. Kim and Alexander L. Zubarev Department of Physics, Purdue University
More information1 Fluctuations of the number of particles in a Bose-Einstein condensate
Exam of Quantum Fluids M1 ICFP 217-218 Alice Sinatra and Alexander Evrard The exam consists of two independant exercises. The duration is 3 hours. 1 Fluctuations of the number of particles in a Bose-Einstein
More informationGreen s Function Approach to B. E. Condensation in the alkali atoms
arxiv:quant-ph/010048v1 8 Oct 00 Green s Function Approach to B. E. Condensation in the alkali atoms Mahendra Sinha Roy October 9, 018 Department of Physics, Presidency College, Calcutta 700 073,India.
More informationBCS Pairing Dynamics. ShengQuan Zhou. Dec.10, 2006, Physics Department, University of Illinois
BCS Pairing Dynamics 1 ShengQuan Zhou Dec.10, 2006, Physics Department, University of Illinois Abstract. Experimental control over inter-atomic interactions by adjusting external parameters is discussed.
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 2 Aug 2004
Ground state energy of a homogeneous Bose-Einstein condensate beyond Bogoliubov Christoph Weiss and André Eckardt Institut für Physik, Carl von Ossietzky Universität, D-6 Oldenburg, Germany (Dated: November
More informationUltracold Fermi and Bose Gases and Spinless Bose Charged Sound Particles
October, 011 PROGRESS IN PHYSICS olume 4 Ultracold Fermi Bose Gases Spinless Bose Charged Sound Particles ahan N. Minasyan alentin N. Samoylov Scientific Center of Applied Research, JINR, Dubna, 141980,
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynamics As a first approximation, we shall assume that in local equilibrium, the density f 1 at each point in space can be represented as in eq.(iii.56), i.e. [ ] p m q, t)) f
More informationInauguration Meeting & Celebration of Lev Pitaevskii s 70 th Birthday. Bogoliubov excitations. with and without an optical lattice.
Inauguration Meeting & Celebration of Lev Pitaevskii s 7 th Birthday Bogoliubov excitations with and without an optical lattice Chiara Menotti OUTLINE OF THE TALK Bogoliubov theory: uniform system harmonic
More informationBose-Einstein Condensation
Bose-Einstein Condensation Kim-Louis Simmoteit June 2, 28 Contents Introduction 2 Condensation of Trapped Ideal Bose Gas 2 2. Trapped Bose Gas........................ 2 2.2 Phase Transition.........................
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationSuperfluidity and superconductivity. IHP, Paris, May 7 and 9, 2007
Superfluidity and superconductivity. IHP, Paris, May 7 and 9, 2007 L.P. Pitaevskii Dipartimento di Fisica, Universita di Trento, INFM BEC CNR,Trento, Italy; Kapitza Institute for Physical Problems, ul.
More informationarxiv: v1 [cond-mat.quant-gas] 18 Sep 2009
Sound propagation in a Bose-Einstein condensate at finite temperatures arxiv:0909.3455v1 [cond-mat.quant-gas] 18 Sep 2009 R. Meppelink, S. B. Koller, and P. van der Straten 1 1 Atom Optics and Ultrafast
More informationarxiv: v1 [cond-mat.stat-mech] 30 Aug 2007
Journal of Low Temperature Physics manuscript No. (will be inserted by the editor) arxiv:0708.414v1 [cond-mat.stat-mech] 30 Aug 007 Luca Salasnich 1 and Flavio Toigo Shell Effects in the First Sound Velocity
More informationPAPER 84 QUANTUM FLUIDS
MATHEMATICAL TRIPOS Part III Wednesday 6 June 2007 9.00 to 11.00 PAPER 84 QUANTUM FLUIDS Attempt TWO questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY REQUIREMENTS
More informationExcitations and dynamics of a two-component Bose-Einstein condensate in 1D
Author: Navarro Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Bruno Juliá Díaz Abstract: We study different solutions and their stability for a two component
More informationGas flow around a longitudinally oscillating plate at arbitrary ratio of collision frequency to oscillation frequency
Gas flow around a longitudinally oscillating plate at arbitrary ratio of collision frequency to oscillation frequency Feli Sharipov and Denize Kalempa Departamento de Física, Universidade Federal do Paraná
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationarxiv: v3 [gr-qc] 30 Mar 2009
THE JEANS MECHANISM AND BULK-VISCOSITY EFFECTS Nakia Carlevaro a, b and Giovanni Montani b, c, d, e a Department of Physics, Polo Scientifico Università degli Studi di Firenze, INFN Section of Florence,
More informationCRITICAL ATOM NUMBER OF A HARMONICALLY TRAPPED 87 Rb BOSE GAS AT DIFFERENT TEMPERATURES
CRITICAL ATOM UMBER OF A HARMOICALLY TRAPPED 87 Rb BOSE GAS AT DIFFERET TEMPERATURES AHMED S. HASSA, AZZA M. EL-BADRY Department of Physics, Faculty of Science, El Minia University, El Minia, Egypt E-mail:
More informationRelativistic magnetotransport in graphene
Relativistic magnetotransport in graphene Markus Müller in collaboration with Lars Fritz (Harvard) Subir Sachdev (Harvard) Jörg Schmalian (Iowa) Landau Memorial Conference June 6, 008 Outline Relativistic
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 25 Feb 2006
1 arxiv:cond-mat/0602600v1 cond-mat.stat-mech 25 Feb 2006 The equation of state for two-dimensional hard-sphere gases: Hard-sphere gases as ideal gases with multi-core boundaries Wu-Sheng Dai 1,2 and Mi
More informationPhysics 112 The Classical Ideal Gas
Physics 112 The Classical Ideal Gas Peter Young (Dated: February 6, 2012) We will obtain the equation of state and other properties, such as energy and entropy, of the classical ideal gas. We will start
More informationHIGH FRICTION LIMIT OF THE KRAMERS EQUATION : THE MULTIPLE TIME SCALE APPROACH. Lydéric Bocquet
HIGH FRICTION LIMIT OF THE KRAMERS EQUATION : THE MULTIPLE TIME SCALE APPROACH Lydéric Bocquet arxiv:cond-mat/9605186v1 30 May 1996 Laboratoire de Physique, Ecole Normale Supérieure de Lyon (URA CNRS 1325),
More informationFrom laser cooling to BEC First experiments of superfluid hydrodynamics
From laser cooling to BEC First experiments of superfluid hydrodynamics Alice Sinatra Quantum Fluids course - Complement 1 2013-2014 Plan 1 COOLING AND TRAPPING 2 CONDENSATION 3 NON-LINEAR PHYSICS AND
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Sep 2005
1 arxiv:cond-mat/0509103v1 [cond-mat.stat-mech] 5 Sep 2005 Normal and Anomalous Averages for Systems with Bose-Einstein Condensate V.I. Yukalov 1,2 and E.P. Yukalova 1,3 1 Institut für Theoretische Physik,
More informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationQuantum Monte Carlo Simulations of Exciton Condensates
Quantum Monte Carlo Simulations of Exciton Condensates J. Shumway a and D. M. Ceperley b a Dept. of Physics and Astronomy, Arizona State University, Tempe, AZ 8583 b Dept. of Physics, University of Illinois,
More information(Super) Fluid Dynamics. Thomas Schaefer, North Carolina State University
(Super) Fluid Dynamics Thomas Schaefer, North Carolina State University Hydrodynamics Hydrodynamics (undergraduate version): Newton s law for continuous, deformable media. Fluids: Gases, liquids, plasmas,...
More informationGaussian fluctuations in an ideal bose-gas a simple model
Gaussian fluctuations in an ideal bose-gas a simple model A Petrova,, O Nedopekin, D Tayurskii and Q A Wang LUNAM Universite, ISMANS, Laboratoire de Physique Statistique et Systeme Complexes, 44, Avenue
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationInteraction between atoms
Interaction between atoms MICHA SCHILLING HAUPTSEMINAR: PHYSIK DER KALTEN GASE INSTITUT FÜR THEORETISCHE PHYSIK III UNIVERSITÄT STUTTGART 23.04.2013 Outline 2 Scattering theory slow particles / s-wave
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationTransport coefficients from Kinetic Theory: Bulk viscosity, Diffusion, Thermal conductivity. Debarati Chatterjee
Transport coefficients from Kinetic Theory: Bulk viscosity, Diffusion, Thermal conductivity Debarati Chatterjee Recap: Hydrodynamics of nearly perfect fluids Hydrodynamics: correlation functions at low
More informationarxiv: v1 [nucl-th] 18 Apr 2019
Prepared for submission to JHEP Emergence of hydrodynamical behavior in expanding quark-gluon plasmas arxiv:1904.08677v1 [nucl-th] 18 Apr 2019 Jean-Paul Blaizot a and Li Yan b a Institut de Physique Théorique,
More informationThe Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by
More informationHolographic zero sound at finite temperature
Holographic zero sound at finite temperature Richard Davison and Andrei Starinets Rudolf Peierls Centre for Theoretical Physics Oxford University arxiv: 1109.XXXX [hep-th] Holographic Duality and Condensed
More informationIntroduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC FERMI GASES
1 INTERNATIONAL SCHOOL OF PHYSICS "ENRICO FERMI" Varenna, July 1st - July 11 th 2008 " QUANTUM COHERENCE IN SOLID STATE SYSTEMS " Introduction to Bose-Einstein condensation 4. STRONGLY INTERACTING ATOMIC
More informationSupplementary Figure 3: Interaction effects in the proposed state preparation with Bloch oscillations. The numerical results are obtained by
Supplementary Figure : Bandstructure of the spin-dependent hexagonal lattice. The lattice depth used here is V 0 = E rec, E rec the single photon recoil energy. In a and b, we choose the spin dependence
More informationSpontaneous Symmetry Breaking in Bose-Einstein Condensates
The 10th US-Japan Joint Seminar Spontaneous Symmetry Breaking in Bose-Einstein Condensates Masahito UEDA Tokyo Institute of Technology, ERATO, JST collaborators Yuki Kawaguchi (Tokyo Institute of Technology)
More informationThe First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
Commun. Theor. Phys. (Beijing, China 35 (2 pp. 42 46 c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationarxiv:hep-ph/ v1 30 Nov 2001
CHAPMAN-ENSKOG EXPANSION OF THE BOLTZMANN EQUATION AND ITS DIAGRAMMATIC INTERPRETATION arxiv:hep-ph/0446v 30 Nov 00 M.E. CARRINGTON A,B, HOU DEFU A,B,C AND R. KOBES B,D a Department of Physics, Brandon
More informationV. FERMI LIQUID THEORY
V. FERMI LIQUID THEORY Markus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus/cours.html (Dated: March 30, 2010)
More informationCondensed matter theory Lecture notes and problem sets 2012/2013
Condensed matter theory Lecture notes and problem sets 2012/2013 Dmitri Ivanov Recommended books and lecture notes: [AM] N. W. Ashcroft and N. D. Mermin, Solid State Physics. [Mar] M. P. Marder, Condensed
More informationPhase transitions and critical phenomena
Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (
More informationSecond Quantization Method for Bosons
Second Quantization Method for Bosons Hartree-Fock-based methods cannot describe the effects of the classical image potential (cf. fig. 1) because HF is a mean-field theory. DFF-LDA is not able either
More informationarxiv:cond-mat/ v1 14 May 1997
Semiclassical wave functions and energy levels of Bose-condensed gases in spherically symmetric traps András Csordás 1, Robert Graham, Péter Szépfalusy 3 1 Research Group for Statistical Physics of the
More informationCHAPMAN-ENSKOG EXPANSION OF THE BOLTZMANN EQUATION AND ITS DIAGRAMMATIC INTERPRETATION
CHAPMAN-ENSKOG EXPANSION OF THE BOLTZMANN EQUATION AND ITS DIAGRAMMATIC INTERPRETATION M.E. CARRINGTON A,B,HOUDEFU A,B,C AND R. KOBES B,D a Department of Physics, Brandon University, Brandon, MB,R7A 6A9
More informationarxiv:cond-mat/ v1 4 Aug 2003
Conductivity of thermally fluctuating superconductors in two dimensions Subir Sachdev arxiv:cond-mat/0308063 v1 4 Aug 2003 Abstract Department of Physics, Yale University, P.O. Box 208120, New Haven CT
More informationin-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators
(by A. A. Shanenko) in-medium wave functions in-medium pair-wave functions and spatial pair particle correlations momentum condensation and ODLRO (off-diagonal long range order) U(1) symmetry breaking
More informationSecond sound and the superfluid fraction in a resonantly interacting Fermi gas
Second sound and the superfluid fraction in a resonantly interacting Fermi gas Meng Khoon Tey Tsinghua University China Workshop on Probing and Understanding Exotic Superconductors and Superfluids Trieste,
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More informationarxiv: v1 [cond-mat.quant-gas] 7 Aug 2016
APS preprint Temperature and interaction dependence of the moment of inertia of a rotating condensate boson gas arxiv:608.059v [cond-mat.quant-gas] 7 Aug 06 Ahmed S. Hassan, Azza M. El-Badry, and Shemi
More informationBose Gases, Bose Einstein Condensation, and the Bogoliubov Approximation
Bose Gases, Bose Einstein Condensation, and the Bogoliubov Approximation Robert Seiringer IST Austria Mathematical Horizons for Quantum Physics IMS Singapore, September 18, 2013 R. Seiringer Bose Gases,
More informationBogoliubov theory of disordered Bose-Einstein condensates
Bogoliubov theory of disordered Bose-Einstein condensates Christopher Gaul Universidad Complutense de Madrid BENASQUE 2012 DISORDER Bogoliubov theory of disordered Bose-Einstein condensates Abstract The
More informationBEC Vortex Matter. Aaron Sup October 6, Advisor: Dr. Charles Hanna, Department of Physics, Boise State University
BEC Vortex Matter Aaron Sup October 6, 006 Advisor: Dr. Charles Hanna, Department of Physics, Boise State University 1 Outline 1. Bosons: what are they?. Bose-Einstein Condensation (BEC) 3. Vortex Formation:
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationAnswer TWO of the three questions. Please indicate on the first page which questions you have answered.
STATISTICAL MECHANICS June 17, 2010 Answer TWO of the three questions. Please indicate on the first page which questions you have answered. Some information: Boltzmann s constant, kb = 1.38 X 10-23 J/K
More informationThe Gross-Pitaevskii Equation A Non-Linear Schrödinger Equation
The Gross-Pitaevskii Equation A Non-Linear Schrödinger Equation Alan Aversa December 29, 2011 Abstract Summary: The Gross-Pitaevskii equation, also called the non-linear Schrödinger equation, describes
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
More informationTheory Seminar Uni Marburg. Bose-Einstein Condensation and correlations in magnon systems
Theory Seminar Uni Marburg 11 November, 2010 Bose-Einstein Condensation and correlations in magnon systems Peter Kopietz, Universität Frankfurt 1.) Bose-Einstein condensation 2.) Interacting magnons in
More information1. Thomas-Fermi method
1. Thomas-Fermi method We consider a system of N electrons in a stationary state, that would obey the stationary Schrödinger equation: h i m + 1 v(r i,r j ) Ψ(r 1,...,r N ) = E i Ψ(r 1,...,r N ). (1.1)
More informationSemiclassical theory of the quasi two-dimensional trapped Bose gas
May 8 EPL, 8 (8) 3 doi:.9/95-575/8/3 www.epljournal.org Semiclassical theory of the quasi two-dimensional trapped Bose gas Markus Holzmann,, Maguelonne Chevallier 3 and Werner Krauth 3(a) LPTMC, Université
More informationSuperfluidity of a 2D Bose gas (arxiv: v1)
Superfluidity of a 2D Bose gas (arxiv:1205.4536v1) Christof Weitenberg, Rémi Desbuquois, Lauriane Chomaz, Tarik Yefsah, Julian Leonard, Jérôme Beugnon, Jean Dalibard Trieste 18.07.2012 Phase transitions
More informationWe can then linearize the Heisenberg equation for in the small quantity obtaining a set of linear coupled equations for and :
Wednesday, April 23, 2014 9:37 PM Excitations in a Bose condensate So far: basic understanding of the ground state wavefunction for a Bose-Einstein condensate; We need to know: elementary excitations in
More informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationThe Ginzburg-Landau Theory
The Ginzburg-Landau Theory A normal metal s electrical conductivity can be pictured with an electron gas with some scattering off phonons, the quanta of lattice vibrations Thermal energy is also carried
More informationarxiv: v1 [nucl-th] 9 Jun 2008
Dissipative effects from transport and viscous hydrodynamics arxiv:0806.1367v1 [nucl-th] 9 Jun 2008 1. Introduction Denes Molnar 1,2 and Pasi Huovinen 1 1 Purdue University, Physics Department, 525 Northwestern
More informationPhysics Nov Bose-Einstein Gases
Physics 3 3-Nov-24 8- Bose-Einstein Gases An amazing thing happens if we consider a gas of non-interacting bosons. For sufficiently low temperatures, essentially all the particles are in the same state
More informationI. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS
I. BASICS OF STATISTICAL MECHANICS AND QUANTUM MECHANICS Marus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris marus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/marus (Dated:
More informationStanding Sound Waves in a Bose-Einstein Condensate
Standing Sound Waves in a Bose-Einstein Condensate G.A.R. van Dalum Debye Institute - Universiteit Utrecht Supervisors: P. van der Straten, A. Groot and P.C. Bons July 17, 2012 Abstract We study two different
More informationErrata 1. p. 5 The third line from the end should read one of the four rows... not one of the three rows.
Errata 1 Front inside cover: e 2 /4πɛ 0 should be 1.44 10 7 ev-cm. h/e 2 should be 25800 Ω. p. 5 The third line from the end should read one of the four rows... not one of the three rows. p. 8 The eigenstate
More informationDrag force and superfluidity in the supersolid striped phase of a spin-orbit-coupled Bose gas
/ 6 Drag force and superfluidity in the supersolid striped phase of a spin-orbit-coupled Bose gas Giovanni Italo Martone with G. V. Shlyapnikov Worhshop on Exploring Nuclear Physics with Ultracold Atoms
More informationThe state of a quantum ideal gas is uniquely specified by the occupancy of singleparticle
Ideal Bose gas The state of a quantum ideal gas is uniquely specified by the occupancy of singleparticle states, the set {n α } where α denotes the quantum numbers of a singleparticles state such as k
More informationThe Superfluid-Insulator transition
The Superfluid-Insulator transition Boson Hubbard model M.P. A. Fisher, P.B. Weichmann, G. Grinstein, and D.S. Fisher, Phys. Rev. B 40, 546 (1989). Superfluid-insulator transition Ultracold 87 Rb atoms
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationThe non-interacting Bose gas
Chapter The non-interacting Bose gas Learning goals What is a Bose-Einstein condensate and why does it form? What determines the critical temperature and the condensate fraction? What changes for trapped
More information