for ultracold dilute non - relativistic bosons in three dimensions , bogoliubov theory gives a successful description of most quantities of interest @xcite .
this approximation breaks down , however , near the critical temperature for the phase transition , as well as for the low temperature phase in lower dimensional systems , due to the importance of fluctuations .
one would therefore like to have a systematic extension beyond the bogoliubov theory , which includes the fluctuation effects beyond the lowest order in a perturbative expansion in the scattering length .
such extensions have encountered obstacles in the form of infrared divergences in various expansions @xcite . only recently , a satisfactory framework has been found to cure these problems @xcite .
functional renormalization @xcite for the average action @xcite systematically copes with the infrared problems by exploring the difficult long range behavior gradually by means of non - perturbative flow equations , which are based on an exact renormalization group equation @xcite . for zero temperature
they lead to a consistent description of bose - einstein condensation and the quantum phase transition for nonrelativistic bosons in arbitrary dimension @xcite .
a key ingredient is the generation of a term in the inverse propagator that is quadratic in the frequency @xmath4 .
it characterizes the long time behavior of the full propagator at long distances , even if no quadratic frequency dependence is present in the microphysical or classical propagator .
the coupling @xmath5 is due to quantum fluctuations and appears during the flow from microphysics to macrophysics .
the extreme infrared limit is governed by a `` relativistic '' model where @xmath5 dominates , with an enhanced ( approximate ) space - time symmetry @xmath6 corresponding to lorentz symmetry . in one or two dimensions
the relativistic term @xmath7 in the inverse propagator is crucial for a correct description of the low temperature behavior .
in contrast , in three dimensions the flow towards the extreme infrared behavior is only logarithmic , such that the quantitative influence remains quite moderate for practical purposes .
nevertheless , the effective coupling @xmath5 is needed in order to avoid the infrared problems .    in this paper , we extend this formalism to a nonvanishing temperature
. we present a quantitative rather accurate picture of bose - einstein condensation in three dimensions and find that the bogoliubov approximation is indeed valid for most quantities .
the same method can be applied for one or two dimensions , such that the present work can also serve as a test of the method .
this is nontrivial , because insufficient truncations of the flow equations can lead to fake dependencies on the microscopic physics .    for dilute non - relativistic bosons in three dimensions
we find an upper bound on the scattering length @xmath0 .
this is similar to the `` triviality bound '' for the higgs scalar in the standard model of elementary particle physics . as a consequence ,
the scattering length is at most of the order of the inverse effective ultraviolet cutoff @xmath8 , which indicates the breakdown of the pointlike approximation for the interaction at short distances .
typically , @xmath8 is of the order of the range of the van der waals interaction . for dilute gases , where the interparticle distance @xmath9 is much larger than @xmath8
, we therefore always find a small concentration @xmath10 .
this provides for a small dimensionless parameter , and perturbation theory in @xmath11 becomes rather accurate for most quantities . for typical experiments with ultracold bosonic alkali atoms
one has @xmath12 , @xmath13 , such that @xmath14 is really quite small .
bosons with pointlike interactions can also be employed for an effective description of many quantum phase transitions at zero temperature , or phase transitions at low temperature @xmath15 . in this case , they correspond to quasi - particles , and their dispersion relation may differ from the one of non - relativistic bosons , @xmath16 . we describe the quantum phase transitions for a general microscopic dispersion relation , where the inverse classical propagator in momentum and frequency space takes the form @xmath17 ( in units where the particle mass @xmath18 is set to @xmath19 ) .
we present the quantum phase diagram at @xmath20 in dependence on the scattering length @xmath0 and a dimensionless parameter @xmath21 , which measures the relative strength of the term quadratic in @xmath22 in @xmath23 . in the limit @xmath24 ( @xmath25 )
our model describes relativistic bosons .
a nonzero temperature @xmath15 introduces a new scale into the problem .
this modifies the effective scaling behavior .
in particular , near the critical temperature @xmath2 , where the bose - einstein condensate dissolves , the long distance physics can be described by classical statistics .
the non - relativistic bosons belong to the universality class of the three - dimensional @xmath26 model , with the associated universal critical exponents .
the flow equations in the classical regime are well studied and lead , for an appropriate level of the truncation , to very accurate results for the critical behavior @xcite .
here we use a rather simple truncation , but we obtain nevertheless a reasonable description of the non - analytical critical behavior . however , the temperature range for the applicability of the universal critical behavior is found to be very small for the small values of @xmath11 encountered for the dilute gases .
the value of the critical temperature for interacting bosons can not be computed within the bogoliubov theory .
we find a temperature shift @xmath3 , compared to the free theory .
in contrast to other estimates , which only take the classical fluctuations into account ( only the zero matsubara frequency ) , we include the full quantum statistics ( all matsubara frequencies ) .
we finally compute the sound velocity at @xmath20 .
this quantity agrees to high precision with the bogoliubov result . for a known scattering length @xmath0
, it can therefore be used as a precise measure of the density .
our paper is organized as follows . as a starting point , we specify our microscopic model in sec .
[ sectmicroscopicmodel ] . in the following sec .
[ sectnonperturbativeflowequations ] , we recall the method of functional renormalization , explain our truncation of the effective average action and also describe the projection of the exact flow equation onto that truncation .
the flow equations in the vacuum are investigated in sec .
[ sectupperboundforthescatteringlength ] and an upper bound for the scattering length is derived . in sec .
[ sectquantumphasediagram ] we explain our method to determine the density and also study the quantum phase diagram at zero temperature .
the effects of a non - vanishing temperature on the condensate are taken into account in sec .
[ secttemperaturedependenceofcondensate ] , while sec .
[ sectcriticaltemperature ] gives our result for the critical temperature in dependence of the interaction strength .
finally , we investigate the sound velocity and draw conclusions in secs .
[ sectsoundvelocity ] and [ sectconclusions ] .
the appendices [ secttruncation ] and [ sectsymmetriesandnoethercurrents ] contain a motivation of our truncation in terms of a systematic derivative expansion and an analysis of symmetries constraining the form of the effective action .
our explicit results for the flow equations of the effective potential and the kinetic coefficients are shown in appendices [ secfloweffectivepotential ] and [ sectflowofkineticcoefficients ] .
appendix [ sectpropdisp ] analyzes the dispersion relation .
our microscopic action describes nonrelativistic bosons , with an effective interaction between two particles given by a contact potential .
it is assumed to be valid on length scales where the microscopic details of the interaction are irrelevant and the scattering length is sufficient to characterize the interaction .
the microscopic action reads @xmath27=\int_x \,{\big \{}\phi^*\,(s\partial_\tau - v\partial_\tau^2-\delta-\sigma)\,\phi\,+\,\frac{1}{2}\lambda(\phi^*\phi)^2{\big \ } } , \label{microscopicaction}\ ] ] with @xmath28 the integration goes over the whole space as well as over the imaginary time @xmath29 , which at finite temperature is integrated on a circle of circumference @xmath30 according to the matsubara formalism .
we use natural units @xmath31 .
we also scale time and energy units with appropriate powers of @xmath32 , with @xmath18 the particle mass . in other words ,
our time units are set such that effectively @xmath33 .
in particular , we use the rescaled chemical potential @xmath34 and @xmath15 stands for the temperature multiplied by @xmath32 . in these units time
has the dimension of length squared . for standard non - relativistic bosons
one has @xmath35 and @xmath36 , but we also consider quasiparticles with a more general dispersion relation described by nonzero @xmath37 .
after fourier transformation , the kinetic term reads @xmath38 with @xmath39 at nonzero temperature , the frequency @xmath40 is discrete , with @xmath41 while at zero temperature this becomes @xmath42 the dispersion relation encoded in eq
. obtains by analytic continuation @xmath43    in this paper , we consider homogeneous situations , i.e. an infinitely large volume without a trapping potential .
many of our results can be translated to the inhomogeneous case in the framework of the local density approximation .
one assumes that the length scale relevant for the quantum and statistical fluctuations is much smaller than the characteristic length scale of the trap . in this case ,
our results can be transferred by taking the chemical potential position dependent in the form @xmath44 , where @xmath45 is the trapping potential .
the microscopic action is invariant under the global @xmath46 symmetry which is associated to the conserved particle number , @xmath47 on the classical level , this symmetry is broken spontaneously when the chemical potential @xmath48 is positive . in this case , the minimum of @xmath49 is situated at @xmath50 .
the ground state of the system is then characterized by a macroscopic field @xmath51 , with @xmath52 .
it singles out a direction in the complex plane and thus breaks the @xmath46 symmetry .
nevertheless , the action itself and all modifications due to quantum and statistical fluctuations respect the symmetry . for @xmath35 and @xmath36 ,
the situation is similar for galilean invariance . at zero temperature
, we can perform an analytic continuation to real time and the microscopic action is then invariant under transformations that correspond to a change of the reference frame in the sense of a galilean boost .
it is easy to see that in the phase with spontaneous @xmath46 symmetry breaking also the galilean symmetry is broken spontaneously : a condensate wave function , that is homogeneous in space and time , would be represented in momentum space by @xmath53 under a galilean boost transformation with a boost velocity @xmath54 , this would transform according to @xmath55 this shows that the ground state is not invariant under such a change of reference frame .
this situation is in contrast to the case of a relativistic bose - einstein condensate , like the higgs boson field after electroweak symmetry breaking .
a relativistic scalar transforms under lorentz boost transformations according to @xmath56 such that a condensate wave function @xmath57 transforms into itself .
we will investigate the implications of galilean symmetry for the form of the effective action in app .
[ secttruncation ] .
an analysis of general coordinate invariance in nonrelativistic field theory can be found in @xcite .
we start with a functional integral representation of the grand canonical partition function @xmath58}. \label{eqpartfunction}\ ] ] in this paper , we work with the formalism of quantum statistics for many particle problems .
in contrast to classical statistics , the fields @xmath59 are parametrized by an euclidean time variable @xmath29 in addition to the space variable @xmath60 .
this euclidean time is wrapped up on a circle of circumference @xmath61 , such that the fields @xmath59 live on a ( generalized ) torus .
the microscopic action @xmath62 $ ] is of a form similar to .
we generalize eq . by introducing a source @xmath63 for the fields @xmath64 and
write @xmath65=e^{w[j]}=\int d \chi \,e^{-s[\chi]+\int_x j\chi}. \label{eqpartfunctwithj}\ ] ] the @xmath1-point correlation functions can now be obtained by functional differentiation of @xmath66 $ ] , while @xmath67 $ ] generates the connected @xmath1-point functions .
for example , the expectation value of @xmath64 obtains from @xmath68}{\delta j(x)}.\ ] ] the thermodynamic potential @xmath69 associated with the grand canonical partition function is given by @xmath70.\ ] ]    the effective action @xmath71 $ ] is defined as @xmath72=\left(-w[j]+\int_x j \phi\right)_{j = j_{\text{ex}}[\phi]}.\ ] ] where @xmath73 is obtained by the inversion of @xmath74}{\delta   j}\bigg{|}_{j = j_\text{ex}}=\phi.\ ] ] it is straightforward to show @xmath75}{\delta \phi}=j$ ] . in the absence of a source @xmath63 , we obtain for the thermodynamic potential of the grand canonical partition function @xmath76\ ] ] where @xmath77 is defined by @xmath78}{\delta \phi}\bigg{|}_{\phi=\phi_{\text{eq}}}=0.\ ] ] the effective action @xmath71 $ ] is the generating functional of the one - point irreducible correlation functions and in a sense , its precise knowledge corresponds to the solution of the theory .
our method determines @xmath71 $ ] with the help of an exact flow equation . for that purpose
, we include an infrared cutoff term @xmath79 $ ] in eq . and define @xmath80}=\int d\chi\ , e^{-s[\chi]-\delta s_k[\chi]+\int j \chi}.\ ] ] in fourier space , the cutoff term reads @xmath81=\int_q r_k(q ) \chi^*(q)\chi(q)\ ] ] and has the properties @xmath82 the effective average action is defined as a modified legendre transform of @xmath83 @xmath84=\left ( -w_k[j]+\int_x j\phi \right)_{j = j_\text{ex}}-\delta s_k[\phi].\ ] ] it has the important property , that it interpolates between the microscopic action @xmath85 $ ] and the full effective action @xmath71 $ ] @xmath86 & \rightarrow s[\phi ] & \quad ( k\rightarrow \infty),\\ \gamma_{k}[\phi ] & \rightarrow \gamma[\phi ] & \quad ( k\rightarrow 0).\end{aligned}\ ] ] in physical terms @xmath87 $ ] is the effective action in the presence of an infrared cutoff at a momentum scale @xmath88 .
only fluctuations with momenta larger than @xmath88 are included .
for example , a finite volume @xmath89 of the system under consideration would lead to an situation that is described by @xmath87 $ ] .
our method to determine @xmath87 $ ] ( and for @xmath90 also @xmath71 $ ] ) relies on the existence of an exact flow equation @xcite @xmath91 here the trace operation includes a momentum integration @xmath92 , as well as a sum over internal indices @xmath93 , according to the two real components in the decomposition @xmath94 . on the r.h.s . of
, @xmath95 is the second functional derivative of @xmath87 $ ] @xmath96)_{ij}(q , p)=\frac{\overset{\rightharpoonup}{\delta}}{\delta\phi_i(-q)}\gamma_k[\phi]\frac{\overset{\leftharpoonup}{\delta}}{\delta \phi_j(p)}.\ ] ] it is therefore a matrix in internal and momentum space .
correspondingly , @xmath97 in eq .
stands for @xmath98 the flow equation describes the evolution of the effective average action with the cutoff scale @xmath88 .    the functional differential equation is hard to be solved exactly . in principle , @xmath99 is described by infinitely many couplings .
we will use here an approximation with only a finite number of couplings .
this is achieved by an ansatz for a specific form of @xmath87 $ ] . in the absence of anomalies
the effective action @xmath71 $ ] is invariant under the same symmetries as the microscopic action @xmath85 $ ] .
this holds also for the effective average action @xmath87 $ ] , provided that the cutoff term @xmath100 $ ] is also invariant .
our ansatz will respect all symmetries of the classical action .
approximate solutions of the exact flow equations obtain from a truncation of the general form of the effective action .
we use here terms with up to two derivatives and truncate @xmath101 with @xmath102 .
this particular form is motivated by a more systematic derivative expansion and an analysis of symmetry constraints ( ward identities ) in appendix [ secttruncation ] .
we introduce the renormalized fields @xmath103 , @xmath104 , the renormalized kinetic coefficients @xmath105 , @xmath106 and we express the effective potential in terms of the renormalized invariant @xmath107 , with @xmath108 this yields @xmath109    for the effective potential , we use an expansion around the @xmath88-dependent minimum @xmath110 of the effective potential and the @xmath88-independent value of the chemical potential @xmath111 that corresponds to the physical particle number density @xmath1 .
we determine @xmath110 and @xmath111 by the requirements @xmath112 more explicitly we take a truncation for @xmath113 of the form @xmath114 in the symmetric phase we have @xmath115 , while in the phase with spontaneous symmetry breaking , we have @xmath116 . in summary , the flow of @xmath99 for fixed @xmath117 is described by four running renormalized couplings @xmath118 , @xmath119 , @xmath120 and @xmath37 .
in addition , we need the anomalous dimension @xmath121 . a computation of @xmath1 requires a flow equation of @xmath122 , which involves the couplings linear in @xmath123 , namely @xmath124 and @xmath125 .
the pressure is calculated by following the @xmath88-dependence of the height of the minimum @xmath126 .
all couplings @xmath118 , @xmath119 , @xmath120 , @xmath37 , @xmath127 , @xmath122 , @xmath128 , @xmath124 , @xmath125 depend on @xmath88 and @xmath15 .
the physical renormalized couplings obtain for @xmath90 .
they specify the thermodynamic potential @xmath129 as well as suitable derivatives of the potential and the correlation function .
the `` initial values '' at the scale @xmath130 are determined by the requirement @xmath131=s[\phi],\ ] ] using the microscopic action @xmath85 $ ] in eq . .
this implies the initial values @xmath132 we remain with the free microscopic couplings @xmath133 , @xmath134 , @xmath135 .
the coupling @xmath133 will be replaced by the scattering length @xmath0 in the next section .
we further choose units for @xmath29 where @xmath136 then our second free coupling is @xmath137 in consequence , besides the thermodynamic variables @xmath15 and @xmath111 our model is characterized by two free parameters , @xmath0 and @xmath5 .
often , we will concentrate on `` standard '' non - relativistic bosons with a linear @xmath29 derivative , such that @xmath138 .
the scattering length @xmath0 remains then the only free parameter . in the vacuum , where @xmath139 , this sets the relevant unit of length .
finally , we have for the infrared cutoff @xmath140 we choose the optimized cutoff function @xcite @xmath141 where we recall , that @xmath116 in the regime with spontaneous symmetry breaking .
it is convenient to work with real fields @xmath142 , @xmath94 , with fourier components @xmath143 the inverse propagator for the fields @xmath144 becomes a @xmath145 matrix in the space of @xmath146 and @xmath147 , given by the second functional derivative of @xmath99 . for a real constant background field @xmath148 ,
@xmath149 the latter becomes diagonal in momentum space @xmath150 for our truncation one has at @xmath117 @xmath151 here , primes denote derivatives with respect to @xmath107 ( not @xmath152 ) . in the phase with spontaneous symmetry breaking , the infrared cutoff in the flow equation
adds to the diagonal term in a piece @xmath153 .
this effectively replaces @xmath154 in eq .
whenever @xmath155 , thus providing for an efficient infrared regularization .
we project the flow equation of the effective average action onto equations for the coupling constants by using appropriate background fields and taking functional derivatives .
the flow equation for the effective potential obtains by using a space- and time - independent background field in eq .
, with @xmath156 @xmath157 the propagator @xmath158 is here determined from @xmath159 with @xmath160 again primes denote a differentiation with respect to @xmath107 .
we switch to renormalized fields by making a change of variables in the differential equation @xmath161 we can now derive the flow equations for the couplings @xmath110 and @xmath162 by appropriate differentiation of with respect to @xmath107 . the flow equation for @xmath163
is given more explicitly in appendix [ secfloweffectivepotential ] .
differentiation with respect to @xmath48 yields the flow of @xmath122 , @xmath124 , @xmath125 .
we use in detail @xmath164 where we recall , that @xmath165 in our truncation . to determine the flow equation of @xmath118
, we use the condition that @xmath166 for all @xmath88 , and therefore @xmath167 we show the flow of @xmath119 and @xmath118 in figs .
[ figflowoflambda ] , [ figflowofrho ] for @xmath168 , @xmath20 and different values of @xmath133 ( with @xmath138 ) .
the change in @xmath118 is rather modest
. this will be different for nonzero temperature .     with the scale parameter @xmath88 for different start values ,
corresponding to @xmath169 ( dotted ) , @xmath170 ( solid ) and @xmath171 ( dashed ) .
the first case is plotted for zero density ( @xmath172 ) only , while the last two cases are plotted also for unit density ( @xmath168 ) .
the curves @xmath172 and @xmath168 are identical within the plot resolution .
the solid and the dashed curve correspond to @xmath173 and @xmath174 , respectively . ]    .
the parameters for the solid and the dashed curves are the same as in figure [ figflowoflambda ] . ]
the flow of @xmath122 is given by @xmath175 and similar for the flow of @xmath124 and @xmath125 , @xmath176 where the last equation holds since @xmath177 in our truncation .
for a derivation of @xmath178 and the flow equations for @xmath120 and @xmath37 , we have to evaluate the flow equation for a background field depending on @xmath179 and @xmath180 .
we use an analytic continuation @xmath181 and obtain the flow equation for @xmath120 from @xmath182 with four - volume @xmath183 .
the projection prescription for @xmath37 is @xmath184 and similar for @xmath127 @xmath185 after the functional differentiation , we evaluate the expressions , , at homogeneous background fields .
these calculations are a little intricate , but standard and straightforward in principle .
more explicit flow equations are given in app .
[ sectflowofkineticcoefficients ] . eventually , it is always possible to perform the matsubara sums and also the spatial momentum integration analytically . in figure [ figflowkinetic ]
we show the flow of @xmath127 , @xmath120 and @xmath37 at zero temperature and for density @xmath168 .
( solid ) , @xmath120 ( dashed ) and @xmath37 for a scattering length @xmath173 , temperature @xmath20 and density @xmath168 . ]
the kinetic coefficient @xmath127 starts on the large scale with @xmath186 , increases a little around @xmath187 and saturates to a constant .
in contrast , the coefficient @xmath120 starts to decrease after a short period of increase with @xmath127 .
for very tiny scales @xmath88 , @xmath120 would finally go to zero .
the frequency dependence of the propagator is then governed by the quadratic frequency coefficient @xmath37 . in three spatial dimensions , however , this decrease of @xmath120 is so slow that it is not relevant on the length scales of experiments .
this is one of the reasons why bogoliubov theory , which neglects the appearance of @xmath37 , describes experiments with ultracold bosonic quantum gases in three dimensions with so much success .
the coefficient @xmath37 is always generated in the phase with spontaneous symmetry breaking @xcite . its @xmath88-dependence is also shown in figure [ figflowkinetic ] .
the vacuum is defined to have zero temperature @xmath20 and vanishing density @xmath172 , which also implies @xmath115 .
the interaction strength @xmath119 at the scale @xmath188 determines the four point vertex at zero momentum .
it is directly related to the scattering length @xmath0 for the scattering of two particles in vacuum , which is experimentally observable .
we therefore want to replace the microscopic coupling @xmath133 by the renormalized coupling @xmath0 . in our units ( @xmath33 ) , one has the relation @xmath189 the vacuum properties can be computed by taking for @xmath20 the limit @xmath190 .
we may also perform an equivalent and technically simple computation in the symmetric phase by choosing @xmath191 such , that @xmath192 .
this guarantees that the boson field @xmath193 is a gap - less propagating degree of freedom .
we first investigate the model with a linear @xmath29-derivative , @xmath194 , @xmath195 .
projecting the flow equation , we find the following equations : @xmath196 the propagator is not renormalized , @xmath197 , @xmath198 , @xmath199 , and one finds @xmath200 . the coupling @xmath125 is running according to @xmath201
since @xmath125 appears only in its own flow equation , it is of no further relevance in the vacuum .
also , no coupling @xmath37 is generated by the flow and we have therefore set @xmath35 on the r.h.s . of eqs . and .    inserting in eq .
the vacuum values @xmath116 and @xmath36 , we find @xmath202 the solution @xmath203 tends to a constant for @xmath204 , @xmath205 .
the dimensionless variable @xmath206 goes to zero , when @xmath88 goes to zero .
this shows the infrared freedom of the theory . for fixed ultraviolet cutoff
, the scattering length @xmath207 as a function of the initial value @xmath133 , has an asymptotic maximum @xmath208 the relation between @xmath0 and @xmath133 is shown in fig .
[ figscatteringbound ] .     in dependence on the microscopic interaction strength @xmath133 ( solid ) .
the asymptotic maximum @xmath209 is also shown ( dashed ) . ]    as a consequence of eq . , the nonrelativistic bosons in @xmath210 are a `` trivial theory '' in the sense that the bosons become noninteracting in the limit @xmath211 , where @xmath212 .
the upper bound has important practical consequences .
it tells us , that whenever the `` macrophysical length scales '' are substantially larger than the microscopic length @xmath8 , we deal with a weakly interacting theory . as an example , consider a boson gas with a typical inter - particle distance substantially larger than @xmath8 .
( for atom gases @xmath8may be associated with the range of the van der waals force . )
we may set the units in terms of the particle density @xmath1 , @xmath168 .
in these units @xmath213 is large , say @xmath214 . this implies a very weak interaction , @xmath215 . in other words
, the scattering length can not be much larger than the microscopic length @xmath8 .
for such systems , perturbation theory will be valid in many circumstances .
we will find that the bogoliubov theory indeed gives a reliable account of many properties .
even for an arbitrary large microphysical coupling @xmath216 , the renormalized physical scattering length @xmath0 remains finite .
let us mention , however , that the weak interaction strength does not guarantee the validity of perturbation theory in all circumstances .
for example , near the critical temperature of the phase transition between the superfluid and the normal state , the running of @xmath162 will be different from the vacuum . as a consequence ,
the coupling will vanish proportional to the inverse correlation length @xmath217 as @xmath15 approaches @xmath2 , @xmath218 .
indeed , the phase transition will be characterized by the non - perturbative critical exponents of the wilson - fisher fixed point . also for lower dimensional systems , the upper bound for @xmath219 is no longer valid - for example the running of @xmath119 is logarithmic for @xmath220 . for our models with @xmath221 ,
the upper bound becomes dependent on @xmath222 .
it increases for @xmath223 . in the limit @xmath224 , it is replaced by the well known `` triviality bound '' of the four dimensional relativistic model , which depends only logarithmically on @xmath213 . finally , for superfluid liquids , as @xmath225 , one has @xmath226 , such that for @xmath227 one finds a large concentration @xmath11 .
the situation for dilute bosons seems to contrast with ultracold fermion gases in the unitary limit of a feshbach resonance , where @xmath0 diverges .
one may also think about a feshbach resonance for bosonic atoms , where one would expect a large scattering length for a tuning close to resonance . in this case , however , the effective action does not remain local .
it is best described by the exchange of molecules .
the scale of nonlocality is then given by the gap for the molecules , @xmath228 .
only for momenta
@xmath229 the effective action becomes approximately local , such that @xmath230 for our approximation .
close to resonance , the effective cutoff is low and again in the vicinity of @xmath231 .
the density sets a crucial scale for our problem .
its precise determination is mandatory for quantitative precision .
we will discuss two different methods for its determination and show that the results agree within our precision .
for @xmath20 , we also find agreement with the ward identity @xmath232 .    the first method is to derive flow equations for the density .
this has the advantage that the occupation numbers for a given momentum @xmath233 are mainly sensitive to running couplings with @xmath234 . in the grand canonical formalism , the density is defined by @xmath235{\big |}_{\phi=\phi_0,\sigma=\sigma_0}\ ] ] we can formally define a @xmath88-dependent density @xmath122 by @xmath236{\big |}_{\phi=\phi_0,\sigma=\sigma_0}=-(\partial_\sigma u)(\rho_0,\sigma_0).\ ] ] the flow equation for @xmath122 is given in eq .
and the physical density obtains for @xmath188 . the term @xmath237 that enters eq . is the derivative of the flow equation for @xmath163 with respect to @xmath48 . to compute it
, we need the @xmath48-dependence of the propagator @xmath158 in the vicinity of @xmath111 . within a systematic derivative expansion ,
we use the expansion of @xmath113 and the kinetic coefficients @xmath238 and @xmath239 to linear order in @xmath240 , as described in appendix [ secttruncation ] . here , @xmath241 and @xmath242 are the coefficient functions of the terms linear in the @xmath29-derivative and linear in @xmath243 , respectively .
no reasonable qualitative behavior is found , if the linear dependence of @xmath238 and @xmath239 on @xmath240 is neglected .
also , the scale dependence of @xmath124 and @xmath125 are quite important .
the flow equations for @xmath124 and @xmath125 can be obtained directly by differentiating the flow equation of the effective potential with respect to @xmath48 and @xmath107 , cf .
the situation is more complicated for the kinetic coefficients @xmath244 and @xmath245 .
their flow equations have to be determined by taking the @xmath48-derivative of the flow equation for @xmath246 and @xmath242 . as discussed in app .
[ secttruncation ] , we use in this paper the approximation @xmath247 .    as a check of both our method and our numerics
, we also use another way to determine the particle density .
this second method is more robust with respect to shortcomings of the truncation , but less adequate for high precision calculations as needed e.g. to determine the condensate depletion .
the second method determines the pressure @xmath248 as a function of the chemical potential @xmath111 . here ,
the effective potential is normalized by @xmath249 at @xmath20 , @xmath172 .
the flow of the pressure can be read of directly from the flow equation of the effective potential and is independent of the couplings @xmath124 and @xmath125 .
we calculate the pressure as a function of @xmath48 and determine the density @xmath250 by taking the @xmath48-derivative numerically .
it turns out that @xmath251 is in very good approximation given by @xmath252 , where the constant @xmath11 can be determined from a numerical fit .
the density is thus linear in @xmath48 .
at zero temperature and for @xmath138 , we can additionally use the ward identities connected to galilean symmetry , which yield @xmath232 .
we compare our methods in figure [ densitycompared ] and find that they give numerically the same result .
we stress again the importance of a reliable method to determine the density , since we often rescale variables by powers of the density to obtain dimensionless variables .    .
we use three different methods : @xmath253 from the flow equation ( triangles ) , @xmath232 as implied by galilean symmetry ( stars ) and @xmath254 , where the pressure @xmath255 ( boxes ) was obtained from the flow equation and phenomenologically fitted by @xmath256 ( solid lines ) .
units are arbitrary and we use @xmath257 , @xmath214 . ]      we want to split the density into a condensate part @xmath258 and a density for uncondensed particles or `` depletion density '' @xmath259 . for our model
the condensate density is given by the `` bare '' order parameter @xmath260 in order to show this , we introduce occupation numbers @xmath261 for the modes with momentum @xmath233 with normalization @xmath262 one formally introduces a @xmath233 dependent chemical potential @xmath263 in the grand canonical partition function @xmath264}=\text{tr}e^{-\beta(h-\omega_3\int_{\vec{p}}\sigma(\vec{p})n(\vec{p}))},\ ] ] with three dimensional volume @xmath265
. then one can define the occupation numbers by @xmath266{\big |}_{\phi=\phi_0,\sigma(\vec{p})=\sigma_0}.\ ] ] this construction allows us to use @xmath88-dependent occupation numbers by the definition @xmath267{\big |}_{\phi=\phi_0,\sigma(\vec{p})=\sigma_0}.\ ] ] one can derive a flow equation for this occupation number @xmath268 @xcite : @xmath269{\big |}_{\phi_0,\sigma_0}(\partial_k\rho_0 ) .
\label{flowofnp}\end{aligned}\ ] ]    we split the density occupation number into a @xmath270-distribution like part and a depletion part , which is regular in the limit @xmath271 @xmath272 one can see from the flow equation for @xmath268 that the only contribution to @xmath273 comes from the second term in equation . within a more detailed analysis @xcite
one finds @xmath274 we therefore identify the condensate density with the bare order parameter @xmath275 correspondingly , we define the @xmath88-dependent quantities @xmath276 and compute @xmath277 by a solution of its flow equation .    even at zero temperature ,
the repulsive interaction connected with a positive scattering length @xmath0 causes a portion of the particle density to be outside the condensate . from dimensional reasons ,
it is clear , that @xmath278 should be a function of @xmath279 .
the prediction of bogoliubov theory or , equivalently , mean field theory , is @xmath280 .
we may determine the condensate depletion from the solution to the flow equation for the particle density , @xmath281 , and @xmath282 .    from galilean invariance for @xmath20 and @xmath138
, it follows that @xmath283 with @xmath284 .
this gives an independent determination of @xmath258 . in figure [ figdepletion ]
we plot the depletion density obtained from the flow of @xmath1 and @xmath285 over several orders of magnitude . apart from some numerical fluctuations for small @xmath279
, we find that our result is in full agreement with the bogoliubov prediction .     as a function of the dimensionless scattering length @xmath286 .
for the solid curve , we vary @xmath0 with fixed @xmath168 , for the dashed curve we vary the density at fixed @xmath174 . the dotted line is the bogoliubov - result @xmath287 for reference .
we find perfect agreement of the three determinations .
the fluctuations in the solid curve for small @xmath286 are due to numerical uncertainties .
their size demonstrates our numerical precision . ]      for @xmath20 a quantum phase transition separates the phases with @xmath115 and @xmath288 . in this section
, we investigate the phase diagram at zero temperature in the cube spanned by the dimensionless parameters @xmath289 , @xmath290 and @xmath291 .
this goes beyond the usual phase transition for nonrelativistic bosons , since we also include a microscopic second @xmath29-derivative @xmath292 , and therefore models with a generalized microscopic dispersion relation . for non - vanishing @xmath5 ( i.e. for a nonzero initial value of @xmath293 with @xmath294 in app .
[ secttruncation ] ) , the galilean invariance at zero temperature is broken explicitly . for large @xmath5
, we expect a crossover to the `` relativistic '' @xmath26 model .
if we send the initial value of the coefficient of the linear @xmath29-derivative @xmath295 to zero , we obtain the limiting case @xmath25 .
the symmetries of the model are now the same as those of the relativistic @xmath26 model in four dimensions .
the space - time - rotations or lorentz symmetry replace galilean symmetry .
it is interesting to study the crossover between the two cases .
since our cutoff explicitly breaks lorentz symmetry , we investigate in this paper only the regime @xmath296 .
detailed investigations of the flow equations for @xmath25 can be found in the literature @xcite .
the phase diagram in the @xmath297 plane with @xmath298 is shown in figure [ figqptsigmava1 ] .
the critical chemical potential first increases linearly with @xmath5 and then saturates to a constant .
the slope in the linear regime as well as the saturation value depend linearly on @xmath299 for @xmath300 .    at @xmath20 ,
the critical exponents are everywhere the mean field ones ( @xmath198 , @xmath301 ) .
this is expected : it is the case for @xmath138 @xcite , and for @xmath302 the theory is equivalent to a relativistic @xmath26 model in @xmath303 dimensions .
this is just the upper critical dimension of the wilson - fisher fixed point @xcite .
-@xmath5 plane for @xmath298 . ]    from section 3 we know that for @xmath138 the parameter @xmath299 is limited to @xmath304 . for @xmath138 and a small scattering length @xmath212
, a second order quantum phase transition divides the phases without spontaneous symmetry breaking for @xmath305 from the phase with a finite order parameter @xmath288 for @xmath306 .
it is an interesting question , whether this quantum phase transition at @xmath307 also occurs for larger scattering length @xmath0 .
we find in our truncation that this is indeed the case for a large range of @xmath0 , but not for @xmath308 . here
, the critical chemical potential suddenly increases to large positive values as shown in fig .
[ figqptsigmaofa ] . for @xmath309
this increase happens even earlier .
( for a truncation with @xmath310 , the phase transition would always occur at @xmath311 . )
we plot the @xmath312 plane of the phase diagram for different values of @xmath5 in figure [ figqptsigmaofa ] .
the form of the critical line can be understood by considering the limits @xmath313 as well as @xmath314 .
-@xmath299 plane for @xmath315 ( dotted ) , @xmath316 ( dashed ) and @xmath138 ( solid ) . ]    for a fixed chemical potential , the order parameter @xmath118 as a function of @xmath0 goes to zero at a critical value @xmath317 as shown in fig .
[ figqptrhoofa ] .
this happens in a continuous way and the phase transition is therefore of second order . for @xmath318 ,
we find @xmath319 .
we emphasize , however , that @xmath317 is of the order of the microscopic distance @xmath8 .
universality may not be realized for such values , and the true phase transition may depend on the microphysics .
for example , beyond a critical value for the repulsive interaction , the system may form a solid .
ultracold atom gases correspond to metastable states which may lose their relevance for @xmath320 .
for @xmath309 and @xmath321 the phase transition occurs for @xmath322 such that universal behavior is expected .    , with @xmath214 .
the density @xmath323 as a function of the scattering length @xmath0 goes to zero at a critical @xmath324 , indicating a second order quantum phase transition at that point . ]
so far , we have only discussed the vacuum and the dense system at zero temperature .
a non vanishing temperature @xmath15 will introduce an additional scale in our problem . for small @xmath325
we expect only small corrections .
however , as @xmath15 increases the superfluid order will be destroyed .
near the phase transition for @xmath326 and for the disordered phase for @xmath327 , the characteristic behavior of the boson gas will be very different from the @xmath328 limit .
for @xmath329 the particle density @xmath1 receives a contribution from a thermal gas of bosonic ( quasi- ) particles .
it is no longer uniquely determined by the superfluid density @xmath118 .
we may write @xmath330 and observe , that @xmath331 is enforced by galilean symmetry only for @xmath332 .
the heat bath singles out a reference frame , such that for @xmath329 galilean symmetry no longer holds . in our formalism
, the thermal contribution @xmath333 appears due to modifications of the flow equations for @xmath334 .
we start for high @xmath88 with the same initial values as for @xmath20 . as long as @xmath335
the flow equations receive only minor modifications . for @xmath336 or smaller , however , the discreteness of the matsubara sum has important effects .
we plot in fig .
[ fignoftemperature ] the density as a function of @xmath15 for fixed @xmath337 .
( solid ) and order parameter @xmath338 ( dashed ) as a function of the temperature @xmath339 .
the units are arbitray with @xmath340 and @xmath214 .
the plot covers only the superfluid phase . for higher temperatures ,
the density is given by the thermal contribution @xmath341 only . ]    in fig .
[ fignofsigma ] we show @xmath342 , similar to fig .
[ densitycompared ] , but now for different @xmath0 and @xmath15 . for @xmath20
the scattering length sets the only scale besides @xmath1 and @xmath48 , such that by dimensional arguments @xmath343 .
bogoliubov theory predicts @xmath344 the first term on the r.h.s . gives the contribution of the ground state , while the second term is added by fluctuation effects . for small scattering length @xmath0 ,
the ground state contribution dominates .
we have then @xmath345 for @xmath168 and @xmath346 can be treated as a small quantity . for @xmath347 and small @xmath0
one expects @xmath348 .
the curves in fig .
[ fignofsigma ] for @xmath349 show that the density , as a function of @xmath48 , is below the curve obtained at @xmath20 .
this is reasonable , since the statistical fluctuations now drive the order parameter @xmath118 to zero . at
very small @xmath48 , the flow enters the symmetric phase .
the density is always positive , but for simplicity , we show the density as a function of @xmath48 in figure [ fignofsigma ] only in those cases , where the flow remains in the phase with spontaneous @xmath46 symmetry breaking .     for different temperatures and scattering length .
we plot @xmath342 in arbitrary units , with @xmath214 , and for a scattering length @xmath350 ( solid and dotted ) , @xmath174 ( dashed and dashed - dotted ) . the temperature is @xmath20 ( solid and dashed ) and @xmath349 ( dotted and dashed - dotted ) . ]    for temperatures above the critical temperature , the order parameter @xmath118 vanishes at the macroscopic scale and so does the condensate density @xmath351 .
the density is now given by a thermal distribution of particles with nonzero momenta .
up to small corrections from the interaction @xmath352 , it is described by a free bose gas , @xmath353 with the `` bose function '' @xmath354    in figure [ figrhooftemperature ] we show the dimensionless order parameter @xmath355 as a function of the dimensionless temperature @xmath356 .     as a function of the dimensionless temperature @xmath357 for scattering length @xmath174 . here , we varied @xmath15 keeping @xmath48 fixed .
numerically , this is equivalent to varying @xmath48 with fixed @xmath15.[figrhooftemperature ] ]    this plot shows the second order phase transition from the phase with spontaneous @xmath46 symmetry breaking at small temperatures to the symmetric phase at higher temperatures .
the critical temperature @xmath2 is determined as the temperature , where the order parameter just vanishes - it is investigated in the next section .
since we find @xmath358 , the condensate fraction @xmath359 as a function of @xmath356 resembles the order parameter @xmath355 .
we plot @xmath127 as a function of @xmath356 in fig .
[ abaroftemperature ] . except for a narrow region around @xmath2 ,
the deviations from one remain indeed small . near @xmath2
the gradient coefficient @xmath127 diverges according to the anomalous dimension , @xmath360 , with @xmath361 the anomalous dimension .
the correlation length @xmath362 diverges with the critical exponent @xmath363 , @xmath364 , such that @xmath365 here , @xmath361 and @xmath363 are the critical exponents for the wilson fisher fixed point of the classical three - dimensional @xmath26 model , @xmath366 , @xmath367 @xcite .
, as a function of the dimensionless temperature @xmath357 , and for scattering length @xmath174 . here , we varied @xmath15 keeping @xmath48 fixed .
numerically , this is equivalent to varying @xmath48 with fixed @xmath15 .
the plot covers only the phase with spontaneous symmetry breaking . for higher temperatures , the symmetric phase has @xmath368 .
the divergence of @xmath127 for @xmath369 reflects the anomalous dimension @xmath361 of the wilson - fisher fixed point . ]    in figure [ figrhoofsigma ] we plot @xmath355 as a function of the chemical potential @xmath48 for different temperatures and scattering lengths .
we find , that @xmath370 is indeed approached in the limit @xmath328 , as required by galilean invariance .
all figures of this section are for @xmath138 .
the modifications for @xmath371 are mainly quantitative , not qualitative .
, as a function of the chemical potential .
we use arbitrary units with @xmath214 .
the curves are given for a scattering length @xmath350 ( solid and dotted ) , @xmath174 ( dashed and dashed - dotted ) and temperature @xmath372 ( solid and dashed ) and @xmath349 ( dotted and dashed - dotted ) .
at zero temperature , galilean invariance implies @xmath323 , which we find within our numerical resolution . ]
the critical temperature @xmath2 for the phase transition between the superfluid phase at low temperature and the disordered or symmetric phase at high temperature depends on the scattering length @xmath0 . by dimensional reasoning
, the temperature shift @xmath373 obeys @xmath374 .
the proportionality coefficient can not be computed in perturbation theory @xcite .
it depends on @xmath5 and we concentrate here on @xmath138 .
monte - carlo simulations in the high temperature limit , where only the lowest matsubara frequency is included , yield @xmath375 @xcite . within the same setting , renormalization group studies @xcite yield a similar result , for composite bosons see @xcite . near @xmath2 , the long wavelength modes with momenta @xmath376 dominate the `` long distance quantities '' .
then a description in terms of a classical three dimensional system becomes valid .
this `` dimensional reduction '' is achieved by `` integrating out '' the nonzero matsubara frequencies .
however , both @xmath377 and @xmath1 are dominated by modes with momenta @xmath378 such that corrections to the classical result may be expected .
we have computed @xmath2 numerically by monitoring the zero of @xmath118 , as shown in fig .
[ figrhooftemperature ] , @xmath379 .
our result is plotted in fig .
[ tcofa ] . in the limit @xmath212
we find for the dimensionless critical temperature @xmath380 , which is in good agreement with the expected result for the free theory @xmath381 .
for the shift in @xmath2 due to the finite interaction strength , we obtain @xmath382 we expect that the result for @xmath383 depends on the truncation and may change somewhat if additional higher order couplings are included .     as a function of the dimensionless scattering length @xmath279 ( points ) .
we also plot the linear fit @xmath384 ( solid line ) . ]
the macroscopic sound velocity @xmath385 is a crucial quantity for the hydrodynamics of the gas or liquid .
it is accessible to experiment . as a thermodynamic observable ,
the adiabatic sound velocity is defined as @xmath386 where @xmath18 is the particle mass ( in our units @xmath387 ) , @xmath251 is the pressure , @xmath1 is the particle density , and @xmath388 is the entropy per particle .
it is related to the isothermal sound velocity @xmath389 by @xmath390 where we use our units @xmath33 .
one needs the `` equation of state '' @xmath391 and @xmath392 by dimensional analysis , one has @xmath393 with @xmath394 ( in bogoliubov theory ) , and @xmath395 ( free theory ) , such that for small @xmath11 @xmath396    at zero temperature the second term in eq .
vanishes , such that @xmath397 . for the isothermal sound velocity one
has @xmath398 we can now use @xmath399 and infer @xmath400    one may also define a microscopic sound velocity @xmath401 , which characterizes the propagation of ( quasi- ) particles . at zero temperature , where we can perform the analytic continuation to real time , we can calculate the microscopic sound velocity from the dispersion relation @xmath402 ( with @xmath403 ) . in turn
, the dispersion relation is obtained from the effective action by setting @xmath404 , where @xmath405 is the full inverse propagator .
we perform the calculation explicitly in app .
[ sectpropdisp ] and find @xmath406    the bogoliubov result for the sound velocity is in our units @xmath407 in three dimensions , the decrease of @xmath120 is very slow and the coupling @xmath37 remains comparatively small even on macroscopic scales , cf .
[ figflowkinetic ] .
we thus do not expect measurable deviations from the bogoliubov result for the sound velocity at @xmath20 . in figure [ figsound ] , we plot our result over several orders of magnitude of the dimensionless scattering length and , indeed , find no deviations from bogoliubovs result .     at zero temperature , as a function of the scattering length @xmath279 . within the plot resolution the curves obtained by varying @xmath0 with fixed @xmath1 , by varying @xmath1 with fixed @xmath0 , and the bogoliubov result , @xmath408 , coincide . ]
we finally show that for @xmath20 the macroscopic and microscopic sound velocities are equal , @xmath409 . for this purpose
, we use @xmath410 from the minimum condition @xmath411 , it follows @xmath412 combining this with the ward identities from app . [ secttruncation ] , namely @xmath413 and @xmath414 , valid at @xmath20 , it follows that the macroscopic sound velocity equals the microscopic sound velocity @xmath415
the use of functional renormalization yields a quantitative description of the whole phase diagram for dilute non - relativistic bosons in three dimensions .
this can describe a gas of ultracold bosonic atoms .
more generally , our results can also be applied to quantum phase transitions or low temperature phase transitions , whenever the most relevant excitations correspond to bosonic quasi - particles . for this reason we deal with a general dispersion relation , involving in the classical propagator terms linear and quadratic in the frequency . as a function of temperature @xmath15 and effective chemical potential @xmath48 ,
we have computed the pressure @xmath251 , the density @xmath1 , condensate @xmath285 , superfluid density @xmath118 , and the sound velocity @xmath401 , in dependence on two system parameters , namely the scattering length @xmath0 and a dimensionless coupling @xmath5 parameterizing the classical dispersion relation . for @xmath20 and @xmath138
we find very good agreement with the bogoliubov theory .
as @xmath15 increases , the condensate melts at a critical temperature @xmath2 , which exceeds the one for the free theory by @xmath416 .
we find a second order phase transition in the universality class of the three - dimensional @xmath26 model , with the associated universal critical exponents .
we have found an upper bound for the scattering length , which is of the order of the microscopic length scale @xmath8 .
this indicates the breakdown of the pointlike approximation . for atom gases
, @xmath8 typically corresponds to the range of the van der waals interaction .
the upper bound is a result of strong fluctuation effects in presence of a large microscopic interaction .
even for an arbitrary strong microscopic interaction , the quantum fluctuations reduce the renormalized coupling , which corresponds to the physical scattering length in vacuum . for non - relativistic atoms with pointlike interactions
, a strong interaction regime can therefore only be realized at high density , @xmath417 .
superfluid liquids , like @xmath225 are good examples for such systems .
first functional renormalization studies for such systems have already been performed @xcite .
our extended treatment overcomes several problems of this early approach and we will apply it to @xmath225 in the future . except for the vicinity of the phase transition , however , the detailed microphysics may play an important role , since all relevant length scales are of the order of the microphysical length scale @xmath8 .
other interesting extensions concern lower dimensional systems .
we look forward to a unified functional renormalization description of bosons with pointlike interactions , for arbitrary dimensions and arbitrary interaction strength .
we use a derivative expansion for the truncation of the effective average action with derivative operators up to four momentum dimensions @xmath418 here , we employ the renormalized fields @xmath419 and coupling functions @xmath163 , @xmath420 , @xmath421 .
we fix the wave function renormalization factor @xmath127 such that @xmath422 .
terms of the form @xmath423 or @xmath424 are not included here , since they are expected to play a sub - leading role . for a systematic derivative expansion they have to be added - the terms with up to two derivatives
can be found in @xcite . in terms of dimensions ,
the operator @xmath425 counts as two space derivatives for the nonrelativistic model with @xmath35 , while it counts as one space dimension for the relativistic model with @xmath426 .
we expand the @xmath88-dependent functions @xmath113 , @xmath241 , @xmath242 , @xmath427,@xmath428 and @xmath429 around the @xmath88-dependent minimum @xmath110 of the effective potential and the @xmath88-independent value of the chemical potential @xmath111 that corresponds to the physical particle number density @xmath1 .
for example , with @xmath430 , one has @xmath431    let us concentrate on the non - relativistic model where @xmath36 , @xmath35 in the microscopic action .
at zero temperature , we can perform an analytic continuation to real time @xmath432 .
the microscopic action is then @xmath433&=&-\int_{-\infty}^{\infty } dt\int d^3x \\ & & { \big \{}\phi^*\,(-i\partial_t-\sigma-\delta)\,\phi\,+\,\frac{1}{2}\lambda(\phi^*\phi){\big
\}}. \label{realtimemicroscopicaction}\end{aligned}\ ] ] in addition to space translations , rotations and time translations , two further symmetries constrain the possible forms of the couplings in @xmath434 . in order to derive these constraints
, we extend eq . to a @xmath435-dependent source @xmath436 .
first , there is a semi - local @xmath46 symmetry of the form @xmath437 this holds since the combination @xmath438 acts as a covariant derivative .
in addition , we have the invariance under galilean boost transformations of the fields @xmath439 while the invariance of the interaction term under this symmetry is obvious , its realization for the kinetic term is more involved . performing the transformation explicitly , one finds @xmath440 such that indeed the combination @xmath441 leads to this invariance . on the other hand , the validity of the galilean symmetry for an effective action guarantees that only the combination or powers of this operator act on @xmath193 .
an operator of the form @xmath442 with @xmath443 would break the symmetry .
( note , that @xmath444 is also invariant . )    both the semi - local @xmath46 symmetry and the galilean symmetry are helpful only at zero temperature . at nonzero temperature ,
the analytic continuation to real time is no longer useful .
an analog version of the semi - local @xmath46 transformation for euclidean time @xmath29 would involve the imaginary part of the chemical potential @xmath48 , which has no physical meaning .
the dependence of physical quantities on @xmath445 is not restricted .
in addition , the galilean symmetry is broken explicitly by the thermal heat bath .    combining semi - local @xmath46 symmetry and galilean symmetry at @xmath20 , we find that the derivative operators @xmath446 , @xmath243 and the chemical potential term @xmath240 are combined to powers of the operator @xmath447 in addition to powers of that operator acting on @xmath193 , only spatial derivatives of terms , that are invariant under @xmath46 transformations , like @xmath448 , may appear . since the symmetry transformations act linearly on the fields , the full effective action @xmath71 $ ] is also invariant .
this also holds for the average action @xmath87 $ ] , provided that the cutoff term @xmath100 $ ] is invariant .
we can write the effective action as an expansion in the operator @xmath449 @xmath450&=&\int_x \bigg{\{}u_0(\rho)\\ \nonumber & & + \frac{1}{2}\tilde{z}(\rho)\left(\phi^*(-i\partial_t-(\sigma-\sigma_0)-\delta)\phi+c.c\right)\\ \nonumber & & + \frac{1}{2}\tilde{v}(\rho)\left(\phi^*(-i\partial_t-(\sigma-\sigma_0)-\delta)^2\phi+c.c\right)\\ & & + \dots\bigg{\ } } \label{eqeffectiveactiongalilean}\end{aligned}\ ] ] performing the wick rotation back to euclidean time , we can compare this to eq .
, and find for @xmath20 the relations @xmath451 and therefore @xmath452    we next compute the inverse propagator in a constant background field by expanding @xmath99 to second order in the fluctuations around this background . for this purpose , it is convenient to decompose @xmath453 the constant condensate field @xmath51 can be chosen to be real without loss of generality .
the fluctuating real fields are the radial mode @xmath454 and the goldstone mode @xmath455 , and @xmath456 .
the truncation of the effective average action reads in that basis    @xmath457 = \int_x{\bigg \{}u(\rho,\sigma ) % \\
% \nonumber % & & & + & \frac{1}{2}z_1(\rho,\sigma ) \left(i\sqrt{2}\phi_0\partial_\tau\phi_2+i\phi_1\partial_\tau
\phi_2-i\phi_2\partial_\tau\phi_1\right)\\ \nonumber & + & \frac{1}{2}z_2(\rho,\sigma ) \left(\sqrt{2}\phi_0(-\delta)\phi_1+\phi_1(-\delta ) \phi_1+\phi_2(-\delta)\phi_2\right)\\ \nonumber & + & \frac{1}{2}v_1(\rho,\sigma )
\left(\sqrt{2}\phi_0(-\partial_\tau^2)\phi_1+\phi_1(-\partial_\tau^2 ) \phi_1+\phi_2(-\partial_\tau^2)\phi_2\right)\\ \nonumber & + & v_2(\rho,\sigma ) \left(i\sqrt{2}\phi_0(\partial_\tau\delta)\phi_2+i\phi_1(\partial_\tau\delta ) \phi_2-i\phi_2(\partial_\tau\delta)\phi_1\right)\\ & + & \frac{1}{2}v_3(\rho,\sigma ) \left(\sqrt{2}\phi_0(-\delta^2)\phi_1+\phi_1(-\delta^2 ) \phi_1+\phi_2(-\delta^2)\phi_2\right){\bigg \}}.\end{aligned}\ ] ]    the inverse propagator @xmath95 can be inferred from an expansion in second order in @xmath454 and @xmath455 .
we keep the linear order in @xmath123 , which will be needed for the flow equation for the density .
this yields @xmath457 & = & \int_x{\bigg \{}u(\rho_0,\sigma_0)+u^{(\sigma)}\,(\sigma-\sigma_0)+\frac{1}{2}(u^\prime+2\rho_0 u^{\prime\prime})\phi_1 ^ 2+\frac{1}{2}u^\prime \phi_2 ^ 2\\ \nonumber & + & \frac{1}{2}\left(z_1+z_1^\prime \rho_0+z_1^{(\sigma)}(\sigma-\sigma_0)\right ) \left(i\phi_1\partial_\tau \phi_2-i\phi_2\partial_\tau\phi_1\right)\\ \nonumber & + & \frac{1}{2}\left(1 + 2 z_2^\prime \rho_0+z_2^{(\sigma)}(\sigma-\sigma_0)\right ) \left(\phi_1(-\delta ) \phi_1\right)\\ \nonumber & + & \frac{1}{2}\left(1+z_2^{(\sigma)}(\sigma-\sigma_0)\right ) \left(\phi_2(-\delta ) \phi_2\right)\\ \nonumber & + & \frac{1}{2}\left(v_1 + 2 v_1^\prime\rho_0+v_1^{(\sigma)}(\sigma-\sigma_0)\right ) \left(\phi_1(-\partial_\tau^2 ) \phi_1\right)\\ \nonumber & + & \frac{1}{2}\left(v_1+v_1^{(\sigma)}(\sigma-\sigma_0)\right ) \left(\phi_2(-\partial_\tau^2 ) \phi_2\right)\\ \nonumber & + & \left(v_2+v_2^\prime\rho_0+v_2^{(\sigma)}(\sigma-\sigma_0)\right ) \left(i\phi_1(\partial_\tau\delta ) \phi_2-i\phi_2(\partial_\tau\delta)\phi_1\right)\\ \nonumber & + & \frac{1}{2}\left(v_3 + 2v_3^\prime \rho_0+v_3^{(\sigma)}(\sigma-\sigma_0)\right ) \left(\phi_1(-\delta^2 ) \phi_1\right)\\ & + & \frac{1}{2}\left(v_3+v_3^{(\sigma)}(\sigma-\sigma_0)\right)\left(\phi_2(-\delta^2 ) \phi_2\right){\bigg \}},\end{aligned}\ ] ]    where we dropped the argument @xmath458 at several places and used the implicit rescaling condition @xmath459 . in our simple truncation , we take at @xmath117 only @xmath460 into account .
we neglect the contribution of the other couplings , i.e. set @xmath461 .
as shown above , it follows from the symmetry requirements at zero temperature , that @xmath462 , @xmath463 and @xmath464 .
the truncation @xmath294 therefore violates the galilean symmetry , as does our choice of the cutoff term @xmath465 . within our approximation ,
it is consistent to set @xmath466 at zero temperature .
also the deviations from this relation at finite temperature are neglected for simplicity in this paper .
this yields eq . .
in the following we discuss the role of continuous symmetries of the microscopic action @xmath85 $ ] .
since all these symmetries are linear in the fields , the full effective action @xmath71 $ ] is also symmetric . in the case
that the cutoff term @xmath100 $ ] is chosen invariant under the symmetry transformation in question , this also holds for the average action @xmath467 $ ] for finite @xmath88 . from noether s theorem
it follows that there exists a conserved current @xmath468 connected with every such symmetry .
if the action is formulated as an integral over the imaginary time @xmath29 the conservation equation implies for the current @xmath469 at zero temperature , we can perform a wick rotation to real time , @xmath470 , and eq .
takes the usual form @xmath471 the noether charge @xmath472 is conserved in time , i.e. @xmath473 .
this holds if @xmath474 falls off sufficiently fast at spatial infinity . at finite temperature
however , the situation is different . a simple analytic continuation to real time is no longer possible , since the configuration space is now a torus with periodicity @xmath475 in the @xmath29-direction
instead , we can integrate eq .
over complex time @xmath29 , giving @xmath476 from the symmetry , it now follows that there exists a solenoidal vector field or three component current @xmath477 .    a global symmetry of an action @xmath71 $ ] ( where @xmath434 could be replaced by @xmath120 or @xmath99 if appropriate ) can be formulated in its infinitesimal form as @xmath478=\gamma[\phi],\ ] ] with @xmath479 independent of @xmath480 .
here @xmath388 is the infinitesimal generator of the symmetry transformation . for a local transformation , where @xmath479 depends on @xmath480 , @xmath481 , we can expand @xmath482=\gamma[\phi]+\int_x \left\ { ( \partial_\mu\epsilon){\cal j}^\mu+(\partial_\mu\partial_\nu\epsilon){\cal k}^{\mu\nu}+\dots\right\}. \label{noetherexpansion}\ ] ]
the global symmetry implies that the expansion on the r.h.s .
of eq . starts with @xmath483 . here and in the following it
is implied that @xmath479 as well as its derivatives are infinitesimal , i.e. we keep only terms that are linear in @xmath479 .
the index @xmath484 goes over @xmath485 , representing @xmath486 in the real time case and @xmath487 for imaginary time .
. implies for arbitrary @xmath488 @xmath489}{\delta\phi}\epsilon s\phi-(\partial_\mu\epsilon){\cal j}^\mu-(\partial_\mu\partial_\nu\epsilon){\cal k}^{\mu\nu}+\dots\right\}=0 . \label{eqnoetherintegral}\ ] ] our notation is for real fields and implies a summation over components , if appropriate . in a complex basis one replaces @xmath490 by @xmath491 .
is valid for all field configurations @xmath193 and not only for those that fulfill the field equation @xmath492}{\delta \phi}=0 $ ] . in consequence ,
the integrand is a total derivative @xmath493}{\delta\phi}\epsilon s\phi-(\partial_\mu\epsilon){\cal j}^\mu-(\partial_\mu\partial_\nu\epsilon){\cal k}^{\mu\nu}+\ldots & \\ & = -\partial_\mu\left(j^\mu \epsilon+\kappa^{\mu\nu}\partial_\nu\epsilon+\ldots\right ) . &
\end{aligned}\ ] ] we can now specialize to @xmath494 and find @xmath495}{\delta\phi } s\phi(x)=-\partial_\mu j^\mu.\ ] ] this defines the noether current @xmath496 . for solutions of the field equation , @xmath492}{\delta\phi}=0 $ ] ,
the current @xmath496 is conserved , @xmath497 .    for a given @xmath480
we can also specialize to @xmath498 which leads to @xmath499 this process can be continued , leading us to a whole tower of identities for the conserved current @xmath496 .
if the action @xmath71 $ ] includes derivatives only up to a finite order @xmath1 , i.e. can be written in the form @xmath72=\int_x{\cal l}(\phi,\partial\phi , \partial\partial\phi , \ldots , \partial^{(n)}\phi),\ ] ] the expansion on the right hand side of only contains terms up to order @xmath500 such that the tower of equations for @xmath496 can be solved . moreover , for homogeneous situations where @xmath492}{\delta\phi}$ ] is solved by a constant @xmath193 , the second term on the r.h.s .
vanishes since it includes a derivative .
we have then @xmath501 .
a convenient way to find the local currents employs parameters @xmath502 that decay sufficiently fast at infinity such that we can partially integrate eq . @xmath503}{\delta\phi } s\phi+\partial_\mu{\cal j}^\mu-\partial_\mu\partial_\nu{\cal k}^{\mu\nu}+\dots\right\}=0.\ ] ] this yields the local identity @xmath78}{\delta \phi}s\phi=-\partial_\mu{\cal j}^\mu+\partial_\mu\partial_\nu { \cal k}^{\mu\nu}-\dots.\ ] ] an expansion of the l.h.s . in derivatives often yields substantial information on @xmath504 etc . by inspection .
our construction yields a unique conserved local current @xmath496 for every generator of a continuous symmetry .
we note , however , that @xmath505 is also a conserved local current if @xmath124 and @xmath506 are independent of @xmath480 .
this remark is important if we want to associate @xmath496 with the current for a physical quantity .
a rotation invariant setting implies @xmath507 , but @xmath508 and @xmath124 may differ from zero .
after these general considerations we now specialize to nonrelativistic real time actions of the form @xmath72=\int_{-\infty}^{\infty}dt \int d^3x{\cal l}(\phi , ( i\partial_t+\delta)\phi,(i\partial_t+\delta)^2\phi,\dots).\ ] ] we assume , that @xmath434 invariant under the same symmetries as the action . from the symmetry under time translations
@xmath509 we find a conserved current @xmath510 .
up to a possible additive constant its @xmath435-component is the energy density , while the spatial components describe energy flux density .
the multiplicative constant @xmath124 gets fixed if we choose the units to measure energy .
the choice @xmath511 corresponds to @xmath512 .
similarly , the invariance under spatial translations @xmath513 implies a conserved current @xmath514 for each spatial direction @xmath515 . up to an additive constant @xmath516
the @xmath435-component is the conserved momentum density , @xmath517 , while the spatial components can be interpreted as a momentum flux density , with the diagonal components @xmath518 describing pressure .    from the global @xmath46 symmetry
@xmath519 we can infer the conservation of the current @xmath520 associated to the conserved particle number . in order to identify the total particle number with the charge of this current , @xmath521 , we need to fix a possible multiplicative constant @xmath124 . for this purpose
, we use the galilean boost invariance , described already in eq . .
it reads in its infinitesimal form @xmath522 and the conserved charge of @xmath523 is the center of mass , again up to an additive constant .
the generator @xmath523 can be decomposed as @xmath524 this implies for the current @xmath525 specializing to the @xmath435-component , identifying the momentum density @xmath526 and reintroducing the particle mass @xmath33 we find @xmath527 from this we can conclude that up to an additive constant @xmath528 is the particle density @xmath529 .    for the effective action we find for @xmath117 and constant @xmath530 the current @xmath531 using the normalization condition @xmath532
, this gives @xmath533 . at zero temperature , this is the particle density and the additive constant @xmath534 vanishes .
at nonzero temperature we can compare to eq . and find @xmath535 .
for completeness we also mention the symmetry under spatial rotations @xmath536 with orthogonal matrix @xmath537 , generators @xmath538 , and @xmath539 the antisymmetric tensor in three dimensions .
the infinitesimal transformation reads @xmath540 the time component of the conserved current @xmath541 is , of course , the angular momentum density .
we derive the flow equation for the effective potential by evaluating the flow equation for the average action for constant fields . inserting a real constant field @xmath542
one finds for @xmath543 the flow at fixed @xmath107 @xmath544 here , @xmath545 is the number of spatial dimensions and we use the abbreviations @xmath546 we dropped the arguments @xmath547 at several places on the right hand side .
primes denote derivatives with respect to @xmath107 . in the phase with spontaneous symmetry breaking , we have @xmath548
.    the matsubara sums over @xmath1 can be carried out by virtue of the formulae @xmath549 with @xmath550 .
this brings us to @xmath551 where we introduced @xmath552 in our simple truncation with @xmath553 , @xmath554 , @xmath461 , and at @xmath117 , the integrand in eq .
becomes mostly independent of the spatial momentum .
the integral can than be carried out and we find @xmath555 with @xmath556 that the momentum integral can be performed analytically is a nice feature of the cutoff .
the limit @xmath328 is obtained by substituting the @xmath557 functions with unity .
the flow of the effective potential contains a subtlety that can be seen in the limit @xmath558 ( @xmath515 ) , where we find @xmath559 the term @xmath560 in the last line is not present if @xmath293 is set to zero from the outset .
if @xmath238 is independent of @xmath107 , this term is independent of @xmath107 and gives only an overall shift of the effective potential .
we show in this appendix our results for the flow equation of the kinetic coefficients @xmath120 , @xmath127 and @xmath37
. we neglect all contributions from momentum dependent vertices . in other words
, we use @xmath107-independent constants @xmath561 , @xmath562 and @xmath554 . in our truncation with @xmath563 , and with the cutoff , we can perform all momentum integrations analytically , leading us to    @xmath564    here , @xmath545 is the number of spatial dimensions , @xmath565 and @xmath566 are as in .
the matsubara sums over @xmath1 can be performed analytically again by using and derivatives thereof .    in the limit @xmath328 ,
the matsubara frequencies are continuous @xmath567 and the sum becomes an integral @xmath568 .
the expressions for @xmath361 and @xmath569 in eq .
agree with those derived in @xcite , while our result for @xmath570 corrects an error in a first version of @xcite ( in equation ( c.8 ) ) .
the inverse propagator is given by the second functional derivative of the effective action @xmath571 and we find from the truncation @xmath572 here we use the abbreviations @xmath573 at zero temperature , we can analytically continue to real time @xmath574 , and find @xmath575    the dispersion relation is found from the on shell condition @xmath576 which yields @xmath577 the solutions for @xmath22 define the dispersion relation .
we find two branches , according to @xmath578 in the phase with spontaneous symmetry breaking , the @xmath579 branch of this solution is an `` optical mode '' , while the @xmath580 branch is a sound mode .
the microscopic sound velocity is @xmath581 . using @xmath582 , @xmath583 , @xmath584 and @xmath585 , we find @xmath586 the `` optical mode '' has at vanishing spatial momentum the frequency @xmath587 which diverges @xmath588 in the limit @xmath589 .    12 n.  n.  bogoliubov , j.  phys .
( moscow ) * 11 * , 23 ( 1947 ) .
s.  t.  beliaev , sov .
jetp * 7 * , 289 ( 1958 ) ; * 7 * , 299 ( 1958 ) ; j.  gavoret , p.  nozires , ann .
( n. y. ) * 28 * , 349 ( 1964 ) ; a.  a.  nepomnyashchii , y.  a.  nepomnyashchii , jetp lett .
* 21 * , 1 ( 1975 ) .
c.  castellani , c.  di castro , f.  pistolesi , g.  c.  strinati , phys .
lett .  * 78 * , 1612 ( 1997 ) ; f.  pistolesi , c.  castellani , c.  d.  castro , g.  c.  strinati , phys .
b * 69 * , 024513 ( 2004 ) . c.  wetterich , phys .
b * 77 * , 064504 ( 2008 ) .
k.  g.  wilson , phys .
b * 4 * , 3174 ( 1971 ) ; k.  g.  wilson , j.  b.  kogut , phys .
* 12 * , 75 ( 1974 ) ; f.  wegner , a.  houghton , phys .
rev .  a * 8 * , 401 ( 1973 ) ; j.  polchinski , nucl .
b * 231 * , 269 ( 1984 ) . c.  wetterich , z.  phys .
c * 48 * , 693 ( 1990 ) ; nucl .
b * 352 * , 529 ( 1991 ) . c.  wetterich , phys .
b * 301 * , 90 ( 1993 ) .
j.  berges , n.  tetradis and c.  wetterich , phys .
rept .   * 363 * , 223 ( 2002 ) ; c.  bagnuls , c.  bervillier , phys .  rept .  *
348 * , 91 ( 2001 ) ; b.  delamotte , arxiv : cond - mat/0702365 . n.  dupuis , k.  sengupta , europhys .
lett .  * 80 * , 50007 ( 2007 ) . j.  o.  andersen and m.  strickland , phys .  rev .
a * 60 * , 1442 ( 1999 ) ; l.  canet , b.  delamotte , d.  mouhanna and j.  vidal , phys .
b * 68 * , 064421 ( 2003 ) ; phys .
d * 67 * , 065004 ( 2003 ) ; c.  bervillier , a.  juttner and d.  f.  litim , nucl .  phys .
b * 783 * , 213 ( 2007 ) .
d.  t.  son and m.  wingate , annals phys .   *
321 * , 197 ( 2006 ) .
d.  f.  litim , phys .
b * 486 * , 92 ( 2000 ) ; j.  m.  pawlowski , annals phys .   *
322 * , 2831 ( 2007 ) . c.  wetterich , arxiv:0712.1926 [ cond-mat.stat-mech ] . t.  papenbrock and c.  wetterich , z.  phys .
c * 65 * , 519 ( 1995 ) .
d.  i.  uzunov , phys .
a * 87 * , 11 ( 1981 ) ; s.  sachdev , quantum phase transitions , cambridge university press , cambridge ( 1999 ) . k.  g.  wilson and m.  e.  fisher , phys .  rev .
* 28 * ( 1972 ) 240 .
a.  pelissetto and e.  vicari , phys .
rept .   * 368 * , 549 ( 2002 ) .
j.  o.  andersen , rev .
phys .   * 76 * , 599 ( 2004 ) .
p.  arnold , g.  moore , phys .
* 87 * , 120401 ( 2001 ) ; v.  a.  kashurnikov , n.  v.  prokofev , b.  v.  svistunov , phys .
lett .  * 87 * , 120402 ( 2001 ) .
g.  baym , j.  -p .
blaizot , m.  holzmann , f.  lalo , d.  vautherin , phys .
lett . * 83 * , 1703 ( 1999 ) ; j.  -p .
blaizot , r.  mndez - galain , n.  wschebor , phys .
e * 74 * , 051116 ( 2006 ) ; j.  -p .
blaizot , r.  mndez - galain , n.  wschebor , phys .
e * 74 * , 051117 ( 2006 ) .
s.  ledowski , n.  hasselmann , p.  kopietz , phys .
a * 69 * , 061601(r ) ( 2004 ) ; n.  hasselmann , s.  ledowski , p.  kopietz , phys .
a * 70 * , 063621 ( 2004 ) .
s.  diehl , h.  gies , j.  m.  pawlowski and c.  wetterich , phys .
a * 76 * , 021602(r ) ( 2007 ) .
t.  gollisch and c.  wetterich , phys .
b * 65 * , 134506 ( 2002 ) .