arxiv: v2 [astro-ph.co] 24 Feb 2012

Transcription

1 BA-TH/65-, CERN-PH-TH/-8, LPTENS-/46 Backreaction on the luminoity-redhift relation from gauge invariant light-cone averaging arxiv:.47v atro-ph.co 4 Feb Ido Ben-Dayan,, Maurizio Gaperini 3,4, Giovanni Marozzi 5, Fabien Nugier 6 and Gabriele Veneziano 5,7 Canadian Intitute for Theoretical Atrophyic, 6 St George, Toronto ON, M5S 3H8 Perimeter Intitute for Theoretical Phyic, Waterloo, Ontario NL Y5, Canada 3 Dipartimento di Fiica, Univerità di Bari, Via G. Amendola 73, 76 Bari, Italy 4 Itituto Nazionale di Fiica Nucleare, Sezione di Bari, Bari, Italy 5 Collège de France, Place M. Berthelot, 755 Pari, France 6 Laboratoire de Phyique Théorique de l École Normale Supérieure, CNRS UMR 8549, 4 Rue Lhomond, 755 Pari, France 7 CERN, Theory Unit, Phyic Department, CH- Geneva 3, Switzerland Uing a recently propoed gauge invariant formulation of light-cone averaging, together with adapted geodeic light-cone coordinate, we how how an induced backreaction effect emerge, in general, from correlated fluctuation in the luminoity ditance and covariant integration meaure. Conidering a realitic tochatic pectrum of inhomogeneitie of primordial inflationary) origin we find that both the induced backreaction on the luminoity-redhift relation and the diperion are larger than naïvely expected. On the other hand the former, at leat to leading order and in the linear perturbative regime, cannot account by itelf for the oberved effect of dark energy at largeredhift. A full econd-order calculation, or even better a reliable etimate of contribution from the non-linear regime, appear to be neceary before firm concluion on the correct interpretation of the data can be drawn. PACS number: 98.8-k, x, 98.8.E I. INTRODUCTION The o-called concordance or ΛCDM) model, baed on a uitable combination of dark matter, dark energy and baryon for an overall critical denity, ha become the reference paradigm for the late i.e. pot-equality epoch evolution of our Univere ee e.g. ). It account equally well for the CMB data, the Large Scale Structure and, even more ignificantly, for the upernovae data in term of a comic acceleration. Strictly peaking thee three tet of the concordance model are not at the ame level of theoretical rigor. While the firt two have to do, by definition, with the inhomogeneitie preent in our Univere, the third i baed on an ideal homogeneou and iotropic Friedmann-Lemaître-Roberton-Walker FLRW) geometry. It i clear that a better treatment of comic acceleration hould take inhomogeneitie into account, at leat in an average tatitical ene. Only when thi i done we can etablih in a convincing way whether ΛCDM give a imultaneou conitent decription of the above-mentioned body of comological data. Thi realization ha led to a vat literature about averaging comological obervable in realitic inhomogeneou comologie ee e.g. 3 for recent review). The concluion, however, are till rather controverial: according to ome author 4 preent inhomogeneitie might explain, by themelve, comic acceleration without any need for dark-energy contribution; according to other 5 the effect of inhomogeneitie i, intead, completely negligible. The truth may lie omewhere in between, in the ene that a quantitative undertanding of inhomogeneitie effect could be important in order to put precie contraint on dark-energy parameter, uch a the critical fraction of dark-energy denity, Ω Λ, and the time evolution of it effective equation of tate, w Λ z). In the firt paper tudying the dynamical effect of averaging, the problem wa approached mainly following Buchert precription 6, namely averaging inhomogeneitie over pacelike hyperurface and computing the enuing backreaction on the averaged geometry. Nonethele it i clear that a proper treatment of comic acceleration which i indeed revealed through the experimental tudy of the luminoity-ditance to redhift relation need to conider the backreaction of averaged inhomogeneitie along our pat light cone, which i a null three-dimenional hyperurface. While a number of paper have emphaized the importance of light-cone averaging ee e.g. 7 9), a fully covariant

2 and gauge invariant formulation of uch a procedure wa only given recently, by generalizing to null hyperurface an analogou precription previouly derived for pacelike hyperurface,. Actually, it turn out that the phyically meaningful covariant and gauge invariant) average of a calar over a null hyperurface reduce to averaging over an appropriate two-dimenional urface a calar object which i non local, a the integrand itelf contain integral along lightlike geodeic curve lying on the given null hyperurface. Nonethele, we hall imply refer to thi procedure a light-cone averaging. The aim of thi paper i to apply uch an averaging procedure to the luminoity ditance of a light ource lying on our pat light-cone, in order to compute the poible backreaction of inhomogeneitie on the luminoity-redhift relation. We will how the emergence of an effect called hereafter induced backreaction which arie from a generic correlation between the inhomogeneitie preent in the variable we want to average e.g. the luminoity ditance) and thoe appearing in the covariant integration meaure. We tre immediately that our induced backreaction account only for a part of the total effect of inhomogeneitie. A we hall dicu in detail below, the leading-order induced backreaction can be computed in term of linear perturbation theory, while a complete leading-order calculation left to future work) would require perturbation theory up to econd order. We will ue a imple phenomenological model of inhomogeneou geometry baed on a patially flat FLRW metric which include, to firt order, calar perturbation of primordial inflationary) origin. Such perturbation are often conveniently parametrized in the longitudinal or Newtonian) gauge 3. However, we will take advantage of the gauge invariance of our formalim to ue an adapted ytem of coordinate defining the o-called geodeic light-cone GLC) gauge introduced in. In that gauge the light-cone average of luminoity ditance and redhift take very imple, exact expreion that keep all the required degree of freedom for being applicable to general geometrie. By further taking the enemble average of the inhomogeneou term, and uing the tochatic propertie of our model of perturbation, we will compute to leading order the induced backreaction on the luminoity ditance, and the correponding variance. If we limit ourelve to firt order computation in the range of cale where linear perturbation theory i reliable, we find that the induced econd-order) correction to the luminoity-redhift relation of the homogeneou CDM model, although much bigger than one could have naïvely expected, are not large enough to mimic a izable fraction of Ω Λ. However, our formalim ugget that other econd-order correction could be even larger, thu confirming the importance of performing a full econd-order calculation or, even better, a reliable calculation in the non-linear hort cale) regime. The paper i organized a follow. In Sect. II we recall definition and baic propertie of GLC coordinate and give exact, non-perturbative expreion for the redhift, the luminoity ditance and the light-cone average of the latter at contant redhift. In Sect. III we give the linear order tranformation expreing the metric of a perturbed FLRW geometry in GLC form, tarting from the more commonly ued Newtonian gauge. We then expre the luminoity ditance a a function of the redhift and of the angular poition of the ource, to firt order in the given perturbed geometry. In Sect. IV we take the light-cone average of the computed luminoity ditance and we how how backreaction effect automatically emerge, in general, from correlation between the averaged variable and the covariant integration meaure. We explicitly define uch an induced backreaction and the comic variance, to leading order, taking into account the tochatic propertie of our model of perturbation and thu including the appropriate enemble average of the inhomogeneou correction. In Sect. V we implement analytical calculation of the induced backreaction term and of the diperion, and in Sect. VI we preent the correponding numerical reult for a pecific pectrum of primordial calar perturbation and for a given model of tranfer function. Our main concluion are finally ummarized in Sect. VII. We alo preent, in Appendix A, a ample of the analytic manipulation needed for reducing our light-cone plu tochatic) average to the form of a one-dimenional integral over the wavenumber k labelling the calar perturbation pectrum. II. THE GLC GAUGE AND THE AVERAGED LUMINOSITY-REDSHIFT RELATION A. Reminder of the GLC gauge and of it main propertie Our purpoe, in thi paper, i to compute the light-cone average of the luminoity ditance a a function of the redhift by applying the gauge invariant procedure introduced in. One of the main virtue of uing a gauge invariant formalim i the freedom of chooing a gauge particularly adapted to the problem at hand. In the cae of pacelike averaging, for intance, a convenient coordinate ytem correpond to a gauge where the averaging A dicued in, thi integration meaure will alo induce backreaction term of the type uually dicued in the literature 6, namely term that arie from generalized) commutation rule between differential operator and averaging integral.

