Light Propagation in the Averaged Universe. arxiv:

Light Propagation in the Averaged Universe arxiv: 1404.2185 Samae Bagheri Dominik Schwarz Bielefeld University Cosmology Conference, Centre de Ciencias de Benasque Pedro Pascual, 11.Aug, 2014

Outline 1 Motivation Introduction The Averaging Problem 2 Averaging procedures Two Approaches 3 Light Propagation Geodesic Equation Two Sets of Equations 4 Summary and Outlook

Introduction We observe the homogeneous matter distribution of the Universe, on scales larger than 100Mpc.... [D. Hogg et al., 2001 (SDSS analysis)] [C. Blake et al., 2007 (WiggleZ Dark Energy Survey)]

On small scales we see the lumpy universe with inhomogeneous distribution of matter.

Effect of Inhomogeneities Matter distribution of real lumpy universe (on small scale) T µν (local) = T µν (discrete) R µν (local) 1 2 g µν (local) R (local) = kt µν (local) Matter distribution of homogeneous universe (on large scale) T µν = T (discrete) µν T (discrete) µν T (discrete) µν R µν 1 2 g µν g αβ R αβ + Q µν = k T (discrete) µν where Q µν is the backreaction term.

Structure in simulation on the top pannel, and the averaging problem in the bottom row. [C. Clarkson et. al. 2010]

The Averaging Problem The gravitational field equations on cosmological scales are obtained by averaging the Einstein field equations of general relativity. Averaging can lead to some effects in GR, because of non-comutativity with time evolution [ t, ] 0 non-linearity of Einstein s equations G µν ( g µν ) G µν (g µν )

procedures The averaging involves integration of tensor operation. Changing the coordinates will change the results in an arbitrary way. Two procejures to handle this problem: 1- Defining a covriant averaging of tensors, via bilocal operators. Proposed by Zalaletdinov (1993) 2- Using only field equations involving averaged scalars. Proposed by Buchert (2000)

Approach 1 Zalaletdinov s macroscopic gravity: Arbitrary tensors can be averaged here. P a (x) = 1 d 4 x g V A a a (x, x )p a (x ) Σ Where A a a (x, x ) = e a I (x)ei a (x ) is a bilocal operator. Σ The averaged Einstein equations reduce to Macroscopic Field Equations: g βɛ M ɛβ 1 2 δɛ γ g µν M µν = 8πG T ɛ γ + C ɛ γ Tγ ɛ is the averaged energy-momentum tensor, Cγ ɛ is the correction term, given by Cγ ɛ = g µν (Zµνγ ɛ 1 2 δɛ γzµνα) α Where Z α µνβ = 2Z α ɛ µɛ νβ Details

Approach 2 Buchert s spatial averaging of scalars: The metric can be written in the synchronous gauge ds 2 = dt 2 (3) g µν dx µ dx ν where (3) g µν is the metric on hypersurface of constant t. The spatial average of a scalar quantity f f (t, x) = 1 d 3 xf (t, x) g V D The obtained Bucherts equations are D 3äD a D = 4πG N ρ + Q 3ȧ2 D a 2 D = 8πG N ρ 1 2 R 1 2 Q Details

Brief Overview: Problem of covarinat averaging Can t average tensor fields in covariant way (coordinate dependent results) Can use bitensors for curvature and matter, but not for metric itself, and leads to complex set of equations (Zalaletdinov). A less ambitious road: averaging scalars only, It seems to retain the main features, gives modified Friedmann equations (Buchert).

Introduction Effects of inhomogeneities on the cosmological expansion: 1. Effect on the dynamics through backreaction or correlation terms. Applying Zalaletdinov or Buchert s approach on Einstein equations and Friedmann equations A. Coley et al. (2002), S. Räsänen (2004), E. Kolb (2005), A. Paranjape, T. Singh (2006), N.Li, D. Schwarz (2007) 2. Effect on light propagation. Investigating geodesic equation of motion in an averaged universe S. Räsänen (2008), C. Clarkson (2008), A. Coley (2009)

