Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the cepheid period - luminosity ( pl ) relation is a classical tool widely used to estimate the distance to local group galaxies and to external galaxies with hubble space telescope ( hst ) observations , as well as , through the calibration of secondary distance indicators , to even more distant stellar systems . this explains the dominant role of these variables on the route to solve cosmological problems such as the evaluation of the hubble constant @xmath8 and , in turn , of the age of the universe . the basic physics underlying the cepheid variability suggests that the pulsation period @xmath11 depends on the mass , luminosity and effective temperature of the pulsator . from stellar evolution theory one expects a close relation between mass and luminosity , and the natural result of these theoretical prescriptions is a period - luminosity - color ( plc ) relation , where the pulsator absolute magnitude @xmath12 in each given photometric passband @xmath13 is a linear function of the period and color index [ @xmath14 , as given by    @xmath15\eqno(1)\ ] ]    however , since the pulsation occurs within a finite zone of the hr diagram , the color term in the plc relation is often neglected and cepheid distances are usually estimated from the most favourite pl relation    @xmath16    where @xmath17 is now the average of the cepheid magnitudes for each given period . it should be noticed that whereas the plc relation holds for any individual pulsator , pl is a `` statistical '' solution which depends both on the pulsation boundaries and on the distribution of pulsators within the instability strip . this explains why distance determinations based on pl relations require statistically significant numbers of cepheids in order to reduce the effects of deviations from the ridge line . the pl relation in bolometric magnitude is traditionally assumed to be metal - insensitive ( see , e.g. , iben & renzini 1984 , freedman & madore 1990 ) and cepheid distances are generally derived by adopting _ universal _ pl relations at different wavelengths , having the slope provided by the cepheids in the lmc and the zero - point referenced to the lmc distance , as obtained with independent methods ( see freedman 1988 ; kennicutt et al . 1998 ; walker 1999 , and reference therein ) , or to calibrating galactic cepheids ( see feast & catchpole 1997 ; lanoix et al . 1999 ) .    in the last years , we have deeply investigated the cepheid pulsational behavior through the computations of nonlinear , nonlocal and time - dependent convective pulsating models which take into account the coupling between pulsation and convection . with respect to linear - nonadiabatic models computed by different authors ( e.g. , chiosi et al . 1993 ; saio & gautschy 1998 ; alibert et al . 1999 ) , our theoretical approach allows reliable predictions on the temperature of both the blue and red edges of the instability strip , as well as about the amplitude and morphology of the light - curve ( see bono et al . 1999a [ paper i ] , 2000a [ paper iii ] , 2000b [ paper vi ] ) . concerning the assumptions on the input physics , computing procedures and the adopted mass - luminosity relation , the reader is referred to paper i and paper iii , which contain also the detailed comparison of our models with the linear results in the literature . here we only remark that from pulsating models computed with three different chemical abundances ( @xmath1=0.004 , 0.008 and 0.02 ) we derived that the _ bolometric _ magnitude of metal - rich variables is , on average , fainter than that of metal - poor stars with the same period ( bono et al . 1999b [ paper ii ] ) . moreover , we found that both the slope and zero - point of synthetic pl relations at different wavelengths depend on the pulsator metallicity , with the amplitude of the metallicity effect decreasing from visual to near - infrared magnitudes ( see paper ii ; caputo et al . 2000 [ paper v ] ) . also the predicted plc relations at the different wavelengths turned out to be , in various degrees , metallicity dependent : as an example , for a given period and @xmath18 color , metal - rich pulsators are brighter than metal - poor ones , whereas they are fainter if the @xmath19 color is adopted ( see paper v ) . on these grounds , it has been shown that cepheid observations in three filters ( @xmath20 or @xmath21 ) allow to simultaneously constrain distance , reddening and metallicity of variables in the magellanic clouds ( caputo et al . 1999 [ paper iv ] ) and in the milky way ( caputo et al . , in preparation ) . in this paper we will apply our theoretical scenario to the cepheids observed within two hst surveys : the `` extragalactic distance scale key project '' ( hereafter kp , see freedman et al . 1994a ) and the `` type ia supernova calibration '' ( hereafter snp , see saha et al . the adopted procedure is briefly presented in sect . 2 and the predicted reddenings and distances are given in sect . the metallicity effects on the cepheid distance scale and on the value of the hubble constant @xmath8 are discussed in sect . the concluding remarks close the paper . the results presented in paper v predict that the intrinsic dispersion of the synthetic pl relations with fixed metallicity decreases with increasing the filter wavelength , in close agreement with the observed trend ( see madore & freedman 1991 ; tanvir 1999 ) . moreover , we found that the predicted plc relations are almost independent of the pulsator distribution within the instability strip and reproduce the tight correlation between pulsation period and photometric parameters of individual pulsators , as earlier suggested by sandage ( 1958 ) , sandage & gratton ( 1963 ) and sandage & tammann ( 1968 ) . since hst observations are in the two passbands @xmath22 and @xmath23 , we adopt the predicted pl(@xmath23 ) and plc(@xmath0 ) relations given in paper v and reported here in table 1 and table 2 , respectively , to get the average reddening and true distance modulus of the cepheid samples from the observed intensity - weigthed mean magnitudes @xmath24 and @xmath25 , for each assumption about the cepheid metallicity . with respect to the classical method based on pl(@xmath22 ) and pl(@xmath23 ) , our procedure aims at reducing the spurious effects caused by both the large dispersion of the observed pl(@xmath22 ) relation and the uncertainty on the predicted pl(@xmath22 ) relation as due to different pulsator distributions ( see paper ii and paper v ) . the offset of each cepheid from the predicted pl(@xmath23 ) and plc(@xmath0 ) relations gives the apparent distance moduli @xmath26=@xmath3+@xmath27 and @xmath28=@xmath3+@xmath29 , respectively . adopting @xmath30=3.1 and @xmath31=1.88 from cardelli et al . ( 1989 ) , the reddening @xmath2 and true distance modulus @xmath3 of each cepheid are derived . eventually , by averaging over the sample of variables , we provide the predicted mean values , for the three adopted metallicities . cccc z & @xmath32 & @xmath33 & @xmath34 + & & & +   + & & & + 0.004&-2.00@xmath350.01&-3.00@xmath350.01&0.14 + 0.008&-2.03@xmath350.01&-2.88@xmath350.01&0.14 + 0.02&-2.18@xmath350.01&-2.53@xmath350.01&0.12 + & & & +    ccccc z & @xmath36 & @xmath37 & @xmath38 & @xmath34 + & & & & +   + & & & & + 0.004 & -3.55 @xmath350.03 & -3.58 @xmath350.03 & 3.75 @xmath350.07 & 0.03 + 0.008 & -3.54 @xmath350.03 & -3.59 @xmath350.02 & 3.74 @xmath350.06 & 0.03 + 0.02 & -3.61 @xmath350.03 & -3.59 @xmath350.04 & 3.85 @xmath350.09 & 0.03 + & & & & + table 3 lists the kp galaxies together with the color excess @xmath2 and true distance modulus @xmath3 , as given by the authors ( see references in the last column ) , and the oxygen - to - hydrogen abundance ratios ( o / h ) of hii regions , as listed in ferrarese et al . for some of these galaxies the cepheid data have been re - analyzed by the kp team and the new reddening and distance values ( ferrarese et al . 1999 ) are given in italics in the same table . the galaxies marked with an asterisk , originally studied by the snp group , were extracted from the hst archive and re - processed by the kp team ( gibson et al . 2000 ) . for them we list in table 3 the kp values , while the snp original distances are reported in table 4 . galaxy & @xmath2 & @xmath39 & 12+log(o / h ) & reference +   + ngc 925 & 0.13 @xmath35 0.02 & 29.84 @xmath350.16 & 8.55@xmath350.15 & silbermann et al . 1996 +   + ngc 1326a & 0.05 @xmath35 0.09 & 31.27 @xmath350.21 & 8.50@xmath350.15 & prosser et al . 1999 + & _ 0.00 @xmath35 0.01 _ & _ 31.43 @xmath35 0.17 _ & & ferrarese et al . 1999 + ngc 1365 & 0.12 @xmath350.06 & 31.31 @xmath350.20 & 8.96@xmath350.20 & silbermann et al . 1999 + & _ 0.15 @xmath35 0.03 _ & _ 31.39 @xmath35 0.19 _ & & ferrarese et al . 1999 + ngc 1425 & 0.07 @xmath350.03 & 31.73 @xmath350.23 & 9.00@xmath350.15 & mould et al . 2000 + & _ 0.07 @xmath35 0.03 _ & _ 31.81 @xmath35 0.17 _ & & ferrarese et al . 1999 +   + ngc 3031 & 0.03 @xmath350.05 & 27.80 @xmath350.20 & 8.75@xmath350.15 & freedman et al . 1994a +   + ngc 3198 & 0.05 @xmath350.03 & 30.80 @xmath350.23 & 8.60@xmath350.15 & kelson et al . 1999 + ngc 3319 & 0.06 @xmath350.08 & 30.78 @xmath350.17 & 8.38@xmath350.15 & sakai et al . 1999 +   + ngc 3351 & 0.12 @xmath35 & 30.01 @xmath350.19 & 9.24@xmath350.20 & graham et al . 1997 + ngc 3368 & 0.11 @xmath35 & 30.20 @xmath350.19 & 9.20@xmath350.20 & gibson et al . 2000 + ngc 3627 * & 0.11 @xmath350.04 & 30.06 @xmath350.23 & 9.25@xmath350.15 & gibson et al . 2000 +   + ngc 4414 & 0.01 @xmath350.05 & 31.41 @xmath350.19 & 9.20@xmath350.15 & turner et al . 1998 +   + ngc 4725 & 0.19 @xmath350.03 & 30.50 @xmath350.16 & 8.92@xmath350.15 & gibson et al . 1999 + & _ 0.19 @xmath35 0.03 _ & _ 30.57 @xmath35 0.18 _ & & ferrarese et al . 1999 +   + ngc 4321 & 0.10 @xmath350.06 & 31.04 @xmath350.17 & 9.13@xmath350.20 & freedman et al . 1994b + ngc 4548 & 0.09 @xmath350.03 & 31.01 @xmath350.28 & 9.34@xmath350.15 & graham et al . 1999 +   + ngc 4496a * & 0.03 @xmath350.05 & 31.02 @xmath350.17 & 8.77@xmath350.15 & gibson et al . 2000 + ngc 4535 & 0.08 @xmath350.02 & 31.02 @xmath350.26 & 9.20@xmath350.15 & macri et al . 1999 + & _ 0.10 @xmath35 0.03 _ & _ 31.10 @xmath35 0.17 _ & & ferrarese et al . 1999 + ngc 4536 * & 0.08 @xmath350.03 & 30.95 @xmath350.17 & 8.85@xmath350.15 & gibson et al . 2000 +   + ngc 4639 * & 0.04 @xmath350.03 & 31.80 @xmath350.18 & 9.00@xmath350.15 & gibson et al . 2000 +   + ngc 5253 * & 0.10 @xmath35 & 27.61 @xmath350.19 & 8.15@xmath350.15 & gibson et al . 2000 +   + ngc 5457 & 0.00 @xmath350.04 & 29.34 @xmath350.17 & 8.50@xmath350.15 & kelson et al . 1996 +   + ngc 7331 & 0.15 @xmath350.05 & 30.89 @xmath350.14 & 8.67@xmath350.15 & hughes et al . 1998 +   + ic 4182 * & -0.02 @xmath35 0.03 & 28.36 @xmath350.18 & 8.40@xmath350.20 & gibson et al . 2000 + ngc 2090 & 0.07 @xmath350.02 & 30.45 @xmath350.16 & 8.80@xmath350.15 & phelps et al . 1998 + [ -4pt ] ngc 2541 & 0.08 @xmath350.05 & 30.47 @xmath350.11 & 8.50@xmath350.15 & ferrarese et al . 1998 + ngc 3621 & 0.23 @xmath350.03 & 29.13 @xmath350.18 & 8.75@xmath350.15 & rawson et al . 1997 +    .empirical cepheid distance to snp galaxies , as published in the original papers by the snp team . [ cols="^,^,^ " , ]     by applying the theoretical relations in table 1 and table 2 to the full samples of cepheids provided by kp and snp studies ( references in table 3 and table 4 ) , we derive the predicted reddenings and true distance moduli listed in table 5 and table 6 , respectively , for each of the three selected metal abundances . to avoid any bias , the predicted true distance modulus includes the correction for negative reddening , if present . the total errors to the predicted values are due to the observed scatter in the pl and plc fits added in quadrature to the intrinsic dispersion of the theoretical relations .    considering that both kp and snp distances adopt pl(@xmath22 ) and pl(@xmath23 ) relations based on the cepheids in the lmc , and since the average metallicity of lmc cepheids is @xmath400.008 ( luck et al . 1998 ) , let us first compare the predicted reddening and distance with @xmath1=0.008 to the lmc - based empirical values . it seems worth mentioning that we have already found a fair agreement between the empirical pl(@xmath22 ) and pl(@xmath23 ) relations and our predictions with @xmath1=0.008 ( see paper v ) .    as for the galaxies in table 3 , we obtain that empirical and predicted results with @xmath1=0.008 ( table 5 ) are in mutual agreement , with a difference ( empirical minus predicted ) of @xmath410.05 mag and @xmath420.05 mag , in the mean . we show in fig . 1 that the discrepancy is for each galaxy well within a total uncertainty on the distance of the order of 10% ( dashed lines ) . this is a reassuring result as to the reliability of the adopted method of analysis and our pulsating models with @xmath1=0.008 . moreover , one has to consider that both kp and snp studies assume @xmath43=18.50 mag and @xmath44=0.10 mag , and that the above differences of @xmath410.05 mag and @xmath420.05 mag with respect to the predictions with @xmath1=0.008 could be easily removed by slightly increasing ( decreasing ) the adopted lmc distance ( reddening ) . on the other hand , the theoretical results depend on assumptions such as the adopted mass - luminosity ratio , bolometric corrections and temperature - color relations . thus , revisions of input physics and/or atmosphere models might modify the results in table 5 and table 6 , but with marginal effects on the relative distances and reddenings . as a fact , by applying our method to the lmc cepheids ( madore & freedman 1991 ) we get @xmath4518.6 mag and @xmath460.03 mag . accounting for these results , the actual discrepancy between empirical and predicted results is @xmath470.05 mag and @xmath480.02 mag , in the mean .    on the contrary , it seems important to notice that different distributions of cepheids within the instability strip may produce non - systematic effects on both distance and reddening . as a matter of example , let us consider a cepheid with absolute magnitudes @xmath49 and @xmath50 , true distance modulus @xmath51 and reddening @xmath52 , and let us assume that @xmath53=@xmath54 , where @xmath55 is given by the pl(@xmath13 ) relation . from the pl(@xmath22)+pl(@xmath23 ) procedure one derives @xmath56 , while following our pl(@xmath23)+plc(@xmath0 ) method one has @xmath57 . since cepheids close to the blue edge of the instability strip have positive values of @xmath58 , with @xmath59 as a consequence of the fact that the width of the instability strip decreases towards longer wavelengths ( see paper v ) , the result is that the actual reddening of cepheids near the pulsation blue edge is underestimated , whereas the opposite occurs for variables close the red edge . concerning the distance , from the pl(@xmath22)+pl(@xmath23 ) procedure one has @xmath60 , while the pl(@xmath23)+plc(@xmath0 ) method yields @xmath61 . with @xmath62 0.35 mag and @xmath63 0.25 mag for cepheids close to the pulsation blue edge 1.0 - 1.5 . ] , the derived reddening and true distance modulus are @xmath10 0.10 mag lower than the actual values , whatever is the adopted method of analysis . on these grounds , one understands that statistically significant samples of variables are required in order to get accurate distance and reddening determinations . otherwise , selection effects on distant cepheids ( the variables close to the blue edge are brighter than those near the red edge , for fixed period ) lead to underestimate both the distance and reddening . as already mentioned , for some of the kp galaxies in table 3 the cepheid data have been reviewed by the same team , leading to a slight systematic increase on the distance . from the data listed in italics in table 3 one finds that the average upward revision is of the order of @xmath100.09 mag . as discussed in ferrarese et al . ( 1999 ) , this is a consequence of the adopted lower cut - off in period , as the empirical distance modulus increases when the lower period cut - off is moved from 10 to 20 - 25 days . on the theoretical side , we show in table 5 ( see values in italics ) that using the same selection as in ferrarese et al . ( 1999 ) causes a systematic increase of the predicted true distance modulus by only @xmath100.05 mag . this suggests that the cut - off in period has minor influence upon our pl(@xmath23)+plc(@xmath0 ) procedure , supporting the belief that robust distance determinations should adopt only plc relations , for which cepheid observations in at least three filters are needed . concerning the snp galaxies , we give in table 6 the predicted distances as obtained from the full set of cepheid data provided in the original papers , while the values in italics refer to the cepheid samples as selected by the snp team . one notices that the adopted selection criteria yield rather discordant effects on the predicted distance . with @xmath1=0.008 , the true distance modulus of four galaxies increases ( from 0.05 mag for ngc 4496a to 0.44 mag for ngc 3627 ) , whereas it decreases by @xmath640.11 mag for ic 4182 . moreover , limiting the comparison to the snp empirical distances obtained from the selected samples ( table 4 ) , the difference with respect to the predicted distance modulus with @xmath1=0.008 ( in italics in table 7 ) ranges from @xmath640.30 mag ( ic 4182 ) to + 0.22 mag ( ngc 4536 ) , which is somehow larger than that ( from @xmath640.12 mag for ic 4182 to + 0.18 mag for ngc 3627 ) found by adopting kp results for the same galaxies ( see table 3 and table 5 ) . it is known that the disagreement between snp and kp distances to the same galaxy ( up to 0.47 mag for ngc 5253 ) deals with different pipeline processing procedures , selection criteria and data reduction strategies , and the interested reader can find a detailed discussion in gibson et al . ( 2000 ) and tammann et al . ( 2000 ) . here we have to face up to the evidence that by using the cepheid data selected by snp studies the overall discrepancy between empirical and predicted distances with @xmath1=0.008 is somehow larger ( see fig . 2 ) than that found with kp data ( filled circles in fig . 1 ) . for this reason , the rest of our analysis will deal exclusively with the results of kp studies .    throughout the above discussion we used the cepheid predicted distance with @xmath1=0.008 , neglecting the evidence that the oxygen abundance of the galaxies in table 3 varies from [ o / h]@xmath650.50 to @xmath10 + 0.40=log(o / h)-log(o / h)@xmath66 , with log(o / h)@xmath66=-3.10 . ] . since for the lmc [ o / h]@xmath650.40 ( pagel et al . 1978 ) , one has that the hii metallicity of the cepheid fields observed with hst may be up to 0.8 dex higher than the `` reference '' sample of lmc variables . on the other hand , the data in table 5 show that the predicted cepheid reddening and true distance modulus decrease as the adopted metallicity @xmath1 increases . this is a consequence of the fact that our metal - rich pulsating models are , on average , fainter than metal - poor ones . resuming the discussion in sect . 3 , let us assume that for our cepheid holds @xmath670 with respect to the pl relations with @xmath1=0.008 . this yields @xmath68=@xmath69 and @xmath70=@xmath51 . however , the variable is brighter with respect to the pl relations with @xmath1=0.02 , by @xmath710.25 mag and @xmath720.20 at log@xmath731.0 ( see paper v ) . as a consequence , @xmath740.06 mag and @xmath750.11 mag . the decrease of the reddening with increasing the adopted metallicity yields that almost all the galaxies show unpleasant negative values of @xmath2 with @xmath1=0.02 . however , as already mentioned , our theoretical scenario suggests @xmath760.03 mag for the cepheids in the lmc , while both the kp and snp studies assume @xmath44=0.10 mag . thus , were our predictions calibrated on this galaxy , then all the reddening values in table 5 and table 6 should be increased by 0.07 mag .    since our purpose is the evaluation of the metallicity effects on the distance scale , in the mean we derive @xmath770.27@xmath350.04 mag dex@xmath6 , at least within the explored range 0.004@xmath78 0.02 . the sense of this result is that when universal relations are used , the cepheids more metal - rich than the calibrators will appear spuriously _ brighter_. on this basis , the predicted correction ( in magnitude ) to empirical lmc - calibrated distance moduli is given by    @xmath79    where @xmath7 is the difference between the cepheid metallicity and the lmc value of @xmath1=0.008 .    in the following section we investigate into the effects of this metallicity correction to kp distances by assuming that the cepheid metal content scales with the hii region metallicity , i.e. adopting @xmath7=@xmath80[o / h ] , the difference between the galaxy [ o / h ] abundance and that of the lmc ( [ o / h]=@xmath640.40 ) . observational pros and cons for this assumption will be discussed . figure 3 shows the difference between the empirical kp data and the predicted true distance modulus when the oxygen abundance of the cepheid fields is taken into account . as expected , the almost flat distribution in fig.1 now shows a clear correlation to the hii region metallicity , with the discrepancy between lmc - calibrated and theoretical distances increasing when moving from metal - poor to metal - rich galaxies . however , one notices that for very few galaxies the discrepancy overcomes the threshold of 10% ( dashed line ) . of particular interest is the behaviour of the eight galaxies ( the six snp galaxies plus ngc 3368 and ngc 4414 ; filled circles ) which give cepheid - calibrated snia luminosities ( see later ) .    as a first straightforward test to the actual occurrence of a metallicity effect , we consider the kp distance to galaxies members of groups or clusters . the lower panel of fig . 4 shows the residuals @xmath81 and @xmath82[o / h ] of each galaxy from the distance modulus and o / h metallicity as averaged over the galaxies of the same group or cluster . a correlation between the metallicity and distance deviations from the mean values can be detected , best - fitted by the relation ( dotted line )    @xmath83,\eqno(4)\ ] ]    in the sense that galaxies whose o / h metallicity is larger than the average appear to have a _ larger _ distance . we believe that depth - effects within a given group or cluster can not be invoked since there is no reason for which the metal - richest galaxies are also the most distant ones . this unexpected result , which appears surprisingly in agreement with our predicted correction [ eq . ( 3 ) ] , disagrees with earlier observational clues . it is known that studies of different fields in m31 ( freedman & madore 1990 ) and m101 ( kennicutt et al . 1998 ) suggested an opposite metallicity correction on distance ( see also kochanek 1997 ; sasselov et al . 1997 ) . following kennicutt et al . ( 1998 ) , the correction ( in magnitude ) to the true distance modulus is given by    @xmath84\eqno(5)\ ] ]    we show in the upper panel of the same fig . 4 that adopting such an empirical correction would imply an even stronger correlation between distance and metallicity , with a slope of + 0.53 ( dotted line ) hard to accept .    applying the predicted metallicity correction to the absolute distance moduli [ eq . ( 3 ) ] , obviously results in a correction to the @xmath8-values based on lmc - calibrated distances . assuming @xmath85 , one has @xmath86[o / h]@xmath10 + 0.124 dex@xmath6 , i.e. an increase of @xmath10 6% in the lmc - based @xmath8 value of any snia calibrator whose metallicity is 0.5 dex larger than that of the lmc . this is a not dramatic variation , nevertheless it seems interesting to settle whether it works to decrease or increase the present dispersion of @xmath8 values . within this context , it is worth noticing that the well known disagreement between the `` high '' and `` low '' @xmath8 values claimed by kp and snp studies , respectively , are not due entirely to the already mentioned different distances to snia calibrating galaxies . as discussed by gibson et al . ( 2000 ) , the methodology adopted by the snp group with kp distances leads to @xmath87 63 km s@xmath6 mpc@xmath6 dex@xmath6 , which is 9% higher than the snp average value ( @xmath87 58 km s@xmath6 mpc@xmath6 dex@xmath6 ) . moreover , the calibration of snia luminosities presents a variety of different approaches and the same kp distances , coupled with the suntzeff et al . ( 1999 ) procedure , would give @xmath87 67 km s@xmath6 mpc@xmath6 dex@xmath6 . however , the absolute value of the hubble constant is out the purpose of this paper . we aim only at determining the metallicity - correction to lmc - based @xmath8 values , with the hope of providing new elements for reducing the present uncertainty of @xmath1010% to a 5% level . taking at the face value the @xmath88 estimates given by gibson et al . ( 2000 , see their table 6 ) we show in the lower panel of fig . 5 that they agree to each other to within 10% ( dashed lines ) , but with a mild tendency to increase as the oxygen abundance of the host galaxy decreases . the linear regression to the points ( dotted line ) is    @xmath89,\eqno(6)\ ] ]    suggesting a metallicity correction as @xmath86 [ o / h]@xmath10 + 0.069 dex@xmath6 which is roughly half the amount predicted on the basis of eq . ( 3 ) , but it runs towards the same direction . on the contrary , applying the kennicutt et al . ( 1998 ) correction to the true distance moduli [ eq . ( 5 ) ] , results in a correction as @xmath90/@xmath80[o / h]@xmath650.111 dex@xmath6 to the lmc - based @xmath8 values . the upper panel of fig . 5 shows that this would produce a more evident correlation between @xmath8 and [ o / h ] , with two estimates overcoming the discrepancy of 10% from the average value ( dashed lines ) . on the other hand , the lower panel of fig . 6 shows that our predicted metallicity correction @xmath90/@xmath80[o / h]@xmath100.124 dex@xmath6 would yield @xmath8 values which agree to within 10% , but slightly increasing as the [ o / h ] abundance increases . eventually , in order to remove any dependence of @xmath8 on the galaxy metallicity ( see upper panel in fig . 6 ) , we suggest that _ at least _ the empirical evidence in fig . 5a should be taken into consideration , i.e. , @xmath90/@xmath80[o / h]@xmath100.069 dex@xmath6 , leading to a slight upward revision of the kp unweighted mean of @xmath91=67.0@xmath354.3 km s@xmath6 mpc@xmath6 to @xmath91=68.6@xmath353.9 km s@xmath6 mpc@xmath6 .    as said before , the determination of @xmath8 deals with several factors apart from the distance to the snia calibrating galaxies . thus , we are giving the predicted value of the hubble constant , but only the result of the predicted metallicity - correction to the kp @xmath8 estimates , still holding their adopted approach in the treatment of snia data and distance to lmc . hst galaxies with published cepheid distances have been studied to the light of a theoretical pulsational scenario based on nonlinear , nonlocal and time - dependent convective models with three selected metallicities ( @xmath1=0.004 , 0.008 , 0.02 ) . theoretical pl and plc relations in the @xmath0 passbands have been used to predict cepheid distance and reddening for each adopted metal content . the comparison of the predictions with @xmath1=0.008 ( the average metallicity of the lmc ) to the lmc - calibrated empirical results provided by the kp group shows a mutual agreement , better than that found with snp results . the predicted pl and plc relations suggest that the cepheid distance should decrease with increasing the adopted metallicity and that lmc - based distance moduli should be corrected according to @xmath920.27@xmath7 mag dex@xmath6 , where @xmath7 is the difference between the cepheid metallicity and the lmc mean value of @xmath1=0.008 . theoretical studies which adopt a linear pulsational approach , being unable to get firm predictions on the red edge of the instability strip , can not disprove this result .    on the observational front , earlier evidences suggest an opposite metallicity correction as given by @xmath93=+0.24@xmath80[o / h ] mag dex@xmath6 ( kennicutt et al . 1998 ) , where @xmath80[o / h ] is the o / h metallicity difference between the cepheid field and the lmc mean value ( [ o / h]@xmath650.40 ) . however , we show that the kp true distance modulus to galaxies members of clusters or groups seems to be correlated with the galaxy o / h metallicity in good agreement with our predictions . * the comparison between baade wesselink radii for galactic and magellanic clouds cepheids ( laney 1998 , 1999 , 2000 ; feast 1999 ; storm 2000 ) suggests that metal - rich variables are fainter than the metal - poor ones . * the purely `` geometric '' distance to the galaxy ngc 4258 , as obtained from maser sources rotating around the central black hole , is 7.2@xmath350.5 mpc ( herrnstein et al . 1999 ) , whereas the value derived from lmc - calibrated pl relations is 8.1@xmath350.4 mpc ( maoz et al . 1999 ) . since for this galaxy one has [ o / h]@xmath100.12 ( zaritsky et al . 1994 ) , the kennicutt et al . ( 1998 ) correction would increase the distance to @xmath10 8.6 mpc , whereas our predicted correction gives @xmath107.6 mpc . * a carefully selected sample of 236 cepheids from the hipparcos catalogue suggests that the slope of the galactic pl relations could be shallower than that observed for lmc variables ( groenewegen & oudmaijer 2000 ) , thus suggesting some caution in the adoption of universal pl relations .    as a whole , both the predictions and the empirical evidences seem to confirm that the metallicity effect on the cepheid distance scale is not dramatic , with a correction factor of @xmath650.14 mag to the lmc - based distance modulus of any galaxy whose metallicity is 0.5 dex larger than that of the lmc . the consequent predicted metallicity correction to the hubble constant is @xmath90/@xmath80[o / h]@xmath10 + 0.124 dex@xmath6 , which means an increase of @xmath10 6.2% in the lmc - based @xmath8 values of any snia calibrator whose metallicity is 0.5 dex larger than that of the lmc . taking the kp @xmath8 estimates at their face value , we show that the observed correlation with the [ o / h ] metallicity of the cepheid fields suggests @xmath90/@xmath80[o / h]@xmath10 + 0.069 dex@xmath6 , which is more consistent with our prediction rather than the metallicity correction based on the kennicutt et al . ( 1998 ) empirical results ( @xmath90/@xmath80[o / h]@xmath650.111 dex@xmath6 . ) on the basis of this , we eventually suggest that the average value provided by the kp team @xmath967 km s@xmath6 mpc@xmath6 should be increased _ at least _ up to @xmath9 69 km s@xmath6 mpc@xmath6 . the adoption of a metallicity correction to lmc - based empirical estimates of @xmath8 might hopefully reduce the present uncertainty of 10% on the average value to a 5% level . unfortunately , the variety of the current methods of calibrating snia luminosities give very different results for the absolute value of the hubble constant . as an example ( see gibson et al . 2000 ) , the suntzeff et al . ( 1999 ) methodology yields a mean value which is from 3.7% to 8.3% larger than that found with the saha et al . ( 1997 ) procedure , even adopting the same distance to the calibrators . in our belief this is a key problem on the route to settle the absolute value of the hubble constant . we are greatly indebted to our referee for several helpful comments . we thank v. castellani and g. bono for discussion and critical suggestions . financial support for this work was provided by the ministero delluniversit e della ricerca scientifica e tecnologica ( murst ) under the scientific project `` stellar evolution '' ( vittorio castellani , coordinator ) .
Abstract:

Answer: we use theoretical period - luminosity and period - luminosity - color relations in the @xmath0 passbands , as based on nonlinear , nonlocal and time - dependent convective pulsating models , to predict the reddening and true distance modulus of distant cepheids observed with the hubble space telescope . by relying on the pulsating models with metal content @xmath1=0.008 , we find that the theoretical predictions agree to the values obtained by the extragalactic distance scale key project on the basis of empirical period - luminosity relations referenced to lmc variables . in the meantime , from the theoretical relations with @xmath1=0.004 and 0.02 we find that the predicted @xmath2 and @xmath3 decrease as the adopted metal content increases . this suggests a metallicity correction to lmc - based distances as given by @xmath4/@xmath50.27 mag dex@xmath6 , where @xmath7 is the difference between the metallicity of the cepheids whose distance we are estimating and the lmc value @xmath1=0.008 . such a theoretical correction appears supported by an existing , although weak , correlation between the cepheid distance and the [ o / h ] metallicity of galaxies within a given group or cluster , as well as by a similar correlation between the @xmath8 estimate and the [ o / h ] metallicity of the galaxies which calibrate the snia luminosity . on the contrary , the metallicity correction earlier suggested on empirical grounds seem to be excluded . eventually we suggest that the average value @xmath9 67 km s@xmath6 mpc@xmath6 provided by the key project team should increase _ at least _ up to @xmath10 69 km s@xmath6 mpc@xmath6 . further observational evidences in support of the predicted scenario are finally presented .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: systems interacting through long - range forces behave very differently from those in which particles interact through short - range potentials . for systems with short - range forces , for arbitrary initial condition , the final stationary state corresponds to the thermodynamic equilibrium and can be described equivalently by either microcanonical , canonical , or grand - canonical ensembles . on the other hand , for systems with unscreened long - range interactions , equivalence between ensembles breaks down  @xcite . often these systems are characterized by a negative specific heat  @xcite in the microcanonical ensemble and a broken ergodicity  @xcite . in the infinite particle limit , @xmath2 , these systems never reach the thermodynamic equilibrium and become trapped in a stationary out of equilibrium state ( ss )  @xcite . unlike normal thermodynamic equilibrium , the ss does not have maxwell - boltzmann velocity distribution . for finite @xmath0 , relaxation to equilibrium proceeds in two steps . first , the system relaxes to a quasi - stationary state ( qss ) , in which it stays for time @xmath3 , after which it crosses over to the normal thermodynamic equilibrium with the maxwell - boltzmann ( mb ) velocity distribution  @xcite . in the limit @xmath2 , the life time of qss diverges , @xmath4 , and the thermodynamic equilibrium is never reached . unlike the equilibrium state , which only depends on the global invariants such as the total energy and momentum and is independent of the specifics of the initial particle distribution , the ss explicitly depends on the initial condition . this is the case for self - gravitating systems  @xcite , confined one component plasmas  @xcite , geophysical systems  @xcite , vortex dynamics  @xcite , etc  @xcite , for which the ss state often has a peculiar core - halo structure  @xcite . in the thermodynamic limit , none of these systems can be described by the usual equilibrium statistical mechanics , and new methods must be developed .    in this paper we will restrict our attention to self - gravitating systems . unfortunately , it is very hard to study these systems in 3d  @xcite . the reason for this is that the 3d newton potential is not confining . some particles can gain enough energy to completely escape from the gravitational cluster , going all the way to infinity . in the thermodynamic limit , one must then consider three distinct populations : particles which will relax to form the central core , particles which will form the halo , and particles which will completely evaporate . existence of three distinct classes of particles , makes the study of 3d systems particularly difficult . on the other hand , the interaction potential in 2d is logarithmic , so that all the particles remain gravitationally bound . similar to magnetically confined plasmas the stationary state of a 2d gravitational system should , therefore , have a core - halo structure  @xcite . we thus expect that the insights gained from the study of confined plasmas might prove to be useful to understand the 2d gravitational systems . our system consists of @xmath0 particles with the total mass @xmath5 in a two dimensional space . at @xmath6 the particles are distributed over the phase space with the initial distribution @xmath7 , and then allowed to relax . our goal is to calculate the one particle distribution function @xmath8 , once the relaxation process has been completed and the stationary state has been established . for now , we will restrict our attention to the azimuthally symmetric systems .    the mean gravitational potential at @xmath9 at time @xmath10 satisfies the poisson equation @xmath11 where @xmath12 , @xmath13 is the particle number density , and @xmath14 is the gravitational constant . it is convenient to define dimensionless variables by scaling lengths , velocities , potential , and energy to @xmath15 ( arbitrary length scale ) , @xmath16 , @xmath17 and @xmath18 , respectively . in 3d space , our system corresponds to infinitely long parallel rods of line density @xmath19 interacting through a pair potential @xmath20 . if the system has a finite number of particles , after sufficiently large time @xmath3 , it will relax to thermodynamics equilibrium with the mb distribution function , given _ exactly _ by @xmath21 where @xmath22 is a normalization constant , @xmath23 is the lagrange multiplier used to conserve the total energy , and @xmath24 is the potential of mean force  @xcite . for a gravitational system of mass @xmath5 , the correlations between the particles vanish as @xmath0 becomes large , so that @xmath25 . substituting eq . ( [ eqmb ] ) in eq . ( [ poisson0 ] ) , we obtain the classical poisson - boltzmann equation in its adimensional form @xmath26 the solution of this equation is  @xcite @xmath27 for large @xmath9 this potential must grow as @xmath28=0 \,,\ ] ] which requires that @xmath29 and @xmath30 . with these values , the distribution function eq . ( [ eqmb ] ) automatically satisfies the constraint @xmath31 while the value of @xmath32 is obtained from the conservation of energy @xmath33 f({\bf r},{\bf v})={\cal e}_0   \,,\end{aligned}\ ] ] where @xmath34 is the renormalized initial energy , see appendix  [ energy ] . in this paper we will restrict our attention to initial distributions of the water - bag form , @xmath35 where @xmath36 is the heaviside step function and @xmath37 . for simplicity , from now on we will measure all lengths in units of @xmath38 , so that @xmath39 . the renormalized energy , appendix  [ energy ] , then reduces to @xmath40 performing the integral in eq . ( [ conserv2 ] ) with @xmath41 given by eq . ( [ pot pb ] ) , we obtain @xmath42 this provides a complete solution to the equilibrium thermodynamics of 2d self - gravitating system in a large ( but finite ) @xmath0 limit . we next compare the analytical solution presented above with the full @xmath0-body molecular dynamics simulation  @xcite . to do this we calculate the number density of particles between @xmath43 $ ] @xmath44 and the number density of particles with velocity between @xmath45 $ ] , @xmath46    . solid line is the theoretical prediction obtained using the mb distribution function , eq . ( [ eqmb ] ) , and the points are the result of molecular dynamics simulation with @xmath47 particles.,width=340 ]    fig . [ fig-1 ] shows an excellent agreement between the theory and the simulations . it is important to stress , however , that to reach the mb equilibrium distribution , has required a week of cpu time , ( a million dynamical times for @xmath47 particles , see section [ core - halo ] ) . up to the crossover time @xmath3 , the system remained trapped in a quasi - stationary state , with the one particle distribution very different from the equilibrium one . we now turn to the discussion of this non - equilibrium quasi - stationary state . for systems interacting through short - range potentials , the final stationary state corresponds to the thermodynamic equilibrium and is exactly described by the mb distribution . in spite of a popular believe that in the thermodynamic limit for systems with long - range unscreened interactions the mean - field becomes exact , this is not quite true . or rather , this is true mathematically , but is irrelevant for real physical systems , since when @xmath2 it takes an infinite time for such a system to relax to the thermodynamic equilibrium . what is correct , is that in the thermodynamic limit the dynamical evolution of a system with long - range interactions is governed _ exactly _ by the collisionless boltzmann ( vlasov ) equation  @xcite @xmath48 where @xmath49 is the one particle distribution function and @xmath50 is the mean force felt by a particle at position @xmath9 . the mb distribution , together with the poisson equation for the mean - field potential , is a stationary solution of the vlasov equation . thus , if we start with this distribution it is guaranteed to be preserved by the vlasov dynamics . however , unlike for the collisional boltzmann equation , mb distribution is not a global attractor of the vlasov dynamics  an arbitrary ( non - stationary ) initial distribution will not evolve to the mb distribution . thus , the collisionless relaxation described by the vlasov equation is much more complex than the collisional relaxation governed by the boltzmann equation for systems with short - range interactions . the final stationary state of vlasov dynamics will depends explicitly on the initial particle distribution . vlasov equation has an infinite number of conserved quantities , called casimir invariants . any local functional of the distribution function is a casimir invariant of the vlasov dynamics . in particular , if we discretize the initial distribution function into surface levels with values @xmath51 , the hypervolume @xmath52 of each level will be preserved by the vlasov flow . the evolution of the distribution function corresponds to the process of filamentation and proceeds _ ad infinitum _ from large to small length scale . thus on a fine - grain scale , the evolution never stops and the stationary state is never reached . however , in practice , there is always a limit to the maximum resolution , and only a coarse - grained distribution function is available in simulations or in experiments . it is this coarse - grained distribution which appears in practice as the stationary state of a collisionless relaxation dynamics .    numerical solution of the vlasov equation is a very difficult task . since 1960 s there has been a tremendous effort to find an alternative way to predict the final stationary state without having to explicitly solve the vlasov equation  @xcite . one of the first statistical approaches was proposed by lynden - bell and has become know as the _ violent relaxation theory _ . this theory is based on the assumption that there exists an efficient phase space mixing during the dynamical evolution . this assumption is similar to the ergodicity of the boltzmann - gibbs statistical mechanics . for systems with short - range interactions , ergodicity  although very difficult to prove explicitly  is almost always found to be satisfied in practice . this , however , is not the case for the efficient mixing hypothesis for systems with long - range interactions . in fact it was found that for most initial conditions , the phase space mixing is very poor . for magnetically confined plasmas , efficient mixing was found to exist only for very special initial conditions , and in general these systems relax to a stationary state very different from the one predicted by the lynden - bell theory . similarly , for 3d gravitational systems , violent relaxation theory was found to work only if the initial distribution satisfied the , so called , virial condition  @xcite . otherwise strong _ particle - density wave _ interactions broke the ergodicity and resulted in a core - halo phase separation . we first briefly review the violent relaxation theory . the basic assumption of this theory is that during the temporal evolution , the system is able to efficiently explore the whole of phase space . to obtain the stationary ( coarse - grained ) distribution @xmath53 , the initial distribution @xmath54 is discretized into the @xmath55 levels , and the phase space is divided into macrocells of volume @xmath56 , which are in turn subdivided into @xmath57 microcells , each of volume @xmath58 , for a d - dimensional system . since vlasov dynamics is incompressible , @xmath59 , each microcell can contain at most one discretized level @xmath60 . the number density of the level @xmath61 inside a _ macrocell _ at @xmath62  number of microcells occupied by the level @xmath61 divided by @xmath57  will be denoted by @xmath63 . ( 2.5,2.5 ) ( 2.,2.)(0,0)(4,4 ) ( 4.,4.1)(a ) ( 1.,1.)(3.,3 . ) ( -0.1,0)(-0.1,1 ) ( -0.43,0.43)d*r * ( 0,-0.1)(1,-0.1 ) ( 0.43,-0.43)d*v * ( 1.15,2.15)(0.48,2.48)(0.15,1.48)(-1,2.5 )    ( -1.1,2.6)@xmath64    ( 2.5,2.5 ) ( 2.,2.)(0,0)(4,4 ) ( 4.,4.1)(b ) ( 0,0)(0.34,0.34 ) ( 0.34,0.68)(0.68,1 . ) ( 0.68,0.34)(1.,0.68 ) ( 1,0)(1.34,0.34 ) ( 1.34,3.68)(1.68,4 . ) ( 1.0,3.)(1.34,3.34 ) ( 1.68,3.)(2.,3.34 ) ( 1.34,2.68)(1.68,3 . ) ( 3,3)(3.34,3.34 ) ( 2.34,2.34)(2.68,2.68 ) ( 1.68,1.68)(2.,2 . ) ( 1.34,1.34)(1.68,1.68 ) ( 1.68,1.)(2.,1.34 ) ( 2.34,1.68)(2.68,2 . ) ( 3.,1.68)(3.34,2 . ) ( 3.34,1.34)(3.68,1.68 ) ( 1.34,.68)(1.68,1 . ) ( 2.68,2)(3.,2.34 ) ( 2,1.34)(2.34,1.68 ) ( 1,2)(1.34,2.34 ) ( 0.68,1.68)(1.,2 . ) ( 0.34,2.68)(0.68,3 . ) ( 0.34,3.34)(0.68,3.68 ) ( 0,3.68)(0.34,4 . ) ( 2,3.34)(2.34,3.68 ) ( 2.68,3.68)(3.,4 . ) ( 2.34,3)(2.68,3.34 ) ( 3.34,3.68)(3.68,4 . ) ( 3.68,2)(4.,2.34 ) ( 3.34,2.34)(3.68,2.68 ) ( 3.68,2.68)(4.,3 . ) ( 2.34,0.34)(2.68,0.68 ) ( 2,0)(2.34,0.34 ) ( 3.,0.68)(3.34,1 . ) ( 3.34,0.34)(3.68,0.68 ) ( 3.68,0)(4.,0.34 )    ( -0.1,0)(-0.1,1 ) ( -0.43,0.43)d*r * ( 0,-0.1)(1,-0.1 ) ( 0.43,-0.43)d*v *    ( -0.1,3.68)(-0.1,4 ) ( -0.4,3.75)_h _ ( 4,3)a ( 4,4)b ( 5.2,3.4)macro - cell ( 5.2,2.9)@xmath65(*r*,*v * )    ( 4,1.68)c ( 4,2.)d ( 5.2,1.75)micro - cell ( 5.2,1.25)@xmath49(*r*,*v * )    ( 1.15,2.15)(0.48,2.48)(0.15,1.48)(-1,2.5 ) ( 2.45,3.15)(1.48,3.48)(1.15,2.48)(-1,2.6 ) ( -1.1,2.45)@xmath64    note that by construction , the total number density of _ all _ levels in a macrocell is restricted to be @xmath66 see fig . [ fig0 ] . using a standard combinatorial procedure @xcite it is then possible to associate a coarse - grained entropy with the distribution of @xmath67 . the entropy is found to be that of a @xmath55-species lattice gas , @xmath68   + [ 1-\sum_{j=1}^p\rho_j   ( { \bf r},{\bf v } ) ] \ln[1-\sum_{j=1}^p\rho_j ( { \bf r},{\bf v } ) ] \right\ } \;,\end{aligned}\ ] ] with the boltzmann constant set to one . if the initial condition is of the water - bag form , eq . ( [ eq waterbag ] ) ( @xmath69 ) , the maximization procedure is particularly simple , yielding a fermi - dirac distribution , @xmath70}+1}\;,\end{aligned}\ ] ] where @xmath71 is the mean particle energy , and @xmath72 and @xmath73 are the two lagrange multipliers required by the conservation of the total number of particles and the total energy eqs . ( [ eq conserv1 ] ) and ( [ conserv2 ] ) . for a 2d self - gravitating system the virial theorem requires that @xmath74 , in a stationary state appendix  [ virial ] . if the initial distribution does not satisfy this condition , the system will undergo strong oscillations before relaxing into the final stationary state in which the virial theorem is satisfied . for a water - bag initial distribution the virial condition reduces to the requirement that @xmath75 . for future convenience , we will define the virial number for water - bag distributions to be @xmath76 , so that @xmath77 , when the initial distribution satisfies the virial condition . if @xmath78 , the envelope radius , defined as @xmath79 , will vary with time until a stationary state is achieved . note that with the above definition , @xmath80 , as it should be . it is possible to show that the temporal evolution of @xmath81 satisfies @xmath82 where @xmath83 $ ] . the derivation of this equation is given in appendix  [ envelope ] . for an initial water - bag distribution , @xmath84 and @xmath85 so that , if the initial distribution satisfies the virial condition , @xmath86 , and the large envelope oscillations are suppressed . as was already noted for magnetically confined plasmas and 3d self - gravitating systems  @xcite , we expect the violent relaxation theory to work well when the initial distribution satisfies the virial condition and there are no macroscopic envelope oscillations . to check this , we compare the predictions of the theory with the full @xmath0-particle molecular dynamics simulations . at @xmath6 , particles are distributed over the phase space in accordance with the water - bag distribution  ( [ eq waterbag ] ) which satisfies the virial condition , @xmath77 .    ) , while the points are the result of molecular dynamics simulation with @xmath47 particles.,width=340 ]    we then numerically solve the poisson equation ( [ poisson0 ] ) , with the distribution function given by equation ( [ eq lynbell ] ) and compare the results with the molecular dynamics simulations . as can be seen from fig . [ fig1 ] there is a reasonably good agreement between the theory and the simulations .    . solid line is the theoretical prediction obtained using the distribution function of eq . ( [ eq lynbell ] ) , while the points are the result of molecular dynamics simulation with @xmath47 particles.,width=226 ]    however , if the virial condition is not met exactly , one notices a deviation in the tail region of the particle distribution , see inset of fig . [ fig2 ] . for initial distributions with @xmath73 significantly different from @xmath87 , there is a clear qualitative change in the ss distribution function . in this case the original homogeneous cluster , separates into a high density core region surrounded by a diffuse halo , fig . [ fig3 ]  the violent relaxation theory fails completely and a new approach must be developed  @xcite .    . solid line is the prediction of the violent relaxation theory eq . ( [ eq lynbell ] ) , while points are the result of molecular dynamics simulation with @xmath47 particles.,width=340 ] the failure of the violent relaxation theory is a consequence of the inapplicability of the efficient mixing hypothesis to strongly oscillating gravitational systems , figs . [ fig2 ] and  [ fig3 ] . density oscillations excite parametric resonances which favor some particles to gain a lot of energy at the expense of the rest . the resulting particle - wave interactions are a form of non - linear landau damping which allows some particles to escape from the main cluster to form a diffuse halo . the process of evaporation will continue as long as the oscillations of the core persist . oscillations will only stop when the core exhaust all of its free energy , and its effective temperature drops to @xmath88 , @xmath89 in eq . ( [ eq lynbell ] ) . note that because of the incompressibility restriction imposed by the vlasov dynamics eq . ( [ eq vol ] ) , the core can not freeze  collapse to the minimum of the potential energy . instead , the distribution function of the core particles progressively approaches that of a fully degenerate fermi gas  @xcite , @xmath90 where @xmath91 is the effective fermi energy . the final stationary state of the cluster will then correspond to a cold core surrounded by a high energy diffuse halo , @xmath92 where @xmath93 is the energy of the one particle resonance . the parameter @xmath94 and the effective fermi energy @xmath91 are determined using the conservation of particle number and energy . the extent and the location of the parametric resonance can be calculated using the canonical perturbation theory  @xcite . in fig . [ fig4 ] we show the poincar section of a test particle @xmath95 moving under the action of an oscillating potential calculated using the envelope equation  ( [ envelope0 ] ) , @xmath96 with , @xmath97 where @xmath98 is the modulus of the test particle angular momentum , conserved by the dynamics , and @xmath99 is fixed at its initial value @xmath100 .    ) with 10 different initial conditions , plotted when the envelope is at its minimum . ( b ) n - particle simulation for a non - virial system with @xmath101 . an excellent agreement is found between the extent of the halo in n - particle simulation and the one particle resonant orbit shown in the poincar plot . , width=340 ]    the resonant orbit is the outermost curve of the poincar plot fig . the first resonant particles move in an almost simple harmonic motion with energy @xmath102 , where @xmath103 is the intersection of the resonant trajectory with the @xmath104 axis . empirically we find that the location of the one particle resonance @xmath103 for values of @xmath105 is very well approximated by a simple expression  @xcite @xmath106 as the relaxation proceeds , the oscillating core becomes progressively colder , while a halo of highly energetic particles is formed . as more and more particles are ejected from the core , their motion becomes chaotic , and a halo distribution becomes smeared out . similar to what happens for magnetically confined plasmas , we find that the distribution function of a completely relaxed halo is very well approximated by the heaviside step function @xmath107 . for notational simplicity , from now on , we will drop the over - bar on the distribution function @xmath108 , but it should always be kept in mind that @xmath49 is stationary only within the coarse graining procedure described above . in order to obtain the density and the velocity distribution after the ss state is achieved , we solve the poisson equation  ( [ poisson0 ] ) @xmath109 with the constraints  ( [ eq conserv1 ] ) and ( [ conserv2 ] ) . since the initial mass distribution has the azimuthal symmetry , the potential must have the form @xmath110 substituting this into poisson equation and noting that @xmath111 we obtain @xmath112\\ % \epsilon_f+\frac{\chi}{\eta}(\epsilon_r-\epsilon_f)+c_1\>j_0(r^ * ) \\ \psi_{halo}(r)&=&\epsilon_r+c_2\>j_0(r^{**})+c_3\>y_0(r^{**})\\ \psi_{out}(r)&=&\ln \left(r\right ) \;,\end{aligned}\ ] ] where @xmath113 and @xmath114 are the bessel functions of the first type and of order @xmath115 ; @xmath116 is the core radius ; @xmath117 and @xmath118 . the integration constants @xmath119 and the value of @xmath116 can be determined using the continuity of the potential and of the gravitational field . the parameter @xmath94 can then be obtained using the conservation of energy . once @xmath119 are calculated , see appendix  [ constants ] , we are left with just two equations for @xmath116 and @xmath94 , @xmath120 where prime denotes the derivative with respect to @xmath9 and @xmath121 ^ 2}{4\eta j_0\left(r_c^{*}\right)}\;.\\ % \left(1-\frac{\chi}{\eta}\right)\left(\epsilon_r-\epsilon_f\right ) % \left[\pi^2\chi(\epsilon_r-\epsilon_f)r_c^2-\frac{r_c\>\psi'_{core}(r_c)}{2}\right]+\frac{\epsilon_r}{2}\\ % { \cal n}(r_c,\chi)&=&r_r \psi'_{out}(r_r)\end{aligned}\ ] ]    this completely determines the distribution function of the final stationary state achieved by a self - gravitating system when its initial distribution deviates from the virial condition . in fig . [ fig5 ] , we compare the predictions of the theory with the molecular dynamics simulations . an excellent agreement is found without any adjustable parameters .    . solid line is the theoretical prediction obtained using the distribution function , eq . ( [ eq corehalo ] ) , and the points are the result of molecular dynamics simulation with @xmath47 particles.,width=340 ]    finally , we explore the life time @xmath3 of a qss of a self - gravitating system with a finite number of particles . to do this we define the crossover parameter @xmath122 ^ 2{\rm d}r}\ ] ] where @xmath123 is the number density of particles inside shells located between @xmath9 and @xmath124 at each time of simulation @xmath10 and @xmath125 where @xmath126 is the stationary distribution given by eq . ( [ eq corehalo ] ) . the dynamical time scale is defined as @xmath127 . in fig . [ fig7]a we plot the value of @xmath128 for systems with different number of particles . [ fig7]b shows that if we scale the time with @xmath129 , where @xmath130 , all the curves collapse onto one universal curve , showing the divergence of the crossover time in the thermodynamic limit . it is interesting to note that for a hamiltonian mean - field ( hmf ) model , @xmath3 was found to diverge with the exponent @xmath131  @xcite , while for a virial 3d self - gravitating system the exponent was found to be @xmath132 . unfortunately , at the moment there is no theory which allows us to predict this exponents _ a priory_.     for different number of particles in the system . after relaxing into the qss , the system crosses over to mb distribution after time @xmath3 . inset ( a ) shows that the relaxation to core - halo state takes approximately @xmath133 and does not depend on the number of particles in the system . when the time is scaled with @xmath3 all the data in ( a ) ( for large times ) collapse onto one universal curve ( b).,width=188 ] we have studied the thermodynamics of 2d self - gravitating system in the microcanonical ensemble . it was shown that the gravitational clusters containing finite number of particles relax to the equilibrium state characterized by the mb distribution . prior to achieving the thermodynamic equilibrium , however , these systems become trapped in a quasi - stationary state , where they stay for time @xmath134 , which diverges as @xmath135 for large @xmath0 . thus , in the limit @xmath2 at fixed total mass @xmath5 , thermodynamic equilibrium can not be reached in a finite time . a new approach , based of the conservation properties of the vlasov dynamics and on the theory of parametric resonances , is formulated and allows us to quantitatively predict the one particle distribution function in the non - equilibrium stationary state . finally , it is curious to consider what will happen to a self - gravitating system in a contact with a thermal bath  the canonical ensemble . in appendix [ virial ] , it is shown that for a 2d self - gravitating system a stationary state is possible , if and only if , @xmath136=1/2 , i.e. when the kinetic temperature is @xmath137 . if such system is put into contact with a thermal bath which has @xmath138 there will be a constant heat flux from the reservoir into the system . this heat will be converted into the gravitational potential energy  since the kinetic energy is fixed by the virial condition  making the cluster expand without a limit . conversely if the bath temperature is @xmath139 , the heat flux will be from the system into the bath . again , since the system can only exist in a stationary state if @xmath137 , the energy for the heat flux can come only from the gravitational potential . in this case the gravitational cluster will contract without a limit , concentrating all of its mass at the origin . thus , in the canonical ensemble no thermodynamic equilibrium is possible , unless the reservoir is at exactly @xmath137 . we hope that the present work will also help shed new light on the collisionless relaxation in 3d self - gravitating systems . unfortunately the 3d problem is significantly more difficult , since besides the core - halo formation , one must also account for the particles evaporating to infinity . this work was partially supported by the cnpq , inct - fcx , and by the us - afosr under the grant fa9550 - 09 - 1 - 0283 . the gravitational potential energy @xmath140 of a system is @xmath141 where the integration extends over all space . unlike in 3d , the gravitational potential of a 2d system diverges at infinity . therefore , some care must be taken with the limits . performing the integration by parts we obtain @xmath142 where @xmath143 is the radius of the bounding sphere . from eq . ( [ gausscond ] ) , @xmath144 , and the total energy is given by @xmath145 the divergence in the last term is common to all energy calculations in 2d . for example , if at @xmath6 the system is distributed with a water - bag distribution eq . ( [ eq waterbag ] ) , its energy is @xmath146 the gravitational self - energy in 2d is always divergent . however , since this divergence is always the same , it can be easily renormalized away . we simply add an infinite constant , @xmath147 , to all gravitational self - energies . the renormalized ( finite ) energy @xmath148 of a self - gravitating system is then @xmath149 which for a water - bag distribution eq . ( [ eq waterbag ] ) becomes @xmath150 the hamiltonian of a general self - confined system is given by @xmath151 where @xmath152 is the momentum of particle @xmath95 , and @xmath153 is the interaction potential . the virial function @xmath154 is defined as @xmath155 taking the time derivative , and using hamilton s equations we obtain @xmath156 where , @xmath157 if @xmath158 is an homogeneous function of order @xmath55 @xmath159 euler s theorem requires that @xmath160 for a stationary state @xmath161 , and we obtain the usual result @xmath162 , where @xmath163 is the mean kinetic energy .    in two dimensions , @xmath158 is a sum of logarithms and euler equation does not apply directly . however , we can still derive a 2d virial theorem by writing the inter - particle interaction potential as @xmath164 , which is a logarithm plus an infinite constant . this is a homogeneous function of order p=0 , so we can use the euler theorem to write @xmath165 note that since the right hand side of this equation only contains derivative of the potential , this expression is equally valid for @xmath166 and for the logarithmic potential . substituting eq . ( [ vir ] ) into eq . ( [ eq a1 ] ) , we arrive at 2d virial theorem  @xcite . @xmath167 in the thermodynamic limit , and after rescaling the velocity to put everything into adimensional form , we obtain the result @xmath168 , quoted in section  [ virial ] . we define the @xmath169 radius of the mass distribution as @xmath170 . deriving twice with respect to time we obtain @xmath171 this reduces to @xmath172 where @xmath173 is the _ emittance _ which commonly appears in plasma physics  @xcite . the last term can be simplified using the poisson equation  ( [ poisson0 ] ) in its dimensionless form , @xmath174}\nonumber\\ & = & \lim_{r_0 \rightarrow \infty}{-\frac{1}{2}\left(r{\partial\psi\over\partial r}\right)^2\bigg|_{r = r_0 } } \ , , \label{aprforce}\end{aligned}\ ] ] which can be obtained directly using eq . ( [ gausscond ] ) , @xmath175 for a water bag initial distribution  ( [ eq waterbag ] ) we define the envelope radius as @xmath176 , so that for @xmath6 @xmath80 , and rewrites  ( [ envelope1 ] ) as @xmath177 this is the envelope equation . if initially , @xmath178 then @xmath179 and the envelope will not oscillate . this is precisely the virial condition . integrating the core - halo distribution function over velocity , the dimensionless poisson equation  ( [ eq poisson ] ) becomes @xmath180 we define @xmath181 for @xmath182 , @xmath183 for @xmath184 , and @xmath185 for @xmath186 . changing variables @xmath187 and @xmath188 , we can rewrite  ( [ eq poisson2 ] ) as @xmath189 the solution of the first two of these equations can be written in terms of the bessel functions of first type and of order @xmath115 , @xmath190 the last equation is solved by @xmath191 where @xmath192 are the integration constants . the regularity of solution at the origin and eq.([gausscond ] ) require that @xmath193 , @xmath194 and @xmath195 , respectively . the potential reduces to @xmath112\\ % \epsilon_f+\frac{\chi}{\eta}(\epsilon_r-\epsilon_f)+c_1\>j_0(r^ * ) \\ \psi_{halo}(r)&=&\epsilon_r+c_2\>j_0(r^{**})+c_3\>y_0(r^{**})\\ \psi_{out}(r)&=&\ln \left(r\right ) \;,\end{aligned}\ ] ] where we have defined @xmath116 such that @xmath111 . the others requirements to continuity of the potential and its derivative are , @xmath196 solving these equations yields the integration constants @xmath119 as a function of @xmath197 and the parameters @xmath198 , @xmath199 the remaining equation of continuity of @xmath200 at @xmath116 and the conservation of energy will determine @xmath116 and @xmath94 , see eq.([econ ] ) . 99 j. w. gibbs , _ collected works _ , longmans , green and co. , ny ( 1928 ) . j. barr and d. mukamel , and s. ruffo , phys . lett . * 87 * , 030601 ( 2001 ) . d. lynden - bell , mon . not . . soc . * 136 * , 101 ( 1967 ) . w. thirring , zeitschrift fr physik , * 235 * , issue 4 , pp.339 - 352 ( 1970 ) . d. lynden - bell and r. m. lynden - bell , , mon . not . . soc . * 181 * , 405 ( 1977 ) . d. mukamel , and s. ruffo , and n. schreiber , phys . * 95 * , 240604 ( 2005 ) . a. ramrez - hernndez and h. larralde and f. leyvraz phys . lett . * 100 * , 120601 ( 2008 ) . klein , phys . rev . * 97 * , 1446 ( 1955 ) . f. bouchet and t. dauxois phys . e * 72 * , 045103 ( r ) ( 2005 ) . kavita jain et al j. stat . ( 2007 ) p11008 t. padmanabhan , physics reports * 188 * , 285 ( 1990 ) ; p .- h . chavanis , int . j. mod . b * 20 * , 22 , 3113 - 3198 ( 2006 ) . y. levin , r. pakter , and t. n. teles , phys . lett . * 100 * , 040604 ( 2008 ) . f. b. rizzato , r. pakter , and y. levin phys . e * 80 * , 021109 ( 2009 ) . chavanis , physica d 200 , 257 ( 2005 ) . j. miller , phys . lett . * 65 * , 21372140 , ( 1990 ) . p. h. chavanis and j. sommeria , journal of fluid mechanics , * 314 * , pp 267 - 297 . a. venaille and f. bouchet , phys . lett . * 102 * , 104501 ( 2009 ) . a. campa , t. dauxois , and s. ruffo , phys . 480 , 57 ( 2009 ) . y. levin , r. pakter , and f.b . rizzato phys . e. * 78 * , 021130 ( 2008 ) . m. joyce , b. marcos , f. s. labini , j. stat . theory and experiment 2009 , p04019 ( 2009 ) . y. levin , rep . . phys . * 65 * , 1577 ( 2002 ) . j. ostriker , astrophys . j. 140 , 1056 ( 1964 ) ; j.j . aly , phys . e 49 , 3771 ( 1994 ) ; j.j . aly and j. perez , phys . e 60 5185 ( 1999 ) w. braun and k. hepp , comm . math . phys . * 56 * , 101 ( 1977 ) ; a. antoniazzi , f. califano , d. fanelli , and s. ruffo , phys . rev . lett . * 98*. 150602 ( 2007 ) . m. antoni and s. ruffo , phys . e * 52 * , 2361 ( 1995 ) t. m. rocha filho , m. a. amato , and a. figueiredo j. phys . a : math . theor . * 42 * 165001 ( 2009 ) . teles , r. pakter , and y. levin , applied phys . 95 , 173501 ( 2009 ) . r. l. gluckstern , phys . 73 * , 1247 ( 1994 ) . t. p. wangler , k. r. crandall , r. ryne , and t. s. wang , phys . rev . beams 1 , 084201 ( 1998 ) . m. reiser , _ theory and design of charged particle beams _ , ( wiley - interscience , 1994 ) ; r.c . davidson and h. qin , _ physics of intense charged particle beams in high energy accelerators _ ( world scientific , singapore , 2001 ) . chavanis , phys . e * 73 * , 066103 ( 2006 ) . yamaguchi , j.barr , f.bouchet , t.dauxois , and s.ruffo , physica a ( 2004 ) 337 , 36 - 66 .
Abstract:

Answer: we study , using both theory and molecular dynamics simulations , the relaxation dynamics of a microcanonical two dimensional self - gravitating system . after a sufficiently large time , a gravitational cluster of @xmath0 particles relaxes to the maxwell - boltzmann distribution . the time to reach the thermodynamic equilibrium , however , scales with the number of particles . in the thermodynamic limit , @xmath1 at fixed total mass , equilibrium state is never reached and the system becomes trapped in a non - ergodic stationary state . an analytical theory is presented which allows us to quantitatively described this final stationary state , without any adjustable parameters .