Noname manuscript No. (will be inserted by the editor)

Averaging in LRS class II spacetimes Petr Kaˇ spar† · Otakar Sv´ıtek§

arXiv:1412.7208v1 [gr-qc] 22 Dec 2014

Received: date / Accepted: date

Abstract We generalize Buchert’s averaged equations [Gen. Rel. Grav. 32, 105 (2000); Gen. Rel. Grav. 33, 1381 (2001)] to LRS class II dust model in the sense that all Einstein equations are averaged, not only the trace part. We derive the relevant averaged equations and we investigate backreaction on expansion and shear scalars in an approximate LTB model. Finally we propose a way to close the system of averaged equations. Keywords LRS family; Cosmology; Averaging

1 Introduction Our universe is considered to be homogeneous and isotropic on the large scale leading to the FLRW model. However, if we move to smaller scales, we can observe a strongly inhomogeneous distribution of structures. If we want to deal with inhomogeneity rigorously and at the same time keep a consistent connection with the FLRW geometry we may consider an averaging formalism to smooth out the metric tensor and at the same time average Einstein equations as well. The problem is that Einstein equations are nonlinear and if we average them straightforwardly we do not obtain an averaged metric tensor as a solution of averaged equations. Instead, we should consider an additional term — the so-called correlation term, which can change the evolution of a smooth metric tensor and lead to the so-called backreaction. This term arises due to the nonlinearity of Einstein equations. It does not need to satisfy the usual energy conditions so it can possibly act as dark energy. While building a rigorous averaging scheme we face the problem that the average value of a tensor field is not well defined. There are several different Institute of Theoretical Physics Faculty of Mathematics and Physics, Charles University in Prague V Holeˇsoviˇ ck´ ach 2, 180 00 Praha 8, Czech Republic E-mail: † [email protected], § [email protected]

2

Petr Kaˇspar† , Otakar Sv´ıtek§

approaches to define averages of tensors. One of the most promising ones is the scheme by Zalaletdinov [1], [2] where not only Einstein equations but also Cartan structure equations (and their integrability conditions) are averaged. A theorem about isometric embedding of a 2-sphere into Euclidian space is applied in the averaging method developed by Korzy´ nski [3]. In [4] Weitzenb¨ock connection for parallel transport is used to define the average value of a tensor field. One of the most popular approaches to averaging is the one investigated by Buchert [5], [6], where only the scalar part of Einstein equations is averaged. Wiltshire used this approach to give an alternative explanation of cosmic acceleration [7]. This theory was also applied to the cosmological perturbation theory [8], [9], [10], [11]. For observational issues see e.g. [12]. In this paper, we will generalize Buchert’s equations to the locally rotationally symmetric (LRS) class II dust family of spacetimes. The LRS family was classified in [13], [14] and recently in [15]. It contains e.g. LRS Bianchi cosmologies, KantowskiSachs model or LTB model and its generalizations. We will use the fact that this family is described by scalars to average the complete set of Einstein equations, including constraints. Although the averaged constraints are shown to be preserved during evolution the averaged system of equations is not closed and additional information has to be supplemented. In the past there have been many attempts to apply Buchert’s approach to the LTB model. Papers [16] and [17] are comprehensive studies of Buchert’s formalism applied to generic LTB models. For a treatment of Buchert’s equations inside LTB spacetime see e.g. [18], [19], [20], [21], [22], [23], [24] and [25]. For an application of Buchert’s formalism to the structure formation see e.g. [26], [27] and [28]. Naturally, one can study inhomogeneities perturbatively on a homogeneous background and many important results are based on this approach. However, we should be cautious about relying solely on a linear perturbative analysis when dealing with a nonlinear theory. The effects of the correlation term indicate what kind of effects one might be missing when using a simple approach. In this sense, rigorous averaging of exact inhomogeneous spacetimes leading to standard cosmological models provides a possibility to qualitatively estimate these effects.

The paper is organized as follows. In Section 2 we review the LRS family and its characterizations, then we briefly mention Buchert’s equations. Next, we average equations describing dust LRS class II family. After a short review of LTB metric in Section 5 we investigate the backreaction in the so-called onion model. We proceed by attempting to close the averaged equations and we finish with conclusion.

