Low-redshift effects of local structure on the Hubble parameter in presence of a cosmological constant

In order to estimate the effects of a local structure on the Hubble parameter we calculate the low-redshift expansion for H(z) and δHH\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\delta H}{H}$$\end{document} for an observer at the center of a spherically symmetric matter distribution in the presence of a cosmological constant. We then test the accuracy of the formulas comparing them with fully relativistic non-perturbative numerical calculations for different cases for the density profile. The low-redshift expansion we obtain gives results more precise than perturbation theory since it is based on the use of an exact solution of Einstein’s field equations. For larger density contrasts the low-redshift formulas accuracy improves respect to the perturbation theory accuracy because the latter is based on the assumption of a small density contrast, while the former does not rely on such an assumption. The formulas can be used to take into account the effects on the Hubble expansion parameter due to the monopole component of the local structure. If the H(z) observations will show deviations from the ΛCDM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda CDM$$\end{document} prediction compatible with the formulas we have derived, this could be considered an independent evidence of the existence of a local inhomogeneity, and the formulas could be used to determine the characteristics of this local structure.


Introduction
The standard cosmological model is based on the assumption that the Universe is homogeneous and isotropic on sufficiently large scales. Nevertheless local observations could be strongly affected by local structure as shown for example in [1], and it is important to study its effects. The analysis of luminosity density data [2] has provided some strong experimental evidence supporting the existence of local inhomogeneities, but it would be important to confirm it using a e-mail: aer@physics.uoc.gr another observable such as the baryonic acoustic oscillations (BAO) measurements [3][4][5][6][7][8]. The BAO scale allows in fact to determine the expansion rate of the Universe H (z) independently from the luminosity distance and as such provides an important source of information as regards our Universe. If the H (z) estimations obtained from BAO observations data will show deviations from the C DM predictions this could be considered an independent evidence of the existence of local inhomogeneities. This motivates the calculation of a low redshift formula for H (z), able to take into account the effects of inhomogeneities which cannot be fully modeled with perturbation theory, as some of the inhomogeneities found for example in [2]. A low-redshift expansion based on the use of exact solutions of Einstein's equations is in fact valid also for large values of the density contrast or of the gravitational potential.
The effects of a local inhomogeneity on cosmological observations have been studied already for different cases [1, such as the equation of state of dark energy or the luminosity distance [18,19,23]. It has been shown for example that the value of the cosmological constant could be affected significantly by the presence of local inhomogeneity seeded by primordial curvature perturbations [9], which could also lead to the wrong conclusion of a varying equation of state for dark energy while only a cosmological constant is present [19]. The origin of these effects is that spatial inhomogeneities can change the energy of the propagating photons, contaminating the cosmological redshift due to the Universe expansion and consequently introducing some errors in the estimation of parameters based on cosmological models which ignore the effects of the inhomogeneities.
In this paper we will focus on the low-redshift effects of inhomogeneities on the Hubble expansion parameter, adopting an analytical approach based on the use of an exact solution of Einstein's equations to model the local structure. We first derive the redshift expansion of the geodesics equations and use it to obtain the expansion of H (z). We then compute a formula for δ H H , the relative difference between the C DM and the inhomogeneous case. Finally we compare the formulas with the numerical calculations based on the integration of the Einstein's equation and the geodesics equations, finding good agreement. We also check that the low-redshift expansion formulas are more precise than the perturbative calculation, especially when the density contrast is larger.

Modeling the local Universe
We use the LTB solution to model the monopole component of the local structure [31][32][33][34][35] where R is a function of the time coordinate t and the radial coordinate r , E(r ) is an arbitrary function of r , and R (t, r ) = ∂ r R(t, r ). The Einstein equations imply that where M(r ) is an arbitrary function of r ,Ṙ = ∂ t R(t, r ), and we choose a system of units in which c = 8π G = 1.
To compute H (z) we need to solve the radial null geodesics [36] and then we substitute in the formula for the Hubble parameter in a LTB space [37,38] where The analytical solution can be derived [39,40] introducing a new coordinate η = η(t, r ), and new functions ρ 0 (r ) and k(r ) given by We will adopt, without loss of generality, the coordinate system in which ρ 0 (r ) is a constant, the so called FLRW gauge. We can then express Eq. (2) in the form The coordinate η, which can be considered a generalization of the conformal time in a homogeneous FLRW Universe, is defined implicitly by Eq. (10). The relation between t and η is obtained by integrating Eq. (10) and is given by [26] t where t b (r ) is a functional constant of integration called the bang function, since it corresponds to the fact that in these models the scalar factor can vanish at different times at different locations. We will consider models with t b (r ) = 0. The solution of Eq. (13) can then be written in the form where ℘ (x; g 2 , g 3 ) is the Weierstrass elliptic function and In terms of η and a(η, r ) the radial null geodesics and the Hubble parameter are given by [12] where The function G(η, r ) has an explicit analytical form, making the above geodesics equations particularly suitable for a lowredshift expansion.

Low-redshift expansion of the Hubble parameter H(z)
In order to obtain a low-redshift formula for the Hubble parameter we expand the function k(r ) as We also need an expansion for t (η, r ), which can be obtained from the exact solution for a(η, r ) according to where we defined the function t 0 (r ) by Using the properties of the Weierstrass elliptic functions ℘, ζ , and σ [41] we can compute the integral in Eq. (14), obtaining where ℘ is the derivative of the Weierstrass elliptic function ℘. ℘ −1 , ζ , and σ are defined by the equations The use of the exact expression for t (η, r ) improves the accuracy for the expansion for the geodesics with respect to previous calculations [23], which were based on a perturbative expansion of t 0 (r ), rather than the use of the exact value. Now we can find the low redshift Taylor expansion for the geodesic equations [19], and then we calculate the Hubble parameter. We expand the solution of the geodesic equations according to After substituting the above expansion in the geodesic equations we can map the solution of the system of differential equations into the solution of a system of algebraic equations for the coefficients of the expansion. Here we give the formulas for the case in which k 0 = 0, which is enough to understand qualitatively the effects of the inhomogeneity. The term k 0 corresponds in fact to the homogeneous component of the curvature function, which in the absence of inhomogeneities is simply the curvature of a FLRW model and as such is not associated to any new physical effect not already known from standard cosmology. For the geodesics we get where M , , T 0 , and K n are dimensionless quantities given by [16] We have also used the following definitions: where ζ is the Weierstrass zeta function. As we can see in the above formulas the effects of the inhomogeneity start to show, respectively, at first order for η(z) and second order for r (z).
In order to obtain a formula for the Hubble parameter as a function of the redshift we need to substitute Eqs. (17) and (18) in Eq. (19), After expanding up to second order in z we get 1 We plot t 0 (0) as a function of K 1 . This is the quantity determining the leading order effect for δ H H as shown in Eq. (51) The procedure to reduce the analytical formula to this form is rather complicated since it involves the expression wherever possible of all the intermediate terms of physically meaningful quantities using the properties of the Weierstrass elliptic functions [41]; we give more details in Appendix A.
As we can see from the first order coefficient H 1 , the leading order in t 0 (0) determines the sign of the correction with respect to the homogeneous case, and for this reason we plot t 0 (0) as a function of K 1 in Fig. 1. At second order we have a more complicated dependency for H 2 , which involves also K 2 and t 0 (r ).
We can easily interpret the linear behavior shown Fig. 1, applying the chain rule for the derivative

Relative difference of H(z) with respect to the homogeneous case
For a flat FLRW solution the expansion rate is given by Since we want to compare the inhomogeneous case with the flat FLRW case we define the relative difference as where we denote by H LT B (z) the expansion rate defined in Eq. (42).
We can now expand the above expression at low redshift to get δ H (z) It is straightforward to verify that in the homogeneous limit, when K 0 = K 1 = K 2 = 0, as expected, δ H 1 = δ H 2 = 0, since t 0 (0) = t 0 (0) = 0. From the expression for δ H 1 H the crucial role played by t 0 (0) is clear, which determines the sign of the relative difference at leading order in the redshift, and according to Eqs. (46) and (47) it is proportional to K 1 .
Using Eq. (47) we can also see that for a fixed K 0 the sign of K 1 determines the sign of δ H 1 H at very low redshift when the second order contributions can be neglected. This is in good agreement with Figs. 3 and 4, which show an approximate linear behavior with opposite slope, corresponding to different values of K 1 .

Testing the accuracy of the formulas
In order to test the accuracy of the analytical formulas and compare it with perturbation theory we perform numerical calculations using LTB solutions corresponding to the density contrasts shown in Fig. 2. Using large density contrasts we can test the limit of the perturbation theory results and compare them with the low-redshift expansion.
As seen in Figs. 3 and 4 the analytical formula for the relative difference of H (z) with respect to the homogeneous Fig. 2 The density contrast δ = δρ ρ b is plotted as a function of the redshift for three different inhomogeneities profiles modeled by LTB solutions Fig. 3 The relative difference with respect to the homogeneous case δ H (z) H is plotted as a function of the redshift for three different density contrasts. The colors correspond to the density contrasts in Fig. 2. The solid lines are for the numerical calculation, the dashed lines are for the analytical formulas, and the dot-dashed lines are for the perturbation theory result The colors correspond to the density contrasts in Fig. 2. As can be seen the low-redshift formula is always better, especially for larger density contrasts case is in good agreement with δ H (z) H obtained by integrating numerically the geodesics and background equations and is more accurate than the perturbation theory. From Fig. 4 in particular we can see that the agreement with the exact numerical calculation for the redshift expansion is in general better than that of the perturbation theory result. For larger density contrasts, as expected, the perturbative calculation is increasingly inaccurate while the red-shit expansion is still in good agreement with the exact numerical result, since it does not rely on the assumption of a small density contrast. This is due to the fact that perturbation theory is based on the assumption that δρ/ρ 1, while the low-redshift expansion is based on an exact solution of the Einstein equations.

Conclusions
We have derived a low-redshift analytical formula for the Hubble parameter for an observer at the center of a spherically symmetric matter distribution, using an exact solution of Einstein's field equations. Such a formula is in good agreement with numerical calculations and is more accurate than the perturbation theory result, especially for large density contrasts. This is due to the fact that perturbation theory is based on the assumption that δρ/ρ 1, while the lowredshift expansion is based on an exact solution of the Einstein equations.
If the H (z) observations will show deviations from the C DM predictions compatible with the formulas we have derived, this could be considered an independent evidence of the existence of a local inhomogeneity, and the formulas could be used to determine the characteristics of this local structure. While the expansion for r (z) depends on our choice of radial coordinate, the formulas for the H (z) are independent of it, since both H and the redshift are physically observable quantities and as such are independent of the gauge choice. Since in the derivation of the formulas the inhomogeneity profile is determined by the coefficients of the expansion of the function K (r ), there is still some coordinate dependency left in the way we parameterize the inhomogeneity. While our choice of gauge is quite natural, since in the F RW gauge, corresponding to ρ 0 (r ) = const, the radial coordinates reduces to the radial F RW comoving coordinate in the limit in which the function k(r ) goes to zero, a totally coordinate independent formula would still be preferable. For this reason in the future it will be interesting to derive formulas which are completely independent of the choice of the coordinate system, parameterizing the inhomogeneity in terms of the redshift expansion of the density ρ(z), which is a physical observable, instead of using the expansion of K (r ).