3 hyperurface are identified with contant-time hyperurface. In many application the averaging hyperurface are choen indeed a the one aociated to a cla of geodeic oberver correponding to contant value of the time parameter t of the ynchronou gauge ee e.g. 4). Similarly, for light-cone average, it i convenient to identify the null hyperurface with thoe on which a null coordinate take contant value. For thi reaon we have introduced in an adapted ytem of coordinate defining what we have called a geodeic light-cone GLC) gauge where the averaging precription greatly implifie, while keeping all the required degree of freedom for application to general geometrie. Furthermore we are alo able to identify the timelike coordinate of the GLC gauge with the comic time t of the ynchronou gauge. A a conequence we can eaily introduce a family of geodeic reference oberver which exactly coincide with the tatic one of the ynchronou gauge. Let u firt recall, a dicued in, that the coordinate x µ τ, w, θ a ), a,, pecifying the metric in the GLC gauge correpond to a complete gauge fixing of the o-called obervational coordinate, defined e.g in 7, 5, 6. The GLC metric depend on ix arbitrary function Υ, a two-dimenional vector U a and a ymmetric matrix γ ab ), and it line element take the form d GLC Υ dw Υdwdτ + γ ab d θ a U a dw)d θ b U b dw)..) In matrix form, the metric and it invere are then given by: gµν GLC Υ Υ U b /Υ Υ Υ + U U b, g µν T Ua T GLC Υ γ ab U a ) T /Υ T γ ab 3,.) where, ), U b U, U ), while the matrice γ ab and γ ab lower and raie the two-dimenional indice. Clearly w i a null coordinate i.e. µ w µ w ), and a pat light-cone hyperurface i pecified by the condition w cont. We can alo eaily check that µ τ define a geodeic flow, i.e. that ν τ) ν µ τ) a a conequence of the relation g ττ ). In the limiting cae of a patially flat FLRW geometry, with cale factor a, comic time t, and conformal time parameter η uch that dη dt/a, the tranformation to the GLC coordinate and the meaning of the new metric component are eaily found a follow : τ t, w r + η, Υ at), U a, γ ab dθ a dθ b a t)r dθ + in θdφ )..3) Even though we will be mainly uing the GLC gauge for a perturbed FLRW metric in the Newtonian gauge, it i important to tre that the equality between the coordinate τ and the proper time t of the ynchronou gauge hold at the exact, non perturbative level: it i alway poible, in fact, to chooe the GLC coordinate in uch a way that τ and t are identified like in the above FLRW limit. In order to illutrate thi point let u conider an arbitrary pace-time metric written in the ynchronou gauge with coordinate X µ t, X i ), i, j,, 3), where the line element take the form: d SG dt + h ij dx i dx j..4) Let u impoe the condition t τ, and check whether we run into any contradiction with the exact metric tranformation g ρσ Xρ X σ SG X) x µ x ν gµν GLC x)..5) Uing Eq..) for the GLC metric, and conidering the tranformation for the g tµ condition SG component, we then obtain the g tµ SG {, } τ + Υ w + U a a ) X µ u ν ν X µ dxµ dλ,.6) where u µ µ τ i the four-velocity of the geodeic GLC oberver, and where λ denote an affine parameter along the oberver world-line. So the requirement τ t boil down to the tatement that along the geodeic flow of the vector field u µ the SG coordinate X i are contant. Thi clearly define the coordinate tranformation in a non-perturbative way, and alo how that the geodeic oberver u µ µ τ of the GLC gauge correpond to a tatic and geodeic) oberver in the ynchronou gauge. It follow that the identification t τ can alway be taken for any pace-time

