this work is part of a series of papers devoted to the study of gas adsorption and capillary condensation in silica aerogels@xcite .
our goal is to provide a theoretical interpretation of @xmath0he adsorption experiments that have been performed over the last 15 years in several low - temperature physics groups@xcite .
more generally , we wish to propose a theoretical framework for describing the _ nonequilibrium _ ( hysteretic ) behavior of fluids in disordered porous solids@xcite .
it is worth noting that the experimental studies of helium in aerogel were originally aimed at understanding the effect of quenched random impurities on _ equilibrium _ critical phenomena and phase separation . in particular ,
the study of the liquid - vapor transition of @xmath0he followed from the suggestion that the critical behavior of mixtures or pure fluids confined in porous media or in contact with a gel network could be interpreted in terms of the random field ising model ( rfim)@xcite . in this perspective , the use of aerogel was motivated by its very tenuous structure ( which minimizes the effects of confinement ) and by the fact that the porosity ( and thus the amount of disorder ) can be varied in a controlled way .
however , in spite of the early claim that mesoscopic equilibrium phase separation of @xmath0he could be observed in a @xmath1 porosity aerogel@xcite , leading to an estimation of the order parameter critical exponent @xmath2 , it now seems more and more questionable that true equilibrium behavior is observed in experiments .
in particular , recent measurements of sorption isotherms using a high - precision capacitive technique@xcite show that hysteresis between filling and emptying persists to temperatures very close to the critical temperature of pure helium ( @xmath3 , and gradual adsorption is always observed above @xmath4 , instead of the sharp vertical step expected for two - phase coexistence .
similar observations have been made by another group@xcite , using a low - frequency mechanical pendulum technique and , more recently , coupling pressure and optical measurements .
these results are reminiscent of the behavior of fluids confined in dense porous glasses which is usually interpreted in terms of capillary condensation .
although this moves away the exciting perspective of studying the putative equilibrium vapor - liquid critical point of @xmath0he inside aerogel , the remarkable changes in the hysteretic behavior of the adsorbed fluid with temperature and porosity remain to be explained . in particular
, the hysteresis loops in very light aerogels have a rectangular ( instead of triangular ) shape at very low temperature ( for instance , at @xmath5 in @xmath6 porosity aerogel@xcite ) with a very steep adsorption branch that suggests the presence of a genuine , but _ nonequilibrium _ , first - order phase transition .
this is in agreement with our theoretical description@xcite which predicts the change from a continuous to a discontinuous adsorption isotherm as one decreases the temperature ( at constant porosity ) or increases the porosity ( at constant temperature ) , the jump in the fluid density corresponding to a macroscopic avalanche in the system .
( note that we only focus here on adsorption .
the behavior on desorption is also interesting but is related to other physical processes as discussed in refs.@xcite . )
this scenario , which resembles the one predicted for the field - driven rfim at zero temperature@xcite and observed in some disordered magnetic systems@xcite , suggests the existence of a line of disorder - induced critical points in the temperature - porosity diagram with an associated scaling behavior . in the present work
we further elaborate on this issue by presenting new results concerning the avalanche properties as a function of temperature and porosity and by studying the scaling behavior of the adsorption isotherms in the vicinity of the nonequilibrium critical points .
we are also interested in characterizing the avalanches geometrically and studying their possible relation with the aerogel structure .
the paper is arranged as follows . in section 2
we review the model and the theory . in section 3
we present the results for the statistical properties of avalanches . in section 4
we discuss the scaling behavior of the adsorption isotherms and present the phase diagram of the model in the temperature - porosity plane .
we conclude , in section 5 , with a discussion of the relevance of our results to experiments .
since both the model and the theory were described in detail in ealier work ( see in particular ref.@xcite ) , we will restrict the presentation to some key aspects that are relevant to our present goal .
the gel - fluid system is modelled by a _ coarse - grained _ lattice - gas where each cell of the size of a silica particle ( typically a few nanometers@xcite ) is occupied either by the solid or the fluid . the distribution of silica particles on the lattice is obtained by a diffusion - limited cluster - cluster aggregation ( dlca ) algorithm@xcite , so as to model the intricate microstructure of actual base - catalized silica aerogels . on the other hand
, fluid particles can equilibrate with a reservoir that fixes their chemical potential @xmath7 and the temperature @xmath8 .
multiple occupancy of a site is forbidden and only nearest - neighbors ( nn ) attractive interactions are considered . for reasons that are explained in ref.@xcite ,
we use a body - centered cubic ( bcc ) lattice and periodic boundary conditions ( @xmath9 where @xmath10 is the linear size ) .    for a given configuration of the solid matrix ,
fluid configurations are obtained by minimizing a grand - potential functional in the local mean field approximation@xcite , @xmath11&=k_bt \sum_{i } [ \rho_i\ln \rho_i + ( \eta_i-\rho_i)\ln ( \eta_i-\rho_i ) ]   -w_{ff }   \sum_{<ij > } \rho_i\rho_j    -w_{sf }   \sum_{<ij > } [ \rho_i(1-\eta_j ) + \rho_j(1-\eta_i)]-\mu \sum_{i } \rho_i\end{aligned}\ ] ] where @xmath12 ( @xmath13 ) is the thermally averaged fluid density at site @xmath14 and @xmath15 is the quenched occupancy variable for the gel particles ( @xmath16 if the site is occupied by the gel and @xmath17 otherwise ) . @xmath18 and @xmath19 denote the fluid - fluid and solid - fluid interactions , respectively , and the double summations run over all distinct pairs of nn sites .
the gel porosity is given by @xmath20 and the ratio @xmath21 controls the wettability of the solid surface .
@xmath22 $ ] thus defines a grand - potential landscape whose extrema are obtained by solving the corresponding euler - lagrange equations @xmath23 . at fixed @xmath8 and @xmath24 , and for a given realization of the solid ,
this yields a set of @xmath25 coupled nonlinear equations    @xmath26 ) } } \ , \end{aligned}\ ] ]    where @xmath27 and the sum runs over the nearest neighbors of site @xmath14 .
a crucial feature of these equations is that they may have a large number of solutions at low temperature ( and in a certain range of @xmath24 ) as a consequence of the randomness of the solid matrix@xcite .
we believe that this feature is the key for understanding the physics of adsorption in disordered porous solids . especially important
are the local _ minima _ of @xmath22 $ ] which are obtained by solving the set of eqs .
( 2 ) iteratively@xcite and which correspond to metastable states ( to be distinguished from the global minimum of @xmath22 $ ] that defines the true equilibrium state ) .
note that despite its mean - field character ( by which thermal fluctuations are only crudely described ) , the present approach fully accounts for the disorder - induced fluctuations .
in contrast with the standard global mean - field approximation , our treatment still allows to describe strongly inhomogenous configurations and to predict non - classical critical behavior .
computing the metastable states is however not sufficient for describing the nonequilibrium behavior of the system , and one also needs to specify a dynamics .
our theoretical description@xcite is based on the use of a non - ergodic , deterministic , _ zero - temperature_-like dynamics which neglects all thermally activated processes that would allow the system to cross the ( grand - potential ) barriers : a change in the fluid configurations is thus only due to a change in the landscape that comes itself from a change in the temperature or from a variation of the external field ( here the chemical potential @xmath24 ) . as @xmath24 varies ( for instance , when the pressure in the reservoir is slowly increased at constant @xmath8 ) , the system either follows the minimum in which it was trapped as this minimum deforms gradually , or it falls instantaneoulsy into another minimum when the former reaches its stability limit .
this latter move is a discontinuous and irreversible process , an _ avalanche _ , which is at the origin of the history - dependent behavior of the system .
the avalanche corresponds to some _ collective _ condensation event inside the gel network which manifests itself by a jump in the adsorption isotherm ( whereas the gradual deformation of a minimum only corresponds to a slight swelling of the liquid domains ) .
the size of the avalanche is thus given by the size of the jump , i.e. the discontinuous change in the average fluid density @xmath28 . as was shown qualitatively in ref.@xcite
, these condensation events occur at different length scales inside the gel , depending on chemical potential , temperature , and porosity .
in particular , when @xmath29 is large enough , for instance @xmath30 , there exists a critical point separating a high - temperature regime where all condensation events are of microscopic or mesoscopic size ( and the adsorption isotherm is smooth in the thermodynamic limit ) , and a low - temperature regime where a macroscopic fraction of the gas condenses , which yields a finite jump in the isotherm .
it is one goal of the present work to put these observations on a more quantitative basis .
we shall discuss below in more detail the relevance of this theoretical description to the actual experimental situation .
let us just recall that this picture is expected to be valid when there are widely separated time scales in the problem .
specifically , the observation time must be much larger than the local equilibration time ( i.e. the time needed to relax from a marginally stable state to a nearby metastable state ) but much smaller than the time associated to thermal activation ( i.e. the time that it takes to escape from a local minimum via a nearby saddle point)@xcite .
this requires that energy barriers are much larger than temperature , which is indeed the case in many experimental situations dealing with random systems .
it may also be noted that our definition of avalanches generalizes the one used at @xmath31 ( see e.g. ref.@xcite ) .
we indeed introduce an explicit temperature - dependence in the problem by allowing the ( mean - field ) free - energy landscape itself to change with @xmath8 ( for a similar approach in the context of random magnets , see e.g. ref .
as shown in our previous works , this dependence explains quite convincingly the variations observed in the hysteretic behavior of gases adsorbed in disordered solids at low enough temperature .
as discussed in the preceding papers of this series@xcite , the number , size and shape of the condensation events along an adsorption isotherm depend on porosity and temperature .
they also depend on the value of the chemical potential at which they occur .
the random character of these events ( each finite - size sample has a different isotherm ) requires a statistical study of their properties , an issue that has been investigated in great detail in the context of the field - driven rfim@xcite .
the first problem is to identify each avalanche properly . at @xmath31 , when the microscopic variables can only take discrete values ( such as the ising spins in the rfim ) , this is an easy task because the energy landscape only changes through the disappearance of minima , which induces the avalanches .
one can then vary the external field until a spin becomes unstable ( this spin is unique when the random field distribution is continuous ) and then count the number of spins that flip during the avalanche to determine its size . at finite temperature , on the other hand , one has to deal with continuous variables ( the local fluid densities @xmath12 in the present case ) , and since @xmath32 , the increment in the chemical potential along the adsorption isotherm , is always finite , one has to distinguish between the gradual deformation of a minimum and the jump associated to the passage to another local minimum , both leading to a change in the total density . moreover , if @xmath32 is not small enough , several avalanches may overlap , which biases the statistics of their number or size distribution .
although , in principle , one could precisely locate the beginning of each avalanche by following the evolution of the eigenvalues of the hessian matrix @xmath33 with @xmath24 or @xmath8 ( an avalanche starts when a local spinodal corresponding to the vanishing of a certain eigenvalue is reached@xcite ) , this is a daunting numerical task for a large system ( with @xmath34 ) and it is more advisable to use another procedure .
let us first note that the determination of the location and size of the avalanches depends on the precision in the computation of the metastable states . since these latter are obtained by solving the set of eqs . ( 2 ) iteratively , it must be checked that the results are not significantly modified when the convergence criterium of the iterative procedure is changed .
all the results presented here were computed with a precision @xmath35 ( i.e. , the calculation was stopped at iteration @xmath36 when @xmath37 ) .
it was found that decreasing @xmath38 to @xmath39 did not change the results by more than @xmath40 ( but in turn the computational time was increased by a considerable factor ) .
porosity aerogel at @xmath41 ( @xmath42 is the average fluid density inside the aerogel and @xmath24 the gas chemical potential ) .
calculations were done in a sample of size @xmath43 using the dichotomy procedure explained in the text.,width=302 ]    consider now the variations @xmath44 in the local fluid densities when increasing the chemical potential by a very small ( but finite ) increment @xmath32 ( for instance @xmath45 , with @xmath18 taken as the energy unit ) .
we show in fig .
1 a small portion of a typical low - temperature isotherm computed in a @xmath46 porosity sample of size @xmath43 at @xmath41 ( @xmath47 ) .
the random character of the condensation events results in the fine staircase structure of the curve , with the vertical steps corresponding to the avalanches separated by nearly horizontal portions corresponding to the smooth deformation of the grand - potential local minima ( note the small scale on the vertical axis ) .
the slope of these near - plateaus depends on temperature and its order of magnitude can be estimated from the isothermal compressibility of the pure lattice gas .
it turns out that _ all _ the @xmath12 s vary by a very small amount in these flat portions of the isotherms ( less than @xmath48 ) whereas the vertical steps are due to a significant variation of _ some _ of the @xmath12 s ( a low temperature , this typically corresponds to a gazeous region becoming liquid ) .
we can use this feature to identify the avalanches unambiguously . specifically , we consider that the change @xmath49 in the isotherm is only due to a smooth deformation of a grand - potential minimum when @xmath50 is smaller that some threshold , which we take equal to @xmath51 .
otherwise there is an avalanche , and we associate to it all the sites for which @xmath52 .
the size of the avalanche is then @xmath53 , where the sum is restricted to these sites@xcite ( by taking @xmath54 one would neglect the small increase of @xmath42 which is due to the  smooth " variation of the @xmath12 s on the other sites of the lattice ) .
this criterium thus implies that avalanches for which @xmath55 can not be detected . by choosing a smaller threshold , e.g. @xmath56 , one may find new avalanches , but their size @xmath57 is very small ( always less than @xmath58 ) . as a consequence ,
avalanches with @xmath59 will not be considered in the following . in any experiment
there is also a threshold below which signals can not be detected .    since one must collect the results of many samples to get a good statistics of the avalanches properties , it is impossible in practice to study a whole isotherm with an increment @xmath32 as small as @xmath60 .
it is much more efficient to proceed by dichotomy , comparing the configurations obtained for two values of @xmath24 , say @xmath61 and @xmath62 ( initially , @xmath61 is a large negative value and @xmath63 , the value at saturation ) . if @xmath64 between the two configurations , one then compares the configurations for @xmath61 and @xmath65 on the one hand , and @xmath65 and @xmath62 on the other hand .
one then looks for an avalanche in the two intervals , and the same procedure is repeated until no avalanche is found between two configurations , or the difference in chemical potential becomes smaller than @xmath60 ( this is the method that was used for obtaining the isotherm shown in fig .
1 , which explains that the points are irregularly spaced ) .    finally , to characterize geometrically an avalanche , which is a random object , one may define an effective  radius of gyration " @xmath66 by attributing to each site @xmath14 a  mass " that corresponds to the change in the local density , @xmath67 where the sum is restricted to the sites contributing to the avalanche ( i.e. with @xmath68 ) and @xmath69 is the  center of mass " of the avalanche@xcite .
let us first study how the number of avalanches with a given size varies with temperature and porosity . to illustrate the general trend , we consider isotherms obtained in @xmath46 and @xmath1 porosity aerogels , as shown in figs .
2(a ) and 3(a ) .
the sequence of avalanches was recorded in @xmath70 samples of linear size @xmath43 and @xmath71 samples of linear size @xmath72 , respectively ( this approximately corresponds to the same ratio @xmath73 , where @xmath74 is the correlation length of the aerogel@xcite ) . figs .
2(b ) and 3(b ) show @xmath75 , the ( unnormalized ) avalanche size distribution ( @xmath76 , the total number of avalanches of size larger than @xmath58 ) . as @xmath8 increases
, @xmath77 decreases from @xmath78 to @xmath79 and from @xmath80 to @xmath81 , respectively .
the size of the largest possible avalanche corresponding to a simultaneous change of all the @xmath12 s from @xmath82 to @xmath17 is @xmath83 .
one thus has @xmath84 and @xmath85 in the @xmath46 and @xmath1 aerogels , respectively .
it is readily seen that there exist much larger avalanches in the lighter aerogel .
this is a direct consequence of the larger void space available , as discussed in our previous papers@xcite . in particular , at @xmath86 , some avalanches span the whole system , signaling that a macroscopic avalanche possibly occurs in the thermodynamic limit .
note however that the statistics for very large avalanches is not good , which explains the erratic behavior seen in fig .
3(b ) for large @xmath57 .
moreover , finite - size effects play an important role in the large - porosity / low - temperature regime and they should be taken into account to properly describe the evolution towards the thermodynamic limit ( see ref.@xcite for a careful study in the case of the rfim ) . unfortunately , this would require a considerable amount of numerical work , and we shall therefore only concentrate our discussion on the non - spanning avalanches .
porosity aerogel at @xmath87 , and @xmath88 ( from top to bottom ) ; data correspond to an average over @xmath70 gel samples of size @xmath43 .
( b)corresponding ( unnormalized ) avalanche size distributions @xmath75 .
the solid lines represent the fit @xmath89 ( color on line).,width=377 ]     porosity aerogel at @xmath90 , and @xmath91 ; data correspond to an average over @xmath71 gel samples of size @xmath72 ( color on line).,width=377 ]    figs .
2(b ) and 3(b ) show that the evolution of @xmath75 with temperature is quite similar in the two aerogels : as @xmath8 increases , both large _ and _ small avalanches are suppressed , and , as a consequence , the total number @xmath77 decreases .
only the number of avalanches of intermediate size stays roughly constant with @xmath8 .
if one refers to what happens in the field - driven rfim at @xmath31 , this result may appear surprising .
indeed , the isotherms get smoother as one increases the temperature , like the magnetization curves @xmath92 in the rfim as one increases the width of the random field distribution@xcite .
this seems to imply that the temperature plays the role of an effective disorder in the system .
therefore , one would expect that large avalanches become rare as @xmath8 increases ( which is indeed what happens ) , but also that small avalanches proliferate , so that @xmath77 increases on average@xcite .
it turns out , however , that temperature also plays another role in our system , a role which is well documented in the literature devoted to capillary condensation@xcite : adsorption in a small cavity becomes more gradual as @xmath8 increases .
specifically , mean - field calculations for a single pore of simple geometry ( for instance a slit of width @xmath93 ) predict that adsorption is continuous and reversible above some critical temperature @xmath94 which decreases as @xmath93 decreases .
since small avalanches are associated with gas condensation in small cavities or crevices defined by the neighboring gel strands , this phenomenon explains the decrease of the number of small avalanches ( when @xmath95 , the temperature at which the hysteresis disappears , there are no avalanches at all and the isotherms become reversible@xcite ) .
it is clear that the random - field picture can only be applied at a scale much larger than the correlation length of the solid : schematically , capillary condensation prevails for @xmath96 whereas the disordered character of the aerogel structure plays a major role for @xmath97 .
the two mechanisms compete at the scale @xmath98 which explains that the number of avalanches of intermediate size does not vary significantly with @xmath8 .
if we now concentrate on the disorder - dominated regime , the most remarkable feature , quite visible in the @xmath99 gel at @xmath86 , is that @xmath75 follows an approximate power - law behavior over several decades .
this is the defining signature of criticality , an issue that has been at the focus of intense research activity in recent years , for instance to explain the scaling properties of the barkausen noise in magnetic materials@xcite . in the rfim picture ,
the power - law behavior is related to the existence of a disorder - induced critical point@xcite . in this context
, the avalanche size distribution should be analyzed using the form @xmath100 where @xmath101 is a critical exponent@xcite , @xmath102 is a scaling function , and @xmath103 a cut - off that depends on the distance to the critical point ( @xmath104 diverges as the critical disorder is approached ) .
our results are however too limited to obtain sensible scaling collapses ( although a crude attempt shows that the scaling function is non - monotonic , in agreement with the rfim behavior@xcite ) .
we simply replace the function @xmath102 by an exponential form only describing the cut - off of the power - law at the scale @xmath104 . figs .
2(b ) and 3(b ) show the fits of @xmath75 in the two aerogels according to this approximate treatment .
the power - law regime is very limited in the @xmath46 aerogel even at the lowest temperature , indicating that the system is clearly outside the critical region . on the other hand , in the @xmath1 gel ,
the curves show a linear portion on a log - log scale that increases as @xmath8 decreases : the cutoff @xmath104 is pushed to larger values of @xmath57 ( from @xmath105 to @xmath106 as @xmath107 varies from @xmath91 to @xmath88 ) and reaches the box size for @xmath86 ( our set of data is however too limited to estimate the critical exponent associated to the divergence of @xmath104 ) .
the exponent @xmath101 increases and becomes approximately equal to @xmath108 in the fully developed power - law regime at @xmath86 .
taking into account the fact that our systems are small and that @xmath101 in a finite system depends on disorder@xcite , this value does not appear inconsistent with that predicted for the rfim , @xmath109@xcite . in any case , the gel - fluid system is clearly close to criticality at this temperature .     recorded in the @xmath46 porosity samples as a function of chemical potential @xmath24 at @xmath86 ( a ) and @xmath110 ( b),title="fig:",height=264 ]   recorded in the @xmath46 porosity samples as a function of chemical potential @xmath24 at @xmath86 ( a ) and @xmath110 ( b),title="fig:",height=264 ]     porosity samples at @xmath86 ( a ) and @xmath111 ( b).,title="fig:",height=264 ]   porosity samples at @xmath86 ( a ) and @xmath111 ( b).,title="fig:",height=264 ]    the number and the size of the avalanches along an adsorption isotherm strongly depend on the value of the chemical potential .
we thus also recorded the value of @xmath24 at which each avalanche of size @xmath57 occured . figs .
4 and 5 show the corresponding point clouds in the @xmath46 and @xmath1 gels at @xmath86 and @xmath88 , and @xmath86 and @xmath91 , respectively . as a general rule ,
the size of the largest avalanches increases with @xmath24 , especially in the last stage of the adsorption ( note the logarithmic scale on the vertical axis ) .
but the most striking feature is that avalanches of very different sizes occur in a small interval of @xmath24 .
for instance , in the @xmath46 gel at @xmath86 and @xmath112 , @xmath57 varies from @xmath58 to more than @xmath113 . in the @xmath1 gel at the same temperature and @xmath114 ,
the avalanches span almost all possible sizes between @xmath58 and @xmath115 .
these observations are in complete contradiction with the picture based on the independent pore model which is often used to interpret capillary condensation phenomena@xcite . in this model , only pores of a certain characterisic size @xmath93 fill at a given value of the chemical potential , and the adsorbed fluid density is directly related to the number of pores of this size that can be found in the material .
( one then try to extract from the isotherm the pore size distribution . )
this picture is clearly wrong in the present case because of the very open and intricate character of the gel structure that allows condensation events of very different sizes .
the geometrical characterization of the avalanches is rather difficult . using @xmath116 as
an estimate of the extension of an avalanche is misleading because events with very different shapes are observed ( see for instance figs .
6 and 7 in ref.@xcite ) .
this variety in the geometry is reflected in the values of the radius of gyration @xmath66 , as defined by eq .
the value of @xmath66 is indeed proportional to the maximal extension of the avalanche , but the prefactor depends on its shape . figs .
6 and 7 show how @xmath66 varies with @xmath57 in the @xmath46 and @xmath1 aerogels .
one can see that avalanches with a similar size may have quite different radii of gyration .
for instance , in the @xmath46 aerogel at @xmath86 , avalanches of size @xmath117 have a radius that varies between @xmath118 and @xmath58 ( note that for a compact object , a factor of @xmath119 in @xmath66 would correspond to a factor @xmath120 in @xmath57 ) .
similarly , in the @xmath1 aerogel at the same temperature , @xmath66 varies from @xmath58 to @xmath121 for avalanches of size @xmath122 .
it is therefore impossible to associate a unique characteristic length to an avalanche of a given size .
moreover , by examining the avalanches that occur in a small interval of @xmath24 , we find that they also have very different radii of gyration .    as the temperature increases , the values of @xmath66 for a given size @xmath57 are less scattered .
indeed , avalanches become more compact as they correspond to condensation events that occur in a single cavity of the aerogel .
this evolution is very clear in the @xmath46 aerogel : at @xmath110 , the radius of gyration of most avalanches is between @xmath119 and @xmath118 and there are almost no avalanches with @xmath66 larger than 6 . on the other hand ,
there is still a large domain of variation of @xmath66 in the @xmath1 aerogel at @xmath111 , and one should go to even higher temperatures to observe a significant reduction .
there is in fact a remarkable resemblance between figs .
6(a ) and 7(b ) which suggests that the avalanche properties are similar apart from a scale factor of about @xmath119 .
the adsorption isotherms have also similar shapes , as can be seen in figs .
2(a ) and 3(a ) .
this confirms the observation made in ref.@xcite that the adsorption processes in a low porosity gel at low temperature and in high porosity gel at a higher temperature have similar properties .     and size @xmath57 of the avalanches in the @xmath46 porosity samples at @xmath86 ( a ) and @xmath110 ( b ) .
the inset in ( a ) is a magnification of the point cloud of the spherical avalanches ( see text ) in log - log scale ; the solid line represents the equation @xmath123 ( color on line).,title="fig:",height=264 ]   and size @xmath57 of the avalanches in the @xmath46 porosity samples at @xmath86 ( a ) and @xmath110 ( b ) .
the inset in ( a ) is a magnification of the point cloud of the spherical avalanches ( see text ) in log - log scale ; the solid line represents the equation @xmath123 ( color on line).,title="fig:",height=264 ]     porosity samples at @xmath86 and @xmath91 . in the inset of ( b ) the solid
line represents the equation @xmath124 ( color on line).,title="fig:",height=264 ]   porosity samples at @xmath86 and @xmath91 . in the inset of ( b )
the solid line represents the equation @xmath124 ( color on line).,title="fig:",height=264 ]      in both fig .
6(a ) and fig .
7(b ) , one may notice that a small set of avalanches stands out . in the @xmath46 gel , at @xmath86 ,
these avalanches have a size @xmath57 in the range @xmath125 and a smaller radius of gyration than the other ones .
they are therefore more compact , and a direct observation of some of them shows that their shape is almost perfectly spherical .
this is confirmed by a systematic study of the asphericity parameter which is always smaller than @xmath126@xcite . as shown in the inset of fig .
6(a ) on a log - log scale , the radius of gyration of these avalanches varies approximately like @xmath127 which is indeed close to the expected value for spherical compact objects ( @xmath128^{1/3}\approx 0.48 s^{1/3}$ ] for a sphere of radius @xmath129 ) .
the radius @xmath129 varies between @xmath130 and @xmath118 , this latter value corresponding approximately to @xmath74 which is also more or less the size of the largest voids in the aerogel @xcite .
moreover , as shown in fig .
8 , most of these avalanches occur at the end of the adsorption process : they thus correspond to the condensation of the last gas bubbles that remain in the large cavities .     of spherical avalanches along the adsorption isotherm in the @xmath46 porosity aerogels at @xmath86 .
the average isotherm is shown as a continuous line.,width=226 ]    a similar behavior is observed in the @xmath1 aerogel at @xmath111 ( fig .
the avalanches are again approximately spherical and occur in the last stage of the adsorption process .
however , since the temperature is higher , the boundaries of the gas bubbles are less well defined , and the radius of gyration of varies approximately like @xmath131 .
the largest observed value of @xmath132 is again comparable to the size of the largest cavities in the aerogel ( @xmath133 ) .
although there are also some other remaining bubbles of more complicated shape ( see figs . 7 and 9 in ref.@xcite ) , these spherical avalanches actually represent the only case where one can relate directly the characteristic size of the condensation events to the structure of the aerogel .
we stress again that , in general , the avalanches have a complicated shape so that one can not describe the condensation of the adsorbed gas in terms of  constant curvature bubbles"@xcite .     as a function of the inverse radius of the spherical avalanches in the @xmath46 and @xmath1 aerogels at @xmath86 and @xmath111 , respectively .
( b ) is a magnification of the small point cloud in ( a ) near the origin corresponding to the avalanches in the @xmath1 aerogel .
the solid line with a slope @xmath17 corresponds to the kelvin equation ( see text ) ( color on line).,width=359 ]    the fact that spherical avalanches are clearly identified gives us the opportunity for testing a macroscopic description of the gas condensation in terms of a competition between volume and surface contributions .
such a description underlies the classical interpretation of capillary condensation based on the kelvin equation@xcite that is commonly used to extract from the isotherms information on the characteristic pore size in the solid .
we thus write the change in the grand potential associated to the condensation of a spherical gas bubble of radius @xmath129 at the chemical potential @xmath24 as@xcite @xmath134 \frac{4\pi}{3}r^3-\gamma 4\pi r^2\ ] ] where @xmath135 is the corresponding bulk pressure , @xmath136 is the pressure of the metastable liquid at the same chemical potential ( both pressures are computed from the equation of state of the bulk lattice gas ) , and @xmath137 is the ( planar ) liquid - gas surface tension .
all these terms except @xmath129 are temperature dependent . assuming that condensation occurs at equilibrium , the kelvin equation ( in its laplace form ) corresponds to setting @xmath138 in the equation . expanding as usual @xmath139 and @xmath140 about @xmath141 and to first - order in @xmath142 ( which is valid for small undersaturations ) , this yields @xmath143 $ ] , where @xmath144 and @xmath145 are the densities of the coexisting liquid and gas at this temperature .
since the surface tension on a lattice depends on the orientation of the surface with respect to the lattice vectors , we estimate @xmath137 by taking the average of the three values calculated in the principle lattice planes @xmath146 , @xmath147 and @xmath148@xcite . the corresponding  kelvin plot " is shown in fig .
we see that the  macroscopic " kelvin equation provides the correct trend for the variation of @xmath149 as a function of @xmath150 , in spite of the fact that actual condensation does not occur at equilibrium ( @xmath151@xcite ) and that the adsorbent potential and the curvature corrections are neglected in eq .
( 4 ) ( note that @xmath149 is ten times smaller in the lighter gel , which is due to the effect of temperature , correctly described by the term @xmath152 , and to the change in the average radius @xmath129 by a factor of about @xmath119 ) .
the scatter in the data , however , shows that the kelvin equation is only an approximation .
the existence of a critical point in our system , at which an infinite avalanche occurs for the first time in the thermodynamic limit and the avalanche size distribution has a power - law behavior , suggests the possibility of performing also a scaling collapse of the adsorption isotherms . in the case of the nonequilibrium rfim at @xmath31 where the disorder is controlled by the width @xmath153 of the distribution of the random fields ,
standard renormalization group arguments suggest that the magnetization curve @xmath92 for different values of @xmath153 may be collapsed using the scaling form@xcite    @xmath154    where @xmath155 is the critical magnetization ( i.e. the magnetization at the critical field @xmath156 for the critical disorder @xmath157 ) , @xmath158 and @xmath159 are the reduced disorder and reduced field respectively and @xmath160 is a universal scaling function ( @xmath161 refers to the sign of @xmath162 ) .
the estimated values of the critical exponents in three dimensions are @xmath163 and @xmath164@xcite .
is rescaled by the density at @xmath165 ) . from left to right : @xmath166 at @xmath167 .
the curves result from an average over @xmath70 samples of size @xmath72 ( color on line).,width=340 ]    the situation in our system is more complicated than in the rfim because the hamiltonian ( or , more correctly , the grand - potential functional @xmath168 ) contains three parameters that can be tuned independently ( the porosity @xmath29 , the interaction ratio @xmath169 , and the chemical potential @xmath24 ) , and , moreover , the temperature controls the ruggedness of the free - energy landscape . although a rigorous proof is lacking , it seems however reasonable to assume that the critical behavior is governed by a unique fixed point , which may or may not be the same as for the @xmath31 nonequilibrium rfim . in other words
, there are only two relevant variables from a renormalization - group standpoint : the chemical potential ( that plays the role of the magnetic field @xmath170 in the rfim ) and an effective disorder @xmath153 which is some smooth funtion of @xmath169 , @xmath29 and @xmath8 .
( in addition , one could also have to take into account the microstructure of the aerogel , on top of the dependence on the porosity : in particular , one may suspect that the fractal character of very light aerogels , e.g. @xmath171 , modifies the values of the critical exponents or , at least , induces a different scaling behavior in some intermediate regime ; this possibility is not taken into account in the present study which does not consider aerogels of porosity larger then @xmath1 . ) for fixed @xmath169 and @xmath8 , there is a critical value @xmath172 of the porosity , and a critical value @xmath173 of the chemical potential .
alternatively , for fixed @xmath169 and @xmath29 , there is a temperature @xmath174 and a chemical potential @xmath175 at which the system becomes critical .
if this picture is correct , one expects to observe the same critical behavior whether the critical line @xmath176 ( or @xmath177 ) is approached along a constant-@xmath178 or a constant-@xmath8 trajectory in the temperature - porosity diagram ( for simplicity , we forget the dependence on @xmath169 which is kept constant and equal to @xmath119 in all our computations@xcite ) .
consider for instance the collection of isotherms displayed in fig .
10 which covers a significant range of temperatures and porosities .
this would correspond to @xmath179 for the real @xmath0he - aerogel system , but this must be considered only as a crude estimation : we indeed recall that our treatment is not aimed at a quantitative description of actual systems because of the coarse - grained character of the model and the use of a ( local ) mean - field approximation .
note that the density of the adsorbed fluid is rescaled by its value at @xmath165 so that all the curves go to @xmath17 as @xmath180 ( one has @xmath181 where @xmath182 is the density of the bulk liquid at coexistence ) .
isotherms for @xmath183 and @xmath184 according to eq .
8 where @xmath185/\mid \phi- \phi_c(t)\mid^{\beta \delta}$ ] and @xmath186 with @xmath187 , @xmath188 , and @xmath189 ( color on line).,width=302 ]    let us first study the scaling behavior when varying the porosity at fixed temperature .
by analogy with the rfim , our basic scaling hypothesis is the existence , in the limit @xmath190 , of only two scaling variables @xmath191 and @xmath192 that measure the distance to the critical point @xmath193 . accordingly
, the adsorption isotherms in the vicinity of the critical point @xmath194 should scale as @xmath195 where @xmath196 are the corresponding scaling functions and @xmath197 and @xmath198 are two nonuniversal scale factors that are temperature - dependent ( they actually depend on all the specific features of the model@xcite : @xmath169 , the lattice structure , the range of interaction , etc ... ) .
note that the scaling variables @xmath199 and @xmath200 need not be exactly @xmath201 and @xmath202 , but their dependence on @xmath29 and @xmath24 should be analytic@xcite . unless one uses the proper variables , the curves do not collapse . in the case of the rfim@xcite , it was shown that the proper scaling must be done with replacing @xmath159 in eq .
( 6 ) by the  rotated " variable @xmath203 , where @xmath204 is some non - universal constant ( on the other hand , there is no need to rotate @xmath162 because the corresponding correction is irrelevant when @xmath205 ) . in practice
, the rotation of the x - axis can be avoided by simply lining up the peaks in the derivative of the curves , with no correction terms , as is often done in analyzing experimental data ( by considering the slope @xmath206 instead of @xmath207 we also get rid of the dependence on @xmath208 ) . what to do on the y - axis is a matter of choice and , for simplicity , we just take @xmath209 . (
note that in ref.@xcite the curves @xmath210 in the rfim are plotted using @xmath211 instead of the more standard choice @xmath212 ; alternatively , it has been proposed to expand @xmath162 to second - order in @xmath213@xcite ; these different choices are equivalent asymptotically if no rotation is required , but they may change significantly the amplitude of the critical regions . ) therefore , to summarize , we look for a scaling of the form @xmath214 where @xmath209 , @xmath215 is the value of the chemical potential at the maximum of @xmath206 , @xmath216 is the derivative of @xmath196 , and @xmath217 and @xmath218 account for the dependence on @xmath197 and @xmath198 . since it is important to have smooth isotherms so that one can take the derivative numerically , the average must be performed over a significant number of samples : the curves shown in fig .
10 result from an average over @xmath70 gel realizations of size @xmath72 ( however , there are still some annoying irregularities in the average slopes , especially around @xmath219 ) .
calculations are of course done in finite systems so that one should also take into account finite - size effects@xcite .
these are especially important in the low - disorder regime ( @xmath220 in which there is a discontinuity in the isotherms in the thermodynamic limit .
this introduces an additional difficulty in the analysis and it is a very demanding numerical task to perform a systematic finite - size scaling study of all the curves ( see ref.@xcite for such a study of the @xmath221 isotherm in the @xmath1 porosity aerogel ) . as a consequence , we shall only consider here the isotherms that are clearly in the strong disorder regime , which leads to reasonably size - independent results : this eliminates in particular the @xmath222 and @xmath223 isotherms in the @xmath1 porosity aerogel .    fig .
11 shows the scaling plot of the slope of the @xmath224 isotherms for @xmath225 and @xmath184 obtained with the choice @xmath226 , @xmath187 , and @xmath189 ( the scale factors @xmath227 and @xmath228 are arbitrarily fixed equal to @xmath17 ) .
it is actually very difficult to adjust independently the three unknown parameters @xmath229 , @xmath2 and @xmath176 ; the values of the exponents used here are the ones that allow for a collapse of a whole set of isotherms , as discussed below .
it can be seen that the @xmath230 and @xmath231 isotherms collapse onto a single curve , which is not the case of the @xmath232 isotherm , especially on the left - hand side of the peak .
this part of the isotherm is indeed dominated by the presence of a liquid layer that coats the aerogel strands@xcite , and in the low - porosity gels the remaining available void space is too small for allowing large avalanches to spread out .
this contribution to the adsorption has nothing to do with a critical phenomenon and it is reasonable to eliminate from the scaling plot all the isotherms obtained in the @xmath46 porosity aerogel . more generally , at each temperature , we choose to discard the isotherms that do not belong to the scaling region because the porosity is obviously too far from @xmath176 ( i.e. , the distance from the effective critical disorder is too large ) .
this finally leads to the global scaling plot shown in fig .
12 , obtained by adjusting also the nonuniversal scale factors @xmath217 and @xmath218 ( which , however , always remain of order @xmath17 ) .
the overall quality of the collapse is not perfect , but similar to that found in the rfim@xcite ( it is actually better because we only keep two isotherms at each temperature@xcite ) .
in fact , as discussed in ref.@xcite , the magnetization curves ( or , here , the isotherms ) are quantities that do not collapse well , which makes it difficult to obtain from the scaling plots a precise determination of the critical exponents or of the location of the critical point .
although our estimated values of the exponents @xmath229 and @xmath2 are compatible with those of the rfim , the uncertainty is much too large to conclude that our model belongs to the same universality class . on the other hand ,
it is noteworthy that the value @xmath233 for the exponent @xmath229 is consistent with the value @xmath234 for the finite - size scaling exponent @xmath235 obtained in ref.@xcite , if one uses for @xmath236 the rfim value , @xmath237@xcite .     ) . @xmath238 and @xmath169 are the corresponding scaling variables and scaling functions ( color on line).,width=302 ]     ( or @xmath176 ) below which there is a first - order nonequilibrium transition in the adsorption isotherm .
the squares indicate the isotherms that are included in the scaling plot of fig
. 12 ( color on line).,width=340 ]     and @xmath239 isotherms in the @xmath1 aerogel according to eq .
9 where @xmath185/\mid t / t_c(\phi)-1\mid^{\beta \delta}$ ] and @xmath240 with @xmath187 , @xmath188 , and @xmath241 ( color on line).,width=302 ]    from the values of @xmath176 extracted from the scaling plots at each temperature , we can map out an approximate phase diagram of the model in the temperature - porosity plane .
this is shown in fig .
13 , where the critical line has been forced to extrapolate towards the zero - disorder limit @xmath242 .
( let us recall however that , in a real system , there can not be a clear - cut transition close to @xmath243 because thermal fluctuations become larger and can allow the system to cross the free - energy barriers on the experimental time - scale ; moreover , as noted before , we expect the fractal character of the microstructure in very light aerogels to come into play . )
the present calculations indicate that the critical temperature drops very quickly with the porosity , and the extrapolation to @xmath31 predicts that there is no transition at finite temperature in gels of porosity smaller than @xmath244 .
this is in agreement with finite - size scaling calculations performed in the @xmath46 porosity aerogel that show that the isotherms are continuous at all temperatures in the thermodynamic limit . on the other hand , in the @xmath1 porosity aerogel
, the study of ref.@xcite predicts a transition at @xmath221 , which is indeed below the critical line of fig . 13 .
morever , as shown in the preceding section , the avalanche size - distribution has a power - law behavior at @xmath86 , which is the critical temperature extracted from the scaling plot ( @xmath245 ) .
it turns out that the same critical temperature is also obtained by studying the scaling behavior when varying the temperature at fixed porosity .
the scaling is now performed according to the expression    @xmath246    where @xmath247 and @xmath248 and @xmath249 are porosity - dependent scale factors .
it is clear from fig .
13 that this procedure can only be used with the @xmath1 porosity aerogel .
the scaling plot of the @xmath110 and @xmath239 isotherms is shown in fig .
14 ( again , the @xmath250 isotherm is discarded because of the presence of important finite - size effects ) .
it can be seen that a good scaling is obtained with the same values of the critical exponents as used before . remarkably , the resulting scaling function is also the same and it can be included in fig .
12 by adjusting the scale factors @xmath248 and @xmath249 .
this is a clear indication that the critical behavior of the system is controlled by the same fixed point for all temperatures and porosities .
we now discuss the comparison of all these results with experiments .
in this work and our earlier work@xcite , we have presented a detailed numerical study of the adsorption of @xmath0he in silica aerogel that explains the modification in the shape of the isotherms as the temperature or the porosity is changed .
we have found that the adsorption occurs via a sequence of irreversible condensation events ( avalanches ) of varying size and shape .
the traditional picture of capillary condensation in mesoporous solids in which pores of increasing size fill up successively as the chemical potential or the gas pressure increases does not hold in aerogel , except , approximately , for the special category of spherical avalanches that occur at the end of the adsorption process .
our calculations predict that there exists a line of critical points in the temperature - porosity diagram separating two regimes of continuous and discontinuous adsorption and characterized by a power - law behavior of the avalanche size distribution .
there is a rather large scaling region where adsorption isotherms at different temperatures and different porosities collapse onto a single universal curve .
although the system sizes considered in this work are too small to reach a definite conclusion , the estimated values of the critical exponents seem to be consistent with those of the three - dimensional field - driven rfim at @xmath31@xcite .    the qualitative agreement between the simulated and experimental loops@xcite is a strong argument in favor of the validity of the above scenario .
however , a complete validation still awaits the experimental description of the scaling behavior of the isotherms and the observation of a power - law behavior in the avalanche size distribution .
this latter issue requires some additional comments .
it may be indeed distressing that avalanches have not been observed so far with @xmath0he in aerogel .
in fact , the only exemple of avalanche behavior in adsorption / desorption experiments concerns the draining of superfluid @xmath0he from nuclepore@xcite , a porous material that contains a random spatial distribution of interconnected cylindrical pores . in this case
, superfluidity appears to be a critical factor for the existence and propagation of avalanches via the presence of a superfluid film on the nuclepore substrate that provides a coupling mechanism even between distant pores .
it is likely however that such a mechanism is not necessary for avalanches to occur in aerogels because of very open character of the microstructure .
there are actually several requirements for observing avalanches in sorption experiments , which are all satisfied in the case of the experiments with nuclepore@xcite , but not in the case of the existing aerogel studies@xcite .
let us first recall that , by definition , avalanches correspond to some collective reorganization of the system _ in response _ to a variation of the external driving force . if one is changing the density rather than the chemical potential or the pressure ( by introducing a controlled amount of gas into the sample and then measuring the pressure change ) , the system may follow a different path among the metastable states , and the vertical jumps in the isotherms that are the signature of the avalanches may not be observed@xcite ( this could explain that there are many data points along the adsorption isotherm at @xmath5 in @xmath6 aerogel displayed in ref.@xcite ) .
secondly , as already noted , the observation time must be larger than the time the system takes to respond to the change in the external control variable .
this can be done by continuously changing the chemical potential at a very slow rate , and measuring the amount of gas in the solid regularly ( every @xmath238 seconds ) . in practice , this procedure is only possible if the ( local ) equilibration time is short ( for instance , it is smaller than one second with superfluid helium in nuclepore ) .
unfortunately , the relaxation times in aerogel are much longer ( of the order of two hours with @xmath0he at low temperature in the experiments of ref.@xcite ) , and one has to change @xmath24 or @xmath251 by a finite step and keep it fixed until a new local equilibrium ( i.e. , a new metastable state ) is reached . with this procedure , the issue of having a system with a long equilibration time boils down to stability of the experimental conditions and of course patience ; in principle , one could then observe avalanches . yet ,
if @xmath149 or @xmath252 is not small enough , the avalanches may overlap and the resulting isotherm may appear continuous .
this is probably the situation with the experiments in ref.@xcite where @xmath253 ( to be compared with @xmath254 in the nuclepore experiments@xcite ) .
finally , a last requirement is that the sensitivity of the measurements is large enough to discriminate the jumps in the adsorbed fluid density .
this condition is more easily satisfied when these jumps are large enough , which implies that one must be close to the critical line displayed in fig .
although one must be very cautious in comparing our temperature scale with the true experimental one , it seems that this last condition is also not satisfied in the existing experiments .
f. detcheverry , e. kierlik , m. l. rosinberg , and g. tarjus , phys .
e * 68 * , 61504 ( 2003 ) .
f. detcheverry , e. kierlik , m. l. rosinberg , and g. tarjus , langmuir * 20 * , 8006 ( 2004 ) .
a. p. y. wong and m. h. w. chan , phys .
65 * , 2567 ( 1990 ) .
d. j. tulimieri , j. yoon , and m. h. w. chan , phys .
lett . * 82 * , 121 ( 1999 ) .
j. beamish and t. herman , physica b * 340 * , 329 ( 2003 ) ; j. beamish and t. herman , j. low temp . phys . * 134 * , 339 ( 2004 ) ; t. herman , ph.d .
thesis , university of alberta ( 2005 ) ; t. herman , j. day , and j. beamish , preprint cond - mat/0505430 . c. gabay , p.e .
wolf , and l. puech , physica b * 284 * , 97 ( 2000 ) ; c. gabay , f. despetis , p.e . wolf , and l. puech , j. low temp .
121 * , 585 ( 2000 ) ; c. gabay , ph.d .
thesis , universit joseph fourier , grenoble i ( 2001 ) ; t. lambert , c. gabay , l. puech , and p.e .
wolf , j. low temp .
phys . * 134 * , 293 ( 2004 ) ; t. lambert , ph.d .
thesis , universit joseph fourier , grenoble i ( 2004 ) .
e. kierlik , p. a. monson , m. l. rosinberg , l. sarkisov , and g. tarjus , phys .
lett . * 87 * , 055701 ( 2001 ) ; e. kierlik , p. a. monson , m. l. rosinberg , and g. tarjus , j. phys . : condens . matter .
* 14 * , 9295 ( 2002 ) ; m.l .
rosinberg , e. kierlik , and g. tarjus , europhys . lett . * 62 * , 377 ( 2003 ) ; f. detcheverry , e. kierlik , m. l. rosinberg , and g. tarjus , physica b * 343 * , 303 ( 2004 ) ; f. detcheverry , e. kierlik , m. l. rosinberg , and g. tarjus , adsorption * 11 * , 115 ( 2005 ) .
f. brochard and p.g .
de gennes , j. phys .
lett . * 44 * , l785 ( 1983 ) ; p. g. de gennes , j. phys . chem . * 88 * , 6469 ( 1984 ) .
j. p. sethna , k. dahmen , s. kartha , j. a. krumhansl , b. w. roberts , and j. d. shore , phys .
* 70 * , 3347 ( 1993 ) ; o. perkovic , k. dahnen , and j. p. sethna , phys .
75 * , 4528 ( 1995 ) .
o. perkovic , k. dahnen , and j. p. sethna , cond - mat/9609072 ( unpublished ) ; k. dahmen and j. p. sethna , phys .
b * 53 * , 14872 ( 1996 ) ; o. perkovic , k. dahnen , and j. p. sethna , phys . rev .
b * 59 * , 6106 ( 1999 ) .
a. berger , a. inomata , j. s. jiang , j. e. pearson , and s. d. bader , phys .
rev . lett . * 85 * , 4176 ( 2000 ) .
j. marcos , e. vives , l. maosa , m. acet , e. duman , m. morin , v. novak , and a. planes , phys .
b * 67 * , 224406 ( 2003 ) .
p. meakin , phys .
lett . * 51 * , 1119 ( 1983 ) ; m. kolb , r. botet , and r. jullien , phys .
* 51 * , 1123 ( 1983 ) .
more precisely , we expect the number of solutions to scale exponentially with the system size .
the fact that only minima are obtained by solving iteratively nonlinear equations such as eqs .
( 2 ) is shown by d. ling , d. r. bowman , and k. levin , phys .
b * 28 * , 262 ( 1983 ) . in practice ,
as explained in ref.@xcite , we use a more efficient algorithm than simple iteration .
however , we checked numerically that it yields the same solutions as the simple iteration procedure . in actual sorption experiments , there are also kinetic effects associated to mass and/or heat transfer which may significantly slow down the relaxation to a new metastable state .
in particular , silica aerogels are exceptional thermal insulators and convection in the adsorbed fluid is probably suppressed by the aerogel strands@xcite .
these effects are neglected in the present treatment in which the fluid responds instantaneously to a change in chemical potential or temperature .
mass transfert can be incorporated in an effective grand - canonical treatment by rescaling the relaxation time so as to take into account the constraint of mass conservation ( woo and monson , phys .
e * 67 * , 041207 ( 2003 ) ) . in any case , at low enough temperature , the ( local ) relaxation time remains much smaller than the time associated to thermally activated processes between metastable states , which justifies the use of a zero - temperature dynamics .
j. p. sethna , k. dahmen , and c. r. myers , nature * 410 * , 242 ( 2001 ) . c. m. soukoulis , k. levin , and g. s. grest , phys .
b * 28 * , 1495 ( 1983 ) . f. j. perez - reche and e. vives , phys . rev .
b * 67 * , 134421 ( 2003 ) ; f. j. perez - reche and e. vives , phys . rev .
b * 70 * , 214422 ( 2004 ) .
it is easy to see from eq .
( 1 ) that the elements of the hessian matrix @xmath255 are given by @xmath256 -w_{ff}n_{ij}$ ] where @xmath257 is the kronecker symbol and @xmath258 when @xmath14 and @xmath259 are nearest neighbors and zero otherwise .
@xmath255 can thus be diagonalized by a unitary transformation .
the threshold @xmath260 ensures that the avalanche size is computed with enough precision .
the calculation of @xmath66 is restricted to avalanches that do not span the whole lattice .
r. evans , j. phys .
: condensed matter * 2 * , 8989 ( 1990 ) .
for a recent review on the barkhausen effect , see g. durin and s. zapperi in _ the science of hysteresis _ , edited by g. bertotti and i. mayergoyz , elsevier ( 2004 ) .
since we consider the avalanche size distribution integrated over a whole isotherm , our exponent @xmath101 identifies with the exponent @xmath261 in the notation of ref.@xcite . s. j. gregg and k. s. w. sing , _ adsorption , surface area , and porosity _ , academic press , new york ( 1982 ) .
we measure the asphericity via the parameter @xmath262/2r_g^2 $ ] .
the normalisation is such that @xmath263 for spherically symmetric objects and @xmath17 for  cigar " shaped objects .
note that @xmath264 is not an invariant quantity with this definition since the components of the radius of gyration tensor are not taken along the principal axes ( see j. rudnick and g. gaspari , j. phys .
a. * 19 * , l191 ( 1986 ) ) .
it is likely that , on average , we underestimate the actual deviation from sphericity .
r. salazar and l. d. gelb , mol .
phys . * 102 * , 1015 ( 2004 ) .
note however that @xmath265 in eq .
5 is a small quantity compared to the bulk and surface contributions .
v. privman and m. e. fisher , phys .
b * 30 * , 322 ( 1984 ) . see e.g. j. cardy , _ scaling and renormalization in statistical physics _ , cambridge university press , cambridge ( 1996 ) .
except at @xmath110 where we only take into account the isotherm in the @xmath30 porosity aerogel ..
r. evans , u. marini bettolo marconi , and p. tarazona , j. chem .
, farady trans . * 2 * , 82 ( 1986 ) .
lilly , a. h. wooters , and r. b. hallock , phys .
lett . * 77 * , 4222 ( 1996 ) ; m. p.
lilly and r. b. hallock , phys .
b * 64 * , 024516 ( 2001 ) ; m. p. lilly , a. h. wooters , and r. b. hallock , phys . rev .
b * 65 * , 104503 ( 2002 ) . moreover , as emphasized in ref.@xcite , adding volumetric bursts to the system induces an unknown degree of local heating .