liquid state theory experienced remarkable successes in the understanding of the properties of dense , simple fluids @xcite and more recently , of complex liquids , like colloidal suspensions or globular polymers @xcite .
however a fully quantitative description of phase transitions can not be obtained by standard integral equations or diagrammatic expansions and requires specially devised approaches .
the freezing transition , often dominated by geometrical constraints , is usually understood on the basis of density functional theory @xcite , while fluctuation induced phenomena , like the liquid - vapour transition , need the completely different theoretical framework provided by the renormalization group method @xcite .
the hierarchical reference theory of fluids ( hrt ) was formulated more than 25 years ago as an attempt to reconcile the renormalization group ( rg ) approach , an extremely general and successful theory of second order phase transitions formulated in field - theoretical language , and the microscopic point of view adopted within liquid state theory @xcite .
the main goal of hrt was to develop a unified theory able to reproduce the universal properties at the critical point , as obtained by rg techniques , together with the non universal features of realistic models of fluids , present both in the critical region ( e.g. critical density and temperature , short range correlations ) as well as far from the critical point .
a comprehensive presentation of the hrt formalism , including a detailed derivation of the basic equations and an extensive discussion of its critical properties , can be found in ref .
the earliest applications of hrt to simple fluids and to the ising model are briefly discussed in the same reference .
however , most of the more interesting results and extensions were developed later . by now
it is therefore appropriate to summarize the advances in hrt which took place in the last decade , providing a guidance to the most significant recent applications .
such a critical discussion of the hrt formalism also gives the opportunity to emphasize its relationship with other theoretical tools which have been independently developed in different frameworks .
since the early nineties , sophisticated functional , non perturbative rg techniques ( nprg ) were devised within field theory , leading to a boost of applications to many and diverse physical situations : from particle physics , to quantum many body theory , to classical statistical mechanics @xcite .
only recently it was recognized that , despite the different language and physical context , a deep relation ties hrt and nprg @xcite , whose developments have proceeded quite independently for 20 years in distinct scientific communities .
to introduce the hierarchical reference theory of fluids it is convenient to refer to a system of @xmath0 classical particles confined in a volume @xmath1 , in spatial dimension @xmath2 , interacting via a pairwise potential @xmath3 .
generalizations to classical spin systems , quantum fluids and magnets do not present additional difficulties and will be examined in the next sections . the grand canonical partition function @xmath4 of the model
is written as : @xmath5= \sum_n \frac{z^{n}}{n!}\int { \rm d } { \bf r_1} ... {\rm d } { \bf r_n } \ , e^ { -\beta\sum_{i < j } v({\bf r}_i - { \bf r}_j)-\beta\sum_i u({\bf r}_i ) } \label{gran}\ ] ] where the fugacity @xmath6 is related to the chemical potential @xmath7 and to the thermal wavelength @xmath8 by @xmath9 with @xmath10 .
an external field @xmath11 is also included for convenience and @xmath4 is considered as a functional of @xmath11 .
a detailed derivation of the hrt equations and a thorough discussion of their properties can be found in ref .
@xcite via diagrammatic methods @xcite . here
we provide a brief introduction of the physical content of the method together with a simpler , less technical derivation .    since the pioneering works by van der waals , a large part of liquid state theory is based on the idea that the strongly repulsive short range interaction , distinctive of the interatomic potentials , provides a constraint to the available phase space , determining an important contribution to the entropy of the fluid , while the mainly attractive , regular part of @xmath3 dominates the internal energy .
this is certainly true for the special , although prototypical , case in which @xmath3 has hard core at short distances , but holds also in the general case of soft core potentials .
the different role played by the singular and regular contributions to the interparticle potential is manifest in lattice gas models , where particles live on the sites of a lattice : hard core repulsion between particles gives rise to the single occupancy constraint , which allows to map the lattice gas onto an ising model with a spin - spin interaction proportional to the regular part of the interparticle potential .
mean field approaches developed for specific models , like van der waals equation in simple fluids or weiss theory in magnets , suggest that only the latter interaction triggers the order parameter fluctuations , responsible for the order - disorder transitions .    in fluid models , both the singular and the regular part contributions are included in the interaction @xmath3 which
should then be split in a reference , singular , part @xmath12 and in the regular remainder @xmath13 .
the partition function ( [ gran ] ) can be written in a suggestive form by introducing the  configurational density " @xmath14 whose thermal average is just the density profile @xmath15 , which reduces to the ordinary average density @xmath16 for @xmath17 . in the following
, the properties of the reference system will be considered as known .
the fluctuating order parameter of the liquid - vapour phase transition in fluids is then identified as the configurational density and @xmath4 is written as @xmath5 = \sum_n\frac{(z\,e^{-\frac{1}{2}\phi(0)})^{n}}{n!}\ ,   \int { \rm d } { \bf r_1} ... {\rm d } { \bf r_n } \ , e^{-\beta h_r}\ , e^{\frac{1}{2}\int { \rm d}{{\bf x}}_1 \,{\rm d}{{\bf x}}_2 \phi({{\bf
x}}_1-{{\bf x}}_2 ) \ , \hat\rho({{\bf x}}_1 ) \,\hat\rho({{\bf x}}_2 ) } \label{xi1}\ ] ] where the regular dimensionless interaction has been defined as @xmath18 and the reference hamiltonian is @xmath19 according to the previous remarks , @xmath20 defines a non - trivial integration measure , limiting the phase space available to the particles . by factorizing the reference @xmath0-particle canonical partition function @xmath21 ,
( [ xi1 ] ) can be also expressed as : @xmath5 =   \sum_n\frac{(z\,e^{-\frac{1}{2}\phi(0)})^{n}}{n!}\,e^{-\beta\,a_r(n ) } \,\left \langle    \exp\left ( \frac{1}{2}\int \frac{{\rm d}{{\bf q}}}{(2\pi)^d } \,\tilde\phi({{\bf q}})\ , \hat\rho({{\bf q}})\,\hat\rho(-{{\bf q}})\right ) \right \rangle_r \label{xi2}\ ] ] where the average @xmath22 refers to the gibbs measure of the reference fluid , and the fourier transforms of the potential @xmath23 and of the density fluctuations @xmath24 have been defined .
such a form is particularly suitable as a starting point for a reformulation of wilson s renormalization group picture in hamiltonian systems . in the momentum space rg formulation @xcite ,
a sequence of systems , labeled by a wave - vector @xmath25 , is defined by allowing density fluctuations only at wave - vectors @xmath26 larger than the cut - off @xmath25 .
the evolution of the physical properties of the model is then followed starting from the limit @xmath27 , where the system reduces to the reference one with a shift @xmath28 in the chemical potential , to the final value @xmath29 , where the cut - off is removed and the original partition function @xmath4 is retrieved . expression ( [ xi2 ] ) can easily accommodate the presence of a cut - off @xmath25 : @xmath30 =   \sum_n\frac{(z\,e^{-\frac{1}{2}\phi(0)})^{n}}{n!}\,e^{-\beta\,a_r(n ) } \,\left \langle    \exp\left ( \frac{1}{2}\int \frac{{\rm d}{{\bf q}}}{(2\pi)^d } \,\tilde\phi({{\bf q}})\,\theta ( q - q)\ , \hat\rho({{\bf q}})\,\hat\rho(-{{\bf q}})\right ) \right \rangle_r \label{xiq}\ ] ] where @xmath31 is the heaviside step function .
this definition embodies basic principle of hrt : a cut - off in density fluctuations is equivalent to a cut - off in the regular part of the two body interaction for the evaluation of the fluctuation corrections to the free energy , leading to the definition of a cut - off dependent potential @xmath32 .
it is customary to perform a legendre transform defining a cut - off dependent helmholtz free energy @xmath33 as a functional of the density profile @xmath15 : @xmath34 = -kt\,\log \xi_q[u({{\bf r } } ) ] - \int { \rm d}{{\bf r}}\ , \left ( u({{\bf r}})-kt\,\log z\right ) \,\rho({{\bf r } } ) \label{leg0}\ ] ] finally , it is also convenient to include in the definition of the free energy functional the mean field contribution due to the residual interaction @xmath35 : @xmath36 = a_q[\rho({{\bf r } } ) ]   + \frac{1}{2}\int { \rm d}{{\bf r}}\,{\rm d}{{\bf r}}^\prime \left [ w({{\bf r}}-{{\bf r}}^\prime ) - w_q({{\bf r}}-{{\bf r}}^\prime)\right ] \ , \rho({{\bf r } } ) \,\rho({{\bf r}}^\prime ) \label{leg}\ ] ] in this way , for @xmath37 , @xmath38 reduces to the mean field approximation of the full free energy , which is recovered as @xmath39 .
this way to artificially eliminate the long wave - length fluctuations from the thermodynamics is referred to as  sharp cut - off " .
however it is often more convenient to deal with a  smooth cut - off " procedure , where the @xmath40-function in the definition of @xmath41 is substituted by a smoothed step function , denoted in the following as @xmath42 which is continuous and differentiable for every @xmath43 and tends to the step function in the @xmath44 limit : @xmath45      after having defined the cut - off dependent free energy , @xmath38 via eq .
( [ leg ] ) , the next step is to write the exact evolution ( or flow ) equation governing the change in the physical properties of the model when the cut - off @xmath25 is varied .
this goal is most easily achieved in the smooth cut - off formulation , where an infinitesimal increase of the cut - off @xmath25 leads to an infinitesimal change in the potential @xmath41 . according to standard first order perturbation theory @xcite we get , in a uniform fluid , @xmath46 where @xmath47 is the total correlation function of a fluid interacting through a potential @xmath48 and @xmath49 .
this smooth cut - off hrt equation has been obtained for the first time in ref .
@xcite as an attempt to apply the renormalization group ideas to microscopic models .
note that the total correlation function @xmath50 rather than the radial distribution function @xmath51 appears in the evolution equation ( [ evol1 ] ) because the mean field contribution has been subtracted in the definition of the @xmath25-dependent free energy ( [ leg ] ) .
the total correlation function can be conveniently expressed in terms of the direct correlation function @xmath52 by use of the ornstein - zernike relation : @xmath53 the limit @xmath54 in eq .
( [ wqsmooth ] ) requires some care , as already recognized in the description of the sharp cut - off hrt equation and also discussed in field theoretical approaches @xcite , because of the emergence of a discontinuity in @xmath55 precisely on the same surface @xmath56 where the @xmath57-function singularity develops in @xmath58 .
such a discontinuity originates from the random phase contribution to the exact direct correlation function , being the single term where the regular part of the potential appears outside any integration .
if @xmath59 denotes the direct correlation function of the reference fluid , the first terms of the perturbative expansion read @xmath60 all the remaining contributions , being integrated , do not display discontinuities in the whole @xmath26-space .
this suggests to introduce a modified correlation function @xmath61 which remains continuous also in the sharp cut - off limit @xmath54 : @xmath62 interestingly , this newly defined function is just the second functional derivative of the free energy functional @xmath63 $ ] with respect to density variations , evaluated at uniform density : @xmath64 where the last term is due to the ideal gas contribution , included in our definition of @xmath38 .
we first write the right hand side of eq .
( [ evol1 ] ) by introducing the fourier transforms : @xmath65 } \ , \frac{\partial \tilde\phi_q({{\bf q}})}{\partial q }
\label{smooth}\ ] ] the easiest way to take the sharp cut - off limit of this flow equation is to introduce the differential operator @xmath66 acting only on the q - dependence of @xmath67 ( while @xmath61 is treated as @xmath25-independent ) : @xmath68   \label{evol2}\ ] ] taking the @xmath54 limit and writing @xmath69 as a finite difference , eq . ( [ evol2 ] )
is then written in the form @xcite : @xmath70 \label{sharp}\ ] ] where the integration is limited to the spherical surface ( in dimension @xmath2 ) defined by @xmath56 and the function @xmath71 is defined by the ornstein - zernike relation as : @xmath72 equations ( [ smooth ] ) and ( [ sharp ] ) govern the hrt flow of the free energy in the smooth and sharp formulations respectively . as previously mentioned , the initial condition at @xmath27 coincides in both cases with the known mean field approximation , while for @xmath39 all density fluctuations are included and the exact free energy is in principle recovered if @xmath71 is known .
hrt then provides the exact evolution of the free energy in terms of the two body correlation function @xmath73 ( or equivalently @xmath71 ) whose flow is determined by similar equations involving the three and four particle correlations @xcite .
these equations , which can be derived along the same lines previously sketched , are most easily obtained by functional differentiation of ( [ smooth ] ) and ( [ sharp ] ) with respect to @xmath15 leading to a hierarchy of differential equations describing the rg flow of the free energy and of the @xmath74-particle direct correlation functions @xmath75 , formally defined as the @xmath76 functional derivative of @xmath77 with respect to the density profile @xmath15 evaluated at uniform density @xmath78 .
the ideal gas contribution is included in our definition of @xmath38 leading to an additional term @xmath79 with respect to the standard definition of the many particle direct correlation functions @xcite .
as an example , we explicitly report the flow equations for the previously defined two point function in smooth cut - off formulation @xcite @xmath80 and in sharp cut - off formulation @xcite : @xmath81 where , in the latter equation , @xmath82 is defined as : @xmath83^{-1 } & \text{for $ \vert { { \bf k}}+{{\bf q}}\vert < q$ } \\ { \cal f}_q({{\bf k}}+{{\bf q } } ) & \text{for $ \vert { { \bf k}}+{{\bf q}}\vert > q$ }   \end{cases}\ ] ] note that , in the limit @xmath84 , these flow equations reproduce exactly the second density derivative of the corresponding free energy evolution ( [ smooth],[sharp ] ) , showing that the density dependence of the free energy flow encodes some ( partial ) information contained in the full infinite hierarchy .      since the first derivation of the hrt equations in the sharp cut - off formulation @xcite a close similarity to some previous differential rg formulation was noticed .
in particular it was shown that close to a critical point , in the long wave - length limit the full hrt hierarchy precisely reproduces the flow equations for the one particle irreducible rg generator @xcite .
the analogy between the hierarchical reference theory and functional rg approaches is however even deeper , as recently shown by caillol @xcite . to highlight the connection between the two approaches it is convenient to start from the definition of the partition function @xmath85 ( [ xiq ] ) .
as pointed out in an early attempt to apply rg techniques to microscopic models @xcite , by performing a hubbard - stratonovich transformation , the partition function can be expressed as a functional integral over an auxiliary scalar field @xmath86 governed by an effective action @xmath87 $ ] .
this commonly adopted procedure , however , is rather formal and does not lend a well defined physical meaning for the hubbard stratonovich field @xmath86 .
it is possible to circumvent this problem by first introducing a reference effective action @xmath88 $ ] by : @xmath89 = -\log \left [   \sum_n\frac{z^{n}}{n!}\ ,   \int { \rm d } { \bf r_1} ... {\rm d } { \bf r_n } \ , e^{-\beta\sum_{i < j } v_r({\bf r}_i - { \bf r}_j)}\ , \prod_{{\bf q}}\delta(\hat\rho({{\bf q}})-\psi_{{\bf q } } ) \right ]   \label{sr}\ ] ] this definition clearly forces the effective field @xmath86 to faithfully represent density fluctuations .
the exact partition function is then written as @xmath90 & = & \int { \cal d } \psi \ , e^{-s_q[\psi_{{\bf q } } ] } \\ s_q[\psi_{{\bf q } } ] & = & s_r[\psi_{{\bf q } } ] - \frac{1}{2v}\,\sum_{{\bf q}}\tilde\phi_q({{\bf q}})\ , \psi_{{\bf q}}\,\psi_{-{{\bf q } } } + \frac{1}{v } \sum_{{\bf q}}\left ( \frac{\phi(0)}{2 } + \beta \tilde u({{\bf q}})\right ) \psi_{-{{\bf q } } }   \label{action}\end{aligned}\ ] ] this expression shows that @xmath87 $ ] can be written as the sum of the effective action in the absence of any cut - off @xmath91 $ ] ( i.e. the value attained at @xmath29 ) plus a  momentum dependent mass term " @xcite : @xmath92 = s[\psi_{{\bf q } } ] + \frac{1}{2v}\,\sum_{{\bf q}}r_q({{\bf q } } ) \
, \psi_{{\bf q}}\,\psi_{-{{\bf q}}}\ ] ] with @xmath93   \label{mass}\ ] ] following wetterich we now perform a  modified legendre transform " @xcite introducing the field @xmath94 conjugated to the external potential @xmath95 , which , due to the definition ( [ action ] ) coincides with the statistical average of the field @xmath86 : @xmath96 = -\log\xi_q + \frac{1}{v}\,\sum_{{\bf q}}\left ( \tilde u({{\bf q}})-kt\,\log z\right ) \,\rho_{-{{\bf q}}}-   \frac{1}{2v}\,\sum_{{\bf q}}r_q({{\bf q}})\ , \rho_{{\bf q}}\,\rho_{-{{\bf q } } } \label{leg1}\ ] ] comparing this definition to eq .
( [ leg ] ) it is clear that the average action coincides with the modified free energy defined in hrt : @xmath97 , providing the link between hrt , based on the gradual inclusion of the fourier components of the regular part of the physical interaction , and the  effective average action " rg method built upon the idea of adding a cut - off and momentum dependent mass term @xmath98 to the original action . from eq .
( [ mass ] ) it follows that for @xmath29 the mass term vanishes identically , while at fixed @xmath25 , @xmath98 is large only for @xmath99 providing an effective cut - off for the long wave - length fluctuations , as requested in all the implementations of the momentum space rg ideas @xcite .
a difference in the initial condition however persists between the average action method and hrt : in the rg approach , @xmath98 diverges as @xmath37 , thus providing an infinite mass to all fluctuations and leading to a soluble problem .
instead in hrt we get @xmath100 : for this special but finite value of the mass term the regular part of the interaction cancels exactly in @xmath87 $ ] and the model reduces to the soluble reference system .
the procedure adopted in hrt is reminiscent of some recent development in functional rg methods @xcite .
the solution of the infinite hierarchy of equations previously derived needs some approximation scheme .
two kinds of approximations are customary in this type of problems : truncations and closures .
a truncation of the hierarchy simply corresponds to neglecting @xmath74-particle correlations for @xmath74 larger than a maximum value @xmath101 ( typically @xmath102 ) .
such a procedure has been followed in the early investigations of rg equations but it is rigorously justified only near a critical point in spatial dimension close to @xmath103 @xcite , while it is not expected to provide an accurate description of a system , like a liquid , where short range correlations are important . instead ,
most of the applications of hrt , both in the smooth and sharp cut - off formulations , adopted a non perturbative closure to the first equation of the hierarchy ( [ smooth],[sharp ] ) and only this type of approximation will be discussed here .
non perturbative closures in fact allow to preserve the convexity of the free energy , providing a physically satisfactory description of phase transitions between homogeneous phases , together with a picture of the critical region fully consistent with the scaling hypothesis .
the evolution equation for the free energy involves the knowledge of the two particle direct correlation function in the presence of a cut - off @xmath25 .
a closure is a parametrization of the direct correlation function expressing @xmath104 in terms of the free energy density @xmath105 . as previously noticed , the compressibility relation , which follows directly from eq .
( [ second ] ) , @xmath106 plays a special role in hrt because , when implemented into the free energy evolution equation , allows to include in the very first equation some information contained into the whole hierarchy .
the simplest non perturbative closure of the first hrt equation is inspired by the well known random phase approximation ( rpa ) in liquid state theory @xcite : @xmath107 the exact initial condition at @xmath37 is obtained for @xmath108 , while the compressibility relation ( [ komp ] ) can be enforced at all @xmath25 via a suitable choice of the parameter @xmath109 . when this closure is substituted into eq .
( [ smooth ] ) or ( [ sharp ] ) the resulting partial differential equation for the free energy @xmath38 as a function of density @xmath16 and cut - off @xmath25 can be solved for a given value of the temperature .
interestingly , this approximate closure can be equivalently formulated starting from the evolution equations ( [ smooth ] ) or ( [ sharp ] ) for the full functional @xmath110 $ ] @xcite : by restricting its form through a suitably chosen parametrization , the functional flow equation gives rise to partial differential equations for the parameters .
for instance , the rpa - like closure is obtained by use of a simple weighted density approximation @xcite for the free energy : @xmath111 = a_r[\rho({{\bf r } } ) ] + \int { \rm d}{{\bf r}}\,f_q(\bar\rho({{\bf r } } ) )   \label{wda}\ ] ] where @xmath112 $ ] is the exact free energy functional of the reference system and the function @xmath113 is the parameter whose evolution will be determined by the flow equation . the weighted density in eq . ( [ wda ] ) @xmath114
is defined by the local average : @xmath115 with a normalized weight function given by @xmath116^{1/2}$ ] @xcite . by substituting this form into the evolution equation ( [ smooth ] ) or ( [ sharp ] ) ,
we obtain the same expression found by use of the ansatz ( [ rpa ] ) if we identify @xmath117 . this class of closures , in the nprg framework
, are known as local potential approximation ( lpa ) because of the similarity between eq .
( [ wda ] ) and the usual local density approximation of density functional theory .
clearly , more refined approximations to the direct correlation function , or equivalently to the free energy functional , may be devised , leading to more accurate representation of the exact free energy flow . in fluids with a hard sphere reference system ( in the continuum or on a lattice ) ,
( [ rpa ] ) can be generalized by including a further contribution , non vanishing only inside the hard sphere diameter , so as to satisfy the core condition @xmath118 .
this constraint can be implemented exactly in the smooth cut - off formulation for yukawa potentials @xcite thanks to the exact analytical solution of the ornstein - zernike equation for this model @xcite .
analogous closures might be implemented on the flow equation for the two particle direct correlation function ( [ csmooth ] ) or ( [ csharp ] ) . a closure to the second equation of the hierarchy amounts to expressing the three and four particle correlations in terms of the two point function .
a possible closure has been proposed already in the earliest stage of hrt @xcite but it has not been implemented yet : @xmath119 \label{c4}\end{aligned}\ ] ] this class of approximations has the advantage of allowing for a non trivial anomalous dimension exponent at the critical point @xcite .
a similar closure has been independently proposed in the context of functional rg @xcite and a simplified form has been numerically investigated in the same framework @xcite .    up to now , all the applications of hrt to microscopic models adopted a closure of the form ( [ rpa ] ) with or without the additional optimization of the potential inside the hard sphere diameter to enforce the core condition .
the self consistent ornstein zernike approximation ( scoza ) is a well known liquid state theory which proved successful in the determination of phase diagrams in fluids characterized by a pair interaction written as the sum of a singular hard sphere part plus a yukawa tail @xmath120 @xcite .
originally , scoza has been formulated as a generalization of the mean spherical approximation @xcite for the direct correlation function outside the core : @xmath121 while enforcing the core condition @xmath122 to define @xmath123 inside the hard sphere diameter .
the parameter @xmath124 is determined by requiring thermodynamic consistency between the internal energy and the compressibility routes to thermodynamics : @xmath125 the two exact expressions ( [ scozac ] ) , together with the core condition and the parametrization ( [ scoza ] ) outside the hard sphere diameter , give rise to a partial differential equation for the free energy as a function of @xmath16 and @xmath126 .
remarkably , this equation is exactly reproduced by hrt in the smooth cut - off formulation ( [ smooth0 ] ) by use of a linear turning on of the interaction ( [ wqsmooth ] ) : @xmath127 where @xmath128 is the switching - on parameter @xcite .
clearly this form does not satisfy the requirements stated in section [ sechrt ] for a selective introduction of density fluctuations : while scoza just tunes the amplitude of @xmath129 , hrt changes simultaneously both amplitude and range of the potential by imposing a cut - off @xmath25 on its long wave - length fourier components thereby qualifying as a microscopic implementation of rg ideas .
in fact , scoza has been shown to provide a good description of the non universal properties in the critical region , as well as non classical critical exponents but fails in reproducing the correct scaling picture predicted by rg and hrt @xcite .
an attempt to reconcile hrt and scoza has been formulated in ref .
@xcite , where a further adjustable parameter is added in the rpa - like closure of the hrt equation , to enforce thermodynamic consistency between the compressibility and the internal energy route .
one of the goals which motivated the formulation of hrt was to provide a description of the critical properties of realistic fluid models consistent with the rg picture .
critical exponents , governing the singularities of thermodynamic quantities in the critical region , acquire universal values quoted in table [ table1 ] @xcite . at the critical density ,
the constant volume specific heat diverges as @xmath130 as the reduced temperature @xmath131 vanishes ; the isothermal compressibility as @xmath132 and the correlation length as @xmath133 .
the shape of the coexistence curve is asymptotically given by @xmath134 , while the equation of state on the critical isotherm is given by @xmath135 .
finally , precisely at the critical point , the pair correlation displays a power law tail of the form @xmath136 .    in the sharp cut - off formulation
, universality is indeed a direct consequence of the asymptotic form of the hrt hierarchy close to the critical point @xcite : all the microscopic details of the model disappear from the asymptotic equations and the full hierarchy acquires a universal form equivalent to a previously studied rg differential generator @xcite .
scaling laws and non classical values for the critical exponents immediately follow from this identification . in the smooth cut - off formulation , the universality of critical phenomena is less apparent because the specific form of the chosen cut - off function @xmath137 ( [ wqsmooth],[mass ] ) survives also in the asymptotic critical region .
universality of the critical properties can then be proved only order by order in the @xmath138 expansion and no general proof is available @xcite . here
we will briefly review the critical properties of the approximate non perturbative class of closures previously introduced both in the sharp and in the smooth cut - off formulation ( [ rpa ] ) .
moreover we will sketch the argument leading to the conclusion that already at the level of this simple approximation , hrt always gives a convex free energy , even at coexistence .
the unphysical divergence of the compressibility at the phase boundary predicted by the sharp cut - off closure @xcite is circumvented by the smooth cut - off formulation which then provides a satisfactory description of first order transitions @xcite .
other , more sophisticated closures have been proposed . as previously noticed ,
the second equation of the hrt hierarchy describes the evolution of the two particle direct correlation function in terms of the three and four particle correlations ( [ csmooth],[csharp ] ) . by use of the closure ( [ c3],[c4 ] )
it has been shown that the critical exponents are correctly reproduced to second order in @xmath138 , including a non vanishing anomalous dimension @xmath139 @xcite .
unfortunately this quite promising non perturbative closure has not been thoroughly investigated in microscopic models , probably due to the problems posed by the numerical solution of the resulting non linear partial differential equation .
recently , a similar ( although simpler ) approximation has been studied in the asymptotic regime , in the framework of nprg , leading to encouraging results for the universal quantities @xcite .      in the sharp cut - off
hrt , the free energy evolution equation ( [ sharp ] ) depends on the two point function @xmath71 evaluated  on shell " , i.e. on the surface @xmath56 .
this means that in the late stage of the integration , when @xmath25 is small , only the long wave - length physics contributes to the free energy flow .
moreover , in the critical region , @xmath71 diverges at small @xmath140 and then we can disregard the additive term in eq .
( [ sharp ] ) and the sub - leading @xmath141 contribution .
finally , the evolution equation acquires a universal form @xmath142 where @xmath143 is a geometric factor equal to the ratio between the volume of the unit sphere in @xmath2 dimensions and @xmath144 .
when we insert a closure of the form ( [ rpa ] ) at right hand side , recalling eq .
( [ calf ] ) , we need to evaluate the direct correlation function at wave - vector @xmath39 . by exploiting the analyticity of @xmath145 and @xmath146 around @xmath147 , we can write @xmath148 where @xmath149 determines the curvature of the direct correlation function in this approximation and attains a finite , positive limit at the critical point . recalling the compressibility sum rule ( [ komp ] ) and inserting this asymptotic form into eq .
( [ sharp2 ] ) we obtain a closed partial differential equation expressing the free energy flow at long wave - length and in the critical region . the rg structure of this equation can be highlighted by a suitable rescaling , as shown in refs .
@xcite , which allows the explicit calculation of the universal properties , like critical exponents and scaling functions , via the fixed point analysis .
the analyticity in the momentum dependence of the direct correlation function , implied by the rpa - like ansatz ( [ rpa ] ) , immediately gives the vanishing of the anomalous dimension , or fisher exponent , @xmath150 in any spatial dimension .
this result is exact for @xmath151 , when the critical exponents are classical , while for @xmath152 such a simple non perturbative approximation to the first hrt equation gives critical properties correct to first order in the @xmath138 expansion .
however this closure leads to a qualitatively incorrect prediction in @xmath153 , when the anomalous dimension is crucial in ensuring the existence of a fixed point solution .
few critical exponents in @xmath154 are reported in table [ table1 ] together with the accepted values obtained from extrapolations of the high temperature expansion in the ising model @xcite .
0.2 cm    .hrt estimates of the critical exponents and compressibility amplitude ratio in three dimensions for the sharp and the smooth formulations of hrt compared to the exact values @xcite obtained by extrapolation of high - temperature series expansions .
the smooth cut - off procedure adopted here is detailed in section [ secsmooth ] . [ cols="^,^,^,^,^,^,^,^",options="header " , ]     the ising model on a cubic lattice was the first lattice system investigated by sharp cut - off hrt @xcite by use of a more sophisticated closure , belonging to the class previously discussed , which consistently preserves the spin normalization constraint @xmath155 .
a thorough analysis of lattice models via the smooth cut - off version of hrt has not been attempted yet , although a very similar procedure , put forward in the framework of non perturbative renormalization group , has been recently developed @xcite .
the results shown in @xcite provided a very accurate estimate of the critical temperature of the ising model and compared favourably with the coexistence curve predicted by high temperature expansion .
a numerical study of the same hrt equation aimed at determining the effective parameters of a scalar field theory able to reproduce the correct pre - critical behaviour of the ising model was performed in ref .
interestingly , hrt showed that a @xmath156 representation of the long wave - length spin fluctuations was inadequate for a faithful representation of the non - universal properties of the ising model , while the inclusion of a @xmath157 term allowed to considerably improve the description of the critical region of the model .    in ref .
@xcite the hrt formalism was generalized to the case of an uniaxial spin model where @xmath158 can assume three values : @xmath159 .
this system can be easily mapped into a symmetric binary mixture on a lattice with opposite signs for like and unlike interactions .
if the coupling is limited to nearest neighbours , the blume - capel model is obtained , while for long range interactions decaying as @xmath160 , the lattice restricted primitive model ( lrpm ) is recovered .
the phase diagrams of both models in the symmetry plane ( zero magnetization ) were obtained by integrating the hrt flow equation @xcite .
the topologies of the phase boundaries in the density - temperature plane , shown in fig .
[ figlrpm ] for the lrpm , turn out to be very similar : at high temperatures a critical @xmath124-line separates a charge - ordered fluid phase at low density from a homogeneous plasma phase .
the critical line ends at a tricritical point where it merges into the liquid - vapour coexistence curve .
the comparison with available simulation data for the lrpm shows that the structure of the phase diagram and the universality classes of the transitions are correctly reproduced , although some quantitative discrepancy is present , suggesting that a better representation of pair correlations should be devised .     plane .
hrt results : the tricritical point is marked by an open star .
full dots give the @xmath124-line above the tricritical temperature and the liquid - vapour binodal at lower temperatures .
the scoza data for the @xmath124-line are shown by the dotted line .
simulation data from refs .
@xcite are represented by filled triangles and open squares respectively .
, width=264,height=264 ]      the hrt formalism can be easily generalized to the case of binary mixtures @xcite : the flow equations in the smooth ( [ smooth ] ) and sharp cut - off ( [ sharp ] ) formulation remain valid for mixtures provided the two body interaction @xmath146 and the pair correlations @xmath71 and @xmath61 are substituted by matrices whose indices run over the two species . the ornstein - zernike relation ( [ calf ] ) acquires the form @xmath161 =   \delta^{ik } \label{calfmat}\ ] ] while the compressibility sum rule ( [ komp ] ) becomes a set of three independent equations which can be written as : @xmath162 the derivation of the hrt evolution equations straightforwardly follows the same steps sketched in section [ sec : exact ] . by adopting the smooth cut - off procedure ,
the resulting expression is : @xmath163^{-1 } \
, \frac{\partial \tilde\phi_q({{\bf q}})}{\partial q } \right \ } \label{smoothmat}\ ] ] where matrix products are understood at right hand side and  tr " stands for the trace operation . in the sharp cut - off limit this form reduces to @xmath164 \label{sharpmat}\ ] ] the trace of the logarithm at right hand side of eq .
( [ sharpmat ] ) must be interpreted in matrix sense and coincides with the logarithm of the determinant of the @xmath165 matrix @xmath166 $ ] . by use of a rpa - like closure of the form ( [ rpa ] ) , the flow equation then becomes a partial differential equation in three independent variables : the cut - off @xmath25 and the densities of the two species @xmath167 .
the analysis of criticality in binary mixtures is performed following the procedure developed in simple fluids @xcite : several universality classes are predicted by the hrt equations , but most of the possible fixed points turn out to be unstable and the critical properties at a generic point along a transition line is of the ising type .
actually , two distinct fixed points correspond to critical exponents related to the case of a one component scalar order parameter .
one is weakly unstable and gives rise to strictly ising - like critical properties , while the other is stable and correctly reproduces the so called phenomenological approach by fisher which , starting from the concept of field mixing , predicts the occurrence of some weak exponent renormalization @xcite : the long wavelength limit of a structure factor @xmath168 at the critical density and concentration diverges with an exponent @xmath169 ( with @xmath170 ) , a @xmath171 larger than the corresponding value for one component fluids , while the isothermal compressibility displays a weak divergence governed by the exponent @xmath172 .
hrt allows to identify the thermodynamic perturbation corresponding to the crossover between these two fixed points @xcite , thereby providing a microscopic interpretation of the experimental difficulty in observing the true asymptotic critical behaviour .
special points in the phase diagram might display genuinely different critical indices , like un - renormalized ising or tricritical exponents .
hrt also naturally suggests an unbiased definition of the order parameter at a critical point of binary fluids @xcite , which in general can not be simply identified either as density or concentration , because critical fluctuations are expected to affect both quantities . the critical point in binary fluids
is signaled by the divergence of the leading eigenvalue of the ( symmetric ) matrix of the structure factors @xmath173 for @xmath84 .
the structure factors are a direct measure of density fluctuations and can be expressed as the equilibrium average of a product of density operators @xmath174 . to discriminate between demixing and liquid - vapour transitions in binary mixtures
it is convenient to introduce suitable linear combinations of @xmath173 , the bathia - thornton structure factors @xcite , which properly characterize the density and concentration fluctuations : @xmath175 where @xmath176 , @xmath177 are the number concentrations of the two species and the total density , respectively .
let us denote by @xmath178 leading eigenvalue of the @xmath165 matrix built with the previously defined bathia thornton structure factors in the @xmath84 limit , and @xmath179 be the corresponding normalized eigenvector . according to these definitions
, the diverging eigenvalue provides a measure of the fluctuations of the particular linear combination of density and concentration @xmath180 ( for @xmath84 ) : @xmath181 , suggesting the identification of @xmath182 as the order parameter of the transition .
it is convenient to pictorially represent the order parameter as an arrow in the @xmath183 plane pointing in the direction identified by @xmath40 . remarkably , the  mixing angle " @xmath40 is related to
well defined thermodynamic quantities @xcite : @xmath184_{coex}&=&\rho\tan\theta \quad . \label{termo2}\end{aligned}\ ] ] the first equation is written in terms of the partial molar volumes @xmath185 which attain a finite limit at the critical point .
the second , equivalent , relation expresses @xmath40 in terms of the differences of concentration @xmath186 and molar volumes @xmath187 at coexistence on an isotherm  isobar near the critical point .
clearly , the direction of the order parameter is defined at each point of a critical line and varies continuously along that line .
while in some mixture , the character of the transition remains almost unaltered by varying the thermodynamic state , in other cases it may instead experience considerable changes , passing from demixing to liquid - vapour by reducing the density .
a plot of the mean field critical lines for two lennard - jones binary fluids mimicking two rare gas mixtures , shown in fig .
[ figmix ] , illustrates these two possibilities . for convenience , the total packing fraction @xmath188 and the volume fraction @xmath189
are used in place of total density @xmath16 and number concentration @xmath190 .
defined in the text : an arrow parallel to the density axis indicates a liquid - vapour transition , an arrow along the volume fraction axis , a demixing transition .
, width=264,height=264 ]    the numerical integration of the hrt equation poses some challenge due to the presence of three independent variables : cut - off @xmath25 , density @xmath16 and concentration @xmath190 .
special attention must be paid to the stability requirements and to the possible occurrence of singularities induced by the presence of large coexistence regions in the phase diagram .
a stable numerical code has been developed for two classes of models of binary fluids : the rare gas mixtures @xcite and the symmetric yukawa mixtures @xcite , which include the ising fluids , defined as hard spheres with an internal ising - like internal degree of freedom whose relative orientation specifies the sign of the pair interaction @xcite . in both cases , the hrt equation in the sharp cut - off formulation has been numerically integrated .
here we just review few results for the model of rare gas mixtures , where extensive experimental data are available @xcite .
the argon - krypton system is modeled as a mixture of particles interacting via lennard - jones ( lj ) potentials @xcite .
each lj interaction has been preliminarily mapped into a reference hard sphere potential plus a lj attractive tail .
the effective hard sphere diameter does not coincide with the lj parameter @xmath191 and is determined according to the standard barker - henderson prescription @xcite which gives rise to a weak temperature dependence .
the bare lj parameters adopted in the numerical calculation have been fixed to reproduce the location of the critical points of the two pure components : @xmath192  , @xmath193  , @xmath194 k and @xmath195 k. the strength of the unlike attractive interaction has been set to @xmath196 k , and additivity of the effective hard sphere diameters has been assumed .
plots of the coexistence regions in the packing fraction / molar concentration plane are given for two temperatures in figs .
[ figarkr1-rx]-[figarkr2-rx ] , where the tie - lines , connecting thermodynamic states in equilibrium , are also shown .
the argon - krypton mixture belongs to the topology of phase diagram illustrated in the left panel of fig [ figmix ] , where the critical points of the two species are connected by a line of critical points .    .
tie lines in the two - phase regions are shown . in the coexistence regions
darker colors correspond to higher pressures .
, width=226,height=226 ]    .
tie lines in the two - phase regions are shown . in the coexistence regions
darker colors correspond to higher pressures .
, width=226,height=226 ]    a different view of the coexistence regions at the same temperatures can be seen in the molar concentration / pressure plane in fig .
[ figarkr - xp ] : at the higher temperature pure argon is in the fluid phase , while at the lower temperature the coexistence regions of the two components merge and define a unique , connected domain .
, right panel @xmath197 .
experimental phase boundaries @xcite are shown as red dots .
, width=264,height=264 ]    a good agreement between theory and experiments @xcite can be appreciated for this system . finally , it is instructive to check that , within hrt , the full coexistence region collapses to a single line when intensive thermodynamic quantities are selected , as required by thermodynamics : in fig .
[ figarkr - mup ] the same coexistence regions are plotted in a diagram where the difference between the chemical potentials of the two species and the pressure are considered as independent variables .
, right panel @xmath197 .
note the collapse of the full coexistence region into a single line .
, width=264,height=264 ]    the phase diagram of the neon - krypton mixture instead displays a different topology , belonging to the class depicted in the right panel of fig .
[ figmix ] .
the same lennard - jones model used for the argon - krypton mixture has been adopted with parameters appropriate for the neon - neon ( @xmath198  , @xmath199 k ) and neon - krypton interactions ( @xmath200 k ) .
the coexistence regions of the two pure components remain separated in the whole phase diagram : while pure neon is in the gas phase in a large temperature range , the neon - krypton mixture displays phase coexistence between states of considerably different composition .
the order parameter smoothly changes character signaling a liquid - vapour transition at low neon concentration and a demixing transition at higher pressures as already captured by the mean field calculations of fig [ figmix ] .
the coexistence region in the three planes packing fraction / molar concentration , molar concentration / pressure and chemical potential / pressure are shown in figs .
[ fignekr1-rx ] , [ fignekr2-rx ] , [ fignekr - xp ] and [ fignekr - mup ] respectively .    .
tie lines in the two - phase regions are shown . in the coexistence regions
darker colors correspond to higher pressures .
, width=226,height=226 ]    .
tie lines in the two - phase regions are shown . in the coexistence regions
darker colors correspond to higher pressures .
, width=226,height=226 ]    , right panel @xmath201 .
experimental phase boundaries @xcite are shown as red dots .
, width=264,height=264 ]    the comparison with the experimental data @xcite shows good agreement at not too high pressure , while some discrepancy is present at high pressures , where the repulsive part of the interaction plays an important role . as discussed in ref .
@xcite , it is likely that a better representation of the soft core part of the interatomic potential must be employed for a fully quantitative description of demixing transitions at high pressures .    , right panel @xmath202 . , width=264,height=264 ]
hrt was originally developed in the framework of classical statistical physics but its formalism can be generalized to quantum hamiltonians .
lattice spin models have been considered in ref .
@xcite , where the reference system was identified as free spins in an external magnetic field which couples with the order parameter .
quantum antiferromagnets on an hypercubic lattice in @xmath2 dimensions are described by the isotropic , nearest neighbour heisenberg hamiltonian @xmath203 and require the introduction of a staggered external field @xmath204 . here , @xmath205 are spin @xmath206 operators defined on the sites @xmath207 of the lattice and @xmath208 is the antiferromagnetic vector in @xmath2 dimensions .
the reference system is defined by the hamiltonian ( [ heis ] ) for @xmath209 and gives rise to a solvable , single site problem . following the standard sharp cut - off formulation of hrt ,
the spin - spin interaction is gradually turned on , by including its fourier components starting from the boundary of the antiferromagnetic brillouin zone .
the exact evolution equation for the free energy per site in sharp cut - off formulation was obtained in ref .
@xcite and involves the two body dynamical spin - spin correlations in imaginary time : @xmath210 + \\ & & \log \left [ \left ( 1-{\cal f}_q^{xx}({{\bf k}},\omega_n ) \tilde j({{\bf
k}})\right ) \ , \left ( 1-{\cal f}_q^{xx}({{\bf k}}+{{\bf g}},\omega_n ) \tilde
j({{\bf k}}+{{\bf g}})\right ) +   { \cal f}_q^{xy}({{\bf k}},\omega_n)\,{\cal f}_q^{xy}({{\bf k}}+{{\bf g}},\omega_n)\ , \tilde j({{\bf k}})\tilde j({{\bf k}}+{{\bf g } } ) \right ] \big \ } \nonumber   \label{gian}\end{aligned}\ ] ] where @xmath211 is the fourier transform of the spin - spin nearest negighbour interaction and @xmath25 ranges from @xmath212 , where fluctuations are absent and the mean field approximation is recovered , to @xmath29 , where all fourier components are introduced .
the structure of this equation is rather similar to its classical counterpart ( [ lsharp ] ) with two main differences : the presence of off diagonal two body correlations @xmath213 , due to the non commutativity of the spin components , and the summation over the matsubara frequencies @xmath214 .
the imaginary time variable in fact naturally emerges when quantum perturbation theory is employed to evaluate the effects of an additional infinitesimal momentum shell in the spin - spin interaction .
notice that in this formulation , the cut - off on fluctuations is set only on the spatial fourier components , while the matsubara frequency variable runs on the whole real axis for all @xmath25 s .
the quantum hrt ( qhrt ) equation appropriate for the antiferromagnetic heisenberg model ( [ heis ] ) was studied by employing a simple closure similar in spirit to the classical rpa - like approximation ( [ lrpa ] ) : the structure of the dynamical correlations was evaluated at mean field level and only  mass renormalization " was allowed by imposing thermodynamic consistency between the zero frequency limit of the dynamical correlations at wave - vector @xmath215 and the derivatives of the free energy @xmath38 with respect to the staggered magnetization .
the explicit form of the adopted closure reads : @xmath216 where @xmath217 is the staggered magnetization , and the two cut - off dependent parameters @xmath218 are related to the transverse and longitudinal staggered susceptibilities : @xmath219\\ \frac{\partial^2}{\partial m_s^2 } \left ( \frac{{\cal a}_q}{n}\right ) & = & \mu_\parallel + \tilde j({{\bf g } } ) \end{aligned}\ ] ] in fig .
[ figheis ] the constant field specific heat and the spontaneous magnetization of a heisenberg antiferromagnet on the cubic lattice below the critical temperature are shown for two choices of the spin @xmath206 : the extreme quantum value @xmath220 and the classical limit @xmath221 .
monte carlo simulations for @xmath222 @xcite are also reported .     for spin @xmath223 and @xmath220 .
left panel : qhrt spontaneous magnetization for @xmath220 and @xmath223 ( classical limit ) .
right panel : constant field specific heat at @xmath224 as a function of temperature compared with classical ( full dots ) @xcite and quantum ( full squares ) @xcite monte carlo simulations , respectively . in the @xmath223
limit the magnetization is normalized to @xmath206 and the specific heat to @xmath225 .
, width=264,height=264 ]    a more detailed comparison between theory and simulations can be found in @xcite , showing that the performances of qhrt have generally a lower quality with respect to the previously shown applications to classical models , suggesting that a better representation of static and dynamic spin correlations has to be developed : the different role of longitudinal and transverse correlations is probably not fully grasped by the adopted single mode approximation , which did not prove accurate enough to provide a consistent picture of the crossover between the renormalized classical and the quantum critical regime of the non linear sigma model @xcite .
the effects of easy - axis anisotropy have been also investigated by use of qhrt @xcite , as well as the role of a weak next nearest neighbour interaction @xcite .
however , a qhrt analysis of quantum ( i.e. zero temperature ) phase transitions , or of the occurrence of off - diagonal long range order has not been attempted yet . as shown in the related non perturbative rg approach for the bose - hubbard model @xcite , the investigation of bose condensation by rg methods requires a different cut - off procedure and it may be more convenient to set a cut - off on the kinetic energy term rather than on the two body interaction . according to this scheme , the reference system becomes a purely classical model , while quantum fluctuations are gradually introduced , leading to the appropriate free energy flow equation .
the smooth cut - off formulation of hrt has not been applied to quantum systems yet , although the derivation of the flow equations does not present any additional difficulty .
a closer connection between qhrt and non perturbative renormalization group methods may provide some clue on the structure of more refined closures to the flow equation .
a proper description of first and second order phase transitions within liquid state theory requires a careful account of long wavelength fluctuations .
the hierarchical reference theory of fluids was developed precisely to achieve this goal . while the general framework of the theory was set since the very first formulation of hrt , the implementations to specific physical systems required the development of approximation schemes which may vary according to the individual features of the model . in this review
we tried to summarize the efforts in this direction carried out in the last ten years . among the features shared by all the non - perturbative closures of the hrt equation we may mention the correct scaling form of the thermodynamics and correlations in the critical region , the prediction of non - classical critical exponents and the exact implementation of the maxwell construction at first order transitions , showing that the inclusion of long wave - length fluctuations gives rise to coexistence regions at constant pressure and chemical potential .
the latter feature is particularly valuable in the case of mixtures , where the _ ad hoc _ implementation of the equal tangent construction , required by most liquid state theories , may prove rather cumbersome .
the smooth cut - off formulation of hrt , proposed several years ago , has been quantitatively investigated only recently . among its peculiar properties ,
it is worth mentioning the absence of the spurious divergence of the isothermal compressibility at the phase boundary : a problem which affects the sharp cut - off formulation .
this allows the location of the spinodal loci within a renormalization group approach , as explictly shown for the yukawa fluid .
extensions of the smooth cut - off hrt equations to the case of mixtures or quantum systems are straightforward and represent a fertile research field for future studies .
the approximations to the hrt hierarchy we have been considering so far were limited to closures to the first equation , which provides the flow equation for the free energy of the model .
more sophisticated closures , at the level of the second equation , which contains information on the two particle correlations , have been formulated but not quantitatively investigated yet .
soft matter physics provides a stimulating environment for any microscopic theory able to capture the mechanisms leading to phase transitions . and
hrt was applied to several systems of physical interest in this context .
however , we feel that improved analytical representations of pair correlations in different frameworks have still to be developed to achieve a degree of accuracy comparable to simulation methods . for instance , a common feature of colloidal suspensions is the presence of short range attractive interactions , often represented by the  sticky limit " . a closure of the form ( [ rpa ] ) , adopted in most of the implementations of hrt , implies a linear relation between the direct correlation function and the attractive part of the potential , which is known to provide poor results for short range interactions @xcite .
extensions of hrt to lattice systems , mixtures and quantum models have been briefly covered in this review , but other interesting directions were explored in the literature , as the application of hrt to polymer physics or to disordered media @xcite .
a special mention must be devoted to the relationship between hrt and non perturbative ( or functional ) renormalization group methods .
the two approaches are in fact intimately related , as already noticed by caillol @xcite , and have been independently developed in different physical frameworks .
we hope that the present review , which specifically devoted section [ sec : nprg ] to this problem , will contribute to clarify this issue providing a dictionary between the two , substantially equivalent , formulations .
it is a pleasure to express our gratitude to all the former students and collaborators who helped us in these investigations , most notably d. pini , m. tau and f. lo verso .
we acknowledge financial support by the cimaina foundation ( lr ) and by miur through a prin/08 grant ( ap ) .
the flow equation for the free energy ( [ sharp ] ) , together with the definition ( [ calf ] ) and the parametrization of the direct correlation function ( [ rpa ] ) , supplemented by the compressibility sum rule ( [ komp ] ) form a closed partial differential equation ( pde ) for the free energy as a function of the density @xmath16 and the cut - off @xmath25 . as a practical example here
we report the explicit equations for a model defined by a spherically symmetric interaction written as the sum of a reference potential @xmath226 ( and corresponding free energy @xmath227 ) and an attractive tail @xmath129 . as usual
we set @xmath228 ; fourier transforms are denoted by tildes .
the initial condition of the pde is set at some large value of @xmath229 , such that the pair interaction @xmath230 can be neglected for @xmath231 .
the free energy at @xmath229 is then given by the mean field approximation : @xmath232 here and in the following , the ideal gas contribution is included in the definition of @xmath233 .
the integration of the pde is then performed backwards , down to @xmath29 , where all fluctuations are included and @xmath38 coincides with the free energy of the model .
a mesh is defined on the density axis @xmath234 . here
@xmath235 is some relatively high value ( typically @xmath236 in units of the hard - sphere diameter ) and the free energy @xmath237 is propagated from @xmath25 to @xmath238 by discretizing the pde on the mesh .
it is convenient to write the pde ( [ sharp ] ) in quasi - linear form by introducing the auxiliary quantity @xmath239 \label{defu}\ ] ] the evolution equation for this new variable can be straightforwardly obtained from eq .
( [ sharp ] ) : @xmath240 -\rho\,q^{d-1}\ , \frac{\tilde\phi(q)}{\tilde\phi(0)}\,\frac{d\,k_d}{2 } \,\frac{\partial^2 u_q}{\partial\rho^2 } \label{diff}\ ] ] where as usual @xmath241 in three dimensions ( @xmath154 ) and @xmath242 is the direct correlation function of the reference system .
usually the hard sphere gas is adopted as a reference system and @xmath242 is evaluated in percus - yevick @xcite or verlet - weis approximation @xcite .
this pde has the structure of a non - linear diffusion equation , where @xmath243 plays the role of time - like variable and @xmath16 of space - like variable .
the stability requirement for this class of pde s is given by the positiveness of the effective diffusion coefficient .
this condition is verified by eq .
( [ diff ] ) provided @xmath244 is positive .
boundary conditions at @xmath245 and @xmath246 are required : it follows from the definition that at low density @xmath247 , while @xmath248 is usually approximated by its mean field expression , i.e. by eq .
( [ defu ] ) with @xmath249 $ ] .
different boundary conditions at high density give very similar solutions to the pde , at least when the range of the potential is not too short .
details on the discretization procedure can be found in refs .
we finally remark that at low values of @xmath25 ( e.g. @xmath250 ) it is convenient to adopt a logarithmic scale in the cut - off , by setting @xmath251 with @xmath252 to slow down the inclusion of long wave - length fluctuations .
note that the diagrammatic resummation quoted in eq .
( 4.17 ) of ref .
@xcite is reported incorrectly . the expression
bond " should be substituted with  @xmath255 bond " throughout and the sentence
exactly one path " in item ( ii ) should read  at most one path " .
p.gianinetti , a.parola and l.reatto , in _ advances in quantum many - body theory vol .
5 : 150 years of quantum many - body theory _
bishop , k.a .
gernoth , n.r .
walet eds .
world scientific , singapore , ( 2001 ) .