Averaging in LRS class II spacetimes

3

2 LRS family Locally rotationally symmetric (LRS) dust spacetimes are defined by the following features [13]: In an open neighborhood of each point p, there exists a nondiscrete subgroup of the Lorentz group which leaves invariant the Riemann tensor and its covariant derivatives up to the third order. Therefore, in LRS spacetimes there exists a preferred direction eµ (the axis of symmetry) at every point. The subgroup can be one or three-dimensional. In the latter case, we can rotate the axis of symmetry and spacetimes are everywhere isotropic these are the FLRW models. We will use the covariant 1+3 splitting of spacetime with the timelike vector field uµ normalized by the condition uρ uρ = −1 and the projection tensor hµν = gµν + uµ uν . In this section we will follow the article of van Elst and Ellis [15]. The preferred spacelike vector field eµ satisfies the following conditions: eρ uρ = 0,

eρ eρ = 1.

(1)

Because of the defining property of the LRS spacetime, all covariantly defined spacelike vectors orthogonal to uµ (acceleration u˙ µ , vorticity ω µ , projected gradient of density hσµ ∇σ ρ, pressure hσµ ∇σ p and expansion hσµ ∇σ θ) must be proportional to eµ - if this condition does not hold, spacelike vectors will not be invariant under the rotation about eµ . u˙ µ = ue ˙ µ, hσµ ∇σ ρ = ρ0 eµ ,

ω µ = ωeµ ,

hσµ ∇σ p = p0 eµ ,

(2)

hσµ ∇σ θ = θ0 eµ .

(3)

µ

Dots denote covariant derivative along the flow vector u and primes denote covariant derivative along the vector eµ . We define the magnitude of the spatial rotation k and the magnitude of the spatial divergence a by (4) k := η αβγδ (∇β eγ ) uδ ,
a := hαβ ∇α eβ ,

(5)

where η αβγδ is the totally antisymmetric Levi-Civita pseudotensor (η0123 = √ − −g). A similar rule works also for the spacelike tracefree symmetric tensors orthogonal to uµ . Following [15] we introduce a tensor field eµν defined by eµ eµν :=

1 (3eµ eν − hµν ) . 2

(6)

Then we have the relations for the shear tensor and the electric and magnetic parts of the Weyl tensor 2 σµν = √ σeµν , 3

2 Eµν = √ Eeµν , 3

2 Hµν = √ Heµν . 3

(7)

Petr Kaˇspar† , Otakar Sv´ıtek§

4

Here we can see that the LRS spacetimes are characterized only by a finite set of scalar functions. For simplicity we will restrict our attention to the LRS class II dust models defined by the relation k = ω = 0. One can also show that the magnetic part of the Weyl tensor is equal to zero, H = 0. This family of spacetimes includes the LTB metric and its generalizations based on foliation by spacelike 2-surfaces with negative or zero curvature scalar. The relevant evolution equations are 1 θ˙ = − θ2 − 2σ 2 − 4πρ, 3 1 2 σ˙ = − √ σ 2 − θσ − E, 3 3 √ E˙ = −4πρσ + 3Eσ − θE, ρ˙ = −ρθ, 1 1 a˙ = − aθ + √ aσ, 3 3

(8) (9) (10) (11) (12)

and the constraints 1 2 σ 0 = √ θ0 − aσ, 3 3 4π 3 E 0 = − aE + √ ρ0 , 2 3 2 4 2 1 16π 2 ρ. a0 = θ2 + √ θσ − σ 2 − √ E − a2 − 9 3 2 3 3 3 3

(13) (14) (15)

If we take the time derivative of the constraints, we can prove that they do not change with time.

3 Buchert’s equations In this section we will review an averaging method developed by Buchert [5]. We will consider only the dust case - for generalization to perfect fluid see [6]. This approach uses 1+3 splitting of spacetime, which is well defined by irrotational dust 4-velocity. However, averaging is well defined only for scalars, therefore only scalar part of Einstein equations is averaged. Given a scalar field A, the average value over three-dimensional spacelike domain D is defined by 1 hAiD = VD