4 metric, and that thi imple connection between GLC and ynchronou gauge ha validity far beyond the particular FLRW cae or it perturbed generalization. We alo remark that, in GLC coordinate, the null geodeic connecting ource and oberver are characterized by the imple tangent vector k µ g µν ν w g µw δ µ τ Υ, which mean that photon reach the oberver travelling at contant w and θ a. Thi make the calculation of the area ditance and of the redhift particularly eay. Conider, for intance, a light ray emitted by a tatic geodeic ource at the interection between the pat light-cone of our oberver, w w, and the patial hyperurface τ τ, and received by uch tatic geodeic oberver at τ τ > τ. The aociated redhift z i then given by : + z ) kµ u µ ) k µ u µ ) o µ w µ τ) µ w µ τ) o Υw, τ, θ a ) Υw, τ, θ a ),.7) where the ubcript o and denote, repectively, a quantity evaluated at the oberver and ource pace-time poition. The expreion for the angular ditance will be explicitly derived in Sect. II B. Let u finally recall that, in GLC coordinate, the covariant average of a calar quantity Sτ, w, θ a ) over the compact two-dimenional urface Σ, defined by the interection of our pat light-cone w w with the pacelike hyperurface τ τ, i imply given by : S w,τ Σ d4 x g δw w )δτ τ )Sτ, w, θ a ) µ τ µ w Σ d4 x g δw w )δτ τ ) µ τ µ w d θ γw, τ, θ a ) Sw, τ, θ a ) d θ γw, τ, θ,.8) a ) where γ det γ ab. In the cae of interet for thi paper, namely light-cone average on urface of contant redhift z z, one then obtain d θ γw, τz, w, θ a ), θ b ) Sw, τz, w, θ a ), θ b ) S w,z d θ γw, τz, w, θ a ), θ,.9) b ) where τz, w, θ a ) ha to be determined by olving the redhift equation.7) for τ a a function of w, τ, z and θ a. Thi general reult will now be applied to the cae in which S i identified with the luminoity ditance d L. 4 B. Light-cone average of the luminoity ditance Let u firt recall that the luminoity ditance d L of a ource at redhift z i related in general to the angular ditance d A of the ource a een from the oberver) by the o-called Etherington or reciprocity) law 7: d L + z) d A..) In the particular cae of an unperturbed, patially flat FLRW background, and for a ource with redhift z, the ditance d A i imply given by d F LRW A z ) a r a η η ),.) where a aη ), while η η denote the conformal time interval between the emiion and obervation of the light ignal. For the unperturbed metric, on the other hand, we have + z a /at), and dη dt/a a dz/h, where H dln a)/dt. Hence: η z d F L LRW dz z ) + z )a dη + z ) η Hz) + z / z dz Ω n + z) 3+wn)..) H n In the lat equality we have ued the tandard patially flat) Friedmann equation for H ee e.g. 8), auming that the given homogeneou model ha perfect fluid ource with preent fraction of the critical denity Ω n and barotropic parameter w n. Expanding in the limit z we alo obtain the expreion d F L LRW z ) z + 3 ) w n Ω n z + Oz 3 H 4 ) H z + q )z + Oz) 3,.3) n

5 which how the well known enitivity of the term quadratic in z to the compoition of the comic fluid through the deceleration parameter q. Let u now dicu how thi well known reult for d L i modified when including generic inhomogeneitie. We recall, to thi purpoe, that in a generic metric background the angular ditance d A can be computed by conidering the null vector k µ dx µ /dλ tangent to the null ray connecting ource and oberver i.e. belonging to the congruence of null geodeic forming the oberver pat light-cone). Here λ i an affine parameter along the ray trajectory, choen in uch a way that λ at the oberver and λ λ at the ource poition. The expanion Θ of the congruence of null ray i then given by Θ µ k µ, and the correponding angular ditance d A i defined by the differential equation 5, 9: d dλ ln d A) Θ µk µ..4) Conider now a generic metric in the GLC gauge, and the null vector k µ µ w uch that k µ dx µ /dλ δ µ τ Υ. We have, in thi cae, where γ det γ ab, and the integration of Eq..4) give, in general, µ k µ Υ τ ln γ) d dλ ln γ),.5) d Aλ) c γλ),.6) where c i independent of λ. In order to fix the contant c we may recall the boundary condition required to guarantee that the null ray generated by k µ belong to the pat light cone centered on the world line of our oberver. In the limit d A or λ ), where θ a θ a ee in particular next ection, Eq. 3.8)), uch condition require in particular that 5 γ in θo in θ,.7) lim λ d A where, in the lat equality, we have ued the already mentioned fact that θ a i contant along the null geodeic. Thi clearly fixe c / in θ, and uniquely determine the angular ditance a d A λ) γ /4 λ) in θ ) /..8) Inerting thi reult into Eq..9) we finally arrive at an exact expreion for the light-cone average of the luminoity ditance, a a function of z, in the GLC gauge: d θ γ d L w,z + z ) / w, τz, w, θ a ), θ b )d A w, τz, w, θ a ), θ b ) d θ γ / w, τz, w, θ a ), θ..9) b ) 5 Thi equation, being exact, can be applied in principle to any given highly inhomogeneou comology. Specifically, it can be applied to a Lemaître-Tolman-Bondi LTB) model even in the cae of a generic oberver hifted away from the ymmetry centre of the iotropic geometry. In thi paper, however, we hall be mainly intereted in working out the conequence of light-cone averaging for a particular model of perturbed FLRW Univere. Since we know how to decribe the latter, for example, in the Newtonian longitudinal) gauge, we need in general to connect uch a coordinate ytem to the GLC gauge. III. LUMINOSITY DISTANCE IN A PERTURBED FLRW GEOMETRY A. Firt-order coordinate tranformation from the Newtonian to the GLC gauge Let u conider the calar perturbation of a conformally flat FLRW background to decribe the particular model of inhomogeneou geometry we are intereted in thi paper. Chooing in particular the o-called longitudinal or conformally) Newtonian gauge NG), uing pherical coordinate r, θ a ) r, θ, φ), and going up to the firt order in perturbation theory, it i well known that the model i parametrized by the following invere) metric tenor 3: g µν NG a η) diag + Φ, + Ψ, + Ψ)γ ab ). 3.)

