there are several excellent review articles on trapped bose - condensed atomic gases written at a relatively introductory level @xcite , with emphasis on the dynamics of the condensate at @xmath0 . in the present lectures ,
i will concentrate on the interplay between the condensate and non - condensate components .
this topic requires a more sophisticated analysis based on the concepts and methods of many body theory .
my lectures will attempt to give a very basic introduction to this kind of approach .
the primary audience i have in mind are graduate students and postdocs coming from atomic and laser physics , rather than condensed matter theory .
however , i hope the experts will also learn something from these lectures . in his lectures , fetter @xcite has given a detailed review of excitations in a trapped dilute bose gas at @xmath0 .
he shows that a convenient way of discussing these excitations is to start from the time - dependent gross - pitaevskii ( gp ) equation for the macroscopic wavefunction @xmath1 .
linearizing around the equilibrium value @xmath2 , one finds of the macroscopic wavefunction are given by the well - known bogoliubov coupled equations . for a uniform bose - condensed gas , the excitation frequencies are those first discussed by bogoliubov in 1947 @xcite .    at @xmath0
, one can assume that all the atoms are in the bose condensate described by @xmath1 .
in contrast , in these lectures , my main topic will be to review our understanding of what excitations are in a bose - condensed gas at _
finite _ temperatures , when there is a large number of atoms in the non - condensate ( in trapped gases , the non - condensate is often referred to as the `` thermal cloud '' ) .
i will make some contact with the ideas discussed in the lectures by burnett @xcite .    the first half of these lectures ( sections  [ sec : elementary excitations and density fluctuations in normal systems ] - [ sec : why are excitations in a bose fluid so interesting ] ) deals with excitations whose very existence depends on self - consistent mean fields ( of various kinds ! ) , rather than on the collisions between atoms .
in the standard language developed in condensed matter physics in the 1960 s , this means the excitations are in the `` collisionless region '' .    the second half of the lectures ( sections  [ sec : hydrodynamic oscillations in a trapped bose gas ] - [ sec : first and second sound in a uniform bose - condensed gas ] ) deals with excitations in the collision - dominated hydrodynamic region .
i review the two - fluid hydrodynamic equations such as given by landau @xcite , generalized to include a trap .
i give an explicit microscopic derivation of such two - fluid equations in a trapped bose gas .
this extends recent work of zaremba , griffin and nikuni ( zgn ) @xcite to include the case when the condensate is not yet in local equilibrium with the non - condensate atoms .
the type of questions i want to address in these lectures include :    1 .
what is the difference between an elementary excitation and a collective mode ?
2 .   can we
_ isolate _ the dynamical role of the condensate on the nature of the excitations ? 3 .   at @xmath0 , with a pure condensate , excitations in a gas must be oscillations of the condensate .
in contrast , above @xmath3 , the excitations are not related to a condensate . what happens to an excitation as we go from @xmath0 to @xmath4 .
how does the excitation get rid of its `` condensate '' dressing ?
what is the essential physics behind the different mean - field theories of excitations which have been discussed in the recent literature : gross - pitaevskii , hartree - fock - popov , hartree - fock - bogoluibov ?
what is the physics behind the dreaded beliaev second - order approximation ?
this is the first approximation which includes damping of the elementary excitations even at @xmath0  @xcite .
why are the excitations and collective modes in a bose - condensed system uniquely interesting , compared to all other many body systems ?
the key reason is , of course , that the condensate couples and hybridizes single - particle excitations with density fluctuations . above @xmath3
, these two excitation branches are uncoupled .
once this is understood , one sees why more detailed and systematic experimental studies of excitations in trapped atomic gases at finite temperature are of great importance .    to give a careful discussion of all these questions
would need  10 lectures . in these  3 lectures
, i will give a speeded - up version .
i will often use a _ uniform _ weakly interacting bose gas to illustrate the structure of the theory .
while i will only sketch the math , i will still try to give a flavour of what is involved .
the approach i will use to discuss these questions is based on the field - theoretic formulation of a bose - condensed system of particles .
as i review elsewhere in this volume  @xcite , this powerful formalism was introduced by beliaev  @xcite in 1957 and extensively developed in the golden period : 1958 - 1965 .
i will introduce this formalism in a very schematic manner - but even experimentalists will find it useful to know some of the  language " used in this approach .
there are other methods to deal with collective modes in bose gases but they are not as useful at _ isolating _ the dynamical role of the condensate , or dealing with finite temperatures .
while i will always have trapped atomic bose gases in mind , much of the general theory  @xcite is valid for any bose - condensed fluid ( gas or liquid ) .
thus i will often make references to superfluid @xmath5he , pointing out similarities with bose gases .
we work with quantum field operators : @xmath6 these operators satisfy the usual bose commutation relations , such as @xmath7=\delta ( { \bf r }   - { \bf r'})$ ] . of course
, if we have several hyperfine atomic states , then we have different field operators @xmath8 , where @xmath9 is the hyperfine state label ( the analogue of a spin label ) .
all observables can be written in terms of these quantum field operators .
two important examples are the local density operator : @xmath10 and the hamiltonian @xmath11 { \hat\psi}({\bf r } ) + { 1 \over   2 } g \int d{\bf r } \hat{\psi}^+({\bf r } ) \hat{\psi}^+({\bf r } )   \hat{\psi}({\bf r } ) \hat{\psi}({\bf r } ) . \label{eq : elementary excitations3}\ ] ] in these lectures , the two - particle interaction will always be approximated by a @xmath12-wave scattering length , appropriate to a dilute bose gas at very low temperatures @xcite .
thus @xmath13 , with @xmath14 ( we set @xmath15 throughout this article ) .    what we want to calculate are various kinds of correlation functions involving different local operators : @xmath16 in these tutorial lectures
, i will not bother to distinguish between the various kinds of correlation functions ( time - ordered , retarded , etc ) or even the difference between correlation functions and response functions .
these are discussed in all the standard textbooks on many body theory @xcite but are only important when one is doing systematic calculations .
a convenient summary of the formalism is given in a recent review article on homogeneous weakly interacting bose gases @xcite .
the density response function @xmath17 involves four quantum field operators and is an example of a two - particle green s function @xmath18 .
one can measure @xmath19 by coupling a weak external field to the local density in the system @xmath20 linear response theory gives @xmath21 for the density response ( see chapters 5 and 9 of ref.@xcite and chapter 2 of ref.@xcite ) . for a uniform system with @xmath22
, we have @xmath23 fourier transforming @xmath24 and @xmath25 , the linear response equation in ( [ eq : elementary excitations7 ] ) reduces to @xmath26 if @xmath27 has a pole at @xmath28 , ( for further details , see section  [ sec : density fluctuation spectrum in the mean - field approximation ] ) then when @xmath29 and @xmath30 of the external potential satisfy @xmath28 , @xmath31 as given by ( [ eq : elementary excitations8 ] ) can be very large even though @xmath32 is small .
we note that @xmath33 and thus clearly the pole at @xmath28 is the signature of an oscillating density fluctuation .
the basic correlation function in the field - theoretic approach is given by @xmath34 this single - particle green s function @xmath35 involves two quantum field operators .
it describes creating an atom at @xmath36 , let it propagate through the system to @xmath37 and then destroying the atom .
all other higher - order correlation functions such as @xmath38 can be constructed out of combinations of @xmath39 .
for example , the _ lowest order _ contributions to the density response function in ( [ eq : elementary excitations5 ] ) are @xmath40 as we discuss in section  [ sec : green s function formulation ] , the terms in the last line of ( [ eq : elementary excitations11 ] ) vanish in a normal bose system but are finite for @xmath41 .
the poles of @xmath35 correspond to what are called elementary excitations ( or quasiparticles ) . in a uniform system
, we have @xmath42 one can show @xcite that these single - particle excitations determine the thermodynamic properties of interacting systems . however , it is very difficult to directly measure the spectrum of @xmath39 since one needs an external field which couples to @xmath43 , ie , an atom reservoir .
later we will see that what makes a bose - condensed system unique is that we can easily access @xmath39 as a result of the effects of the bose condensate .
finally we introduce the key idea of a single - particle self - energy through dyson s equation : @xmath44 where @xmath35 is the interacting single - particle green s function , @xmath45 is the non - interacting single - particle green s function and all effects of the two - particle interactions are contained in the self - energy function @xmath46 . in a uniform system , we can fourier transform this dyson equation to give @xmath47 where @xmath48 .
this is easily solved to give @xmath49}.\label{eq : elementary excitations15}\ ] ] thus we see that @xmath35 may have a single - particle pole at the quasiparticle energy @xmath50 in general , @xmath51 , where @xmath52 describes the damping of the single - particle excitations .
field - theoretic calculations @xcite involve a systemmatic ( diagrammatic ) procedure to calculate @xmath53 and from this to obtain @xmath54 .
we note that the self - energy @xmath55 is by definition highly non - perturbative . as an illustration ,
let us consider the self - energy to first order in @xmath56 .
this hartree - fock approximation is shown in fig.1 . in our @xmath12-wave approximation ,
the total self - energy is simply    @xmath57    and therefore ( [ eq : elementary excitations15 ] ) gives @xmath58}. \label{eq : elementary excitations18}\ ] ] hence the normal hf excitation energy has the dispersion relation @xmath59
the simplest approximation for the density response function introduced in section  [ sec : elementary excitations and density fluctuations in normal systems ] is @xcite @xmath60 in a uniform bose system , the fourier transform of this gives @xmath61\over \omega_1 - \omega_2 - \omega } a({\bf k } , \omega_1 ) a({\bf k } - { \bf   q } , \omega_2 ) .
\label{eq : density fluc2}\ ] ] here @xmath62 is the bose distribution function and @xmath63 single - particle spectral density @xmath64 .
if we use @xmath65 as given in ( [ eq : elementary excitations18 ] ) , we find ( [ eq : density fluc2 ] ) reduces to @xmath66 we note that @xmath67 has a continuum of poles given by @xmath68 .
it is easy to understand the physics which gives rise to this `` ideal - gas '' spectrum .
the fourier transform of the local density operator in ( [ eq : elementary excitations2 ] ) is given by ( for a uniform system ) @xmath69 clearly @xmath70 creates a `` particle - hole '' density fluctuation with the following features : @xmath71 the spectrum ( [ eq : density fluc3 ] ) describes a broad incoherent superposition of particle - hole states @xcite and is not a true collective mode , such as we discuss next .
the mean field approximation ( mfa ) for the density response is also called other names : scf ( self - consistent field ) , rpa ( randon phase approximation ) , but all involve the same physics .
the mfa was introduced by bohm and pines in the period 1951 - 1953 , work which had a pivotal effect in our understanding collective effects in all many - particle systems .
we recall the linear response expression in ( [ eq : elementary excitations7 ] ) , where @xmath72 is the _ full _ density response function for interacting bose gas .
if we introduce the self - consistent hartree mean - field : @xmath73
then we can approximate the linear response equation in ( [ eq : elementary excitations7 ] ) by @xmath74 .
\label{eq : density fluc7}\ ] ] the system is assumed to respond as if the atoms propagate independently ( as described by @xmath75 but are moving in an effective field @xmath76 . in a uniform system ,
the fourier transform of ( [ eq : density fluc7 ] ) is @xmath77,\label{eq : density fluc8}\ ] ] which gives the well - known mfa expression for @xmath72 @xcite : @xmath78 this result for @xmath79 may have new poles given by zeros of the denominator , @xmath80 assuming that it is distinguisable from the incoherent ideal gas density fluctuation spectrum given by ( [ eq : density fluc3 ] ) , this pole at @xmath81 is called a zero sound mode ( a plasmon in charged systems ) . this language was first introduced in 1957 by landau in fermi systems @xcite but the concept is generally applicable in any interacting many body system .
physically , it is clear that zero sound is a `` collisionless '' density oscillation arising from dynamic self - consistent mean fields .
we make a few comments about such `` zero sound '' collective modes :    1 .
@xmath82 and @xmath83 ( the poles of @xmath54 ) are _ both _ states of an interacting many - particle system .
however the collective mode disappears if there are no interactions , while the single - particle excitations still exist in a non - interacting gas .
because of the low density , dynamic mean fields in normal systems are too weak to allow the existence of a well - defined ( weakly - damped ) zero sound mode in the bose gases of current interest ( for @xmath84 ) . if @xmath85 , we can then approximate @xmath86 .
the situation is quite different for @xmath87 , when a coherent mean - field due to the condensate is present .
one _ also _ expects a collective pole to appear in @xmath88 in the collision - dominated hydrodynamic region .
however this sound wave pole is induced by rapid collisions producing local equilibrium , as we discuss in sections  [ sec : hydrodynamic oscillations in a trapped bose gas ] - [ sec : first and second sound in a uniform bose - condensed gas ] .
ordinary sound is not described by the mfa discussed above , ie , it is not the result of dynamic mean fields .
all the formalism we have been discussing for @xmath89 can be extended in a natural way to bose - condensed systems @xcite making use of the fundamental decomposition which separates out the condensate and non - condensate parts of the quantum fields : @xmath90 here the average is over a restricted ensemble @xcite consistent with @xmath91 .
the most profound way of doing this is to add a symmetry - breaking field @xmath92 \label{eq : green function2}\ ] ] and work with @xmath93 , taking the limit @xmath94 at the end .
this gives the system a `` hunting license '' @xcite to have finite value of @xmath95 and one finds @xmath96 a key point is that if @xmath97 , then a direct consequence is the existence of `` anomalous '' or `` off - diagonal '' propagators @xcite @xmath98 these describe new condensate - induced correlations between non - condensate atoms at different space - time points @xmath99 and @xmath100 . in a sense
, these anomalous correlation functions are as important as the macroscopic wavefunction @xmath101 .
we also note that @xmath102 is the `` pair '' function that burnett  @xcite    discusses in detail , using a different formalism .
> from now on , we leave the symmetry - breaking label on the averages implicit .    in addition , if @xmath103 we find that correlation functions involving three non - condensate field operators can be finite @xmath104 in particular , @xmath105 .
this describes the condensate - induced coupling of non - condensate density fluctuations @xmath106 and the single - particle field fluctuations ( see section  [ sec : why are excitations in a bose fluid so interesting ] ) .    clearly one has to work with a @xmath107 matrix single - particle propagator @xmath35 when @xmath91 , namely @xmath108 where we have introduced a @xmath107 matrix green s function for the non - condensate atoms @xmath109 beliaev @xcite effectively showed ( in modern matrix notation ) @xmath110 where we use the standard convention that repeated indices @xmath111 are summed over .
this is the famous dyson - beliaev equation , involving a 2 x 2 matrix self - energy @xmath112 .
clearly all components @xmath113 will share the same single - particle excitation spectrum . in a uniform system , we have @xmath114 which is a set of linear algebraic equations which are easy to solve for @xmath115 and @xmath116 .
we will now use this beliaev formalism to discuss various simple approximations for the self - energies @xmath117 . the interaction energy in ( [ eq : elementary excitations3 ] ) , namely @xmath118 splits into various distinct contributions when we use ( [ eq : green function1 ] ) to separate out the condensate parts ( see fig .
2 ) . at @xmath119
we can ignore the @xmath120 and @xmath121 contributions because so few atoms are in non - condensate .
this is the famous bogoliubov approximation @xcite , and it is equivalent to the linearized gp theory . in beliaev language ,
the bogoliubov self - energies are shown in fig .
3 . for a uniform gas ,
these give @xcite @xmath122 } \nonumber\\ { \tilde g}_{12}(p , \omega ) & = & { -n_cg\over \omega^2 -   \left[\epsilon^2_p + 2\epsilon_p n_cg\right]}\ .\label{eq : green   function11}\end{aligned}\ ] ] these bogoliubov single - particle green s functions contain the same physics as discussed in fetter s lectures @xcite .
they clearly have poles at the frequencies @xmath123 , where @xmath124 at low @xmath125 , the single - particle excitation is phonon - like @xmath126 , with the bogoliubov phonon velocity @xmath127 it turns out that this simple @xmath0 bogoliubov theory already exhibits most of the structure which will always arise in bose - condensed systems .
this is why it plays the role of the `` h - atom '' in discussions of bose - condensed fluids @xcite .    in this approximation , as we have stressed , the starting point assumes that all atoms are in the condensate and hence @xmath128 .
however , we can still use @xmath129 in ( [ eq : green function11 ] ) to estimate the depletion .
this makes use of the relation @xmath130 which is an exact formula for @xmath131 ( see ch.4 of ref .
@xcite ) in terms of the so - called single - particle _ spectral density _ [ see ( [ eq : density fluc2 ] ) ] .
the spectral density given by the bogoliubov approximation above is @xmath132 showing the characteristic negative energy pole at @xmath133 this last feature emphasizes that creating an atom in a bose - condensed system involves a coherent weighted combination of creating an excitation _ and _ destroying an excitation . calculating the depletion
, one finds @xmath134 \nonumber\\ & = & \int d{\bf q}\left\{n^0(e_q)u^2_q   -n^0(-e_q)v^2_q\right\}\nonumber\\ & = & \int d{\bf q}\left\{v^2_q + n^0(e_q)(u^2_q +   v^2_q)\right\}\nonumber\\ & = & { 8\over 3 } n_c \left(n_c a^3\over\pi\right)^{1/2 } , \label{eq : green function16}\end{aligned}\ ] ] a result first obtained by bogoliubov @xcite and reviewed in detail by fetter @xcite .
we have made use of the key identity which the bose distribution satisfies , @xmath135 .
\label{eq : green function17}\ ] ] we see from ( [ eq : green function16 ] ) that the @xmath0 depletion of the condensate formally arises from the negative energy pole in @xmath136 .
it is easy to evaluate the density response function @xmath137 in ( [ eq : density fluc1 ] ) using the bogoliubov spectral density in ( [ eq : green function15 ] ) .
the negative energy poles are seen to give rise to new poles at @xmath138 which have finite weight even at @xmath0 .
these are in addition to the `` normal '' particle - hole poles of @xmath67 at @xmath139 , as given by ( [ eq : density fluc3 ] ) .
at finite @xmath140 , a lot of atoms are thermally excited out of the condensate in a dilute bose gas .
these `` normal '' self - energies must be included . in the beliaev self - energy formalism , this was first done by popov in 1965 @xcite .
we have to add the _ ordinary _ hartree - fock self - energies , as shown in fig.4 .
this gives the first - order popov approximation ( the excited atoms are treated as as ideal bose gas , and thus the approximation is not self - consistent ) @xmath141 these self - energies lead to [ compare with ( [ eq : green function11 ] ) ] @xmath142 with @xmath143[eq : beyond3 ] where @xmath144    the `` chemical potential '' @xmath145 as defined in ( [ eq : beyond4 ] ) here was introduced by hugenholtz and pines @xcite at @xmath146 ( and generalized to @xmath147 later @xcite ) .
hp showed that if @xmath148 had a gapless excitation spectrum in the @xmath149 limit , then the true chemical potential @xmath150 must satisfy @xmath151 . in the popov approximation
, calculation gives @xcite @xmath152 and thus @xmath153 , where @xmath154 is the ideal gas condensate density at temperature @xmath140 .
the generalized gp equation determines the chemical potential , which we find convenient to denote by @xmath155 . in the popov approximation , this equation is given by @xmath156 and hence in a uniform gas , we see that @xmath157 . in the popov approximation ,
the single - particle green s functions @xmath158 in ( [ eq : beyond2 ] ) have poles given by @xmath159^{1/2}.\label{eq : beyond7}\ ] ] this is formally the same as the @xmath0 bogoliubov excitation frequency , except that now @xmath154 is temperature - dependent . for a more detailed discussion , see chapter 3 of the recent review by shi and griffin @xcite .
the _ full _ hartree - fock bogoliubov ( hfb ) approximation involves calculating @xmath117 self - consistently using the complete matrix @xmath107 beliaev propagator , as shown in fig .
we note that above @xmath3 , the hfb reduces to the usual self - consistent hartree - fock approximation . moreover , the hfb is the best _ single - particle _
approximation for a bose - condensed system , in a variational sense . within this limitation
, the hfb will give the best single - particle approximation for the thermodynamic properties .    however , there is one  bad " aspect of the full hfb which was noticed in the 1960 s and which has been the subject of a several recent theoretical papers on bose - condensed gases .
namely , the hfb does not obey the hugenholtz - pines ( hp ) theorem .
in a uniform gas , this means that the single - particle excitations will have an energy gap in the long wavelength limit _
_ @xmath160 0 .
it is easy to check this by calculating the two chemical potentials introduced earlier : @xmath161 the static equation of motion for @xmath162 using the hfb approximation gives @xmath163{\sqrt n_{c}}= \mu_{c}\sqrt   n_{c}\rightarrow\mu_{c } =   g(n+\tilde{n}+\tilde{m}).\label{eq : beyond9}\ ] ] one `` solution '' of this problem is to leave out @xmath164 ( which is the source of the problem ) but keep @xmath165 .
this corresponds to the self - consistent hfp(opov ) approximation @xcite which , while approximate , has several nice features .
it gives the correct excitation spectrum for both @xmath166 0 and @xmath84 .
moreover , the spectrum is gapless at all temperatures .
the origin of this  problem " with the hfb is clear .
the hfb keeps all self - energies which are first order in the interaction @xmath56 .
however , by computing these self - energy diagrams using self - consistent @xmath113 propagators , one is clearly bringing in terms to all orders in @xmath56 .
this brings one into dangerous territory ! moreover , one easily can check that @xmath164 must be at least of order @xmath56 .
thus the @xmath167 contribution to @xmath168 is at least of @xmath169 ) .
this suggests that to `` fix - up '' the hfb as a theory of excitations , we have to include all self - energy contributions to at least second order in @xmath56 .
moreover , it is not obvious that it is consistent to calculate @xmath170 self - consistently @xcite and ignore @xmath164 , since the lowest order interaction contributions to both quantities are found to be of the same order .
indeed , one can show for a uniform gas that keeping the lowest order asymptolic correction to the ideal gas result gives @xmath171 , where @xmath172 is the critical density of an ideal gas at temperature @xmath140 @xcite .
however , we note that an expansion around the ideal gas results is not necessarily a good approximation in trapped gases ( see also the discussion of the results in ( [ eq : two - fluid7 ] ) . at @xmath146 ,
the popov approximation is a poor guide for corrections to the bogoliubov approximation since a calculation analogous to ( [ eq : green function16 ] ) gives @xmath173 @xcite .
so , finally , we come to the work of beliaev @xcite who , in 1957 , evaluated the @xmath117 self - energies keeping all contributions up to order @xmath174 .
this was originally done at @xmath146 . for finite temperatures @xmath175 ,
beliaev s work has been recently extended by several authors @xcite .
when one includes the second - order diagrams for @xmath117 , we note that one must be careful to also treat the hfb first - order diagrams correctly to order @xmath174 since there are many terms which cancel .
as expected , the beliaev - type second order calculation cures all the problems of the first - order hfb approximation .
we make a few comments on the results for a uniform bose gas .
the self - energies are found to be    @xmath176    with the renormalized interaction defined by @xmath177\label{eq : beyond11}\ ] ] and @xmath178.\label{eq : beyond12}\ ] ] here the anomalous and normal densities , @xmath179 and @xmath180 , are computed using the hfp approximation for @xmath181 .
@xmath182 is the second - order approximation to the @xmath183-matrix , as shown in fig .
6 . using these results to calculate the chemical potential , we find ( for @xmath184 @xmath185 this last result gives precisely the hfb result for @xmath186 , and hence the hp relation is satisfied . in the last line of ( [ eq : beyond13 ] ) , we have used the key result that for @xmath187 , @xmath188 , as proven at all temperatures by talbot and griffin @xcite .
we note that one can rewrite the expression in ( [ eq : beyond13 ] ) in the alternative form @xmath189,\label{eq : beyond14}\ ] ] where we have introduced a `` renormalized '' anomalous density @xmath190.\label{eq : beyond15}\end{aligned}\ ] ] in contrast to @xmath191 in ( [ eq : beyond12 ] ) , @xmath192 has no ultraviolet divergence from large @xmath193 contributions .
proukakis , burnett and coworkers @xcite have pointed out that the self - consistent ladder diagram approximation for the @xmath194matrix ( see ref.@xcite and ch.4 of ref.@xcite ) can be expressed in terms of @xmath195 . in a uniform gas ,
this relation is simply @xmath196 .
this suggests that the effect of @xmath197 can be partially included by working in terms of the self - consistent @xmath194matrix @xcite .
we have seen in section  [ sec : beyond the bogoliubov approximation ] that beliaev s second - order calculation satisfies the hp theorem .
it has another feature whose significance was only realized in the early 1960 s . using the above @xmath0 beliaev expression ( [ eq : beyond14 ] ) for the chemical potential @xmath150 , one finds at @xmath0 ( a result first found by lee and yang using a different approach ) @xmath198 ,   \label{eq : approximation1}\ ] ] where we have used [ see ( [ eq : green function16 ] ) ] @xmath199
one also finds that the phonon pole of @xmath200 has a velocity given by @xmath201 ,   \label{eq : approximation3}\ ] ] again expressed in terms of the total density @xmath202 .
one may easily check that @xmath203 that is , the beliaev phonon excitation at @xmath0 has a velocity precisely equal to the _ compressional _ sound velocity .
this last result raises the basic question : why does the elementary excitation ( pole of @xmath204 have a velocity corresponding to that of a density fluctuation ( pole of @xmath205 ?
later this strange identity was proven to be correct at @xmath0 to _ all _ orders of perturbation theory by gavoret and nozi@xmath206res @xcite . to be precise , at @xmath207 in any bose - condensed fluid ( liquid or gas ) , gn showed ( for @xmath184 @xmath208 where @xmath209 is given by the compressibility as in ( [ eq : approximation4 ] ) .
we discuss this equivalence in section  [ sec : why are excitations in a bose fluid so interesting ] .
let us consider the full time - dependent hfb equation of motion for @xmath210 , given by @xmath211 \phi({\bf r},t ) .
\label{eq : approximation6}\ ] ] there are various approximations to this key equation which have been used in recent discussions of trapped bose gases :    * _ gp : _ ignore @xmath165 and @xmath164 completely .
this is valid at @xmath146 @xcite . *
_ static hfp : _ keep @xmath212 but set @xmath213 and @xmath214 @xcite . * _ static hfb : _
keep @xmath215 but set @xmath216 and @xmath217 @xcite . as we have noted , this produces an energy gap in the excitation spectrum ( since @xmath218 . * _ dynamic hfb : _ treat all dynamic mean - fields due to @xmath219 and @xmath220 on an equal basis in a generalized mean - field calculation of the density response function @xmath221 .
one finds @xcite that the poles of @xmath72 given by this kind of calculation @xmath222 are _ identical _ to the poles of @xmath181 given by the second - order beliaev approximation @xmath223 ) .
recent work on excitations at @xmath2240 in trapped gases is easily understood in terms of such a generalized linear response calculation of @xmath72 .
for example , minguzzi and tosi @xcite only keep fluctuations in @xmath225 and @xmath226 but @xmath227 is completely ignored ( dynamic hfpopov ) ; in contrast , giorgini @xcite keeps fluctuations in both @xmath226 and @xmath227 but only those induced by the condensate fluctuations @xmath225 .
let us briefly discuss the excitation frequencies in trapped gases at finite temperatures  @xcite . for @xmath228 @xmath3 ,
gp theory @xmath229 predictions @xcite are in excellent agreement with the jila @xcite and mit @xcite experiments . for @xmath230 @xmath3 , however , there are significant differences between most current theoretical results and the measured excitation energies  @xcite .
the problem clearly is that mean - field theories must include the collective dynamics of _ both _ the condensate and the non - condensate .
the beliaev treatment does this , as do the papers @xcite mentioned above , as well as work based on a collisionless kinetic equation @xcite .
these improved treatments can correctly obtain the in - phase rigid oscillation of the condensate and non - condensate equilibrium profiles ( the so - called kohn mode ) , as discussed in section  [ sec : in - phase oscillation of the condensate and non - condensate ] .
we have seen in section  [ sec : physics of the beliaev approximation ] that at @xmath146 , the elementary excitations have the same spectrum as density fluctuations .
in fact , one can prove this equivalence to be a direct consequence of the bose broken symmetry , @xmath231 .
thus , this equivalence is valid at all @xmath232 . of course , as soon as @xmath233 , the phonon velocity is no longer simply given by the compressional sound velocity ( [ eq : approximation4 ] ) as found at @xmath146 .
this fundamental equivalence was only proven and understood in the early 1960 s ( gavoret and nozi@xmath206res @xcite , hohenberg and martin @xcite , bogoliubov @xcite ) .
it was implicitly assumed by landau @xcite in his 1941 paper and later by feynman @xcite in his 195354 papers on excitations .
of course , the clear distinction between elementary excitations and density fluctuations in many - body systems was only clarified in the period starting from 1957 . only then were bose fluids realized to be very special @xcite .
today , we understand this key hidden assumption of the work of landau and of feynman has its microscopic basis in @xmath234 .
let us try to understand the origin of this equivalence of the field fluctuation and density fluctuation spectrum .
the density operator in a bose - condensed system can be decomposed using  ( [ eq : green function1 ] ) as follows : @xmath235 where the first three terms correspond to @xmath236 .
clearly density fluctuations in a bose - condensed system ( ie , @xmath237 ) have a direct coupling to the single - particle field fluctuations , due to the possibility of atoms coming in or out of condensate reservoir . in momentum space
, we recall that ( [ eq : excitations1 ] ) is equivalent to @xmath238 making the same decomposition as in ( [ eq : excitations1 ] ) for the density response function ( [ eq : elementary excitations5 ] ) gives @xmath239 we note that in the bogoliubov approximation , we have @xmath240 and the equivalence of @xmath241 and @xmath181 spectrum is then trivial ( see also @xcite ) @xmath242    to illustrate how the hybridization of single - particle and density fluctuation arises , we discuss the time - dependent hfp approximation on which ref .
@xcite is based . for simplicity ,
we consider the uniform gas case . in the mfa of section  [ sec : density fluctuation spectrum in the mean - field approximation ] , we have @xmath243 , \label{eq : excitations5}\ ] ] and hence for @xmath244 , @xmath245 using this result in the hfp equation of motion  ( [ eq : approximation6 ] ) , ie , with @xmath246 , we find @xmath247 in this approximation , which only keeps the @xmath248 fluctuations induced by the condensate , the equation of motion for @xmath249 reduces to the gp equation , but with @xmath250 .
thus the characteristic poles will be determined by the zeroes of @xmath251 \propto [ \omega^2 - ( \epsilon^2_p + 2 g n_c\epsilon_p ) ] [ 1 - 2g\chi^0 ] - 4g^2 2n_c \epsilon_p \chi^0 .
\label{eq : excitations8}\ ] ] this is the same denominator ( up to a common factor ) exhibited by the response functions obtained by minguzzi and tosi  @xcite .
this calculation illustrates how the bogoliubov single - particle mode and the zero sound density fluctuation are coupled and hybridized .
more generally , one can prove that ( at all @xmath252 : @xmath253 where @xmath254 is a bose broken - symmetry vertex function which vanishes if @xmath255 . using the  dielectric formalism " developed in the early 1970
s ( for a review , see ch.5 of ref .
@xcite ) , one can prove that the self - energies @xmath117 are such that @xmath72 and @xmath181 have the same poles  and that there are none specifically associated with @xmath256 .
the advantage of this diagrammatic formalism @xcite    is that one manifestly sees that the poles of @xmath181 and @xmath72 are identical , within a given approximation for the  building blocks " of the dielectric formalism .
we now see why we can study elementary excitations from the density fluctuation spectrum when @xmath237 . to directly measure @xmath257 , we need an external perturbation which can change the number of particles in the system .
bec indirectly provides us with such a probe !
this equivalence of the single - particle excitations with density fluctuations lies at the heart of the phenomenon of superfluidity .
it essentially _ restricts _ the possible excited states to density fluctuations and thus ensures the stability of superfluid motion .
however , as nozires says in a very lucid review article  @xcite , `` that the @xmath146 superfluid equations merge with those of ordinary hydrodynamics does not alter the fact that real understanding must be based on a microscopic description based on long - range order induced by a condensate '' .
in the previous sections , we have been discussing collective oscillations of the condensate and non - condensate in the mean - field or _ collisionless _ region .
we now turn to the collision - dominated _ hydrodynamic _ region @xmath258 , where @xmath259 is the wavelength of the collective mode and @xmath260 is the collisonal mean - free - path of the elementary excitations @xcite .
experiments in trapped bose gases are just now starting to probe the hydrodynamic two - fluid region but this region is very rich in new physics and deserves careful study .
it is unfortunate that the collisionless region at @xmath0 described by the gross - pitaevskii equation is also commonly referred to as `` hydrodynamic '' theory @xcite .
i guess we are stuck with this terminology , although the terms `` quantum hydrodynamics '' or `` superfluid hydrodynamics '' would be useful compromises .
we will limit our analysis to _ finite _ temperatures @xmath261 , where the number of atoms in the non - condensate is comparable to the condensate @xmath262 .
we will concentrate on the _ derivation _ of two - fluid hydrodynamic equations starting from a simple microscopic model and try to clarify the physics .
these two - fluid hydrodynamic equations are expressed in terms of fluctuations of the two components : @xmath263 and @xmath264 .
we can use them to `` derive '' the two - fluid equations in the conventional landau form , which are written in terms of the fluctuations of local thermodynamic variables ( temperature , entropy , pressure , etc ) .
we mainly consider the linearized version of the zgn - type two - fluid equations , which give the hydrodynamic normal modes of the coupled system of two components . as a concrete example , we discuss the _ in - phase _ rigid motion of the equilibrium density profiles of the condensate and non - condensate .    we need to find equations of motion for both the condensate and the non - condensate . as a first approximation @xcite
, we can use the time - dependent hf popov equation of motion ( see section  [ sec : physics of the beliaev approximation ] ) , for the order parameter @xmath265 @xmath266\phi({\bf r } , t ) , \label{eq : hf popov equation motion1}\ ] ] where @xmath267 and @xmath268 is the local non - condensate density .
it is clear that this equation of motion for @xmath269 is not closed since it involves @xmath270 if we treat the non - condensate in ( [ eq : hf popov equation motion1 ] ) _ statically _ , this means that the condensate moves in the _
static _ hartree - fock mean field of non - condensate given by @xmath271 @xcite more generally , we need to have an equation of motion for the fluctuations in @xmath272 of the non - condensate ( ie , the density of the excited atoms ) , which will be discussed in section  [ sec : dynamics of the non - condensate atoms ] .
we can rewrite equation of motion in ( [ eq : hf popov equation motion1 ] ) in the _ quantum hydrodynamic _
variables for phase and amplitude , defined by @xmath273 where the superfluid velocity is given by @xmath274 using ( [ eq : hf popov equation motion3 ] ) , ( [ eq : hf popov equation motion1 ] ) can be shown to be equivalent to ( all quantities depend on @xmath275 and @xmath183 ) @xmath276 & = &   -\mbox{\boldmath$\nabla$}\mu_c,\label{eq : hf popov equation motion5}\end{aligned}\ ] ] where the _ superfluid chemical potential _ is given by ( in the dynamic hfp approximation ) @xmath277 the first term in ( [ eq : hf popov equation motion6 ] ) is the `` quantum pressure '' term , which is _ ignored _ in the so - called thomas - fermi ( tf ) approximation to the condensate equation of motion @xcite . if we omit the non - condensate term , equations ( [ eq : hf popov equation motion5 ] ) and ( [ eq : hf popov equation motion6 ] ) describe the well - known gross - pitaevskii approximation @xcite    .
following ref.@xcite , we base our discussion on the simplest possible _ kinetic equation _ for the distribution function @xmath278 valid in the semi - classical limit : @xmath279 this distribution function is determined by the quantum boltzmann equation @xmath280 f({\bf r , p } , t ) = \left.{\partial f({\bf r , p } , t)\over\partial   t}\right|_{\mbox{collision } } , \label{eq : semi - classical limit2}\ ] ] where the effective potential @xmath281.\label{eq : semi - classical limit3}\ ] ] includes the dynamic hf field .
once we solve ( [ eq : semi - classical limit2 ] ) for @xmath282 we can find the _ non - condensate _ density from @xmath283 in our simple model , the excited atoms are still particle - like , ie , @xmath284 this last fact leads to many simplifications but restricts the results to finite temperatures close to @xmath3 .
we now give the explicit form for the collision term in ( [ eq : semi - classical limit2 ] ) : @xmath285 + c_{12}[f ] .
\label{eq : semi - classical limit6}\ ] ] the term @xmath286 $ ] describes collisions between excited atoms and is given by @xcite , @xmath287 & = & 2g^2 \int{d{\bf p}_1\over(2\pi)^3}{d{\bf   p}_2\over(2\pi)^3}{d{\bf p}_3\over(2\pi)^3}(2\pi)^3\delta({\bf   p}+{\bf p}_1 - { \bf p}_2 - { \bf p}_3)\nonumber\\ & \times & 2\pi\delta(e+e_1-e_2-e_3 )   \left[(1+f)(1+f_1)f_2f_3-ff_1(1+f_2)(1+f_3)\right ] .
\label{eq : semi - classical limit7}\end{aligned}\ ] ] in contrast , @xmath288 $ ] describes collisions involving one condensate atom .
this kind of collision brings about equilibrium between the _ excited _ atoms and _ condensate _ atoms and it is given by ( these collision terms were ignored in ref.@xcite ) @xmath289&= & 2g^2 \int{d{\bf p}_1\over(2\pi)^3}{d{\bf   p}_2\over(2\pi)^3}{d{\bf p}_3\over(2\pi)^3}(2\pi)^3\delta({\bf p}_1 +   { \bf p}_s - { \bf p}_2 - { \bf p}_3)\nonumber\\ & & \qquad\hbox{}\times   2\pi\delta(e+\epsilon_c - e_2-e_3)(2\pi)^3[\delta({\bf p}-{\bf p}_1 )   -\delta({\bf p}-{\bf p}_2 ) - \delta({\bf p}-{\bf p}_3)]\nonumber\\ & & \qquad\hbox{}\times   \left[n_c(1+f_1)f_2f_3-n_cf_1(1+f_2)(1+f_3)\right ] , \label{eq : semi - classical limit8}\end{aligned}\ ] ] with @xmath290 and @xmath291 .
these two collision terms were first derived by kirkpatrick and dorfman @xcite , who considered a uniform gas and worked in a frame in which @xmath292 .
our present calculation @xcite is for trapped bose gas and works in the lab frame . as a result , ( [ eq : semi - classical limit8 ] ) takes into account that a condensate atom has energy @xmath293 and momentum @xmath294 , where @xmath155 is defined in ( [ eq : hf popov equation motion6 ] ) .
we note that :    * @xmath286 $ ] conserves number , momentum and energy of the colliding excited atoms and therefore @xmath295 = 0 \nonumber\\ & \int & d{\bf p}{\bf p } c_{22}[f ] = 0 \nonumber\\ & \int & d{\bf p } e_p c_{22 } [ f ] = 0 \label{eq : semi - classical   limit9}\end{aligned}\ ] ] * @xmath288 $ ] conserves momentum and energy of the colliding atoms , with the condensate atoms having energy @xmath296 and momentum @xmath297 .
therefore , we have    @xmath298 = 0 \nonumber\\ & \int & d{\bf p } ( e_p - \epsilon_c ) c_{12 } [ f ] = 0 ,   \label{eq : semi - classical limit10}\end{aligned}\ ] ]    but @xmath299 \equiv \gamma_{12}[f ] \neq 0 .
\label{eq : semi - classical limit11}\ ] ] the fact that @xmath300 is finite follows since @xmath301 does _ not _ conserve the number of _ excited _ atoms ; @xmath300 will be referred to as a source term .    as usual in dealing with the collision - dominated hydrodynamic region
, we will assume that the collisions @xmath286 $ ] produce _ local thermal equilibrium _ among the excited atoms .
this is described by the _
bose distribution : @xmath302}-1 } , \label{eq : semi - classical limit12}\ ] ] where now @xmath303 and @xmath304 all depend on @xmath305 .
one can easily checks that @xmath306 $ ] vanishes exactly .
this uses the key identity for bose distribution given by ( [ eq : green function17 ] ) .
in contrast , one finds that substituting @xmath307 in ( [ eq : semi - classical limit8 ] ) gives @xmath308 & = & 2g^2 \int{d{\bf p}_1\over(2\pi)^3}{d{\bf   p}_2\over(2\pi)^3}{d{\bf p}_3\over(2\pi)^3}(2\pi)^3\delta({\bf   p}_1+{\bf p}_s - { \bf p}_2 - { \bf p}_3)\nonumber\\ & \times & 2\pi\delta(e_1+\epsilon_c - e_2-e_3)(2\pi)^3[\delta({\bf   p}-{\bf p}_1 ) -\delta({\bf
p}-{\bf p}_2 ) - \delta({\bf p}-{\bf   p}_3)]\nonumber\\ & \times &   \left[n_c - n_c e^{-\beta({\tilde\mu } ' - \mu_c)}\right](1 +   f_{01 } ) f_{02 } f_{03 } , \label{eq : semi - classical limit13}\end{aligned}\ ] ] where @xmath309 in a _
linearized _ theory .
we note that @xmath310 $ ] in ( [ eq : semi - classical limit13 ] ) vanishes when @xmath311 ie , only if the condensate atoms are in _ diffusive equilibrium _ with the non - condensate atoms .
the same kind of factor arises in the theory of the growth dynamics of a condensate in trapped bose gases developed by gardiner , zoller and coworkers @xcite .
we next derive _ hydrodynamic equations _ for the non - condensate atoms by taking moments of the kinetic equation : @xmath312 f({\bf r , p } , t ) = c_{22}[f ]   + c_{12}[f ] .
\label{eq : semi - classical limit15}\ ] ] we obtain @xcite @xmath313 & \rightarrow & { \partial{\tilde   n}\over\partial t}+\mbox{\boldmath$\nabla$}\cdot{\tilde n}{\bf v}_n =   \gamma_{12}\nonumber\\ \int d{\bf pp } [ \mbox{kinetic eq . } ] & \rightarrow & m{\tilde   n}\left({\partial\over\partial t}+{\bf   v}_n\cdot\mbox{\boldmath$\nabla$}\right){\bf v}_n =   -\mbox{\boldmath$\nabla$}{\tilde p}-{\tilde   n}\mbox{\boldmath$\nabla$}u - m({\bf v}_n - { \bf   v}_s)\gamma_{12}\nonumber\\ \int d{\bf p}p^2 [ \mbox{kinetic eq . } ] & \rightarrow & { \partial { \tilde   p}\over\partial t}+\mbox{\boldmath$\nabla$}\cdot({\tilde p}{\bf v}_n )   =   -{2\over 3}{\tilde p}\mbox{\boldmath$\nabla$}\cdot{\bf v}_n\nonumber\\ & \quad & \quad + { 2\over 3}[\mu_c + { 1\over 2}m({\bf
v}_n - { \bf   v}_s)^2 - u]\gamma_{12}. \label{eq : semi - classical limit16}\end{aligned}\ ] ] here @xmath300 is defined in ( [ eq : semi - classical limit11 ] ) and @xmath314 with the local fugacity @xmath315 .
in sections  [ sec : hydrodynamic oscillations in a trapped bose gas ] and [ sec : dynamics of the non - condensate atoms ] , we derived a set of 7 equations for non - condensate and the 4 equations for the condensate .
this gives a coupled set of equations for the variables : @xmath316.\label{eq : two - fluid1}\ ] ] > from ( [ eq : semi - classical limit13 ] ) , one sees @xmath300 is function of @xmath317 , which in turn depends on @xmath318 and @xmath304 .
thus our system of hydrodynamic equations is closed and can be solved since it involves 11 local variables and 11 equations .    for reasons discussed earlier , the explicit effect of the @xmath319 collisions has disappeared in the hydrodynamic equations ( [ eq : semi - classical limit16 ] ) .
however , the @xmath319 collisions justify and enforce the local equilibrium form ( [ eq : semi - classical limit12 ] ) used for @xmath320 which describes the excited atoms .
in contrast , @xmath301 is still explicitly present in ( [ eq : semi - classical limit16 ] ) through the @xmath300 source term .
this means that the condensate is not in _ diffusive _ equilibrium with non - condensate .
we recall from ( [ eq : semi - classical limit13 ] ) that @xmath321 ,   \label{eq : two - fluid2}\ ] ] and thus @xmath322 $ ] vanishes only when @xmath323 .
one might expect that this condition would arise in the case of strong @xmath301 collisions .
a detailed study @xcite shows that even in this case , one must allow for fluctuations in the variable @xmath324 and thus we can not simply set @xmath325 since @xmath326 in ( [ eq : semi - classical limit16 ] ) , we have to generalize the hf popov description to ensure that [ compare with the result in ( [ eq : hf popov equation motion5 ] ) ] @xmath327 this is needed for the continuity equation to be correct @xmath328 where @xmath329 as discussed in detail in ref.@xcite , a simple microscopic model based on the neglect of anomalous pair correlation functions @xmath330 leads to both ( [ eq : two - fluid3 ] ) and ( [ eq : semi - classical limit16 ] ) .
we can linearize our coupled two - fluid hydrodynamic equations given by ( [ eq : hf popov equation motion5 ] ) and ( [ eq : semi - classical limit16 ] ) around the static equilibrium solution , where @xmath331 .
the fluctuations @xmath332 are precisely the quantities which can be directly measured in trapped gases .
the only static equilibrium functions we need are @xmath333 and @xmath334 , which are given by ( [ eq : semi - classical limit17 ] ) , @xmath335 where @xmath336 and @xmath337 is given by the _
static _ hfp equation ( [ eq : beyond6 ] ) for @xmath338 .
we note that these static thermodynamic properties for a bose - condensed gas are quite different from the expressions first given by lee and yang ( 1958 ) for a uniform gas @xcite .
they calculated only the first - order interaction corrections to the thermodynamic properties of an ideal bose gas , using @xmath339 here @xmath340 is the critical density at the temperature @xmath140 for an ideal uniform bose gas .
the non - perturbative approximation we use for thermodynamic properties [ as illustrated by ( [ eq : two - fluid6 ] ) ] is consistent with the lee - yang results when we expand to first order in @xmath56 @xcite .
we can use our generalized zgn equations to derive two - fluid equations in the traditional form first given by landau ( 1941 ) . in linearized form , these are @xcite @xmath341 where we have included an external potential , @xmath342 is the local entropy density and @xmath343 of course , the landau equations describe the case when the superfluid and normal fluid are in dynamic local equilibrium with each other .
the landau form of two - fluid equations is very natural for superfluid @xmath5he , for which they were developed .
they involve the usual thermodynamic variables : pressure , entropy and temperature fluctuations .
however , in a trapped gas , the more natural variables are those used in a zgn - type formulation @xcite , namely @xmath344 and @xmath345 .    in the generalized zgn equations discussed in ref .
@xcite , the fluctuations around static equilibrium involve fluctuations in @xmath346 , namely [ see ( [ eq : two - fluid2 ] ) ] @xmath347 this leads to an additional equation of motion for @xmath348 , with a characteristic relaxation time @xmath349 .
this gives rise to a new relaxational mode of the condensate .
the zgn limit @xcite corresponds to @xmath350 while the landau two - fluid limit corresponds to @xmath351 while the relaxational mode is not included , the zgn equations appear to give a reasonable first approximation for the other hydrodynamic modes , as shown in the explicit results given in ref.@xcite .
we note that landau two - fluid equations given by ( [ eq : two - fluid8 ] ) have been recently used by shenoy and ho @xcite to work out the hydrodynamic normal modes of a trapped superfluid gas .
the landau two - fluid equations , which assume total local equilibrium , are valid at all temperatures .
our present analysis is restricted to finite @xmath352 , but could be generalized to the case of very low temperatures . at @xmath0 , of course , the distinction between the collisionless and hydrodynamic region is not a useful one ( see , moreover , the comment at the end of section  [ sec : why are excitations in a bose fluid so interesting ] ) . however , at @xmath0 , we have @xmath353 and @xmath354 and then the landau two - fluid equations formally reduce to hydrodynamic equations for _ one _ component : @xmath355 these equations are the basis of recent work by pitaevskii and stringari @xcite . in conjuction with the local density approximation @xcite @xmath356 + u_{\rm   ext}({\bf r } )
,   \label{eq : two - fluid12}\ ] ] where @xmath357 is the chemical potential of a _ uniform _ interacting bose - condensed gas at @xmath0 , such as given by ( [ eq : approximation1 ] ) , one may use ( [ eq : two - fluid11 ] ) to find collective mode frequencies past the mfa ( or gp approxiamation ) .
two interesting collective modes exhibited by the zgn equations involve _ rigid _ centre - of - mass oscillations of the equilibrium density profiles : @xmath358 one finds @xcite an in - phase solution ( kohn mode ) with @xmath359 , with a frequency given by the trap frequency @xmath360 ; and an out - of - phase solution , with @xmath361 and @xmath362 in opposite directions .
the latter is the _ analogue _ of the second sound mode in liquid @xmath5he .
first let us work out the static equilibrium solution @xmath363 of ( [ eq : semi - classical limit15 ] ) , which must satisfy @xmath364 = 0 .
\label{eq : oscillation of the cond2}\ ] ] the solution is easily verified to be @xmath365 where ( within the thomas - fermi approximation for @xmath366 ) @xmath367 .
\label{eq : oscillation of the cond4}\end{aligned}\ ] ] using ( [ eq : oscillation of the cond3 ] ) , we find @xmath368 is given by ( [ eq : two - fluid6 ] ) , with @xmath369    we next consider the in - phase solution described by ( [ eq : oscillation of the cond1 ] ) , with @xmath370 ( not dependent on the position @xmath275 ) .
this is a solution of zgn equations with @xmath371 , where @xmath360 is the parabolic trap frequency in @xmath372 direction @xmath373 .
the proof proceeds as follows ( in the linearized theory ) .
given ( [ eq : oscillation of the cond1 ] ) , we have @xmath374 and hence @xmath375   \nonumber \\ & = & \mbox{\boldmath$\nabla$ } \left({\mbox{\boldmath$\eta$ } } \cdot [ 2 g   \mbox{\boldmath$\nabla$}{\tilde n}_0 + g   \mbox{\boldmath$\nabla$}n_{c0 } ] \right).\label{eq : oscillation of the   cond7}\end{aligned}\ ] ] this result is equivalent to [ using ( [ eq : oscillation of the cond4 ] ) ] @xmath376 with @xmath377 the solution of ( [ eq : oscillation of the cond8 ] ) is clearly @xmath378    in a similar way , one can show that the hydrodynamic equations for the non - condensate lead to the same equation of motion ( [ eq : oscillation of the cond8 ] ) .
the rigid motion of the non - condensate density profile corresponds to a distribution function given by @xmath379 this can be verified to solve the kinetic equation ( [ eq : semi - classical limit15 ] ) , using the fact that @xmath380 = 0 $ ] , and leads to the expression in ( [ eq : oscillation of the cond1 ] ) , namely @xmath381 this in - phase dipole mode was first exhibited by zgn in the hydrodynamic region @xcite .
it is an important test of the proper treatment of the non - condensate dynamics .
this kind of mode is generic , existing in the collisionless region as well @xcite .
the simplicity of the zgn two - fluid equations allow a very transparent discussion of first and second sound in a uniform gas @xcite . we find two sound - wave solutions : @xmath382 , @xmath383 .
first sound has a velocity given by ( to lowest order in @xmath56 ) @xmath384 and mainly involves the non - condensate , with @xmath385 ( but in - phase ) .
second sound has a velocity given by ( to lowest order in @xmath56 ) @xmath386 and mainly involves the oscillation of the condensate , with @xmath387 ( and out - of - phase ) .
it is clear that second sound is the collision - dominated hydrodynamic mode which is the _ analogue _ of the collisionless bogoliubov phonon ( [ eq : green function13 ] ) which we discussed in section  [ sec : green s function formulation ] .
the frequency vanishes ( becomes `` soft '' ) above @xmath3 and it is the expected goldstone mode of the bose broken - symmetry @xcite .
in contrast to liquid @xmath5he , both second and first sound in a bose gas involve density fluctuations ( ie , second sound is not a temperature wave in a gas ! ) .
can we study first and second sound pulse propagation in cigar - shaped traps ? in this volume , ketterlee @xcite reviews the beautiful pulse propagation experiments carried out at mit at low temperatures in the collisionless region @xcite . if we could get into the hydrodynamic region , one should be able to see _ two _ pulses propagate
roughly with the speeds @xmath388 and @xmath389 given above .
this would be very dramatic and a direct confirmation of the superfluid dynamics described by the two - fluid equations .
inside the broad , relatively uniform condensate distribution along the axis of the cigar - shaped trap , one could excite both first and second sound pulses , with the relative intensity of each propagating pulse given by @xcite @xmath390 the velocity of second sound would give [ using ( [ eq : sound velocity2 ] ) ] a direct measurement of the condensate density @xmath391 .
perhaps of more interest , the velocity of the first sound pulse would give a relatively _ direct _ measurement [ using ( [ eq : sound velocity1 ] ) ] of the non - condensate density @xmath392 `` underneath '' the condensate .
in sections [ sec : dynamics of the non - condensate atoms ] and [ sec : the two - fluid hydrodynamic equations ] , we discussed the linearized form of the generalized two - fluid equations @xcite which describe how the condensate relaxes to equilibrium with the excited atoms
. these hydrodynamic equations with @xmath393 ( i.e. @xmath394 ) should also be useful in the study of how a condensate grows and comes into equilibrium with the thermal cloud @xcite . while based on a different formalism , our theory has points of contact with the work by gardiner , zoller and co - workers on the kinetics of condensate growth @xcite .
however , our analysis , while assuming the existence of a condensate , takes the dynamical evolution of the non - condensate fully into account .
the two - fluid dynamics of trapped bose - condensed gases promises to be a very rich subject , as discussed in ref .
the region where @xmath393 is not easily studied in superfluid @xmath5he , largely because local thermal equilibrium is almost too easy to reach in a liquid .
however , some years ago pitaevskii @xcite treated the regime very close to @xmath395 , when the superfluid density @xmath396 was not necessarily in equilibrium with the normal fluid density @xmath397 .
one of the interesting aspects of trapped bose gases is that even if the centre of the trap is at large enough density so that one is in the hydrodynamic region described by dynamic local thermal equilibrium , this description will always break down in the low density non - condensate thermal cloud far enough from the trap centre .
recent theoretical studies @xcite have shown that the cross - over region between the hydrodynamic and collisionless regions in this low density region plays the dominant role as the source of damping of certain kinds of hydrodynamic oscillations .    more generally , i would like to emphasize that the study of collective modes in trapped bose gases finally provides us with the opportunity of making a quantitative test of the complex dynamics of a system with a bose - broken symmetry . as i discuss elsewhere in this volume @xcite , the complexity of dealing with a liquid such as superfluid @xmath5he never allowed one to make a thorough confrontation of the microscopic theory of bose - condensed systems with experiment @xcite .
bose gases finally allow us to do this , in a much cleaner fashion @xcite .
i would hope that future studies on trapped bose gases will also stimulate new interest in the dynamics of superfluid @xmath5he as a bose - condensed liquid .
the work presented here on the two - fluid hydrodynamic equations has been done in close collaboration with eugene zaremba and tetsuro nikuni .
in addition , i have benefitted from stimulating discussions with many of my colleagues , and especially would like to thank k.burnett , m.l .
chiofalo , e. cornell , w. ketterle , l.p .
pitaevskii , g.v .
shlyapnikov , d.m .
stamper - kurn , s.stringari and p. zoller .                                                                                1
.   hartree - fock self - energies in a normal bose gas .
various interaction terms involving condensate ( wiggly line ) and non - condensate ( solid line ) atoms .
self - energy diagrams in the bogoliubov approximation @xmath229 .
self - energy diagrams in the first - order popov approximation .
the propagators are for an ideal bose gas . 5 .
self - consistent hartree - fock - bogoliubov ( hfb ) approximation for self - energies .
ladder diagram approximation for the @xmath183-matrix .