Z D

1 d XJA = VD 3

Z

d3 X

p det gij A,

(16)

D

p where J := detgij , gij is the metric of the spacelike hypersurface, X i are the comoving coordinates and VD is the proper volume of the three-dimensional

Averaging in LRS class II spacetimes

5

domain D. From this definition we can see that time derivative and averaging do not commute. We have a commutation relation · hAiD

  Z Z   d  1 V˙ D 1 3 ˙ + J A˙  = d XJA = − hAiD + d3 X JA dt VD VD VD D D D E = − hθiD hAiD + hAθiD + A˙ , (17) D

where the expansion rate θ is related to the velocity of the fluid uµ according to the definition θ = uµ;µ . Next, in analogy with the FLRW spacetime we introduce the dimensionless scale factor aD and the effective Hubble parameter HD 1  VD 3 , (18) aD = VDi

hθiD =

V˙D a˙D =3 =: 3HD . VD aD

(19)

VDi is the volume of the initial domain which is geodetically evolved into VD . Now we have a formalism for averaging scalars. To obtain a scalar equation from Einstein equations, we have to contract them with available tensors - i.e. g µν , uµ and ∇µ . After contracting we get the Raychaudhuri equation, the Hamiltonian constraint and the continuity equation. Now we perform averaging and use the commutation rule (17) 3



a˙ D aD

a¨D + 4πG hρiD = QD , aD

2 −

(20)

hRiD QD 8πG hρiD + =− , 3 6 6

(21)

a˙D hρiD = 0. aD

(22)

∂t hρiD + 3

hRiD denotes the average value of the spatial Ricci scalar, hρiD means the average density of the fluid and QD that shows possible backreaction (due to inhomogeneity and anisotropy) is defined by QD :=

E

 2D 2 (θ − hθiD ) − 2 σ2 D . 3 D

The scalar σ 2 = 21 σµν σ µν is constructed from the shear tensor.

(23)

Petr Kaˇspar† , Otakar Sv´ıtek§

6

4 Averaging LRS class II dust spacetime Now, we will generalize the above approach to LRS class II dust solutions. Originally, Buchert considered spacetimes with a dust [5] or a perfect fluid [6] source. He did not assume any symmetries or simplifications and his equations can be applied to a large class of metrics. Here we will restrict to spacetimes with the special LRS symmetry. For this family we will generalize Buchert’s equations in the sense that all Einstein equations are averaged consistently. Given a preferred spacelike direction eµ , all the equations describing the LRS metric are scalar. It means that we can perform averaging (which is covariantly defined for scalars). We will define averaging over the spacelike domain D according to (16). In order to obtain averaged equations we need to derive commutation relations for the time and spatial derivatives (with respect to the preferred direction). For the LRS class II spacetime it is possible to express √ derivative along the preferred direction e by the formula e = g11 ∂1 , where g11 is the metric function of a particular solution inside the LRS class II spacetime. The most studied solution inside this class is the LTB model, for which √ R0 . We will show the the square root of the metric function reads g11 = √1+2B basic facts about the LTB metric in the next section. The commutation rule between coordinate derivative and averaging reads     ∂1 J ∂1 J ∂1 hAi − h∂1 Ai = A − hAi . (24) J J For simplicity we will restrict to the class II LRS spacetime with the condition p = 0 ⇔ ρ˙ = −ρθ (dust models) which includes LTB spacetimes and their generalizations. For simpler notation we shall omit symbol D at the averaging bracket in the rest of the paper (however we retain the symbols not defined as direct averages - aD , QD ). If we average the equations (8) - (15) we obtain 

 2  2  1 2 2 · θ − hθi − 2 σ 2 , (25) hθi = − hθi − 4π hρi + 3 3 1 1  2  2 · 2 hσi = − √ hσi − hθi hσi − hEi + √ hσi − σ 2 3 3 3 1 + (hθσi − hθi hσi) , (26) 3 √ · hEi = −4π hρi hσi + 3 hEi hσi − hθi hEi √ −4π (hρσi − hρi hσi) + 3 (hEσi − hEi hσi) , (27) ·

