recent developments in the physics of ultra - cold gases have led to the discovery of bose - einstein condensation ( bec ) in trapped clouds of alkali atoms @xcite and stimulated a tremendous boost in theoretical studies of weakly interacting bose gases . as in previous years , these studies rely on the binary approximation for the interparticle interaction .
the latter is characterized by the 2-body scattering length @xmath0 , which assumes the presence of a small gaseous parameter @xmath1 ( n is the gas density ) . especially intensive
are the attempts to reach beyond the ordinary mean field approach and to develop a theory which can properly describe the behavior of finite temperature elementary excitations of a trapped bose - condensed gas and in particular , explain the jila @xcite and mit @xcite experiments on energy shifts and damping rates of the excitations .    the commonly used mean field theory ( for @xmath2 )
is based on the bogolyubov quasiparticle approach developed originally for a spatially homogeneous bose - condensed gas at @xmath3 @xcite and employed by lee and yang @xcite ( see also @xcite ) at finite temperatures .
the generalization of the bogolyubov method for spatially inhomogeneous systems has been described by de gennes @xcite . in the case of a bose - condensed gas
it should be completed by the equation for the wavefunction of the spatially inhomogeneous condensate , derived by pitaevskii @xcite and gross @xcite .    for spatially homogeneous gases the theory beyond the mean field approach was also developed .
beliaev @xcite constructed the zero - temperature diagram technique which allows one to find corrections to the energies of bogolyubov excitations , proportional to @xmath4 , where @xmath5 is the condensate density .
the corrections are provided by the interaction between the excitations ( in particular , through the condensate ) and contain both real ( energy shift ) and imaginary ( damping rate ) parts . at @xmath6
the latter originates from spontaneous decay of a given excitation ( @xmath7 ) to two other excitations ( @xmath8 and @xmath9 ) , with smaller energies and momenta : @xmath10 .
a universal expression for the chemical potential in terms of the self - energy functions has been found by pines and hugenholtz @xcite .
it should be emphasized that the corrections proportional to @xmath11 already depend on the contribution of 3-body interactions and , hence , can not be obtained within the binary approximation .
the beliaev approach was employed by popov @xcite at finite temperatures . in this case
the corrections to the beliaev self - energies contain infra - red singularities , i.e. they tend to infinity for momenta @xmath12 .
this prompted popov to make a renormalization of the theory , which links the microscopic approach with phenomenological landau hydrodynamics @xcite .
the popov theory eliminates the infra - red singularities and allows one to describe the behavior of low - energy excitations ( phonons ) at temperatures much smaller than the mean field interparticle interaction @xmath13 ( @xmath14 , with @xmath15 being the atom mass ) .
the damping of phonons in this temperature range is determined by the beliaev damping processes and has also been calculated by hohenberg and martin @xcite .
a simplified approach within the dielectric formalism was used by szepfalusy and kondor @xcite for calculating the damping rates of excitations in the phonon branch of the spectrum .
they found that at temperatures @xmath16 the damping rate of a given excitation ( @xmath7 ) originates from the scattering of thermal excitations ( @xmath8 and @xmath9 ) on the excitation @xmath7 through the processes @xmath17 .
since the characteristic energies of the thermal excitations @xmath8 , @xmath9 turn out to be much larger than the energy of the excitation @xmath7 , this damping channel can be treated as landau damping .
it should be noted that the damping rates can be simply found by considering the interaction between the excitations as a small perturbation and using fermi s golden rule .
this allows one to properly take into account the bogolyubov nature of the thermal excitations .
the damping rates of phonons in a spatially homogeneous bose - condensed gas , in particular for the szepfalusy - kondor mechanism , have been calculated in the recent contributions @xcite .    in order to reach beyond the mean field theory at @xmath18 one
should further develop the popov approach .
one can also proceed along the lines of the beliaev theory , since any physical quantity should be determined by combinations of the beliaev self - energies , which do not contain the infrared singularities . in this paper
we choose the latter way and construct the perturbation theory for a bose - condensed gas , which allows us to find the next to leading order terms ( the terms proportional to @xmath4 ) in the energy spectrum of the elementary excitations . as in @xcite , we consider the excitations in the so - called collisionless regime , where their de broglie wavelength is much smaller than the mean free path of the thermal excitations .
we start with the case of a spatially homogeneous bose - condensed gas and find temperature - dependent energy shifts and damping rates for bogolyubov excitations of any energy . at temperatures @xmath16 the small parameter of the theory proves to be @xmath19 in contrast to @xmath20 for @xmath6 .
the appearance of the extra factor ( @xmath21 ) originates from the bose occupation numbers of thermal excitations with energies of order @xmath13 , which are the most important in the perturbation theory .
as shown below , the damping of excitations with energies @xmath22 is determined by both the szepfalusy - kondor ( @xmath23 ) and beliaev ( @xmath24 ) processes , and can no longer be treated as landau damping .    the theory is generalized for the case of excitations in a spatially inhomogeneous ( trapped ) bose - condensed gas .
a new ingredient here is related to the inhomogeneous density profile of the condensate and the discrete structure of the excitation spectrum .
we develop the method of calculating the self - energy functions and derive the equations for finding the wavefunctions and energies of the excitations ( generalized bogolyubov - de gennes equations ) .
the use of the theory is demonstrated by two examples .
the first one concerns quasiclassical low - energy excitations of a trapped bose - condensed gas in the thomas - fermi regime . the term `` low - energy '' assumes that the excitation energy @xmath25 is much smaller than the mean field interparticle interaction @xmath26 ( @xmath27 is the maximum condensate density ) , and the quasiclassical character of the excitations requires the condition @xmath28 , where @xmath29 is the characteristic trap frequency .
we consider anisotropic harmonic traps , where the discrete structure of the excitation spectrum is not important ( see below and in @xcite ) . on the contrary
, the inhomogeneity of the condensate density profile has a crucial consequence for the damping rates and energy shifts of quasiclassical low - energy excitations .
the most important turns out to be the boundary region of the condensate , where @xmath30 @xcite .
therefore , the result for the damping rates and energy shifts is completely different from that in spatially homogeneous gases .
finally , we analyze the frequency shifts and damping of axially propagating sound waves in cylindrical bose condensates .
as found , the nature of damping is similar to that in the case of phonons in spatially homogeneous bose condensates .
we show that the attenuation of axially propagating sound wave packets in the recent mit experiment @xcite can be well explained as a consequence of this damping .
we consider a weakly interacting bose - condensed gas confined in an external potential @xmath31 .
the grand canonical hamiltonian of the gas can be written as @xmath32 , where ( hereinafter @xmath33 ) @xmath34 and the term @xmath35 assumes a point interaction between atoms .
the field operator of atoms @xmath36 can be represented as the sum of the above - condensate part @xmath37 and the condensate wavefunction @xmath38 which is a _ c_-number . as the interparticle interaction @xmath39 contains both terms conserving the number of above - condensate particles and terms transferring two above - condensate particles to the condensate ( or two condensate particles to the above - condensate part ) , the diagram technique should include both the normal green function @xmath40 and the anomalous green function @xmath41 ( see , e.g. @xcite ) .
the sums of the contributions of all irreducible diagrams will be represented by the normal ( @xmath42 ) and anomalous ( @xmath43 ) self - energies ( see fig.1 and fig.2 ) .
the former corresponds to the processes conserving the number of above - condensate particles , and the latter describes absorption ( or emission ) of two particles to ( out of ) the condensate .
the green function and self - energy operators satisfy beliaev - dyson equations @xcite @xmath44 where the green functions @xmath45 and @xmath46 describe forward and backward propagation of a particle characterized by the hamiltonian @xmath47 .
we confine ourselves to the case of repulsive interaction between the atoms ( @xmath48 ) . to develop the finite temperature perturbation theory for calculating dynamic properties and finding the excitation spectrum of a weakly interacting bose - condensed gas we will use the non - equilibrium generalization @xcite of the matsubara diagram technique . in eqs .
( [ boperg]),([boperf ] ) we perform an analytical continuation of the matsubara frequencies @xmath49 ( @xmath50 is an integer number ) to the upper half - plane , which corresponds to the replacement @xmath51 . then , multiplying both sides of eqs .
( [ boperg ] ) and ( [ boperf ] ) by @xmath52 and @xmath53 , respectively , we arrive at the system of equations in the frequency - coordinate representation : @xmath54 { \cal g}% ( \varepsilon ; { \bf r},{\bf r}^{\prime } ) + \sigma _ { a}(\varepsilon ) { \cal f}% ( \varepsilon ; { \bf r},{\bf r}^{\prime } ) + \delta ( { \bf r}-{\bf r}^{\prime } ) , \label{bcoordg } \\
-\varepsilon { \cal f}(\varepsilon ; { \bf r},{\bf r}^{\prime } ) & = & \left [ -% \frac{\delta } { 2}+v({\bf r})-\mu + \sigma ^{+}(\varepsilon ) \right ] { \cal f}% ( \varepsilon ; { \bf r},{\bf r}^{\prime } ) + \sigma _ { a}(\varepsilon ) { \cal g}% ( \varepsilon ; { \bf r},{\bf r}^{\prime } ) .   \label{bcoordf}\end{aligned}\ ] ] here the action of the integral self - energy operators on the green functions is written in a compact form @xmath55 ( and similar relations for the other combinations ) .
( [ bcoordg]),([bcoordf ] ) should be completed by a generalized gross - pitaevskii equation for the condensate wavefunction : @xmath56 \psi _ { 0}({\bf r})=0   \label{gpeq}\ ] ] and by the normalization condition @xmath57 where @xmath58 is the condensate density , @xmath59 is the density of above - condensate particles , and @xmath60 the total number of particles in the gas .
( [ bcoordg]),([bcoordf ] ) can be solved by using the bogolyubov transformation for the green functions : @xmath61 where the index @xmath7 stands for the set of quantum numbers , and the functions @xmath62 , @xmath63 satisfy generalized bogolyubov - de gennes equations @xmath64 \!\!u_{\nu}({\bf r } ) \!-\!\!\sigma _ { a}(\!\varepsilon_{\nu}\!)v_{\nu}({\bf r } ) , \!\!\ !
\label{bdggenu } \\ \!-\varepsilon_{\nu}v_{\nu}({\bf r})\!&\!=\!&\!\left[\!-\frac{\delta } { 2}\ ! + \!v({\bf r})\!-\!\mu\!+\!\!\sigma ^{+}\!(\!\varepsilon_{\nu}\!)\!\right]\ ! \!v_{\nu}({\bf r})\!-\ ! \!\sigma_{a}(\!\varepsilon_{\nu}\!)u_{\nu}({\bf r}).\!\!\ !   \label{bdggenv}\end{aligned}\ ] ]    in the bogolyubov - de gennes approach only the terms bilinear in @xmath65 operators are retained in the interaction hamiltonian @xmath39 , which assumes that the condensate density is much larger than the density of above - condensate particles .
then , the self - energy operators take the form @xmath66 the result of their action on the condensate wavefunction @xmath67 and the functions @xmath68 , @xmath69 is reduced to @xmath70 and similar relations for the other combinations .
then , eq.([gpeq ] ) becomes the ordinary gross - pitaevskii equation @xmath71 \psi _ { 0}({\bf r})=0 ,   \label{gpeq0}\ ] ] and eqs .
( [ bcoordg]),([bcoordf ] ) are transformed to the ordinary bogolyubov - de gennes equations @xmath72
u_{\nu}({\bf r})-n_{0}({\bf r } ) \widetilde{u}v_{\nu}({\bf r } ) ,   \label{bdgu } \\
-\varepsilon_{\nu}v_{\nu}({\bf r } ) & = & \left [ -\frac{\delta } { 2}+v({\bf r})- \mu+2n_{0}({\bf r})\widetilde{u}\right ] v_{\nu}({\bf r})-n_{0}({\bf r } ) \widetilde{u}u_{\nu}({\bf r } ) .
\label{bdgv}\end{aligned}\ ] ] taking into account eq.([gpeq0 ] ) , in terms of the functions @xmath73 these equations can be rewritten as @xmath74    for a trapped bose - condensed gas in the thomas - fermi regime , where @xmath75 is much larger then the spacing between the trap levels , the kinetic energy term in eq .
( [ gpeq0 ] ) can be omitted and one has @xcite
@xmath76 if the argument of the square root is positive and zero otherwise . for the low - energy excitations ( @xmath77 ) of thomas - fermi condensates eqs .
( [ fplus0]),([fmin0 ] ) can be reduced to hydrodynamic equations obtained by stringari @xcite and solved in the case of spherically symmetric harmonic potential @xmath78 and for some excitations in a cylindrically symmetric potential .
an analytical method of solving eqs .
( [ fplus0]),([fmin0 ] ) ( or the corresponding hydrodynamic equations ) for the low - energy excitations of thomas - fermi condensates in an anisotropic harmonic potential @xmath31 has been developed in @xcite .
for a spatially homogeneous gas the generalized gross - pitaevskii equation ( [ gpeq ] ) is equivalent to the pines - hugenholtz identity @xcite . in the bogolyubov approach it simply gives @xmath79 , and eqs .
( [ bdgu]),([bdgv ] ) lead to the bogolyubov spectrum @xmath80 where @xmath81 is the momentum of the excitation .    under the condition @xmath82 , for which the bogolyubov approach was originally developed
, one can simply put @xmath83 equal to the total density @xmath84 in eq.([bogspectrum ] ) . for @xmath85 , which can be the case at @xmath86 ,
the dispersion law becomes essentially temperature dependent @xcite . in a spatially homogeneous gas
the temperature dependence predominantly originates just from the presence of above condensate particles , with @xmath87 where @xmath88 is the bec transition temperature .
this leads to the replacement @xmath89 in eq.([sigmas0 ] ) and gives @xmath90 .
the dispersion law will be still given by eq.([bogspectrum ] ) in which the condensate density is now temperature dependent : @xmath91 $ ] .
in this section we present the results for the damping rates and energy shifts of elementary excitations in an infinitely large spatially homogeneous bose - condensed gas .
as one can see from eqs.([bdgu]),([bdgv ] ) , for finding the energy spectrum and wavefunctions of the excitations it is sufficient to calculate the self - energies @xmath42 , @xmath92 and @xmath43 .
we will perform the calculations in the frequency - momentum representation and for physical transparency consider temperatures @xmath93 ( the opposite limiting case has been discussed by popov @xcite with regard to the phonon branch of the excitation spectrum ) . in the zero order approximation in the parameter @xmath94
we have the well - known mean field result : @xmath95 , @xmath96 , with @xmath97 ( see above ) . in this approach
we obtain the bogolyubov quasiparticle excitations with the spectrum ( [ bogspectrum ] ) , which we use in order to calculate the next order in @xmath98 .
the latter is determined by the contribution of diagrams containing one quasiparticle loop @xcite ( see figures 1 and 2 ) .
actually in this approach we represent the hamiltonian as the sum of the ( diagonalized ) bogolyubov hamiltonian and the perturbation @xmath99 originating from @xmath100 ( [ hpert ] ) and containing the terms proportional to @xmath101 and @xmath102 : @xmath103.\ ] ] retaining only the temperature - dependent contributions , after laborious calculations for the normal self - energy we obtain @xmath104 , where @xmath105 here @xmath106 , @xmath107 , @xmath108 , @xmath109 is given by eq.([bogspectrum ] ) , @xmath110 is the equilibrium occupation number , @xmath111 and @xmath112 .
similarly , the correction to the anomalous self - energy is given by @xmath113 the resonant parts @xmath114 originate from the terms where one of the intermediate quasiparticles is created and another one annihilated , and the non - resonant parts @xmath115 from the terms where both intermediate quasiparticles are created ( annihilated ) .
temperature independent terms in the non - resonant parts , found by beliaev @xcite , are omitted in eqs.([sigman])-([sigmaar ] ) .
each of the self - energies ( [ sigman])-([sigmaar ] ) is singular at @xmath116 and at least for small momenta the corrections become larger than the mean field values ( [ sigmas0 ] ) . nevertheless , keeping in mind that any physical quantity is determined by the combinations of the self - energies , which do not contain the infra - red singularities , we will still treat @xmath117 and @xmath118 as perturbations .    for a spatially homogeneous gas the pines - hugenholtz identity @xmath119
gives the first order correction to the chemical potential @xmath120 and the relation between @xmath5 and the chemical potential , @xmath121 , coincides with that found by popov @xcite .
the @xmath122 functions in generalized bogolyubov - de gennes equations ( [ bdggenu ] ) , ( [ bdggenv ] ) can be written as @xmath123 and @xmath124 , and in terms of the functions @xmath125 these equations take the form @xmath126 where @xmath127 considering the terms @xmath128 , @xmath129 in eqs .
( [ form1]),([form2 ] ) as small perturbations we put @xmath130 in the expressions for these quantities , following from eqs .
( [ sigman])-([sigmaar ] ) .
then , solving eqs .
( [ form1]),([form2 ] ) , for the excitation energy we obtain @xmath131 , where @xmath132 _ { p\rightarrow ( \varepsilon _ { p},{\bf p})}\ ! ,   \label{cor}\ ] ]    as @xmath133 and @xmath134 are complex , the correction to the excitation energy has both a real and an imaginary part : @xmath135 .
the former gives the energy shift , and the latter is responsible for damping of the excitations . under the condition @xmath136 a straightforward calculation of eq.([cor ] ) on the basis of eqs .
( [ sigman])-([sigmaar ] ) and ( [ s+]),([s- ] ) , for the phonon branch of the spectrum ( @xmath137 ) yields @xmath138 it is important to emphasize that in this case the energy shift is determined by both resonant and non - resonant terms , whereas the damping rate @xmath139 is described solely by the resonant contributions .
this type of damping originates from quasi - resonant scattering of thermal excitations from a given excitation ( landau damping ) and is absent at @xmath6 .
both the energy shift and the damping rate are determined by the interaction of a given excitation with intermediate quasiparticles having energies @xmath140 .
the damping rate @xmath139 ( [ reshomim ] ) coincides with that found in recent contributions @xcite and contains a slight numerical difference from the earlier szepfalusy - kondor result @xcite .
the energy shift for the phonon branch of the spectrum was also calculated in @xcite . in the latter work
the expansion of the self - energy functions near the point @xmath141 was used and formally divergent integrals were canceling each other in the final expression for the energy shift , which have led to the result by approximately factor 6 smaller then the shift ( [ reshomre ] ) obtained by the exact integration .
eqs.([reshomre]),([reshomim ] ) clearly show that the small parameter of the theory is @xmath142 ( see eq.([small ] ) ) , whereas in the zero temperature approach @xcite the small parameter is @xmath143 .
the presence of an additional large factor @xmath144 at finite temperatures @xmath145 originates from the bose enhancement diagrams containing one quasiparticle loop : compared to the zero - temperature case the contribution of each of these diagrams is multiplied by the bose factor @xmath146^{-1}$ ] ( or @xmath147 ) . as the most important is the contribution of intermediate quasiparticles with energies @xmath148 , for @xmath145 the bose factor @xmath149 .
the criterion similar to eq.([small ] ) was found by popov @xcite as the condition which allows one to use the mean field approach at finite temperatures and to renormalize the theory for reaching beyond this approach .
remarkably , the criterion ( [ small ] ) is fulfilled even at temperatures very close to @xmath150 .
for @xmath151 we have @xmath152 , and eq.([small ] ) gives @xmath153 , which coincides with the well known ginzburg criterion @xcite for the absence of critical fluctuations .
the criterion ( [ small ] ) also ensures that the main contribution to the damping rate originates from the interaction of a given excitation with thermal excitations through the condensate , i.e. , from the first term in @xmath154 ( [ hint ] ) .    at energies @xmath155 ( but still @xmath156 ) the results for the energy shift and damping rate , following from eq.([cor ] ) ,
can be represented in the form @xmath157 the functions @xmath158 and @xmath159 have been calculated numerically .
for @xmath160 we have @xmath161 , @xmath162 and eqs .
( [ homre ] ) , ( [ homim ] ) coincide with eqs .
( [ reshomre ] ) , ( [ reshomim ] ) .
the dependence of @xmath163 and @xmath164 on the excitation energy is presented in fig.3 . for @xmath165 the damping rate increases with @xmath166 and reaches its maximum at @xmath167 .
further increase of @xmath166 leads to slowly decreasing @xmath164 . for single - particle excitations ( @xmath168 )
we obtain @xmath169 .
accordingly , the damping rate can be written as @xmath170 , where @xmath171 is the elastic cross section , and @xmath172 the particle velocity .
this damping rate exceeds the beliaev temperature - independent term @xmath173 even at @xmath174 close to @xmath175 , if @xmath176 .
in contrast to the phonon branch of the spectrum , for @xmath155 the damping is provided by both the szepfalusy - kondor ( @xmath17 ) and beliaev ( @xmath177 ) processes and , hence , can no longer be treated as landau damping .
the small parameter of the theory is still given by eq.([small ] ) , since even at @xmath178 the energy of at least one of the thermal excitations is of order @xmath179 .
the energy shift for @xmath180 is negative .
the modulus of the shift increases with @xmath166 and reaches its maximum at @xmath181 .
the further increase of @xmath166 decreases @xmath182 .
the latter is equal to zero for @xmath183 , and becomes positive at larger @xmath166 .
the above results for the damping rate and energy shift of a given excitation are obtained in the so called collisionless regime : we assume that the de broglie wavelength of the excitation , @xmath184 , is much larger than the mean free path of the thermal quasiparticles with energies @xmath185 , which are mostly responsible for the damping and shifts .
it is also assumed that the excitation energy @xmath166 greatly exceeds the damping rate of these thermal excitations .
the latter is of order @xmath186 ( see fig.3 ) , and for @xmath187 the two requirements of the collisionless regime are well satisfied under condition ( [ small ] ) . in the phonon branch of the excitation spectrum ( @xmath188 ) these requirements
are equivalent to each other , and the collisionless criterion can be simply written as @xmath189 as clearly seen , in the phonon branch one can always find excitations which do not satisfy eq.([cols ] ) and , hence , require a hydrodynamic description with regard to their damping rates and energy shifts .
the collisionless criterion ( [ cols ] ) provides an additional argument on support of the above used perturbative approach for solving eqs .
( [ form1 ] ) , ( [ form2 ] ) . under condition ( [ cols ] ) the term @xmath190 , the term @xmath191 , and the term @xmath192 .
the non - mean - field shift @xmath163 is actually the shift of the excitation energy @xmath109 at a given condensate density @xmath5 . on the other hand , @xmath166 is determined by the bogolyubov dispersion law ( [ bogspectrum ] ) , with the temperature - dependent condensate density @xmath193 , and , hence , is temperature - dependent by itself .
therefore , at a given @xmath174 one will also have the mean - field temperature - dependent energy shift @xmath194 . as the condensate density decreases with increasing temperature , @xmath195 is always negative .
for @xmath145 it greatly exceeds the above calculated shift @xmath163 at any @xmath196 .
the ratio @xmath197 decreases with temperature , but even for @xmath198 one has @xmath199 and @xmath200 .
we now generalize the above obtained results for the case of elementary excitations in a spatially inhomogeneous ( trapped ) bose - condensed gas . as already mentioned in the introduction ,
a new ingredient here is related to the inhomogeneous density profile of the condensate and the discrete structure of the excitation spectrum .
this requires us to develop a new method of calculating the self - energy functions in generalized bogolyubov - de gennes equations ( [ bdggenu ] ) , ( [ bdggenv ] ) .
the self - energy operators in these equations are the sums of the zeroth and first order terms : @xmath201 and a similar relation for @xmath202 . at temperatures
@xmath203 the zero order value of the above - condensate density in the condensate spatial region is coordinate independent and equal to the above - condensate density in the ideal gas approach : @xmath204 . on the contrary , the self - energies @xmath133 , @xmath205
depend explicitly on the condensate density .
due to the discrete structure of the energy spectrum of excitations the expressions for these quantities should be written in the form of sums over the discrete states of intermediate quasiparticles @xmath8,@xmath9 . in the frequency - coordinate representation
we have @xmath206    as we saw in the previous section , in the spatially homogeneous case all physical quantities are determined by the contribution to the self - energy functions @xmath133 , @xmath205 from intermediate quasiparticles with energies of order the mean field interaction between particles .
the same holds for a spatially inhomogeneous ( trapped ) bose - condensed gas in the thomas - fermi regime , where the mean field interaction @xmath26 greatly exceeds the level spacing in the trapping potential .
the intermediate quasiparticles with energies of order @xmath26 are essentially quasiclassical . with regard to the integral operator @xmath207 in the generalized gross - pitaevskii equation ( [ gpeq ] ) , which is solely determined by non - resonant contributions
, this immediately allows one to replace the summation over the discrete intermediate states by integration .
the kernel of this integral operator varies at distances @xmath208 of order the correlation length @xmath209 which is much smaller than the characteristic size of the condensate .
therefore , the result of the operator action on the condensate wavefunction can be written in the local density approximation and , hence , should rely on eq.([mucorr ] ) with coordinate - dependent condensate density @xmath210 : @xmath211_{\varepsilon\rightarrow 0 } \psi_0({\bf r})= \int d^3r'\left[\sigma^{(1)}(\varepsilon,{\bf r},{\bf r'})- \sigma_a^{(1)}(\varepsilon,{\bf r},{\bf r'})\right]_{\varepsilon \rightarrow 0}\psi_0({\bf r'})=-\beta\psi_0 ^ 2({\bf r}).\ ] ] this result can be easily obtained from eqs .
( [ tnsn ] ) , ( [ trsa ] ) , where one should put @xmath212 , neglect the difference between @xmath213 and @xmath214 , and make a summation over @xmath9 . replacing the summation over @xmath8 by integration one
should also take into account that for quasiclassical excitations the functions @xmath215 can be represented in the form @xmath216 where @xmath217 is the ratio of the local to total density of states for bogolyubov quasiparticles of a given symmetry , described by the classical hamiltonian @xmath218    on the basis of eq.([missed ] ) we obtain the generalized gross - pitaevskii equation in the form @xmath219 where @xmath220 is coordinate independent .
compared to the ordinary gross - pitaevskii equation ( [ gpeq0 ] ) , eq.([gpeq1 ] ) contains an extra term @xmath221\psi_0 $ ] in the lhs .
one can easily check that eq.([gpeq1 ] ) coincides with the equation @xmath222\psi_0=0\ ] ] obtained by averaging the non - linear schrdinger equation for the field operator . as mentioned above , for @xmath145 the above - condensate density @xmath223 in the condensate spatial region
is mainly determined by the coordinate - independent ( ideal gas ) contribution @xmath224 .
the coordinate - dependent correction to the above - condensate density , @xmath225 , turns out to be equal to the anomalous average @xmath226 : @xmath227 accordingly , the quantity @xmath228\psi_0= \widetilde{u}(2\langle\hat\psi^{\prime\dagger}\hat\psi^{\prime } \rangle + \langle\psi^{\prime}\psi^{\prime}\rangle)$ ] .    taking advantage of eqs  ( [ 0 ] ) , ( [ 00 ] ) and
( [ gpeq1 ] ) , the generalized bogolyubov - de gennes equations ( [ bdggenu ] ) , ( [ bdggenv ] ) are reduced to @xmath229 where the quantities @xmath230 , @xmath231 are given by eqs .
( [ s+ ] ) , ( [ s- ] ) , and @xmath25 is the exact value of the excitation energy . eqs .
( [ gpeq1 ] ) , ( [ fplus ] ) and ( [ fmin ] ) represent a complete set of equations for finding the energy shifts and damping rates of the elementary excitations .    a precise calculation of the self - energy functions in eqs .
( [ fplus ] ) , ( [ fmin ] ) depends on the value of @xmath25 and on the trapping geometry . in this section
we will make general statements on how the calculation can be performed . in most of the cases ( except the case of the lowest excitations with zero orbital angular momentum in spherically symmetric traps )
the characteristic time scale in the self - energy operators , @xmath232 , is much smaller than the inverse level spacing in the trap .
therefore , the summation over the discrete intermediate states can be replaced by integration .
this is a direct consequence of the general statement that the time - dependent discrete fourier sum can be replaced by its integral representation at times much smaller than the inverse frequency spacing ( see e.g. , @xcite ) .
the kernels of the non - resonant parts of the self - energy operators , ( [ tnsn ] ) and ( [ tnsa ] ) , vary at distances @xmath208 which do not exceed the correlation length @xmath209 .
as @xmath233 is much smaller than the characteristic size of the condensate , the non - resonant parts of the self - energies can be calculated in the local density approximation .
one can easily find from eqs .
( [ tnsn ] ) , ( [ tnsa ] ) that they lead to the result of eqs .
( [ sigman ] ) , ( [ sigmaan ] ) ( with @xmath5 replaced by the coordinate - dependent density @xmath210 ) , multiplied by @xmath234 . for the lowest excitations @xmath7 one should put @xmath235 in eqs .
( [ sigman]),([sigmaan ] ) , and for quasiclassical excitations @xmath7 take @xmath196 from the bogolyubov dispersion law @xmath236 .
the calculation of the resonant contributions to the self - energies is more subtle .
using eq.([form ] ) for the functions @xmath237 in eqs .
( [ trsn ] ) , ( [ trsa ] ) , one can see that all resonant contributions contain the quantity @xmath238 writing @xmath239 as the integral over time @xmath240 , we obtain @xmath241 where the quantum - mechanical correlation function @xmath242 we will turn from the integration over the quantum states @xmath9 of the quasiclassical thermal excitations to the integration along the classical trajectories of motion of bogolyubov - type quasiparticles in the trap . following a general method ( see @xcite ) , employed in @xcite for the damping of low - energy excitations , we obtain @xmath243 where @xmath244 is the coordinate of the classical trajectory with initial momentum @xmath81 and coordinate @xmath245 .
eq.([k2 ] ) will be used in the next sections where we demonstrate the facilities of the theory .    concluding this section
, we emphasize the key role of harmonicity of the trapping potential for temperature - dependent energy shifts of the excitations .
as mentioned in the previous section , in the spatially homogeneous case at a given temperature the non - mean - field shift is much smaller than the shift @xmath246 appearing in the mean field approach simply due to the temperature dependence of the condensate density in the bogolyubov dispersion law ( [ bogspectrum ] ) . for the thomas - fermi condensate in a harmonic confining potential the situation is different . in this case
the spectrum of low - energy ( @xmath77 ) excitations is independent of the mean field interparticle interaction @xmath26 ( chemical potential ) and the condensate density profile@xcite .
hence , the temperature - dependent energy shifts can only appear due to non - thomas - fermi corrections . for finding these corrections one should use the mean field self - energies @xmath247 , @xmath248 , where the only difference from the @xmath6 case is related to the presence of above - condensate particles in the condensate spatial region at finite @xmath174 through the coordinate - independent term @xmath249 in @xmath250
. then eqs .
( [ fplus ] ) , ( [ fmin ] ) take the form of ordinary bogolyubov - de gennes equations ( [ fplus0]),([fmin0 ] ) , and eq.([gpeq1 ] ) becomes the ordinary gross - pitaevskii equation ( [ gpeq0 ] ) , with the chemical potential @xmath251 replaced by @xmath252 .
the latter circumstance changes the condensate wavefunction compared to that at @xmath6 and ensures the temperature dependence of @xmath253 .
accordingly , the excitation energies @xmath25 in eqs .
( [ fplus0]),([fmin0 ] ) also become temperature dependent .
this type of approach , which for a spatially homogeneous gas would immediately lead to the result of lee and yang @xcite , has been used in recent numerical calculations of the energy shifts of the lowest quadrupole excitations in spherically symmetric @xcite and cylindrically symmetric @xcite harmonic traps .
the presence of the coordinate - dependent part of the above - condensate density , @xmath254 , in these calculations is not adequate , since the anomalous average equal to this part was omitted and equations for the excitations did not contain the corrections to the self - energies , also proportional to @xmath4 .
however , at @xmath145 , where @xmath255 , the coordinate - dependent part @xmath254 as itself should not significantly influence the result , and the calculations @xcite should actually demonstrate how important are the mean field non - thomas - fermi effects .
the results of @xcite show the absence of energy shifts of the excitations at temperatures @xmath256 in the jila experiment @xcite and in this sense agree with the experimental data , but do not describe the upward and downward shifts of the excitation energies , observed experimentally at higher temperatures ( in this respect it is worth mentioning that the calculations @xcite performed for the thermal cloud in the hydrodynamic regime agree surprisingly well with the experiment @xcite ) . on the other hand ,
the calculation @xcite shows a downward shift of the energy of the lowest quadrupole excitation with increasing temperature in the conditions of the mit experiment @xcite .
this is consistent with the experimental data and indicates that for not very small thomas - fermi parameter @xmath257 the mean field non - thomas - fermi effects can be important for temperature - dependent shifts of the lowest excitations .
below we will assume a sufficiently small thomas - fermi parameter @xmath257 and demonstrate the use of the theory by the examples where the influence of non - thomas - fermi effects on the energy shifts of the excitations is not important .
we will discuss the thomas - fermi condensates in a harmonic confining potential on the basis of eqs .
( [ gpeq1 ] ) , ( [ fplus ] ) , ( [ fmin ] ) . neglecting the kinetic energy term in eq.([gpeq1 ] ) , we arrive at a quadratic equation for @xmath253 . expanding the solution of this equation in powers of @xmath258 and retaining only the terms independent of @xmath258 and the terms linear in @xmath258 , for the condensate density
we obtain @xmath259 where @xmath260 is the density of the thomas - fermi condensate in the ordinary mean field approach .
we first consider the damping and energy shifts of quasiclassical ( @xmath261 ) low - energy excitations of a trapped thomas - fermi condensate , i.e. , the quasiclassical excitations with energies much smaller than the mean field interaction between particles @xmath26 . in this case
the terms in eqs .
( [ fplus ] ) and ( [ fmin ] ) , originating from the kinetic energy of the condensate , can be omitted from the very beginning @xcite .
then , using eq.([cd1 ] ) and treating the terms containing @xmath231 and @xmath230 as perturbations , we obtain @xmath262 , where @xmath263 is the excitation energy in the mean field approach , and the correction to the excitation energy @xmath264 is given by the relation @xmath265 here @xmath266 are the zero - order wavefunctions of the excitations , determined by the ordinary bogolyubov - de gennes equations  ( [ fplus0 ] ) , ( [ fmin0 ] ) , with @xmath267 .    in the case of quasiclassical excitations also the kernels of resonant parts of integral operators in eq.([pertenergy ] ) vary on a distance scale @xmath208 which does not exceed the correlation length @xmath233 .
this can be already seen from eqs .
( [ q ] ) , ( [ k2 ] ) : the characteristic time scale @xmath232 in eq.([q ] ) is much shorter than @xmath268 and important is only a small part of the classical trajectory , where the condensate density is practically constant and @xmath269 , with @xmath270 .
the correlation length @xmath233 is not only much smaller than the size of the condensate , but also smaller than the width of the boundary region of the condensate , where @xmath271 .
therefore , the action of all integral operators on the functions @xmath266 in eq.([pertenergy ] ) can be calculated in the local density approximation .
accordingly , for each of these operators one can use the quantity following from eqs .
( [ sigman])-([sigmaar ] ) , with @xmath272 and @xmath196 from the bogolyubov dispersion law @xmath236 .
then , using eqs .
( [ form ] ) we can express the energy shift @xmath273 and the damping rate @xmath274 through the energy shift @xmath275 and damping rate @xmath276 of the excitation of energy @xmath25 in a spatially homogeneous bose - condensed gas with the condensate density equal to @xmath277 : @xmath278 the second term in the integrand of eq.([ecorghoms ] ) originates from the temperature dependence of the shape of the condensate wavefunction .
for any ratio @xmath279 this positive term dominates over the negative term @xmath275 .
the latter circumstance can be easily established from the results for @xmath280 in fig.3 .
thus , for quasiclassical low - energy excitations the energy shift @xmath273 will be always positive , irrespective of the trapping geometry and the symmetry of the excitation .
we confine ourselves to the case of cylindrical symmetry , where for the states with zero angular momentum one finds @xmath281 with @xmath282 , @xmath283 being the characteristic size of the condensate in the radial and axial direction , @xmath284 , @xmath285 the radial and axial frequencies , and @xmath286 the radial coordinate .
the main contribution to the integral in eq.([ecorqhom ] ) comes from the boundary region of the condensate , where @xmath287 . from eqs .
( [ fplus]),([fmin ] ) one can easily see that in this region the possibility to omit the non - thomas - fermi effects originating from the kinetic energy of the condensate requires the condition @xmath288 .
this condition ensures that the characteristic width of the boundary region greatly exceeds the excitation wavelength , and we arrive at the following relations for the energy shifts and damping rates of the excitations : @xmath289    it is important to emphasize that in the boundary region of the condensate , responsible for the energy shifts and damping rates of the quasiclassical excitations , the quantities @xmath290 , @xmath291 , and @xmath205 are determined by the contribution of intermediate quasiparticles which have energies comparable with @xmath25 . moreover ,
in this spatial region the quasiparticle energies are of order the local mean field interparticle interaction . as a consequence , the energy shift @xmath273 ( [ qcshift ] ) and the damping rate @xmath274 ( [ qcdamp ] ) are practically independent of the condensate density profile .
for the same reason the damping rate is determined by both the szepfalusy - kondor and beliaev damping processes .
therefore , similarly to the damping of excitations with energies @xmath292 in a spatially homogeneous gas , the damping of quasiclassical low - energy excitations of a trapped bose - condensed gas can no longer be treated as landau damping .
the derivation of eqs .
( [ qcshift ] ) , ( [ qcdamp ] ) assumes that the motion of the excitation @xmath7 is quasiclassical for all degrees of freedom .
we now turn to the condensate excitations in cigar - shaped cylindrical traps , which are quasiclassical only in the axial direction and correspond to the lowest modes of the radial motion .
we will consider low - energy excitations ( @xmath77 ) , i.e. , the excitations with the axial wavelength much larger than the correlation length @xmath233 .
in the recent mit experiment @xcite localized excitations of this type were created in the center of the trap by modifying the trapping potential using the dipole force of a focused off - resonant laser beam .
then , a wave packet traveling along the axis of the cylindrical trap ( axially propagating sound wave ) was observed . in the mean field approach
the sound waves propagating in an infinitely long ( axially homogeneous ) cylindrical bose condensate have been discussed in @xcite .    for revealing the key features of the non - mean - field effects ( damping and the change of the sound velocity ) we confine ourselves to the same trapping geometry .
with regard to realistic cylindrical traps this will be a good approach if the mean free path of sound waves is smaller than the characteristic axial size of the sample . as found in @xcite , for axially propagating sound waves radial oscillations of the condensate
are absent , and the wavefunctions @xmath293 , with @xmath294 and @xmath295 being the axial momentum .
the dispersion law @xmath296 is characterized by the sound velocity equal to @xmath297 , where @xmath298 is the maximum density of the thomas - fermi condensate in the ordinary mean field approach .
it should be noted from the very beginning that , according to eq.([cd1 ] ) , @xmath299 is related to the corrected value of the maximum condensate density @xmath27 as @xmath300 .
therefore , being interested in the sound velocity at a given value of the maximum condensate density , one should substitute this expression to eq.([disp ] ) .
this immediately changes the sound velocity to @xmath301 in the leading term ( [ disp ] ) of the dispersion law and provides a contribution to the frequency shift of the sound wave @xmath302    the damping rate and other contributions to the frequency shift can be found directly from eq.([pertenergy ] ) by using the wavefunctions @xmath303 ( [ fk ] ) .
the intermediate quasiparticles giving the main contribution to the damping rate and frequency shift have energies @xmath304 , i.e. much larger than the frequency of the considered sound wave , @xmath305 ( see below ) .
therefore , similarly to the case of phonons in a spatially homogeneous condensate , the non - resonant terms ( [ tnsn ] ) , ( [ tnsa ] ) contribute only to the frequency shift . as already mentioned above ,
the characteristic distance scale @xmath306 in the kernels of the self - energies ( [ tnsn ] ) , ( [ tnsa ] ) is of order the correlation length @xmath233 , and the sum of their contributions to the frequency shift , @xmath307 , can be calculated by using the local density approximation for the action of the self - energy operators on the functions @xmath308 . as a result , we express @xmath307 through the non - resonant part of the energy shift @xmath309 in the spatially homogeneous condensate with the condensate density @xmath310 : @xmath311 for @xmath312 one can directly find from eqs .
( [ sigman ] ) , ( [ sigmaan ] ) that @xmath313 .
then we obtain @xmath314    the resonant terms ( [ trsn ] ) , ( [ trsa ] ) contribute to both the frequency shift and damping rate .
this means that the latter is determined by the szepfalusy - kondor scattering processes and , since the characteristic energies of intermediate quasiparticles are much larger than @xmath305 , can be treated as landau damping .
the resonant contributions to the frequency shift and damping rate can not be found in the local density approximation , as the characteristic distance scale @xmath306 in the kernels of the self - energy operators in eq.([pertenergy ] ) is of order the radial size of the condensate . for finding these contributions
one has to substitute the resonant parts of the self - energies , ( [ tnsn])-([trsa ] ) , to eq.([pertenergy ] ) and , by using eqs.([q])-([k2 ] ) , turn from summation over quasiclassical states @xmath8 , @xmath9 of intermediate quasiparticles to the integration along classical trajectories of their motion .
then , a direct calculation yields ( cf .
@xcite ) @xmath315 where @xmath316 is the classical trajectory starting at the phase space points @xmath317 on the ( hyper)surface of constant energy @xmath213 , @xmath318 , and @xmath319    generally speaking , the integration in eq.([deltak1 ] ) is a tedious task as it requires a full knowledge of the classical trajectories on a time scale @xmath320 .
this is also the case in the idealized cylindrical trap , because of coupling between the radial and axial degrees of freedom .
we will rely on the approach which assumes a fast radial motion of quasiparticles compared to their motion in the axial direction and , hence , requires the frequency of the sound wave , @xmath305 , significantly smaller than the radial frequency @xmath284 . then on a time scale @xmath320 the quasiparticles with energies @xmath321 ( which are the most important for the energy shifts and damping of the sound wave ) oscillate many times in the radial direction ,
whereas their axial variables @xmath322 , @xmath323 only slightly change and , hence , can be adiabatically separated from the fast radial variables @xmath324 , @xmath325 . in this case
it is convenient to integrate eq,([deltak1 ] ) over @xmath326 and , using eq.([fk ] ) , represent it in the form @xmath327 where the integration is performed over the entire classically accessible region of the phase space .    since @xmath322 is close to @xmath328 , in the exponent of the integrand we can write @xmath329 , where the axial velocity @xmath330 is obtained from the exact hamiltonian equations of motion by averaging over the fast radial variables : @xmath331 . for the classical radial motion ( @xmath332 ) the averaging procedure simply reduces to integration over @xmath333 under the condition @xmath334 at fixed values of @xmath213 , @xmath335 and @xmath328 , with the weight proportional to the local density of states for the radial motion : @xmath336 where @xmath337 . finally , averaging the function @xmath338 over the fast radial variables and integrating over @xmath339 in eq.([interm ] ) , we obtain @xmath340 the resonant contribution to the frequency shift , given by the real part of eq.([dkf ] ) , after the integration proves to be @xmath341 .
the sum of this quantity with the non - resonant term ( [ dnonres ] ) and @xmath342 ( [ s1 ] ) leads to the frequency shift of the sound wave @xmath343 the imaginary part of eq.([dkf ] ) gives the damping rate @xmath344    except for the numerical coefficients , eqs,([sndshift ] ) , ( [ sd ] ) are similar to eq.([reshomre ] ) , ( [ reshomim ] ) for the damping rate and energy shift of phonons in a spatially homogeneous bose condensate .
this a consequence of the fact that the condensate boundary region practically does not contribute to the damping rate and frequency shift of axially propagating sound waves , in contrast to the case of excitations quasiclassical for both axial and radial degrees of freedom .    in the mit experiment @xcite the characteristic spatial size of created localized excitations was @xmath345 @xmath251 m and , accordingly , so was the initial wavelength of propagating sound .
according to the experimental data , the propagating pulse died out during 25 ms , and after that only the lowest quadrupole excitation characterized by a much loner damping time ( @xmath346 ms ) was observed .
we believe that the attenuation of axially propagating sound in the mit experiment @xcite on the time scale of 25 ms can be well explained as a consequence of damping .
the characteristic frequency of the waves in the packet can be estimated as @xmath347 .
then eq.([sd ] ) gives the damping rate independent of the condensate density @xmath299 . in the mit experiment
the temperature @xmath3480.5 @xmath251k was roughly only twice as large as @xmath26 , which decreases the damping rate by approximately @xmath349 compared to that given by eq.([sd ] ) . in these conditions
we obtain a characteristic damping time of @xmath350 ms , relatively close to the measured value .
the relative change of the sound velocity , @xmath351 , increases with decreasing condensate density @xmath27 .
however , even at the lowest densities of the mit experiment @xcite ( @xmath352 @xmath353 ) the quantity @xmath354 does not exceed @xmath355 and is practically invisible .
we acknowledge fruitful discussions with m.w .
reynolds , m.a . baranov and a. griffin . this work was supported by the dutch foundation fom and by the russian foundation for basic studies .