one of the most important problems today in cosmology is the nature of dark matter and dark energy that must dominate the matter content of the universe @xcite .
the existence of dark matter is suggested by the anomalies in the dynamics of galaxies and clusters of galaxies @xcite .
dark energy seems to be an inevitable consequence of the present acceleration of the universe as indicated by the type ia supernovae ( sne ia ) data , which asks for a fluid of negative pressure which does not agglomerate at small scales @xcite . while in general the nature of dark matter
can be connected with relics of fundamental theories , like axions @xcite , there are two main candidates to represent dark energy : cosmological constant @xcite and quintessence which is a self interacting scalar field @xcite@xcite .
both models suffers of many drawbacks , mainly linked with fine tuning of parameters , as in the quintessence model @xcite , or disagreement with the theoretically predicted values , like in the cosmological constant case @xcite .
the traditional chaplygin gas model ( cgm ) and the  generalized chaplygin gas model ( gcgm ) have been widely considered as alternatives to the cosmological constant and to quintessence as the dark energy that drives the present acceleration of the universe @xcite@xcite .
the appealing of the gcgm comes from , among other reasons , the fact that it can unify the description of dark matter and dark energy .
the generalized chaplygin gas ( gcg ) is characterized by the equation of state @xmath14 where @xmath15 and @xmath16 are constants that in general are considered as positive numbers .
the conservation equation leads to the relation @xmath17 some of the main interesting features of eq .
( [ densidade ] ) are the following . initially , the gcgm behaves as a dust fluid , with @xmath18 , while at late times the gcgm behaves as a cosmological constant term , @xmath19 .
hence , the gcgm interpolates a matter dominated phase ( where the formation of structure can happen ) and a de sitter phase . at the same time
, the sound velocity of the gcgm is positive which assures that , in spite of exhibiting a negative pressure , it is stable against perturbations of the background configuration @xcite .    a lot of effort has been done in order to constrain the free parameters of the gcgm .
type ia supernovae data @xcite@xcite , the spectra of anisotropy of cosmic microwave background radiation @xcite , the mass power spectrum @xcite , gravitational lenses @xcite , x - ray data @xcite and age estimates of high-@xmath20 objects @xcite have extensively been used in this sense .
this is essential for many obvious reasons .
one of them is that if the parameter @xmath16 is close to zero , the gcgm reduces essentially to the cosmological constant model .
if observations indicate that this is the preferred value , the appealing for the cosmological constant as dark energy is reinforced . in a previous work @xcite
, we have analysed the case of the traditional chaplygin gas model , where @xmath21 , using the sne ia data .
one important aspect of this analysis is the use of the bayesian statistics , for one side , and also the fact that all free parameters ( hubble constant , the curvature term , the value of the sound velocity and the proportion of ordinary matter and chaplygin gas ) of the model were taken into account in the treatment of the problem .
this contrasts with most of the analyses previously found in the literature in the sense that they employ a more simplified statistical analysis or some input parameters were fixed a priori @xcite@xcite . in that work @xcite , it was concluded that the unification scenario ( i.e. , without pressureless matter ) , and a sound velocity near ( but not equal ) to the velocity of light were preferred .
moreover , a closed universe was clearly favoured .
the aim of the present work is to apply the same analysis , using the bayesian statistics and taking into account all free parameters , to the gcgm . in order to keep contact with the previous work
, we will use the same supernovae sample .
this sample is much smaller than those available today ( see , for example , ref .
@xcite@xcite ) .
but , the quality of the data is extremely good , which leads to a smaller value of @xmath22 ( the quantity that measures the quality of the fitting ) : with this restricted sample the best - fitting gives @xmath22 around @xmath23 , while for the complete `` gold '' sample @xcite , the value mounts to unity . in the gcgm
there are five free parameters instead of the four free parameters of the traditional chaplygin gas model ( @xmath21 ) : the hubble constant @xmath24 , the parameter @xmath0 related to the sound velocity of the chaplygin gas , the curvature density parameter @xmath25 , the ordinary matter parameter @xmath26 ( or alternatively , the chaplygin gas parameter @xmath27 ) and the equation of state parameter @xmath16 .
one important point in performing this generalized analysis is the range of validity of the parameter @xmath16 .
usually , it has been assumed that @xmath28 .
the reason is twofold : @xmath16 greater than one could lead , even if in the future , to a sound velocity greater than the velocity of the light ; @xmath16 negative should lead to an imaginary sound velocity , that is , instability .
however , an analysis of the supernovae data may taken into account values of @xmath16 greater than one or negative .
in particular , in reference @xcite , using a different statistical approach with respect to the method employed here , a value of @xmath29 near @xmath30 was favoured .
in fact , values of @xmath16 outside that limited range can be considered .
the reason lies in the fact that the gcgm with the equation of state ( [ eos ] ) may be a phenomenological manifestation of a theory expressed in terms , for example , of scalar field with a potential , or even of tachyon condensate model @xcite . in this case
, the notion of sound velocity can not be identified with the ordinary expression used in the perfect fluid approach , and no violation of causality or instability may occur if @xmath16 is greater than one or negative .
hence , in our analysis the parameter @xmath16 will vary from negative to highly positive values .
the prediction for each parameter ( @xmath16 , @xmath24 , @xmath25 , @xmath26 , @xmath27 and @xmath0 ) is taken marginalizing , i.e. , integrating over the other parameters .
we stress that such marginalization procedure gives predictions that can be quite different from a two dimensional analysis , for example , where one searches for the simultaneous preferred values for any two parameters .
this marginalization is in fact essential in order to have the correct prediction for a given free parameter .
we will consider the following cases : all five parameters free ; spatial curvature zero ( four free parameters ) ; pressureless matter given only by the barionic component ( four free parameters ) ; pressureless matter absent ( four free parameters ) ; spatial curvature zero and no pressureless matter component ( three free parameters ) ; spatial curvature zero and pressureless matter component given by the barionic component ( three free parameteres ) .
for the more general case with five free parameters the predictions we obtain are the following : @xmath31 , @xmath32 , @xmath5 , @xmath33 , @xmath34 , @xmath35 .
moreover the predicted age of the universe and the value of de deceleration parameter today are @xmath36 and @xmath37 , respectively .
the results are consistent with , for example , the age of globular clusters @xcite .
two crucial parameters are @xmath16 and @xmath38 , since if @xmath39 and/or @xmath40 are favoured , the gcgm becomes equivalent to the cosmological constant model ( @xmath41cdm ) .
the value of the parameter @xmath0 is in fact peaked at unity .
but for reasons discussed later , this does not imply that this is the most favoured value ; indeed it is just expected to be near unity . on the other hand , the parameter @xmath16 is peaked in a negative value , but the spread is very large . in particular , the traditional chaplygin gas with @xmath42 can not be excluded or even disfavoured : in most of the cases this value of @xmath16 is favoured with a probability around @xmath43 .
hence , our analysis does not exclude the traditional chaplygin gas model .
this paper is organized as follows . in next section ,
the basic equations and quantities are exhibited .
section @xmath44 explains the type ia supernovae fitting .
in section @xmath45 , we display the results of the bayesian analysis of type ia supernovae data with the gcgm .
we present our conclusions in section @xmath46 .
regardless the possibility of unification of dark matter and dark energy using the gcgm , we will consider a model for the universe today with two fluids : pressureless matter and the gcg .
in doing so , our aim is to verify whether the unified model is the case favoured by observations .
moreover , even if in the gcgm there is no need of a dark matter component , there is baryonic matter in the universe which is represented also by a pressureless matter .
hence , the dynamics of the universe is driven by the friedmann s equation @xmath47 where @xmath48 and @xmath49 stand for the pressureless matter and generalized chaplygin gas component .
dot means derivative with respect to the cosmic time @xmath50 . as usual
, @xmath51 indicates a flat , closed and open spatial section . in the conservation equation for the generalized chaplygin gas component
( [ be5 ] ) , the equation of state ( [ eos ] ) has been introduced .
the equations expressing the conservation law for each fluid ( [ be4])-([be5 ] ) lead to @xmath52 henceforth , we will parametrize the scale factor such that its value today is equal to unity , @xmath53 . hence
, @xmath54 and @xmath55 are the pressureless matter and gcg densities today . eliminating from the last relation the parameter @xmath56 , the gcg density at any time can be reexpressed as @xmath57 where @xmath58 .
this parameter @xmath0 is connected with the sound velocity for the gcg today by the relation @xmath59 .
the luminosity distance is given by the expression @xmath60 @xmath61 being the co - moving coordinate of the source . using the expression for the propagation of light @xmath62 and the friedmann s equation ( [ be3 ] ) , we can reexpress the luminosity distance as @xmath63\quad , \ ] ] where @xmath64 and @xmath65^{1/(1+\alpha ) } -\omega _ { k0}(z^{\prime } + 1)^{2}\}^{1/2}}\quad , \ ] ] with the definitions @xmath66 where @xmath67 .
the final equations have been also expressed in terms of the redshift @xmath68 .
the age of the universe today can be expressed as @xmath69^{1/(1+\alpha ) } -\omega _ { k0}(z^{\prime } + 1)^{2}\}^{1/2}}\quad , \ ] ] where @xmath70 , so @xmath71 has units of years .
the value of the decelerated parameter today , @xmath72 , is given by @xmath73    the value @xmath74  of the scale factor that shows the beginning of the recent accelerating phase of the universe is given by the roots of the equation @xmath75 and is related to the redshift value @xmath76 such that @xmath77 or @xmath78 .
we proceed by fitting the sne ia data using the gcg model described above .
essentially , we compute the quantity distance moduli , @xmath79 and compare the same distance moduli as obtained from observations .
the quality of the fitting is characterized by the @xmath80 parameter of the least - squares statistic , as defined in ref .
@xcite , @xmath81 in this expression , @xmath82 is the measured value , @xmath83 is the value calculated through the model described above , @xmath84 is the measurement error , @xmath85 is the dispersion in the distance modulus due to the dispersion in galaxy redshift due to peculiar velocities .
this quantity we will taken as @xmath86 where , following ref .
@xcite , @xmath87 .    in this article , we have used the same @xmath88 type ia supernovae in table @xmath2 of ref .
@xcite , so we do not need to analyse and  compare again the cgm with the @xmath41cdm ( this work has been done in ref .
@xcite ) , but instead we can concentrate our efforts on the gcgm . although nowadays there are approximately @xmath89 sne ia available , the choice of these @xmath88 are confirmed to be of excellent quality by analysing the low values of @xmath90 ( the estimated errors for degree of freedom , i.e. , @xmath80 divided by @xmath88 , the number of sne ia ) , shown in table [ tablebestfitgcg ] .
@xmath91 values range between @xmath92 to @xmath93 , while other larger sne ia samples have been used in many cosmological models @xcite , typically resulting @xmath91 near the unity .
.the best - fitting parameters , i.e. , when @xmath94 is minimum , for each type of spatial section and matter content of the generalized chaplygin gas model .
@xmath24 is given in @xmath95 , @xmath96 in units of @xmath97 , @xmath71 in @xmath98 and @xmath74 in units of @xmath99 .
@xmath100 and @xmath101 are the constants of the relation @xmath102 . [
cols="^,^,^,^,^,^,^ " , ]     table [ tableparestgcg ] is very comprehensive and resumes our results , listing all the bayesian parameter estimations for six different cases of spatial section and matter content within the context of generalized chaplygin gas model .
the following sub - sections will analyse in detail this table and accompanying figures .
more emphasis is spent into the analyses of the generalized chaplygin gas parameter @xmath16 , because it is the subject of constant debate .    in figures
[ figsalpha ] , [ figsa ] and [ figsomegas ] , the one - dimensional pdf shapes are matched against half - gaussians to estimate the probabilities outside the left and right edges , in this way increasing the precision of the parameter estimations without excessively enlarging the parameter space .
the parameter estimation of @xmath16 is strongly related to the @xmath0 parameter , i.e. , the joint pdfs for the parameter space of @xmath104 do not have a gaussian shape , so the marginalization process changes the peak values and credible regions .
table [ tableparestgcg ] , figures [ figsalphaa ] and [ figsalpha ] are carefully analysed below with respect to @xmath16 . in brief , the estimations for @xmath16 are quite spread ( i.e. , they have large credible regions ) and they differ a lot depending on if two or one - dimensional parameter space is used .     for the generalized chaplygin gas model .
the joint pdf peak is shown by the large dot , the credible regions of @xmath105 ( @xmath106 ) by the red dotted line , the @xmath107 ( @xmath108 ) in blue dashed line and the @xmath109 ( @xmath110 ) in green dashed - dotted line .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the joint pdf peak is shown by the large dot , the credible regions of @xmath105 ( @xmath106 ) by the red dotted line , the @xmath107 ( @xmath108 ) in blue dashed line and the @xmath109 ( @xmath110 ) in green dashed - dotted line .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the joint pdf peak is shown by the large dot , the credible regions of @xmath105 ( @xmath106 ) by the red dotted line , the @xmath107 ( @xmath108 ) in blue dashed line and the @xmath109 ( @xmath110 ) in green dashed - dotted line .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the joint pdf peak is shown by the large dot , the credible regions of @xmath105 ( @xmath106 ) by the red dotted line , the @xmath107 ( @xmath108 ) in blue dashed line and the @xmath109 ( @xmath110 ) in green dashed - dotted line .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]      in figures [ figsalphaa ] , the gcgm with @xmath46 free parameters has a joint pdf peak of @xmath113 at @xmath114 with @xmath115 ( @xmath116 ) , @xmath117 ( @xmath118 ) and @xmath119 ( @xmath120 ) credible region contours with pdf values of @xmath121 , @xmath122 and @xmath123 , respectively .
the two - dimensional credible regions are not gaussian - shaped , and their extremes give the estimations @xmath124 at @xmath115 , @xmath125 at @xmath117 , and @xmath126 at @xmath127 confidence levels .
the symbol ( @xmath128 ) in the upper credible value for @xmath16 is due to the fact that the graphics is limited to @xmath129 .
the marginalization ( integration ) of @xmath130 over the @xmath0 space yields @xmath131 , shown in figure [ figsalpha ] , where the credible intervals are @xmath132 at @xmath115 and @xmath133 at @xmath117 ( the @xmath127 pdf level can not be seen with @xmath134 ) .
the pdf peak of @xmath135 at @xmath136 is greater than the pdf value of @xmath137 for @xmath21 ( ordinary chaplygin gas model ) , so the region with pdf greater than @xmath137 has a cdf of @xmath138 , i.e. , we can say that @xmath21 is ruled out at a @xmath139 ( @xmath140 )  confidence level ( or it is preferred at @xmath141 ) .
table [ tableparestgcg ] shows that positive values of @xmath16 are more likely at a @xmath142 ( @xmath143)confidence level .
the @xmath45-parameter case with @xmath144 has a joint pdf peak of @xmath145 at @xmath146 with @xmath115 , @xmath117 and @xmath119 credible region contour levels of @xmath147 , @xmath148 and @xmath149 , respectively .
the two - dimensional credible regions give the estimations @xmath150 at @xmath151 , @xmath152 at @xmath117 , and @xmath153 at @xmath154 confidence levels .
the estimation based on @xmath131 yields @xmath155 at @xmath115 and @xmath156 at @xmath117 confidence levels . the pdf peak of @xmath157 at @xmath158 is greater than the pdf value of @xmath159 for @xmath21 such that @xmath160 is ruled out at a @xmath161 ( @xmath162 )  confidence level .
positive values of @xmath16 are preferred at a @xmath163 ( @xmath164 )  confidence level .     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines and the @xmath167 ( @xmath168 ) credible regions are given by blue dashed lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines and the @xmath167 ( @xmath168 ) credible regions are given by blue dashed lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines and the @xmath167 ( @xmath168 ) credible regions are given by blue dashed lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines and the @xmath167 ( @xmath168 ) credible regions are given by blue dashed lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]
the joint pdf has a peak value of @xmath169 at @xmath170 with @xmath115 , @xmath117 and @xmath119 credible region contour levels of @xmath171 , @xmath172 and @xmath149 , respectively .
these levels indicate that the joint two - dimensional pdf is strongly peaked .
the two - dimensional credible regions give the estimations @xmath173 at @xmath115 , @xmath174 at @xmath117 , and @xmath175 at @xmath127 confidence levels .
the estimations based on one - dimensional pdf yield the same peak value ( due to the highly peaked joint pdf ) , with @xmath176 at @xmath151 and @xmath177 at @xmath117 confidence levels .
the pdf peak of @xmath178 at @xmath179 is greater than the pdf value of @xmath180 for @xmath21 , and we can say that @xmath21 is preferred at a @xmath181 ( @xmath182 )  confidence level .
positive values of @xmath16 are preferred at a @xmath183 ( @xmath184)confidence level .
the joint pdf has a peak value of @xmath185 at @xmath186 with @xmath115 , @xmath117 and @xmath119 credible region contour levels of @xmath187 , @xmath188 and @xmath189 , respectively , which show a highly peaked behaviour .
the two - dimensional credible regions give the estimations @xmath190 at @xmath151 , @xmath191 at @xmath117 , and @xmath192 at @xmath127 confidence levels .
the estimations based on one - dimensional pdf give @xmath193 at @xmath115 and @xmath194 at @xmath117 confidence levels .
the pdf peak of @xmath195 at @xmath196 is greater than the pdf value of @xmath197 for @xmath21 , and we can say that @xmath21 is preferred at a @xmath198 ( @xmath199 )  confidence level .
positive values of @xmath16 are preferred at a @xmath200 ( @xmath201 )  confidence level .
the joint pdf has a peak value of @xmath202 at @xmath203 with @xmath115 , @xmath117 and @xmath119 credible region contour levels of @xmath204 , @xmath205 and @xmath149 , respectively .
the two - dimensional credible regions give the estimations @xmath206 at @xmath115 , @xmath207 at @xmath117 , and @xmath208 at @xmath127 confidence levels .
the estimations based on one - dimensional pdf give @xmath209 at @xmath115 and @xmath210 at @xmath117 confidence levels .
the pdf peak of @xmath211 at @xmath212 is greater than the pdf value of @xmath213 for @xmath21 , and we can say that @xmath214 is preferred at a @xmath215 ( @xmath216 )  confidence level .
positive values of @xmath16 are preferred at a @xmath217 ( @xmath218 )  confidence level .
the joint pdf has a peak value of @xmath219 at @xmath220 with @xmath115 , @xmath117 and @xmath119 credible region contour levels of @xmath221 , @xmath222 and @xmath149 , respectively .
the two - dimensional credible regions give the estimations @xmath223 at @xmath115 , @xmath224 at @xmath117 , and @xmath225 at @xmath127 confidence levels .
the estimations based on one - dimensional pdf give @xmath226 at @xmath115 and @xmath227 at @xmath117 confidence levels .
the pdf peak of @xmath228 at @xmath229 is greater than the pdf value of @xmath230 for @xmath21 , and we can say that @xmath21 is preferred at a @xmath161 ( @xmath162 )  confidence level .
positive values of @xmath16 are preferred at a @xmath231 ( @xmath232 )  confidence level .
the non - gaussian shape of the credible region contours in the @xmath233 two - dimensional parameter space is the culprit of this difference for the estimations of @xmath16 , so it is crucial to emphasize how the @xmath16 parameter is estimated : from a joint two - dimensional pdf or an one - dimensional pdf ( after marginalization of the joint two - dimension pdf ) .
the most robust parameter estimation is based on the one - dimensional pdf , therefore giving small importance to regions in the joint two - dimensions which have a small joint pdf volume ( cdf ) .
comparing the estimations , the ones based on two - dimensional parameter space always give larger values for @xmath29 , after marginalization the pdf is peaked at lower values of @xmath234 .
we also expect that a large sample of sne ia will narrow the credible regions of @xmath16 , so the issue of a ordinary chaplygin gas model ( cgm ) being preferred or ruled out will be probably better addressed .
the generalized chaplygin gas model  ( gcgm ) estimations for @xmath0 are given by table [ tableparestgcg ] , which if compared with the results of the cgm ( table @xmath46 of ref .
@xcite ) clearly indicate that the pdf peak is closer to @xmath40 and the dispersion is increased when the gcgm is considered .
this behaviour is due to the shape of the joint pdf in the parameter space @xmath235 , figures [ figsalphaa ] explicitly show that as @xmath16 increases there is a narrower region of high pdf values closer to higher values of @xmath0 .
but the estimations for @xmath0 shown in table [ tableparestgcg ] and figures [ figsa ]  have to be interpreted .
they do not indicate that the @xmath236 limit ( @xmath40 ) is favoured because the resolution for the @xmath38 sampling is @xmath237 , i.e , a pdf  peak at @xmath238 means the peak happens for @xmath239 indeed .
therefore , in figures [ figsa ] , the graphics for the pdf of @xmath0 which are peaked at @xmath240 have the shape of the graphics for the case @xmath112 , with the pdf peak closer to @xmath40 .
@xcite for a detailed comparison between the cgm and @xmath241 cases .     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath242 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath242 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath242 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]     for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath165 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath242 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines .
the cases for @xmath111 are not shown here because they are similar to the ones with @xmath112 . ]
the analyses of @xmath16 presented above have considered the joint pdfs for the two - dimensional parameter space of @xmath235 , and the extremes of the credible regions have given estimations for @xmath235 at @xmath151 , @xmath117 , and @xmath154 confidence levels .
these estimations for @xmath0 are almost always different from the estimations based on the one - dimensional pdf of @xmath0 , obtained by marginalizing @xmath245 over the @xmath16 space .
we stress again that the most robust estimations are based on the one - dimensional pdfs .
the general behaviour of table [ tableparestgcg ] is that the pdf  peak is at @xmath246 and the @xmath117 credible interval begins at @xmath247 and finishes at @xmath248 .      the parameter estimation of @xmath24 shown in table [ tableparestgcg ] is almost the same for all cases , around @xmath249 with a narrow credible region of @xmath250 .
this estimation , like the one for the cgm in table @xmath46 of ref .
@xcite , is compatible with other sne ia analyses .
table [ tableparestgcg ] and figure [ figsomegas ] show that the curvature density parameter @xmath25 estimation for the gcgm is almost like for the cgm of ref .
@xcite , i.e. , it also strongly depends on the cold dark matter content of the gcgm .    in the @xmath46-parameter case
a closed universe ( @xmath251 ) is preferred at @xmath252 ( @xmath253 ) and a spatially flat universe is ruled out at @xmath254 ( @xmath255 ) confidence level , i.e. , the pdf of @xmath256 is smaller than the pdf of @xmath257 and this region has a cdf of @xmath254 ( @xmath255 ) .
we can also see the figure [ figomegam0c0 ] that the integral of the joint pdf over the upper region gives the same high likelihood of a closed universe , although the joint pdf peak is located in the open universe region .     for the generalized chaplygin gas model , where @xmath258 is a integral of @xmath259 over the @xmath260 parameter space .
it shows a different shape for the confidence regions with respect to the @xmath241 but similar to the cgm of ref .
the joint pdf peak has the value @xmath261 for @xmath262 ( shown by the large dot ) , the credible regions of @xmath263 ( @xmath106 , shown in red dotted line ) , @xmath264 ( @xmath108 , in blue dashed line ) and @xmath109 ( @xmath265 , in green dashed - dotted line ) have pdf levels of @xmath266 , @xmath267and @xmath268 , respectively .
as @xmath269 , the probability for a spatially flat universe is on the line @xmath270 , above it we have the region for a closed universe ( @xmath271 , @xmath251 ) , and below , the region for an open universe ( @xmath272 , @xmath273 ) . ]    but assuming the hypothesis that @xmath112 , a @xmath45-parameter case , the behaviour is opposite , a closed universe ( @xmath251 ) is preferred at @xmath274 ( @xmath275 ) and a spatially flat universe is preferred at @xmath276 ( @xmath277 ) confidence level , i.e. , the pdf of @xmath256 is smaller than the pdf of @xmath278 and this region has a small cdf of @xmath279 .
therefore , if we impose a low value for @xmath280 then the spatially flat universe is favoured .    ,
@xmath280 and @xmath27 for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath281 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath282 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines . as @xmath283 , for @xmath111 we have @xmath284 , for @xmath112 then @xmath285 and for @xmath286 we also have @xmath287 . ]    , @xmath280 and @xmath27 for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath281 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath282 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines . as @xmath283 , for @xmath111
we have @xmath284 , for @xmath112 then @xmath285 and for @xmath286 we also have @xmath287 . ]    , @xmath280 and @xmath27 for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath281 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath282 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines . as @xmath283 , for @xmath111 we have @xmath284 , for @xmath112 then @xmath285 and for @xmath286 we also have @xmath287 . ]    , @xmath280 and @xmath27 for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath281 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath282 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines . as @xmath283 , for @xmath111 we have @xmath284 , for @xmath112 then @xmath285 and for @xmath286 we also have @xmath287 . ]    , @xmath280 and @xmath27 for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath281 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath282 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines . as @xmath283 , for @xmath111
we have @xmath284 , for @xmath112 then @xmath285 and for @xmath286 we also have @xmath287 . ]    , @xmath280 and @xmath27 for the generalized chaplygin gas model .
the solid lines are the pdf , the @xmath281 ( @xmath166 ) regions are delimited by red dotted lines , the @xmath282 ( @xmath168 ) credible regions are given by blue dashed lines and the @xmath243 ( @xmath244 ) regions are delimited by green dashed - dotted lines . as @xmath283 , for @xmath111
we have @xmath284 , for @xmath112 then @xmath285 and for @xmath286 we also have @xmath287 . ]
the generalized chaplygin gas model favours small values of @xmath26 like the ordinary cgm , see table [ tableparestgcg ] and table @xmath46 of ref .
the graphics for @xmath26 in figure [ figsomegas ] , for the @xmath46-parameter case and @xmath45-parameter case with @xmath144 , shows a estimation for @xmath26 peaked at @xmath288 at a high level of confidence , with @xmath289 ( for the @xmath46-parameter case ) and @xmath290 ( for @xmath45-parameter case with @xmath144 ) at @xmath115 ( @xmath166 ) confidence level .
as @xmath26 was sampled with bin size @xmath291 , then the peaks can be somewhere between @xmath292 and @xmath237 .
the joint pdf for the parameter space of @xmath293 , i.e , @xmath294 , is shown in figure [ figomegam0c0 ] .
the pdf peak is now located at @xmath295 , with the same the central value of @xmath26 after bayesian marginalization .
but , in this two - parameter space with simultaneously values of @xmath296 , the credible regions would be given by @xmath297 at @xmath115 ( @xmath106 ) and @xmath298 at @xmath117 ( @xmath299 ) confidence levels .
of course , the integral ( marginalization ) of @xmath300 over the @xmath27 space yields @xmath301 , shown in figure [ figsomegas ] .
the conclusion is that the marginalization decreases the confidence regions , even if it does not change the peak value of probability .
the generalized chaplygin gas density parameter @xmath27 estimation ( table [ tableparestgcg ] ) is like the estimation in the case of cgm shown by table @xmath46 of ref .
@xcite , and is quite spread , i.e. , it has a large credible region , for example @xmath302 in the case of @xmath46-parameter space .
figure [ figsomegas ] show the pdf of @xmath303 for the @xmath46-parameter space ( showing @xmath304 at @xmath115 confidence level )  and the @xmath45-parameter space with @xmath144 ( when @xmath305 ) .
as discussed above for the gcgm , the joint pdf for the parameter space of @xmath306 in figure [ figomegam0c0 ] shows  that the credible regions would be given by @xmath297 at @xmath115 ( @xmath106 ) and @xmath298 at @xmath117 ( @xmath299 ) confidence level .
so , the marginalization of @xmath307 over the @xmath26 space gives @xmath308 , yielding a different peak value as well as different ( narrower ) confidence regions .
the estimates for the dynamical age of the universe today , @xmath71 , given in table [ tableparestgcg ] , show that the generalized  chaplygin gas model increases the @xmath117 credible regions ( e.g. , @xmath309 ) , when compared to the estimations of the cgm in table @xmath46 of ref .
@xcite ( e.g. , @xmath310 ) .
the peak values of @xmath311 range between @xmath312 to @xmath313 , while the cgm features @xmath71 peaked between @xmath314 to @xmath315 .
we expect that a large sample of sne ia will narrow the credible regions and allow , by means of independent estimations of @xmath71 from different observations ( globular agglomerates , etc ) , to verify which cosmological model is preferred : gcgm , cgm and @xmath316cdm .
table [ tablebestfitgcg ] shows the @xmath71 for the best - fitting parameters without any marginalization ( i.e. , the @xmath101-parameter model has its @xmath101 parameters taken simultaneously yielding a minimum @xmath80 ) : the @xmath71 ranges from @xmath317 to @xmath318 , well inside the estimated credible regions for @xmath71 .
all the estimations based on the generalized chaplygin gas model ( table [ tableparestgcg ] ) suggest an accelerating universe today , i.e. , @xmath320 , with a high confidence level ranging from @xmath277 ( @xmath321 ) to @xmath322 ( @xmath323 ) , depending on the curvature and cold dark matter densities of the universe . the peak values of @xmath319 range between @xmath324 to @xmath325 . the data of cgm in table @xmath46 of ref .
@xcite show a similar behaviour .    although the @xmath319  values of table [ tablebestfitgcg ] , from @xmath326 to @xmath327 , are inside the estimated credible regions for @xmath319 , the estimated peak values of @xmath319 differs from the @xmath319 for the best - fitting parameters .
we have also estimated the value of the scale factor that shows the beginning of the recent accelerating phase of the universe .
table [ tableparestgcg ] shows that the most probable value for @xmath74 ranges between @xmath328 to @xmath329 , where @xmath99 is the value of scale factor today , with usually an small upper credible region and a larger lower credible region .
the probability of @xmath330 is theoretically the same having an accelerating universe today ( @xmath320 ) , and these independent probability calculations agree almost exactly as table [ tableparestgcg ]  shows , indicating that the bayesian probability analyses of this work are accurate and reliable for the given sne ia data .
we can obviously express the scale factor value @xmath74 in terms of redshift value @xmath76 , @xmath331 or @xmath332 .
for example , @xmath333 corresponds to @xmath334 .
we expect that a large sample of sne ia with large redshift values will help improve this type of estimation of @xmath76 or equivalently @xmath74 .
the @xmath74 values are almost equal in table [ tablebestfitgcg ] for the best - fitting parameters , between @xmath335 and @xmath336 , and this range fits well inside the the estimated credible regions for @xmath74 of table [ tableparestgcg ] .
the probability of having a eternally expanding universe , given here by demanding @xmath337 during the future evolution of the universe , is estimated to be exactly one or almost one ( at almost a @xmath338 level ) , see table [ tableparestgcg ] .
the estimations for the cgm and @xmath339cdm cases ( tables @xmath45 and @xmath46 , ref .
@xcite ) were miscalculated , they are also exactly one or almost one .
in this paper we investigated some constraints on the generalized chaplygin gas model ( gcgm ) using type ia supernovae data and bayesian statistics .
the model is built with the chaplygin gas and a pressureless fluid as the matter content of the universe , whose density parameters today are @xmath27 and @xmath26 respectively , such that the curvature term is @xmath340 .
the parameter @xmath0 related to the sound velocity of the chaplygin gas , the hubble parameter @xmath24 and the @xmath29 parameter are also independent free variables .
we allowed the parameter @xmath16 to vary between negative and large positive values . in principle , this could lead to instability problems at perturbative level ( @xmath29 negative ) , or superluminal sound velocity ( @xmath16 greater than one ) , but both potential problems can be resolved using a fundamental description for the gcg .
when the five independent parameters are free , the bayesian estimation yields @xmath341 , @xmath342 , @xmath343 , @xmath344 , @xmath302 , @xmath345 , @xmath36 , @xmath346 and @xmath333 where @xmath311 is the age of the universe , @xmath319 is the value of the deceleration parameter today and @xmath74 is the value of the scale factor that marks the start of the recent accelerating phase of the universe .
we have also estimated the cosmological parameters fixing the pressureless matter ( through the constraints of nucleosynthesis ) or the curvature ( imposing the flat condition ) , or even both .
the estimations of the remaining free parameters do not change dramatically in each case , with some exceptions .
for example , the prediction for the curvature parameter depends on whenever the pressureless matter is fixed or not : a spatially flat universe is favoured only if @xmath26 is fixed at low values .
in particular , we verify that the @xmath16 parameter is strongly related to the @xmath0 parameter . in this manner ,
the best value of @xmath16 have the peak values and the credible regions changed in function of the marginalization process used .
the @xmath16 estimations are very spread and differ a lot if one or two - dimensional parameter space is used .
many authors @xcite@xcite , prefer to estimate @xmath16 and @xmath96 simultaneously using two - dimensional credible regions and they do not marginalize the pdf to obtain the more statistically robust one - dimensional pdfs for @xmath16 and @xmath0 , so any comparison should be taken carefully as the credible region estimations on two - dimensional parameter space are usually different from the one - dimensional estimations .    here
we take the marginalization approach to obtain one - dimensional pdfs , and one important point arises in this analysis : despite the dependency on the @xmath0 parameter , the value @xmath21 , that represents the ordinary cgm , is not ruled out since this particular case is preferred at @xmath347 to @xmath348 , depending on the spatial section and matter content .
we calculated the value of the scale factor when the dark energy starts to prevail , @xmath349 , for the case of five free parameters .
in general , the value for @xmath74 ranges between @xmath350 and @xmath351 , where @xmath99 is the value of the scale factor today .
this results are coherent with the suposition that the recent acceleration of the universe could be started with the redshift parameter @xmath20 between @xmath288 and @xmath2 obtained by the type ia supernovae data .    in the case of the ordinary matter
parameter @xmath26 our analyses show that the preferred values are small , around zero .
this situation is the same of the case of the ordinary chaplygin gas .
the result reinforces the idea of the quartessence model @xcite , that unifies dark matter and dark energy into a single fluid , the chaplygin gas .
the predicted age of the universe is well above the age of globular clusters . at same time , the deceleration parameter is consistent with the estimations from the wmap data .
further developments of the present work include the use of a large sample of type ia supernovae .
however , the sample of @xmath88 sne ia used here has an excellent quality with respect to the low values of @xmath91 .
moreover , the dispersion of the parameter estimations obtained here has the same order of estimations using larger samples . in some cases
our estimation are even less spread , as our preliminary investigation on larger sne ia indicates .
we also plan to apply the bayesian analysis to the gcgm using large scale structures data ( 2dfrgs ) and to the anisotropy of the cosmic microwave background . in this way
, we may consistently cross the estimations of different observational data and constrain more sharply the predictions to each free parameter of the model .    *
acknowledgments * we would like to thank orfeu bertolami , ioav waga and martn makler for many discussions about the generalized chaplygin gas model .
the authors are also grateful by the received financial support during this work , from cnpq ( j. c. f. and s. v. b. g. ) and facitec / pmv ( r. c. jr . ) .