hρi = − hρi hθi , 1 1 2 · hai = − hai hθi + √ hai hσi + (haθi − hai hθi) 3 3 3 1 + √ (haσi − hai hσi) , 3     1 2 ∂1 J ∂1 J √ 0 0 hσi = √ hθi − hai hσi + g11 σ − hσi 3 J J 3

(28)

(29) (30)

Averaging in LRS class II spacetimes

r

  1 − + hθi g11 √ aσ − hai hσi , g11     3 4π ∂1 J ∂1 J √ 0 0 hEi = − hai hEi + √ hρi + g11 E − hEi (31) 2 J J 3 r         ∂1 J ∂1 J 3 √ 1 g11 ρ − hρi − − 4π g11 √ aE − hai hEi , 3 J J 2 g11 0

hai = + + − −

g11 3



∂1 J θ J





7

∂1 J J



2 − 3



√



2 2 4 2 1 16π 2 2 2 hθi + √ hθi hσi − hσi − √ hEi − hai − hρi (32) 9 3 2 3 3 3 3         ∂1 J ∂1 J 2 √ 1 √ 2 g11 a − hai + g11 √ θ2 − hθi J J 9 g11         2 1 1 4 √ √ 2 2 √ g11 √ θσ − hθi hσi − g11 √ σ − hσi g11 3 g11 3 3         2 √ 1 1 √ 1 2 2 √ g11 √ E − hEi − g11 √ a − hai g11 2 g11 3     16π √ 1 g11 √ ρ − hρi . 3 g11

The underlined parts of the equations denote additional terms due to averaging. We can recognize the well-known Buchert’s equation (25) with the kinematical backreaction term and the mass conservation equation (28). If we want to restrict the above equations to the LTB model (the line element is given in the next section) we can substitute the following expressions for the magnitude of the spatial divergence a, itsp average value haiD and the square root of the spatial part of the metric J = detgij √ 4πR2 2 1 + 2B R0 R2 sin θ , haiD = . (33) a = hνµ eµ ;ν = , J= √ R VD 1 + 2B One can show that the averaged constraint equations (30) - (32) are preserved in time. The key role in the calculation is played by the equation [15]  0 2 1 · (34) (f 0 ) = f˙ − √ σf 0 − θf 0 , 3 3 and its averaged version. Now, we can take time derivative of the constraints (30) - (32). Using commutation rules and equations (8) - (15), a slow but straightforward computation will show that the constraints do not evolve in time. The explicit computation for equation (30) is shown in the Appendix. All of Einstein equations are averaged now. It means that we can investigate not only backreaction on the expansion rate but also on shear scalar or electric part of the Weyl scalar. The problem is that the equations are not closed. We need additional

 relations to close the system because for example hθi is independent of θ2 . In the next chapters we will give some suggestions for closing the system of equations.

Petr Kaˇspar† , Otakar Sv´ıtek§

8

So far we have not seen an analog of the Hamiltonian constraint. For example Sussman in [16] used three-dimensional curvature R instead of the function a. The relation between three-dimensional curvature R and the magnitude of the spatial divergence a reads [15]   2 R = − 2a0 + a2 − 2K , (35) 3 where K is the Gaussian curvature of the 2-D spacelike group orbits orthogonal to eµ and uµ . The Hamiltonian constraint has the form 2 R = 16πρ − θ2 + 2σ 2 . 3

(36)

The averaged Hamiltonian constraint is Buchert’s equation (21) 2 2 2 hRi = 16π hρi − hθi + 2 hσi 3     2  2  2 2 − θ − hθi + 2 σ 2 − hσi . 3

(37)

5 LTB metric The most important representative of the dust LRS class II family is the LTB spacetime. In this section we will briefly review its properties. The Lemaˆıtre-Tolman-Bondi (LTB) metric [29], [30], [31] is a spherically symmetric exact solution of Einstein equations. It corresponds to an inhomogeneous dust with the stress energy tensor Tµν = ρuµ uν ,