6 6 Here γ ab diag r, r in θ ), 3.) and Φ, Ψ are the uual gauge invariant Bardeen potential 3, general function of η, r, θ and φ. We will aume, for implicity, that the matter ource have vanihing or negligible) aniotropic tre, o that Φ Ψ. For the ubequent computation in particular, for the application of the imple average precription given previouly we need to re-expre thi metric in GLC form, finding the tranformation from the NG coordinate y µ η, r, θ, φ) to the GLC coordinate x ν τ, w, θ, θ ), and computing the reparametrized metric a g ρσ xρ x σ GLC x) y µ y ν gµν NG y). 3.3) To thi purpoe we will introduce the ueful zeroth-order) light-cone variable η ± η ± r, uch that with correponding partial derivative η η + + η ), r η + η ), 3.4) η + +, r +, ± η ± η ± r ). 3.5) Uing thee variable we olve the three differential equation obtained from Eq. 3.3) for the component gglc ττ, gglc ww, gwa GLC, by impoing the boundary condition that i) the tranformation i non ingular around the oberver poition at r, and ii) that the two-dimenional patial ection r cont are locally parametrized at the oberver poition by tandard pherical coordinate, i.e. θ a ) θ a θ, φ). The ought for tranformation can then be written, to firt order in Ψ, a follow, τ η w η + + θ a θ a + η in dη aη ) + Ψη, r, θ a ), 3.6) η η + dx ˆΨη +, x, θ a ), 3.7) η η + dx ˆγ ab η +, x, θ a ) x η + dy b ˆΨη+, y, θ a ), 3.8) where ˆΨη +, η, θ a ) Ψη, r, θ a ), ˆγ ab η +, η, θ a ) γ ab η, r, θ a ) and η in repreent an early enough time when the perturbation or better the integrand) wa negligible. We can eaily check that, to zeroth order in Ψ, we recover the homogeneou tranformation.3) a expected. To firt order in Ψ we can ue again Eq. 3.3) to compute the non-trivial entrie of the GLC metric.), and obtain: Υ aη) U a ˆγab + η γ ab aη) + ˆΨη +, η +, θ a ) η η η + dx + ˆΨη+, x, θ a ) dx b ˆΨη+, x, θ a ) η + aη) γab dx + η + { ˆγ ab η +, x, θ a ) x + Ψη, r, θ a ) γ ab + η η + dη aη ) r Ψη, r, θ a ); η in 3.9) dη aη ) b Ψη, r, θ a ) η in dy b ˆΨη+, y, θ a ) η + lim x η + ˆγ ac η η + dx c ˆγ bd η +, x, θ a ) ˆγ ab η +, x, θ a ) x x η + dy d ˆΨη+, y, θ a ) dy b ˆΨη+, y, θ a ) ; 3.) η + ) } + a b.3.) The term ˆΨη +, η +, θ a ), appearing in the general equation for Υ, denote at any η) the value of the perturbation potential evaluated at the tip of the light-cone connecting the origin r ) to the point y µ η, r, θ a ) namely, ˆΨη +, η +, θ a ) Ψη + r,, θ a )).

7 The above tranformation can be immediately applied to obtain an explicit expreion for the redhift parameter z. By inerting the reult 3.9) into Eq..7) and conidering that the ource i located on the pat light-cone of the oberver identified by the equation w w η ), we obtain: where J I + I r, and where: + z aη ) + Jz, θ a ), 3.) aη ) 7 I + I r η η+ η η dx + ˆΨη +, x, θ a ) Ψ Ψ o η in dη aη ) aη ) rψη, η η, θ a ) η η dη r Ψη, η η, θ a ), 3.3) η in dη aη ) aη ) rψη,, θ a ). 3.4) We have defined η ± η ± r, Ψ Ψη, η η, θ a ), Ψ o Ψη,, θ a ) and we have ued the zeroth-order light-cone condition η + η + η. It hould be treed, however, that while the integral appearing in I r are evaluated at contant r namely along timelike geodeic), all the other integral are evaluated at fixed η + i.e. along null geodeic on the oberver pat light-cone). The contribution aociated to I r can alo be rewritten a I r v v o ) ˆn, 3.5) where ˆn i the unit tangent vector along the null geodeic connecting ource and oberver, and where v,o η,o η in dη aη ) aη, ) Ψη, r, θ a ) 3.6) are the peculiar velocitie of ource and oberver aociated to a geodeic configuration perturbed up to firt order in the NG gauge. In a realitic ituation one hould add to Eq.3.) imilar term taking into account a poible intrinic non-perturbative) motion of ource and oberver, unle our theoretical prediction are compared with data already corrected for thee latter Doppler contribution. Note alo that our firt-order expreion for the redhift 3.), valid in general for any given cale factor aη), i in full agreement with the expreion obtained in Eq.38) of, where z i computed for the particular cae of a CDM-dominated model. B. The luminoity ditance including firt-order calar perturbation We now apply the above coordinate tranformation to find, in the Newtonian gauge and to firt order in perturbation theory, the other relevant quantitie for thi paper. Let u tart with the determinant γ appearing in the angular ditance.8). For a ource emitting light at time η and radial ditance r we obtain from Eq. 3.), to firt order, γ λ ) det γ ab λ ) a r in θ) + 4Ψ + 4 J z, θ a ), 3.7) where: J 4 η η + x dx ˆγ ab η+, x, θ a ) dy a b ˆΨη +, y, θ a ) η+ η η the latter equality follow upon a imple integration by part). Hence, from.8): η η dη η η η η θ + in θ) φ Ψη, η η, θ a ) 3.8) d A λ ) a r Ψ J in θ ) /. 3.9) in θ Note that, to firt order, we have γ γ γ ), and that, for thee diagonal matrix element, the operator b commute with γ db.

8 The lat factor can be eaily computed, to firt order, by uing Eq. 3.8) and the fact that the θ a are contant along the null geodeic. It i eay to check that it amount to a redefinition of J, and that the above angular ditance become where: J η η η η d A λ ) a r Ψ J z, θ a ), 3.) dη η η θ + cot θ θ + in θ) φ Ψη, η η, θ a ) η η η η η η 8 dη η η η η Ψ 3.) here i the two-dimenional Laplacian operator on the unit two-phere). For the full explicit expreion of the luminoity ditance d L at contant redhift what we need, at thi point, i the firt-order expanion of the factor a r aη )r appearing in Eq. 3.). To thi purpoe we tart from Eq. 3.), conidering z a a contant parameter localiing the given light ource on the pat light-cone w w ) of our oberver, and we look for approximate olution for η η z, θ a ). Let u firt define the zero-order olution η ) through the exact relation aη ) ) a + z, 3.) where a aη ). Expanding 3.) with repect to the parameter δη η η ) we then find: aη) + H δη Jz, θ a ) + H δη Jz, θ a ), 3.3) + z a + z ) where H dln aη ) ))/dη ), o that: H δη Jz, θ a ). 3.4) On the other hand, by applying Eq. 3.7) to the light-cone w w at the ource poition, we readily obtain w η + ηψ av η, 3.5) where we have introduced the zero-order quantity η η η ), and denoted by Ψ av the average value of Ψ along the unperturbed) null geodeic connecting ource and oberver: η η + dx ˆΨη +, x, θ a ) η η dη Ψη, η η, θ a ) ηψ av. 3.6) Combining thi reult with Eq. 3.4) we can then determine, to firt order, the value of the radial coordinate r z, θ a ) correponding to the given redhift z : r z, θ a ) w η ) z ) δη + ηψ av η Jz, θ a ) H η + Ψ av. 3.7) Proceeding in the ame way for a we obtain, from Eq. 3.), 3.), a z, θ a ) aη ) ) + Jz, θ a ), 3.8) which, in turn, allow u to compute the value of a r on the contant-z -urface: a r z, θ a ) aη ) ) η + Ψ av + ) Jz, θ a ). 3.9) H η The angular ditance 3.9), for a ource at redhift z, can now be written a 3 d A z, θ a ) aη ) ) η + Ψ av + ) J Ψη, η η, θ a ) J. 3.3) H η 3 Let u note that in the firt order term we can alway afely identify η with it unperturbed value η ).

