Open Access

The luminosity-redshift relation in brane-worlds: I. Analytical results

  • Zoltán Keresztes1,
  • László Á Gergely1,
  • Botond Nagy1 and
  • Gyula M Szabó1
PMC Physics A20071:4

DOI: 10.1186/1754-0410-1-4

Received: 26 June 2007

Accepted: 2 October 2007

Published: 2 October 2007

Abstract

The luminosity distance – redshift relation is analytically given for generalized Randall-Sundrum type II brane-world models containing Weyl fluid either as dark radiation or as a radiation field from the brane. The derived expressions contain both elementary functions and elliptic integrals of the first and second kind. First we derive the relation for models with the Randall-Sundrum fine-tuning. Then we generalize the method for models with cosmological constant. The most interesting models contain small amounts of Weyl fluid, expected to be in good accordance with supernova data. The derived analytical results are suitable for testing brane-world models with Weyl fluid when future supernova data at higher redshifts will be available.

PACS Codes: 98.62.Py, 98.80.Jk, 11.25.-w

1 Introduction

At present the Universe is considered a general relativistic Friedmann space-time with flat spatial sections, containing more than 70% dark energy and at about 25% of dark matter. Dark energy could be simply a cosmological constant Λ, or quintessence or something entirely different. There is no widely accepted explanations for the nature of any of the dark matter or dark energy (even the existence of the cosmological constant remains unexplained).

An alternative to introducing dark matter would be to modify the law of gravitation, like in MOND [13] and its relativistic generalization [4, 5]. These theories are compatible with the Large scale structure of the Universe [68]. However in spite of the successes, certain problems were signaled on smaller scales [913].

Quite remarkably, supernova data, which in the traditional interpretation yield to the existence of dark energy, can be explained by certain f(R) [14, 15] or inverse curvature gravity models [16]. However the parameter range, in which the latter is in goood agrement with the supernova data, also presents stability problems [17, 18].

Modifications of the gravitational interaction could also occur by enriching the space-time with extra dimensions. Originally pioneered by Kaluza and Klein, such theories contained compact extra dimensions. The so-called brane-world models, motivated by string/M-theory, containing our observable 4-dimensional universe (the brane) as a hypersurface, were introduced in [1921] and [22], the latter model allowing for a non-compact extra dimension.