(38)

where uµ is 4-velocity of the dust with density ρ. For a recent review of LTB metric, see e.g. [32], [33]. The line element reads ds2 = −dt2 +

(R0 )2 dr2 + R2 (t, r)[dθ2 + sin2 (θ)dφ2 ], 1 + 2B(r)

(39)

where B(r) is an arbitrary function and the prime denotes partial derivative with respect to r. Function R(t, r) is a solution of Einstein equation 2 R,t = 2B +

2M Λ + R2 . R 3

(40)

M = M (r) is another arbitrary function of integration. The energy density ρ is determined by the equation 4πρ =

M0 . R0 R2

(41)

Averaging in LRS class II spacetimes

9

0.00025

·

0.00020

0.00015

0.00010

r 2

4

6

8

10

Fig. 1 Density profile at the time t = 20.

The function B(r) is related to the quasi-local spatial curvature [16] and M (r) is the gravitational mass contained within a comoving spherical shell at a given r. Equation (40) can be integrated to yield ZR

˜ dR q

0

2B +

2M ˜ R

˜2 + 13 ΛR

= t − tB (r).

(42)

Here tB (r) is a third free function of r (called the bang time function). In the LTB model, in general, the Big Bang is not simultaneous as in the FRW case, but it depends on the radial coordinate r. The given formulas are invariant under a transformation r˜ = g(r). We can use this freedom to choose one of the functions B(r), M (r) and tB (r). For Λ = 0 the above equation can be solved explicitly. The evolution can be elliptic (B < 0), parabolic (B = 0) or hyperbolic (B > 0).

6 Backreaction inside the LTB onion model As an example of backreaction computation we consider an approximate LTB model (the so-called onion model) investigated in [34] by Biswas, Mansouri and Notari, who computed corrections to the luminosity distance–redshift relation. It represents a spacetime with radial shells of overdense and underdense regions. The function B(r) is nonzero (B(r) > 0), so the evolution of the LTB model is hyperbolic. The metric function R(t, r) reads !  1/3  1/3   6 81 1 1 R(t, r) := t2/3 r 1 + t2/3 sin2 πr . (43) π 4000π 2 2π r The function B(r) is given as follows B(r) =

r sin2 πr. 2π

(44)

Petr Kaˇspar† , Otakar Sv´ıtek§

10 2

4

-0.0001

6

l

8

-0.0002 -0.0003 -0.0004 -0.0005

Fig. 2 Backreaction term

2 3



 

 θ2 − hθi2 − 2 σ 2 in the evolution equation for expansion

depending on the averaging scale l and normalized by hθi· .

The density profile at the time t = 20 can be seen in Figure 1. The density is computed using the formula (41). The coordinates were chosen so that the function M (r) is given by M (r) = 4/3πr3 . First, we investigate the backreaction term in Buchert’s equation (25). We numerically integrate the underlined part of equation (25) depending on the averaging scale l. As one can see from Figure 2, the backreaction normalized · by hθi is negative. It has a peak for the averaging scale l ≈ 0.8 . The value of · backreaction normalized by hθi is of the order of 10−4 . The backreaction term (without normalization) is positive and it leads to an increase of expansion. We can investigate also the backreaction terms in other equations which do not appear in the Buchert framework and which can supplement his equations. For example here we will show the result for backreaction in the averaged evolution equation for shear (26) (specifically the whole underlined part of the equation is considered). All results depend on the averaging scale l. As we can see from Figure 3 - for small scales, the contribution of all backreaction terms · in the evolution equation for shear normalized by hσi is negative with a peak around l ≈ 0.9 . For larger scales the contribution is smaller and positive. The turning point is for l ≈ 1.2 . To be more precise, there exist regions where the backreaction changes the sign twice for a very small increase of l. If we compare the backreaction with the time derivative of the shear scalar, we can see that their ratio is of the order of 10−4 − 10−3 . It means that the backreaction plays a more important role in the averaged equation for shear than in the averaged equation for expansion. Note that we investigated only an approximate LTB model. Due to nonlinearity it is not clear if the backreaction behavior shown above will be similar for exact solutions (even when they are close to the onion model in some specific sense). We used this non-exact model because it has suitable properties for investigation of averaging and backreaction.

