Effects of Nonlinear Inhomogeneity on the Cosmic Expansion with Numerical Relativity
Total Page:16
File Type:pdf, Size:1020Kb
week ending PRL 116, 251302 (2016) PHYSICAL REVIEW LETTERS 24 JUNE 2016 Effects of Nonlinear Inhomogeneity on the Cosmic Expansion with Numerical Relativity Eloisa Bentivegna1,2,* and Marco Bruni3 1Dipartimento di Fisica e Astronomia, Università degli Studi di Catania, Via Santa Sofia 64, 95123 Catania, Italy 2INFN, Sezione di Catania, Via Santa Sofia 64, 95123 Catania, Italy 3Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth PO1 3FX, United Kingdom (Received 7 December 2015; revised manuscript received 8 March 2016; published 24 June 2016) We construct a three-dimensional, fully relativistic numerical model of a universe filled with an inhomogeneous pressureless fluid, starting from initial data that represent a perturbation of the Einstein–de Sitter model. We then measure the departure of the average expansion rate with respect to this homogeneous and isotropic reference model, comparing local quantities to the predictions of linear perturbation theory. We find that collapsing perturbations reach the turnaround point much earlier than expected from the reference spherical top-hat collapse model and that the local deviation of the expansion rate from the homogeneous one can be as high as 28% at an underdensity, for an initial density contrast −2 of 10 . We then study, for the first time, the exact behavior of the backreaction term QD. We find that, for small values of the initial perturbations, this term exhibits a 1=a scaling, and that it is negative with a linearly growing absolute value for larger perturbation amplitudes, thereby contributing to an overall deceleration of the expansion. Its magnitude, on the other hand, remains very small even for relatively large perturbations. DOI: 10.1103/PhysRevLett.116.251302 Cosmology as a physical theory of the Universe was one wishes to correctly interpret the data that will be born soon after the formulation of general relativity one produced by the upcoming precision surveys [8,9].While hundred years ago [1], yet the extent to which relativistic some approaches have been introduced to estimate the role of nonlinearity may affect structure formation remains largely relativistic corrections in N-body simulations [10–15],the unexplored. With the increasing volume of cosmological only viable avenue to an exact computation of the systematic data and their precision, more sophisticated modeling is errors resulting from the omission of these effects is the required, and thus it is becoming timely to quantify these direct numerical integration of Einstein’s equation in the relativistic effects. The current theoretical framework for corresponding scenarios. Integrating the equations of general cosmology is based on three main ingredients: a homo- relativity, possibly coupled to stress-energy sources, is the geneous and isotropic Friedmann-Lemaître-Robertson- field of numerical relativity, a framework strongly motivated Walker (FLRW) background, relativistic perturbation theory by gravitational-wave-source modeling, but which has, over to describe fluctuations in the early Universe and at very large the years, developed in a number of parallel areas such as scale, and Newtonian methods, notably N-body simulations, cosmology, mathematical relativity, and modified gravity to study the evolution of fluctuations into the nonlinear [16]. Some of this work has already been aimed at studying regime of structure formation. Reconciling this framework inhomogeneous cosmologies [17–21]. While these numeri- with the observations requires the existence of dark compo- cal-relativity studies do not yet aspire to the level of realism nents, cold dark matter (CDM) and a cosmological constant achieved by N-body simulations [22,23], they are useful Λ or some other form of dark energy. The resulting standard test beds to quantify the relativistic effects of nonlinear cosmological model, ΛCDM, satisfies a vast class of inhomogeneity on the cosmic expansion. observational constraints, in particular the high-precision In this Letter, we integrate Einstein’s equation coupled measurements of the cosmic microwave background anisot- to an inhomogeneous irrotational pressureless fluid (dust) ropies [2]. However, the existence and nature of these dark with a three-dimensional density profile and no continuous constituents are one of the most debated topics not only symmetries. We choose initial data corresponding to a in modern cosmology, but also in theoretical physics. perturbed Einstein–de Sitter (EdS) model, i.e., a flat FLRW One aspect that has been the subject of intense debate is model with dust, with the aim of measuring, with no the question whether nonlinear relativistic “backreaction” approximations, the departures of the fully nonlinear effects due to formation of structures may play an important numerical solution from the idealized FLRW background role in the average cosmic expansion [3–7]. and its perturbations. On the numerically generated space- Quantifying the systematic errors involved in the differ- times, we measure a number of local and average properties ent modeling approximations, such as the use of Newtonian of cosmological interest, such as the growth of over- gravity for structure formation, is a crucial undertaking if densities and the formation of voids, the inhomogeneous 0031-9007=16=116(25)=251302(6) 251302-1 © 2016 American Physical Society week ending PRL 116, 251302 (2016) PHYSICAL REVIEW LETTERS 24 JUNE 2016 2 and average expansion rate, and the backreaction term γ¯ij ¼ aðtÞ δij, where the scale factor aðtÞ is a solution of defined in the averaging framework [3]. The main results of Friedmann’s equations this study are that (i) those perturbations that are large a_ 2 8πρ¯ ä 4π enough to collapse stop partaking in the cosmic expansion ¼ ¼ − ρ¯; ð Þ (i.e., reach the “turnaround” point) much earlier than a2 3 a 3 1 expected from a spherical top-hat collapse model with the same initial density contrast, (ii) locally, the effects of where the dot represents a time derivative, and we denote nonlinear inhomogeneities can be substantial, leading to a the EdS-background quantities with an overbar. The matter ρ¯ ∼ a−3 departure from the average expansion rate of over 28% at continuity equation gives for the background density. ρ the underdensities, and (iii) the average expansion rate is Starting from the inhomogeneous density , one can define δ ¼ðρ − ρ¯Þ=ρ¯ hardly affected by the inhomogeneities, with a backreaction the density contrast ; its growth in the term which is never larger than 10−8. synchronous-comoving gauge is governed, at first order, by Method.—We integrate Einstein’s equation and the fluid 3 3 δ00 þ δ0 − δ ¼ 0 ð Þ conservation equation using a variant of the Baumgarte- 2a 2a2 . 2 Shapiro-Shibata-Nakamura formulation [24–26], along with the Wilson formulation for the hydrodynamical The system of (1) and (2) is then solved by system [27] and the conformal transverse traceless formu- t 2=3 aðtÞ¼a ; ð Þ lation for the Einstein constraints [28,29], an approach i t 3 already used in cosmological settings [30–32]. We choose i −3=2 to represent the spacetime in the synchronous-comoving δðtÞ¼δþaðtÞþδ−aðtÞ ; ð4Þ gauge [33], popular in cosmological perturbation theory [22], which corresponds to the Lagrangian coordinates of where δþ and δ− are the so-called growing and decaying the observers at rest with the matter. modes. We will use these expressions below as a consistency To integrate this system, we use the Einstein Toolkit [34], check in the small-perturbation regime. a free, open-source community infrastructure for numerical Another useful framework is that of cosmological aver- relativity. In particular, we use the MCLACHLAN code [35] aging [3], where Einstein’s equation is reduced from a set of for the evolution of the gravitational variables, the CARPET partial differential equations for the fields to a set of ordinary [36] package for handling adaptive mesh refinement, and differential equations in time for some of their averages the multigrid elliptic solver CT_MULTILEVEL [37] to gen- over a given spatial region D. Defining its volume as erate initial data; this is then coupled to a new module Z 3 pffiffiffi 3 which evolves the hydrodynamical equations. All equa- aD ¼ γd x; ð5Þ tions are discretized using fourth-order finite differencing. D The Einstein Toolkit is routinely used for simulations in where γ is the determinant of the spatial metric γij, one finds relativistic astrophysics, and passes a variety of tests [34]. that the average scale factor aD satisfies a system similar to Likewise, as will be presented elsewhere, the new module Friedmann’s (1), and in particular that correctly reproduces several exact cosmological models ä 4π M Q with varying degrees of inhomogeneity. All results pre- D ¼ − D þ D ; ð Þ 3 6 sented are convergent at the correct rate as the grid spacing aD 3 aD 3 is decreased, and we use this fact to extrapolate the where continuum solution of the evolution system and to estimate the error bars resulting from its numerical integration at Z pffiffiffi 3 finite resolution. These are the quantities that appear in MD ¼ γρd x ð7Þ all plots. D Perturbations and averaging.—We recall two appro- 2 Q ¼ ðhK2i − hKi2 Þ − 2hA2i : ð Þ aches commonly used to solve the evolution system D 3 D D D 8 approximately, so that we can compare our solution to K K ≡ −γ_ =2 these schemes and check that we obtain the correct Here is the trace of the extrinsic curvature ij ij , 2 ij behavior in the appropriate regime. A ¼ AijA =2, Aij is the traceless part of Kij, and h·iD For irrotational dust in the synchronous-comoving gauge, denotes the average of a field over D. Note that −K represents the line element can be written (with no loss of generality the local expansion rate, and in the FLRW background H ¼ 2 2 i j ¯ [33])asds ¼ −dt þ γijdx dx , where γij is the spatial −K=3 is the Hubble parameter.