The curved generalizations of the model presented in [22] have evolved into a 5-dimensional alternative to general relativity, in which gravity has more degrees of freedom. In contrast with standard model fields, these evolve in the whole 5-dimensional bulk. In this generalized Randall-Sundrum type II (RS) theory, the brane has a tension λ and gravitational dynamics is governed by the 5-dimensional Einstein equation. Its projections to our observable 4-dimensional universe (the brane) are the twice contracted Gauss equation, the Codazzi equation and an effective Einstein equation, the latter being obtained by employing the junction conditions across the brane [23]. The effective Einstein equation (for the case of symmetric embedding and no other contribution to the bulk-energy-momentum than a bulk cosmological constant) was first given in a covariant form in [24]. Supplementing this by the pull-back to the brane of the bulk energy momentum tensor Π ˜ a b MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuqHGoaugaacamaaBaaaleaacqWGHbqycqWGIbGyaeqaaaaa@30FD@ , which is
P a b = 2 κ ˜ 2 3 ( g a c g b d Π ˜ c d ) T F MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqqGqbaudaWgaaWcbaGaemyyaeMaemOyaigabeaakiabg2da9maalaaabaGaeGOmaidcciGaf8NUdSMbaGaadaahaaWcbeqaaiabikdaYaaaaOqaaiabiodaZaaadaqadaqaaiabdEgaNnaaDaaaleaacqWGHbqyaeaacqWGJbWyaaGccqWGNbWzdaqhaaWcbaGaemOyaigabaGaemizaqgaaOGafuiOdaLbaGaadaWgaaWcbaGaem4yamMaemizaqgabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaemivaqLaemOrayeaaaaa@4741@
(1)
(with κ ˜ 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacamaaCaaaleqabaGaeGOmaidaaaaa@2F93@ the bulk coupling constant and g ab the induced metric on the brane) the effective Einstein equation reads [23]:
G a b = Λ g a b + κ 2 T a b + κ ˜ 4 S a b a b + P a b . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGhbWrdaWgaaWcbaGaemyyaeMaemOyaigabeaakiabg2da9iabgkHiTiabfU5amjabdEgaNnaaBaaaleaacqWGHbqycqWGIbGyaeqaaOGaey4kaSccciGae8NUdS2aaWbaaSqabeaacqaIYaGmaaGccqWGubavdaWgaaWcbaGaemyyaeMaemOyaigabeaakiabgUcaRiqb=P7aRzaaiaWaaWbaaSqabeaacqaI0aanaaGccqWGtbWudaWgaaWcbaGaemyyaeMaemOyaigabeaakiabgkHiTmrtHrhAL1wy0L2yHvtyaeHbnfgDOvwBHrxAJfwnaGabaiab+btifnaaBaaaleaacqWGHbqycqWGIbGyaeqaaOGaey4kaSIaeeiuaa1aaSbaaSqaaiabdggaHjabdkgaIbqabaGccqGGUaGlaaa@5BF6@
(2)
Here κ2 is the brane coupling constant, related to the bulk coupling constant and the brane tension λ as 6κ2 = κ ˜ 4 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWF6oWAgaacamaaCaaaleqabaGaeGinaqdaaaaa@2F97@ λ, and
Λ = κ 2 2 λ κ ˜ 2 2 n c n d Π ˜ c d MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqqHBoatcqGH9aqpdaWcaaqaaGGaciab=P7aRnaaCaaaleqabaGaeGOmaidaaaGcbaGaeGOmaidaaiab=T7aSjabgkHiTmaalaaabaGaf8NUdSMbaGaadaahaaWcbeqaaiabikdaYaaaaOqaaiabikdaYaaacqWGUbGBdaahaaWcbeqaaiabdogaJbaakiabd6gaUnaaCaaaleqabaGaemizaqgaaOGafuiOdaLbaGaadaWgaaWcbaGaem4yamMaemizaqgabeaaaaa@43BF@
(3)
represents a cosmological "constant" which possibly varies due to the normal projection of the bulk energy-momentum tensor (this includes the contribution – Λ ˜ g a b MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuqHBoatgaacaiabdEgaNnaaBaaaleaacqWGHbqycqWGIbGyaeqaaaaa@324B@ due to the bulk cosmological constant Λ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuqHBoatgaacaaaa@2E30@ ). The source term S ab is quadratic in the brane energy-momentum tensor T ab :
S a b = 1 4 [ T a c T b c + 1 3 T T a b g a b 2 ( T c d T c d + 1 3 T 2 ) ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGtbWudaWgaaWcbaGaemyyaeMaemOyaigabeaakiabg2da9maalaaabaGaeGymaedabaGaeGinaqdaamaadmaabaGaeyOeI0Iaemivaq1aaSbaaSqaaiabdggaHjabdogaJbqabaGccqWGubavdaqhaaWcbaGaemOyaigabaGaem4yamgaaOGaey4kaSYaaSaaaeaacqaIXaqmaeaacqaIZaWmaaGaemivaqLaemivaq1aaSbaaSqaaiabdggaHjabdkgaIbqabaGccqGHsisldaWcaaqaaiabdEgaNnaaBaaaleaacqWGHbqycqWGIbGyaeqaaaGcbaGaeGOmaidaamaabmaabaGaeyOeI0Iaemivaq1aaSbaaSqaaiabdogaJjabdsgaKbqabaGccqWGubavdaahaaWcbeqaaiabdogaJjabdsgaKbaakiabgUcaRmaalaaabaGaeGymaedabaGaeG4mamdaaiabdsfaunaaCaaaleqabaGaeGOmaidaaaGccaGLOaGaayzkaaaacaGLBbGaayzxaaGaeiilaWcaaa@5D38@
(4)
and a b MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFWesrdaWgaaWcbaGaemyyaeMaemOyaigabeaaaaa@3A48@ is the electric part of the bulk Weyl tensor C ˜ a b c d MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGdbWqgaacamaaBaaaleaacqWGHbqycqWGIbGycqWGJbWycqWGKbazaeqaaaaa@332E@ , given as
a c = C ˜ a b c d n b n d . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFWesrdaWgaaWcbaGaemyyaeMaem4yamgabeaakiabg2da9iqbdoeadzaaiaWaaSbaaSqaaiabdggaHjabdkgaIjabdogaJjabdsgaKbqabaGccqWGUbGBdaahaaWcbeqaaiabdkgaIbaakiabd6gaUnaaCaaaleqabaGaemizaqgaaOGaeiOla4caaa@48A0@
(5)

In a cosmological context and suppressing any energy exchange between the brane and the bulk, this latter term generates the so-called dark radiation. Otherwise it can be called a Weyl fluid.

A review of many aspects related to the theories described by the effective Einstein equation (2) can be found in [25]. Both early cosmology [26] and gravitational collapse [2732] are essentially modified in these theories. There is also possible to replace dark matter with geometric effects in the interpretation of galactic rotation curves, weak lensing and galaxy cluster dynamics [3336].

The possible modifications of gravitational dynamics are even more versatile in the so-called induced gravity models. These can be regarded as brane-world models enhanced with the first quantum-correction arising from the interaction of the brane matter with bulk gravity. The induced gravity correction couples to the 5-dimensional Einstein-Hilbert action with the coupling constant γ κ ˜ 2 / κ 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFZoWzcuWF6oWAgaacamaaCaaaleqabaGaeGOmaidaaOGaei4la8Iae8NUdS2aaWbaaSqabeaacqaIYaGmaaaaaa@34F1@ . The simplest of such models, the DGP model was introduced in [37]. This model however suffers from linear instabilities (ghost modes in the perturbations), as shown for de Sitter branes [3840]. The ghost modes withstand even the introduction of a second brane [41]. Generalizations of the DGP model are discussed covariantly in [42] and [43] when the embedding is symmetric, and in [44] when it is asymmetric. In these models the role of the effective Einstein equation (2) is taken by a more complicated equation (see for example Eq. (29) of [44]), which contains the square of the Einstein tensor G ab . This implies that in certain sense the degree of nonlinearity of the theory is squared. In a cosmological setup the square root of this equation can be taken, leading to a set of modified Friedmann and Raychaudhuri equations, which however contain a sign ambiguity ε = ±1 due to the involved square root. These are called the BRANE1 [DGP(-)] branch for ε = -1 and BRANE2 [DGP(+)] for ε = 1 in the terminology of [42] [or [45], respectively]. Both the original Randall-Sundrum type II model and the DGP model are contained as special subcases. Notably, the BRANE2 branch contains cosmological models which self-accelerate at late-times. We give in Fig 1 a diagram containing a classification of these theories and how they emerge as different limits from each other.
https://static-content.springer.com/image/art%3A10.1186%2F1754-0410-1-4/MediaObjects/13066_2007_Article_4_Fig1_HTML.jpg
Figure 1

A diagram presenting various brane-world models and their inter-relations. LWRS is the generalized Randall-Sundrum model with cosmological constant and a Weyl fluid reflecting a brane radiating into the bulk during nowadays or at least until recent cosmological times.

In this paper we discuss analytically the luminosity distance – redshift relation in various generalized Randall-Sundrum type II brane-world models described by Eq. (2). Our analytical approach can enhance the confrontation of these models with current and most notably, with future supernova observations. We note that recently analytical results have been given in Ref. [46] for a wide class of phantom Friedmann cosmologies too, in terms of elementary and Weierstrass elliptic functions.

In section 2 we review the notion of luminosity distance, its relation with the redshift and how these can be measured independently. This section was included mainly for didactical purposes.

In section 3 we review the modification of this relation in the Randall-Sundrum type II brane-world scenario. These include the introduction of the parameters Ω λ and Ω d which can be traced back to the source terms S ab and a b MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFWesrdaWgaaWcbaGaemyyaeMaemOyaigabeaaaaa@3A48@ of the modified Einstein equation (2). The other cosmological parameters are Ω ρ , representing (baryonic and dark) matter and ΩΛ. We do not include bulk sources in the analysis, with the notable exception of a bulk cosmological constant.

Section 4 contains the derivation of the analytic expression for the luminosity distance – redshift relation for the brane-worlds which are closest to the original Randall-Sundrum scenario [22], thus with no cosmological constant (Randall-Sundrum fine-tuning). The generic expression (35) of the luminosity distance derived here is given in terms of elementary functions and elliptic integrals of the first and second kind. From this most generic case we take the subsequent limits: Ω d = 0 (subsection 4.2), Ω λ = 0 (subsection 4.3); and both Ω d = Ω λ = 0, this being the general relativistic Einstein-de Sitter case (subsection 4.4).

Such models however could not allow for late-time acceleration, therefore in section 5 we discuss the luminosity distance – redshift relation for brane-worlds with Λ. First we present in subsection 5.1 a class of models, for which the luminosity distance can be given in terms of elementary functions alone. These models are characterized by an extremely low value of the brane tension, thus are in conflict with various constraints on brane-world models.

Next, in subsection 5.2 we discuss brane-worlds for which the brane-characteristic contributions Ω λ and Ω d represent small perturbations. This is a good assumption as observational evidences suggest that general relativity is a sufficiently accurate theory of the universe, and as such the deviations from it could not be very high, at least at late-times. We give analytical expressions in terms of both elementary functions and elliptic integrals of the first and second kind for the luminosity distance, to first order accuracy in the chosen small parameters of the model. Some of the most lengthy computations needed in order to achieve the result are presented in the Appendix.

Section 6 contains the concluding remarks.

Throughout the paper c = 1 was employed.

2 The luminosity-redshift relation

The Friedmann-Lemaître-Robertson-Walker (FLRW) metric
d s F L R W 2 = d τ 2 + a 2 ( τ ) [ d r 2 1 k r 2 + r 2 ( d θ 2 + sin 2 θ d ϕ 2 ) ] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazcqWGZbWCdaqhaaWcbaGaemOrayKaemitaWKaemOuaiLaem4vaCfabaGaeGOmaidaaOGaeyypa0JaeyOeI0IaemizaqgcciGae8hXdq3aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkcqWGHbqydaahaaWcbeqaaiabikdaYaaakiabcIcaOiab=r8a0jabcMcaPmaadmaabaWaaSaaaeaacqWGKbazcqWGYbGCdaahaaWcbeqaaiabikdaYaaaaOqaaiabigdaXiabgkHiTiabdUgaRjabdkhaYnaaCaaaleqabaGaeGOmaidaaaaakiabgUcaRiabdkhaYnaaCaaaleqabaGaeGOmaidaaOGaeiikaGIaemizaqMae8hUde3aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkcyGGZbWCcqGGPbqAcqGGUbGBdaahaaWcbeqaaiabikdaYaaakiab=H7aXjabdsgaKjab=v9aQnaaCaaaleqabaGaeGOmaidaaOGaeiykaKcacaGLBbGaayzxaaaaaa@636B@
(6)
describes a homogeneous and isotropic universe. Here τ is cosmological time, (r, θ, ϕ) are comoving coordinates, a is the scale factor and k = 0, ±1 the curvature index. The proper radial distance is defined as ar. A useful alternative form of the FLRW metric is
d s F L R W 2 = d τ 2 + a 2 ( τ ) [ d χ 2 + r 2 ( χ ; k ) ( d θ 2 + sin 2 θ d ϕ 2 ) ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazcqWGZbWCdaqhaaWcbaGaemOrayKaemitaWKaemOuaiLaem4vaCfabaGaeGOmaidaaOGaeyypa0JaeyOeI0IaemizaqgcciGae8hXdq3aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkcqWGHbqydaahaaWcbeqaaiabikdaYaaakiabcIcaOiab=r8a0jabcMcaPiabcUfaBjabdsgaKjab=D8aJnaaCaaaleqabaGaeGOmaidaaOGaey4kaSIaemOCai3aaWbaaSqabeaacqaIYaGmaaGccqGGOaakcqWFhpWycqGG7aWocqWGRbWAcqGGPaqkcqGGOaakcqWGKbazcqWF4oqCdaahaaWcbeqaaiabikdaYaaakiabgUcaRiGbcohaZjabcMgaPjabc6gaUnaaCaaaleqabaGaeGOmaidaaOGae8hUdeNaemizaqMae8x1dO2aaWbaaSqabeaacqaIYaGmaaGccqGGPaqkcqGGDbqxcqGGSaalaaa@64FD@
(7)
with
r = r ( χ ; k ) = { sin χ , k = 1 , χ , k = 0 , sinh χ , k = 1. , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGYbGCcqGH9aqpcqWGYbGCcqGGOaakiiGacqWFhpWycqGG7aWocqWGRbWAcqGGPaqkcqGH9aqpdaGabeqaauaabeqadmaaaeaacyGGZbWCcqGGPbqAcqGGUbGBcqWFhpWyaeaacqGGSaalaeaacqWGRbWAcqGH9aqpcqaIXaqmcqGGSaalaeaacqWFhpWyaeaacqGGSaalaeaacqWGRbWAcqGH9aqpcqaIWaamcqGGSaalaeaacyGGZbWCcqGGPbqAcqGGUbGBcqGGObaAcqWFhpWyaeaacqGGSaalaeaacqWGRbWAcqGH9aqpcqGHsislcqaIXaqmcqGGUaGlaaaacaGL7baacqGGSaalaaa@5869@
(8)

χ being an other comoving radial coordinate.

If a photon stream emitted by an astrophysical light source travel without collisions, the number of photons dN γ from a comoving elementary volume of the 6-dimensional phase space ( x MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWG4baEgaWcaaaa@2E37@ , p MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGWbaCgaWcaaaa@2E27@ ) is conserved [47]. Thus the phase space density
f ( t , x , p ) = d N d 3 x d 3 p = d N ω 2 d τ d A d ω d Ω MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGMbGzcqGGOaakcqWG0baDcqGGSaalcuWG4baEgaWcaiabcYcaSiqbdchaWzaalaGaeiykaKIaeyypa0ZaaSaaaeaacqWGKbazcqWGobGtaeaacqWGKbazdaahaaWcbeqaaiabiodaZaaakiqbdIha4zaalaGaemizaq2aaWbaaSqabeaacqaIZaWmaaGccuWGWbaCgaWcaaaacqGH9aqpdaWcaaqaaiabdsgaKjabd6eaobqaaGGaciab=L8a3naaCaaaleqabaGaeGOmaidaaOGaemizaqMae8hXdqNaeeiiaaIaemizaqMaemyqaeKaeeiiaaIaemizaqMae8xYdCNaeeiiaaIaemizaqMaeuyQdCfaaaaa@55B7@
(9)
of a photon stream is constant in time. Here ω denotes the frequency of the photons, dA and dΩ stand for the elementary area normal to the direction of propagation and for the elementary solid angle around the direction of propagation, respectively (see Fig 2). Eq. (9) holds true for any kind of cosmological evolution, provided d3 x MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWG4baEgaWcaaaa@2E37@ dτdA and d3 p MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGWbaCgaWcaaaa@2E27@ ω2dωdΩ are valid for the photons [47]. The luminosity of the source is https://static-content.springer.com/image/art%3A10.1186%2F1754-0410-1-4/MediaObjects/13066_2007_Article_4_IEq11_HTML.gif = dE em /dt em (total energy produced in unit time; the suffix em refers to emission).
https://static-content.springer.com/image/art%3A10.1186%2F1754-0410-1-4/MediaObjects/13066_2007_Article_4_Fig2_HTML.jpg
Figure 2

A schematic representation of the propagation in the curved space-time of the light emitted by a supernova explosion in a distant galaxy and collected on the telescope mirror. A dimensional magnification (à la Wheeler) shows the elementary area dA normal to the direction of propagation n and the elementary solid angle dΩ around n.

A telescope detects the photon flux MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFXeIraaa@3786@ = dE rec / rec /A M (the suffix rec refers to reception). This is the energy detected during unit time on the telescope mirror surface A M . (The surface A M is understood to be perpendicular to the incident light stream.)

From their definition, one can easily find a relation between MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFXeIraaa@3786@ and https://static-content.springer.com/image/art%3A10.1186%2F1754-0410-1-4/MediaObjects/13066_2007_Article_4_IEq11_HTML.gif :
A M = d E r e c / d τ r e c d E e m / d τ e m . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaamrtHrhAL1wy0L2yHvtyaeHbnfgDOvwBHrxAJfwnaGabaiab=ftigjabdgeabnaaBaaaleaacqWGnbqtaeqaaaGcbaGae8NeHWeaaiabg2da9maalaaabaGaemizaqMaemyrau0aaSbaaSqaaiabdkhaYjabdwgaLjabdogaJbqabaGccqGGVaWlcqWGKbaziiGacqGFepaDdaWgaaWcbaGaemOCaiNaemyzauMaem4yamgabeaaaOqaaiabdsgaKjabdweafnaaBaaaleaacqWGLbqzcqWGTbqBaeqaaOGaei4la8IaemizaqMae4hXdq3aaSbaaSqaaiabdwgaLjabd2gaTbqabaaaaOGaeiOla4caaa@581A@
(10)
As the energy of the photon stream in the comoving elementary phase space volume is dE = ħω dN, from Eq. (9) we find
d E r e c d E e m = A M A t o t ( ω r e c ω e m ) 3 d ω r e c d ω e m d A r e c d A e m d τ r e c d τ e m . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaiabdsgaKjabdweafnaaBaaaleaacqWGYbGCcqWGLbqzcqWGJbWyaeqaaaGcbaGaemizaqMaemyrau0aaSbaaSqaaiabdwgaLjabd2gaTbqabaaaaOGaeyypa0ZaaSaaaeaacqWGbbqqdaWgaaWcbaGaemyta0eabeaaaOqaaiabdgeabnaaBaaaleaacqWG0baDcqWGVbWBcqWG0baDaeqaaaaakmaabmaabaWaaSaaaeaaiiGacqWFjpWDdaWgaaWcbaGaemOCaiNaemyzauMaem4yamgabeaaaOqaaiab=L8a3naaBaaaleaacqWGLbqzcqWGTbqBaeqaaaaaaOGaayjkaiaawMcaamaaCaaaleqabaGaeG4mamdaaOWaaSaaaeaacqWGKbazcqWFjpWDdaWgaaWcbaGaemOCaiNaemyzauMaem4yamgabeaaaOqaaiabdsgaKjab=L8a3naaBaaaleaacqWGLbqzcqWGTbqBaeqaaaaakmaalaaabaGaemizaqMaemyqae0aaSbaaSqaaiabdkhaYjabdwgaLjabdogaJbqabaaakeaacqWGKbazcqWGbbqqdaWgaaWcbaGaemyzauMaemyBa0gabeaaaaGcdaWcaaqaaiabdsgaKjab=r8a0naaBaaaleaacqWGYbGCcqWGLbqzcqWGJbWyaeqaaaGcbaGaemizaqMae8hXdq3aaSbaaSqaaiabdwgaLjabd2gaTbqabaaaaOGaeiOla4caaa@76FD@
(11)

Here we have used that from the isotropy of the FLRW universe dΩ rec = dΩ em and we integrate the first to the solid angle encompassing the mirror surface, the second to the whole solid angle (cf. the definitions of E rec , E em ). In Eq. (11) A tot represents the proper area of a sphere centered in the light source and containing the reception point on its surface, at the time of reception.

Due to cosmological evolution the elementary area dA changes as a2 and the frequency of the light is redshifted during cosmic expansion, ω 1/a [47]. In the cosmological evolution of the comoving elementary phase space volume element changes accordingly: 1/a. Therefore
= 1 A t o t ( a a 0 ) 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaamrtHrhAL1wy0L2yHvtyaeHbnfgDOvwBHrxAJfwnaGabaiab=ftigbqaaiab=jrimbaacqGH9aqpdaWcaaqaaiabigdaXaqaaiabdgeabnaaBaaaleaacqWG0baDcqWGVbWBcqWG0baDaeqaaaaakmaabmaabaWaaSaaaeaacqWGHbqyaeaacqWGHbqydaWgaaWcbaGaeGimaadabeaaaaaakiaawIcacaGLPaaadaahaaWcbeqaaiabikdaYaaakiabcYcaSaaa@4786@
(12)
where a0 is the present value of the scale factor, and a is understood to be the scale factor at emission time. In the FLRW universe the proper area of a sphere with comoving radius r em is A t o t = 4 π a 0 2 r e m 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGbbqqdaWgaaWcbaGaemiDaqNaem4Ba8MaemiDaqhabeaakiabg2da9iabisda0GGaciab=b8aWjabdggaHnaaDaaaleaacqaIWaamaeaacqaIYaGmaaGccqWGYbGCdaqhaaWcbaGaemyzauMaemyBa0gabaGaeGOmaidaaOGaeiilaWcaaa@3F84@ and the redshift z is defined as
1 + z = a 0 a e m . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqaIXaqmcqGHRaWkcqWG6bGEcqGH9aqpdaWcaaqaaiabdggaHnaaBaaaleaacqaIWaamaeqaaaGcbaGaemyyae2aaSbaaSqaaiabdwgaLjabd2gaTbqabaaaaOGaeiOla4caaa@389B@
(13)
The luminosity distance d L is defined as in Euclidean geometry:
d L ( z ) : = ( 4 π ) 1 / 2 = a 0 r e m ( 1 + z ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaWgaaWcbaGaemitaWeabeaakiabcIcaOiabdQha6jabcMcaPiabcQda6iabg2da9maabmaabaWaaSaaaeaat0uy0HwzTfgDPnwy1egaryqtHrhAL1wy0L2yHvdaiqaacqWFsectaeaacqaI0aaniiGacqGFapaCcqWFXeIraaaacaGLOaGaayzkaaWaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaIYaGmaaGccqGH9aqpcqWGHbqydaWgaaWcbaGaeGimaadabeaakiabdkhaYnaaBaaaleaacqWGLbqzcqWGTbqBaeqaaOGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeiOla4caaa@5567@
(14)

This definition is rigorous as long as we are dealing with the (homogeneous and isotropic) FLRW universe (irrespective of the value of the curvature index k) and the radius of a sphere is measured in the proper distance ra (the FLRW metric (6) guarantees that the surface of a sphere with radius ra is 4πa2r2).

According to Eq. (7) the comoving coordinate r em can be written in terms of an other radial comoving coordinate χ em (representing the location of the source):
https://static-content.springer.com/image/art%3A10.1186%2F1754-0410-1-4/MediaObjects/13066_2007_Article_4_Equ15_HTML.gif
(15)
d L (z) ≡ a0 (1 + z) (χ em ; k).
Disregarding possible deflections by perturbations of the FLRW universe, a light ray follows radial null geodesics of the FLRW metric, characterized by = /a(τ) = da/a2H. Here H = a ˙ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuWGHbqygaGaaaaa@2E00@ /a is the Hubble parameter. Then
χ e m = χ ( a e m ) = a e m a 0 d a a a ˙ = a e m a 0 d a a 2 H ( a ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaWgaaWcbaGaemyzauMaemyBa0gabeaakiabg2da9iab=D8aJjabcIcaOiabdggaHnaaBaaaleaacqWGLbqzcqWGTbqBaeqaaOGaeiykaKIaeyypa0Zaa8qmaeaadaWcaaqaaiabdsgaKjabdggaHbqaaiabdggaHjqbdggaHzaacaaaaaWcbaGaemyyae2aaSbaaWqaaiabdwgaLjabd2gaTbqabaaaleaacqWGHbqydaWgaaadbaGaeGimaadabeaaa0Gaey4kIipakiabg2da9maapedabaWaaSaaaeaacqWGKbazcqWGHbqyaeaacqWGHbqydaahaaWcbeqaaiabikdaYaaakiabdIeaijabcIcaOiabdggaHjabcMcaPaaaaSqaaiabdggaHnaaBaaameaacqWGLbqzcqWGTbqBaeqaaaWcbaGaemyyae2aaSbaaWqaaiabicdaWaqabaaaniabgUIiYdGccqGGUaGlaaa@5D0C@
(16)
By employing Eq. (13) the radial variable χ can also be expressed in terms of an integral over the redshift as
χ e m ( z ) = 1 a 0 0 z d z H ( z ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaWgaaWcbaGaemyzauMaemyBa0gabeaakiabcIcaOiabdQha6jabcMcaPiabg2da9maalaaabaGaeGymaedabaGaemyyae2aaSbaaSqaaiabicdaWaqabaaaaOWaa8qmaeaadaWcaaqaaiabdsgaKjqbdQha6zaafaaabaGaemisaGKaeiikaGIafmOEaONbauaacqGGPaqkaaaaleaacqaIWaamaeaacqWG6bGEa0Gaey4kIipakiabc6caUaaa@45BA@
(17)

which completes the definition (15) of the luminosity distance d L in terms of the redshift z.

Differentiating Eq. (15) with χ given by Eq. (17) with respect to z gives
1 H ( z ) = [ 1 k d L 2 ( z ) a 0 2 ( 1 + z ) 2 ] 1 / 2 d d z [ d L ( z ) 1 + z ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaiabigdaXaqaaiabdIeaijabcIcaOiabdQha6jabcMcaPaaacqGH9aqpdaWadaqaaiabigdaXiabgkHiTmaalaaabaGaem4AaSMaemizaq2aa0baaSqaaiabdYeambqaaiabikdaYaaakiabcIcaOiabdQha6jabcMcaPaqaaiabdggaHnaaDaaaleaacqaIWaamaeaacqaIYaGmaaGccqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkdaahaaWcbeqaaiabikdaYaaaaaaakiaawUfacaGLDbaadaahaaWcbeqaaiabgkHiTiabigdaXiabc+caViabikdaYaaakmaalaaabaGaemizaqgabaGaemizaqMaemOEaOhaamaadmaabaWaaSaaaeaacqWGKbazdaWgaaWcbaGaemitaWeabeaakiabcIcaOiabdQha6jabcMcaPaqaaiabigdaXiabgUcaRiabdQha6baaaiaawUfacaGLDbaacqGGSaalaaa@5CB1@
(18)

therefore if independent measurements of d L and z are available for a set of light sources, the Hubble-parameter H(z) and in consequence the cosmological dynamics can be determined.

From the combined measurements of the large-scale structure of the Universe [48, 49] and of the structure of the cosmic microwave background [50] the conclusion was reached that the space geometry has flat spatial sections. Therefore in what follows we consider k = 0. Then the luminosity distance-redshift relation becomes
d L ( z ) = ( 1 + z ) 0 z d z H ( z ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaWgaaWcbaGaemitaWeabeaakiabcIcaOiabdQha6jabcMcaPiabg2da9iabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPmaapedabaWaaSaaaeaacqWGKbazcuWG6bGEgaqbaaqaaiabdIeaijabcIcaOiqbdQha6zaafaGaeiykaKcaaaWcbaGaeGimaadabaGaemOEaOhaniabgUIiYdGccqGGUaGlaaa@454A@
(19)

In practice, the function d L (z) is conveniently measured with distant supernovae of type Ia. The luminosity is evaluated by photometry, while the redshift from spectroscopic analysis of the host galaxy.

Each cosmological model has its own prediction for the shape of the function d L (z) [see Eq. (15) with χ given by Eq. (17) for generic k, or Eq. (19) for k = 0]. This is how the measured d L (z) data turn into a cosmological test.

3 The luminosity-redshift relation in Randall-Sundrum type II brane-worlds

We consider FLRW branes with k = 0 and brane cosmological constant Λ, embedded symmetrically. The bulk is the Vaidya-anti de Sitter space-time with cosmological constant Λ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuqHBoatgaacaaaa@2E30@ , and it contains bulk black holes with masses m on both sides of the brane. The black hole masses can change if the brane radiates into the bulk. An ansatz comparable with structure formation has been advanced in [51] for the Weyl fluid m/a4 for the case when the brane radiates, m = m0a α , where m0 is a constant and α = 2, 3. For α = 0 the Weyl fluid is known as dark radiation and then the bulk space-time becomes Schwarzschild-anti de Sitter. The brane tension and the two cosmological constants are inter-related as
2 Λ = κ 2 λ + κ ˜ 2 Λ ˜ . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqaIYaGmcqqHBoatcqGH9aqpiiGacqWF6oWAdaahaaWcbeqaaiabikdaYaaakiab=T7aSjabgUcaRiqb=P7aRzaaiaWaaWbaaSqabeaacqaIYaGmaaGccuqHBoatgaacaiabc6caUaaa@3AD9@
(20)
The Friedmann equation gives the Hubble parameter to Λ, m, the scale factor a and the matter energy density ρ on the brane:
H 2 = Λ 3 + κ 2 ρ 3 ( 1 + ρ 2 λ ) + 2 m 0 a 4 α . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGibasdaahaaWcbeqaaiabikdaYaaakiabg2da9maalaaabaGaeu4MdWeabaGaeG4mamdaaiabgUcaRmaalaaabaacciGae8NUdS2aaWbaaSqabeaacqaIYaGmaaGccqWFbpGCaeaacqaIZaWmaaWaaeWaaeaacqaIXaqmcqGHRaWkdaWcaaqaaiab=f8aYbqaaiabikdaYiab=T7aSbaaaiaawIcacaGLPaaacqGHRaWkdaWcaaqaaiabikdaYiabd2gaTnaaBaaaleaacqaIWaamaeqaaaGcbaGaemyyae2aaWbaaSqabeaacqaI0aancqGHsislcqWFXoqyaaaaaOGaeiOla4caaa@4B05@
(21)
In the matter dominated era the brane is dominated by dust, obeying the continuity equation
ρ ˙ + 3 H ρ = 0 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacuWFbpGCgaGaaiabgUcaRiabiodaZiabdIeaijab=f8aYjabg2da9iabicdaWiabcYcaSaaa@35FA@
(22)

which gives ρ ~ a-3. We introduce the following dimensionless quantities:

Ω tot = ΩΛ + Ω ρ + Ω λ + Ω d ,

Ω ρ = κ 2 ρ 0 3 H 0 2 , Ω λ = κ 2 ρ 0 2 6 λ H 0 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeqacaaabaGaeuyQdC1aaSbaaSqaaGGaciab=f8aYbqabaGccqGH9aqpdaWcaaqaaiab=P7aRnaaCaaaleqabaGaeGOmaidaaOGae8xWdi3aaSbaaSqaaiabicdaWaqabaaakeaacqaIZaWmcqWGibasdaqhaaWcbaGaeGimaadabaGaeGOmaidaaaaakiabcYcaSaqaaiabfM6axnaaBaaaleaacqWF7oaBaeqaaOGaeyypa0ZaaSaaaeaacqWF6oWAdaahaaWcbeqaaiabikdaYaaakiab=f8aYnaaDaaaleaacqaIWaamaeaacqaIYaGmaaaakeaacqaI2aGncqWF7oaBcqWGibasdaqhaaWcbaGaeGimaadabaGaeGOmaidaaaaakiabcYcaSaaaaaa@4DFD@
(24)
Ω d = 2 m 0 a 0 4 α H 0 2 , Ω Λ = Λ 3 H 0 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeqacaaabaGaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqGH9aqpdaWcaaqaaiabikdaYiabd2gaTnaaBaaaleaacqaIWaamaeqaaaGcbaGaemyyae2aa0baaSqaaiabicdaWaqaaiabisda0iabgkHiTGGaciab=f7aHbaakiabdIeainaaDaaaleaacqaIWaamaeaacqaIYaGmaaaaaOGaeiilaWcabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqGH9aqpdaWcaaqaaiabfU5ambqaaiabiodaZiabdIeainaaDaaaleaacqaIWaamaeaacqaIYaGmaaaaaOGaeiOla4caaaaa@4932@
(25)
The subscript 0 denotes the present value of the respective quantities. In terms of these notations the Friedmann equation becomes
H 2 = H 0 2 [ Ω Λ + Ω ρ a 0 3 a 3 + Ω d a 0 4 α a 4 α + Ω λ a 0 6 a 6 ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGibasdaahaaWcbeqaaiabikdaYaaakiabg2da9iabdIeainaaDaaaleaacqaIWaamaeaacqaIYaGmaaGcdaWadaqaaiabfM6axnaaBaaaleaacqqHBoataeqaaOGaey4kaSIaeuyQdC1aaSbaaSqaaGGaciab=f8aYbqabaGcdaWcaaqaaiabdggaHnaaDaaaleaacqaIWaamaeaacqaIZaWmaaaakeaacqWGHbqydaahaaWcbeqaaiabiodaZaaaaaGccqGHRaWkcqqHPoWvdaWgaaWcbaGaemizaqgabeaakmaalaaabaGaemyyae2aa0baaSqaaiabicdaWaqaaiabisda0iabgkHiTiab=f7aHbaaaOqaaiabdggaHnaaCaaaleqabaGaeGinaqJaeyOeI0Iae8xSdegaaaaakiabgUcaRiabfM6axnaaBaaaleaacqWF7oaBaeqaaOWaaSaaaeaacqWGHbqydaqhaaWcbaGaeGimaadabaGaeGOnaydaaaGcbaGaemyyae2aaWbaaSqabeaacqaI2aGnaaaaaaGccaGLBbGaayzxaaGaeiOla4caaa@5CC8@
(26)
In particular at present time this gives Ω tot = 1. Then the radial coordinate (16) becomes
χ e m = 1 H 0 a e m a 0 a d a [ Ω Λ a 6 + Ω ρ a 0 3 a 3 + Ω d a 0 4 α a α + 2 + Ω λ a 0 6 ] 1 / 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaWgaaWcbaGaemyzauMaemyBa0gabeaakiabg2da9maalaaabaGaeGymaedabaGaemisaG0aaSbaaSqaaiabicdaWaqabaaaaOWaa8qmaeaadaWcaaqaaiabdggaHjabdsgaKjabdggaHbqaamaadmaabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqWGHbqydaahaaWcbeqaaiabiAda2aaakiabgUcaRiabfM6axnaaBaaaleaacqWFbpGCaeqaaOGaemyyae2aa0baaSqaaiabicdaWaqaaiabiodaZaaakiabdggaHnaaCaaaleqabaGaeG4mamdaaOGaey4kaSIaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqWGHbqydaqhaaWcbaGaeGimaadabaGaeGinaqJaeyOeI0Iae8xSdegaaOGaemyyae2aaWbaaSqabeaacqWFXoqycqGHRaWkcqaIYaGmaaGccqGHRaWkcqqHPoWvdaWgaaWcbaGae83UdWgabeaakiabdggaHnaaDaaaleaacqaIWaamaeaacqaI2aGnaaaakiaawUfacaGLDbaadaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaaaaaabaGaemyyae2aaSbaaWqaaiabdwgaLjabd2gaTbqabaaaleaacqWGHbqydaWgaaadbaGaeGimaadabeaaa0Gaey4kIipakiabc6caUaaa@6EA5@
(27)

This is a complicated integral, which cannot be computed analytically in the majority of cases. In what follows we will analyze various specific cases of the above integral, when an analytic solution is possible. The cases α = 2, 3 represent the Weyl fluid compatible with structure formation, while α = 0 represents the dark radiation.

4 Branes with Randall-Sundrum fine-tuning

In the original Randall-Sundrum scenario the bulk cosmological constant Λ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuqHBoatgaacaaaa@2E30@ is fine-tuned with the brane tension λ such that cf. Eq. (20) the brane cosmological constant vanishes. For simplicity we also assume throughout this section α = 0. By imposing a vanishing cosmological constant on the brane, ΩΛ = 0 such that the polynomial of rank 6 in the denominator of the integrand in Eq. (27) shrinks to a polynomial of rank 3. Therefore its roots can be found analytically. Following general procedures, the luminosity distance – redshift relation can be then given analytically in terms of elliptic functions. This is done in the following subsection. In the second and third subsections of this chapter we discuss the limits Ω d → 0 (when the bulk is anti de Sitter) and the late-time universe limit ρ/λ → 0. The general relativistic (Einstein-deSitter) limit is found in the fourth subsection, when further Ω λ → 0 is taken.

4.1 Schwarzschild-AdS bulk

With no brane cosmological constant, Eq. (27) becomes:
χ e m = 1 a 0 3 / 2 H 0 Ω ρ 1 / 2 a e m a 0 a d a [ a 3 + Ω d Ω ρ a 0 a 2 + Ω λ Ω ρ a 0 3 ] 1 / 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaWgaaWcbaGaemyzauMaemyBa0gabeaakiabg2da9maalaaabaGaeGymaedabaGaemyyae2aa0baaSqaaiabicdaWaqaaiabiodaZiabc+caViabikdaYaaakiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabigdaXiabc+caViabikdaYaaaaaGcdaWdXaqaamaalaaabaGaemyyaeMaemizaqMaemyyaegabaWaamWaaeaacqWGHbqydaahaaWcbeqaaiabiodaZaaakiabgUcaRmaalaaabaGaeuyQdC1aaSbaaSqaaiabdsgaKbqabaaakeaacqqHPoWvdaWgaaWcbaGae8xWdihabeaaaaGccqWGHbqydaWgaaWcbaGaeGimaadabeaakiabdggaHnaaCaaaleqabaGaeGOmaidaaOGaey4kaSYaaSaaaeaacqqHPoWvdaWgaaWcbaGae83UdWgabeaaaOqaaiabfM6axnaaBaaaleaacqWFbpGCaeqaaaaakiabdggaHnaaDaaaleaacqaIWaamaeaacqaIZaWmaaaakiaawUfacaGLDbaadaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaaaaaabaGaemyyae2aaSbaaWqaaiabdwgaLjabd2gaTbqabaaaleaacqWGHbqydaWgaaadbaGaeGimaadabeaaa0Gaey4kIipakiabc6caUaaa@6DCD@
(28)
Following the method given in [52] we find the following roots of the denominator:
α = 2 Ω d a 0 3 Ω ρ ( 1 + 2 cosh Ψ 3 ) , β = Ω d a 0 3 Ω ρ ( 1 + cosh Ψ 3 + i 3 sinh Ψ 3 ) ,
(29)
and its complex conjugate β*. The auxiliary quantity Ψ is defined as
cosh Ψ = 1 + 27 Ω λ Ω ρ 2 2 Ω d 3 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacyGGJbWycqGGVbWBcqGGZbWCcqGGObaAcqqHOoqwcqGH9aqpcqaIXaqmcqGHRaWkdaWcaaqaaiabikdaYiabiEda3iabfM6axnaaBaaaleaaiiGacqWF7oaBaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabikdaYaaaaOqaaiabikdaYiabfM6axnaaDaaaleaacqWGKbazaeaacqaIZaWmaaaaaOGaeiOla4caaa@465E@
(30)
We introduce the following real combinations of the complex roots
b 1 = β + β 2 = Ω d a 0 3 Ω ρ ( 1 + cosh Ψ 3 ) , a 1 = β β 2 i = Ω d a 0 3 Ω ρ sinh Ψ 3 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaGaemOyai2aaSbaaSqaaiabigdaXaqabaGccqGH9aqpdaWcaaqaaGGaciab=j7aIjabgUcaRiab=j7aInaaCaaaleqabaGaey4fIOcaaaGcbaGaeGOmaidaaiabg2da9maalaaabaGaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqWGHbqydaWgaaWcbaGaeGimaadabeaaaOqaaiabiodaZiabfM6axnaaBaaaleaacqWFbpGCaeqaaaaakmaabmaabaGaeyOeI0IaeGymaeJaey4kaSIagi4yamMaei4Ba8Maei4CamNaeiiAaG2aaSaaaeaacqqHOoqwaeaacqaIZaWmaaaacaGLOaGaayzkaaGaeiilaWcabaGaemyyae2aaSbaaSqaaiabigdaXaqabaGccqGH9aqpdaWcaaqaaiab=j7aIjabgkHiTiab=j7aInaaCaaaleqabaGaey4fIOcaaaGcbaGaeGOmaiJaemyAaKgaaiabg2da9maalaaabaGaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqWGHbqydaWgaaWcbaGaeGimaadabeaaaOqaamaakaaabaGaeG4mamdaleqaaOGaeuyQdC1aaSbaaSqaaiab=f8aYbqabaaaaOGagi4CamNaeiyAaKMaeiOBa4MaeiiAaG2aaSaaaeaacqqHOoqwaeaacqaIZaWmaaGaeiOla4caaaaa@6E11@
(31)
Then Eq. (28) is written conveniently as
χ e m = 1 a 0 3 / 2 H 0 Ω ρ 1 / 2 a e m a 0 a d a ( a α ) [ ( a b 1 ) 2 + a 1 2 ] 1 / 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWFhpWydaWgaaWcbaGaemyzauMaemyBa0gabeaakiabg2da9maalaaabaGaeGymaedabaGaemyyae2aa0baaSqaaiabicdaWaqaaiabiodaZiabc+caViabikdaYaaakiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiabeg8aYbqaaiabigdaXiabc+caViabikdaYaaaaaGcdaWdXaqaamaalaaabaGaemyyaeMaemizaqMaemyyaegabaGaeiikaGIaemyyaeMaeyOeI0Iae8xSdeMaeiykaKYaamWaaeaacqGGOaakcqWGHbqycqGHsislcqWGIbGydaWgaaWcbaGaeGymaedabeaakiabcMcaPmaaCaaaleqabaGaeGOmaidaaOGaey4kaSIaemyyae2aa0baaSqaaiabigdaXaqaaiabikdaYaaaaOGaay5waiaaw2faamaaCaaaleqabaGaeGymaeJaei4la8IaeGOmaidaaaaaaeaacqWGHbqydaWgaaadbaGaemyzauMaemyBa0gabeaaaSqaaiabdggaHnaaBaaameaacqaIWaamaeqaaaqdcqGHRiI8aOGaeiOla4caaa@64FB@
(32)
The integration can be carried out by employing the formulae (239.07) and (341.53) of Ref. [53]. We obtain
χ e m = 1 a 0 3 / 2 H 0 Ω ρ 1 / 2 { g ( α + A ) [ F ( ϕ 0 , ε ) F ( ϕ e m , ε ) ] 2 g A [ E ( ϕ 0 , ε ) E ( ϕ e m , ε ) + 2 g A [ sin ϕ 0 1 ε 2 sin 2 ϕ 0 1 + cos ϕ 0 sin ϕ e m 1 ε 2 sin 2 ϕ e m 1 + cos ϕ e m ] } ,
(33)
where F (ϕ, ε) is the elliptic integral of the first kind; E (ϕ, ε) is the elliptic integral of the second kind (with variable ϕand argument ε); and we have introduced the following standard notations, cf. Ref. [53]:
A 2 = ( b 1 α ) 2 + a 1 2 , ε 2 = A + b 1 α 2 A , g = 1 A , ϕ = arccos ( A + α a A α + a ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeabbaaaaeaacqWGbbqqdaahaaWcbeqaaiabikdaYaaakiabg2da9iabcIcaOiabdkgaInaaBaaaleaacqaIXaqmaeqaaOGaeyOeI0ccciGae8xSdeMaeiykaKYaaWbaaSqabeaacqaIYaGmaaGccqGHRaWkcqWGHbqydaqhaaWcbaGaeGymaedabaGaeGOmaidaaOGaeiilaWcabaGae8xTdu2aaWbaaSqabeaacqaIYaGmaaGccqGH9aqpdaWcaaqaaiabdgeabjabgUcaRiabdkgaInaaBaaaleaacqaIXaqmaeqaaOGaeyOeI0Iae8xSdegabaGaeGOmaiJaemyqaeeaaiabcYcaSaqaaiabdEgaNjabg2da9maalaaabaGaeGymaedabaWaaOaaaeaacqWGbbqqaSqabaaaaOGaeiilaWcabaGae8x1dOMaeyypa0JagiyyaeMaeiOCaiNaei4yamMaei4yamMaei4Ba8Maei4Cam3aaeWaaeaadaWcaaqaaiabdgeabjabgUcaRiab=f7aHjabgkHiTiabdggaHbqaaiabdgeabjabgkHiTiab=f7aHjabgUcaRiabdggaHbaaaiaawIcacaGLPaaacqGGUaGlaaaaaa@68FF@
(34)
By employing Eqs. (13), (17) and (19), after a lengthy, but straightforward calculation, the luminosity distance-redshift relation emerges:
d L λ d = ( 1 + z ) Ω d 1 / 2 H 0 Ω ρ { ( B 1 + B 2 ) [ F ( ϕ 0 , ε ) F ( ϕ e m , ε ) ] + 2 B 2 [ E ( ϕ e m , ε ) E ( ϕ 0 , ε ) + sin ϕ 0 1 ε 2 sin 2 ϕ 0 1 + cos ϕ 0 sin ϕ e m 1 ε 2 sin 2 ϕ e m 1 + cos ϕ e m ] } ,
(35)
with
B 1 = 3 1 / 4 ( 1 + 2 cosh Ψ 3 ) 3 ( 4 cosh 2 Ψ 3 1 ) 1 / 4 , B 2 = ( 4 cosh 2 Ψ 3 1 3 ) 1 / 4 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaGaemOqai0aaSbaaSqaaiabigdaXaqabaGccqGH9aqpdaWcaaqaaiabgkHiTiabiodaZmaaCaaaleqabaGaeGymaeJaei4la8IaeGinaqdaaOWaaeWaaeaacqaIXaqmcqGHRaWkcqaIYaGmcyGGJbWycqGGVbWBcqGGZbWCcqGGObaAdaWcaaqaaiabfI6azbqaaiabiodaZaaaaiaawIcacaGLPaaaaeaacqaIZaWmdaqadaqaaiabisda0iGbcogaJjabc+gaVjabcohaZjabcIgaOnaaCaaaleqabaGaeGOmaidaaOWaaSaaaeaacqqHOoqwaeaacqaIZaWmaaGaeyOeI0IaeGymaedacaGLOaGaayzkaaWaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaI0aanaaaaaOGaeiilaWcabaGaemOqai0aaSbaaSqaaiabikdaYaqabaGccqGH9aqpdaqadaqaamaalaaabaGaeGinaqJagi4yamMaei4Ba8Maei4CamNaeiiAaG2aaWbaaSqabeaacqaIYaGmaaGcdaWcaaqaaiabfI6azbqaaiabiodaZaaacqGHsislcqaIXaqmaeaacqaIZaWmaaaacaGLOaGaayzkaaWaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaI0aanaaGccqGGUaGlaaaaaa@6966@
(36)
and
ε 2 = 1 2 + 3 cosh Ψ 3 2 4 cosh 2 Ψ 3 1 , ϕ = arccos { ( 1 + z ) Ω d [ 12 cosh 2 Ψ 3 3 1 2 cosh Ψ 3 ] 3 Ω ρ ( 1 + z ) Ω d [ 12 cosh 2 Ψ 3 3 + 1 + 2 cosh Ψ 3 ] + 3 Ω ρ } . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaacciGae8xTdu2aaWbaaSqabeaacqaIYaGmaaGccqGH9aqpdaWcaaqaaiabigdaXaqaaiabikdaYaaacqGHRaWkdaWcaaqaamaakaaabaGaeG4mamdaleqaaOGagi4yamMaei4Ba8Maei4CamNaeiiAaG2aaSaaaeaacqqHOoqwaeaacqaIZaWmaaaabaGaeGOmaiZaaOaaaeaacqaI0aancyGGJbWycqGGVbWBcqGGZbWCcqGGObaAdaahaaWcbeqaaiabikdaYaaakmaalaaabaGaeuiQdKfabaGaeG4mamdaaiabgkHiTiabigdaXaWcbeaaaaGccqGGSaalaeaacqWFvpGAcqGH9aqpcyGGHbqycqGGYbGCcqGGJbWycqGGJbWycqGGVbWBcqGGZbWCdaGadeqaamaalaaabaGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGcdaWadaqaamaakaaabaGaeGymaeJaeGOmaiJagi4yamMaei4Ba8Maei4CamNaeiiAaG2aaWbaaSqabeaacqaIYaGmaaGcdaWcaaqaaiabfI6azbqaaiabiodaZaaacqGHsislcqaIZaWmaSqabaGccqGHsislcqaIXaqmcqGHsislcqaIYaGmcyGGJbWycqGGVbWBcqGGZbWCcqGGObaAdaWcaaqaaiabfI6azbqaaiabiodaZaaaaiaawUfacaGLDbaacqGHsislcqaIZaWmcqqHPoWvdaWgaaWcbaGae8xWdihabeaaaOqaaiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPiabfM6axnaaBaaaleaacqWGKbazaeqaaOWaamWaaeaadaGcaaqaaiabigdaXiabikdaYiGbcogaJjabc+gaVjabcohaZjabcIgaOnaaCaaaleqabaGaeGOmaidaaOWaaSaaaeaacqqHOoqwaeaacqaIZaWmaaGaeyOeI0IaeG4mamdaleqaaOGaey4kaSIaeGymaeJaey4kaSIaeGOmaiJagi4yamMaei4Ba8Maei4CamNaeiiAaG2aaSaaaeaacqqHOoqwaeaacqaIZaWmaaaacaGLBbGaayzxaaGaey4kaSIaeG4mamJaeuyQdC1aaSbaaSqaaiab=f8aYbqabaaaaaGccaGL7bGaayzFaaGaeiOla4caaaaa@A943@
(37)
Here ϕ runs in the range 0..π/2. In computing for other values of ϕ, we can use the following addition rules for the elliptic integrals:
F ( m π ± ϕ , ε ) = 2 m K ± F ( ϕ , ε ) , E ( m π ± ϕ , ε ) = 2 m E ± E ( ϕ , ε ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeGabaaabaGaemOrayKaeiikaGIaemyBa0gcciGae8hWdaNaeyySaeRae8x1dOMaeiilaWIae8xTduMaeiykaKIaeyypa0JaeGOmaiJaemyBa0Maem4saSKaeyySaeRaemOrayKaeiikaGIae8x1dOMaeiilaWIae8xTduMaeiykaKIaeiilaWcabaGaemyrauKaeiikaGIaemyBa0Mae8hWdaNaeyySaeRae8x1dOMaeiilaWIae8xTduMaeiykaKIaeyypa0JaeGOmaiJaemyBa0MaemyrauKaeyySaeRaemyrauKaeiikaGIae8x1dOMaeiilaWIae8xTduMaeiykaKIaeiilaWcaaaaa@6197@
(38)

where K and E are the complete elliptic integrals of the first and second kind. Eqs. (30) and (35)–(37) represent the analytical expression of the luminosity distance-redshift relation for FLRW branes with Randall-Sundrum fine-tuning. They are given in terms of the well-known elliptic integrals of first and second kind, and the cosmological parameters Ω ρ , Ω λ and Ω d .

4.2 Limit of no black hole in the bulk

In this subsection we consider the case Ω d = 0. The derivation follows closely the steps of the previous subsection, however the formulae are simpler. The auxiliary expression (30) for Ψ is well defined only for Ω d ≠ 0 and we have to address the question how to obtain suitable limits of the results derived for Ω d ≠ 0. For any Ω d 1
cosh Ψ 27 Ω λ Ω ρ 2 2 Ω d 3 1. MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacyGGJbWycqGGVbWBcqGGZbWCcqGGObaAcqqHOoqwcqGHijYUdaWcaaqaaiabikdaYiabiEda3iabfM6axnaaBaaaleaaiiGacqWF7oaBaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabikdaYaaaaOqaaiabikdaYiabfM6axnaaDaaaleaacqWGKbazaeaacqaIZaWmaaaaaOGaeS4AI8JaeGymaeJaeiOla4caaa@4784@
(39)
But
cosh Ψ = cosh Ψ 3 ( 4 cosh 2 Ψ 3 3 ) 4 cosh 3 Ψ 3 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacyGGJbWycqGGVbWBcqGGZbWCcqGGObaAcqqHOoqwcqGH9aqpcyGGJbWycqGGVbWBcqGGZbWCcqGGObaAdaWcaaqaaiabfI6azbqaaiabiodaZaaadaqadaqaaiabisda0iGbcogaJjabc+gaVjabcohaZjabcIgaOnaaCaaaleqabaGaeGOmaidaaOWaaSaaaeaacqqHOoqwaeaacqaIZaWmaaGaeyOeI0IaeG4mamdacaGLOaGaayzkaaGaeyisISRaeGinaqJagi4yamMaei4Ba8Maei4CamNaeiiAaG2aaWbaaSqabeaacqaIZaWmaaGcdaWcaaqaaiabfI6azbqaaiabiodaZaaacqGGSaalaaa@5725@
(40)
as cosh (Ψ/3) 1 also holds. Thus
cosh Ψ 3 ± sinh Ψ 3 3 Ω λ 1 / 3 Ω ρ 2 / 3 2 Ω d . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacyGGJbWycqGGVbWBcqGGZbWCcqGGObaAdaWcaaqaaiabfI6azbqaaiabiodaZaaacqGHijYUcqGHXcqScyGGZbWCcqGGPbqAcqGGUbGBcqGGObaAdaWcaaqaaiabfI6azbqaaiabiodaZaaacqGHijYUdaWcaaqaaiabiodaZiabfM6axnaaDaaaleaaiiGacqWF7oaBaeaacqaIXaqmcqGGVaWlcqaIZaWmaaGccqqHPoWvdaqhaaWcbaGae8xWdihabaGaeGOmaiJaei4la8IaeG4mamdaaaGcbaGaeGOmaiJaeuyQdC1aaSbaaSqaaiabdsgaKbqabaaaaOGaeiOla4caaa@54A9@
(41)
By employing Eq. (41) in the generic expressions derived in the preceding subsection, we obtain the luminosity distance-redshift relation in a very similar form to Eq. (35), but with different coefficients:
d L λ = 2 3 4 ( 1 + z ) Ω λ 1 / 6 H 0 Ω ρ 2 / 3 { 1 2 ( 1 3 3 ) [ F ( ϕ 0 , ε ) F ( ϕ e m , ε ) ] + E ( ϕ e m , ε ) E ( ϕ 0 , ε ) + sin ϕ 0 1 ε 2 sin 2 ϕ 0 1 + cos ϕ 0 sin ϕ e m 1 ε 2 sin 2 ϕ e m 1 + cos ϕ e m } , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqadeWabaaabaGaemizaq2aa0baaSqaaiabdYeambqaaGGaciab=T7aSbaakiabg2da9maalaaabaGaeGOmaiZaaOqaaeaacqaIZaWmaSqaaiabisda0aaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPiabfM6axnaaDaaaleaacqWF7oaBaeaacqaIXaqmcqGGVaWlcqaI2aGnaaaakeaacqWGibasdaWgaaWcbaGaeGimaadabeaakiabfM6axnaaDaaaleaacqWFbpGCaeaacqaIYaGmcqGGVaWlcqaIZaWmaaaaaOWaaiqabeaadaWcaaqaaiabigdaXaqaaiabikdaYaaadaqadaqaaiabigdaXiabgkHiTmaalaaabaWaaOaaaeaacqaIZaWmaSqabaaakeaacqaIZaWmaaaacaGLOaGaayzkaaGaei4waSLaemOrayKaeiikaGIae8x1dO2aaSbaaSqaaiabicdaWaqabaGccqGGSaalcqWF1oqzcqGGPaqkcqGHsislcqWGgbGrcqGGOaakcqWFvpGAdaWgaaWcbaGaemyzauMaemyBa0gabeaakiabcYcaSiab=v7aLjabcMcaPiabc2faDbGaay5EaaaabaGaey4kaSIaemyrauKaeiikaGIae8x1dO2aaSbaaSqaaiabdwgaLjabd2gaTbqabaGccqGGSaalcqWF1oqzcqGGPaqkcqGHsislcqWGfbqrcqGGOaakcqWFvpGAdaWgaaWcbaGaeGimaadabeaakiabcYcaSiab=v7aLjabcMcaPaqaamaaciqabaGaey4kaSYaaSaaaeaacyGGZbWCcqGGPbqAcqGGUbGBcqWFvpGAdaWgaaWcbaGaeGimaadabeaakmaakaaabaGaeGymaeJaeyOeI0Iae8xTdu2aaWbaaSqabeaacqaIYaGmaaGccyGGZbWCcqGGPbqAcqGGUbGBdaahaaWcbeqaaiabikdaYaaakiab=v9aQnaaBaaaleaacqaIWaamaeqaaaqabaaakeaacqaIXaqmcqGHRaWkcyGGJbWycqGGVbWBcqGGZbWCcqWFvpGAdaWgaaWcbaGaeGimaadabeaaaaGccqGHsisldaWcaaqaaiGbcohaZjabcMgaPjabc6gaUjab=v9aQnaaBaaaleaacqWGLbqzcqWGTbqBaeqaaOWaaOaaaeaacqaIXaqmcqGHsislcqWF1oqzdaahaaWcbeqaaiabikdaYaaakiGbcohaZjabcMgaPjabc6gaUnaaCaaaleqabaGaeGOmaidaaOGae8x1dO2aaSbaaSqaaiabdwgaLjabd2gaTbqabaaabeaaaOqaaiabigdaXiabgUcaRiGbcogaJjabc+gaVjabcohaZjab=v9aQnaaBaaaleaacqWGLbqzcqWGTbqBaeqaaaaaaOGaayzFaaGaeiilaWcaaaaa@BE5B@
(42)
where
ε 2 + 1 2 + 3 4 , ϕ = arccos ( 3 1 ) Ω λ 1 / 3 ( 1 + z ) Ω ρ 1 / 3 ( 3 + 1 ) Ω λ 1 / 3 ( 1 + z ) + Ω ρ 1 / 3 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaacciGae8xTdu2aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkdaWcaaqaaiabigdaXaqaaiabikdaYaaacqGHRaWkdaWcaaqaamaakaaabaGaeG4mamdaleqaaaGcbaGaeGinaqdaaiabcYcaSaqaaiab=v9aQjabg2da9iGbcggaHjabckhaYjabcogaJjabcogaJjabc+gaVjabcohaZnaalaaabaGaeiikaGYaaOaaaeaacqaIZaWmaSqabaGccqGHsislcqaIXaqmcqGGPaqkcqqHPoWvdaqhaaWcbaGae83UdWgabaGaeGymaeJaei4la8IaeG4mamdaaOGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeyOeI0IaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabigdaXiabc+caViabiodaZaaaaOqaaiabcIcaOmaakaaabaGaeG4mamdaleqaaOGaey4kaSIaeGymaeJaeiykaKIaeuyQdC1aa0baaSqaaiab=T7aSbqaaiabigdaXiabc+caViabiodaZaaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPiabgUcaRiabfM6axnaaDaaaleaacqWFbpGCaeaacqaIXaqmcqGGVaWlcqaIZaWmaaaaaOGaeiOla4caaaaa@7053@
(43)

Again, ϕfor this case emerges in the limit Ω d → 0 from the generic expression Eq. (37), by employing Eq. (41) as in the limiting process expressions of the type ∞ × 0 appear.

4.3 Late-time universe limit

In the late-time universe ρ λ and in consequence Ω λ = 0 can be safely assumed. We keep however the dark radiation in the model. Eq. (27) simplifies considerably, and a straightforward integration gives the luminosity distance – redshift relation
d L d ( z ) = 2 1 + z H 0 Ω ρ ( ( Ω ρ + Ω d ) ( 1 + z ) Ω ρ + Ω d ( 1 + z ) ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaemizaqgaaOGaeiikaGIaemOEaONaeiykaKIaeyypa0ZaaSaaaeaacqaIYaGmdaGcaaqaaiabigdaXiabgUcaRiabdQha6bWcbeaaaOqaaiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aaSbaaSqaaGGaciab=f8aYbqabaaaaOWaaeWaaeaadaGcaaqaaiabcIcaOiabfM6axnaaBaaaleaacqWFbpGCaeqaaOGaey4kaSIaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqGGPaqkcqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkaSqabaGccqGHsisldaGcaaqaaiabfM6axnaaBaaaleaacqWFbpGCaeqaaOGaey4kaSIaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkaSqabaaakiaawIcacaGLPaaacqGGUaGlaaa@5D5D@
(44)
We can also prove that this result emerges as the Ω λ → 0 limit from the generic results, Eqs. (30) and (35)–(37). When Ω λ → 0 Eq. (30) gives Ψ → 0. Then Eqs. (37) and (36) give
ε 2 = 1 ϕ = arccos { Ω ρ 2 ( 1 + z ) Ω d + Ω ρ } , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaacciGae8xTdu2aaWbaaSqabeaacqaIYaGmaaGccqGH9aqpcqaIXaqmaeaacqWFvpGAcqGH9aqpcyGGHbqycqGGYbGCcqGGJbWycqGGJbWycqGGVbWBcqGGZbWCdaGadeqaamaalaaabaGaeyOeI0IaeuyQdC1aaSbaaSqaaiab=f8aYbqabaaakeaacqaIYaGmcqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkcqqHPoWvdaWgaaWcbaGaemizaqgabeaakiabgUcaRiabfM6axnaaBaaaleaacqWFbpGCaeqaaaaaaOGaay5Eaiaaw2haaiabcYcaSaaaaaa@5171@
(45)

B2 = 1 = -B1.

By noting that E (ϕ, 1) = sin ϕ, we obtain from Eq. (35):
d L d = 2 ( 1 + z ) Ω d 1 / 2 H 0 Ω ρ [ sin ϕ 0 1 + cos ϕ 0 sin ϕ e m 1 + cos ϕ e m ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaemizaqgaaOGaeyypa0ZaaSaaaeaacqaIYaGmcqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkcqqHPoWvdaqhaaWcbaGaemizaqgabaGaeGymaeJaei4la8IaeGOmaidaaaGcbaGaemisaG0aaSbaaSqaaiabicdaWaqabaGccqqHPoWvdaWgaaWcbaGaeqyWdihabeaaaaGcdaWadaqaamaalaaabaGagi4CamNaeiyAaKMaeiOBa4gcciGae8x1dO2aaSbaaSqaaiabicdaWaqabaaakeaacqaIXaqmcqGHRaWkcyGGJbWycqGGVbWBcqGGZbWCcqWFvpGAdaWgaaWcbaGaeGimaadabeaaaaGccqGHsisldaWcaaqaaiGbcohaZjabcMgaPjabc6gaUjab=v9aQnaaBaaaleaacqWGLbqzcqWGTbqBaeqaaaGcbaGaeGymaeJaey4kaSIagi4yamMaei4Ba8Maei4CamNae8x1dO2aaSbaaSqaaiabdwgaLjabd2gaTbqabaaaaaGccaGLBbGaayzxaaGaeiilaWcaaa@6ABD@
(47)

By inserting the values ϕ em = ϕ(z) and ϕ0 = ϕ(0), we recover the luminosity distance – redshift relation (44).

4.4 General relativistic (Einstein-de Sitter) limit

The general relativistic limit of the luminosity distance – redshift relation for dust matter and k = 0 = Λ (Einstein-de Sitter model) can be obtained by direct integration of Eq. (27):
d L G R ( z ) = 2 1 + z H 0 Ω ρ 1 / 2 [ 1 + z 1 ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaem4raCKaemOuaifaaOGaeiikaGIaemOEaONaeiykaKIaeyypa0ZaaSaaaeaacqaIYaGmdaGcaaqaaiabigdaXiabgUcaRiabdQha6bWcbeaaaOqaaiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaGGaciab=f8aYbqaaiabigdaXiabc+caViabikdaYaaaaaGccqGGBbWwdaGcaaqaaiabigdaXiabgUcaRiabdQha6bWcbeaakiabgkHiTiabigdaXiabc2faDjabcYcaSaaa@4B86@
(48)

It is straightforward to check that the above result stems out from Eq. (44) by simply switching off the dark radiation.

The general relativistic limit of the luminosity distance – redshift relation should also emerge in the limit Ω λ → 0 of Eq. (42). To see this, we note that when Ω λ → 0, both ϕπ and ϕ0π. Therefore the elliptic integrals of the first and second kind both tend to finite values, thus the differences evaluated at ϕand ϕ0 vanish. Then the only terms which should be carefully investigated are the last two terms of Eq. (42), which are of the type 0/0. By employing Eq. (43), for the last term we obtain:
lim Ω λ 0 Ω λ 1 / 6 sin ϕ e m 1 ε 2 sin 2 ϕ e m 1 + cos ϕ e m = Ω ρ 1 / 6 3 1 / 4 1 + z . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWfqaqaaiGbcYgaSjabcMgaPjabc2gaTbWcbaGaeuyQdC1aaSbaaWqaaGGaciab=T7aSbqabaWccqGHsgIRcqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=T7aSbqaaiabigdaXiabc+caViabiAda2aaakmaalaaabaGagi4CamNaeiyAaKMaeiOBa4Mae8x1dO2aaSbaaSqaaiabdwgaLjabd2gaTbqabaGcdaGcaaqaaiabigdaXiabgkHiTiab=v7aLnaaCaaaleqabaGaeGOmaidaaOGagi4CamNaeiyAaKMaeiOBa42aaWbaaSqabeaacqaIYaGmaaGccqWFvpGAdaWgaaWcbaGaemyzauMaemyBa0gabeaaaeqaaaGcbaGaeGymaeJaey4kaSIagi4yamMaei4Ba8Maei4CamNae8x1dO2aaSbaaSqaaiabdwgaLjabd2gaTbqabaaaaOGaeyypa0ZaaSaaaeaacqqHPoWvdaqhaaWcbaGae8xWdihabaGaeGymaeJaei4la8IaeGOnaydaaaGcbaGaeG4mamZaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaI0aanaaGcdaGcaaqaaiabigdaXiabgUcaRiabdQha6bWcbeaaaaGccqGGUaGlaaa@6FBE@
(49)
Accordingly, the second to last term gives
lim Ω λ 0 Ω λ 1 / 6 sin ϕ 0 1 ε 2 sin 2 ϕ 0 1 + cos ϕ 0 = Ω ρ 1 / 6 3 1 / 4 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWfqaqaaiGbcYgaSjabcMgaPjabc2gaTbWcbaGaeuyQdC1aaSbaaWqaaGGaciab=T7aSbqabaWccqGHsgIRcqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=T7aSbqaaiabigdaXiabc+caViabiAda2aaakmaalaaabaGagi4CamNaeiyAaKMaeiOBa4Mae8x1dO2aaSbaaSqaaiabicdaWaqabaGcdaGcaaqaaiabigdaXiabgkHiTiab=v7aLnaaCaaaleqabaGaeGOmaidaaOGagi4CamNaeiyAaKMaeiOBa42aaWbaaSqabeaacqaIYaGmaaGccqWFvpGAdaWgaaWcbaGaeGimaadabeaaaeqaaaGcbaGaeGymaeJaey4kaSIagi4yamMaei4Ba8Maei4CamNae8x1dO2aaSbaaSqaaiabicdaWaqabaaaaOGaeyypa0ZaaSaaaeaacqqHPoWvdaqhaaWcbaGae8xWdihabaGaeGymaeJaei4la8IaeGOnaydaaaGcbaGaeG4mamZaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaI0aanaaaaaOGaeiOla4caaa@66F2@
(50)

Adding everything together, we recover the general relativistic result (48).

5 Branes with Λ

In this section we discuss certain cases of Randall-Sundrum type brane-worlds with cosmological constant, for which analytical expressions for the luminosity-redshift relation can be found.

5.1 A brane with analytically integrable luminosity distance-redshift relation

If we do not impose the Randall-Sundrum fine-tuning in Eq. (20) and we keep the brane cosmological constant Λ, the polynomial in the denominator of the integrand in Eq. (27) can be simplified for certain values of the dimensionless Ω-s. In particular, if we choose
Ω d = 0 and 4 Ω λ Ω Λ = Ω ρ 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeqadaaabaGaeuyQdC1aaSbaaSqaaiabdsgaKbqabaGccqGH9aqpcqaIWaamaeaacqqGHbqycqqGUbGBcqqGKbazaeaacqaI0aancqqHPoWvdaWgaaWcbaacciGae83UdWgabeaakiabfM6axnaaBaaaleaacqqHBoataeqaaOGaeyypa0JaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabikdaYaaakiabcYcaSaaaaaa@43C5@
(51)
the expression under the square root of denominator becomes a quadratic expression, and the integral can be given in terms of elementary functions [54]:
d L Λ = κ 2 λ / 2 = 2 1 / 3 ( 1 + z ) 6 H 0 Ω ρ 1 / 3 Ω Λ 1 / 6 { ln ( 1 h + h 2 ) [ 1 + h ( 1 + z ) ] 2 ( 1 + h ) 2 [ 1 h ( 1 + z ) h 2 ( 1 + z ) 2 ] + 2 3 arctan 2 h 3 h 2 3 arctan 2 ( 1 + z ) h 3 h } , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeGabaaabaGaemizaq2aa0baaSqaaiabdYeambqaaiabfU5amjabg2da9GGaciab=P7aRnaaCaaameqabaGaeGOmaidaaSGae83UdWMaei4la8IaeGOmaidaaOGaeyypa0ZaaSaaaeaacqaIYaGmdaahaaWcbeqaaiabigdaXiabc+caViabiodaZaaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPaqaaiabiAda2iabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabigdaXiabc+caViabiodaZaaakiabfM6axnaaDaaaleaacqqHBoataeaacqaIXaqmcqGGVaWlcqaI2aGnaaaaaOWaaiqabeaacyGGSbaBcqGGUbGBdaWcaaqaaiabcIcaOiabigdaXiabgkHiTiabdIgaOjabgUcaRiabdIgaOnaaCaaaleqabaGaeGOmaidaaOGaeiykaKIaei4waSLaeGymaeJaey4kaSIaemiAaGMaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeiyxa01aaWbaaSqabeaacqaIYaGmaaaakeaacqGGOaakcqaIXaqmcqGHRaWkcqWGObaAcqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabcUfaBjabigdaXiabgkHiTiabdIgaOjabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPiabgkHiTiabdIgaOnaaCaaaleqabaGaeGOmaidaaOGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKYaaWbaaSqabeaacqaIYaGmaaGccqGGDbqxaaaacaGL7baaaeaadaGaceqaaiabgUcaRiabikdaYmaakaaabaGaeG4mamdaleqaaOGagiyyaeMaeiOCaiNaei4yamMaeiiDaqNaeiyyaeMaeiOBa42aaSaaaeaacqaIYaGmcqGHsislcqWGObaAaeaadaGcaaqaaiabiodaZaWcbeaakiabdIgaObaacqGHsislcqaIYaGmdaGcaaqaaiabiodaZaWcbeaakiGbcggaHjabckhaYjabcogaJjabcsha0jabcggaHjabc6gaUnaalaaabaGaeGOmaiJaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeyOeI0IaemiAaGgabaWaaOaaaeaacqaIZaWmaSqabaGccqWGObaAaaaacaGL9baacqGGSaalaaaaaa@ACE5@
(52)

with h = (Ω ρ /2ΩΛ)1/3.

The first condition (51) merely simplifies the bulk to an anti de Sitter space-time. The second condition (51) by contrast, yields to a much more serious constraint:

κ2λ = 2Λ (53)

The second condition (51), together with the constraint (23) leads to a quadratic equation for Ω λ . For Ω ρ = 0.27 this has two solutions [54]:

ΩΛ = 0.704,   Ω λ = 0.026 (54)

corresponding to the brane tension‡ λ1 = 38.375 × 10-60TeV4 and

ΩΛ = 0.026,   Ω λ = 0.704. (55)

corresponding to the brane tension λ2 = 1.4173 × 10-60TeV4.

It is interesting to note that while solution (55) is ruled out by the recent supernova data, solution (54) is quite close to the present observational value of ΩΛ [55]. From a brane point of view, however the value of the brane tension in the model (54) is far too small, thus it does not describe our physical world. Indeed, all lower limits set for λ are much higher than λ2.

In the two-brane model of Ref. [21] the minimal brane tension depends on the value of the Planck mass M P and on the characteristic curvature scale of the bulk l as λ min = 3 M P 2 / 4 π l 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaaiiGacqWF7oaBdaWgaaWcbaGagiyBa0MaeiyAaKMaeiOBa4gabeaakiabg2da9iabiodaZiabd2eannaaDaaaleaacqWGqbauaeaacqaIYaGmaaGccqGGVaWlcqaI0aancqWFapaCcqWGSbaBdaahaaWcbeqaaiabikdaYaaaaaa@3E42@ [21]. Table-top experiments [5658] on possible deviations from Newton's law currently probe gravity at sub-millimeter scales. As a result they constrain the characteristic curvature scale of the bulk to l ≤ 44 μm. The brane tension therefore (in units c = 1 = ħ) is constrained as λ > 715.887 TeV4. (For a detailed discussion see section 6 of [59], where a slightly lower bound for the brane tension was derived, based on the previously available estimate l ≤ 0.1 mm for the characteristic curvature scale of the bulk.) Big Bang Nucleosynthesis constraints give a much milder lower limit, λ 1 MeV4 [60]. An astrophysical limit λ > 5 × 108 MeV4 (depending on the equation of state of a neutron star) has also been derived [61]. This latter value of λmin is in between the two previous lower limits.

The interpretation of the model (54) is the following. The condition (53) on the models with small brane tension implies Λ ˜ MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacuqHBoatgaacaaaa@2E30@ = 0, thus the bulk becomes flat. As such, it has no effect on dynamics and the fifth dimension becomes superfluous. In fact what we face here is a GR model with stiff fluid scaling as a-6.

5.2 Branes with Ω d 1 and Ω λ 1

In this subsection we assume that both Ω λ and Ω d are small, however we allow for arbitrary values of ΩΛ. These assumptions are motivated by observational evidence that at present our universe is extremely close to a ΛCDM model. A Taylor series expansion of Eq. (27) gives, to leading order in the small parameters:
d L Λ λ d = d L Λ CDM + Ω λ I λ + Ω d I d , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaeu4MdWecciGae83UdWMaemizaqgaaOGaeyypa0Jaemizaq2aa0baaSqaaiabdYeambqaaiabfU5amjabboeadjabbseaejabb2eanbaakiabgUcaRiabfM6axnaaBaaaleaacqWF7oaBaeqaaOGaemysaK0aaSbaaSqaaiab=T7aSbqabaGccqGHRaWkcqqHPoWvdaWgaaWcbaGaemizaqgabeaakiabdMeajnaaBaaaleaacqWGKbazaeqaaOGaeiilaWcaaa@4B05@
(56)
with
d L Λ CDM = a 0 ( 1 + z ) H 0 a e m a 0 d a a 1 / 2 [ Ω Λ a 3 + Ω ρ a 0 3 ] 1 / 2 , I λ = a 0 7 ( 1 + z ) 2 H 0 a e m a 0 d a a 7 / 2 [ Ω Λ a 3 + Ω ρ a 0 3 ] 3 / 2 , I d ( α ) = a 0 5 α ( 1 + z ) 2 H 0 a e m a 0 a α 3 / 2 d a [ Ω Λ a 3 + Ω ρ a 0 3 ] 3 / 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeGabaa8BuaabaqadmaaaeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaeu4MdWKaee4qamKaeeiraqKaeeyta0eaaaGcbaGaeyypa0dabaWaaSaaaeaacqWGHbqydaWgaaWcbaGaeGimaadabeaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPaqaaiabdIeainaaBaaaleaacqaIWaamaeqaaaaakmaapedabaWaaSaaaeaacqWGKbazcqWGHbqyaeaacqWGHbqydaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaakmaadmaabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqWGHbqydaahaaWcbeqaaiabiodaZaaakiabgUcaRiabfM6axnaaBaaaleaaiiGacqWFbpGCaeqaaOGaemyyae2aa0baaSqaaiabicdaWaqaaiabiodaZaaaaOGaay5waiaaw2faamaaCaaaleqabaGaeGymaeJaei4la8IaeGOmaidaaaaakiabcYcaSaWcbaGaemyyae2aaSbaaWqaaiabdwgaLjabd2gaTbqabaaaleaacqWGHbqydaWgaaadbaGaeGimaadabeaaa0Gaey4kIipaaOqaaiabdMeajnaaBaaaleaacqWF7oaBaeqaaaGcbaGaeyypa0dabaWaaSaaaeaacqWGHbqydaqhaaWcbaGaeGimaadabaGaeG4naCdaaOGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKcabaGaeGOmaiJaemisaG0aaSbaaSqaaiabicdaWaqabaaaaOWaa8qmaeaadaWcaaqaaiabdsgaKjabdggaHbqaaiabdggaHnaaCaaaleqabaGaeG4naCJaei4la8IaeGOmaidaaOWaamWaaeaacqqHPoWvdaWgaaWcbaGaeu4MdWeabeaakiabdggaHnaaCaaaleqabaGaeG4mamdaaOGaey4kaSIaeuyQdC1aaSbaaSqaaiab=f8aYbqabaGccqWGHbqydaqhaaWcbaGaeGimaadabaGaeG4mamdaaaGccaGLBbGaayzxaaWaaWbaaSqabeaacqaIZaWmcqGGVaWlcqaIYaGmaaaaaOGaeiilaWcaleaacqWGHbqydaWgaaadbaGaemyzauMaemyBa0gabeaaaSqaaiabdggaHnaaBaaameaacqaIWaamaeqaaaqdcqGHRiI8aaGcbaGaemysaK0aa0baaSqaaiabdsgaKbqaaiabcIcaOiab=f7aHjabcMcaPaaaaOqaaiabg2da9aqaamaalaaabaGaemyyae2aa0baaSqaaiabicdaWaqaaiabiwda1iabgkHiTiab=f7aHbaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPaqaaiabikdaYiabdIeainaaBaaaleaacqaIWaamaeqaaaaakmaapedabaWaaSaaaeaacqWGHbqydaahaaWcbeqaaiab=f7aHjabgkHiTiabiodaZiabc+caViabikdaYaaakiabdsgaKjabdggaHbqaamaadmaabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqWGHbqydaahaaWcbeqaaiabiodaZaaakiabgUcaRiabfM6axnaaBaaaleaacqWFbpGCaeqaaOGaemyyae2aa0baaSqaaiabicdaWaqaaiabiodaZaaaaOGaay5waiaaw2faamaaCaaaleqabaGaeG4mamJaei4la8IaeGOmaidaaaaakiabc6caUaWcbaGaemyyae2aaSbaaWqaaiabdwgaLjabd2gaTbqabaaaleaacqWGHbqydaWgaaadbaGaeGimaadabeaaa0Gaey4kIipaaaaaaa@D090@
(57)

The first expression is the general relativistic luminosity distance – redshift relation in the presence of a cosmological constant (in the ΛCDM model). The next two integrals represent the correction functions scaling the small coefficients Ω λ and Ω d .

All integrands have the same expression Ω Λ a 3 + Ω ρ a 0 3 in the denominator. The roots of this cubic polynomial are:
α = a 0 ( Ω ρ Ω Λ ) 1 / 3 , β = a 0 2 ( Ω ρ Ω Λ ) 1 / 3 ( 1 + i 3 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaacciGae8xSdeMaeyypa0JaeyOeI0Iaemyyae2aaSbaaSqaaiabicdaWaqabaGcdaqadaqaamaalaaabaGaeuyQdC1aaSbaaSqaaiab=f8aYbqabaaakeaacqqHPoWvdaWgaaWcbaGaeu4MdWeabeaaaaaakiaawIcacaGLPaaadaahaaWcbeqaaiabigdaXiabc+caViabiodaZaaakiabcYcaSaqaaiab=j7aIjabg2da9maalaaabaGaemyyae2aaSbaaSqaaiabicdaWaqabaaakeaacqaIYaGmaaWaaeWaaeaadaWcaaqaaiabfM6axnaaBaaaleaacqWFbpGCaeqaaaGcbaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaaaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaIZaWmaaGcdaqadaqaaiabigdaXiabgUcaRiabdMgaPnaakaaabaGaeG4mamdaleqaaaGccaGLOaGaayzkaaaaaaaa@5626@
(58)
and β*. Then d L Λ CDM MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaeu4MdWKaee4qamKaeeiraqKaeeyta0eaaaaa@33FD@ can be rewritten as:
d L Λ CDM = a 0 ( 1 + z ) H 0 Ω Λ 1 / 2 a a 0 d a a 1 / 2 ( a α ) 1 / 2 ( a β ) 1 / 2 ( a β ) 1 / 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaeu4MdWKaee4qamKaeeiraqKaeeyta0eaaOGaeyypa0ZaaSaaaeaacqWGHbqydaWgaaWcbaGaeGimaadabeaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPaqaaiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiabfU5ambqaaiabigdaXiabc+caViabikdaYaaaaaGcdaWdXaqaamaalaaabaGaemizaqMaemyyaegabaGaemyyae2aaWbaaSqabeaacqaIXaqmcqGGVaWlcqaIYaGmaaGccqGGOaakcqWGHbqycqGHsisliiGacqWFXoqycqGGPaqkdaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaakiabcIcaOiabdggaHjabgkHiTiab=j7aIjabcMcaPmaaCaaaleqabaGaeGymaeJaei4la8IaeGOmaidaaOGaeiikaGIaemyyaeMaeyOeI0Iae8NSdi2aaWbaaSqabeaacqGHxiIkaaGccqGGPaqkdaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaaaaGccqGGUaGlaSqaaiabdggaHbqaaiabdggaHnaaBaaameaacqaIWaamaeqaaaqdcqGHRiI8aaaa@6D3F@
(59)
The integration can be carried out by employing Eq. (260.00) of [53] and we obtain the result:
d L Λ CDM ( z ) = ( 1 + z ) [ F ( ϕ 0 , ε ) F ( ϕ , ε ) ] 3 1 / 4 H 0 Ω ρ 1 / 3 Ω Λ 1 / 6 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaqhaaWcbaGaemitaWeabaGaeu4MdWKaee4qamKaeeiraqKaeeyta0eaaOGaeiikaGIaemOEaONaeiykaKIaeyypa0ZaaSaaaeaacqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkcqGGBbWwcqWGgbGrcqGGOaakiiGacqWFvpGAdaWgaaWcbaGaeGimaadabeaakiabcYcaSiab=v7aLjabcMcaPiabgkHiTiabdAeagjabcIcaOiab=v9aQjabcYcaSiab=v7aLjabcMcaPaqaaiabiodaZmaaCaaaleqabaGaeGymaeJaei4la8IaeGinaqdaaOGaemisaG0aaSbaaSqaaiabicdaWaqabaGccqqHPoWvdaqhaaWcbaGae8xWdihabaGaeGymaeJaei4la8IaeG4mamdaaOGaeuyQdC1aa0baaSqaaiabfU5ambqaaiabigdaXiabc+caViabiAda2aaaaaGccqGGSaalaaa@6230@
(60)
with the variable ϕand argument ε of the elliptic integral of the first kind F (ϕ, ε) given by
ε 2 + 1 2 + 3 4 , ϕ = arccos ( 1 3 ) Ω Λ 1 / 3 + ( 1 + z ) Ω ρ 1 / 3 ( 1 + 3 ) Ω Λ 1 / 3 + ( 1 + z ) Ω ρ 1 / 3 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqaaeGabaaabaacciGae8xTdu2aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkdaWcaaqaaiabigdaXaqaaiabikdaYaaacqGHRaWkdaWcaaqaamaakaaabaGaeG4mamdaleqaaaGcbaGaeGinaqdaaiabcYcaSaqaaiab=v9aQjabg2da9iGbcggaHjabckhaYjabcogaJjabcogaJjabc+gaVjabcohaZnaalaaabaGaeiikaGIaeGymaeJaeyOeI0YaaOaaaeaacqaIZaWmaSqabaGccqGGPaqkcqqHPoWvdaqhaaWcbaGaeu4MdWeabaGaeGymaeJaei4la8IaeG4mamdaaOGaey4kaSIaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabigdaXiabc+caViabiodaZaaaaOqaaiabcIcaOiabigdaXiabgUcaRmaakaaabaGaeG4mamdaleqaaOGaeiykaKIaeuyQdC1aa0baaSqaaiabfU5ambqaaiabigdaXiabc+caViabiodaZaaakiabgUcaRiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPiabfM6axnaaDaaaleaacqWFbpGCaeaacqaIXaqmcqGGVaWlcqaIZaWmaaaaaOGaeiOla4caaaaa@6FD4@
(61)

(Note that ε2 is the same as in the case ΩΛ = 0 = Ω d , while ϕis different. Here 0 ≤ ϕπ/2 while for other values of ϕ, we use Eqs. (38).)

It is relatively easy to integrate the contribution of the term linear in Ω λ in terms of the variable t = a3/4. After a partial integration meant to reduce the powers in the denominator we employ
a a 0 d a a 1 / 2 [ a 3 + Ω ρ Ω Λ a 0 3 ] 1 / 2 = Ω Λ 1 / 3 3 1 / 4 a 0 Ω ρ 1 / 3 [ F ( ϕ 0 , ε ) F ( ϕ , ε ) ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWdXaqaamaalaaabaGaemizaqMaemyyaegabaGaemyyae2aaWbaaSqabeaacqaIXaqmcqGGVaWlcqaIYaGmaaGcdaWadaqaaiabdggaHnaaCaaaleqabaGaeG4mamdaaOGaey4kaSYaaSaaaeaacqqHPoWvdaWgaaWcbaacciGae8xWdihabeaaaOqaaiabfM6axnaaBaaaleaacqqHBoataeqaaaaakiabdggaHnaaDaaaleaacqaIWaamaeaacqaIZaWmaaaakiaawUfacaGLDbaadaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaaaaGccqGH9aqpdaWcaaqaaiabfM6axnaaDaaaleaacqqHBoataeaacqaIXaqmcqGGVaWlcqaIZaWmaaaakeaacqaIZaWmdaahaaWcbeqaaiabigdaXiabc+caViabisda0aaakiabdggaHnaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabigdaXiabc+caViabiodaZaaaaaGcdaWadaqaaiabdAeagjabcIcaOiab=v9aQnaaBaaaleaaieaacqGFWaamaeqaaOGaeiilaWIae8xTduMaeiykaKIaeyOeI0IaemOrayKaeiikaGIae8x1dOMaeiilaWIae8xTduMaeiykaKcacaGLBbGaayzxaaaaleaacqWGHbqyaeaacqWGHbqydaWgaaadbaGaeGimaadabeaaa0Gaey4kIipakiabbccaGiabcYcaSaaa@739A@
(62)
and obtain
I λ = 1 + z 15 H 0 Ω ρ 2 { 8 Ω Λ + 3 Ω ρ ( Ω Λ + Ω ρ ) 1 / 2 ( 1 + z ) 8 Ω Λ + 3 Ω ρ ( 1 + z ) 3 [ Ω Λ + Ω ρ ( 1 + z ) 3 ] 1 / 2 } + 8 Ω Λ 5 / 6 ( 1 + z ) 15 3 4 H 0 Ω ρ 7 / 3 [ F ( ϕ 0 , ε ) F ( ϕ , ε ) ] , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqadeGabaaabaGaemysaK0aaSbaaSqaaGGaciab=T7aSbqabaGccqGH9aqpdaWcaaqaaiabigdaXiabgUcaRiabdQha6bqaaiabigdaXiabiwda1iabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabikdaYaaaaaGcdaGadeqaamaalaaabaGaeGioaGJaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqGHRaWkcqaIZaWmcqqHPoWvdaWgaaWcbaGae8xWdihabeaaaOqaaiabcIcaOiabfM6axnaaBaaaleaacqqHBoataeqaaOGaey4kaSIaeuyQdC1aaSbaaSqaaiab=f8aYbqabaGccqGGPaqkdaahaaWcbeqaaiabigdaXiabc+caViabikdaYaaaaaGccqGHsislcqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkdaWcaaqaaiabiIda4iabfM6axnaaBaaaleaacqqHBoataeqaaOGaey4kaSIaeG4mamJaeuyQdC1aaSbaaSqaaiab=f8aYbqabaGccqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkdaahaaWcbeqaaiabiodaZaaaaOqaamaadmaabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqGHRaWkcqqHPoWvdaWgaaWcbaGae8xWdihabeaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPmaaCaaaleqabaGaeG4mamdaaaGccaGLBbGaayzxaaWaaWbaaSqabeaacqaIXaqmcqGGVaWlcqaIYaGmaaaaaaGccaGL7bGaayzFaaaabaGaey4kaSYaaSaaaeaacqaI4aaocqqHPoWvdaqhaaWcbaGaeu4MdWeabaGaeGynauJaei4la8IaeGOnaydaaOGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKcabaGaeGymaeJaeGynauZaaOqaaeaacqaIZaWmaSqaaiabisda0aaakiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabiEda3iabc+caViabiodaZaaaaaGcdaWadaqaaiabdAeagjabcIcaOiab=v9aQnaaBaaaleaacqaIWaamaeqaaOGaeiilaWIae8xTduMaeiykaKIaeyOeI0IaemOrayKaeiikaGIae8x1dOMaeiilaWIae8xTduMaeiykaKcacaGLBbGaayzxaaGaeiilaWcaaaaa@A9D9@
(63)

with the variable ϕ and argument ε given in Eq. (61).

The last term of Eq. (56) is much more complicated to evaluate. For α = 1 and 4 the source term Ω d merely contribute to Ω ρ and ΩΛ, respectively. The more interesting cases are for α = 0, 2, 3. The last term of Eq. (55) for α = 2 consits of elementary functions:
I d ( 2 ) = 1 + z 3 Ω ρ Ω Λ + Ω ρ ( 1 + z ) 3 1 + z 3 Ω ρ Ω Λ + Ω ρ , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaqhaaWcbaGaemizaqgabaGaeiikaGIaeGOmaiJaeiykaKcaaOGaeyypa0ZaaSaaaeaacqaIXaqmcqGHRaWkcqWG6bGEaeaacqaIZaWmcqqHPoWvdaWgaaWcbaacciGae8xWdihabeaakmaakaaabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqGHRaWkcqqHPoWvdaWgaaWcbaGae8xWdihabeaakiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPmaaCaaaleqabaGaeG4mamdaaaqabaaaaOGaeyOeI0YaaSaaaeaacqaIXaqmcqGHRaWkcqWG6bGEaeaacqaIZaWmcqqHPoWvdaWgaaWcbaGae8xWdihabeaakmaakaaabaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaGccqGHRaWkcqqHPoWvdaWgaaWcbaGae8xWdihabeaaaeqaaaaakiabcYcaSaaa@59F6@
(64)
while I d ( 0 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaqhaaWcbaGaemizaqgabaGaeiikaGIaeGimaaJaeiykaKcaaaaa@31E5@ and I d ( 3 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGjbqsdaqhaaWcbaGaemizaqgabaGaeiikaGIaeG4mamJaeiykaKcaaaaa@31EB@ are more complicated to evalute, and we give details of the derivation in the Appendix. By passing to the variable z instead of a, we obtain:
I d ( 0 ) = ( 1 + z ) 3 H 0 Ω ρ 2 [ 4 Ω Λ + 3 Ω ρ Ω Λ + Ω ρ 4 Ω Λ + 3 ( 1 + z ) 3 Ω ρ ( 1 + z ) Ω Λ + ( 1 + z ) 3 Ω ρ ] 8 ( 1 + z ) Ω Λ 1 / 6 3 H 0 Ω ρ 5 / 3 I d ( ϕ , ε ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqadeGabaaabaGaemysaK0aa0baaSqaaiabdsgaKbqaaiabcIcaOiabicdaWiabcMcaPaaakiabg2da9maalaaabaGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKcabaGaeG4mamJaemisaG0aaSbaaSqaaiabicdaWaqabaGccqqHPoWvdaqhaaWcbaacciGae8xWdihabaGaeGOmaidaaaaakmaadmaabaWaaSaaaeaacqaI0aancqqHPoWvdaWgaaWcbaGaeu4MdWeabeaakiabgUcaRiabiodaZiabfM6axnaaBaaaleaacqWFbpGCaeqaaaGcbaWaaOaaaeaacqqHPoWvdaWgaaWcbaGaeu4MdWeabeaakiabgUcaRiabfM6axnaaBaaaleaacqWFbpGCaeqaaaqabaaaaOGaeyOeI0YaaSaaaeaacqaI0aancqqHPoWvdaWgaaWcbaGaeu4MdWeabeaakiabgUcaRiabiodaZiabcIcaOiabigdaXiabgUcaRiabdQha6jabcMcaPmaaCaaaleqabaGaeG4mamdaaOGaeuyQdC1aaSbaaSqaaiab=f8aYbqabaaakeaacqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkdaGcaaqaaiabfM6axnaaBaaaleaacqqHBoataeqaaOGaey4kaSIaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKYaaWbaaSqabeaacqaIZaWmaaGccqqHPoWvdaWgaaWcbaGae8xWdihabeaaaeqaaaaaaOGaay5waiaaw2faaaqaaiabgkHiTmaalaaabaGaeGioaGJaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaeuyQdC1aa0baaSqaaiabfU5ambqaaiabigdaXiabc+caViabiAda2aaaaOqaaiabiodaZiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiab=f8aYbqaaiabiwda1iabc+caViabiodaZaaaaaGccqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabcIcaOiab=v9aQjabcYcaSiab=v7aLjabcMcaPiabc6caUaaaaaa@9682@
(65)
and
I d ( 3 ) = ( 1 + z ) 3 H 0 Ω ρ [ 1 Ω Λ + Ω ρ 1 ( 1 + z ) Ω Λ + ( 1 + z ) 3 Ω ρ ] + 2 ( 1 + z ) 3 H 0 Ω Λ 5 / 6 Ω ρ 2 / 3 I d ( ϕ , ε ) .
(66)
where
I d ( ϕ , ε ) = ( 2 + 3 ) 3 4 ( 1 + 3 ) 3 { F ( ϕ 0 , ε ) F ( ϕ , ε ) ( 3 + 3 ) [ E ( ϕ 0 , ε ) E ( ϕ , ε ) ] + ( 2 + 3 ) ( 3 + 3 ) × [ sin ϕ 0 1 ε 2 sin 2 ϕ 0 1 + ( 2 + 3 ) cos ϕ 0 sin ϕ 1 ε 2 sin 2 ϕ 1 + ( 2 + 3 ) cos ϕ ] } . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqadeWabaaabaGaemysaK0aaSbaaSqaaiabdsgaKbqabaGccqGGOaakiiGacqWFvpGAcqGGSaalcqWF1oqzcqGGPaqkcqGH9aqpdaWcaaqaaiabcIcaOiabikdaYiabgUcaRmaakaaabaGaeG4mamdaleqaaOGaeiykaKcabaWaaOqaaeaacqaIZaWmaSqaaiabisda0aaakiabcIcaOiabigdaXiabgUcaRmaakaaabaGaeG4mamdaleqaaOGaeiykaKYaaWbaaSqabeaacqaIZaWmaaaaaOWaaiqaaeaafaqabeWabaaabaaabaGaemOrayKaeiikaGIae8x1dO2aaSbaaSqaaiabicdaWaqabaGccqGGSaalcqWF1oqzcqGGPaqkcqGHsislcqWGgbGrcqGGOaakcqWFvpGAcqGGSaalcqWF1oqzcqGGPaqkaeaaaaaacaGL7baaaeaacqGHsisldaqadaqaaiabiodaZiabgUcaRmaakaaabaGaeG4mamdaleqaaaGccaGLOaGaayzkaaGaei4waSLaemyrauKaeiikaGIae8x1dO2aaSbaaSqaaiabicdaWaqabaGccqGGSaalcqWF1oqzcqGGPaqkcqGHsislcqWGfbqrcqGGOaakcqWFvpGAcqGGSaalcqWF1oqzcqGGPaqkcqGGDbqxcqGHRaWkdaqadaqaaiabikdaYiabgUcaRmaakaaabaGaeG4mamdaleqaaaGccaGLOaGaayzkaaWaaeWaaeaacqaIZaWmcqGHRaWkdaGcaaqaaiabiodaZaWcbeaaaOGaayjkaiaawMcaaaqaaiabgEna0oaaciaabaWaamWaaeaadaWcaaqaaiGbcohaZjabcMgaPjabc6gaUjab=v9aQnaaBaaaleaacqaIWaamaeqaaOWaaOaaaeaacqaIXaqmcqGHsislcqWF1oqzdaahaaWcbeqaaiabikdaYaaakiGbcohaZjabcMgaPjabc6gaUnaaCaaaleqabaGaeGOmaidaaOGae8x1dO2aaSbaaSqaaiabicdaWaqabaaabeaaaOqaaiabigdaXiabgUcaRiabcIcaOiabikdaYiabgUcaRmaakaaabaGaeG4mamdaleqaaOGaeiykaKIagi4yamMaei4Ba8Maei4CamNae8x1dO2aaSbaaSqaaiabicdaWaqabaaaaOGaeyOeI0YaaSaaaeaacyGGZbWCcqGGPbqAcqGGUbGBcqWFvpGAdaGcaaqaaiabigdaXiabgkHiTiab=v7aLnaaCaaaleqabaGaeGOmaidaaOGagi4CamNaeiyAaKMaeiOBa42aaWbaaSqabeaacqaIYaGmaaGccqWFvpGAaSqabaaakeaacqaIXaqmcqGHRaWkcqGGOaakcqaIYaGmcqGHRaWkdaGcaaqaaiabiodaZaWcbeaakiabcMcaPiGbcogaJjabc+gaVjabcohaZjab=v9aQbaaaiaawUfacaGLDbaaaiaaw2haaiabc6caUaaaaaa@BD14@
(67)

Thus, the analytic expression of the generic luminosity distance – redshift relation on branes with cosmological constant and small values of Ω λ and Ω d is given to first order accuracy in these small parameters by Eqs. (56), (60), (63)–(67).

6 Concluding remarks

The main purpose of this paper was to present the analytical formulation of the luminosity distance – redshift relation in the generalized Randall-Sundrum type II brane-world models containing a Weyl fluid either in the form of dark radiation or as radiation leaving the brane and feeding the bulk black holes. We have given the luminosity distance in terms of elementary functions and elliptical integrals of first and second type and we have also shown how the different limits arise from the generic result. Our results hold for:

(a) Models with Randall-Sundrum fine-tuning (Λ = 0), with or without dark radiation from the bulk and with or without considerable contribution from the energy-momentum squared source terms, discussed in section 4.

(b) The models discussed in subsection 5.1, obeying Λ = κ2/2, integrable in terms of elementary functions and

(c) Models with a brane cosmological constant, discussed to first order accuracy in both the Weyl fluid and energy-momentum squared sources.

This last class of models, presented in subsection 5.2 in the latest times of the cosmological evolution are only slightly different from the Λ CDM model, as they have Ω d 1 and Ω λ 1. The derived modifications in the luminosity distance – redshift formula then represent corrections to the corresponding formula of the ΛCDM model.

While the focus of the present paper is the integrability of the luminosity distance – redshift relation in various brane-world models, in a forthcoming paper [62] we will discuss how well the presently available supernova data support the brane-world models with a small amount of Weyl fluid.

Appendix A. The evaluation of the integral I d

We can integrate the last term of Eq. (55) for α = 0, 3, as follows. First we pass to the variable t = a3/2 and we perform a partial integration in order to reduce the powers in the denominator of the integrand
I d ( 0 ) = ( 1 + z ) a 0 5 3 H 0 Ω Λ 3 / 2 a 3 / 2 a 0 3 / 2 d t t 4 / 3 [ t 2 + Ω ρ a 0 3 Ω Λ ] 3 / 2 = ( 1 + z ) Ω Λ 1 / 2 3 H 0 a 0 Ω ρ 2 [ 4 t 2 + 3 Ω ρ a 0 3 Ω Λ t 1 / 3 t 2 + Ω ρ a 0 3 Ω Λ ] a 3 / 2 a 0 3 / 2 8 ( 1 + z ) Ω Λ 1 / 2 9 H 0 a 0 Ω ρ 2 a 3 / 2 a 0 3 / 2 t 2 / 3 d t t 2 + Ω ρ a 0 3 Ω Λ , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaafaqabeWabaaabaGaemysaK0aa0baaSqaaiabdsgaKbqaaiabcIcaOiabicdaWiabcMcaPaaakiabg2da9iabgkHiTmaalaaabaGaeiikaGIaeGymaeJaey4kaSIaemOEaONaeiykaKIaemyyae2aa0baaSqaaiabicdaWaqaaiabiwda1aaaaOqaaiabiodaZiabdIeainaaBaaaleaacqaIWaamaeqaaOGaeuyQdC1aa0baaSqaaiabfU5ambqaaiabiodaZiabc+caViabikdaYaaaaaGcdaWdXaqaamaalaaabaGaemizaqMaemiDaqhabaGaemiDaq3aaWbaaSqabeaacqaI0aancqGGVaWlcqaIZaWmaaGcdaWadaqaaiabdsha0naaCaaaleqabaGaeGOmaidaaOGaey4kaSYaaSaaaeaacqqHPoWvdaWgaaWcbaacciGae8xWdihabeaakiabdggaHnaaDaaaleaacqaIWaamaeaacqaIZaWmaaaakeaacqqHPoWvdaWgaaWcbaGaeu4MdWeabeaaaaaakiaawUfacaGLDbaadaahaaWcbeqaaiabiodaZiabc+caViabikdaYaaaaaaabaGaemyyae2aaWbaaWqabeaacqaIZaWmcqGGVaWlcqaIYaGmaaaaleaacqWGHbqydaqhaaadbaGaeGimaadabaGaeG4mamJaei4la8IaeGOmaidaaaqdcqGHRiI8aaGcbaGaeyypa0ZaaSaaaeaacqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkcqqHPoWvdaqhaaWcbaGaeu4MdWeabaGaeGymaeJaei4la8IaeGOmaidaaaGcbaGaeG4mamJaemisaG0aaSbaaSqaaiabicdaWaqabaGccqWGHbqydaWgaaWcbaGaeGimaadabeaakiabfM6axnaaDaaaleaacqWFbpGCaeaacqaIYaGmaaaaaOWaamWaaeaadaWcaaqaaiabisda0iabdsha0naaCaaaleqabaGaeGOmaidaaOGaey4kaSIaeG4mamZaaSaaaeaacqqHPoWvdaWgaaWcbaGae8xWdihabeaakiabdggaHnaaDaaaleaacqaIWaamaeaacqaIZaWmaaaakeaacqqHPoWvdaWgaaWcbaGaeu4MdWeabeaaaaaakeaacqWG0baDdaahaaWcbeqaaiabigdaXiabc+caViabiodaZaaakmaakaaabaGaemiDaq3aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkdaWcaaqaaiabfM6axnaaBaaaleaacqWFbpGCaeqaaOGaemyyae2aa0baaSqaaiabicdaWaqaaiabiodaZaaaaOqaaiabfM6axnaaBaaaleaacqqHBoataeqaaaaaaeqaaaaaaOGaay5waiaaw2faamaaDaaaleaacqWGHbqydaahaaadbeqaaiabiodaZiabc+caViabikdaYaaaaSqaaiabdggaHnaaDaaameaacqaIWaamaeaacqaIZaWmcqGGVaWlcqaIYaGmaaaaaaGcbaGaeyOeI0YaaSaaaeaacqaI4aaocqGGOaakcqaIXaqmcqGHRaWkcqWG6bGEcqGGPaqkcqqHPoWvdaqhaaWcbaGaeu4MdWeabaGaeGymaeJaei4la8IaeGOmaidaaaGcbaGaeGyoaKJaemisaG0aaSbaaSqaaiabicdaWaqabaGccqWGHbqydaWgaaWcbaGaeGimaadabeaakiabfM6axnaaDaaaleaacqWFbpGCaeaacqaIYaGmaaaaaOWaa8qmaeaadaWcaaqaaiabdsha0naaCaaaleqabaGaeGOmaiJaei4la8IaeG4mamdaaOGaemizaqMaemiDaqhabaWaaOaaaeaacqWG0baDdaahaaWcbeqaaiabikdaYaaakiabgUcaRmaalaaabaGaeuyQdC1aaSbaaSqaaiab=f8aYbqabaGccqWGHbqydaqhaaWcbaGaeGimaadabaGaeG4mamdaaaGcbaGaeuyQdC1aaSbaaSqaaiabfU5ambqabaaaaaqabaaaaOGaeiilaWcaleaacqWGHbqydaahaaadbeqaaiabiodaZiabc+caViabikdaYaaaaSqaaiabdggaHnaaDaaameaacqaIWaamaeaacqaIZaWmcqGGVaWlcqaIYaGmaaaaniabgUIiYdaaaaaa@E77E@
(A.1)
and
I d ( 3 ) = ( 1 + z ) a 0 2 3 H 0 Ω Λ 3 / 2 a 3 / 2 a 0 3 / 2 t 2 / 3 d t [ t 2 + Ω ρ a 0 3 Ω Λ ] 3 / 2 = ( 1 + z ) Ω Λ 1 / 2 3 H 0 a 0 Ω Λ 1 / 2 Ω ρ [ t 2 t 1 / 3 t 2 + Ω ρ a 0 3 Ω Λ ] a 3 / 2 a 0 3 / 2 + 2 ( 1 + z ) 9 H 0 a 0 Ω Λ 1 / 2 Ω ρ a 3 / 2 a 0 3 / 2 t 2 / 3 d t t 2 + Ω ρ a 0 3 Ω Λ .
(A.2)
By a change of the integration variable to x = Ω Λ 1 / 3 t 3 / 2 / Ω ρ 1 / 3 a 0 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWG4baEcqGH9aqpcqqHPoWvdaqhaaWcbaGaeu4MdWeabaGaeGymaeJaei4la8IaeG4mamdaaOGaemiDaq3aaWbaaSqabeaacqaIZaWmcqGGVaWlcqaIYaGmaaGccqGGVaWlcqqHPoWvdaqhaaWcbaacciGae8xWdihabaGaeGymaeJaei4la8IaeG4mamdaaOGaemyyae2aaSbaaSqaaiabicdaWaqabaaaaa@4344@ we can employ Eq. (260.52) of [53] in order to evaluate the remaining integral:
a 3 / 2 a 0 3 / 2 t 2 / 3 d t t 2 + Ω ρ a 0 3 Ω Λ = 3 a 0 Ω ρ 1 / 3 2 Ω Λ 1 / 3 Ω Λ 1 / 3 a Ω ρ 1 / 3 a 0 Ω Λ 1 / 3 Ω ρ 1 / 3 x 2 d x x 3 + 1 = 3 a 0 Ω ρ 1 / 3 Ω Λ 1 / 3 I d ( ϕ , ε ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeGadiaihiaKjiaZjuaabeqacmaaaeaadaWdXaqaamaalaaabaGaemiDaq3aaWbaaSqabeaacqaIYaGmcqGGVaWlcqaIZaWmaaGccqWGKbazcqWG0baDaeaadaGcaaqaaiabdsha0naaCaaaleqabaGaeGOmaidaaOGaey4kaSYaaSaaaeaacqqHPoWvdaWgaaWcbaacciGae8xWdihabeaakiabdggaHnaaDaaaleaacqaIWaamaeaacqaIZaWmaaaakeaacqqHPoWvdaWgaaWcbaGaeu4MdWeabeaaaaaabeaaaaaabaGaemyyae2aaWbaaWqabeaacqaIZaWmcqGGVaWlcqaIYaGmaaaaleaacqWGHbqydaqhaaadbaGaeGimaadabaGaeG4mamJaei4la8IaeGOmaidaaaqdcqGHRiI8aaGcbaGaeyypa0dabaWaaSaaaeaacqaIZaWmcqWGHbqydaWgaaWcbaGaeGimaadabeaakiabfM6axnaaDaaaleaacqaHbpGCaeaacqaIXaqmcqGGVaWlcqaIZaWmaaaakeaacqaIYaGmcqqHPoWvdaqhaaWcbaGaeu4MdWeabaGaeGymaeJaei4la8IaeG4mamdaaaaakmaapedabaWaaSaaaeaacqWG4baEdaahaaWcbeqaaiabikdaYaaakiabdsgaKjabdIha4bqaamaakaaabaGaemiEaG3aaWbaaSqabeaacqaIZaWmaaGccqGHRaWkcqaIXaqmaSqabaaaaaqaamaalaaabaGaeuyQdC1aa0baaWqaaiabfU5ambqaaiabigdaXiabc+caViabiodaZaaaliabdggaHbqaaiabfM6axnaaDaaameaacqWFbpGCaeaacqaIXaqmcqGGVaWlcqaIZaWmaaWccqWGHbqydaWgaaadbaGaeGimaadabeaaaaaaleaadaWcaaqaaiabfM6axnaaDaaameaacqqHBoataeaacqaIXaqmcqGGVaWlcqaIZaWmaaaaleaacqqHPoWvdaqhaaadbaGae8xWdihabaGaeGymaeJaei4la8IaeG4mamdaaaaaa0Gaey4kIipaaOqaaaqaaiabg2da9aqaamaalaaabaGaeG4mamJaemyyae2aaSbaaSqaaiabicdaWaqabaGccqqHPoWvdaqhaaWcbaGae8xWdihabaGaeGymaeJaei4la8IaeG4mamdaaaGcbaGaeuyQdC1aa0baaSqaaiabfU5ambqaaiabigdaXiabc+caViabiodaZaaaaaGccqWGjbqsdaWgaaWcbaGaemizaqgabeaakiabcIcaOiab=v9aQjabcYcaSiab=v7aLjabcMcaPiabc6caUaaaaaa@A3BA@
(A.3)

Here I d (ϕ, ε) is given by Eqs. (67).

Note

‡ All values of the brane tension given in this subsection are in units c = 1 = ħ.

Declarations

Acknowledgements

This work was supported by OTKA grants no. T046939, 69036 and T042509. 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.

Authors’ Affiliations

(1)
Departments of Theoretical and Experimental Physics, University of Szeged

References

  1. Milgrom M: A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. Astrophys J. 1983, 270: 365-10.1086/161130.View ArticleADSGoogle Scholar
  2. Milgrom M: A modification of the Newtonian dynamics – Implications for galaxies. Astrophys J. 1983, 270: 371-10.1086/161131.View ArticleADSGoogle Scholar
  3. Milgrom M: A Modification of the Newtonian Dynamics – Implications for Galaxy Systems. Astrophys J. 1983, 270: 384-10.1086/161132.View ArticleADSGoogle Scholar
  4. Bekenstein JD: Relativistic gravitation theory for the MOND paradigm. Phys Rev. 2004, D70: 083509-Erratum – 2005 ibid. D 71 069901ADSGoogle Scholar
  5. Mannheim PD: Alternatives to Dark Matter and Dark Energy. Prog Part Nucl Phys. 2006, 56: 340-10.1016/j.ppnp.2005.08.001.View ArticleADSGoogle Scholar
  6. Skordis C, Mota DF, Ferreira PG, Boehm C: Large Scale Structure in Bekenstein's theory of relativistic Modified Newtonian Dynamics. Phys Rev Lett. 2006, 96: 011301-10.1103/PhysRevLett.96.011301.View ArticleADSGoogle Scholar
  7. Slosar A, Melchiorri A, Silk J: Did Boomerang hit MOND?. Phys Rev. 2005, D72: 101301-ADSGoogle Scholar
  8. McGaugh S: Confrontation of MOND Predictions with WMAP First Year Data. Astrophys J. 2004, 611: 26-10.1086/421895.View ArticleADSGoogle Scholar
  9. Corbelli E, Salucci P: Testing MOND with Local Group spiral galaxies. 2006, astro-ph/0610618Google Scholar
  10. Lokas EL, Mamon GA, Prada F: The Draco dwarf in CDM and MOND. EAS Publ Ser. 2006, 20: 113-118. 10.1051/eas:2006056. astro-ph/0508668View ArticleGoogle Scholar
  11. Zhao HS: A review on success and problem of MOND on globular cluster scale. Edited by: Jerjen H, Binggeli B. 2005, 198: astro-ph/0508635, to appear in IAU ColGoogle Scholar
  12. Nusser A, Pointecouteau E: Modeling the formation of galaxy clusters in MOND. Mon Not Roy Astron Soc. 2006, 366: 969-976.View ArticleADSGoogle Scholar
  13. Aguirre A, Schaye J, Quataert E: Problems for MOND in Clusters and the Ly-alpha Forest. Astrophys J. 2001, 561: 550-10.1086/323376.View ArticleADSGoogle Scholar
  14. Amarzguioui M, Elgaroy O, Mota DF, Multamaki T: Cosmological constraints on f(R) gravity theories within the Palatini approach. Astron Astrophys. 2006, 454: 707-10.1051/0004-6361:20064994.View ArticleADSMATHGoogle Scholar
  15. Borowiec A, Godlowski W, Szydlowski M: Accelerated Cosmological Models in Modified Gravity tested by distant Supernovae SNIa data. Phys Rev. 2006, D74: 043502-ADSGoogle Scholar
  16. Mena O, Santiago J, Weller J: Constraining Inverse Curvature Gravity with Supernovae. Phys Rev Lett. 2006, 96: 041103-10.1103/PhysRevLett.96.041103.View ArticleADSGoogle Scholar
  17. De Felice A, Hindmarsh M, Trodden M: Ghosts, Instabilities, and Superluminal Propagation in Modified Gravity Models. JCAP. 2006, 06 (08): 005-View ArticleMathSciNetGoogle Scholar
  18. De Felice A, Mukherjee P, Wang Y: Observational Bounds on Modified Gravity Models. 2007, arxiv:0706.1197Google Scholar
  19. Arkani-Hamed N, Dimopoulos S, Dvali G: The Hierarchy Problem and New Dimensions at a Millimeter. Phys Lett B. 1998, 429: 263-10.1016/S0370-2693(98)00466-3.View ArticleADSGoogle Scholar
  20. Arkani-Hamed N, Dimopoulos S, Dvali G: Phenomenology, Astrophysics and Cosmology of Theories with Sub-Millimeter Dimensions and TeV Scale Quantum Gravity. Phys Rev. 1999, D59: 086004-ADSGoogle Scholar
  21. Randall L, Sundrum R: Large mass hierarchy from a small extra dimension. Phys Rev Lett. 1999, 83: 3370-10.1103/PhysRevLett.83.3370.View ArticleADSMathSciNetMATHGoogle Scholar
  22. Randall L, Sundrum R: An Alternative to Compactification. Phys Rev Lett. 1999, 83: 4690-10.1103/PhysRevLett.83.4690.View ArticleADSMathSciNetMATHGoogle Scholar
  23. Gergely LÁ: Generalized Friedmann branes. Phys Rev. 2003, D68: 124011-ADSMathSciNetGoogle Scholar
  24. Shiromizu T, Maeda K, Sasaki M: The Einstein Equations on the 3-Brane World. Phys Rev. 2000, D62: 024012-ADSMathSciNetGoogle Scholar
  25. Maartens R: Brane-world Gravity. Living Rev Rel. 2004, 7: 1-Google Scholar
  26. Binétruy P, Deffayet C, Ellwanger U, Langlois D: Brane cosmological evolution in a bulk with cosmological constant. Phys Lett B. 2000, 477: 285-10.1016/S0370-2693(00)00204-5.View ArticleADSGoogle Scholar
  27. Bruni M, Germani C, Maartens R: Gravitational Collapse on the Brane: A No-Go Theorem. Phys Rev Lett. 2001, 87: 231302-10.1103/PhysRevLett.87.231302.View ArticleADSMathSciNetGoogle Scholar
  28. Dadhich N, Ghosh SG: Gravitational collapse of null fluid on the brane. Phys Lett B. 2001, 518: 1-10.1016/S0370-2693(01)01057-7.View ArticleADSMathSciNetMATHGoogle Scholar
  29. Dadhich N, Ghost SG: Collapsing sphere on the brane radiates. Phys Lett B. 2002, 538: 233-10.1016/S0370-2693(02)01996-2.View ArticleADSMathSciNetMATHGoogle Scholar
  30. Casadio R, Germani C: Gravitational collapse and black hole evolution: do holographic black holes eventually "anti-evaporate"?. Prog Theor Phys. 2005, 114: 23-10.1143/PTP.114.23.View ArticleADSMathSciNetMATHGoogle Scholar
  31. Gergely LÁ: Black holes and dark energy from gravitational collapse on the brane. 2006, hep-th/0603254Google Scholar
  32. Pal S: Braneworld gravitational collapse from a radiative bulk. Phys Rev. 2006, D74: 124019-ADSGoogle Scholar
  33. Mak MK, Harko T: Can the galactic rotation curves be explained in brane world models?. Phys Rev. 2004, D70: 024010-ADSMathSciNetGoogle Scholar
  34. Pal S, Bharadway S, Kar S: Can extra dimensional effects replace dark matter?. Phys Lett B. 2005, 609: 194-10.1016/j.physletb.2005.01.043.View ArticleADSGoogle Scholar
  35. Harko T, Cheng KS: Galactic metric, dark radiation, dark pressure and gravitational lensing in brane world models. Astrophys J. 2006, 636: 8-10.1086/498141.View ArticleADSGoogle Scholar
  36. Harko T, Cheng KS: The virial theorem and the dynamics of clusters of galaxies in the brane world models. Phys Rev. 2007, D76: 044013-ADSGoogle Scholar
  37. Dvali G, Gabadadze G, Porrati M: 4D Gravity on a Brane in 5D Minkowski Space. Phys Lett B. 2000, 485: 208-10.1016/S0370-2693(00)00669-9.View ArticleADSMathSciNetMATHGoogle Scholar
  38. Koyama K: Are there ghosts in the self-accelerating brane universe?. Phys Rev. 2005, D72: 123511-ADSGoogle Scholar
  39. Gorbunov D, Koyama K, Sibiryakov S: More on ghosts in DGP model. Phys Rev. 2006, D73: 044016-ADSGoogle Scholar
  40. Charmousis C, Gregory R, Kaloper N, Padilla A: DGP Specteroscopy. 2006, hep-th/0604086Google Scholar
  41. Izumi K, Koyama K, Tanaka T: Unexorcized ghost in DGP brane world. 2006, hep-th/0610282Google Scholar
  42. Sahni V, Shtanov Y: Braneworld models of dark energy. JCAP. 2003, 03 (11): 014-View ArticleMathSciNetGoogle Scholar
  43. Maeda K, Mizuno S, Torii T: Effective Gravitational Equations on Brane World with Induced Gravity. Phys Rev. 2003, D68: 024033-ADSMathSciNetGoogle Scholar
  44. Gergely LÁ, Maartens R: Asymmetric brane-worlds with induced gravity. Phys Rev. 2005, D71: 024032-ADSMathSciNetGoogle Scholar
  45. Lazkoz R, Maartens R, Majerotto E: Observational constraints on phantom-like braneworld cosmologies. Phys Rev. 2006, D74: 083510-ADSMathSciNetGoogle Scholar
  46. Dabrowski MP, Stachowiak T: Phantom Friedmann cosmologies and higher-order characteristics of expansion. Annals of Physics. 2006, 321: 771-10.1016/j.aop.2005.10.006.View ArticleADSMathSciNetGoogle Scholar
  47. Padmanabhan T: Theoretical Astrophysics. 2002, Cambridge Univ. Press, Cambridge UK, 3: 168-171.View ArticleGoogle Scholar
  48. Doroshkevich A, Tucker DL, Allam S, Way MJ: Large scale structure in the SDSS galaxy survey. Astronomy and Astrophysics. 2004, 418: 7-10.1051/0004-6361:20031780.View ArticleADSMATHGoogle Scholar
  49. Eisenstein DJ, Zehavi I, Hogg DW, et al: Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS Luminous Red Galaxies. Astrophys J. 2005, 633: 560-10.1086/466512.View ArticleADSGoogle Scholar
  50. Tegmark M, Strauss MA, Blanton MR, et al: Cosmological parameters from SDSS and WMAP. Phys Rev. 2004, D69: 103501-ADSGoogle Scholar
  51. Pal S: Structure formation on the brane: A mimicry. Phys Rev. 2006, D74: 024005-ADSGoogle Scholar
  52. Bronstein IN, Semendyayev KA, Musiol G, Mühlig H: Handbook of Mathematics. 2004, Springer-Verlag Berlin Heidelberg New York, 4View ArticleGoogle Scholar
  53. Byrd PF, Friedman MD: Handbook of Elliptic Integrals for Engineers and Scientists. 1971, Springer-Verlag Berlin Heidelberg New YorkView ArticleGoogle Scholar
  54. Nagy B, Keresztes Z: On the luminosity-redshift relation in brane-worlds with cosmological constant. 2006, Publications of the Astronomy Department of the Eötvös University (PADEU), 17: 221-227. astro-ph/0606662Google Scholar
  55. Liddle A: An Introduction to Modern Cosmology. 2003, John Wiley & SonsGoogle Scholar
  56. Long JC, et al: New experimental limits on macroscopic forces below 100 microns. Nature. 2003, 421: 922-10.1038/nature01432.View ArticleADSGoogle Scholar
  57. Gundlach JH, Schlamminger S, Spitzer CD, et al: Laboratory Test of Newton's Second Law for Small Acceleration. Phys Rev Lett. 2007, 98: 150801-10.1103/PhysRevLett.98.150801.View ArticleADSGoogle Scholar
  58. Kapner DJ, Cook TS, Adelberger EG: Tests of the Gravitational Inverse-Square Law below the Dark-Energy Lenght Scale. Phys Rev Lett. 2007, 98: 021101-10.1103/PhysRevLett.98.021101.View ArticleADSGoogle Scholar
  59. Gergely LÁ, Keresztes Z: Irradiated asymmetric Friedmann branes. JCAP. 2006, 06 (01): 022-View ArticleMathSciNetGoogle Scholar
  60. Maartens R, Wands D, Bassett BA, Heard IPC: Chaotic Inflation on the brane. Phys Rev. 2000, D62: 041301 (R)-ADSGoogle Scholar
  61. Germani C, Maartens R: Stars in the braneworld. Phys Rev. 2001, D64: 124010-ADSMathSciNetGoogle Scholar
  62. Szabó Gy M, Gergely LÁ, Keresztes Z: The luminosity-redshift relation in brane-worlds: II. Confrontation with experimental data. 2007, astro-ph/0702610Google Scholar

Copyright

© Keresztes et al. 2007

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.