9 The firt-order, non-homogeneou, non-averaged) expreion of d L in our perturbed background, referred to the unperturbed value.), i finally given by: d L z, θ a ) + z )a η d Lz, θ a ) d F L LRW z ) Ψη, η η, θ a ) + Ψ av + ) J J. 3.3) H η If we apply thi general reult to the particular cae of a CDM-dominated Univere we find almot full agreement with the reult for the luminoity ditance at contant redhift computed in. After everal manipulation, in fact, it turn out that our d L i equivalent to the one found in, modulo a term which can be written a v ˆn in the notation of Eq. 3.5)). Such a term give a ubleading contribution to the backreaction and can be neglected with no impact on our final reult. It hould be noted, however, that by expanding the above expreion in the mall z limit, and comparing the reult with the analogou expanion of the homogeneou ditance.3), we could eaily introduce a redefined value of H, ay H ren, uch that it invere correpond to the coefficient of the linear term in z for the perturbed relation. With uch a renormalized Hubble parameter we have that H ren d L tend moothly to H d F L LRW for z, and we recover full agreement with provided the ame renormalization i applied there. Such a renormalization of H i alo uggeted by and cloely related to) the firt-order computation of the calar expanion factor µ u µ for the flow wordline u µ µ τ of a local geodeic oberver. At thi point we could go on by taking the average of H ren d L z, θ a ), computing the aociated backreaction, and evaluating the correction to the tandard homogeneou relation given by H d F L LRW z ). We have performed that exercie, but we have found that the contribution of the renormalization term give large z -independent contribution to the variance. Furthermore, we think that renormalizing H at z i phyically incorrect ince, at very mall-z, the backreaction i dominated by hort-cale inhomogeneitie which are deeply inide the non-linear regime where even the concept of a Hubble flow become inappropriate). We could intead try to renormalize H at ome mall but finite z, e.g. at a redhift correponding to the cloet ued upernovae ay z.5, ee e.g. ), but then the reult although much better behaved) would depend on the choice of the particular renormalization point. Thu, it eem bet to conider jut the full unrenormalized expreion 3.3) in a limited region of z where one can trut the approximation made, and ue that expreion for a phenomenological parametrization of the backreaction effect which could poibly include a redefinition of H. 9 IV. COMBINING SPACE-TIME AND ENSEMBLE AVERAGES In the following ection the inhomogeneou deviation from the tandard FLRW quantitie are ourced by a tochatic background of primordial perturbation, atifying Ψ, Ψ, where the bar denote tatitical or enemble) average ee Sec. V). Hence, if we limit ourelve to a firt-order computation of d L, we immediately obtain. Non-trivial effect can only be obtained from quadratic and higher-order perturbative correction, or from the pectrum of the two-point correlation function d L z, θ a )d L z, θ a ) dicued in detail in ). In thi paper we will conider the enemble average not of d L but of d L, where the angular bracket refer to the light-cone average defined in Eq..9) ee e.g. 3, 4 for previou attempt of combining enemble average with average over pacelike hyperurface). We will ee that the light-cone average automatically induce quadratic and higher-order) backreaction term, due to the inhomogeneitie preent both in d L and in the covariant integration d L d F LRW L meaure, and ince the enemble average of uch term i non-vanihing we obtain, in general, d L d F L LRW. We will tart thi ection with ome general conideration on how to combine pace-time and enemble averaging, and how to iolate thoe term in d L that we may genuinely call backreaction effect, i.e. effect on averaged quantitie due to inhomogeneitie. We hall alo dicu how to etimate the variance around mean value due to uch fluctuation 4. Many of thee conideration can be certainly found elewhere, but are nonethele preented here for the ake of being elf-contained. Let u conider a typical average over the compact urface Σ topologically equivalent to a two-phere) embedded on the pat light-cone w w at contant z. We imply denote uch an average by: S Σ Σ d µ S Σ d µ, 4.) 4 The importance of the comic variance for a precie meaurement of the comological parameter, taking into account backreaction effect from averaging on domain embedded in a patial hyperurface according to 6), ha been recently pointed out alo in 5.

