the concept of solvent mediated interaction among colloidal particles provides an invaluable tool for the physical understanding of the behavior of complex fluids . when two hard are immersed in a molecular solvent , the simple excluded volume constraint is enough to cause an effective attraction between them driven by purely entropic effects , as shown in the seminal work by asakura and oosawa ( ao ) @xcite .
the depletion mechanism has been extensively studied , mostly in hard sphere fluids @xcite , and represents a very successful paradigm for the interpretation of colloidal aggregation . according to the ao prediction , which neglects correlations ,
the strength of attraction is proportional to the fluid density and its range coincides with the diameter of the solvent molecules .
however , both numerical simulations and liquid state theories proved that the effective interaction acquires a remarkable structure when the density of solvent molecules is increased , due to the presence of coordination shells in liquids .
attractive interactions between solvent molecules reduce particle accumulation near the walls giving rise to smoother solvent mediated forces @xcite . in the neighborhoods of the critical point of the fluid the critical casimir effect @xcite is expected to give rise to long range tails in the solvent mediated attraction between colloidal particles : the effective force between two bodies acquires the range of the correlation length of the underlying critical solvent and is described by a universal scaling law which depends on a very limited number of features ( e.g. the geometry of the bodies , the boundary conditions at the colloidal surfaces , the system dimensionality ) .
this remarkable phenomenon has been indeed experimentally observed by use of advanced optical techniques in a lutidine solvent close to its critical point @xcite .
numerical simulations have been performed for two representative models of solvent : a fluid and a fluid of hard spheres with three attractive patches . in both cases
the emergence of a long range attractive tail in the effective interaction between two large spheres immersed in the solvent is clearly visible @xcite .
an enhanced rate of aggregation among colloidal particles is expected to occur as a result of critical casimir effect , and indeed also this consequence has been directly observed experimentally @xcite .
these investigations suggest a novel way to control colloidal aggregation by tuning the solvent temperature .
critical casimir effect is due to the confinement of order parameter fluctuations induced by the presence of external boundaries . as such
, it governs the long range properties of the effective interaction , while at short range excluded volume effects are still present , providing the essential ingredient for the occurrence of the depletion phenomenon .
the interplay between these two quite different physical mechanisms which act at different has not been deeply investigated .
numerical simulations showed that , on approaching the critical point , the effective force indeed loses the short range oscillations due to the liquid layering near the surface of the bodies and acquires a smooth attractive form @xcite but no theoretical study of this effect has been attempted yet .
more generally , a thorough investigation of the form of solvent mediated interactions in the whole phase diagram of a correlated fluid has not been performed .
such a study would clarify the emergence of the long range tails in the effective interactions on approaching the critical region , as well as the coexistence of the depletion mechanism , acting at short distances , and the critical casimir effect , governing the physics at large distances .    in this paper
we tackle this problem by examining the solvent mediated forces between two hard walls immersed in a hard core yukawa fluid .
the theoretical investigation is performed in the region in order to exclude wetting phenomena , which deserve a separate analysis .
preliminarily , we must choose the most appropriate theoretical tool for such a study : a natural choice is density functional theory ( dft ) , which proved accurate in the description of confined fluids . however , no specific implementation of dft has been investigated in the critical region and most of the liquid state theories , which provide the theoretical basis for the practical formulation of dft , give a poor description of the critical regime .
in section ii we propose a weighted density approximation ( wda ) especially designed for being accurate also near the critical point .
this theory is validated against other dft prescriptions as well as available numerical simulations in several states . in section iii
this approach is applied to the critical region of a yukawa fluid : the density profiles and the solvent mediated force between two walls are determined by numerical minimization of the density functional .
the development of long range tails in both the density profile and the effective force is quantitatively investigated , confirming the trends already shown in simulations . in section
iv the universal properties of the effective interaction in the critical region are discussed .
the emergence of scaling laws is investigated in the critical and regime .
we also show that , within our approximate dft formulation , the critical casimir scaling function is deeply related to the universal bulk properties of the critical fluid .
a simple prescription for the theoretical evaluation of the critical and casimir scaling function is compared to the numerical results and to available simulations in ising systems .
section v offers some final comments ad perspectives .
the purpose of this work is the investigation of the solvent mediated interaction between two planar parallel walls immersed in a classical fluid in a wide range of temperatures and densities , including the critical region . in order to achieve this goal
we need an accurate description of the properties of a confined fluid .
the most successful theoretical approach to study inhomogeneous systems is density functional theory @xcite .
although alternative techniques , like integral equations or scaled particle theory , have been proposed @xcite , dft is generally considered to be the most accurate and versatile tool for dealing with inhomogeneous systems and has been applied in several frameworks : from the study of fluids in nanopores , to the structure of the liquid - vapor interface , to the theory of freezing ( see e.g. @xcite ) .      according to mermin
s extension @xcite of the hohenberg and kohn density functional theorem @xcite , the equilibrium density profile @xmath0 of a confined fluid can be found by minimizing a suitably defined density functional @xmath1 $ ] at fixed chemical potential @xmath2 and temperature @xmath3 .
@xmath1 $ ] can be conveniently expressed in terms of the external potential @xmath4 and the intrinsic free energy functional @xmath5 $ ] as : @xmath6=\mathcal{f}[n(\bm{r})]-\int \mathrm{d}\bm{r}\ , n(\bm{r } ) \bigl ( \mu -\phi(\bm{r})\bigr ) .
\label{eq : omega}\ ] ] the intrinsic free energy functional is exactly known only in the ideal gas limit @xmath7= \int \mathrm{d}\bm{r}\ , n(\bm{r})\bigl(\log \left({\lambda}^3 n(\bm{r})\right)-1 \bigr)\ ] ] ( here @xmath8 and @xmath9 is the thermal de broglie wavelength ) , while , in an interacting system , it is customary to separate this contribution , splitting @xmath5 $ ] as the sum of the ideal and the excess term @xmath10 $ ] : @xmath11=\mathcal{f}^{\mathrm{id } } [ n(\bm{r})]+\mathcal{f}^{\mathrm{ex}}[n(\bm{r})].\ ] ] several approximations for the excess part have been proposed over the years for dealing with specific problems .
hard sphere fluids are successfully described by rosenfeld s fundamental measure theory ( fmt ) @xcite also at bulk densities close to the solid transition .
even if the fmt is widely used and it its implementation is straightforward , at least for the planar geometry , rosenfeld s approximation for the excess free energy holds in principle only for fluids of purely hard particles of any shape .
attractive contributions in the interparticle potential are generally added as a mean field perturbation to the reference hard sphere excess free energy functional .
however , such an approach gives only qualitative predictions , particularly in the critical region , because it does not take into account correlations arising from the attractive tail of the potential . also the weighted density approximations ( wda ) by tarazona @xcite and curtin and ashcroft @xcite , whose prediction of the density profile of the hard sphere fluid are rather accurate , at least at moderate densities ,
includes the effects of an attractive tail in the potential only at the mean field level .
furthermore , approaches based on a truncation of the gradient expansion of the excess free energy functional , such as the square - gradient approximation , can describe only the slowly varying part of the density profile but fail in dealing with the short range structure close to the wall .    in order to describe the microscopic properties of the solvent mediated forces in the whole phase diagram of a fluid , including the liquid - vapor transition
, we need an approximation of the excess free energy functional which , in the uniform limit , provides an accurate description of the homogeneous fluid both in the dense and in the critical regime .
unfortunately such an implementation of dft has not been devised yet .
a novel wda - based approximation for the excess free energy functional based on the hierachical reference theory of fluids @xcite will be now introduced .    according to the weighted density approximation , the excess free energy functional
is expressed in terms of a weighted density @xmath12 as @xmath13=\int \mathrm{d}\bm{r}\ ,   n(\bm{r}){\psi}^{\mathrm{ex } } \bigl(\bar{n}(\bm{r } ) \bigr ) , \label{eq : wda_excfreeenergy}\ ] ] where @xmath14 is the excess free energy per particle of the homogeneous system evaluated at bulk density @xmath15 and the weighted density @xmath12 is written as a local average of the density profile , in terms of an isotropic weight function @xmath16 : @xmath17 the weight function can be generally dependent on the local value of an auxiliary reference density @xmath18 and has to satisfy the normalization requirement @xmath19 to ensure that in the homogeneous limit the weighted density coincides with the actual density of the fluid .    in order to get a reliable approximation for the excess free energy functional we have to choose properly the two key ingredients which characterize our wda ansatz , namely the homogeneous free energy @xmath20 and the weight function @xmath21 .
the only available microscopic liquid state theory which is able to account both for non critical and critical properties of a homogeneous fluid is the hierarchical reference theory ( hrt ) @xcite , which will be therefore adopted in this work for the evaluation of the excess free energy @xmath22 of the uniform fluid .
although the general formalism of hrt can be applied to fluids and mixtures with arbitrary pair interactions , quantitative results for specific models require the closure of the exact hrt equations by introducing some approximation .
a closure which proved remarkably accurate has been implemented in the case of a hard core yukawa ( hcy ) fluid @xcite , because the resulting hrt equations considerably simplify by use of the known solution of the ornstein - zernike equation available for this interaction .
the hcy potential is defined as the sum of a pure hard core term of diameter @xmath23 and an attractive yukawa tail of inverse range
@xmath24 : @xmath25 where the parameter @xmath26 , which defines the energy scale , is positive . in the following
we will investigate this model taking @xmath23 and @xmath27 as the units of length and temperature respectively .
having established the form of our excess functional in the homogeneous limit , we proceed with the explicit definition of @xmath10 $ ] for general density profiles .
we first focus our attention on the effects of the attractive part of the potential .
it is well known that its main contribution to the internal energy is given by the hartree term : this circumstance has been extensively recognized in the previous treatments , where attractive interactions were included just through such a contribution .
we are therefore led to isolate this term in the excess free energy , by writing : @xmath13=\mathcal{f}^{\mathrm{ex}}_{\mathrm{r}}[n(\bm{r})]+ \mathcal{f}^{\mathrm{ex}}_{\mathrm{h}}[n(\bm{r } ) ] , \label{eq : split}\ ] ] where the hartree contribution is given by @xmath28= \frac{1}{2}\int\mathrm{d}\bm{r}^\prime\int \mathrm{d}\bm{r}^{\prime\prime } \ ,
n(\bm{r}^\prime )   v_{\mathrm{y}}(|\bm{r}^\prime-\bm{r}^{\prime\prime}| ) n(\bm{r}^{\prime\prime}),\ ] ] whereas the reference term @xmath29 $ ] , defined by eq .
( [ eq : split ] ) , contains both the entropic contribution to the free energy , arising from hard core repulsion , and the correlations induced by the attractive interaction .
our choice is then to use wda to represent only the part of the intrinsic free energy functional , retaining the exact description of the hartree energy : @xmath30=\int \mathrm{d}\bm{r}\ ,   n(\bm{r}){\psi}^{\mathrm{ex}}_{\mathrm{r}}\bigl(\bar{n}(\bm{r } ) \bigr),\ ] ] where @xmath31 the form of the weight function can be determined following the strategy put forward by tarazona @xcite , by requiring that the two particle direct correlation function reduces , in the homogeneous limit , to the known form of the underlying bulk liquid state theory which , in our case , is the direct correlation function @xmath32 predicted by hrt : @xmath33}{\delta n(\bm{r } )   \delta n(\bm{r}^\prime)}\right|_{n(\bm{r})=\rho_b }   \hskip -5 mm & = & c(|\bm{r}-\bm{r}^\prime|,{\rho}_b)+ \beta \ , v_{\mathrm{y}}(|\bm{r}-\bm{r}^\prime| ) \nonumber \\ & \equiv & c_{\mathrm{r}}(|\bm{r}-\bm{r}^\prime|,{\rho}_b ) , \label{eq_constraint_wda_r}\end{aligned}\ ] ] where the last equality is a definition .
this constraint can be fulfilled only by a density dependent weight function . in a previous work curtin and
ashcroft @xcite proposed that the weight function should depend on the local value of the weighted density itself , i.e. @xmath34 . unfortunately , even for a purely repulsive hard sphere system , this hypothesis leads to difficulties in the solution of the nonlinear differential equation for the weight function which must be imposed in order to implement the constraint ( [ eq_constraint_wda_r ] ) .
this problem can be overcome choosing a position independent auxiliary density @xmath18 , as suggested by leidl and wagner @xcite . in our application of dft to the description of confined systems , we assume that the auxiliary density actually coincides with the homogeneous density at the same chemical potential ( also referred to as _ bulk density _ ) , i.e. we set @xmath35 it is straightforward to obtain from eq .
( [ eq_constraint_wda_r ] ) an algebraic equation for the fourier transform of the weight function : @xmath36 which , at least in the cases examined in this work , always admits real solutions .
the physical root can be determined enforcing the normalization condition ( [ eq : normalization_w ] ) , recalling that the compressibility sum rule , satisfied by the hrt direct correlation function , requires @xmath37    , @xmath38 ) is very close to the critical point .
, width=264 ]    as in other implementations of wda , the range of the weight function always remains comparable to the size of the molecules .
this feature is preserved also in the critical region , as shown in fig .
[ fig : weightfunction ] .
once the intrinsic density functional has been determined , the approximate equilibrium properties , such as the density profile @xmath0 and the grand canonical potential @xmath39 $ ] , can be obtained minimizing the functional ( [ eq : omega ] ) at fixed temperature and chemical potential . in a hcy fluid ,
the value of the chemical potential is related to the bulk density by : @xmath40    a central quantity for the present investigation is the force acting on the two planar hard walls in the hcy fluid . in this geometry , symmetry requires that all the local properties may depend on the single coordinate @xmath41 , orthogonal to the two plates , placed at @xmath42 and @xmath43 respectively .
remarkably , if the wall separation @xmath44 is greater than @xmath23 , the force per unit surface @xmath45 acting on the plates , sometimes called _ solvation force _ , can be expressed as a pressure difference @xcite : @xmath46}{\sigma } \right)\right|_{t,\mu}-p_b(t,\mu ) , \label{eq : solvation_derivata}\ ] ] where @xmath47/\sigma$ ] is the grand canonical potential per unit surface of the fluid confined in the region @xmath48 $ ] , determined by the minimization of the approximated grand canonical functional at fixed @xmath49 and @xmath50 is the pressure of the fluid at the same values of temperature @xmath3 and chemical potential @xmath51 . on the other hand , when @xmath52 there are no particles between the walls and the attractive force per unit surface acting on the walls arises uniquely from the presence of the fluid in the regions @xmath53 and @xmath54 : the first term in eq .
( [ eq : solvation_derivata ] ) vanishes and the force is given by @xmath55 the solvation force @xmath45 defined above is a difference of pressures and goes to zero in the limit @xmath56 . by means of standard functional identities it is possible to express , without any further approximation ,
the derivative of the grand potential per unit surface in terms of the contact density : @xmath57 and the solvation force can be finally written as @xcite : @xmath58 we also remark that this version of wda exactly satisfies the contact value theorem @xcite @xmath59 leading to the more suggestive identity @xmath60    the above results show that , in order to obtain the force acting on the walls we just need to perform the minimization of the functional in the region @xmath48 $ ] in order to evaluate the contact density @xmath61 .
the minimization has been be carried out by a simple iterative ( picard ) method , taking advantage of the exact implicit relation for the approximate equilibrium density profile @xmath62 given by : @xmath63 } ,   \label{eq : iteration_minimization}\ ] ] where the potential of mean force @xmath64 is defined as : @xmath65}{\delta n(\bm{r})}\bigg|_{\rho(z ) } + \int \mathrm{d}\bm{r}^\prime \rho(z^\prime)v_{\mathrm{y}}(|\bm{r}-\bm{r}^\prime|).\ ] ] a spatial @xmath66 has been generally used in the numerical minimization , while for bulk densities @xmath67 and close to the critical point @xmath68 has been reduced up to two orders of magnitude .
typically , up to few thousand picard iterations were necessary to achieve a precision of one part in @xmath69 for the density profile and one part in @xmath70 in the grand potential .
the minimization of the previously defined grand canonical functional allows to evaluate the equilibrium properties of the confined fluid . within this approach
, the relevant quantities can be found at every temperature and bulk density of interest , also in the vicinity and below the critical point of the hcy fluid .
most of the calculations refer to a yukawa fluid with range @xmath71 , where several simulation results are available .
we performed minimizations of the functional for values of the temperature @xmath3 above the critical point ( @xmath72 ) and bulk reduced densities @xmath73 up to @xmath74 .
in the high temperature limit our model reduces to a hard sphere fluid , whose properties have been extensively investigated by numerical simulations .
figure [ fig : hs_simulation_fmt ] shows the density profiles @xmath62 of a hard sphere fluid near a hard wall at two different values of @xmath75 .
the agreement of the wda prediction with the monte carlo ( mc ) data of ref .
@xcite is very good up to reduced densities of the order of @xmath76 , while at higher values the phase of oscillations in the density profile are correctly captured , although a slight underestimation of the peak value is observed .
the comparison of our density profiles with those predicted by the `` white bear '' version of the fmt @xcite shows , as expected , that rosenfeld s theory gives more accurate estimates of the oscillation peaks , particularly at high density .     for a hard sphere fluid between two hard walls ( distance @xmath77 ) at different values of the bulk reduced density @xmath73 .
mc simulation data ( symbols ) are taken from ref .
we have obtained the fmt density profiles ( lines ) from the minimization of the approximated fmt - wb functional @xcite . to enhance visual clarity the density profiles at @xmath78 and @xmath79
are shifted upward by @xmath80 and @xmath81 respectively.,width=264 ]    when the temperature is decreased the contribution of the yukawa tail to the density profile becomes relevant .
we compared the density profile obtained within our dft approximation with the mc simulation data for @xmath71 at temperature @xmath82 and for @xmath83 at @xmath84 .
figure [ fig : t2_diversirho ] shows that the wda estimate is remarkably accurate at reduces densities @xmath80 and @xmath85 ; small deviations from the mc simulation data appear at reduced density @xmath86 .
we note that at @xmath87 the contact reduced density is overestimated of about @xmath88 with respect to simulation data , even if the contact theorem is verified with a relative error of the order of @xmath89 .
this disagreement in the contact value is due to different estimates of the grand canonical potential per unit volume of the homogeneous fluid and it is compatible with the the spread in the values of the bulk pressure obtained within different simulation techniques @xcite .     of a yukawa hard sphere fluid ( @xmath90 ) at reduced temperature @xmath82 confined between two hard walls at different values of bulk reduced density @xmath73 .
lines represent the predictions of the present wda .
points are mc data from ref .
@xcite ( reduced density @xmath86 ) and ref .
@xcite ( reduced density @xmath80 and @xmath85 ) .
the distance between the walls is @xmath91.,width=264 ]    in panel a ) of figure [ fig : profili_forze_louis ] we compare our results for the density profile of the hcy fluid characterized by @xmath83 with the recent mc simulations from ref .
this figure shows that at relatively low densities the agreement between our approximation and mc predictions is remarkable also for attractive potentials of shorter range .
in particular we predict accurately the kink in the density profile at @xmath92 , which is only qualitatively reproduced within mean field approximation @xcite .     of a hcy fluid at a single hard wall ( @xmath83 ) for bulk reduced density @xmath93 and temperature @xmath84 .
lines represent the predictions of the present wda .
points are mc data from ref .
panel b ) : force per unit surface acting between two infinite parallel hard walls immersed in a hcy fluid ( @xmath83 ) in the same thermodynamic state of panel a ) .
lines are the predictions of the wda from eq .
( [ eq : solvation_rho ] ) , whereas points represent the md simulation of ref .
@xcite.,width=264 ]    the minimization of the grand canonical functional provides both the value of the contact density @xmath94 and of the grand free energy .
it is therefore possible to obtain the depletion force either by calculating the derivative w.r.t .
@xmath95 of the grand potential , as in eq .
( [ eq : solvation_derivata ] ) , or by making use of eq .
( [ eq : solvation_rho ] ) .
the consistency between the two estimates is a good check for the accuracy of the numerical procedure . to give an example , for a hard sphere fluid
, the relative difference between the two results is less than @xmath96 if the absolute value of the force per unit surface and @xmath97 is larger than @xmath98 . nonetheless at smaller values of the force
the result obtained by differentiation of the grand free energy is less stable , due to errors introduced by the discretization .
the solvent mediated force acting on two parallel hard walls immersed in a hard sphere fluid , is compared with the monte carlo data of wertheim et al .
@xcite as well as the predictions based on fmt in fig .
[ fig : forza_hs_piani ] , showing a nice agreement also at relatively high densities . at reduced density @xmath99
the force maximum per unit surface is of the order of @xmath100 , and the oscillations due to the packing of the hard spheres are damped within two or three diameters .
furthermore , at this value of the reduced density , a small deviation between wda and fmt is present only at the first minimum . at reduced density @xmath76 the force at distances of the order of the hard sphere diameter @xmath23 is hundred times larger than at @xmath101 .
even this feature is well reproduced by both wda and fmt , as can be seen in the inset .
the strong oscillating behavior of the mc data is captured by wda with a correct phase , even if the peak values are a little underestimated , whereas fmt behaves considerably better .     between two infinite planar hard walls immersed in a fluid of hard spheres of diameter @xmath23 for two values of the bulk reduced density .
the value of the force within the wda approximation has been obtained via eq .
( [ eq : solvation_rho ] ) .
the fmt - wb result comes from the minimization of the `` white bear '' version of the fmt functional @xcite via eq .
( [ eq : solvation_rho ] ) .
the mc data is taken from ref .
the inset highlights the behavior of the force at @xmath102 at small distances.,width=264 ]    the solvent mediated force per unit area between two walls in a hcy fluid has not been extensively investigated by numerical simulations . in panel b ) of fig .
[ fig : profili_forze_louis ] we show a comparison between our results and the monte carlo data of ref .
@xcite for a hcy fluid of inverse range @xmath83 at @xmath84 and @xmath103 .
even if the net force is quite small , our prediction agrees very well with the numerical simulations at all values of the wall separation .
we stress that , particularly at small distances , the wda force is much more accurate than any mean field perturbation method ( see ref .
@xcite ) .
the detailed comparisons of our novel dft with both numerical simulations and state of the art theories allow to conclude that in the high temperature limit the present wda is able to correctly reproduce the density profile and the effective interactions between the hard walls with a very satisfactory accuracy up to reduced densities of about @xmath85 .
moreover , this formulation of wda appears to be the best available dft for a hcy fluid at finite temperature .
we performed the minimization of the wda density functional at several values of temperatures and reduced bulk densities for a hcy fluid of inverse range @xmath71 confined between two hard walls . at high temperatures the system behaves like a hard sphere fluid , whereas when the temperature is decreased the contribution of the yukawa tail becomes more and more relevant , and the strongly oscillating character of both the density profiles and the solvent mediated force is lost .
figure [ fig : profili_05 ] shows the dependence of the density profile on temperature at fixed bulk reduced density @xmath104 . at reduced temperature @xmath105
the system behaves like a hard sphere fluid .
as the temperature is lowered towards its critical value , the density profile gradually becomes monotonic losing the oscillating features typical of hard spheres and the density at contact assumes values four times lower than the bulk density .
the range of the perturbation produced by the wall extends at larger and larger distances as the temperature approaches @xmath106 , giving rise to a region where a kind of drying of the wall can be observed .    ) at a single hard wall at bulk reduced density @xmath107 and different values of the temperature .
the inset shows the same data , also including the density profile at temperature @xmath105 using a different scale .
, width=264 ]    the attractive tail in the pair interaction of the hcy fluid smoothes the density profile reducing the layering of particles . as a consequence ,
the effective force between the two walls loses the strongly repulsive peak present at @xmath108 when the interaction between the fluid particles is purely hard sphere ( see fig .
[ fig : forza_hs_piani ] ) . fig .
[ fig : forze_rhocritico ] shows the force per unit surface between the walls for different values of the temperature at the critical bulk reduced density @xmath109 . the repulsive contribution to the interaction force
, present at @xmath105 gradually disappears at lower temperatures and the force becomes purely attractive and monotonic , confirming the findings of the numerical simulations in a different model @xcite . by approaching the critical temperature ( @xmath110 )
the effective force becomes weaker and weaker at short distance , as can be seen in the right panel of fig .
[ fig : forze_rhocritico ] : its amplitude reduces almost by a factor two due to a @xmath111 change in temperature
. however , a closer look to the long distance tail of the solvent mediated force shows that its range indeed increases close to the critical temperatures , as expected on the basis of scaling arguments .
however this occurs at very large separations ( @xmath112 for the data shown in the figure ) .    )
obtained from the minimization of the wda functional using eq .
( [ eq : solvation_rho ] ) at @xmath113 and different values of the temperature.,width=264 ]    fig .
[ fig : forze_tcritico ] shows the force per unit surface between the walls in different density regimes when the value of the temperature is close to @xmath106 .
we note that when the bulk density is higher than @xmath114 the force is an order of magnitude larger than for @xmath115 .
the force is monotonic and purely attractive for all values of the reduced density and its range grows near @xmath114 , as expected .    )
obtained from the minimization of the wda functional using eq .
( [ eq : solvation_rho ] ) at @xmath116 and different values of the reduced density.,width=264 ]    in lattice fluid models the coexistence curve is symmetric about the critical temperature and the critical isochore coincides with the locus @xmath117 of the maxima of the isothermal susceptibility . for such systems , according to the renormalization group terminology , the path to the critical point orthogonal to the relevant odd operator coincides with the critical isochore . in a real fluid
the coexistence curve is asymmetric about the critical isochore . in this case
a good approximation for the same path is given by the line @xmath118 in the phase diagram defined as the locus of the points @xmath119 such as : @xmath120 where @xmath121 is the isothermal compressibility .
[ fig : fit_profili_rhocritico ] shows the density profile of the hcy fluid at a hard wall along the line @xmath118 .
its behavior at distances larger than the bulk correlation length is well fitted by an exponential of the form : @xmath122 where @xmath75 is the bulk density of the fluid , @xmath123 is a negative amplitude factor and @xmath124 is the bulk correlation length .
this exponential decay of the density profile is observed also if the system is far from the critical region and is probably related to the location in the bulk phase diagram of the point we investigated with respect to the line @xcite . following the argument of ref .
@xcite , we expect that the exponential long range behavior of the density profile reflects in an analogous exponential decay of the force between the two walls .
provided we do not cross the line , this decay should be present both far from the critical point and in the critical region , where it agrees with the predictions of the theory of the critical casimir effect @xcite : @xmath125 the exponential decay of the solvation force is indeed confirmed by the exact solution of a @xcite ising slab under symmetry breaking boundary conditions and was observed in monte carlo simulations of three dimensional simple fluids @xcite .
[ fig : fit_forze_rhocritico ] shows the long distance exponential decay of the force per unit surface and @xmath97 between two planar hard walls mediated by a hcy fluid along the critical line @xmath118 .
the force obtained with our approach is very well fitted by eq .
( [ eq : fit_forza ] ) at large distances ( @xmath126 ) , whereas at short distances , where depletion effects become relevant , the solvent mediated force , always attractive , displays a plateau ( see inset ) .    ) at a single hard wall in the critical region . here
points : fits of the density profiles according to eq .
( [ eq : fit_profilo ] ) performed at distances larger than three times the correlation length of the homogeneous hcy fluid at the same temperature and bulk density .
@xmath75 , @xmath123 and @xmath124 are free fitting parameters and the results obtained for @xmath75 and @xmath124 agree well with the bulk values of density and correlation length respectively ( the accuracy is better than @xmath128 for the bulk density and @xmath129 for the correlation lenght ) .
the bulk correlation lengths obtained from the fitting procedure are @xmath130 , from the lowest to the highest reduced temperature .
the inset shows a magnification of the same density profiles at short distances . ,
width=264 ]    ) at a single hard wall along the line @xmath118 . here
points : fits of the force according to eq .
( [ eq : fit_forza ] ) performed at distances larger than four times the bulk correlation length of the homogeneous hcy fluid at the same temperature and bulk reduced density .
@xmath131 , and @xmath124 are free parameters in the fitting procedure .
the result obtained for @xmath124 agrees well with the bulk value of the correlation length ( the accuracy of the fit is better than @xmath132 ) .
the bulk correlation lengths obtained from the fitting procedure are @xmath133 , from the lowest to the highest reduced temperature .
the inset shows the force profile at short distance.,width=264 ]      the same wda formalism previously introduced may be generalized to other interesting geometries , in addition to the planar one .
most importantly , it can be used to evaluate the effective interaction between two spherical particles in a solvent , with obvious applications to the study of aggregation in colloidal suspensions .
the direct minimization of the dft , although numerically feasible , represents a task considerably more complex than in planar geometry .
therefore , in this first application of the formalism , we have chosen to resort to the simple but effective derjaguin approximation @xcite , which allows to express the interaction between two convex objects starting from the knowledge of the force between two planar walls , independently on the physical origin of the force . according to this approximation
, the force @xmath134 between two spheres of radius @xmath135 can be written as : @xmath136 where @xmath95 is the minimal surface to surface distance between the sphere and @xmath137 is the force per unit surface between the two walls at distance @xmath41 .
this approximation gives accurate results provided @xmath138 and if the interaction potential between the two walls decays rapidly at large distances .
when the force between the walls is mediated by a hard sphere fluid with particles of diameter @xmath23 it is possible to show that derjaguin s expression is the best approximation of the true depletion interaction without taking in account curvature effects @xcite and it is accurate in the limit of @xmath139 , where @xmath140 is the size ratio .
the limits of derjaguin approximation when applied to depletion interactions is a debated issue .
particularly , it is a matter of discussion the size ratio at which derjaguin s theory starts to fail , and how its accuracy depends on the concentration of depletant . in figure
[ fig : depletion_sferedure ] we compare the prediction for the depletion potential @xmath141 between two big hard spheres in a fluid of smaller hard spheres obtained both by derjaguin approximation and mc simulations at two different values of the bulk density of the smaller particles .
the predictions for the depletion potential obtained by derjaguin approximation are accurate at @xmath142 only at values of the packing fraction of the small spheres lower than @xmath143 . at @xmath144 derjaguin approximation overestimates by about @xmath145 both the contact value and the first repulsive peak in the potential , while the oscillations at larger values of distance are underestimated when compared to mc data .
the rather poor performance of derjaguin approximation is probably due to the presence of a strong repulsive peak in the solvation force between the two walls at @xmath146 ( see fig .
[ fig : forza_hs_piani ] ) which appears to be a peculiarity of the slab geometry .
we expect that , for smoother inter - wall effective interactions , the agreement would be considerably better .
between two hard spheres in a fluid of smaller hard spheres ( size ratio @xmath142 ) at bulk packing fraction @xmath147 equal to @xmath88 ( left panel ) and @xmath148 ( right panel ) .
mc ( points ) and fmt predictions obtained with the insertion trick ( purple lines ) are taken from ref .
the red line at @xmath149 is obtained by use of derjaguin approximation starting from the solvation force between two planar hard walls evaluated from the fmt - wb @xcite approximation .
the green lines represent the depletion force obtained within derjaguin approximation when the solvation force between the walls is given by the present wda approximation.,width=264 ]    the results presented above show that derjaguin approximation can be safely adopted is the size ratio between the depletant and the colloid is sufficiently small and up to moderate densities of depletant ( i.e. @xmath150 and @xmath151 ) , and we expect that similar considerations apply when the depletant is a yukawa hard core fluid . within the limits of validity of the derjaguin approximation , we can determine the solvent mediated potential @xmath152 between two hard spheres of radius @xmath135 immersed in a yhc fluid in order to examine the phase stability of such a colloidal suspension . according to noro - frenkel extended law of corresponding states @xcite , fluids characterized by short ranged interaction potentials obey the same equation of state , when expressed in terms of reduced variables .
in particular , it was observed that the dimensionless second virial coefficient @xmath153 where @xmath154 is the second virial coefficient of a hard sphere system with particles of radius @xmath135 , assumes a value of about @xmath155 at the critical point independently on the particular form of the interaction and that its value remains constant in a relatively large density range across the critical point @xcite .
it is therefore possible to estimate the the gas - liquid spinodal line for a system of hard sphere colloidal particles dispersed in a hcy fluid , by evaluating their reduced second virial coefficient .    ) in the @xmath156 plane .
the hrt results for the coexistence curve @xcite are shown by a black line .
green dots represent the mc data from ref .
blue points show the aggregation boundary of two big hard spheres ( size ratio @xmath157 ) predicted on the basis of our wda plus derjaguin approximation and noro - frenkel criterion .
the dashed line connecting the points is a guide to the eye : aggregation takes place on the right of this boundary .
the inset shows the same phase diagram in a wider density interval.,width=264 ]    the blue points in figure [ fig : aggregation ] identify the phase separation line of a hs fluid induced by the presence of a depletant modeled as a hcy fluid .
the size ratio between the depletant and the guest hs particles is @xmath142 . at high values of the reduced temperature
the phase separation occurs at reduced density of about @xmath80 , as expected in the limit of hs depletant .
when the temperature decreases , the concentration of depletant needed to induce phase separation decreases monotonically . at the depletant critical temperature
we observe phase separation at depletant concentrations @xmath158 much lower than the critical one .
this implies that , at this value of the size ratio @xmath159 , the phase separation is not related to the presence of long range tails in the effective force , which characterizes the critical region of the solvent , but is still mainly due to the short range attraction generated by the depletion mechanism .
we note that whenever a direct short range repulsion is present between the colloidal particles , as for the case of charged systems , the strong attraction due to depletion is severely weakened and particle aggregation takes place at considerably larger solvent densities . instead , in the critical region , the long range tails of the solvent mediated ( casimir ) force is not effectively contrasted by the additional short range repulsion . in extreme circumstances ( i.e. when the direct repulsion between particles is sufficiently strong ) , ordinary depletion may be fully screened and phase separation inhibited .
however , aggregation is generally expected in a small pocket within the critical region , due to the emergence of long range casimir forces . for repulsive wall - solvent interactions ,
this pocket will be centered at solvent densities larger than the critical one , due to the strong asymmetry of the critical casimir forces ( see e.g. figs . [
fig : forze_tcritico ] and [ fig : scala_micro_tcritica ] ) .
the aim of this section is to evaluate the solvent mediated interaction induced between two walls when the depletant is in the critical regime .
thermal fluctuations in fluids occur on a characteristic range , determined by the correlation length @xmath124 , which usually is comparable with the molecular diameter or with the range of the interactions . when approaching a second order phase transition , the range of the fluctuations of the order parameter ( the particle density in the case of a simple fluid ) grows larger up to diverging at the critical point . in this regime
the effective force between two bodies immersed in the critical fluid acquires a universal form and obeys scaling laws , as many physical properties near criticality do . in 1978
@xcite fisher and de gennes first recognized that a confinement of the critical fluctuations of the order parameter gives rise to an universal , long ranged interaction which they named critical casimir force . according to the scaling approach , the universal contribution to the force per unit surface @xmath160 acting between two infinite plates confining a critical fluid can be written as @xcite : @xmath161 where @xmath95 is the distance between the two walls .
the upper sign refers to the temperature ( while the lower to the one ) and the scaling variables ( @xmath162 ) @xmath163 are defined in terms of the two scaling fields @xmath164 here @xmath165 is the bulk correlation length at @xmath166 , @xmath167 is a non - universal metric factor and @xmath168 , @xmath169 and @xmath170 are the usual critical exponents .
the function @xmath171 is usually referred to as the scaling function of the critical casimir force in planar geometry .
this function is universal in sense that it depends only on the bulk universality class , on the boundary conditions imposed at the confining surfaces and on the geometry of the system ( which in this case is @xmath172 ) .
we remark that there is no extra metric factor associated with @xmath173 and that there is a dependence on the sign of the field @xmath44 because boundary conditions at the walls break the bulk symmetry @xmath174 .
according to the theory of scaling , eq .
( [ eq : casimir_film ] ) represents the asymptotic decay of the solvent mediated force as @xmath175 and for @xmath176 .
the bulk yukawa fluid under investigation belongs to the 3d ising universality class and the boundary conditions are determined by the affinity of the wall surfaces with the fluid particles : if the contact density is less than the bulk density the boundary condition is of type @xmath177 , otherwise of type @xmath178 . in this work we only deal with temperatures ( @xmath179 ) and with symmetric @xmath180 boundary conditions , which arise for purely repulsive interactions between the fluid particles and two identical confining hard walls .
the analysis of the critical casimir force and of the related universal scaling function @xmath181 is a rather difficult task both experimentally , for the small forces involved , and theoretically , for the lack of an accurate description of critical fluids in confined geometries .
most of the results present in the literature deal with the temperature dependence of the casimir fluctuation induced interaction at zero magnetic field .
an indirect estimate of the scaling function for the film geometry and the 3d ising universality class , at @xmath166 under @xmath182 and @xmath183 boundary conditions was given in ref .
@xcite monitoring the thickness of a binary fluid film at different temperatures near @xmath106 .
the first direct evaluation of the critical casimir force was performed in 2008 @xcite for a system consisting of a colloidal particle close to a wall immersed in a binary mixture ( geometry ) of water and lutidine at different compositions and for both the relevant boundary conditions .
this experiment , however , allowed to probe only the exponential tail of the scaling function .
a mc study of the solvent mediated potential between two spherical particles in a simple fluid along the critical isochore has been performed in refs .
@xcite with different boundary conditions .
however , the determination of the full casimir scaling function could not be obtained in the temperature range examined in the simulations . along the symmetry line ( @xmath166 )
, more precise estimates of the universal casimir scaling functions for the 3d ising universality class and film geometry have been obtained via mc simulations of the ising model @xcite .
few theoretical approaches were devised to address this problem : in addition to the mean field results @xcite , its is worthwhile mentioning the extended local functional method @xcite and a long wavelength analysis of density functional theory @xcite .
the latter investigations have been also extended away from the symmetry line ( @xmath184 ) providing predictions on the shape of the critical casimir scaling function in the case @xcite .
monte carlo simulations at @xmath185 were recently performed in ref .
@xcite .
the wda approach developed above allows the study of this problem starting from the microscopic hcy fluid model confined between two walls . according to the scaling hypothesis , the effective force per unit surface between the two walls
@xmath160 should depend on the physical control parameters @xmath186 only through the combination ( [ eq : casimir_film ] ) , implying the collapse of different data sets onto the same universal curve .
figure [ fig : scala_micro_rhocritica ] shows the scaling function obtained from independent calculations at different temperatures along the previously defined critical line @xmath187 ( i.e. @xmath188 ) .
note that , even at reduced temperature @xmath127 as low as @xmath189 , our estimates show a marked temperature dependence , and the data along different isotherms do not collapse as we expected . at the lowest temperature we investigated , a significant difference between our prediction and the mc simulations of ref .
@xcite suggests the presence of strong corrections to scaling .
we also remark that the curves at the lowest temperatures develop a kink at small values of @xmath173 , due to the singular behavior of the scaling function at @xmath190 .
in fact , at any given reduced temperature @xmath191 , the quantity @xmath192 tends to zero as @xmath193 , forcing the estimate of the scaling function to vanish .    ) from the microscopic force obtained within the present wda approximation .
the corresponding bulk correlation lengths are @xmath194 , from the lowest to the highest reduced temperature .
the lines connecting the points are a guide to the eye .
the mc data are taken from ref . @xcite and refer to @xmath180 boundary conditions and the different sets correspond to different estimates of the corrections to scaling.,width=264 ]    figure [ fig : scala_micro_tcritica ] shows the scaling function at fixed temperature near @xmath106 for different values of the scaling variable @xmath195 , corresponding to different bulk reduced densities .
the scaling function is always negative and a strong asymmetry is evident between the curves at density above and below @xmath114 . for positive values of the scaling field @xmath44
( i.e. @xmath196 in our case ) , the magnitude of the force becomes larger and larger and the peak is shifted towards small values of @xmath173 .     and for different values of the scaling variable @xmath197 .
the lines connecting the points ar a guide to the eye .
the largest correlation length , @xmath198 , corresponds to the critical reduced density @xmath113.,width=264 ]      although the direct numerical evaluation of the critical casimir scaling function predicted by this class of dft is not conclusive , due to severe corrections to scaling , an accurate estimate of the asymptotic behavior can be obtained by a long wavelength ( lw ) analysis of the dft equations .
in fact , following ref .
@xcite we note that :    * the density profile @xmath62 displays reflection symmetry about @xmath199 , limiting the range of interesto to @xmath200 $ ] ; * when the walls are far apart ( @xmath201 ) , the difference between the density profile corresponding to a wall to wall distance @xmath95 and its single wall limit , reached for @xmath202 , is significant only for @xmath203 ; * as a consequence , the effective force per unit surface @xmath204 ( hence the casimir scaling function ) just depends on the long distance tail of the density profile , which is expected to be a slowly varying function of the coordinate @xmath41 .
therefore our wda intrinsic free energy functional can be approximated by keeping only the lowest term in a gradient expansion about the bulk density @xmath75 : @xmath205}{\sigma}=l\,\varphi(\rho_b ) + \nonumber \\ & & \,\int { \rm d } z \ , \bigg [ \frac{b}{2 } \left ( \frac{{\rm d}\delta n(z ) } { { \rm d } z}\right)^2 +   \varphi(\rho_b+\delta n(z ) ) -\varphi(\rho_b ) \bigg ] ,   \label{eq : longw}\end{aligned}\ ] ] where @xmath206 is the free energy density in the bulk , times @xmath169 .
this expression coincides with the long wavelength limit of our wda functional , the stiffness @xmath207 being related to the range of the direct correlation function in the homogeneous system @xmath208 : @xmath209 in the presence of short range interactions , the direct correlation function is analytic in @xmath210 away from the critical point , where it displays a @xmath211 singularity .
however , within our approximate closure of the hrt equations , the critical exponent @xmath212 and analyticity is preserved also at criticality @xcite , keeping the stiffness @xmath207 finite in the whole phase diagram .
this implies that the long wavelength limit of the structure factor of the homogeneous fluid follows the orstein - zernike ansatz : @xmath213 with @xmath214^{-1}\ ] ] and @xmath215 .
close to the critical point , the hrt bulk free energy density @xmath206 acquires a scaling form : @xmath216 where @xmath114 is the critical density and @xmath217 is the chemical potential ( the temperature dependence of these quantities is understood ) , while @xmath218 and @xmath219 are non universal metric factors and @xmath127 is the reduced temperature . in the following
it will be convenient to express the universal quantities in terms of the scaling field @xmath220 instead of the previously defined variable @xmath197 . within our hrt closure ,
the critical exponents are @xmath221 , @xmath222 , @xmath223 in @xmath224 , which agree within @xmath111 with the accepted values .
the metric factors appearing in the scaling function are implicitly defined by the requirement that @xmath225 has the following expansion at small @xmath226 @xcite : @xmath227 in fig .
[ fig : scalapsi_g4 ] the asymptotic hrt scaling function @xmath225 is shown together with a parameterization of the exact result for the 3d ising universality class .
although the two curves are indistinguishable on this scale , calculations at different reduced temperatures , also shown , suggest the presence of important corrections to scaling .     as predicted by hrt
compared with a parameterization of the exact result from ref .
@xcite ( black curve ) .
the red line , showing the hrt asymptotic result , is identical to the exact result on the scale of the figure .
a few rescaled free energies obtained from the integration of the hrt equations at different reduced temperatures are also shown .
inset : universal amplitude ratio @xmath228 for a hcy fluid at the critical density as a function of the reduced inverse correlation length ( green points ) .
the usually accepted value @xcite is shown by a blue triangle whereas the red square represents the asymptotic hrt value . , width=264 ]    the minimization of the long wavelength functional ( [ eq : longw ] ) in slab geometry gives rise to a differential equation whose solution allows to evaluate the asymptotic decay of the effective force between two hard walls in a critical fluid .
the derivation , already detailed in ref .
@xcite and not repeated here , provides a closed form for the critical casimir scaling function in terms of two universal quantities : the bulk free energy scaling function @xmath225 and the universal amplitude ratio @xmath228 . defining the auxiliary quantity @xmath229 by the implicit relations : @xmath230 the critical casimir scaling function in three dimensions is given by : @xmath231
the universal amplitude ratio @xmath228 is expressed in terms of the non universal metric factors previously introduced as @xmath232 again , the evaluation of @xmath228 from the hrt equations displays severe correction to scaling in a hcy fluid , as shown in the inset of fig . [
fig : scalapsi_g4 ] .
more importantly , the usually quoted `` exact '' value @xcite @xmath233 turns out to differ significantly from the hrt prediction @xmath234 .
the asymptotic study of the dft equations allows to extract the critical casimir scaling function just from bulk quantities via eqs .
( [ sigma1]-[eq : scalingc ] ) .
it is then instructive to contrast these predictions with the outcome of the direct minimization of the hrt functional , already shown in fig .
[ fig : scala_micro_rhocritica ] .
such a comparison can be found in fig .
[ fig : micro_vs_long ] , where the scaling functions obtained from the microscopic dft at a few reduced temperatures @xmath235 in the critical region are shown to agree remarkably well with the predictions of the long wavelength analysis , provided both the scaling function for the free energy @xmath225 and the universal amplitude ratio @xmath228 are consistently evaluated at the same reduced temperature @xmath235 .
and @xmath228 evaluated at the same reduced temperature as the dft calculation .
black line : asymptotic limit of the critical casimir scaling function .
points : prediction of @xmath181 from monte carlo simulations @xcite .
, width=264 ]    the presence of strong corrections to scaling in both @xmath225 and @xmath228 , already highlighted , induces strong effects in the critical casimir scaling function which , at reduced temperatures lower than @xmath189 , is still quite far from its asymptotic limit .
the main effect is due to the growth of the amplitude ratio , which , as shown in fig .
[ fig : scalapsi_g4 ] , appears to reach its universal value only extremely close to the critical point , according to the prediction of hrt for the model of critical fluid investigated here .
we remark that , due to the already quoted difference between the hrt estimate of the universal amplitude ratio @xmath228 and the value obtained via series expansions , the critical casimir scaling function predicted by our dft significantly differs from the one obtained in mc simulations , as can be seen in fig .
[ fig : micro_vs_long ] .
we presented a novel density functional , based on the weighted density paradigm , able to describe classical inhomogeneous fluids in a large portion of their phase diagram , critical point included .
this is the first attempt to describe the effects of correlations induced by attractive interactions in confined fluids .
the theory is based on the description of the uniform system provided by the hierarchical reference theory , one of the few liquid state approaches able to cope with long range density fluctuations .
this technique , applied to the evaluation of the effective interaction between two hard walls in a fluid , allowed for the investigation of the crossover between a depletion - like mechanism at high temperatures and the critical casimir effect emerging near the critical point .
our method does not rely on a long wavelength approximation and provides a complete picture of the solvent mediated force for any wall separation , displaying the presence of important non universal contributions in the effective interaction at short distances , even in the critical region .
we believe that this dft will be useful in investigating other correlated systems , when density fluctuations are expected to play an important role .
we showed that , at large separations , the solvent mediated force per unit surface between the walls decays exponentially on the scale of the correlation length in the whole portion of the phase diagram to the left of the line .
such a behavior can not be considered as a signature of the onset of critical casimir effect : only the product between the _ amplitude _ of the long range exponential tail and the cube of the correlation length is a genuine universal quantity .
our microscopic approach allows for the determination of the universal quantities characterizing the critical casimir effect , namely the scaling function @xmath236 , both along the critical isochore ( @xmath188 ) and in the regime .
strong corrections to scaling have been observed in the hcy fluid we investigated : the universal features appear to emerge only in a narrow neighborhood of the critical point , at least in the model we examined .
it would be useful to compare this prediction with numerical simulations for the hcy fluid model as well as with theoretical investigations of other systems , like the ising model , where the correction to scaling may be weaker .
these studies will hopefully clarify the origin of the discrepancy between the hrt estimate of the universal amplitude ratio @xmath228 and the commonly accepted value .
the approach presented in this work allows for further improvements . in our density functional ,
the ideal gas term and the hartree contribution to the internal energy have been treated exactly , while the remaining term has been approximated by use of a weighted density trick .
the next step will be to treat the hard sphere term by the fundamental measure theory , known to be very accurate in dealing with excluded volume effects , limiting the weighted density contribution only for the residual correlation term .
this adjustment is expected to increase the accuracy of the theory at high density , without however modifying the description of the universal properties of the critical casimir effect .    in this first application we just considered a planar geometry ,
whose implications for the phenomenon of colloidal aggregation depend upon further assumptions , namely the derjaguin approximation , which however turns out to be rather inaccurate when the two external bodies are not very close .
a natural further step will be to perform the functional minimization in cylindrical geometry , appropriate for dealing with two spherical particles thereby avoiding any mapping .
finally , a very interesting application of this formalism will be the investigation of temperatures , where wetting phenomena are expected close to the first order transition boundary .
the hierarchical reference theory of fluids provides a consistent description of the full transition line and then it appears to be the natural starting point for the development of a microscopic theory , beyond mean field , for the study of phase coexistence near a wall .