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: galaxy clusters are dark matter haloes which form mainly through the process of gravitational collapse , and so they are expected to be self - similar : small clusters or galaxy groups must be scaled - down versions of the more massive systems . self - similarity predicts a series of relationships between different cluster observables such as x - ray luminosity @xmath9 , mass @xmath12 , temperature @xmath13 and galaxy velocity dispersion @xmath10 , which are the same in all systems ( kaiser 1991 ) . the radial velocity dispersion @xmath10 of the galaxies in a cluster probes the depth and shape of the potential well , assuming that the luminous matter traces reasonably well the dark matter in clusters . moreover , the intracluster gas emits x rays through a process of thermal bremsstrahlung and its bolometric luminosity @xmath9 is found to be strongly correlated with @xmath10 , as both gas and galaxies share a common potential well ( solinger & tucker 1972 ) . also , since the x - ray emission can be well modelled by thermal emission of a hot , optically thin plasma , @xmath9 and the plasma temperature @xmath13 must be correlated . finally , because galaxies are embedded in the intracluster medium , a correlation is also expected between the gas temperature @xmath13 and the galaxies velocity dispersion @xmath10 . self - similar models consider that the only energy source in the cluster comes from the gravitational collapse , predicting the following scaling relations : @xmath14 . whereas there seems to be a general agreement between different measurements that @xmath15 , the measurement of @xmath0 has given somewhat contradictory results so far . some authors have found that @xmath16 indeed , although with quite large measurement errors or rather small data samples , while others find slopes larger than 4 ( see table [ literature ] ) . part of the differences in the results could come from different ways of selecting the sample , with a preference for regular clusters in some of these surveys . it has also been suggested that clusters and groups do not follow the same @xmath17 scaling relation , the latter being flatter than the former ( e.g. mahdavi et al . 2000 , xue & wu 2000 , dellantonio et al . but there are other measurements contrary to that conclusion ( mahdavi & geller 2001 , mulchaey & zabludoff 1998 ) .    for more distant clusters ( @xmath18 between @xmath19 and @xmath20 ) there is some evidence that the slope is also @xmath21 ( borgani et al . 1999 , girardi & mezetti 2001 ) , although only small samples are available at the moment , and more data are needed to reduce the error bars . there is a clear need of new measurements performed on large samples ( more than 100 clusters ) , selected from a homogeneous dataset . the largest cluster ( not group ) samples available up to date , on which this kind of study has been performed , are literature compilations combining data from different authors and/or instruments ( for example xue & wu 2000 , mahdavi & geller 2001 ) . the present work comes to partially fill this gap , as we have a large homogeneous subsample of clusters selected from reflex , a flux - limited x - ray selected cluster catalogue built under well defined selection criteria . the structure of this paper is as follows . in section [ thedata ] we give a brief introduction to the reflex catalogue and the way in which the subsample in this paper has been selected . section [ measurement ] illustrates in detail the process of the @xmath9 and @xmath10 measurements , which lead to the the fitting of the @xmath1 in section [ lxsigma ] . in section [ nosubst ] we propose the removal of those clusters in the sample which are multiple systems . in section [ discussion ] we analyse possible biases due to the nature of the selected sample . a study on the @xmath22 and @xmath23 relation follows in section [ lxt ] , and we finish with a summary in section [ summary ] . we are assuming h@xmath24 km s@xmath3 mpc@xmath3 in a flat universe with @xmath25 . the reflex ( rosat - eso flux - limited x - ray ) cluster survey ( bhringer et al 2001 , bhringer et al 2003 in preparation , guzzo et al . 2003 in preparation ) has identified and measured successfully the redshift for all southern galaxy clusters ( at galactic latitude @xmath26 ) down to a flux limit @xmath27 erg s@xmath3 @xmath28 in the rosat all - sky survey . this sample has already provided us with results concerning the measurement of the cluster x - ray luminosity function at low redshift ( bhringer et al 2002 ) , the power spectrum ( schuecker et al . 2001a ) , the spatial correlation function ( collins et al . 2000 ) , the measurement of @xmath29 and @xmath30 ( schuecker et al . 2003a ) or the measurement of the equation of state parameter @xmath31 of the dark energy ( schuecker et al . 2003b ) , among others . a detailed description of the sample construction can be found in bhringer et al . ( 2001 ) . reflex is an x - ray selected sample and the first compilation of the reflex catalogue includes 452 clusters down to a flux limit of @xmath32 erg s@xmath3 @xmath28 in the rosat energy band from 0.1 to 2.4 kev . the survey covers a total area in the sky of 4.24 ster . two regions corresponding to the magellanic clouds and the galactic plane have been excluded because of the difficulties in identifying clusters of galaxies in these densely crowded regions , as well as some strips of the sky which were not observed by rosat due to power - off of the detector when the satellite was crossing the earth s radiation belts ( see figure 2 in bhringer et al 2001 for a survey exposure map ) . optical follow - up has been performed through long - slit and multi - object spectroscopy at the 1.5 m , 2.2 m and 3.6 m eso telescopes at la silla ( chile ) during a series of observing campaigns , which started back in 1991 . we complement our sample of galaxy redshifts by looking into the literature ( using the nasa / ipac extragalactic database , ned ) and adding those redshifts already published for the previously known clusters . the cluster membership of each individual galaxy has been decided based on the search for velocity peaks in the velocity distribution of the galaxies observed in the field . all galaxies were assigned the same weight in the computation of the cluster s median redshift because in most cases not enough photometric material is available to give proper weights ( i.e. to identify a real cd galaxy ) . anyway , bright cluster galaxies ( if present ) were used to decide on the main redshift of the system when it showed multiple velocity peaks in the redshift histogram along the line of sight . in other cases the choice among multiple velocity peaks was decided taking into account the correlation between the peak position and the x - ray emission centre and also the strength of the peak . finally , other clusters ended up being split into different systems when multiple x - ray emission centres associated to them were detected .    we have selected from reflex those clusters with more than seven measured galaxy redshifts ( see fig . [ z_per_clu ] ) , including only galaxies at a maximum projected distance from the cluster center of 0.5 h@xmath33 mpc , due to the lack of data at larger distances from the cluster centre . outliers in the velocity space were rejected by applying the 3 sigma clipping method : galaxies with velocity differences with the cluster centre velocity larger than three times the measured velocity dispersion were not used in the subsequent iterations to avoid possible contamination by foreground / background objects . by following this procedure we end up with 171 clusters out of the 452 in the whole catalogue . it is a low - redshift sample , the median redshift being 0.076 . only 14 clusters have @xmath18 above 0.2 ( see fig . [ z_histo ] ) .     is measured . ] counts in the rosat images were measured using the growth curve analysis ( gca ) and then converted into fluxes , as described in bhringer et al . luminosities have first been transformed into the rosat rest - frame band ( 0.1 - 2.4 kev ) ( see bhringer et al . 2002 for details ) and afterwards into bolometric luminosities . we assumed a mekal spectral model with an icm metallicity of 0.3 z@xmath34 ( anders & grevesse 1989 ) . we used xspec ( arnaud 1996 ) to obtain the relation between the luminosity in the 0.1 - 2.4 kev band and the bolometric one , as a function of the cluster temperature(see bhringer et al . 2003 for a detailed conversion table ) , which was computed from the velocity dispersions assuming the empirical relation between @xmath10 and t measured by girardi et al . ( 1996 ) . the radial velocity dispersions @xmath10 have been computed by means of the biweight estimator of central location and scale ( within rostat by beers , flynn & gebhardt 1990 ) . this estimator was first suggested by tuckey ( in andrews et al . 1927 ) as a better way to study non - gaussian or contaminated normal distributions than the gaussian estimators ( mean and standard deviation ) . the biweight location estimator belongs to the family of estimators known as @xmath12 estimators of location . it works by minimizing a function of the deviations of each observation from the estimate of location and has the advantage that it does not assume an underlying gaussian distribution , which in our case means to have a perfectly relaxed cluster and might not be the case in most of the sample . it is less affected by points in the wings of the distribution and is robust for a broad range of non - gaussian underlying populations . we used the biweight estimator of central location and the gapper estimator of scale ( based on the gaps between order statistics , wainer & thissen , 1976 ) when the number of redshifts available in the cluster was @xmath35 . for @xmath36 the biweight estimator was chosen for both location and scale . errors were obtained in all cases by jacknifing of the biweight . the choice of the different estimators has been done on the basis of the tests carried out by beers et al . ( 1990 ) for samples with different number of data points . the usual cosmological correction and the correction for velocity errors ( danese , de zotti & tullio , 1980 ) were also applied . the whole procedure was performed iteratively until the results converged . the clusters with the larger errors in the velocity dispersion estimate are among the ones for which @xmath37 . one source of concern is the fact that we are computing the @xmath10 using galaxies which are quite close to the cluster centre , at distances less than 0.5 h@xmath38 mpc ( median distance is 0.04@xmath38 mpc ) . girardi et al . ( 1996 ) find that the value of @xmath10 is a function of the distance to the cluster centre , and it reaches a stable value when galaxies at a projected distance of @xmath39 1 h@xmath38 mpc are considered . at smaller distances , the @xmath10 does not show a constant behaviour : sometimes it raises up , sometimes it goes down , and in some cases it remains stable . so , in a large sample like ours we expect that this effect will cancel out . from our current data , it is not possible to determine the @xmath10 at clustercentric distances larger than 0.5 h@xmath38 mpc . a power law was fit to the @xmath40 relation by means of the orthogonal distance regression method ( odr , boggs et al . 1992 ) , which takes into account errors on both variables . we took the data in units of @xmath41 erg s@xmath3 ( @xmath42 ) and velocity dispersions in units of @xmath43 km s@xmath3 ( @xmath44 ) , that is , we put the origin of coordinates more or less at the center of the datapoint cloud to help the model find the right y - axis intercept . we will use these units throughout the paper , but for table [ literature ] , where we quote @xmath9 in erg s@xmath3 and @xmath10 in km s@xmath3 for comparison purposes with the works referenced there . the best fit to @xmath45 is found to be ( see figure [ odr_data_all ] )    @xmath46    the spearman rank correlation coefficient for this sample is @xmath47 and the probability that a random distribution had this value of @xmath48 or larger by chance is @xmath49 .     vs. radial velocity dispersion @xmath10 for a reflex sample of 171 clusters . the linear fit gives a log - slope of @xmath50 . and the dashed lines are the 1-sigma errors . filled circles correspond to clusters at redshift @xmath51 , open triangles are clusters with @xmath52 and open circles are clusters at @xmath53 . ] the distribution of errors in the sample has been assumed to be normal in the odr fitting . to test for this assumption , we have performed a maximum likelihood analysis with sigma clipping to find which is the best fitting gaussian to the actual error distribution in x. as errors in x are far larger than errors in y , we can neglect the effect of the latter . we find that the best gaussian has central location @xmath54 and standard deviation @xmath55 ( figure [ gaussianity_x_7]).the error distribution along the x axis is thus reasonably gaussian .     and standard deviation @xmath55 . ] we can also perform a crude estimate of the intrinsic dispersion of the sample , @xmath56 , @xmath57 in error units . @xmath58 is the best fitting gaussian @xmath59 parameter as found above for the current error distribution , that is @xmath60 . @xmath61 is the @xmath59 parameter of the error distribution in the ideal case in which only the measurement errors are responsible for the dispersion in the data , which would be a gaussian with @xmath62 . we find @xmath63 in error units , that is , the intrinsic dispersion of the sample is slightly larger than the dispersion due to the errors in the data . to see what this means in terms of physical units , let us assume that a typical value of the error in the velocity dispersion is the median of the error distribution , @xmath64 km s@xmath3 . the intrinsic dispersion is then    @xmath65    finally , to test for the stability of the fitting , we repeated it by removing those points with an error in the velocity dispersion larger than 30 percent of their value . we find @xmath66 thus showing the robustness of the fitting against the presence of large errors in @xmath10 . we may also worry about the small number of redshifts ( 7 at minimum ) with which we are computing the @xmath10 in some clusters . we have repeated the fit , this time taking only clusters for which @xmath10 had been determined from at least 30 individual redshifts . in this way , we are left with a sample of only 57 clusters . the best fit relation is @xmath67 which is in agreement at the one sigma confidence level with the fit obtained for the whole sample . figure [ odr_data_30 ] shows the best fit found . the best fitting gaussian to the error distribution in x is centered at @xmath68 and has standard deviation @xmath55 . this gives a value for the intrinsic dispersion of @xmath69 km s@xmath3 . vs. radial velocity dispersion @xmath10 for a reflex sample of 57 clusters for which more than 30 galaxy redshifts entered the @xmath10 computation . the linear fit gives a log - slope of @xmath70 . and the dashed lines are the 1-sigma errors . filled circles correspond to clusters at redshift @xmath51 , open triangles are clusters with @xmath52 and open circles are clusters at @xmath53 . ] self - similarity in the purely gravitational hierarchical scenario of structure formation implies that for galaxy clusters @xmath14 . although several authors have found a good agreement between this prediction and their measurements for the @xmath17 relation , others have found steeper values ( see table [ literature ] ) .    [ cols="^,^,^,^,^,^",options="header " , ]     our value of @xmath71 is in good agreement with the slope measured by @xcite,@xcite , @xcite and @xcite . as for the intercepts , our result is compatible with @xcite and @xcite at about the one - sigma confidence level .    from fig . [ sigmas_histo ] one can see that our sample is mostly populated by clusters and groups are practically absent , if we consider as `` groups '' those systems with @xmath72 340 km s@xmath3 ( following the criterium used by mahdavi et al . more specifically , only 11 out of the 171 systems are groups . some authors ( e.g. helsdon & ponman 2000 , mahdavi et al . 2000 , xue & wu 2000 , dellantonio et al . 1994 ) have found that the scaling laws are different for clusters and groups . in particular , they find that clusters show a steeper slope in the @xmath17 relation . unfortunately , e can not derive a meaningful @xmath17 relation for groups from our sample , due to the small number of groups observed . the complexity involved in the measurements of @xmath10 and @xmath9 ( because of the presence of substructures , cooling flows , and other physical phenomena ) , the largely unknown intrinsic dispersion of their relationship , and the use of different fitting methods could explain the different results reached by the different works . for example , @xcite illustrate the different results obtained by using odr regression and ordinary least squares ( ols ) regression in the analysis . substructures in clusters constitute a problem in the radial velocity dispersion estimates , as the different components show different @xmath10 and an estimate of a single @xmath10 for the cluster may not be meaningful . we can not make a reliable study of the substructures in our clusters as in many cases we do not have enough redshifts to perform it . the sensitivity for substructure detection depends strongly on the number of galaxies with measured redshift available for each cluster . in our sample , though , we have large differences in the number of galaxies per cluster and trying to detect substructure is not feasible .    in a paper by schuecker et al ( 2001b ) they identify the clusters in the reflex sample which show multiple components in the x - ray emission . when we take those clusters out from our sample , we end up with 123 clusters . in this case we have    @xmath73    which is again in very good agreement with the one obtained for the whole sample . this result shows that our sample is not largely affected by the problems of substructure , probably because we are selecting galaxies very close to the x - ray peak emission , that is , to the central potential well , where one would expect that the strength of the gravitational potential will erase any substructures which might be present . even in the case of redshifts taken from the literature ( see section [ in_redshifts ] ) we have not considered galaxies farther than 0.5 h@xmath38 mpc , to be consistent with the characteristics of our own data . the galaxy redshifts used in this work come from three different sources . when a cluster was already known and well studied , redshifts were taken from the literature , possibly from different authors too . when a cluster was already known but had little data available in the literature , we took our own redshifts in order to secure and/or improve the cluster redshift . finally , in the case of clusters discovered in this project , a complete optical follow - up has been carried out . it is thus not unreasonable to think that such a variety of redshift sources might introduce some undesirable effect particularly on the @xmath10 estimate and therefore in the final @xmath0 relation . to check whether this is the case , we have restricted ourselves to a sample of 38 clusters for which more than 7 redshifts per cluster were available , all of them being redshifts measured by us from our own data . this is then to be considered a completely homogeneous sample in terms of both luminosity and redshift , as both quantities have been measured in the same way for all 38 clusters . the result of the odr fitting is @xmath74 . in this case the fit is quite poor , probably because this sample suffers from important selection effects . first , as we have now a small number of data points we should worry about the problems of measuring @xmath10 with galaxies at small projected distance from the cluster centre , as explained in section [ sigma ] .    the minimum number of galaxy redshifts set to measure the @xmath10 may have a role too , specially in this case where not too many clusters are available for the fitting . when we consider the whole sample , clusters with smaller number of individual galaxies are given much less weight into the fit than those with a larger number of observed galaxies . that is not the case here , as all of the clusters show small numbers ( only 3 have more than 15 measured @xmath18 ) , so the differences in the weights are not so dramatic .      of course , we would like to use all 452 clusters in the reflex catalogue to estimate the @xmath0 relation . this can not be done , unfortunately , as the survey was targeted to measure cluster redshifts and not velocity dispersions in these clusters . while multi - object spectroscopy observations were used for a few distant , compact clusters ( due to the small field of view , 5 arcmin , of the efosc spectrograph on the eso 3.6 m telescope , see e.g. guzzo et al . 1999 ) , providing around 15 redshifts per cluster , most of the new reflex cluster redshifts come from ( multiple ) single slit observations which do not deliver enough galaxy redshifts to compute a velocity dispersion , although the main goal was always to have at least 5 spectra per target . we have compared the @xmath75 distribution for the whole catalogue with the one corresponding to the 171 cluster subsample used in this work . figure [ lx_z ] shows as filled circles the clusters from the subsample , and as open circles the ones from the whole reflex catalogue . one can see that both samples are approximately equally distributed . the results from a kolmogorov - smirnov test are compatible with this conclusion : we find a maximum distance value between the two samples ( the entire catalogue and the 171 cluster sample ) of 0.11 , the probability of this value being larger if coming from a random sample being 0.22 . distribution for all clusters in the reflex catalogue ( open circles ) and those from the 81 cluster subsample used in this work ( filled circles ) . ]      a bias is expected when measuring the @xmath0 relation from a flux - limited sample . specifically , at the flux cut - off level only the brightest clusters will be considered , for a given velocity dispersion . as our sample has actually been drawn from a flux - limited sample , we may worry about whether we are introducing such a bias towards larger luminosities . this bias can be avoided by considering a volume - limited subsample . therefore , we have repeated the analysis for the largest volume - limited sample that we could construct from the 171 cluster sample . it was built by taking all clusters with luminosity @xmath76 h@xmath77 erg s@xmath3 , which corresponds to a redshift limit of @xmath78 . to this redshift the comoving volume surveyed by reflex is @xmath79 h@xmath80 mpc@xmath81    we have a total of 51 `` useful '' clusters ( according to the selection criteria outlined in section [ thedata ] ) out of the 88 reflex clusters in this volume . the odr fitting gives @xmath82 , only compatible at the three - sigma error level with the one obtained from the whole sample , which is an indication of some bias present in the total sample . another important scaling relation in galaxy clusters is the x - ray luminosity - temperature @xmath83 . again , measurements show that it does not follow the scaling relations predicted by pure gravitational collapse models ( @xmath85 ) . the observed @xmath22 relation is closer to @xmath86 ( e.g. novicki et al . 2002 , fairley et al . 2000 , xue & wu 2000 , arnaud & evrard 1999 , reichart et al . 1999 ) , which seems to indicate the presence of additional non - gravitational heating sources in the icm . @xmath23 has also been studied in the literature . there seems to be a consensus in that @xmath87 ( see for example xue & wu 2000 , girardi et al . 1996 , bird , mushotzky & metzler 1995 , lubin & bahcall 1993 ) , again not in agreement with what would be expected in a pure gravitational collapse scenario , @xmath88 . we have studied these relations within our cluster sample to see to which degree our data confirm or otherwise those results . asca temperatures from two different catalogues ( horner 2001 and ikebe et al . 2002 ) were used for 54 clusters included in our 171 cluster sample . ikebe et al.s measurements were preferred whenever both works had data for the same cluster , as they used a two - temperature model ( instead of horner s one component model ) for the isothermal plasma , allowing a multiphase intra - cluster medium , thus obtaining more accurate t estimates . again , we performed the odr fitting in the log - log space , obtaining the following best - fit relations ( see also figures [ lx_t ] and [ vdisp_t ] ) :    @xmath89    @xmath90     for a subsample of 54 clusters drawn from the total sample , for which asca temperatures were available . the solid line is the odr best - fit model , with one - sigma errors marked as dashed lines . ]     for a subsample of 54 clusters drawn from the total sample , for which asca temperatures were available . the solid line is the odr best - fit model , with one - sigma errors marked as dashed lines . ] the @xmath22 estimate that we find is in good agreement ( at one - sigma confidence level ) with the slope and intercepts previously measured by fairley et al . 2000 and arnaud & evrard 1999 .    for the @xmath23 relation , we are finding here a larger slope than found by other authors ( which is @xmath91 ) . to check that we are not biased for including clusters with substructure in the sample , we select a new set of clusters from the sample in section [ nosubst ] for which temperature data are available . we have 33 systems , for which we find    @xmath92    in agreement with the previous result . furthemore , we can select only clusters with the most reliable @xmath10 determination , that is , those with @xmath10 measured from more than 30 galaxy redshifts . in doing so , we end up with a sample of 16 clusters , giving    @xmath93      it is interesting to probe these relations also in a volume - limited sample . from the one we constructed in sect . [ vlted ] , we select those clusters for which temperature data are available , ending up with a 28 cluster sample . on this we find that    @xmath94    @xmath95    the @xmath96 is still the same , but the @xmath97 is flatter than what is found in the 54 cluster sample , steeper than what is expected from self - similar models ( at one - sigma confidence level ) , and in good agreement with what is found by other authors ( xue & wu 2000 , girardi et al . 1996 and lubin & bahcall 1993 , for example ) .    in conclusion , we have derived an @xmath22 similar to what is found by previous studies , and which is in contradiction with self - similarity . but the @xmath97 we are measuring is not reliable , as it has shown to vary depending on the sample s selection criteria . the scatter and the large error bars in the velocity dispersions make it impossible to obtain a consistent relation . still , we can draw a useful conclusion , which is that self - similarity is ruled out in any case . we present an estimate of the @xmath98 relation computed with a new large and homogeneous sample of clusters from the reflex survey . all studies on the @xmath0 relation in large samples ( with more than @xmath99 clusters ) have been performed so far on compilations from the literature requiring the use of models to combine different data from different authors and/or instruments ( e.g. xue & wu 2000 , mahdavi & geller 2001 ) .    from the reflex catalogue we construct a sample of 171 clusters for which we can derive reliable estimates of @xmath10 , and then fit the @xmath100 relation , finding    @xmath101    in agreement with what is expected from self - similar models of cluster formation . the large size and homogeneity of the sample allows us to make an estimate of the intrinsic dispersion of the @xmath0 relation . we find that @xmath102 km s@xmath3 . measurements of @xmath9 and @xmath10 in clusters with multiple components may not be reliable . the identification in the literature of those clusters allows us to build a new sample where we remove those with significant substructure in their x - ray emission in schuecker et al . ( 2001b ) . this procedure leaves us with a sample of 123 clusters , on which we fit again the @xmath100 relation , finding    @xmath103    compatible at one - sigma confidence level with the slope and intercept that we find when using the 171 cluster sample , suggesting that the sample is not noticeably affected by clusters with multiple components . the fact that we are considering only galaxies close to the cluster centre may be a reason for that . we have also investigated the @xmath104 relation in a volume - limited sample , finding @xmath82 . this significantly flatter slope is an indication that some bias may be present on the whole sample due to the flux limit imposed on the data . the @xmath22 and @xmath97 scaling relations have also been explored , finding that @xmath7 and @xmath8 , in a subsample of 54 clusters for which @xmath13 was available from the literature . these results are consistent ( within the error bars ) with the one obtained for the @xmath0 , as we find @xmath105 when using these @xmath22 and @xmath23 to derive @xmath0 . the slope found in the @xmath106 relationship is steeper than the value predicted by a purely gravitational collapse model of cluster formation , and is in good agreement with previous measurements performed on other cluster samples . as for the @xmath107 , the slope we find is only compatible with the self - similar value at the three - sigma confidence level , and slightly larger than the value found by other authors . also , if we consider a volume - limited sample , the slope found is then in good agreement with previous studies , yet incompatible with self - similarity . the fact that this model fails to reproduce the measured slope is attributed to the contribution to the icm energy budget from other non - gravitational heating sources ( see e.g. bialek , evrard & mohr 2002 ) , like agns , supernovae , gas cooling ( pearce et al . 2000 ) , the presence of cool cores ( allen & fabian 1998 ) , or the effect of shocks in cluster formation ( cavaliere et al . 1997 ) . here we are faced with a paradoxical result : the slope of @xmath108 supports the self - similar scenario , whereas @xmath106 and @xmath23 do not , and all three relations are consistent with each other within one sigma confidence level . clearly , more data are required to clarify the situation and reduce the error bars . this is a reflection of the general situation in this topic in the literature : so far there is a general agrement on @xmath106 and @xmath107 , but there is not such an agreement on the @xmath0 . we thank c. altucci as her laurea thesis served as a starting point and guideline for this work . we are grateful to the whole reflex team which helped in the compilation of this catalogue . we are also indebted to mrs . yu xing , who provided us with the data needed to compute the bolometric luminosity out of the luminosity in the rosat band . we also thank m. girardi and s. borgani for stimulating discussions and suggestions . we are indebted to t. beers who gently provided us with the rostat routines , as well as to the authors of the odr package , p. boggs and j. rogers . ps acknowledges support by dlr under the grant no.50or970835 . this research has made use of the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration . this research has made use of nasa s astrophysics data system .
Abstract:

Answer: we present an estimate of the bolometric x - ray luminosity - velocity dispersion ( @xmath0 ) relation measured from a new , large and homogeneous sample of 171 low redshift , x - ray selected galaxy clusters . the linear fitting of @xmath1 gives @xmath2 erg s@xmath3 h@xmath4 . furthermore , a study of 54 clusters , for which the x - ray temperature of the intracluster medium t is available , allows us to explore two other scaling relations , @xmath5 and @xmath6 . from this sample we obtain @xmath7 and @xmath8 , which are fully consistent with the above result for the @xmath9-@xmath10 . the slopes of @xmath5 and @xmath6 are incompatible with the values predicted by self - similarity ( @xmath11 ) , thus suggesting the presence of non - gravitational energy sources heating up the intracluster medium , in addition to the gravitational collapse , in the early stages of cluster formation . on the other hand , the result on @xmath1 supports the self - similar model . [ firstpage ]    cosmology : observations  galaxies : clusters : general  galaxies : fundamental parameters


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 crossing of observational data from high redshift supernovae of type ia ( sne ia ) @xcite , cosmic microwave background ( cmb ) radiation @xcite , matter power spectra @xcite , x - rays from clusters of galaxies @xcite , and weak gravitational lensing @xcite strongly suggests that the present universe is dynamically dominated by a dark sector which is responsible for about @xmath0 of its total energy content . it is usually assumed that this dark sector has two different components : ( i ) dark matter , which is supposed to consist of weakly interacting massive particles ( wimps ) with zero effective pressure and ( ii ) dark energy , a mysterious entity which is equipped with a negative pressure . dark matter candidates include axions ( a particle present in the multiplet of grand unified theories ) and neutralinos ( light particles present in broken supersymmetric models ) , but none of these particles could be detected until now . the most natural dark energy candidate is a cosmological constant which arises as the result of a combination of quantum field theory and general relativity . however , its theoretical value is between 60 - 120 orders of magnitude greater than the observed value for the dark energy . an alternative to the cosmological constant is a self - interacting scalar field , known as quintesssence . for a brief but enlightening review of these and other proposals , see @xcite and references therein . there exists another route of investigations in which dark matter and dark energy are described within a one - component model . according to this idea , dark matter and dark energy are just  different faces " of a single , exotic fluid . to the best of our knowledge , the first proposal along this line was the chaplygin gas in its original and modified forms @xcite . however , this unified description of dark energy and dark matter , in spite of many attractive features , seemed to suffer from a major drawback : it predicted strong small scale oscillations or instabilities in the matter power spectrum , in complete disagreement with the observational data @xcite . ( on the other hand , the observed matter power spectrum corresponds to the baryonic matter distribution which does not exhibit strong oscillations @xcite , so that this point is still controversial . ) the apparently unrealistic predictions of the unified chaplygin gas type models are the result of an adiabatic perturbation analysis . it has been suggested that non - adiabatic perturbations may alleviate or even avoid this problem@xcite . this paper explores to what extent a viscous fluid can provide a unified description of the dark sector of the cosmic medium . the general influence of shear and bulk viscosity on the character of cosmological evolution has been studied , e.g. in @xcite , in the context of bianchi type i models . under the conditions of spatial homogeneity and isotropy , a scalar bulk viscous pressure is the only admissible dissipative phenomenon . the cosmological relevance of bulk viscous media has been investigated in some detail for an inflationary phase in the early universe ( see @xcite and references therein ) . however , as was argued in @xcite , an effective bulk viscous pressure can also play the role of an agent that drives the present acceleration of the universe . ( notice that the possibility of a viscosity dominated late epoch of the universe with accelerated expansion was already mentioned in @xcite ) . for a homogeneous and isotropic universe , the @xmath1cdm model and the ( generalized ) chaplygin gas models can be reproduced as special cases of this imperfect fluid description @xcite . moreover , in a gas dynamical model the existence of an effective bulk pressure can be traced back to a non - standard self - interacting force on the particles of the gas . while these investigations were performed for the homogeneous and isotropic background dynamics ( a study of the background dynamics which is similar to the setup of the present paper was recently performed in @xcite ) , a first perturbation theoretical analysis for a unifying viscous fluid description of the dark sector was performed in @xcite .    the bulk viscous pressure @xmath2 will be described by eckart s expression @xcite @xmath3 , where the ( non - negative ) quantity @xmath4 is the ( generally not constant ) bulk viscosity coefficient and @xmath5 is the fluid expansion scalar which in the homogeneous and isotropic background reduces to @xmath6 , where @xmath7 is the hubble parameter and @xmath8 is the scale factor of the robertson - walker metric . by this assumption we ignore all the problems inherent in eckart s approach which have been discussed and resolved within the israel - stewart theory @xcite ( see also @xcite and references therein ) . we expect that for the applications we have in mind here , the differences are of minor importance . it is obvious that the bulk viscosity contributes with a negative term to the total pressure and hence a dissipative fluid seems to be a potential dark energy candidate . however , a cautionary remark is necessary here . in traditional non - equilibrium thermodynamics the viscous pressure represents a ( small ) correction to the ( positive ) equilibrium pressure this is true both for the eckart and for the israel - stewart theories . here we shall admit the viscous pressure to be the dominating part of the pressure . this is clearly beyond the established range of validity of conventional non - equilibrium thermodynamics . as already mentioned , non - standard interactions are required to support such type of approach @xcite . of course , this reflects the circumstance that dark energy is anything but a  standard " fluid . ( we mention that viscosity has also been suggested to have its origin in string landscape @xcite ) . to successfully describe the transition to a phase of accelerated expansion , preceded by a phase of decelerated expansion in which structures can form , it is necessary that the viscous pressure is negligible at high redshifts but becomes dominant later on .    in ref . @xcite a one - component bulk viscous model ( bv model ) of the cosmic medium was investigated in which baryons were not taken into account . as far as the sne ia data are concerned , the results of this model were similar to those obtained for a generalized chaplygin gas ( gcg ) model . but while gcg models predict small scale instabilities and oscillations at the perturbative level , the corresponding matter power spectrum of the bv model turned out to be well behaved . of course , the results of a model that does not include baryons can not be seen as conclusive ; after all , the observed power spectrum describes the distribution of baryonic matter . in the present paper , we perform an advanced analysis which properly includes a separately conserved baryon component . thus we establish a more realistic viscous fluid scenario of the cosmic substratum .    for the quantitative calculations will use a constant coefficient of bulk viscosity . a qualitative analysis with the help of dynamical system methods is applied to visualize the space of more general cosmological background scenarios . the corresponding phase space reveals interesting new features compared with those generally found in similar models ( see , for example , @xcite and references therein ) . in the present approach there is an entire singular axis . solutions of the desired type are generated , i.e. , solutions for which an initial subluminal expansion is followed by superluminal expansion . these solutions are characterized by a bulk viscosity coefficient with a power law dependence on the energy density , @xmath9 , where @xmath10const and @xmath11 . we test the results of our model both against sne ia data and the observed matter power spectrum . both the results for type ia supernovae and for the matter power spectrum show that the bv model is competitive with the @xmath1cdm model as well as with quintessence and different chaplygin gas type models . in particular , our minimum @xmath12 value for the sne ia data is similar to the @xmath12 values of those models . furthermore , the matter power spectrum represents a good fit to the corresponding observational data . we argue that the absence of oscillations and instabilities is a consequence of the fact that the pressure perturbations in our model are intrinsically non - adiabatic . the paper is organized as follows . in section [ background ] we define our model and specify the homogeneous and isotropic background dynamics . section [ dynamical system ] performs a qualitative analysis using dynamical system methods . in section [ sn ] , the sne ia data are used to restrict the values of the physically relevant free parameters . in section [ matter power spectrum ] the matter power spectrum is determined and compared with large scale structure observations . in section [ conclusions ] we present our conclusions . a bulk viscous fluid is characterized by an energy density @xmath13 and a pressure @xmath14 which has a conventional component @xmath15 and a bulk viscosity component @xmath16 , such that @xmath17 on thermodynamical grounds the bulk viscosity coefficient @xmath18 is positive , assuring that the viscosity pushes the effective pressure towards negative values . in fact , the expression ( [ bulk ] ) is the original proposition for a relativistic dissipative process @xcite . as already mentioned , it follows from the more general israel - stewart theory @xcite in the limit of a vanishing relaxation time . we shall assume this approximation to be valid throughout the paper . let us consider the cosmic medium to consist of a viscous fluid of the type ( [ bulk ] ) , which is supposed to characterize the dark sector , and of a pressureless fluid that describes the baryon component . hence , the relevant set of equations is @xmath19 the ( super)subscripts @xmath20 and @xmath21 indicate the viscous and the ( baryonic ) matter components , respectively . since the matter is pressureless , the total pressure @xmath14 of the cosmic medium coincides with the pressure @xmath22 of the viscous fluid . for the homogeneous and isotropic background dynamics we shall restrict ourselves to the flat friedmann - lematre - robertson - walker ( frlw ) metric , @xmath23 favored by the cmb anisotropy spectrum @xcite . the dynamic equations then are : @xmath24 where the bulk viscosity coefficient was assumed to have a power law dependence on the energy density @xmath25 according to @xmath26 with @xmath27 const . the dot denotes differentiation with respect to the cosmic time . of course , not all these equations are independent . equation ( [ em3 ] ) is decoupled and leads to @xmath28 here and in the following , quantities with a subscript 0 refer to the present epoch and we have used @xmath29 . from now on we will set @xmath30 , thus assuming the dissipative pressure to be the dominating contribution .    with a change of variables , @xmath31 where the prime means derivative with respect to the scale factor @xmath8 , and using relation ( [ em1 ] ) , the conservation equation for the viscous component becomes @xmath32 this equation can be recast in an integral form : @xmath33 we will be mainly interested in the case @xmath34 . there is a simple , direct relation between @xmath25 and the scale factor @xmath8 for @xmath35 . in this particular case , we find , @xmath36^{2}-1\biggl\}\quad , \qquad \quad\left(\nu = 0\right)\ , \label{viscous - dens}\ ] ] where @xmath37 is an integration constant . it can easily be verified that @xmath38 when @xmath39 ( pressureless matter ) and @xmath40 cte ( cosmological constant ) when @xmath41 . it is expedient to notice that for @xmath42 the total energy density @xmath43 coincides with the energy density of a specific gcg @xcite . generally , a gcg is characterized by an equation of state ( @xmath44 const @xmath45 ) @xmath46 which corresponds to an energy density @xmath47^{{1}/{\left(1+\alpha\right)}}\ , \label{rhogcg}\ ] ] where @xmath48 is another ( non - negative ) constant . with the identifications @xmath49 and @xmath50 it becomes obvious that for @xmath51 the energy densities @xmath43 from ( [ viscous - dens ] ) and @xmath52 in ( [ rhogcg ] ) coincide . the analogy between dissipative and gcg models was also pointed out in @xcite . the structure of the total energy density @xmath53 allows us to perform a decomposition of the viscous - baryon system into three non - interacting components plus the baryonic fluid : @xmath54 with @xmath55 and @xmath56 the individual equations of state of the different components are : @xmath57 consequently , our system is equivalent to a mixture of two dark energy type fluids , one dark matter fluid and the baryon component . while the bv model is a unified description of dark matter / energy , the option of a decomposition may be useful in comparing the scenario with some observational data . the x - ray measurement of galactic clusters @xcite is an example which requires a separation of a matter component of the cosmic medium . for a chaplygin gas , different decompositions have been proposed in the literature @xcite . moreover , the decomposition ( [ decomp ] ) with ( [ split ] ) - ( [ decompp ] ) reveals , that the background dynamics of our model differs from that of the @xmath1cdm model ( including a baryon component ) by the existence of the second dark energy component with energy density @xmath58 . the relevant cosmological quantities are the present density parameters @xmath59 and @xmath60 which represent the fractions of the viscous fluid and of the pressureless matter , respectively , with respect to the total density ( remember that we restrict ourselves to the flat case ) , the hubble parameter today @xmath61 , the present value of the deceleration parameter @xmath62 , and the age of the universe @xmath63 . the apparently new parameter of the present model , the viscosity coefficient , can be expressed in terms of the other quantities . in fact , the relation ( [ viscous - dens ] ) can be rewritten in a more convenient way . using friedmann s equation today , we have @xmath64 combining ( [ viscous - dens ] ) with ( [ def1 ] ) and ( [ def2 ] ) we obtain the relation @xmath65 now , if we use equation ( [ em5 ] ) with @xmath66 and the definition of the deceleration parameter @xmath67 , we find @xmath68 from the definition of @xmath69 in equation ( [ defk ] ) we have @xmath70 using equation ( [ em5 ] ) , evaluated today , it follows that @xmath71 which leads to @xmath72 hence , the viscous fluid energy density becomes , @xmath73\quad . \label{def3}\ ] ] all these relations will be useful for defining the observables that will be constrained later on by the sne ia and matter power spectrum data . the relations of the previous section that will be used to compare the theoretical predictions of our model with observational data , are valid under the restriction @xmath74 . it is important to have an idea of how serious this restriction is . if the case @xmath74 is very particular , our results will not be conclusive . therefore it is desirable to have at least a qualitative analysis of the general case ( [ em1])-([em4 ] ) which , because of its complexity , can not be solved analytically . such an analysis could reveal whether or not the case @xmath74 is , in a sense , typical . to this purpose we shall use dynamical system techniques for a system of non - linear differential equations @xcite . the system of equation ( [ em1])-([em4 ] ) can be recast in a more convenient form by using geometric unities @xmath75 and fixing the scale so that @xmath76 . defining @xmath77 and @xmath78 , the following system of equations is obtained : @xmath79 since the density @xmath80 of the baryonic matter must be non - negative , the solutions that have physical meaning are those which satisfy @xmath81 . ( upper pane ) and for @xmath82 ( lower panel ) . for @xmath83 the directions of the arrows in the lower panel have to be reversed . ]     ( upper pane ) and for @xmath82 ( lower panel ) . for @xmath83 the directions of the arrows in the lower panel have to be reversed . ]    in a first step we have to identify the critical points @xmath84 for which @xmath85 . these points are either attractors or repellers or saddle points . the nature of a critical point is determined by the eigenvalues of the matrix    @xmath86 ,    where @xmath87 , @xmath88 and @xmath89 . in general , there are two eigenvalues @xmath90 for the characteristic equation of this matrix . if both eigenvalues are positive , than the critical point is repeller ; if both are negative , the critical point is an attractor ; if one is positive and the other one is negative , than the critical point is a saddle . with the help of a suitable transformation @xcite , the critical points at infinity can be studied . a projection onto the poincar sphere @xcite will then allow us to represent the entire phase space in a finite region in the @xmath91 plane . it turns out that the entire @xmath92 axis is a stationary solution for the system , representing a static space - time . a singular point , denoted by @xmath93 , is @xmath94 where @xmath34 . the characteristic equation for this critical point is given by @xmath95 where @xmath96 . the eigenvalues are @xmath97 and @xmath98 . the phase diagram depends on the value of the parameter @xmath99 . if @xmath11 , the critical point @xmath100 is an attractor : all solutions converge to it . this point represents the de sitter phase . hence , all solutions approach a de sitter phase asymptotically . if @xmath101 , the point @xmath100 becomes a saddle point and all solutions approach a minkowski space - time in the origin . if @xmath102 , @xmath100 is still a saddle point , but the directions of the curves change : all solutions leave the minkowski space - time . these different behaviors are shown in figure ( [ fig1 ] ) . for the point @xmath100 we have @xmath103 , which corresponds to a vanishing baryon density . ( recall that the range @xmath104 is unphysical . ) the deceleration parameter @xmath105 can be written as @xmath106 , equivalent to @xmath107 . combination with eq . ( [ dotx ] ) shows that there is accelerated expansion @xmath108 for @xmath109 . the special case @xmath74 has @xmath110 and accelerated expansion for @xmath111 . this solution is of the type of the curves that approach the point @xmath100 from right below in the left panel of fig . notice that there is a projection effect due to the circumstance that infinity is mapped onto one single point on the @xmath112 axis . the curve that points from the origin towards @xmath100 in the left part of fig . [ fig1 ] corresponds to a phantom scenario with @xmath113 which separates the physical and unphysical regions . we are not interested here in this type of models . the type ia supernovae allow us to obtain informations about the universe at high redshifts . today , observers have detected about 300 of these objects with redshifts up to almost @xmath114 . in what follows , we will use the more restricted sample of 182 sne ia of the gold06 data set which consists of objects for which the observational data are of very high quality @xcite . the relevant quantity for our analysis is the moduli distance @xmath115 , which is obtained from the luminosity distance @xmath116 by the relation @xmath117 .the best - fitting parameters , i.e. , when @xmath118 is minimum , for the bv model , the gcg model , the traditional chaplygin gas model ( @xmath119 ) and the @xmath120cdm model ( @xmath121 ) for a flat universe ( @xmath122 ) . @xmath61 is given in @xmath123 , @xmath124 in units of @xmath125 ( @xmath126 - speed of light ) and @xmath63 in @xmath127 . [ cols="^,^,^,^,^ " , ]     in order to compare the theoretical results with the observational data , the first step is to compute the quality of the fitting through the least square fitting quantity @xmath12 . we adopt the hst ( hubble space telescope ) prior @xcite for @xmath61 , as well as the cosmic nucleosynthesis prior for the baryonic density parameter @xmath128 ( where @xmath129 is @xmath61 divided by @xmath130 km / s / mpc ) , which leads to @xmath131 the quantities @xmath132 are the distance moduli , observationally measured for each supernova of the @xmath133 gold06 sne ia dataset @xcite and the @xmath134 are the corresponding theoretical values . the @xmath135 represent the measurement errors and include the dispersion in the distance moduli , resulting from the galaxy redshift dispersion which is due to the peculiar velocities of the objects ( cf . @xcite ) . the smaller the @xmath12 , the better the agreement between the theoretical model and the observational data . table [ tablebestfitsneia ] shows that the @xmath12 value for bv model is competitive with the gcg model ( cf . ( [ rhogcg ] ) ) , the traditional chaplygin gas model ( cf . ( [ rhogcg ] ) with @xmath119 ) and the @xmath120cdm model . compared with the @xmath120cdm and the ( generalized ) chaplygin gas models , the age of the universe is greater for the bulk viscous model and its deceleration parameter is less negative , i.e. , the universe is less accelerated . using bayesian statistics , a probability distribution can be constructed from the @xmath136 parameter @xcite . in the present case , there are three free parameters : the hubble parameter today @xmath137 , the pressureless baryonic density parameter @xmath138 and the auxiliary quantity @xmath37 which is connected to the present value of the deceleration parameter ( cf . ( [ bq0 ] ) ) . marginalization over one or two of these parameters will lead to corresponding two- or one - dimensional representations . the details of the bayesian statistics and computational analysis , relying on betocs ( * * b**ay**e**sian * * t**ools for * * o**bservational * * c**osmology using * * s**ne ia ) , are given in ref . @xcite . the main feature of the parameter estimation for the viscous fluid model is the gaussian - like distribution around the preferred value .    , the baryon density and the hubble constant . the joint pdf peak is shown by the large dot , the confidence regions of @xmath139 ( @xmath140 ) by the red dotted line , the @xmath141 ( @xmath142 ) in blue dashed line and the @xmath143 ( @xmath144 ) in green dashed - dotted line . ]    , the baryon density and the hubble constant . the joint pdf peak is shown by the large dot , the confidence regions of @xmath139 ( @xmath140 ) by the red dotted line , the @xmath141 ( @xmath142 ) in blue dashed line and the @xmath143 ( @xmath144 ) in green dashed - dotted line . ] this can easily be seen both in the table [ tableparestsneia ] ( almost symmetric error bars ) and in the two and one - dimensional diagrams of figures [ figs2dsneia ] and [ figs1dsneia ] . in figure [ figs2dsneia ] we show the two - dimensional probability distribution for the parameters @xmath37 , @xmath60 and @xmath61 at one , two and three sigma , corresponding , respectively , to @xmath145 , @xmath146 and @xmath147 confidence levels . the corresponding one - dimensional probability distributions are shown in figure [ figs1dsneia ] . the main features of the model are : the age of the universe is considerably larger than for the @xmath1cdm and gcg models , around @xmath148 gy , with a small dispersion ; the baryon density is around @xmath149 , a little larger than predicted by nucleosynthesis , but there is a marginal agreement ; the hubble constant is around @xmath150 km / s / mpc , i.e. much smaller than the one predicted by the cmb , which is around @xmath151 km / s / mpc , but in agreement with the @xmath1cdm and gcg models when only the data from sne ia are used ; the deceleration parameter is bigger ( smaller absolute value ) than in those models . for the parameter estimations in the @xmath1cdm and in the gcg models , see ref . @xcite and references therein .    , the baryonic density , the hubble constant and the deceleration parameter . the joint pdf peak is shown by the large dot , the confidence regions of @xmath139 ( @xmath140 ) by the red dotted line , the @xmath141 ( @xmath142 ) in blue dashed line and the @xmath143 ( @xmath144 ) in green dashed - dotted line . ]    , the baryonic density , the hubble constant and the deceleration parameter . the joint pdf peak is shown by the large dot , the confidence regions of @xmath139 ( @xmath140 ) by the red dotted line , the @xmath141 ( @xmath142 ) in blue dashed line and the @xmath143 ( @xmath144 ) in green dashed - dotted line . ]    , the baryonic density , the hubble constant and the deceleration parameter . the joint pdf peak is shown by the large dot , the confidence regions of @xmath139 ( @xmath140 ) by the red dotted line , the @xmath141 ( @xmath142 ) in blue dashed line and the @xmath143 ( @xmath144 ) in green dashed - dotted line . ]    , the baryonic density , the hubble constant and the deceleration parameter . the joint pdf peak is shown by the large dot , the confidence regions of @xmath139 ( @xmath140 ) by the red dotted line , the @xmath141 ( @xmath142 ) in blue dashed line and the @xmath143 ( @xmath144 ) in green dashed - dotted line . ] the gaussian nature of the distribution is responsible for the fact that the best fitting set of values is very close to the marginalized parameter estimations , with a small dispersion . generally , this is neither true for the @xmath1cdm model nor for the gcg model . again : the most remarkable differences between the viscous model and the @xmath1cdm model are the larger age of the universe and the smaller absolute value of the deceleration parameter . in order to write the first order perturbed equations , we re - express the field equations as    @xmath152    we introduce the quantities , @xmath153 where @xmath154 , @xmath13 and @xmath155 are the known background solutions for the metric , the energy density and the four velocity , respectively , and @xmath156 , @xmath157 and @xmath158 are the corresponding perturbations . we will perform the calculations using the synchronous coordinate condition , @xmath159 since we will be interested mainly in perturbations that are inside the horizon , the choice of the gauge is not really essential . the perturbed ricci tensor takes the form , @xmath160 where @xmath161 \quad , \ ] ] are the perturbations of the christoffel symbols . the relevant non - vanishing components are : @xmath162 for the components of the perturbed energy - momentum tensor we have @xmath163 + using the definitions , @xmath164 and performing a plane wave decomposition of all perturbed functions according to @xmath165 where @xmath166 is the wave vector , we find , after a long but straightforward calculation , the following set of first order perturbed equations : @xmath167 replacing @xmath168 by using the last relation , we obtain @xmath169 with the help of the constraint ( [ last ] ) , this system can be further reduced to a set of three coupled differential equations : @xmath170 + now it is convenient to introduce the background equations and the quantities in eq . ( [ def3 ] ) from which we find @xmath171 \quad , \label{background2}\\ \frac{12\pi g\xi_0}{h_0 } & = & \frac{1}{2}(1 - 2q_0 ) \quad . \label{background3}\end{aligned}\ ] ] the following relations will be useful later on : @xmath172 + in a next step we divide the set of equations ( [ 3system1 ] ) - ( [ 3system3 ] ) by @xmath173 and redefine @xmath174 . again , after a fairly long calculation , the perturbed equations reduce to the following system of coupled linear differential equations : @xmath175\biggl\}\dot\delta_b & = &   \nonumber \\ ( 1 + g)\delta\dot p_v + \biggr\{2\frac{\dot a}{a}(1 + g ) + \frac{1 - 2q_0}{2}\biggr[\frac{f^2}{f^2 - 1 } + g\biggl ] & &   \nonumber \\ + \biggr(k\ , l_h\biggl)^2\frac{(1 - 2 q_0)a}{9\omega_{b0}(f^2 - 1)}\biggl\}\delta p_v & & \quad , \end{aligned}\ ] ] where @xmath176 . now , we re - express the perturbed equations in terms of the variable @xmath8 according to ( [ prime ] ) , which implies @xmath177 we also need the following relations , which result from the background equations : @xmath178 with the redefinition @xmath179 , the final form for the set of perturbed equations is : @xmath180      evidently , equations ( [ f - pe1])-([f - pe3 ] ) are too complicated to admit analytical solutions . hence , we proceed integrating these equations numerically . in order to do so , an important problem is to fix the initial conditions . a scale invariant primordial spectrum ( as predicted by the inflationary scenario ) is assumed . but , since we are interested in the power spectrum today , we have to follow the evolution of this spectrum until the present phase . the corresponding procedure for the @xmath1cdm model has been carried out in terms of the bbks transfer function @xcite . this function is characterized by @xmath181 where @xmath182 and @xmath183 are the density parameters of dark matter and baryonic matter of the @xmath1cdm model . the baryonic transfer function is approximated by by the numerical fit @xcite : @xmath184^{-1/4}\,.\end{aligned}\ ] ] the quantity of interest is the power spectrum for the baryonic matter today . this spectrum is given by @xmath185 where @xmath186 is the total density parameter of the @xmath1cdm model . the function @xmath187 in ( [ powers ] ) is the growth function @xmath188^{-1}\,\end{aligned}\ ] ] in which @xmath189 represents the fraction of the total energy , contributed by the cosmological constant . the normalization coefficient @xmath190 in ( [ powers ] ) can be fixed by using the cobe measurements of the cmb anisotropy spectrum . this coefficient is connected to the quadrupole momentum @xcite @xmath191 of this spectrum by the relation @xmath192 where @xmath193 is the present cmb temperature . the quadrupole anisotropy is taken as @xmath194 this value is obtained from the cobe normalization and is consistent with the more recent results of the wmap measurements which use the prior of a scale invariant spectrum @xcite . taking all these estimates into account , one may fix @xmath195 which corresponds to the initial vacuum state for the density perturbations used in the bbks transfer function . we recall that all the relations ( [ tq])-([a ] ) are valid for the @xmath1cdm model . strictly speaking , we would have to redo all the calculations to obtain the analogue of the expression ( [ powers ] ) for our bulk viscous model . however , since we are interested in the shape of the spectrum , there exists a simpler way which relies on the circumstance that the growth function @xmath187 in ( [ cpt ] ) is a pure background quantity . use of the @xmath1cdm growth function and the @xmath1cdm initial values also for our model will result in an overall shift of the spectrum @xmath196 without affecting the structure of the latter . the resulting data will then be a fit to a  wrong " , i.e. unphysical @xmath1cdm model , in which the composition of the dark sector differs from the standard one with roughly @xmath197 of dark matter and @xmath198 of a cosmological constant . however , this procedure provides us with the correct structure of the spectrum for our bulk viscous model . keep in mind that the background dynamics of our model differs from that of the @xmath1cdm model by the existence of a second dark energy component ( component 2 ) in eq . ( [ split ] ) . well inside the matter dominated phase ( say , @xmath199 ) , we use the @xmath1cdm initial conditions , insert them in equations ( [ f - pe1])-([f - pe3 ] ) , and let this system evolve . then we compute the power spectrum @xmath200 at @xmath201 , equivalent to @xmath202 . finally , we compare the theoretical results with the power spectrum obtained by the 2dfgrs observational program . the spectrum depends essentially on two parameters : @xmath138 and @xmath62 . from now on , relying on the results of the type ia sne analysis and the primordial nucleosynthesis constraint , we fix @xmath203 ( according to the type ia sne analysis of the previous section , this is marginally satisfied at least at @xmath204 ) .    in figures [ figpower1 ] and [ figpower2 ] we present the best fitting for @xmath205 . to have a good agreement , the ratio between the quantity of dark matter / dark energy in the  wrong " @xmath1cdm reference model has to be in the range from @xmath206 for @xmath207 , to about @xmath208 for @xmath209 . it is seen that the baryonic matter power spectrum of the viscous fluid model is in good agreement with the observational data for values of @xmath62 that are close to those predicted by the type ia sne analysis , that is , @xmath210 . we remark _ en passant _ that a variation of the value of @xmath60 within the interval predicted by the type ia sne at @xmath204 does not change substantially the results presented in figures [ figpower1 ] and [ figpower2 ] . in any case , there is no instability in the computed power spectrum . this is true not only for the baryonic fluid power spectrum but also for dark viscous fluid power spectrum . we note , however , that as @xmath62 approaches zero , there is a power depression at large scales ( small @xmath69 ) compared with the @xmath1cdm reference model . ( left ) and @xmath211 ( right ) . the initial conditions using the bbks transfer function correspond to  wrong " ratios ( see text ) of dark matter and dark energy : @xmath212 and @xmath213 ( left ) and @xmath214 and @xmath215 ( right ) . ]     ( left ) and @xmath211 ( right ) . the initial conditions using the bbks transfer function correspond to  wrong " ratios ( see text ) of dark matter and dark energy : @xmath212 and @xmath213 ( left ) and @xmath214 and @xmath215 ( right ) . ]     ( left ) and @xmath216 ( right ) . the initial conditions using the bbks transfer function correspond to  wrong " ratios of dark matter and dark energy : @xmath217 and @xmath218 ( left ) and @xmath219 and @xmath215 ( right ) . ]     ( left ) and @xmath216 ( right ) . the initial conditions using the bbks transfer function correspond to  wrong " ratios of dark matter and dark energy : @xmath217 and @xmath218 ( left ) and @xmath219 and @xmath215 ( right ) . ]      in the homogeneous and isotropic background , the bv model is equivalent to a gcg with @xmath220 . the qualitative differences between both models at the perturbative level can be traced back to a difference in the pressures that characterize the cosmic medium in each of these cases . in the background both pressures coincide . but while the perturbations for the gcg are entirely adiabatic , this is not the case for our present model .    in the background we have @xmath221 for the viscous fluid and @xmath222 for the gcg with @xmath223 ( cf . ( [ pgcg ] ) ) . use of friedmann s equation @xmath224 ( for the spatially flat case ) reveals that in both cases we have an equation of state @xmath225 , where the constants are related by @xmath226 . the first order pressure perturbations for the gcg are @xmath227 where @xmath228 is the adiabatic sound speed square . for @xmath229 the sound speed square is positive , giving rise to a ( non - observed ) oscillatory behavior . for @xmath230 ( our case ) this square is negative which results in instabilities which are not observed either .    for the viscous fluid the pressure perturbations are qualitatively different . in particular , the propagation of perturbations is no longer given by the adiabatic sound speed . the perturbations are non - adiabatic . the non - adiabaticity is characterized by @xmath231 where @xmath232 . in general , the right hand side of ( [ nad ] ) does not vanish ( for a similar situation see @xcite ) . it is this non - adiabatic character of the perturbations which makes the power spectrum well behaved on small scales in contrast to the gcg case . a more detailed analysis of this feature will be given elsewhere . we presented a two - component cosmological model in which one component represents the dark sector , the other one pressureless baryons . the dark sector is described by a single bulk viscous fluid with a constant bulk viscosity coefficient @xmath233constant . a dynamical system analysis was used to show that the assumption of a constant @xmath4 does not seem to be too restrictive and that this model embodies quite general features and does not represent a very particular configuration . in the homogeneous and isotropic background the total energy density ( i.e. , including the baryons ) is equivalent to that of a gcg with @xmath234 . but while the perturbation dynamics of a gcg is entirely adiabatic , the present model exhibits non - adiabatic features . as a consequence , it can avoid the problems that usually plague such type of unified dark matter / dark energy models , namely the appearance of ( non - observed ) small scale oscillations or instabilities . the predictions of our model were compared with sne ia data and with the 2dfgrs results . for the comparison with sne ia data we employed the so - called gold sample of 182 good quality high redshift supernovae . our bayesian statistics analysis leads to the following parameter evaluation at @xmath235 level : @xmath236 , @xmath237 , @xmath238 km / s / mpc and @xmath239 gy . the main differences in comparison with a similar analysis for the @xmath1cdm and the gcg models @xcite are : the age of the universe is considerably larger ; the deceleration parameter is larger ( smaller in absolute value ) ; the baryon density is only marginally compatible with the primordial nucleosynthesis result . a remarkable feature is the ( almost ) perfect gaussian distribution for all parameters with a quite small dispersion .    to compare the results of our perturbation analysis with the observed matter power spectrum we have used the bbks transfer function to fix the initial condition at a very high redshift ( @xmath199 ) , deep in the matter dominated phase . since the bbks transfer function has been conceived for a two - component description of the dark sector , dark matter and a cosmological constant , this must be seen with some caution . we have argued that , in spite of this problem , the resulting power spectrum preserves its shape and can be compared with observations . a good fitting can be obtained for values of @xmath62 which are in agreement with the type ia sne analysis . the best results are found for @xmath240 . qualitatively , our analysis confirms previous results for a simplified model of the cosmic medium , in which the baryon component was neglected @xcite . r. maartens 1997 causal thermodynamics in relativity in _ proceedings of the hanno rund conference on relativity and thermodynamics _ ed s d maharaj , university of natal , durban pp 10 - 44 . ( astro - ph/9609119 ) . r. colistete jr , j. c. fabris , s.v.b . gonalves and p.e . de souza , int . d13 * , 669 ( 2004 ) ; r. colistete jr . , j. c. fabris and s.v.b . gonalves , int . j. mod d14 * , 775 ( 2005 ) ; r. colistete jr . and j. c. fabris , class . * 22 * , 2813 ( 2005 ) ; r. colistete jr . and r. giostri , _ betocs using the 157 gold sne ia data : hubble is not humble _ , arxiv : astro - ph/0610916 .
Abstract:

Answer: we propose a scenario in which the dark components of the universe are manifestations of a single bulk viscous fluid . using dynamical system methods , a qualitative study of the homogeneous , isotropic background scenario is performed in order to determine the phase space of all possible solutions . the specific model which we investigate shares similarities with a generalized chaplygin gas in the background but is characterized by non - adiabatic pressure perturbations . this model is tested against supernova type ia and matter power spectrum data . different from other unified descriptions of dark matter and dark energy , the matter power spectrum is well behaved , i.e. , there are no instabilities or oscillations on small perturbation scales . the model is competitive in comparison with the currently most popular proposals for the description of the cosmological dark sector .