The luminosity distance – redshift relation for a wide class of generalized Randall-Sundrum type II brane-world models with Weyl fluid is compared to the presently available supernova data. We find that there is a class of spacially flat models with different amounts of matter Ωρ and Weyl fluid Ωd, which have a very similar fitting quality. The best-fit models are equally likely and can be regarded as extensions of the ΛCDM model, which is also included. We examine three models with different evolutionary history of the Weyl fluid, characterized by a parameter α = 0, 2 and 3. The first model describes a brane which had radiated energy into the bulk some time ago, but in recent times this energy exchange has ceased and only a dark radiation (α = 0) is left. In the other two models the Weyl-fluid describes a radiating brane throughout the cosmological evolution, up to our days. We find that the throughout of the fitting surface extends over a wider Ωd-range with increasing α, but the linear correlation of Ωd and Ωρ holds all over the examined Ωd range.
PACS codes: 11.25-w, 98.62.Py, 98.80.Es
Current observational data [1–3] suggest that the cosmological model of a Universe with only baryonic matter has to be modified. In the easiest way, the model can be reconciled with observations by the introduction of a cosmological constant Λ and of considerable amount of dark matter (ΛCDM model). Because the energy densities of both baryonic and dark matter decrease during cosmological evolution, the cosmological constant will dominate the late-time evolution. This process was first suggested for the explanation of Ia supernovae data, which suggest that our Universe has reached an accelerating phase. In a Λ-dominated universe, the luminosity distance increases faster with redshift than in the model without Λ , exactly as required by the supernova data.
Generally, the agreement with experiments can be achieved by introducing a dark energy component of the Universe, which replaces Λ. Such a dark energy in general does not clump. A recent analysis  shows that a dark energy model with varying dark energy density going through a transition from an accelerating to a decelerating phase at redshift 0.45 fits well the observational data. Based on observations, the dark energy equation of state w = p/ρ is within about -1 ± 0.1 .
It has been expected for some time that alternative gravitational theories, motivated by string/M-theory could replace dark matter and dark energy by geometric effects. The curved generalizations (see for example the review ) of the original Randall-Sundrum type II model  consist of a hypersurface with tension λ (the brane), representing our observable universe, embedded in a 5-dimensional space-time (the bulk). Gravitational dynamics on the brane is governed by an effective Einstein equation [9, 10]. The sources of gravity in the effective Einstein equation include terms due to the asymmetric embedding of the brane into the bulk , nonstandard model fields in the bulk, and even quantum corrections approximated as induced gravity effects [12–15].
The most relevant source term for early cosmology is a quadratic source term in the energy-momentum tensor . This term dominates over the linear term before the Big Bang Nucleosynthesis (BBN). In the simplest case of cosmological symmetries and suppression of the energy exchange between the brane and the bulk and whenever the bulk contains a static black hole, the Weyl curvature of the bulk generates a so-called Weyl fluid effect on the brane. In Fig 1 of our companion paper  (to be referred in what follows as paper I) we classify the different brane-world theories and their inter-relations. They are divided into two branches, one containing the original Randall-Sundrum type II model (BRANE1) and the other the flat DGP model (BRANE2). The model with Weyl fluid belongs to the BRANE 2 branch.
Supernova data were confronted with the induced gravity models [17–20]. When they are combined with the Sloan Digital Sky Survey (SDSS) baryonic peak, these seem to rule out the flat DGP models [17, 18]. However it was argued in  that the Cosmic Microwave Background (CMB) shift parameter can over-turn this conclusion. Structure formation and CMB were also considered in the DGP models in Ref. .
Most recently, the authors of  tested the accelerating phase of the universe's expansion with a comparison of the models and the supernova data. They have tested the ΛCDM model, the DGP model and three wCDM models with equations of state where w(a) (i) was constant with scale factor a, (ii) varied as w(a) = w0 + wa(1 - a) for redshifts probed by the supernovae but fixed at -1 for earlier epochs, and (iii) varied as w0 + Wa(1 - a) since the recombination. Their main conclusion is that all the five examined models explain equally well the acceleration, and none of them could be selected as a preferred model, based on the Ia-type supernova data.
The authors of Ref.  have compared the predictions of the flat BRANE1 and BRANE2 models to the Gold  and Supernova Legacy Survey (SNLS)  supernova data sets, incorporating the baryon acoustic peaks into the analysis. These brane-world models in certain parameter range (when their induced gravity parameter Ωl is small; the flat DGP models falling outside this range) are satisfied by both supernova data sets. The BRANE1 models fit better to the SNLS data, while the BRANE2 models fit better to the Gold data set. Since the analysis depends very weakly on the bulk cosmological constant , the value of was fixed at zero. With this modification, the BRANE1 model fits better to the SNLS data than the ΛCDM model and fits comparably well to the Gold data. The same conclusion holds for the BRANE2 model. In the analysis of  the dark radiation dimensionless parameter Ωd is switched off.
Using two recent supernova data sets, the CMB shift parameter, and the baryon oscillation peaks, the authors of Ref.  have found that the LDGP model (a subclass of the BRANE1 models with the effective energy density having a phantom-like behavior due to extra-dimensional effects, see Fig 1 of Paper I) fits the observations if it is very close to the ΛCDM model. The modification of the LDGP model with respect to the ΛCDM model appears in the form of a linear term in the Friedmann equation, H/rc, where H is the Hubble parameter and rc a crossover scale. This model includes a cosmological constant, possibly screened by the modified gravity, however the comparison with observations sets strong constraints on the screening.
The first comprehensive study of the generalized Randall-Sundrum type II (RS) brane-worlds tested against astronomical data was presented in Ref. . Agreement with earlier supernova data has been established in the presence of a cosmological constant. In this analysis the dark radiation from the bulk was switched off (Ωd = 0) and the energy-momentum squared term was kept. Under these assumptions, for flat spatial sections and matter parameter Ωρ = 0.3 the maximum likelihood method gave Ωλ = 0.004 ± 0.016 for the parameter characterizing the source term quadratic in the energy-momentum. This in turn implies a tiny value of the brane tension, which is disfavored by generic brane-world arguments. Moreover, much lower values for Ωλ emerge from both CMB and BBN.
In contrast to Refs.  and  the analysis of  keeps both Ωλ and Ωd, the latter obeying |Ωd| < 0.01. The best fit is obtained at Ωρ = 0.15, ΩΛ = 0.80, Ωλ = 0.026 and Ωd = 0.008. With the high value of the brane tension set by either (a) the value of the 4-dimensional Planck constant and sub-millimeter tests  on possible deviations from Newton's law (in units c = 1 = ħ these give TeV4 see [29, 7]), (b) astrophysical considerations MeV4 or (c) BBN constraints MeV4 , the quadratic source term barely counts at late-times in the cosmological evolution.
Given the high limits for the values of λ, in any realistic model Ωλ can be safely ignored. This is a crucial difference of our forthcoming analysis as compared to the one presented in Refs.  and , where the corresponding cosmological parameter Ωλ was kept.
The next question is whether the source term arising from the Weyl curvature of the bulk may be kept, in other worlds, whether Ωd ≠ 0. The Weyl curvature of the bulk gives an energy density ρd = 6m/κ2a4, where κ2 = 8πG is the gravitational coupling constant. In Ref.  it was shown that the BBN limits constrained the dark radiation component as -1.23 ≤ ρd (zBBN)/ργ (zBBN) = 0.11. Combining this with CMB constraints reduces this range to -0.41 ≤ ρd (zBBN)/ργ (zBBN) ≤ 0.105. Here ργ is the energy density of the background photons. Another constraint for the value of the dark radiation at BBN was derived in  as -1 <ρd (zBBN)/ρν (zBBN) < 0.5 where ρν is the energy density contributed by a single, two-component massless neutrino. This constraint was derived for high values of the 5-dimensional Plank mass.
In the simplest case the Weyl source term evolves as a radiation, thus its present value is obviously tiny. This is the reason why all mentioned references  and  comparing RS brane-worlds with observations disregard dark radiation. But is this a necessary assumption? Formulating the question the other way around: if we include even a small component of dark radiation into the late-time universe model we face a serious problem. Due to the fact that the energy density of dark radiation decreases as a-4 (compared to that of matter which is a-3), even an amount of dark radiation of the same order as the amount of baryonic matter nowadays implies dark radiation dominance in the past, for example during structure formation. This conclusion is contradicted by numerical simulations, which favorize cold dark matter as the dominant component of the Universe during structure formation .
However we can generalize the validity of the model by lifting the requirement of a constant mass m in the dark radiation energy density. A constant m implies a static Schwarzschild-anti de Sitter bulk and no energy exchange between the brane and the bulk. Therefore dark radiation is a manifestation of an equilibrium configuration with a static bulk, and it may be well possible that such a situation is reached only at the latest stages of the evolution of the brane-world Universe. Whenever m depends on a certain, non-zero power of a, the evolution of the energy density of the Weyl source term evolves in a non-standard way, allowing to escape from the argument of a small Weyl fluid left nowadays.
We propose here the LWRS (Lambda-Weyl fluid-Randall-Sundrum) model, a specific RS model with i) cosmological constant, ii) the brane radiating away energy during various stages of the cosmological evolution, characterized by the index a and iii) a Weyl fluid depending on the actiul value of a during the latest stage of cosmological evolution, which can be tested by supernova observations. For the inclusion of the LWRS model in the classification of brane-world models, see Fig 1 of paper I.
The LWRS model takes into account the possibility of an energy exchange between the brane and the bulk. This idea is not new. Indeed, it was already proposed that during an inflationary phase on the brane radiation is emitted and black holes thermally nucleate in the bulk . Later on, but still in the high energy regime, the brane radiates such that the mass function of the bulk black hole increases with a4 . This means that the Weyl source term becomes a constant in this era. The brane continues to radiate away energy during structure formation , a process leading to a bulk black hole mass function m ∝ aα, with 1 ≤ α ≤ 4. (Other models with the brane radiating energy into the bulk are also known .)
For α = 0 the Weyl fluid is known as dark radiation, for α = 2, 3 it gives the correct growth factor during structure formation. For α = 1, 4 it is indistinguishable from dark matter and a cosmological constant, respectively. Therefore pure dark radiation can emerge only in the low-z limit, while at earlier times a dynamic bulk – brane interaction governed by energy exchange should be present.
2. Confronting the models with the selected supernova data
As type Ia supernovae result from the explosion of white dwarf stars with identical mass, they show remarkable similarities. By employing well established calibration methods, one can calculate the maximal luminosity of the object (in the reference system of the explosion). This is done by analyzing the time-dependent variation of the emitted luminosity and the spectrum, a method known as the Multi-Color Light Curve analysis [39, 23]. In this process the observed parameters, the shape of the light curve and the spectral distribution of the emission have to be converted into the reference system of the host galaxy. For distant supernovae this translates to take into account the time dilation and the so-called K-correction . While these methods depend on z, they are independent on the specific cosmological model. After performing these corrections, we have well-calibrated maximal luminosities for the supernovae of type Ia and in consequence they are considered as standard candles.
In 2003 a list of dL-z data pairs were published for 230 supernovae of type Ia , and 60 of them had low absoprtion (i.e. AV.1) and cosmological redshift (z > 0.01). To give result which can easily be compared to earlier works, we also involve this selected low-absorption data set for the most of the examinations. The basic supernova data we use here is the improved Gold set , which was released in 2006.
We confront with supernova observations several models from paper I. In Fig 1 we represent graphically on both logarithmic and linear scales their luminosity distance -redshift relations up to z = 2.5. The plots are for k = 0 and Ωρ = 0.27 (according to the combined analysis of the SDSS and WMAP 1-year data in Ref ). In particular, the luminosity distance – redshift relation is shown for the following models:
The LWRS model (the perturbative solution given by Eqs. (56), (60)–(61), (63), (65) and (67) of paper I, with Ωλ= 0 and for α = 0) for the two values of the late-time dark radiation Ωd= -0.05 and Ωd= 0.05 (the curves 1 and 3, respectively). The latter models contain a brane which radiates energy at early times (for Ωd> 0) and during structure formation, such that a bulk black hole is formed and its mass increases continuously. As this process slows down, the Weyl curvature of the bulk induces the late-time dark radiation on the brane.
The ΛCDM model, given by Eqs. (60)–(61) of paper I (curve 2).
The solutions with brane tension λ = 2Λ/κ2and no dark radiation (given by Eq. (52) of paper I) for both admissible values for this model, at ΩΛ = 0.704 (curve 4) and ΩΛ = 0.026 (curve 6). The former is similar to the class of models discussed in .
The late-time universe ΩΛ = 0 limit of the RS model with Randall-Sundrum fine-tuning, containing a huge amount of dark radiation Ωd= 0.73, given by Eq. (44) of paper I (curve 5).
In Fig. 1 we plot these models in a comparison to low-absorption supernova data from Ref.  (red triangles) together with the Gold set  (black dots). The error bars are indicated in the respective colors. The diagrams with linear scale are more instructive, as they emphasize the difference among the predictions of the chosen models and how they fit data, while the logarithmic scale better disseminate between the low z points.
The models represented by the curves 1, 3 and 4 by eye seem to compare as well with the supernova observations as the ΛCDM model (curve 2). By contrast, the models represented by the curves 5 and 6 seem to be not supported by observations. The model with no cosmological constant and significant dark radiation Ωd = 0.73, Ωλ = 0 (curve 6) and the model with Λ = κ2λ/2 and ΩΛ = 0.025 (curve 5) are significantly inconsistent with the observations, as they give χ2 = 213 and 395, respectively‡. All other models shown on Fig 1 are comparable with the supernova observations, as it was expected by a simple glance.
The χ2 = 50 value found for the Λ = κ2λ/2-model with ΩΛ = 0.74 (curve 4) is slightly better than χ2 found the ΛCDM model. However, as mentioned earlier, the tiny brane tension λ = 38.375 × 10-60 TeV4, several order of magnitudes lower than all existing lower limits rules out this model as well.
The best fitting models are the models with brane cosmological constant; a high value of the brane tension (leading to ΩΛ ≈ 0) and a small contribution of dark radiation, Ωd = ±0.05 (the curves 1 and 3). For Ωd = -0.05 we find χ2 = 65, which is still acceptable. For Ωd = 0.05 we get χ2 = 49.
Values of Ωd between these limits are also admissible. It is likely that by increasing Ωd towards higher positive values, χ2 remains compatible, however the accuracy of the perturbative solution is deteriorated with increasing Ωd, therefore higher orders in the expansion would be necessary to take into account.
3. The Gold2006 set of supernovae
More recently, Riess et al.  have published a new set of 182 gold supernovae, including new HST observations and recalibrations of the previous measurements. It is an interesting question how this recalibration influenced the above conclusions for the well-fitting models with dark radiation.
We applied the same tests to the Gold2006 data set as described in the previous section. First we assumed that Ωρ = 0.27, as before, cf. Ref . In this case the critical value of χ2 is 197 at 80% level and 209 at 90% confidence level. Then the models represented by the curves 1–4 of Fig 2 behave as follows. The model with a small amount of negative dark radiation is disfavored at 80% confidence (χ2 = 204). The models with λ = 2Λ/k2 and ΩΛ = 0.704 are ruled out at 90% confidence level, too (as χ2 = 221). As expected from the previous analysis, the ΛCDM model (χ2 = 192) and the LWRS model with Ωd = 0.05 (giving χ2 = 194) compete closely. We also remark that varying Ωd between -0.03 and 0.07, the χ2 remains under the critical value.
For gaining a deeper insight we have then calculated the predictions of the models between Ωd = -0.10 ÷ 0.10 with a stepsize of 0.01 in Ωd, with Ωρ allowed to freely vary in the domain 0.15 ÷ 0.35 and z in the range 0 ÷ 3. Then we looked for the best fit of the Gold2006 set in the Ωd - Ωρ space. This is represented on Fig 3. The global minimum of the surface is at Ωd = 0.040, Ωρ = 0.225 (χ2 = 190.52), which suggests an interesting opportunity for a Universe with less baryonic density and with dark radiation, compatible with the Gold2006 supernova data. The 1-σ confidence interval is centered about this value. The ΛCDM model (where Ωd is exactly 0) has the local minimum of Ωρ = 0.275 (χ2 = 195.8), but this is outside the 1-σ confidence interval.
Similar conclusions emerge from the plot in the ΩΛ - Ωρ plane, Fig 4. Here the global minimum of the surface is at ΩΛ = 0.735, Ωρ = 0.225. The local minimum of the CDM model is at ΩΛ = 0.725.
We note that there is a forbidden parameter range in both planes Ωd - Ωρ and ΩΛ - Ωρ, represented by white regions on Figs 3 and 4. This is because the Friedmann equation for these brane-world models
combined with ΩΛ + Ωρ + Ωd = 1 gives the constraints
in the Ωd - Ωρ plane and
in the ΩΛ - Ωρ plane.
The forbidden region increases in both cases with z. If we would like to extend the limits to z → ∝, we obtain the limiting curves in the Ωd - Ωρ plane and in the ΩΛ - Ωρ plane. However the LWRS model being valid only for low values of z, we represent on the graphs only the forbidden range for z = 3.
4. The compatibility of the LWRS model Ωd = 0.04 and α = 0 with cosmological evolution
The energy density of dark radiation decreases too fast during cosmological evolution to result in a considerable amount nowadays. We discuss this problem and its possible remedy in detail here. First we comment on the problem, then we show how an energy exchange between the brane and the bulk can leave a considerable amount of dark radiation.
The constraint derived in  for the energy density of the dark radiation:
where is the energy density of the background photons at the beginning of BBN. The coefficient
contains [11, 43] the effective number g* of relativistic degrees of freedom, which depends on the temperature. According to g* = 10.75 at the beginning of BBN, when TBBN = 1.16 × 1010 K. Thus ργ (zBBN) = 7.37 × 1025 J m-3 emerges, giving the constraint
Note, that the domain of allowable negative values is larger than the one for positive values.
As for today the background photons have cooled to T0 = 2.725 K and for such low temperatures g* = 3.36 [43, 44] their energy density ργ = ργ (z = 0) is
With the value Km s-1 Mpc-1 of the Hubble constant , cf. Eq. (23) of paper I the present day cosmological parameters ρ and Ω (both for background and dark radiation) relate as
Thus the present value of Ωγ is
which is quite negligible. If the Weyl source term were to evolve as radiation, its value would be even smaller, cf. Eq. (4). Indeed Eqs. (6) and (8) imply
|Ωd| is of the same order of magnitude or smaller as Ωγ.
However if the brane is radiating during structure formation, the mass parameter m becomes a function of the scale factor m ∝ aα, with 1 ≤ α ≤ 4 . Then the energy density scales as a4-α.
Now let us suppose that the brane is in an equilibrium (non-radiating) configuration with α = 0 in the domain 0 ≤ z ≤ z1. In a preceding era z1 <z ≤ z* the brane radiates such that α ≠ 0, finally right after the beginning of BBN, at z* <z ≤ zBBN there is equilibrium once more (α = 0). Here zBBN = (TBBN/T0) - 1 = 4.26 × 109. According to this evolution
Inserting this in Eq. (6) and employing Eq. (8) we obtain:
In the particular case α = 0 we recover the constraint (10) set on pure dark radiation. However for any α > 0 we get
Let us specify this result for the best fit value Ωd = 0.04. Depending on α we obtain the following numerical relations between the redshifts characterizing the switching on and off of the radiation leaving the brane:
It is evident that the value of z* increases with z1 (this dependence becoming an approximate scaling for higher values of z1) and decreases with α.
The lower limit in the LWRS model is z1 = 3. Then
For the higher values of α the duration of the radiative brane regime necessary to produce a high value of Ωd today is quite short.
5. LWRS models with α = 2, 3 confronted with supernova data
In the absence of a known mechanism for changing a, we examine here the cases when α = 2 and α = 3 hold throughout the cosmological evolution, up to nowadays. For this we confront these models with the Gold 2006 set exactly as described before. To preserve the validity of the perturbative solution, the range of Ωd was selected to be -0.1–0.1, and we probed the range 0.15–0.35 of Ωρ. The assumption for flatness was kept, too.
The results are qualitatively similar to the α = 0 case. The remarkable difference is that the peak of the minimum turned into a "trough", which lies aslope in the Ωρ - Ωd space. This means that instead of a district solution, a complete model family exists in both cases, which can equally well explain the supernova data. The Ωρ dependence of Ωd is less in the α = 2 model as compared to α = 0, and is very small if α = 3. The Ωd = 0 case is the ΛCDM model where these model intersect. The steeper slope of the minimum trough thus allows a much lower range of Ωρ in the α = 2, and especially in the α = 3 models, with a value close to 0.3. On the other hand, the range of Ωd gets more and more wide with increasing α, which results in the conclusion that the presence of Ωd is mathematically plausible, and they have to be accounted for in RS cosmology.
Due to the higher slopes of the 1-σ and 2-σ contours in these α = 2, 3 models, Ωρ is much less affected by the Weyl fluid, while Ωd can have various values in the detriment of ΩΛ. Therefore the Weyl fluid can explain some of the dark energy.
The luminosity distance given in paper I as function of redshift in terms of elementary functions and elliptical integrals of first and second type for various brane-world models with Weyl fluid was confronted with the available supernova data sets, including the Gold2006 data . The tested models were:
The models with Randall-Sundrum fine-tuning, discussed in section 4 of paper I, with a considerable amount of dark radiation as a bulk effect, and a high value of the brane tension.
The two models discussed in subsection 5.1 of paper I, which obey Λ = κ2λ/2, have no dark radiation and were integrable in terms of elementary functions.
The LWRS models (subsection 5.2 of paper I), with a brane cosmological constant, for which the luminosity distance could be given analytically as function of redshift to first order accuracy in the dark radiation. (Due to its smallness, the source term Ωλ quadratic in the energy density was suppressed in the perturbative models of paper I.)
The brane-world models (A) although interesting for historical reasons, do not comply with observations. Even if we introduce an extremely high amount of dark radiation Ωd = 0.73, tentatively replacing the cosmological constant in the energy balance ΩΛ + Ωρ + Ωd + Ωλ = 1, these models are quickly outruled by supernova data (curve 5 of Fig 1). Dark radiation is not capable to replace the cosmological constant in producing a late-time acceleration, since it scales as usual radiation. The more we go back in the past, the higher becomes its domination over matter. Therefore a cosmological constant or dark energy is still needed in the generalized Randall-Sundrum type II models.
Our analysis has also dismissed immediately the model (B) with ΩΛ = 0.026. Surprisingly, the other toy model (B) with ΩΛ = 0.704 was in good agreement with the Gold2006 data, but ruled out by its low value of the brane tension, similarly as the models discussed in Ref. . A low brane tension is in disagreement with various upper limits set by cosmological and astrophysical tests.
The perturbative approach of subsection 5.2 of paper I can be considered valid for a Weyl fluid with -0.1 < Ωd < 0.1. In this range the LWRS brane-world models (C) were confronted with supernova data and for α = 0 the dark radiation with significant negative energy density ruled out. The fact that a positive dark radiation (corresponding to a bulk black hole rather than to a bulk naked singularity) is favoured by the presently available best supernova data is in accordance with the early behavior of the RS model with late-time dark radiation, where the brane radiating away energy in early times leads to a black hole, which can further grow during structure formation.
The remaining LWRS brane-world models with α = 0 and Ωd between -0.03 and 0.07 (and Ωρ changed accordingly) turned out to be excellent candidates for describing our universe, as they show remarkable agreement with the Gold2006 supernova data sets. If Ωρ is allowed to vary in the range (0.15, 0.35), the preferred values are Ωd = 0.040, Ωρ = 0.225, ΩΛ = 0.735.
The preferred cosmological parameters determined by comparing the LWRS model with α = 0 with supernova data alone are in perfect accordance with the WMAP 3-year data. Indeed according to Ref.  and from which emerge. The value of Ωρ determined by comparing the LWRS model with the supernova data alone is well in the middle of the domain allowed by the WMAP 3 year data.
We have then proved that the preferred value of Ωd = 0.04 is compatible with the known history of the Universe if the brane radiates away energy into the bulk during a relatively short period of the cosmological evolution. Such a process occurring between z = 24 and z = 3 could increase the amount of dark energy today with a factor of 103 as compared to the non-radiating brane, exactly as required by the LWRS model with α = 0.
The LWRS models with α = 1 (α = 4) are identical with the ΛCDM model with the only difference that some fraction of the dark matter (of the cosmological constant) has geometric origin.
Finally, the LWRS models with α = 2 and α = 3 do not present a sharp minimum, but rather an elongated trought shape in the parameter space, with the slope increasing with the value of α and Ωρ ≈ 0.3. This means that in this class of models a wide range of values for Ωd (with a slight preference for negative values) and corresponding values for ΩΛ are fitting to the supernova data.
We must note that the reliability of these values is somehow deteriorated by the relatively small number of high-z supernova and by the inherent difficulties in the calibration of the available data. An obvious source of error is that data from the Gold2006 set is a combination of measurements taken on different instruments  and in fact it has been already signaled that the Gold2006 data set is not statistically homogeneous .
The conclusion of this paper is somewhat similar to that of Ref. : the presently available supernova data are not enough to discern among several cosmological models. However the difference between the predictions of the acceptable models of our analysis (the ΛCDM model, the LWRS brane-world with α = 0 and Ωd = 0.04 and the models with Weyl fluid and α = 2, 3) are increasing with z. One may reasonably hope that the very far (z > 2) supernovae, which will be discovered for sure in the following decade, will improve their comparison.
We thank Botond Nagy for interactions in the early stages of this research. This work was supported by OTKA grants no. 46939, 44665, 42509 and 69036. LÁG and GyMSz were further supported by the János Bolyai Grant of the Hungarian Academy of Sciences and GyMSz by the Magyary Zoltán Higher Educational Public Foundation.
Department of Experimental Physics, University of Szeged
Department of Theoretical Physics, University of Szeged
Department of Applied Science, London South Bank University
Doroshkevich A, Tucker DL, Allam S, Way MJ: Large scale structure in the SDSS galaxy survey.Astron Astrophys 2004, 418:7.MATHView ArticleADS
Tegmark M, Strauss MA, Blanton MR, et al.: Cosmological parameters from SDSS and WMAP.Phys Rev D 2004, 69:103501.View ArticleADS
Spergel DN, Bean R, Doré O, et al.: Wilkinson Microwave Anisotropy Probe (WMAP) Three Year Results: Implications for Cosmology. arXiv:astro-ph/0603449 2006.
Sahni V, Starobinski A: The Case for a Positive Cosmological Lambda-term.Int J Mod Phys D 2000, 9:373.ADS
Fay S, Tavakol R: A model-independent dark energy reconstraction schem using the geometrical form of the luminosity-distance relation.Phys Rev D 2006, 74:083513.View ArticleADS
Krauss LM, Jones-Smith K, Huterer D: Dark Energy: A Cosmological Constant, and Type Ia Supernovae. arXiv:astro-ph/0701692 2007.
Keresztes Z, Gergely LÁ, Nagy B, Szabó Gy M: The luminosity-redshift relation in brane-worlds: I.Analytical results 2007. In Eq (15) ?? should read r
Fairbairn M, Goobar A: Supernova limits on brane world cosmology. arXiv:astro-ph/0511029 2005.
Alam U, Sahni V: Confronting Braneworld Cosmology with Supernova data and Baryon Oscillations.Phys Rev D 2006, 73:084024.View ArticleADS
Alam U, Sahni V: Supernova Constraints on Braneworld Dark Energy. arXiv:astro-ph/0209443 2002.
Maartens R, Majerotto E: Observational constraints on self-accelerating cosmology.Phys Rev D 2006, 74:023004.View ArticleMathSciNetADS
Koyama K, Maartens R: Structure formation in the DGP cosmological model.J Cosmol Astroparticle Phys JCAP 2006,06(01):016.View ArticleMathSciNet
Barger V, Gao Y, Marfatia D: Accelerating cosmologies tested by distance measures.Phys Let B 2007, 648:127.View ArticleADS
Ries AG, Strolger L-G, Tonry J, et al.: Type Ia Supernova Discoveries atz> 1 from the Hubble Space Telescope: Evidence for Past Deceleration and Constraints on Dark Energy Evolution.Asrophys J 2005, 607:665.View ArticleADS
Astier P, Guy J, Regnault JN, et al.: The Supernova Legacy Survey: Measurement of ΩM, ΩΛand w from the First Year Data Set.Astron Astrophys 2006, 447:31.View ArticleADS
Lazkoz R, Maartens R, Majerotto E: Observational constraints on phantom-like braneworld cosmologies.Phys Rev D 2006, 74:083510.View ArticleMathSciNetADS
Dąbrowski MP, Godlowski W, Szydlowski M: Brane universes tested against astronomical data.Int J Mod Phys D 2004, 13:1669.MATHView ArticleADS
Langlois D, Sorbo L, Rodríguez-Martínez M: Cosmology of a brane radiating gravitons into the extra dimension.Phys Rev Lett 2002, 89:171301.View ArticleMathSciNetADS
Pal S: Structure formation on the brane: A mimicry.Phys Rev D 2006, 74:024005.View ArticleADS
Hebecker A, March-Russell J: Randall-Sundrum II cosmology, AdS/CFT, and the bulk black hole.Nuclear Physics B 2001, 608:375. Leeper E, Maartens R and Sopuerta C, Dynamics of radiating braneworlds 2004, Class Quant Grav. 21 1125. Gergely L Á, Leeper E and Maartens R, Asymmetric radiating brane-world 2004 Phys. Rev D70 104025. Jennings D and Vernon I R, Graviton emmission into non-Z2 symmetric brane world spacetimes 2005 J. Cosmol. Astroparticle Phys. JCAP 05 (07) 011. Langlois D, Is our Universe brany? 2006 Progress of Theoretical Physics Supplement163 258.MATHView ArticleMathSciNetADS
Riess AG, Filippenko AV, Challis P, et al.: Observational evidence from Supernovae for an Accelerating Universe and a Cosmological Constant.Astrophys J 1998, 116:1009.
Schmidt B: 2005 in: The new cosmology, Proceedings of the 16th Int. In Summer School. Singapore World Scientific; 2003.
Tonry JL, Schmidt BP, Barris B, et al.: Cosmological Results from High-z Supernovae.Astrophys J 2003, 594:1.View ArticleADS
Riess AG, Strolger L-G, Casertano S, et al.: New Hubble Space Telescope Discoveries of Type Ia Supernovae atz> 1: Narrowing Constraints on the Early Behavior of Dark Energy.to appear in Astrophys J 2007., 656:e-printatarXiv.org