10 where d µ i the appropriate meaure provided by our gauge invariant precription ee Eq..9)) and S i the poibly non local) calar obervable d L in our cae). We can conveniently extract, from both d µ and from S, a zeroth-order homogeneou contribution by defining: d µ d µ) ) + µ), S S ) + σ), 4.) and ue the poibility of recaling both integral in Eq. 4.) by the ame contant, in order to normalize d µ) ). We then eaily get: S d µ) ) + µ) + σ) + µ) + σ), 4.3) S ) d µ) ) + µ) + µ where we have dropped, for implicity, the ubcript of the averaging region Σ, and where we have defined by... average with repect to the unperturbed meaure d µ) ) we hall drop the ubcript hereafter). Let u now perform the enemble average, denoted by an overbar, paying attention to the fact that enemble average do not factorize, i.e. AB A B. A imple calculation lead to: S/S ) + σ + µσ ) + µ ). 4.4) Thi lat equation i uppoedly exact but, a uch, pretty uele. It become an intereting equation, though, if we can expand the quantitie µ and σ in a perturbative erie: µ i µ i, σ i σ i, 4.5) and we further aume that the firt order quantitie µ, σ have vanihing enemble average a it i the cae for typical comological perturbation coming from inflation). In that cae we can eaily expand the reult and obtain, for intance: where S/S ) + σ + IBR + σ 3 + IBR ) IBR µ σ µ σ, 4.7) IBR 3 µ σ µ σ + µ σ µ σ µ µ σ + µ µ σ, 4.8) and where we have ued again the non-factorization property, i.e. µ σ µ σ, and o on. We ee that the reult contain, to a given order, both term that depend on expanding S to that order but not on the precie averaging precription), and induced backreaction IBR) term that depend on correlation between the fluctuation of S and thoe in the meaure. Thee latter term only depend on lower-order perturbation of S and the meaure eparately. In particular, our firt-order calculation provide the full econd-order IBR effect that come from the above interplay of µ and σ, although the full econd-order reult need alo the harder computation of σ but not of µ ). Note alo that, whenever the fluctuation of S and d µ are uncorrelated, all IBR effect drop out. We will ee below how to apply the above general reaoning to the particular cae of d L. Let u now dicu intead the iue of the variance, i.e. of how broad i the ditribution of value for S/S ) around it mean value S/S ). Thi diperion i due to both the fluctuation on the averaging urface and to thoe due to enemble fluctuation. Let u thu define: ) VarS/S ) S/S ) S/S ) S/S ) ) S/S ) ). 4.9) Inerting the definition 4.) we get, after a little algebra: ) VarS/S ) σ + µ) + µ ) σ + µ) + µ ). 4.) If we now make the ame aumption a before on expanding µ and σ, we find that the econd term in 4.) i at leat of fourth order. Therefore, for the leading term in the variance we find the amazingly imple reult ee alo 5): VarS/S ) σ. 4.)

11 A in the cae of IBR we only need to know the firt-order perturbation, but thi time the effect i completely independent of the averaging meaure. A we hall ee, the diperion which i actually the quare root of the variance) turn out to be larger than the averaging correction due to IBR. We may note, at thi point, that we could have alo conidered the diperion of the angular average S/S ) due to the tochatic fluctuation, namely: Var S/S ) which, after calculation imilar to the one above, give S/S ) S/S ) ) S/S ) ) S/S ) ), 4.) Var S/S ) σ. 4.3) Such a quantity i much maller than the previou one, a can be inferred from Section VI, indicating that the main reaon for the diperion lie in the angular catter of the data rather than in their tochatic ditribution due to the enemble. Let u finally identify the quantitie appearing in IBR for our particular cae. To thi purpoe let u tre that, to decribe the real impact of inhomogenetie on the obervational data, we hould conider the calar S correponding to the true oberved quantity. In the cae of the upernovae data thi quantity hould be the received flux of radiation, which i proportional to d L. On the other hand, if we conider the average of S d L intead of the average of the flux, the difference i only of econd order in our perturbative expanion 5. Since in thi paper we will not evaluate uch genuine econd-order contribution to the average we can limit ourelve, for the ake of implicity, to conider hereafter S d L z, θ a ). In uch a cae, from Eq. 3.3) we immediately find σ A + A + A 3 + A 4 + A 5, 4.4) where we have defined: A Ψ ; A Ψ av ; A 3 ) I + ; H η A 4 ) I r ; A 5 J. 4.5) H η On the other hand, in our particular cae, the general meaure d µ will be given by Eq..9): d µ d θ γw, τz, w, θ a ), θ a ). 4.6) We need to tranform the integral in d θ over the -phere Σ) to the tandard polar coordinate d θ dθdφ of the Newtonian gauge. To firt order, it i eay to check by taking into account the Jacobian determinant of the tranformation θ a θ a, ee Eq. 3.8)) that: d θ γ dφdθ in θ a r z, θ a )) Ψ ). 4.7) The normalized unperturbed meaure i then given by d µ) ) d Ω/4π, where d Ω dφdθ in θ. Proceeding a in the previou ection, and uing in particular Eq. 3.9), we alo obtain: µ Ψ + 4Ψ av + ) Jz, θ a ) A + A + A 3 + A 4 ). 4.8) H η The econd-order induced backreaction IBR, and the variance of d L /d FLRW L, can now be obtained by inerting the reult 4.4), 4.8) into Eq. 4.6) and 4.), repectively. Their explicit evaluation, firt analytic and then numerical, will be preented in the next two ection. 5 Since the flux i proportional to d L z, θa ) the ditance modulu ee Sec. VI) hould be computed a a function of d L ) / rather than of d L. Up to the econd order the difference between the two i hown in the following expreion: d L ) / d F L LRW + σ + IBR 3 ) ) σ ; d L d F L LRW + σ + IBR.

12 V. INDUCED BACKREACTION AND DISPERSION: ANALYTIC CONSIDERATIONS Thi ection will be devoted to a ytematic evaluation of the variou term contributing to the induced backreaction IBR defined in the previou ection, a well a to the variance aociated to the angle and enemble)-averaged value of d L z, θ a ). Let u tart by noting that the implet way to implement the enemble average of our tochatic background of calar perturbation Ψ i to conider their Fourier decompoition in the form: Ψη, x) π) 3/ d 3 k e i k x Ψ k η)e k), 5.) where auming that the fluctuation are tatitically homogeneou E i a unit random variable atifying E k) E k) a well a the enemble-average condition: E k )E k ) δ k + k ). 5.) According to the reult of Sec. IV, all correction we need to compute are bilinear term in the potential Ψ alway occurring in the combination A i A j or A i A j, where the quantitie A i, A j are defined in Eq. 4.5). Let u illutrate a typical computation with one of the implet contribution aociated to A Ψ. Following the notation of Sect. IV we will denote with the angular bracket the integration over the two-urface Σ embedded on the light-cone at w w, z z ). The meaure of integration i the unperturbed normalized one, i.e. d Ω/4π. We then obtain: d 3 k d 3 k Ψ Ψ π) 3 E k)e k d Ω ) Ψ k η )e ir k ˆx Ψ k η )e ir k ˆx 4π rη η rη η d 3 k π) 3 Ψ kη ) dco θ) e ik η co θ e ik η co θ d 3 k π) 3 Ψ kη ) dk k P Ψk, η ), 5.3) d 3 k d 3 k Ψ Ψ π) 3 E k)e k d Ω ) 4π Ψ kη )e ir k ˆx d Ω rη η 4π Ψ k η )e ir k ˆx rη η d 3 k π) 3 Ψ kη ) dco θ) e ik η co θ dco θ ) e ik η co θ d 3 ) k ink η) ) π) 3 Ψ kη ) dk ink η) k η k P Ψk, η ), 5.4) k η where in the econd line of both term we made ue of iotropy i.e. Ψ k only dependent on k k ), and defined θ and θ a the angle between k and x rˆx and between k and x rˆx. We have alo introduced the o-called dimenionle) power pectrum of Ψ, defined in general by: P Ψ k, η) k3 π Ψ kη). 5.5) More complicated example, that contain almot all the ubtletie of thee computation, will be preented in Appendix A. In general we have many contribution like the above one, appearing in both the induced backreaction and the variance, and generated by all the A i term of Eq. 4.5). In thi paper we will conider the particularly imple cae of a CDM-dominated background geometry, with a time-independent pectral ditribution of ub-horizon calar perturbation, η Ψ k. In uch a cae all the relevant contribution can be parameterized in the form A i A j A i A j dk k P Ψk)C ij k, η, η ), 5.6) dk k P Ψk)C i k, η, η ) C j k, η, η ) 5.7)