Geodesic Equation Considering geodesic equation for a vector field k µ k µ,νk ν + Γ µ νρk ν k ρ = 0 (k µ,νg µα + k µ g µα,ν + k µ,αg µν )k ν = 0 By averaging the space time, we notice g µα,ν g µα,ν In the averaged form we have (k,ν g µ µα + k,α g µ µν + k µ g µα,ν )k ν = ( g µα,ν g µα,ν )k µ k ν The modified term: T µαν = g µα,ν g µα,ν T µαν k µ k ν = I α

Geodesic Equation Property (1): T sym µαν = (T µαν + T ναµ )/2 is a symmetric tensor. Property (2): I α k α = 0. Assume Flat Friedmann-Lemaître Model: - averaged geometry g µν = ḡ µν - four-velocity ū µ = ( 1, 0) - and frequency ω = ū µk µ (k 0 = ω and k i = ωe i /a)

The general algebraic structure for T sym µαν: T sym µαν = f 1 2 (ḡ µαū ν + ḡ να ū µ ) + f 2 ū µ ū α ū ν + f 3 ḡ µν ū α f 1, f 2 and f 3 are functions of cosmic time. The contraction I α = k µ T sym µαν k ν = f 1 ( ω)ḡ µα k µ + f 2 ω 2 ū α I α k α = f 2 ω 3 must vanish as shown before. Thus f 2 0. Thus the inhomogeneity of the light propagation equation: I α = f 1 ωḡ µα k µ

Geodesic Equation Geodesic equation (background level) (k µ,νg µα + k µ g µα,ν + k µ,αg µν )k ν = 0 Geodesic equation (averaged form) (k µ,ν g µα + k µ,α g µν + k µ g µα,ν )k ν = T µν k µ k ν I α By considering isotropic homogeneous spatially flat FL universe as the background, we have two sets of equations for α = 0 and α = i.

First Set of Equations (α = 0) ( 2a 2 k 0,0 2a2 Hk 0 )k 0 + ( a 2 k 0,j + a2 δ ij k i,0 )kj = 0 ( 2a 2 k 0,0 2a2 Hk 0 )k 0 + ( a 2 k 0,j + a2 δ ij k i,0 )kj = T 000k 0 k 0 + T i0j k i k j where H = ȧ a. Averaged geodesic equation (time component) ( ω)( ω + ω,i e i /a + Hω) = I 0 = f 1 ω 2 The equation in terms of the affine parameter dω ω dλ = kµ x : µ dω dλ + H eff ω 2 = 0 H eff H + f 1 is the effective Hubble rate.

First Set of Equations (α = 0) ( 2a 2 k 0,0 2a2 Hk 0 )k 0 + ( a 2 k 0,j + a2 δ ij k i,0 )kj = 0 ( 2a 2 k 0,0 2a2 Hk 0 )k 0 + ( a 2 k 0,j + a2 δ ij k i,0 )kj = T 000k 0 k 0 + T i0j k i k j where H = ȧ a. Averaged geodesic equation (time component) ( ω)( ω + ω,i e i /a + Hω) = I 0 = f 1 ω 2 The equation in terms of the affine parameter dω ω dλ = kµ x : µ dω dλ + H eff ω 2 = 0 H eff H + f 1 is the effective Hubble rate.

Second Set of Equations (α = i) k j,0 k0 ḡ ij + k j,l kl ḡ ij + k 0,i k0 ḡ 00 + k j,i kl ḡ jl + k j k 0 ḡ ij,0 + k j k l ḡ ij,l = 0and 0 a 2 (k j,0 k0 δ ij + k j,l kl δ ij + k j,i kl δ jl + 2Hk j k 0 δ ij k 0,i k0 ) = 2T ij0 k 0 k j Modified light propagation equation (spatial components) aγ ij e j ω( ω + ei a ω,i + Hω) + aω 2 γ ij (ė j + e j k ek ) = I i = ω 2 aγ ij e j f 1 where ḡ ij = a 2 γ ij, and denotes a covariant derivative with respect to the 3-metric γ ij. ė j + e j k ek = 0 De i Dλ = 0 i.e. light rays propagate alonge straight lines.