growing body of experimental information on fixed - target nucleus - nucleus collisions from ags and cern - sps accelerators provides a reliable basis for the systematic investigation of strongly interacting hadronic matter under extreme conditions .
particularly , production of particles containing strange quarks is expected to reflect the reaction dynamics at the early stage of collisions . in this context ,
@xmath0 mesons , particles consisting mainly of @xmath4 pairs , are of a special interest .
since the interaction of @xmath5 mesons with non - strange hadronic matter is suppressed according to ozi rule , @xmath5 mesons are expected to decouple easily from the hadronic fireball .
several calculations have been done for @xmath0 meson properties in dense hadronic environment ( mainly at higher baryon densities ) @xcite .
it has been found that modifications of the @xmath0 width and mass are sensitive to the strangeness content of the surrounding medium .
recently , new experimental data from cern - sps on @xmath0 meson production in central pb+pb collisions at 158gev / a beam energy became available .
advantage of the rich experimental program at cern is that it allows to study the @xmath0 production via different @xmath0 decay channels .
results of na49 collaboration  @xcite are based on the @xmath0 meson identification via its hadronic decay @xmath6 , while na50 collaboration has recently reported @xcite about preliminary analysis of @xmath5 mesons identified via the dileptonic decay channel @xmath7 .
the purpose of this paper is to show , that the @xmath5 meson production spectra reconstructed via @xmath2 and @xmath1 decay channels can be significantly different up to the level of present experimental observations .
after introductory considerations in section  ii we derive in section  iii expressions for momentum distributions of @xmath0 mesons detected via @xmath8 and @xmath9 channels .
we take into account rescattering of decay kaons and in - medium modification of @xmath0 meson properties discussed in section  iv .
numerical results are presented in sections  v and vi and conclusions are drawn in section  vii .
for our estimates of the observed @xmath5 meson spectra we assume a simple two - stage picture of central heavy ion collisions , which is close to those considered within cascade - transport @xcite , hydrodynamical @xcite and thermodynamical approaches @xcite .
initial stage of a collision is characterized by temperature close to the qcd phase transition @xmath10  mev .
then the system expands up to the point , when numbers of different kinds of particles freeze in - a chemical freeze - out .
thermodynamical parameters of this stage could be obtained by fitting the final total hadron multiplicities .
typical temperature is found to be @xmath11  mev . during the second stage of the expansion ,
elastic scatterings change momentum distributions of hadrons until they cease and distributions freeze in .
the freeze - out temperature can be extracted from a simultaneous fit to the single - particle @xmath12-spectra of different particles supplemented by the analysis of particle correlation data . according to refs .
@xcite one has @xmath13  mev .
in our considerations we assume that the fireball created in heavy - ion collision consists , mainly , of pions , kaons and excited mesonic resonances in central rapidity region .
mean free path @xmath14 of @xmath0 mesons in hadron gas is estimated in ref .
comparison with mean free paths of pions and kaons , @xmath15 , from ref .
@xcite gives @xmath16 for temperatures @xmath17 .
hence , one can expect that @xmath0 mesons decouple from the pion - kaon subsystem at some earlier stage between the chemical and thermal freeze - out . upon this stage @xmath0
mesons stream freely out from the fireball .
pions and kaons , on the other hand , may still participate in mutual secondary interactions up to the stage of total thermal freeze - out .
therefore , if a @xmath0 meson decays inside a fireball via hadronic channel we have to take into account possible interaction of its decay products with surrounding hadronic environment .    as a result of secondary interactions ,
the reconstructed invariant mass of a given pair of @xmath0 daughter kaons falls out from the original @xmath0 meson peak into the region identified as a combinatorical background .
therefore , the @xmath0 mesons decaying in medium , can be partially unrecognized in experimental analysis of @xmath2 pairs .
together with negligible final state interaction of secondary dimuons originating from @xmath18 decays this behaviour results in the relative suppression of @xmath0 meson yield observed via hadronic @xmath8 channel .
such mechanism has been quantitatively studied in ref .
@xcite where suppression at the level 4060% has been obtained from simulation using rqmd code .    in this paper
we discuss effects which may enhance the suppression of observed @xmath5 mesons identified via kaon channel .
possible increase of a @xmath0 meson width in medium will enlarge the probability of @xmath5 decay inside a fireball enhancing thus consequences of mechanism studied in ref .
alternatively , we consider a possibility that the kaon decay channel of a @xmath0 meson becomes kinematically quenched in the medium .
additionally we argue that substantial relative difference in properties of @xmath19 and @xmath20 meson in hadronic environment caused by the isospin asymmetry and/or by the large baryonic admixture , would also prevent the reconstruction of @xmath0 mesons decaying into kaon pairs inside the medium .
let us denote the phase - space distribution of @xmath0 mesons in the center - of - mass system of two colliding nuclei at freeze - out as @xmath21 .
then the primary momentum distribution of @xmath0 mesons is given by [ integr ] _ 0(p ) = _
d^3 ^p _ f _ ( _ , ) , where integration goes over the fireball volume within a freeze - out hyper - surface @xmath22 ( surface normal vector @xmath23 contracted with a @xmath0 meson momentum @xmath24 ) @xcite . here
we assume the case of a time - like freeze - out hyper - surface , @xmath25 , which is relevant for applications below ( for discussions of alternative cases cf .
@xcite ) . in the absence of any in - medium modifications of decay products and final state rescattering
the shape of _ observed _ momentum distributions @xmath26 of muon pairs @xmath27 and @xmath28 of kaon pairs @xmath29 would be the same . here
@xmath30 is @xmath0 meson total width and @xmath31 , @xmath32 are the @xmath5 partial decay widths in kaon and muon decay channels . accordingly , in the experimental analysis
, @xmath0 meson distribution is reconstructed from momentum distribution of decay products multiplied by the corresponding inverse branching ratio .    in this section
we derive expressions for the apparent momentum distributions of @xmath0 mesons reconstructed via kaon and dimuon decay channels taking into account possible modifications of meson properties in medium and consequences of kaon rescattering .
let us consider @xmath5 meson suffered the last interaction at position @xmath33 inside the hadronic fireball .
the probability , that @xmath0 meson lives for a time @xmath34 is [ pbx ] d_(t)&= & , where @xmath35 is total width of a moving meson with energy @xmath36 .
asterisk * denotes in - medium values of the quantities .    after traveling for a given time @xmath34 with velocity @xmath37
the @xmath5 meson decays at the position @xmath38 with a probability @xmath39 into two kaons . in - medium values of widths
@xmath40 are determined by the current local temperature and density of the system .
daughter kaons from @xmath5 meson decay have momenta @xmath41 in the rest frame of the @xmath0 meson .
value of @xmath42 follows from equation @xmath43 where @xmath44 and @xmath45 are in - medium spectra of kaons and anti - kaons . in the center - of - mass system of two colliding nuclei ( cms )
the kaon momenta are equal to [ kmoment ] _
k^&= & , + & = & p_kk \ { _ k+ _ ( _ -1 ) ( _ k _ ) } , where @xmath46 . the unit vector @xmath47 is uniformly distributed in the @xmath0 rest frame , direction of @xmath48 is determined by @xmath0 meson momentum @xmath49 in cms .    for a successful identification of @xmath5 mesons in invariant mass spectrum of observed @xmath2 pairs
it is essential that momenta of daughter kaons do not change while leaving the hadronic fireball .
here we investigate two mechanisms which may change momenta of daughter kaons : rescattering in surrounding hadronic environment and change of momentum due to in - medium @xmath50 meson potential .
probability that secondary kaon and anti - kaon leave a fireball without rescattering is determined by their mean free paths @xmath51 , the time which kaons need to reach the fireball border @xmath52 , and their velocities @xmath53 ( in cms ) as follows : [ pl]p_(t ) & = & .thus , the probability to register @xmath0 meson in the kaon channel when it decays in medium can be expressed as : [ probx ] p_i(1-p_(t))+p_ii_(t ) , where @xmath54 and @xmath55 are probabilities to identify a @xmath0 meson from rescattered ( @xmath54 ) and non - rescattered ( @xmath55 ) kaons .
single rescattering of one of the kaons changes momentum of the kaon pair by vector @xmath56 , with @xmath57 being a average thermal momentum of pions .
correspondingly the invariant mass changes by @xmath58 . at high temperatures ( @xmath59 )
average thermal momentum of pions is @xmath60  mev and hence a single rescattering may shift @xmath61 far from the @xmath5 meson mass . in this case
we put @xmath62 .
neglecting in - medium effects , one takes @xmath63 , assuming that without a hard rescattering all kaon pairs from @xmath0 decays can be identified .
however the modification of particle properties can provide another mechanism preventing @xmath0 identification
. in medium kaons feel an effective mean - field potential inducing a change of their spectra @xmath64 .
leaving the fireball , kaons have to come back to their vacuum mass shell .
this changes the invariant mass and momentum of the pair . going out of medium
, kaon stays on the same energy level , and its momentum outside the fireball becomes @xmath65 ( @xmath66 ) , where the kaon energy @xmath67 is evaluated at the moment of @xmath0 decay . as a result ,
momentum of the kaon pair is changed by @xmath68 , where we put @xmath69 and @xmath70 .
the invariant mass shift is then given by @xmath71 .
in the nucleon free , isospin symmetrical meson gas , in - medium spectra of kaons and anti - kaons are identical @xmath72 and the invariant mass shift reduces to @xmath73 $ ] .
this expression can be interpolated between the limit cases @xmath74 and @xmath75 as @xmath76 . therefore in the symmetrical case @xmath77 and
@xmath78 vanish for small @xmath79 .
moreover in the isospin symmetrical case kaons suffer only a small mass modification ( @xmath8030  mev for @xmath59 ) .
hence , the release of kaons from a small potential well formed inside a fireball does not affect the momentum and invariant mass of a pair strongly enough to prevent a @xmath0 meson reconstruction .
for this case we take @xmath81 .    in the case of rather strong isospin asymmetry and/or significant baryonic admixture @xmath19 and @xmath20 mesons
can have quite different spectra in medium  @xcite .
this may lead to a wide spread of the kaon pair invariant masses , making the reconstruction of @xmath0 mesons decaying inside fireball impossible ( @xmath82 ) .
note , that in this case our results are insensitive to the details of kaon propagation in medium .
they depend only on the total @xmath0 meson width via function @xmath83 .
finally the probability that a @xmath0 meson created at position @xmath84 will be detected via the kaon channel is : & & _
0^ t^*_kk(t ) d_(t ) p_(t ) p_ii + [ kaonprob ] & & = _ ( _ ) + _ 0^ _ d t^*_kk(t ) d_(t ) p_(t ) p_ii , where @xmath85 is the time of flight of a @xmath0 meson through the fireball , for @xmath86 we have @xmath87 and @xmath88)$ ] .
momentum distribution of the kaon pairs from @xmath0 decays is then obtained by averaging expression ( [ kaonprob ] ) over all @xmath47 directions in the @xmath0 meson rest frame , and integrating over the fireball volume with a primary @xmath0 distribution - see eq.([integr ] ) .
then , multiplying by the inverse branching ratio we express `` reconstructed '' @xmath0 distribution in kaon channel as : _
k(p ) & = & _ ( _ ) + _ 0^_t d_(t ) _
kk^*(t ) p_ii , + [ ratrescat ] & & = _ ( t ) . here and the brackets denote integration :  =
_ d^3 ^p _ f _ ( _ , ) (  ) .
note that according to this definition @xmath89 .
consideration of the di - muon pair momentum distribution is more straightforward . with a change of total @xmath0 width in medium the branching ratio of the @xmath90 decay @xmath91 changes too .
therefore , the apparent @xmath0 distribution reconstructed in the muonic channel is : [ ratmu ] _ (
p ) & = & _ ( _ ) + _ tot^0 _ 0^ _ dt d_(t ) .
note that for @xmath92 we have @xmath93 , for @xmath94 we have @xmath95 and if @xmath96 then @xmath97 .    in our numerical estimates we will consider ratio of the reconstructed @xmath5 meson distributions eq.([ratrescat],[ratmu ] ) which does not depend on the overall normalization : [ eta ] ( p)=_k(p)/_(p ) .    for a direct comparison with experimental results on @xmath12 distributions of identified @xmath5 mesons , one has to integrate over the rapidity interval accessible to the experiments , utilizing @xmath98 .
dependence of the ratio ( [ eta ] ) on @xmath12 reads [ etamt ] ( m_t)=<_k(p)>_y/<_(p)>_y .
we also define ratio of the apparent and primary @xmath0 momentum distributions for @xmath2 and @xmath1 decay channels as a function of @xmath12 : _
k,(m_t)=<_k,(p)>_y/<_0(p ) > _ y
main hadronic decay channels of @xmath0 meson are @xmath100 and @xmath101 .
let us now consider the change of a @xmath0 decay width in a hot meson gas due to the modification of kaon , pion and @xmath102-meson properties . in - medium properties of pions can be effectively incorporated by a small mass shift @xmath103 with @xmath104 .
similar behaviour is expected for kaons .
spectral function of @xmath102 meson in medium was extensively investigated in the context of di - lepton production in heavy - ion collisions at sps energies @xcite . in baryonic matter , coupling of @xmath102 mesons to resonance
nucleon - hole modes  @xcite together with the modification of pions  @xcite play dominant role .
the @xmath102 meson becomes very broad in medium and its spectral density strength is driven effectively to lower energies in analogy to brown - rho - scaling picture  @xcite . in purely mesonic systems
the @xmath102 mass is found to be almost independent of the temperature due to cancellation of @xmath105- and @xmath106-loop contributions  @xcite .
the @xmath102 width , on the other hand , is expected to increase in meson gas considerably e.g. by 80  mev at @xmath107  mev and by 160  mev at @xmath108  mev  @xcite .
the partial width of the @xmath109 decay in medium depends on the kaon mass @xmath110 as : [ gkk]_kk^*= _
kk^0 p_cm^3(m_^2,m_k^*,m_k^*)/ p_cm^3(m_^2,m_k , m_k ) , where @xmath111 is a kaon momentum in the rest frame of @xmath5 meson decay , obeying the equation : @xmath112 . assuming isospin symmetry we have @xmath113  mev .
then vacuum width is equal to @xmath114  mev .
note that for @xmath115 the in - medium width @xmath116 decreases fast and vanishes for @xmath117  mev .
the second hadronic decay channel @xmath101 is effected by the decrease of the @xmath102 meson mass , increase of the @xmath102 meson width and the bose - einstein enhancement factor for pions . for the @xmath0 meson at rest we write [ grho]_^ * = ^0__(m_^*,_^*,t)/ _ ( m_^0,_^0,0 ) , and [ grho ] _ ( m_,,t ) & = & _ 2m_^m_-m _ [ ( -m_)^2-m_^2]^32 + & & , where @xmath118 and @xmath119 are in - medium @xmath102 meson mass and width , respectively . vacuum values are @xmath120  mev , @xmath121  mev , and @xmath122  mev .
we use the constant width approximation for the @xmath102 meson spectral density and denote @xmath123 .
the bose - einstein distribution of pions with temperature @xmath124 is @xmath125 . in the zero - width limit
we obviously have @xmath126 .
numerical evaluation of eq .
( [ grho ] ) leads to the approximated relation @xmath127\big ) $ ] valid for @xmath128  mev , @xmath129  mev , and @xmath130 .
the shift of a pion mass produces a minor effect and it is , therefore , neglected here .
finally , the total width of @xmath0 is given by @xmath131 . here and below we do not consider the @xmath0 meson mass shift , which is small and cancels as soon as we consider the ratio of the momentum spectra , cf .
eq.([eta],[etamt ] ) .
for completeness we remark that the di - lepton decay channel @xmath132 also suffers a modification in medium , because the vector - meson  photon coupling is suppressed by meson fluctuations  @xcite . for the @xmath133 coupling the suppression factor is determined only by kaon fluctuations @xmath134 . here
@xmath135^{-1}\,$ ] , @xmath136 and @xmath137  mev is the pion decay constant . because of the large kaon mass this correction is negligibly small : @xmath138 .
finally we , specify , how the in - medium properties of kaons and @xmath102 mesons relax to their vacuum values during the fireball expansion .
assuming that @xmath139 and @xmath140 are proportional to the density of the system , we have m_k,(t ) = m_k,^0 , _ ( t ) = _ ^0 , where @xmath141 and @xmath142 are input parameters .
before finishing this section we would like to point out that the estimations done here are valid only for an almost baryon free fireball . in the presence of baryons the calculation of a phi self - energy
becomes more elaborated , cf .
. however , in this case the spectra of kaons and anti - kaons differ dramatically , inhibiting thereby the @xmath0 reconstruction even more .
it corresponds to the case @xmath143 when the reconstructed @xmath0 distributions ( [ ratrescat],[ratmu ] ) are the functions of a total width only . to simulate effectively the in - medium modification of @xmath144 we will use eq .
( [ gkk ] ) and vary the kaon mass .
after the considerations above let us now specify the model of the fireball expansion , which we will use in our numerical calculations . assume a simple homogeneous spherical fireball with the constant density and temperature profiles .
the @xmath0 meson momentum distribution [ dist ] f_(,)=,is determined by the temperature @xmath145 , flow velocity profile @xmath146 and radius @xmath147 .
time of flight of a @xmath0 meson and its daughter kaons through the medium contained in eqs .
( [ pl ] ) and ( [ ratrescat ] ) can be expressed as @xmath148 and @xmath149 , where @xmath150 stands for a time , during which a particle with velocity @xmath151 passes a distance from position @xmath152 to the border of a sphere with the radius @xmath153 .
since the fireball is expanding with radial velocity @xmath154 , this time satisfies the equation @xmath155 .
this implies : _ r ( , ) & = & ( + [ tau ] & -&(-v_fr ) ) ( ^2-v_f^2)^-1 . solution ( [ tau ] ) is valid for @xmath156 and @xmath157 . in the case
@xmath158 we put @xmath159 .    during the expansion @xmath160
, fireball density drops as @xmath161 and the temperature decreases as @xmath162 as expected for relativistic pion gas , cf .
the kaon mean free path is @xmath163 and therefore _
k(t ) & = & _ k^0 r^3(t)/r_0 ^ 3 .to incorporate the freeze - out effect we will assume that as soon as @xmath164 kaons become free and @xmath165 as well as @xmath166 . the freeze - out time is then given by [ taofo ] _
f.o.= ( -1 ) .
parameters @xmath167 , @xmath168 , @xmath154 , and @xmath169 serve as input for numerical evaluations below .
the parameter @xmath170 is used for an effective parameterization of the freeze - out time @xmath171 .
it should not be considered as a true freeze - out temperature , since our estimation is based on the simplified hydrodynamical description of a fireball .
the model set up here is a rather crude approximation .
however , the final results are found to be rather insensitive to the details of hydrodynamical evolution of a fireball , being determined mainly by the values of @xmath172 , @xmath173 , and @xmath154 .
we shall vary the input parameters within a broad range to illustrate different possibilities .
in this section we perform numerical evaluation of our expressions for @xmath5 meson yields reconstructed via @xmath2 and @xmath1 channels in central pb+pb collisions at 158gev / n sps energy .
first we investigate to what extend rescattering of secondary kaons enhanced by the in - medium modification of a @xmath0 meson width can suppress experimentally observed yield of @xmath5 mesons identified via @xmath2 channel .
we remind that results of ref .
@xcite give maximal suppression factor 40% for the @xmath5 meson observation in the kaon decay channel .
in our evaluation we use several combinations of input parameters .
freeze - out temperatures @xmath145 of @xmath5 mesons distributed according to eq.([dist ] ) vary between @xmath174 and @xmath170 : ( @xmath175 ) @xmath176  mev , ( @xmath177 ) @xmath178  mev , ( @xmath179 ) @xmath180  mev .
size of the fireball @xmath147 at the stage of the @xmath0 freeze - out has to be comparable with @xmath0 meson mean free path @xmath181 at given temperature @xmath182 with @xmath183 . for different temperature parameters above
we take according to ref .
@xcite : @xmath184  fm , @xmath185  fm , and @xmath186  fm .
first we evaluate suppression factor ( [ etamt ] ) without any modifications of particle properties in medium , i.e. , @xmath187 . in this case
we have @xmath188 and @xmath189 .
flow velocities corresponding to selected temperatures @xmath145 are adjusted to reproduce the slope of @xmath5 meson @xmath12 distribution measured by the na50 collaboration : @xmath190  mev : @xmath191 , @xmath192 , @xmath193 .
we take @xmath194 and vary the mean free path of kaons @xmath195 within the interval @xmath196 , where @xmath197 follows from estimations of ref .
@xcite : @xmath198  fm , @xmath199  fm , and @xmath200  fm .
we vary also @xmath170 between 100  mev and 80  mev in agreement with analysis  @xcite .
this translates into the interval of freeze - out time values @xmath201 .
all these variations produce the upper grey areas shown in fig .
[ fig : nomed ]    ) as a function of @xmath12 calculated for three sets of parameters ( @xmath202 without inclusion of in - medium modifications of kaons and @xmath102 mesons .
the upper grey area corresponds to variations of kaon mean free paths @xmath203 and freeze - out temperature @xmath170 as described in the text .
the lower grey areas are obtained for the same set of parameters but with the common expansion velocity @xmath204.,width=188 ]    for all three sets of parameters @xmath205 we observe that @xmath206 does not fall below @xmath207 significantly .
this is related to the large expansion velocity of the fireball . in this case @xmath208 and @xmath5 mesons decay after the thermal freeze - out . to illustrate this effect we recalculate ratio @xmath206 for the same three cases fixing @xmath204 what corresponds to @xmath209 .
results obtained ( @xmath210 ) are shown as lower gray areas in fig .
[ fig : nomed ] .    to reproduce results of rqmd calculations described in ref .
@xcite we take somewhat larger size of a fireball with @xmath211 , freeze - out temperature @xmath212  mev and @xmath213  fm .
this corresponds to the lowest limit allowed by the analysis  @xcite .
the results are shown in fig .
[ fig : nomedr15 ] by solid lines .
the limiting scenario considered in @xcite , when the freeze - out volume is determined by the last kaon interactions , can be reproduced with @xmath214  mev .
this case is shown by dash lines in fig .  [
fig : nomedr15 ] .    ) calculated for cases @xmath215@xmath216 with @xmath211 , @xmath217  fm and @xmath213  fm .
solid lines are calculated with @xmath212  mev and dash lines correspond to @xmath214  mev .
no medium effects are included.,width=188 ]    we take solid lines in fig .
[ fig : nomedr15 ] as a reference point for our further investigation of in - medium effects .
first we consider modifications of @xmath102 meson properties .
we choose @xmath218 mass shift to be @xmath219  mev for our three parameter sets .
the @xmath218 meson width depends on the temperature .
relying on ref .
@xcite we take : @xmath220  mev , @xmath221  mev , and @xmath222  mev . comparison with fig .
[ fig : nomedr15 ] ( solid lines ) shows slight decrease of ratio @xmath206 which corresponds to small increase of the total @xmath0 meson width by 40% due to @xmath223 channel .
[ fig : med ] shows results obtained taking into account modification of @xmath50 meson properties in - medium .     calculated for in - medium modification of @xmath0 meson properties .
results for different parameters sets @xmath215@xmath216 are depicted by solid , dash and dotted lines , respectively .
the left plane corresponds to the case when the kaon mass increases in medium @xmath224  mev , whereas the right plane shows results for a decreasing kaon mass @xmath225  mev .
thin lines show suppression factors in the muon @xmath226 ( b ) and kaon @xmath227 ( c ) channels calculated for parameter set ( @xmath175 ) .
upper plots are calculated for the case @xmath228 , lower plots correspond to @xmath143 . in all cases modification of @xmath102 meson properties in medium
is taken into account.,width=302 ]    first we investigate the case when @xmath100 channel is closed initially ( @xmath224  mev ) and it opens only during the fireball expansion .
results are shown in the left part of fig .
[ fig : med ] , where the upper plot is calculated for @xmath228 and the lower one corresponds to a strong suppression of the kaon channel in medium with @xmath143 . comparing ratios @xmath206 in fig .
[ fig : nomedr15 ] and in fig .
[ fig : med ] shown by thick solid line ( a ) calculated for parameter set ( @xmath175 ) we observe that kinematical quenching of @xmath229 channel by increasing kaon mass decreases ratio @xmath206 very slightly for both values of @xmath55 .
this happens because @xmath0 meson width becomes very small in this case and therefore probability of @xmath5 meson decay inside the expanding fireball and consequently also probability for rescattering of daughter kaons are small .
compare lines ( c ) , calculated for parameter set ( @xmath175 ) , with a corresponding line in fig .  [ fig : nomedr15 ] .
lines ( b ) in fig .
[ fig : med ] ( left side ) show that @xmath226 becomes larger than one for increasing kaon mass .
since increase of @xmath230 by 10%20% does not change considerably the effective @xmath12 slope of @xmath231 distribution , we do not need to readjust the flow velocity parameter .
this increase of @xmath230 leads at the end to a small decrease of @xmath232 .
let us now consider the case when the @xmath0 meson width increases strongly in hadronic medium due to the increase of @xmath116 .
we simulate this effect by decrease of the kaon mass in medium , which can result e.g. from rescattering of kaons on pions through @xmath233 and heavier kaonic resonances @xcite . here
we restrict ourselves to rather conservative modification of kaon masses @xmath234  mev@xmath235 which corresponds to the @xmath0 width @xmath236  mev . in this case @xmath237 and at small @xmath238 region @xmath230 can be suppressed up to 4060% .
thus , for a given freeze - out temperature @xmath145 and total width @xmath144 we readjust flow velocity @xmath239 to reproduce the slope of the @xmath12 distribution measured in di - muon channel by na50 @xcite .
for our three sets of parameters ( @xmath240 ) we obtain new flow velocities : ( @xmath175 ) @xmath241 , ( @xmath177 ) @xmath242 , ( @xmath179 ) @xmath243 .
right part of fig .
[ fig : med ] shows our results obtained for @xmath225  mev .
corresponding partial width is @xmath244  mev , and the total width is @xmath245  mev .
this leads to @xmath246 , which provides a strong suppression of @xmath247 as it is shown for the parameter set ( @xmath175 ) in fig .
[ fig : med ] ( right side , lines c ) .
we find @xmath248 for @xmath228 and @xmath249 for @xmath143 .
however , since the ratio @xmath250 also shown in fig .  [ fig : med ] ( lines b ) is also suppressed the resulting ratio @xmath206 remains on the level @xmath251 for small @xmath238 , provided we put @xmath143 and @xmath252 for @xmath228 .
taking even larger values of the total decay width in the most preferable case ( @xmath175 ) we obtain @xmath253 for @xmath254  mev , @xmath242 , and @xmath255 for @xmath256  mev , @xmath257 ( both for @xmath258 ) .
it is instructive to investigate ratio ( [ eta ] ) as a function of the @xmath0 momentum . averaging over rapidity mixes momenta within a broad interval , e.g. , for @xmath259 momenta @xmath260  mev what partially washes out the final suppression effect . to illustrate this point we show ratio @xmath261 as a function of a @xmath0 momentum in fig .
[ fig : pdep ] .    ) as a function of @xmath0 meson momentum .
curve styles correspond to those in fig .
[ fig : med ] with readjusted flow velocities .
curves a are calculated for @xmath262 , whereas sets b are calculated with @xmath62 .
, height=188 ]    for increasing kaon mass and correspondingly vanishing @xmath31 we obtain a decrease of @xmath261 at small momenta by 4050% ( left plane ) .
this is again a direct consequence of a rapid fireball expansion which brings @xmath263 in the integral ( [ ratrescat ] ) fast to its vacuum value . in the extreme case for @xmath264 and @xmath265
we get @xmath266 . for decreasing kaon mass
the reduction is even stronger @xmath26770% .
momentum dependence of @xmath232 at small momenta @xmath268 is completely washed out by rapidity averaging .    in fig .
[ fig : ydist ]     ( full circles ) , @xmath269 ( open circles ) and @xmath270 ( open diamonds ) calculated for the parameter set ( @xmath175 ) with @xmath225  mev , the flow velocity @xmath241 and @xmath62 . ,
width=188 ]    we compare the original rapidity distribution of @xmath0 mesons @xmath271 , with the distributions which can be reconstructed via kaonic @xmath269 , and muonic @xmath270 decay channels .
calculations are done for the parameter set ( @xmath175 ) with @xmath272  mev , flow velocity @xmath241 and @xmath143 .
distributions are normalized to give @xmath273 .
we observe a considerable broadening of the rapidity distribution measured in the kaon channel .
we have studied distributions of @xmath5 mesons in heavy - ion collisions at sps energies reconstructed via hadronic @xmath2 and dilepton @xmath274 decay channels .
the analysis of @xmath5 meson mean free path allows to suppose that @xmath0 mesons decouple from the hadronic system at somewhat earlier stage before the common breakup of the hadronic fireball .
therefore , kaon pairs originated from the @xmath0 decays inside a fireball can be rescattered or absorbed .
such kaon pairs will not contribute to a @xmath0 meson reconstruction , whereas the leptonic probes can leave a fireball freely .
we derive the expressions ( [ ratrescat ] ) and ( [ ratmu ] ) for the apparent momentum distribution of @xmath0 mesons in kaonic and muonic channels respectively .    within a simple model of spherically expanding fireball
we investigate dependence of a relative suppression factor of the hadronic channel with respect to the dileptonic one on parameters of the system and on the @xmath0 meson in - medium properties . for a vacuum @xmath0-meson width @xmath275  mev the maximal suppression 0.60.8
is obtained for the fireball size and expansion time @xmath276  fm .
these values are in agreement with results of rqmd simulations @xcite .
the crucial parameter is the fast expansion of a fireball with @xmath2770.5 corresponding to the @xmath0 freeze - out temperature range @xmath278170  mev .
width of hadronic @xmath0 meson decay channels ( @xmath100 and @xmath279 ) can be modified in medium due to changes of the meson properties .
we have found that increase of the @xmath280 channel width due to the broadening of @xmath102 meson and decrease of @xmath218 meson mass leads alone to a tiny increase of the suppression .
other possibility is kinematical quenching of the kaon decay channel , which we simulate by simultaneous increase of @xmath19 and @xmath20 masses . since
total width @xmath281 of @xmath5 meson in medium becomes small ( increase of the @xmath280 channel width is not strong enough ) @xmath0 mesons decay mainly outside the fireball , where vacuum properties of @xmath5 mesons are restored and rescattering of daughter kaons is negligible .
together with relative amplification of the muon decay channel by 20% the resulting suppression factor found for a quenched kaon decay channel was @xmath252 .
the increase of the @xmath0 meson width in medium provides , on the other hand , a mechanism for strong suppression ( @xmath249 ) of the kaonic detection channel due to the enhancement of @xmath5 decay probability inside a fireball increasing thus rescattering of daughter kaons .
however , increase of the @xmath0 total width reduces simultaneously the branching ratio of @xmath282 decay and suppresses the spectrum of @xmath5 mesons reconstructed via @xmath1 decay channel .
obtained suppression of the muonic decay channel at the level @xmath28360% requires readjustment of the flow velocity to be compatible with experimental slope of na50 . adjusted flow velocity @xmath2840.4 for @xmath285  mev and @xmath278170  mev ( compare to @xmath2770.5 obtained for a vacuum @xmath0 width ) gives final net relative suppression factor of kaon channel to muon channel @xmath251
this value is close to experimental observations at cern sps @xcite .
strong increase of @xmath5 meson in - medium width can take place if the kaon mass decreases in medium by 30  mev .
the mechanism for such kaon mass modification can be similar to that studied in ref .
@xcite .
finally we suppose that to improve understanding of experimental results on @xmath5 meson production in heavy ion collisions at cern sps @xcite further detailed investigations taking into account in - medium effects within transport or hydrodynamical models are necessary .
the authors would like to thank yu.b .
ivanov , k.  redlich , v.d .
toneev , and d.n .
voskresensky for discussions .
work of e.k . was supported by the grant of ect@xmath286 , p.f .
was partially supported by slovak grant agency for sciences grant no . 2/5085/00
.    99 c.m .
ko et al . , phys . rev .
* c45 * , 1400 ( 1992 ) . m.  asakawa and c.m .
ko , nucl .
* a572 * , 732 ( 1994 ) .
ko and d.  seibert , phys .
* c49 * , 2198 ( 1994 ) ; phys .
rev .  * 50 * , r559 ( 1994 ) .
k.  haglin , nucl .
* a584 * , 719 ( 1995 ) .
c.  song , phys .
* b388 * , 141 ( 1996 ) .
w.  smith and k.l .
haglin , phys .
* c57 * , 1449 ( 1998 ) .
f.  klingl , t.  waas , and w.  weise , phys .
v.  friese for na49 , j. phys .
* g23 * , 1837 ( 1997 ) ; f.  sikler for na49 , nucl . phys .
* a661 * , 45c ( 1999 ) ; s.v .
afanasiev et al . , to be published .
n.  willis for na50 , nucl
a661 * , 534c ( 1999 ) ; l.  villatte for na50 , in proc .
xxxviii international winter meeting on nuclear physics , bormio , italy , january 24 - 29 ( 2000 ) .
bass et al .
, prog .  part .
* 42 * , 313 ( 1999 ) .
bratkovskaya , w.  cassing , nucl .
* a619 * , 413 ( 1997 ) .
bravina et al . , phys .
rev .   * c60 * , 044905 ( 1999 ) .
j.  sollfrank et al . , j.  phys .
* g25 * , 363 ( 1999 ) .
p.  filip , e - print : hep - ex/9605001 , nucl - th/9801058 ; p.  filip and j.  pisut , phys .
lett .   * b411 * , 187 ( 1997 ) .
j.  cleymans , k. redlich , and d.k .
srivastava , phys . rev . *
c55 * , 1431 ( 1997 ) . yu.b .
ivanov , e.g. nikonov , and v.d .
toneev , privat communications .
peressounko and yu.e .
pokrovsky , e - print : hep - ph/0009025 .
p.  braun - munzinger , i.  heppe , and j.  stachel , phys .
* b465 * , 15 ( 1999 ) ; i.  heppe , ( diploma thesis , uni .
heidelberg ) .
b.  kmpfer , e - print : hep - ph/9612336 ; b.  kmpfer et al .
, j. phys . * g23 * , 2001 ( 1997 ) .
b.  tomik , u.a .  wiedemann and u.
heinz , e - print : nucl - th/9907096 .
t.  peitzmann , e - print : nucl - th/0006025 .
t.  alber et al .
( na49 coll . )
j. * c2 * , 643 ( 1998 ) .
f.  cooper and g.  frye , phys .
* d10 * , 186 ( 1974 ) .
bugaev , nucl .
* a606 * , 559 ( 1996 ) ; d.  molnar and m.  gyulassy , phys .
rev .  * c62 * , 054907 ( 2000 ) .
haglin and s.  pratt , phys . lett . *
b328 * , 255 ( 1994 ) . h.a .
weldon , z. phys . * c54 * , 431 ( 1992 ) . s.c
.  johnson , b.v .
jacak , and a.  drees , eur .
j. * c18 * , 645 ( 2001 ) . g.e .
brown et al . , phys .
c * 43 * , 1881 ( 1991 ) ; g.e .
brown , c.m .
ko and k.  kubodera , z.  phys .
* a341 * , 301 ( 1992 ) ; m.  lutz , a.  steiner and w.  weise , nucl .
* a574 * , 755 ( 1994 ) ; e.e .
kolomeitsev , b.  kmpfer and d.n .
voskresensky , nucl .
* a588 * , 889 ( 1995 ) ; e.e .
kolomeitsev , b.  kmpfer and d.n .
voskresensky , int .
j.  mod .
* e5 * , 313 ( 1996 ) ; c.m .
ko , j.  phys .
* g27 * , 327 ( 2001 ) ; m.f.m .
lutz and c.l .
korpa , eprint : nucl - th/0105067 ; b.  friman and h.j .
pirner , nucl.phys .
* a617 * , 496 ( 1997 ) ; m.  lutz , b.  friman , and g.  wolf , nucl .
phys .  * a661 * , 526 ( 1999 ) ; m.  post , s.  leupold , and u.  mosel , eprint : nucl - th/0008027 .
m.  asakawa et al . , phys .
rev .  * c46 * , 1159 ( 1992 ) .
g.  chanfray and p.  schuck , nucl .
* a555 * , 329 ( 1993 ) ; m.  herrmann , b.  friman , and w.  nrenberg , nucl.phys .
* a560 * ( 1993 ) 411 ; m.  urban et al .
* a641 * , 433 ( 1998 ) .
brown and m.  rho , phys .
66 * , 2720 ( 1991 ) ; g.q .
li , c.m .
ko , and g.e .
brown , nucl .
a606 * , 568 ( 1996 ) ; y.  kim et al . , phys .
* c62 * , 015202 ( 2000 ) .
titov , t.i .
gulamov , and b.  kmpfer , phys .
* d53 * , 3770 ( 1996 ) ; c.  song , phys .
* d53 * , 3962 ( 1996 ) ; h.  van  hees and j.  knoll , eprint : hep - ph/0007070 .
r.  rapp and c.  gale , phys .
rev .  * c60 * , 024903 .
lee , c.  song , and h.  yabu , phys .
lett . , 407 ( 1995 ) . c.  song , s.h .
lee and c.m .
ko , phys .
* c52 * , 476 ( 1995 ) .
voskresensky , jetp * 78 * , 793 ( 1994 ) .
shuryak , phys . rev . *
d42 * , 1764 ( 1990 ) ; nucl . phys . *
a533 * , 761 ( 1991 ) ; e.v .
shuryak and v.  thorsson , nucl .
* a536 * , 739 ( 1992 ) .