13 where C ij, C i are contant pectral coefficient), valid for any given model of time-independent) calar perturbation pectrum. With uch parametrization, the leading-order induced backreaction, called IBR in the expanion 4.6) of d L /d F L LRW, can be written a 4 5 dk 4 5 IBR A i A j A i A j k P Ψk) C ij k, η, η ) C i k, η, η ) C j k, η, η ). 5.8) i j The diperion from Eq.4.)) take intead the form: ) / d L 5 5 Var d F L LRW σ A i A j i j / i j dk k P Ψk) 5 i j 5 C ij k, η, η ) / ) We have analytically computed all the required coefficient C ij, C i, and their final value are given in Table I for C ij k, η, η ), and in Table II for C i k, η, η ). In Table I we alo how the mall-k limit k η ) of C ij and of the product C i C j. For notational convenience we have introduced in the table the dimenionle variable l k η, and we have ued the definition: SinIntl) l dx x in x. 5.) TABLE I: The pectral coefficient C ijk, η, η ) for the A ia j term defined by Eq. 5.6). We alo give the k η limit up to leading order in k η) of C ij and of the product C ic j for the coefficient defined in Table II. A ia j C ijk, η, η ) C ij for k η C i C j for k η A A l 3 A A SinIntl) + l + 4 l ) 9) 9 l ) A A 3 inl) l l l 6 6 ) ) ) A A 4 f inl) co l f l η l η f 3 η H η A A 5 inl) l l 3 8 A A + co l + lsinintl) 4 l 4 4 l 9 9 l ) A A 3 H ) ) η ) A A 4 f +f inl) f +f l f η l η 3 η A A 5 3l l ) co l + l in l + l 3 SinIntl) l ) 3 ) A 3A 3 inl) l l 3 ) ) A 3A 4 f f co l inl) f f l f η l η 3 η ) ) A 3A 5 inl) l l 3 ) A 4A 4 f +f l f f ) ) ) η 3 η co l + + l ) in l f f l f l η 3 η ) ) A 4A 5 f +3f co l + f f in l + f +f ) +lsinintl)) f f l η η l η η 3 A 5A 5 5l 4 + l + 4 l + l 4 ) col) +l 6 + l ) inl) + l 5 SinIntl) l 3 ) l 4 ) l 3 It hould be noted that the computation of ome of the above term require the explicit expreion of the cale factor aη). In thoe cae we have aumed a dut-dominated cale factor with aη) aη )η/η ), and we have defined f, η, η in dη aη) aη, ) η3, η 3 in 3η, l 3 3 l 36 ) l 4 8 ) l 4 3 η, 5.) recall that η in atifie, by definition, η in η, ). Uing uch a CDM-dominated model alo impoe the following zeroth-order) relation between z and η: η + z ) /. 5.) a H 9

14 4 TABLE II: The pectral coefficient C i for the A i A j term defined by Eq. 5.7). A i C ik, η, η ) A A A 3 A 4 in l l SinIntl) l ) ) in l l ) f η co l in l l A 5 ) Thi relation will be ued in ection VI when performing the numerical integration over k a a function of z. In order to proceed with ome qualitative and quantitative conideration it i important, at thi point, to pecify the propertie of the power pectrum P Ψ k). Limiting ourelve to ub-horizon perturbation, and conidering the tandard CDM model, we can imply obtain Ψ k by applying an appropriate, time-independent tranfer function to the primordial inflationary) pectral ditribution ee e.g. 6). The power pectrum of the Bardeen potential i then given by P Ψ k) ) 3 5 RT k), ) n k R A k, 5.3) where T k) i a contant tranfer function which take into account the ub-horizon evolution of the mode re-entering during the radiation-dominated era, and R i the primordial power pectrum of curvature perturbation, amplified by inflation, outide the horizon. The typical parameter of uch a pectrum, namely the amplitude A, the pectral index n and the cale k, are determined by the reult of recent WMAP obervation 7. In our computation we will ue, in particular, the following approximate value: A.45 9, n.96, k /a. Mpc. 5.4) Finally, ince we are mainly intereted in the overall magnitude of the tranfer function, it will be enough for our need to approximate T k) with the effective hape of the tranfer function for denity perturbation 6 with no baryon ocillation, namely T k) T k), where 8: T q) L L + q C q), L 73 q) lne +.8q), C q) q, q k, 5.5) 3.4k eq and where k eq i the cale correponding to matter-radiation equality. We can eaily ee that the above tranfer function goe to for k k eq, while it fall like k log k for k k eq. In the next ection we will numerically evaluate the coefficient of Table I and II, and we will dicu the correponding backreaction effect on d L, together with it diperion, for the calar perturbation pectrum decribed by Eq. 5.3) and 5.5). VI. INDUCED BACKREACTION AND DISPERSION: NUMERICAL RESULTS AND DISCUSSION In the following computation we will et a, Ω m, and we will ue h H / km Mpc ).7. In that cae we obtain 8 k eq.36 Mpc, and we can more preciely define the aymptotic regime of our tranfer function a T for k < 3 Mpc, and T k log k for k >.5 Mpc. Thi econd cale i already deep inide the o-called non-linear regime, roughly correponding to k >. h Mpc ee e.g. 9). Let u tart by comparing the behaviour of the power pectrum with the behaviour of the coefficient in Table I and II. A hown in Table I, all the coefficient C ij and C i C j will go at mot a O) while their combination C ij C i C j will vanih at leat a Ok η ) in the IR limit k η. A a conequence, the infrared part of the term in the induced backreaction and in the diperion will give a ubleading contribution, and we can afely fix our infrared cut-off to be k H i.e. we can limit ourelve to ub-horizon mode without affecting the final reult. 6 The relation between denity perturbation and metric perturbation i of coure fully under control in the linear perturbative regime, but could be unreliable in the high-k non-linear regime.