Estimation for f 1 (t) ds 2 = e 2φ dt 2 + a 2 (t)e 2ψ γ ij dx i dx j g µν = ḡ µν + δg µν, by construction δg µν = 0 φ 0 (ISW effect) e 2φ φ 0. T sym 000 = δg 00,0 = 2 e 2φ φ f 1 = 2 e 2φ φ H eff = H + 2 e φ φ Following δ(r, t) ρ(r,t) ρ(t) ρ(t) For an over-dense region (δ > 0) = φ < 0 H eff < H. For an underdense region (δ < 0) = φ > 0 H eff > H.

Summary A cosmological solution obtained by averaging underlying inhomogeneous Universe. Buchert s and Zalaletdinov s averaging approaches, leads to space time corrections to the standard cosmological equations. The aim of our work was to investigate the effects of averaging on light propagation. Our central result is a modification of null geodesic equation, which is a fully covariant vector equation. By applying this light propagation equation to a flat, spatially isotropic and homogeneous averaged metric (compatible with the standard model of cosmology), This modification can be expressed as an effective Hubble rate So far we studied a single light ray. The next logical step is to study the equation of geodesic deviation to find a modification to the luminosity and angular diameter distances.

Macroscopic Gravity (MG) Zalaletdinov s macroscopic gravity: Arbitrary tensors can be averaged here. P a (x) = 1 d 4 x g A a a V (x, x )p a (x ) Σ Σ Where A a a (x, x ) = e a I (x)ei a (x ) is a bilocal operator. - The averaged of the connection Γ µ ρν on M is by definition taken to be the connection Γµ ρν of the averaged manifold M. - The averaged metric g µν can be chosen as the metric on averaged manifold.

The macroscopic Riemann tensor is M µ ναβ M µ ναβ = α Γµ νβ β Γ µ να + Γµ σα Γσ νβ Γµ σβ Γσ να M µ ναβ Rµ ναβ where Rµ ναβ is the average of the microscopic Riemann tensor. connection correlation tensor defined by Z α µ β[γ νσ] = Γα β[γ Γµ νσ] Γα β[γ Γµ νσ] The averaged Einstein equations reduce to Macroscopic Field Equations: g βɛ M ɛβ 1 2 δɛ γ g µν M µν = 8πG T ɛ γ + C ɛ γ T ɛ γ is the averaged energy-momentum tensor, C ɛ γ is the correction term, given by C ɛ γ = g µν (Z ɛ µνγ 1 2 δɛ γ Z α µνα ) Where Z α µνβ = 2Z α ɛ µɛ νβ Back

Buchert s spatial averaging approach Buchert s spatial averaging of scalars: The metric can be written in the synchronous gauge ds 2 = dt 2 (3) g µνdx µ dx ν where (3) g µν is the metric on hypersurface of constant t. The spatial average of a scalar quantity f f (t, x) = 1 d 3 xf (t, x) g V D D where V D = D d 3 x g. The starting point is expansion and averaging do not commut: where Θ is a expansion rate. t f D tf D = Θf D Θ D f D Θ µ ν = 1 2 g µρ tg ρν; Θ = Θ µ µ

1 The averaged scale factor is defined as a D V 3 D, which ensures Θ D = 3 tln(a D ) = 3H D. The obtained Bucherts equations are 3 äd = 4πG N ρ + Q a D where Q is the backreaction term 3 ȧ2 D ad 2 = 8πG N ρ 1 2 R 1 2 Q Q = 2 3 [ Θ2 D Θ 2 D ] 2 σ2 D and σ is a shear scalar σ µ ν = Θµ ν 1 3 δµ ν Θ; σ2 = 1 2 σµ ν σν µ. Back