Averaging in LRS class II spacetimes

11

0.0010 0.0005

2

4

6

8

l -0.0005 -0.0010 -0.0015 -0.0020 -0.0025

Fig. 3 Backreaction term

√1 3



 hσi2 − σ 2 +

1 3

(hθσi − hθi hσi) in the evolution equation

for shear depending on the averaging scale l and normalized by hσi· .

7 Averaged LRS dust class II equations One of the most important equations in cosmology is the evolution equation for the expansion scalar. independent



 In the averaged equation (25) we

have  variables θ2 and σ 2 . To obtain evolution equation for θ2 we multiply (8) by 2θ. Then we perform averaging and we obtain the equation 

2 · 1  3  2 3 2 3 θ − hθi (45) θ = − hθi − 4 hθi hσi − 8π hρi hθi + 3 3    



 3 2 + hθi − hθi θ2 − 4 θσ 2 − hθi hσi − 4π (hρθi − hρi hθi).

 In a similar way we derive an evolution equation for σ 2

σ2

·

 4 2 1 3 2 2 = − √ hσi − hθi hσi − 2 hEi hσi + hθi hσi − θσ 2 3 3 3     2 4  2 2 3 hθi σ − hθi hσi − √ σ 3 − hσi − 3 3 −2 (hEσi − hEi hσi).

(46)



  Now, we also need equations for e.g. θ3 , θσ 2 or σ 3 (and of course an evolution equation for hρi, hEi and hσi given in Section 4). We could obtain these evolution equations by the same procedure. Thus we have an infinite number of equations for the correlation terms. Here we need to adopt an ansatz. For example we can consider a reasonable assumption that for a given order the correlation terms are negligibly small and we can truncate the hierarchy to obtain a finite set of equations. We can also assume that some terms are proportional to each other. In this approach the inhomogeneities are modeled by different relations for correlation functions. The question is what kind of spacetime may correspond to the given set of averaged equations. We have started with scalar equations characterizing a LRS class II dust spacetime (containing the LTB metric and its generalizations). By averaging we can not leave this class, instead we may end up

12

Petr Kaˇspar† , Otakar Sv´ıtek§

in a special subclass of LRS class II dust models. We performed averaging of the evolution equation for expansion (25) and of the evolution equations for different products of expansion, shear, density and electric part of the Weyl tensor. From this construction we can see that the averaged equations contain an averaged LTB model, but generally not e.g. the homogeneous LRS Bianchi cosmologies. In the above described approach we have evolution equations for averages of different powers and products of the expansion, shear, density and electric part of the Weyl tensor. The problem is that if we derive an evolution equation for the averaged nonlinear terms, then more complicated terms appear in the relevant correlation terms (as is evident in equations (45) and (46)). If we want to close the system of equations we need to effectively eliminate these higherorder terms. In general, we may express the ”unwanted terms” as a suitable function of the lower-order averaged terms whose evolution equation is known. This kind of ansatz does not need to make all higher-order correlation terms necessarily vanish, it only serves to close the system of averaged equations through the selected relations - these may be for example expressed as products of the averaged terms of lower order.

8 Conclusion We generalized Buchert’s equations for the LRS class II dust model. We used the property that this family is characterized only by scalars and we employed a similar technique for averaging. However, the system of averaged equations is not closed. Buchert considered the so-called scaling solutions [35] to close the system of equations. In our work, we first investigated the influence of backreaction on the expansion and shear scalars for an approximate LTB model which describes fluctuating radial inhomogeneities. Then we proposed an infinite system of equations which supplement the averaged equations for expansion. In this approach inhomogeneities are modeled by the form of the correlation terms. Finally, we discussed how to close the system of averaged equations.