15 Furthermore, a can be een in Table I and II, all the coefficient C ij and C i C j will go at mot a Ok 3 η 3 ) in the UV limit k η. A a conequence, all the integral involved in the induced backreaction and in the diperion will be UV-convergent and their main contribution will come from the range / η k.5 Mpc. In thi range, ome of the backreaction coefficient grow a poitive power of k while the tranfer function i not yet decreaing a log k/k. In particular, there are only two leading contribution correponding to the integral controlled by the coefficient C 44 and C 55 thi i the reaon why the diperion of the angular average of d L, 4.3), i maller than the one of d L itelf 4.)). For uch coefficient we have, in fact, the following behaviour for k η : C 44 ) f + f ) k H η 3, 6.) C 55 k3 η 3 SinIntk η). 6.) 5 From the above expreion it i eay to undertand that A 4 A 4 will give the larget contribution for z ee Fig. ), while A 5 A 5 will give the larget contribution for z ee Fig. ). It i alo clear that both the induced backreaction IBR and the diperion depend in principle on the UV cut-off k UV eventually ued to evaluate the integral. On the other hand, when k UV i taken inide the regime where the pectrum goe like log k) /k 4, the dependence on the cut-off will not be too trong. In particular the leading contribution to IBR, given by A 4 A 4, depend very weakly on the particular value of k UV. The leading contribution to the diperion 5.9), controlled by A 5 A 5, ha a omewhat tronger dependence on k UV becaue of the extra power of k in the integrand of A 5 A 5 ). The numerical integration of A 4 A 4 and A 5 A 5 are preented in Fig. and, where we illutrate the magnitude of the backreaction effect a a function of the redhift in the range of value relevant to upernovae obervation), and a a function of it dependence on the cut-off ranging from k UV. Mpc to k UV + ). We hould emphaize that the two abovementioned contribution have a clear and ditinct phyical meaning in the Newtonian gauge. Going back to their explicit expreion it appear that A 4 A 4 repreent a Doppler effect while A 5 A 5 i aociated with lening, both in qualitative agreement with previou claim in the literature 9,. Note that the latter term can only appear in the variance, becaue it contribute to σ, while A 4 A 4 appear alo in IBR ince A 4 i preent in both µ and σ A4 A FIG. : The reult of the numerical integration of Eq. 5.6) for A 4A 4 i plotted a a function of z for three different value of the UV cut-off: k. Mpc thick blue line), k Mpc thin purple line), k + dahed red line).

16 6..8 A5 A FIG. : The reult of the numerical integration of Eq. 5.6) for A 5A 5 i plotted a a function of z for three different value of the UV cut-off: k. Mpc thick blue line), k Mpc thin purple line), k + dahed red line). It i important to tre that, although all conidered backreaction contribution are at any redhift) UV-finite, they can induce relatively big effect at the cale relevant for upernovae obervation ee Fig. and ), provided the given pectrum i extrapolated to ufficiently large value of k. Let u now um up all contribution to the induced backreaction 5.8) and to the diperion 5.9), and compare the reult for d L ± d CDM L σ with the homogeneou luminoity-ditance of a pure CDM model and of a ucceful ΛCDM model. We will include into d L only the IBR contribution 4.7), namely we will et d L d CDM L +IBR ). It i clear that a full computation hould include additional contribution, of the ame order a IBR or even larger, ariing from econd-order perturbation of d L we are referring to the term called σ in Eq. 4.6)). Nonethele, we believe that modulo cancellation) our computation may etimate a reliable lower limit on the trength of the poible correction to the luminoity-redhift relation in the context of our inhomogeneou geometry. The comparion between the homogeneou and inhomogeneou averaged) value of d L can be conveniently illutrated by plotting the o-called ditance modulu m M 5 log d L or, even better, by plotting the difference between the ditance modulu of the conidered model and that of a flat, linearly expanding Milne-type geometry, ued a reference value ee e.g. 8). In uch a cae we can plot, in particular, the following quantity: + z )z m M) 5 log d L 5 log. 6.3) H The reult are illutrated in Fig. 3 for the cae of a cut-off k UV. Mpc, and in Fig. 4 for k UV Mpc. The averaged luminoity-redhift relation of our inhomogeneou model i compared, in particular, with that of a pure CDM model and with that of a ΛCDM model with Ω Λ. and with Ω Λ.73. We have alo explicitly hown the expected diperion around the averaged reult, by plotting the curve correponding to d L ±d CDM L the coloured area appearing in the figure). σ bounding Let u tre again that the choice of the cut-off may affect even if not dramatically) the final reult when the value of k UV are varying in the range. ) Mpc, while the precie choice become le important at higher value of k UV. We hould recall, however, that the inhomogeneou model adopted in thi paper i fully under control only in the linear perturbative regime, and that the pectrum cannot be extrapolated at cale higher than about k Mpc without taking into account the complicated effect of it non-linear dynamical evolution.

17 7.5 m M FIG. 3: The ditance-modulu difference of Eq. 6.3) i plotted for a pure CDM model thin line), for a CDM model including the contribution of IBR dahed blue line) plu/minu the diperion coloured region), and for a ΛCDM model with Ω Λ.73 thick line) and Ω Λ. dahed-dot thick line). We have ued for all backreaction integral the cut-off k. Mpc...5 m M FIG. 4: The ditance-modulu difference of Eq. 6.3) i plotted for a pure CDM model thin line), for a CDM model including the contribution of IBR dahed blue line) plu/minu the diperion coloured region), and for a ΛCDM model with Ω Λ.73 thick line) and Ω Λ. dahed-dot thick line). We have ued for all backreaction integral the cut-off k Mpc.