Acknowledgments We would like to thank R. Sussman for a useful discussion. We would also like to express our gratitude to the referee for correcting and improving our paper. P.K. was supported by grants GAUK 398911 and SVV-267301. O.S. ˇ 14-37086G. acknowledges the support of grant GACR

Appendix In this appendix we will show the computations demonstrating that the averaged constraint equations (30) - (32) are preserved in time. We start with an unaveraged constraint (13)

Averaging in LRS class II spacetimes

13

and perform straightforward averaging without using the commutation rules which gives us

0 1  2 σ = √ θ0 − haσi . 3 3

(A.1)

Now we take time derivative of (A.1) and use the commutation rule (17). We obtain the following expression







  1  0 ·  (σ 0 )· − hθi σ 0 + θσ 0 = √ (θ ) − hθi θ0 + θθ0 3 2 − [h(aσ)· i − hθi haσi + hθaσi] . 3

(A.2)

Now we need to commute prime and dot derivatives. This is done applying the commutation rule (34)



   1 0 2

θσ − hθi σ 0 + θσ 0 = (A.3) (σ) ˙ 0 − √ σσ 0 − 3 3 D  E





  1 ˙ 0 − √2 σθ0 − 1 θθ0 − hθi θ0 + θθ0 − 2 [h(aσ)· i − hθi haσi + hθaσi] . √ (θ) 3 3 3 3



Next, we apply the unaveraged evolution equations for θ, σ and a ((8), (9) and (12)) and we obtain the following expression:  2 0  2 0 0  1 0



 2

2

θσ − θσ − E − √ σσ 0 − θσ − hθi σ 0 + θσ 0 = − √ σσ 0 − 3 3 3 3 3  



  1 0

  1 2 0 2

√ − θθ − 4 σσ 0 − 4π ρ0 − √ σθ0 − θθ − hθi θ0 + θθ0 (A.4) 3 3 3 3     2 1

2 1 1

− − haθσi + √ aσ 2 − √ aσ 2 − haθσi − haEi − hθi haσi + hθaσi . 3 3 3 3 3 From the above expression we can see that several terms cancel each other. Moreover, using constraint equations (13) and (14) to further simplify the above equation (A.4) it is straightforward to see that the left hand side is equal to the right hand side. This means that the constraint equation (30) does not change in time. In the same way it can be shown that the constraint equations (31) and (32) are preserved in time, too.

References 1. Zalaletdinov, R.M.,: Averaging out the Einstein equations, Gen. Rel. Grav. 24, 1015 (1992). 2. Zalaletdinov, R.M.,: Towards a theory of macroscopic gravity, Gen. Rel. Grav. 25, 673 (1993). 3. Korzy´ nski, M.: Covariant coarse graining of inhomogeneous dust flow in general relativity, Class. Quantum Grav. 27, 105015 (2010). 4. Brannlund, J., van den Hoogen, R., Coley, A.: Averaging geometrical objects on a differentiable manifold, Int.J.Mod.Phys. D19 1915-1923 (2010) 5. Buchert, T.: On Average Properties of Inhomogeneous Fluids in General Relativity: Dust Cosmologies, Gen. Rel. Grav. 32, 105 (2000). 6. Buchert, T.: On Average Properties of Inhomogeneous Fluids in General Relativity: Perfect Fluid Cosmologies , Gen. Rel. Grav. 33, 1381 (2001). 7. Wiltshire, D.L.: Cosmic clocks, cosmic variance and cosmic averages , New J. Phys., 9, 377 (2007). 8. Li, N., Schwarz, D.: Onset of cosmological backreaction, Phys. Rev. D 76, 083011 (2007).

14

Petr Kaˇspar† , Otakar Sv´ıtek§

9. Li, N., Schwarz, D.: Scale dependence of cosmological backreaction, Phys. Rev. D, 78, 083531 (2008). 10. Behrend, J., Brown, I. A., Robbers, G.: Cosmological backreaction from perturbations, JCAP 01, 013 (2008). 11. Clarkson, Ch., Ananda, K., Larena, J.: Influence of structure formation on the cosmic expansion, Phys. Rev. D 80, 083525 (2009). 12. Larena, J., Alimi, J. M., Buchert, T., Kunz, M., Corasaniti, P. S.: Testing backreaction effects with observations, Phys. Rev. D 79, 083011 (2009). 13. Ellis, G.F.R.: Dynamics of Pressure-Free Matter in General Relativity, J. Math. Phys. 8, 1171 (1967). 14. Stewart, J.M., Ellis, G.F.R.: Solutions of Einstein’s Equations for a Fluid Which Exhibit Local Rotational Symmetry , J. Math. Phys. 9, 1072 (1968). 15. van Elst, H., Ellis, G.F.R.: The covariant approach to LRS perfect fluid spacetime geometries, Class. Quantum Grav. 13 1099 (1996). 16. Sussman, R.: Back-reaction and effective acceleration in generic LTB dust models, Class. Quant. Grav. 28, 235002 (2011). 17. Chuang, C. H., Gu, J. A., Hwang, W. Y. P.: Inhomogeneity-induced cosmic acceleration in a dust universe, Class. Quantum Grav. 25, 175001 (2008). 18. Paranjape, A., Singh, T. P.: The possibility of cosmic acceleration via spatial averaging in Lemaˆıtre-Tolman-Bondi models, Class. Quantum Grav. 23, 6955 (2006). 19. R¨ as¨ anen S.: Backreaction in the Lemaˆıtre-Tolman-Bondi model, JCAP 0411, 010 (2004). 20. Bolejko, K., Andersson, L.: Apparent and average accelerations of the Universe, JCAP 10, 003 (2008). 21. Sussman, R.: On spatial volume averaging in Lemaˆıtre-Tolman-Bondi dust models. Part I: back reaction, spacial curvature and binding energy, arXiv:0807.1145 [gr-qc], (2008). 22. Sussman, R.: Quasi-local variables and scalar averaging in LTB dust models, AIP Conf. Proc. 1241, 1146-55 (2010). 23. Mattsson, T., Mattsson, M.: Exploiting scale dependence in cosmological averaging, JCAP02, 004 (2008). 24. Mattsson, T., Mattsson, M.: On the role of shear in cosmological averaging, JCAP 1010, 021 (2010). 25. Mattsson, T., Mattsson, M.: On the role of shear in cosmological averaging II: large voids, non-empty voids and a network of different voids, JCAP 1105, 003 (2011). 26. R¨ as¨ anen S.: Cosmological acceleration from structure formation , Int.J.Mod.Phys. D15, 2141-2146 (2006). 27. R¨ as¨ anen S.: Evaluating backreaction with the peak model of structure formation, JCAP 04, 026 (2008). 28. Paranjape, A., Singh, T. P.: Structure formation, backreaction and weak gravitational fields, JCAP 03, 023 (2008). 29. Lemaˆıtre, G.: L’Universe en Expansion, Ann. Soc. Sci. Bruxelles A53, 51-85 (1933). 30. Tolman, R.C.: Effect of Inhomogeneity on Cosmological Models, Proc. Nat. Acad. Sci. U.S.A. 20, 169-76 (1934). 31. Bondi, H.: Spherically Symmetric Models in General Relativity, Mon. Not. Roy. Astron. Soc. 107, 410 (1947). 32. Bolejko, K., Krasi´ nski, A., Hellaby, C., C´ el´ erier, M.N.: Structures in the Universe by exact methods: formation, evolution, interactions, Cambridge: Cambridge University Press, (2009). 33. Hellaby, C.: Modelling Inhomogeneity in the Universe, 5th International School on Field Theory and Gravitation, Cuiaba, Brazil, 20-24 April 2009, Proc. Sci. PoS(ISFTG) 005 (2009). 34. Biswas, T., Mansouri, R., Notari, A.: Non-linear structure formation and ’apparent’ acceleration: an investigation, JCAP 12, 0712 (2007). 35. Buchert, T., Larena, J., Alimi, J.M.: Correspondence between kinematical backreaction and scalar field cosmologies - the ‘morphon field’, Class. Quant. Grav. 23, 6379 (2006). 320, 1 (1997).