# On Backreaction in Newtonian cosmology

On Backreaction in Newtonian cosmology Abstract We clarify that a result recently stated by Kaiser is contained in a theorem of Buchert and Ehlers that is widely known for its main result: that there is no global kinematical backreaction in Newtonian cosmology. Kaiser cites this paper, re-derives parts of the theorem, but incompletely restates its content. He makes further claims, which cannot be proven beyond the limited context of Newtonian cosmology. We also discuss recent papers of Rácz et al. and Roukema, who claim the existence of global backreaction within the Newtonian framework. cosmological parameters, cosmology: theory, dark energy 1 AVERAGE PROPERTIES OF NEWTONIAN MODELS In a recent paper Kaiser (2017b, hereafter K17) considers, as do Buchert & Ehlers (1997, hereafter BE – see also Ehlers & Buchert 1997), the Euler–Poisson system in the fluid approximation for a ‘dust’ matter model in the mean field approximation of Newtonian gravity, equations (2) and (3) in K17 [cf. point (B) below]. Buchert and Ehlers performed spatial averaging of the kinematical scalars of the system to obtain the following general expansion law, given in BE, equation (B4):   \begin{eqnarray} {3{{\ddot{a}}_{\cal D} \over a_{\cal D}} + 4 {\pi} G {M_{\cal D}\over a_{\cal D}^3 } - \Lambda = :{\cal Q}_{\cal D} } \nonumber \\ {\quad=\langle {\nabla }\cdot [{\boldsymbol u} ({\nabla } \cdot {\boldsymbol u}) - ({\boldsymbol u} \cdot {\nabla }) {\boldsymbol u} ]\rangle _{\cal D} - {2\over 3} \langle {\nabla } \cdot {\boldsymbol u} \rangle _{\cal D}^2 }\nonumber \\ {\quad\quad+2\left(\Omega ^2 - \Sigma ^2 \right) + 2\left(\Omega _{ij}\langle \hat{\boldsymbol \omega }_{ij}\rangle _{\cal D} - \Sigma _{ij}\langle \hat{\boldsymbol \sigma }_{ij}\rangle _{\cal D}\right), } \end{eqnarray} (1)where $$a_{\cal D} \propto V_{\cal D}^{1/3}$$ denotes the volume scalefactor of any compact averaging domain $$\cal D$$ with volume $$V_{\cal D}$$, $$\langle \cdot \rangle _{\cal D}: = 1/ V_{\cal D} \int _{\cal D}{\mathrm{d}}V\cdot$$ the averaging operator, $$\boldsymbol u$$ the peculiar-velocity field (with its shear $$\hat{\boldsymbol \sigma }_{ij}$$ and vorticity $$\hat{\boldsymbol \omega }_{ij}$$), defined with respect to a global reference background flow, usually taken to be a Hubble flow (and Kaiser recalls correctly that this reference flow is a priori not fixed to be a Hubble flow). Omitting global shear Σ and global vorticity Ω of the model universe, we arrive at the result of Kaiser (i.e. the kinematical backreaction term, equation (9), in K17).1 Kaiser (K17, appendix A) re-derives the result of BE stating that ‘we show how this may be obtained directly’. Inspection of his derivation shows, however, that it is essentially identical to the line of reasoning of BE, but less straightforward. Kaiser obtains a linear, background-dependent term that (as he now appreciates, see footnote 1) is zero. The direct derivation of BE does not lead to such a zero term. Follow-up papers and reviews also provided complementary derivations and thorough discussions of the properties of the kinematical backreaction term (e.g. Buchert 2000a; Buchert 2008, section 3; Buchert & Räsänen 2012, section 2.5). The result (1) is valid for all compact domains from an infinitesimal domain up to any scale, e.g. up to the boundary of a compact model universe, and the second term in the second line of (1) is there on all scales. Kaiser’s paper titled Why there is no Newtonian backreaction is an improvement over the title There is no kinematic backreaction of (Kaiser 2017a), but still overstates and misleadingly interprets the content of the BE theorem for the following reasons: What BE showed, and what Kaiser re-derives, is that there is no global kinematical backreaction in Newtonian cosmology for boundary-free compact model universes.2 There is, however, kinematical backreaction, describing cosmic variance of the deviations from the assumed reference background, as Kaiser correctly states, and which is explored in the paper by Buchert et al. (2000). Thus, the suggestion that there is ‘no Newtonian backreaction’ per se is not supported by the text. Kinematical backreaction is present in Newtonian models in the interior of, e.g. a three-torus model3 (and the expectation values of the peculiar-velocity invariants are in general non-zero on regional domains). The term $${\cal Q}_{\cal D}$$ in (1) can be expressed in terms of Minkowski Functionals (Buchert 2000a, section 1.11), (Buchert 2008, section 3.1.2). These functionals describe morphological properties of the density distribution and depend on all correlation functions of a density field arising out of smoothing a point-like distribution. Consequently, even small-amplitude interior backreaction terms are important in characterizing the cosmic web. Kaiser, in K17, section 2, leaves the impression that a more elementary approach using discrete particles leads to the same result as that obtained from the Euler–Poisson system for a dust fluid in the mean field approximation. This impression comes from the comparison of the Newtonian equation of motion for particles with the Euler equation (while K17 re-derives the main result, equation 9, within a fluid mechanical approach, section 3 and appendix A). Such an approach has been discussed by Ellis & Gibbons (2014), who show4 that a set of central configuration constraints has to be also satisfied. It should be clear that the general link between a particle approach and the governing fluid equations is more subtle: Starting from discrete particles, coarse-graining the Klimontovich density in phase space and taking velocity moments leads to extra anisotropic multistream stresses modifying the Euler equation, and deviations from mean field gravity modifying the Poisson equation (see e.g. Buchert & Domínguez 2005) not done in K17 and BE. This also leads to extra backreaction terms.5 There is no ‘significant advance over the approach followed e.g. in BE’ (K17): purely space-dependent arguments already assume the vanishing of the above terms.6 Only including these terms would render the result more general.7 Kaiser omits a discussion of the uniqueness of solutions to the Poisson equation. This problem, discussed in BE, is crucial to Newtonian cosmology since we have to specify boundary conditions. BE considered various possibilities, including cases where kinematical backreaction does not globally vanish (Charlier-type models), among them the possibility of introducing a background and deviations thereof defined on a three-torus. BE also discuss why the latter is the only relevant possibility for Newtonian cosmological models (apart from Charlier-type models without a homogeneity scale). To show this, we recall that the Poisson equation does not uniquely define the potential. It is unique if we require square integrability of the potential in the case of $$\mathbb {R}^3$$. [The sum in K17, equation (5), is in general divergent, a problem of which Newton and Einstein were well aware (see e.g. Ellis & Gibbons 2014, section 1.1, and references therein).] The possibility of a square integrable potential over $$\mathbb {R}^3$$ is, however, irrelevant for cosmology. The potential is also unique if the spatial average of the source vanishes. Introducing a background and deviations thereof allows one to consider periodic peculiar-velocity fields and density contrasts with a vanishing average of the source of Poisson’s equation for a periodic peculiar-potential (a three-torus model). This architecture, being identical to that of Newtonian N-body simulations, implies that the peculiar-potential exists and is unique apart from the addition of a solution of the Laplace equation, whose only periodic solutions are spatially constant. This constant term can be removed by the translation invariance of Newton’s equations and set to zero without loss of generality. In general relativity (GR), these problems do not arise and a model universe (that is, in general, background-free) is not restricted to obey this global constraint. In addition, the kinematical backreaction variable couples to the averaged intrinsic curvature (section 3.2.2 ff in Buchert 2008), removing the ‘Newtonian anchor’ (section 3.2.3 in Buchert 2008; section 2.5 in Buchert 2000a), which enforces the vanishing of structure averages globally on an assumed background. This is the reason why the globally vanishing backreaction in Newtonian cosmology is ‘by construction’: Kaiser misinterprets this insightful remark and notes that ‘... $${\cal Q}_{\cal D}$$ tends to zero very rapidly in the limit of large volumes regardless of whether the structure is assumed to be periodic’. As we pointed out above, periodicity is a necessary element of the architecture of Newtonian models applied to cosmology; by the very definition of $${\cal Q}_{\cal D}$$, the term decays rapidly with volume, but it is crucial that it has to vanish exactly on the periodicity scale. We note the importance of a non-vanishing but small $${\cal Q}_{\cal D}$$ for the evolution of averaged scalar curvature in the GR context, see the recent paper by Bolejko (2017a), who illustrates this by employing exact GR solutions. Scale-dependent deviations from an assumed background are unavoidable.8 2 CONCLUDING REMARKS AND DISCUSSION We conclude that there is no new result in K17. Known for 20 yr, the result of BE is key to Newtonian cosmology, and is uncontroversial among researchers working on backreaction problems. However, Kaiser’s re-statement that Newtonian models cannot lead to non-vanishing backreaction is correct globally, and relevant for the rebuttal of recent work by Rácz et al. (2017) on backreaction calculated within the Newtonian framework. We therefore discuss this latter work, as well as a recent complementary work by Roukema (2017), and we point out why the conclusions drawn from these models on the possibility of global backreaction contradicts the result of BE. Since kinematical backreaction is non-vanishing in the interior of a Newtonian simulation, and attains substantial values by going to smaller but still cosmologically relevant scales (Buchert et al. 2000), it is unsurprising that Rácz et al. (2017) re-discover expansion variance by implementing a multiscale volume partitioning in N-body simulations. These toy-simulations may be interpreted to illustrate the effect on differential expansion properties, and only on a restricted range of scales. It is clear, however, and a consequence of the theorem (1), that there cannot be backreaction in these simulations globally. The key-issue to understand this fact is the non-local nature of the gravitational interaction. Rácz et al. (2017) adopt assumptions on individual domains of the global model universe. In their paper, the assumption of vanishing backreaction within a spherical idealization of domains is adopted. The non-locality of gravitation that is reflected by global constraints (boundary conditions in Newtonian cosmology and propagating spatial constraints in GR 3+1–foliated models) assures, however, that any property of the regional domain is ‘known’ to the rest of the model universe. Thus, any approximate assumption made must be rendered consistent with the whole model universe. Otherwise, an error is introduced that – as in the case of Rácz et al. (2017) – results in a fictitious global contribution of backreaction. To avoid this error, the generic situation has to be implemented, accounting in particular for backreaction on regional domains (i.e. deviations from spherical symmetry). It is of equal importance that the domains of the volume partition are be joined properly to obtain consistency with the global constraint. Furthermore, although an evolution of the background away from a homogeneous Hubble flow is a generic property of GR cosmologies on any scale,9 it lacks a physical basis within the Newtonian framework. The separate universe conjecture (Dai, Pajer & Schmidt 2015, and references therein), i.e. making restrictive assumptions on regional domains and considering the region to evolve separated from the environment, has to be put under careful scrutiny, especially in non-linear modelling attempts. Another example for an uncontrolled implementation of this conjecture is the recent paper by Roukema (2017) that suggests an alternative explanation for structure-emerging global backreaction within the Newtonian framework. While Roukema (2017) uses generic properties on subdomains, in particular including backreaction on regional domains,10 he implements the assumption of ‘silent virialization’11 by stabilizing collapsed domains with the consequence of obtaining an excess global expansion over the assumed background expansion. However, virialization in collisionless Newtonian systems can be modelled exactly by including velocity dispersion, as is explained in point (B) of Section 1. This cannot lead to a global backreaction, as explained in footnote 7. Furthermore, it does not help to appeal to the algebraic similarity between the Newtonian backreaction model of Buchert et al. (2000) and the corresponding relativistic model of Buchert, Nayet & Wiegand (2013) through the claim that ‘virialization’ induces an excess of expansion in a ‘relativistically realistic’ (Roukema 2017) situation, while the actual realization of the model is still confined to the Newtonian setting. Also in relativistic N-body simulations (that include virialization), the spatial GR constraints and their propagation has to be controlled, representing what in Newtonian theory are the boundary conditions. These simulations may also face the constraints discussed here for Newtonian models, if periodic boundary conditions on deviations from a fixed background model are implemented. A promising route to background-free investigations of the average properties of GR universe models is furnished by multiscale volume partitions that have been investigated during the last decade (a comprehensive list of works may be found in recent reviews, e.g. in Buchert & Räsänen 2012, Wiltshire 2014). In particular, the exact volume partitioning formulas for GR models given in Buchert & Carfora (2008) and Wiegand & Buchert (2010) can be profitably used to demonstrate the consequences of any regional assumption made. (The present debate motivated us to explicitly demonstrate these consequences for an exact volume partition of Newtonian models in a forthcoming paper.) It may well be possible to model regional domains of the real Universe using Newtonian equations, but this does not provide justification to model the entire observable Universe by a single Newtonian solution. Kaiser’s remarks appear to assume that small regions combine into a single cosmology in a trivial way (i.e. that the ‘background’ of all small regions should be identical). A kaleidoscopic modelling of the Universe by consecutive Newtonian domains restricts the generality substantially, for example, by calculating the average Ricci curvature of the model universe. Kaiser’s discussion of the relevance of the BE theorem for GR constitutes an opinion. More efforts to justify that the Universe is everywhere close to the same background model have been invested by Green & Wald (2014). However, their assumptions are too restrictive to apply to cosmological backreaction: see (Buchert et al. 2015) and references therein, which provide an extensive discussion of the physics of cosmological backreaction. Finally, quoting statements of early papers that pioneered the subject risks being anachronistic if developments in the literature over the subsequent 15 yr are not taken into account. (The remarks cited by Kaiser are still agreed upon today, but nuanced.) Acknowledgements This work is part of a project that, in the final stage, has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement ERC advanced grant 740021–ARTHUS, PI: T. Buchert). Thanks to Jürgen Ehlers in memoriam with whom I share the results of BE, and thanks to Nick Kaiser, Martin Kerscher, Pierre Mourier, Jan Ostrowski, Boud Roukema, Syksy Räsänen and David Wiltshire for correspondence, discussions and feedback, as well as to an anonymous referee for an insightful and very careful report. 1 The three corrections (1a), (1b) and (1c) in Buchert (2017a) of Kaiser’s re-derivation of part of BE’s results in Kaiser (2017a) have now been included (apart from the erroneous factor 3/2 in K17, equation (9), which should read 2/3) after consulting (Buchert 2017a) (N. Kaiser, private communication), except that Kaiser still claims that his derivation is more general. 2 We shall explain this in point (C). Note that kinematical backreaction also vanishes for spherical regions embedded into a reference background Buchert, Kerscher & Sicka 2000, Buchert 2011, section 7.2. This latter provides a compact proof of the Newton iron sphere theorem. Kaiser quotes instead the work of Einstein & Straus (1945), which is out of context. That backreaction vanishes for a single spherical region is well known. 3 Kaiser confuses the scale-dependent volume expansion rate $$H_{\cal D}$$ with the background expansion rate H. It is trivial that ‘There is no Newtonian backreaction on a(t) from structure’ (K17), but there is Newtonian backreaction on $$a_{\cal D}(t)$$ from structure, except on the global scale where $$a_{\cal D}(t)$$ reduces to a(t). Kaiser discusses how the equation for a(t) can be chosen to recover the usual Newtonian equations for point particles, but this does not imply that the average motion of the point particles follows the same a(t). This dismisses the non-local aspect of averaging contained in $$a_{\cal D}(t)$$. 4 Kaiser now cites this paper without pointing out the limitations explained by Ellis & Gibbons (2014). 5 In general relativity, the averaged equations have been given by including isotropic stresses (Buchert 2001). 6 Recall that in situations of multiple streams, space-dependent variables cease to be functions and we have to move to a description in phase space. 7 Keeping the mean field assumption, it is straightforward to include velocity dispersion: Taking velocity moments of the phase-space density, the Euler–Jeans equation replaces the Euler equation featuring multistream stresses Πij. The additional source, $$-\frac{1}{\varrho }\frac{{\partial} }{{\partial} x_j}\Pi _{ij} =: - \Psi _i$$, results in a divergence term, $$- \langle {\partial} \Psi _k / {\partial} x_k \rangle _{\cal D}$$, adding to the kinematical backreaction term in (1), and leaving the main result of BE on the global vanishing of the kinematical backreaction term unchanged. There is no such term in Kaiser’s derivation. 8 In GR, this holds globally and can be traced back to first principles like the non-conservation of intrinsic curvature (equation (13) in Buchert 2000b, Buchert (2008) and Buchert & Carfora 2008). 9 The quantitative importance of the backreaction effect in GR is a result of the observation that the average model evolves away from a predefined FLRW (Friedmann–Lemaître–Robertson–Walker) background due to the non-conservation of curvature (see footnote 8). A scale-dependent ‘background’ (affected by structure formation) has been investigated in the GR context (Roy & Buchert 2012), together with global instability properties of the FLRW model (Roy et al. 2011). 10 Roukema (2017) employs the analytical backreaction model of Buchert et al. (2000) that is based on exact average properties, and an analytical model for the evolution of inhomogeneities that has earlier been well-tested against N-body simulations. 11 The word ‘silent’ here refers to the separate universe conjecture rather than to so-called ‘silent model universes’ in GR, (cf. Bolejko 2017b, and references therein) that have vanishing magnetic part of the Weyl tensor. Note that even ‘silent model universes’ do not obey the separate universe conjecture, since the spatial GR constraints and their consistent propagation have to be obeyed, which again reflects the non-locality of the gravitational interaction. REFERENCES Bolejko K., 2017a, J. Cosmology Astropart. Phys. , 06, 025 CrossRef Search ADS   Bolejko K., 2017b, preprint (arXiv:1707.01800v3) Buchert T., 2000a, JGRG , 9, 306, (arXiv:gr-qc/0001056) Buchert T., 2000b, Gen. Rel. Grav. , 32, 105 CrossRef Search ADS   Buchert T., 2001, Gen. Rel. Grav. , 33, 1381 CrossRef Search ADS   Buchert T., 2008, Gen. Rel. Grav. , 40, 467 CrossRef Search ADS   Buchert T., 2011, Class. Quant Gravity , 28, 164007 CrossRef Search ADS   Buchert T., 2017a, preprint (arXiv:1704.00703v1) Buchert T., Carfora M., 2008, Class. Quant Gravity , 25, 195001 CrossRef Search ADS   Buchert T., Domínguez A., 2005, A&A , 438, 443 CrossRef Search ADS   Buchert T., Ehlers J., 1997, A&A , 320, 1 (BE) Buchert T., Räsänen S., 2012, Ann. Rev. Nucl. Part. Sci. , 62, 57 CrossRef Search ADS   Buchert T., Kerscher M., Sicka C., 2000, Phys. Rev. D , 62, 043525 CrossRef Search ADS   Buchert T., Nayet C., Wiegand A., 2013, Phys. Rev. D , 87, 123503 CrossRef Search ADS   Buchert T.et al.  , 2015, Class. Quant Gravity , 32, 215021 CrossRef Search ADS   Dai L., Pajer E., Schmidt F., 2015, J. Cosmology Astropart. Phys. , 10, 059 CrossRef Search ADS   Ehlers J., Buchert T., 1997, Gen. Rel. Grav. , 29, 733 CrossRef Search ADS   Einstein A., Straus E. G., 1945, Rev. Mod. Phys , 17, 120 CrossRef Search ADS   Ellis G. F. R., Gibbons G. W., 2014, Class. Quant Gravity , 31, 025003 CrossRef Search ADS   Green S., Wald R. M., 2014, Class. Quant Gravity , 31, 234003 CrossRef Search ADS   Kaiser N., 2017a, preprint (arXiv:1703.08809v1) Kaiser N., 2017b, MNRAS , 469, 744 (K17) Rácz G., Dobos Lszl., Beck Rbert., Szapudi Istvn., Csabai Istv., 2017, MNRAS , 469, L1 CrossRef Search ADS   Roukema B. F., 2017, preprint (arXiv:1706.06179v1) Roy X., Buchert T., 2012, Class. Quant Gravity , 29, 115004 CrossRef Search ADS   Roy X.et al.  , 2011, Class. Quant Gravity , 28, 165004 CrossRef Search ADS   Wiegand A., Buchert T., 2010, Phys. Rev. D , 82, 023523 CrossRef Search ADS   Wiltshire D. L., 2014, in Novello M., Perez Bergliaffa S.E., eds., Cosmology and Gravitation: Proc. XVth Brazilian School , Cambridge Scientific Publishers, p. 203 © 2017 The Author Published by Oxford University Press on behalf of the Royal Astronomical Society http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Monthly Notices of the Royal Astronomical Society: Letters Oxford University Press

# On Backreaction in Newtonian cosmology

, Volume 473 (1) – Jan 1, 2018
4 pages

/lp/ou_press/on-backreaction-in-newtonian-cosmology-btm054eqWf
Publisher
journal_eissn:11745-3933
ISSN
1745-3925
eISSN
1745-3933
D.O.I.
10.1093/mnrasl/slx160
Publisher site
See Article on Publisher Site

### Abstract

Abstract We clarify that a result recently stated by Kaiser is contained in a theorem of Buchert and Ehlers that is widely known for its main result: that there is no global kinematical backreaction in Newtonian cosmology. Kaiser cites this paper, re-derives parts of the theorem, but incompletely restates its content. He makes further claims, which cannot be proven beyond the limited context of Newtonian cosmology. We also discuss recent papers of Rácz et al. and Roukema, who claim the existence of global backreaction within the Newtonian framework. cosmological parameters, cosmology: theory, dark energy 1 AVERAGE PROPERTIES OF NEWTONIAN MODELS In a recent paper Kaiser (2017b, hereafter K17) considers, as do Buchert & Ehlers (1997, hereafter BE – see also Ehlers & Buchert 1997), the Euler–Poisson system in the fluid approximation for a ‘dust’ matter model in the mean field approximation of Newtonian gravity, equations (2) and (3) in K17 [cf. point (B) below]. Buchert and Ehlers performed spatial averaging of the kinematical scalars of the system to obtain the following general expansion law, given in BE, equation (B4):   \begin{eqnarray} {3{{\ddot{a}}_{\cal D} \over a_{\cal D}} + 4 {\pi} G {M_{\cal D}\over a_{\cal D}^3 } - \Lambda = :{\cal Q}_{\cal D} } \nonumber \\ {\quad=\langle {\nabla }\cdot [{\boldsymbol u} ({\nabla } \cdot {\boldsymbol u}) - ({\boldsymbol u} \cdot {\nabla }) {\boldsymbol u} ]\rangle _{\cal D} - {2\over 3} \langle {\nabla } \cdot {\boldsymbol u} \rangle _{\cal D}^2 }\nonumber \\ {\quad\quad+2\left(\Omega ^2 - \Sigma ^2 \right) + 2\left(\Omega _{ij}\langle \hat{\boldsymbol \omega }_{ij}\rangle _{\cal D} - \Sigma _{ij}\langle \hat{\boldsymbol \sigma }_{ij}\rangle _{\cal D}\right), } \end{eqnarray} (1)where $$a_{\cal D} \propto V_{\cal D}^{1/3}$$ denotes the volume scalefactor of any compact averaging domain $$\cal D$$ with volume $$V_{\cal D}$$, $$\langle \cdot \rangle _{\cal D}: = 1/ V_{\cal D} \int _{\cal D}{\mathrm{d}}V\cdot$$ the averaging operator, $$\boldsymbol u$$ the peculiar-velocity field (with its shear $$\hat{\boldsymbol \sigma }_{ij}$$ and vorticity $$\hat{\boldsymbol \omega }_{ij}$$), defined with respect to a global reference background flow, usually taken to be a Hubble flow (and Kaiser recalls correctly that this reference flow is a priori not fixed to be a Hubble flow). Omitting global shear Σ and global vorticity Ω of the model universe, we arrive at the result of Kaiser (i.e. the kinematical backreaction term, equation (9), in K17).1 Kaiser (K17, appendix A) re-derives the result of BE stating that ‘we show how this may be obtained directly’. Inspection of his derivation shows, however, that it is essentially identical to the line of reasoning of BE, but less straightforward. Kaiser obtains a linear, background-dependent term that (as he now appreciates, see footnote 1) is zero. The direct derivation of BE does not lead to such a zero term. Follow-up papers and reviews also provided complementary derivations and thorough discussions of the properties of the kinematical backreaction term (e.g. Buchert 2000a; Buchert 2008, section 3; Buchert & Räsänen 2012, section 2.5). The result (1) is valid for all compact domains from an infinitesimal domain up to any scale, e.g. up to the boundary of a compact model universe, and the second term in the second line of (1) is there on all scales. Kaiser’s paper titled Why there is no Newtonian backreaction is an improvement over the title There is no kinematic backreaction of (Kaiser 2017a), but still overstates and misleadingly interprets the content of the BE theorem for the following reasons: What BE showed, and what Kaiser re-derives, is that there is no global kinematical backreaction in Newtonian cosmology for boundary-free compact model universes.2 There is, however, kinematical backreaction, describing cosmic variance of the deviations from the assumed reference background, as Kaiser correctly states, and which is explored in the paper by Buchert et al. (2000). Thus, the suggestion that there is ‘no Newtonian backreaction’ per se is not supported by the text. Kinematical backreaction is present in Newtonian models in the interior of, e.g. a three-torus model3 (and the expectation values of the peculiar-velocity invariants are in general non-zero on regional domains). The term $${\cal Q}_{\cal D}$$ in (1) can be expressed in terms of Minkowski Functionals (Buchert 2000a, section 1.11), (Buchert 2008, section 3.1.2). These functionals describe morphological properties of the density distribution and depend on all correlation functions of a density field arising out of smoothing a point-like distribution. Consequently, even small-amplitude interior backreaction terms are important in characterizing the cosmic web. Kaiser, in K17, section 2, leaves the impression that a more elementary approach using discrete particles leads to the same result as that obtained from the Euler–Poisson system for a dust fluid in the mean field approximation. This impression comes from the comparison of the Newtonian equation of motion for particles with the Euler equation (while K17 re-derives the main result, equation 9, within a fluid mechanical approach, section 3 and appendix A). Such an approach has been discussed by Ellis & Gibbons (2014), who show4 that a set of central configuration constraints has to be also satisfied. It should be clear that the general link between a particle approach and the governing fluid equations is more subtle: Starting from discrete particles, coarse-graining the Klimontovich density in phase space and taking velocity moments leads to extra anisotropic multistream stresses modifying the Euler equation, and deviations from mean field gravity modifying the Poisson equation (see e.g. Buchert & Domínguez 2005) not done in K17 and BE. This also leads to extra backreaction terms.5 There is no ‘significant advance over the approach followed e.g. in BE’ (K17): purely space-dependent arguments already assume the vanishing of the above terms.6 Only including these terms would render the result more general.7 Kaiser omits a discussion of the uniqueness of solutions to the Poisson equation. This problem, discussed in BE, is crucial to Newtonian cosmology since we have to specify boundary conditions. BE considered various possibilities, including cases where kinematical backreaction does not globally vanish (Charlier-type models), among them the possibility of introducing a background and deviations thereof defined on a three-torus. BE also discuss why the latter is the only relevant possibility for Newtonian cosmological models (apart from Charlier-type models without a homogeneity scale). To show this, we recall that the Poisson equation does not uniquely define the potential. It is unique if we require square integrability of the potential in the case of $$\mathbb {R}^3$$. [The sum in K17, equation (5), is in general divergent, a problem of which Newton and Einstein were well aware (see e.g. Ellis & Gibbons 2014, section 1.1, and references therein).] The possibility of a square integrable potential over $$\mathbb {R}^3$$ is, however, irrelevant for cosmology. The potential is also unique if the spatial average of the source vanishes. Introducing a background and deviations thereof allows one to consider periodic peculiar-velocity fields and density contrasts with a vanishing average of the source of Poisson’s equation for a periodic peculiar-potential (a three-torus model). This architecture, being identical to that of Newtonian N-body simulations, implies that the peculiar-potential exists and is unique apart from the addition of a solution of the Laplace equation, whose only periodic solutions are spatially constant. This constant term can be removed by the translation invariance of Newton’s equations and set to zero without loss of generality. In general relativity (GR), these problems do not arise and a model universe (that is, in general, background-free) is not restricted to obey this global constraint. In addition, the kinematical backreaction variable couples to the averaged intrinsic curvature (section 3.2.2 ff in Buchert 2008), removing the ‘Newtonian anchor’ (section 3.2.3 in Buchert 2008; section 2.5 in Buchert 2000a), which enforces the vanishing of structure averages globally on an assumed background. This is the reason why the globally vanishing backreaction in Newtonian cosmology is ‘by construction’: Kaiser misinterprets this insightful remark and notes that ‘... $${\cal Q}_{\cal D}$$ tends to zero very rapidly in the limit of large volumes regardless of whether the structure is assumed to be periodic’. As we pointed out above, periodicity is a necessary element of the architecture of Newtonian models applied to cosmology; by the very definition of $${\cal Q}_{\cal D}$$, the term decays rapidly with volume, but it is crucial that it has to vanish exactly on the periodicity scale. We note the importance of a non-vanishing but small $${\cal Q}_{\cal D}$$ for the evolution of averaged scalar curvature in the GR context, see the recent paper by Bolejko (2017a), who illustrates this by employing exact GR solutions. Scale-dependent deviations from an assumed background are unavoidable.8 2 CONCLUDING REMARKS AND DISCUSSION We conclude that there is no new result in K17. Known for 20 yr, the result of BE is key to Newtonian cosmology, and is uncontroversial among researchers working on backreaction problems. However, Kaiser’s re-statement that Newtonian models cannot lead to non-vanishing backreaction is correct globally, and relevant for the rebuttal of recent work by Rácz et al. (2017) on backreaction calculated within the Newtonian framework. We therefore discuss this latter work, as well as a recent complementary work by Roukema (2017), and we point out why the conclusions drawn from these models on the possibility of global backreaction contradicts the result of BE. Since kinematical backreaction is non-vanishing in the interior of a Newtonian simulation, and attains substantial values by going to smaller but still cosmologically relevant scales (Buchert et al. 2000), it is unsurprising that Rácz et al. (2017) re-discover expansion variance by implementing a multiscale volume partitioning in N-body simulations. These toy-simulations may be interpreted to illustrate the effect on differential expansion properties, and only on a restricted range of scales. It is clear, however, and a consequence of the theorem (1), that there cannot be backreaction in these simulations globally. The key-issue to understand this fact is the non-local nature of the gravitational interaction. Rácz et al. (2017) adopt assumptions on individual domains of the global model universe. In their paper, the assumption of vanishing backreaction within a spherical idealization of domains is adopted. The non-locality of gravitation that is reflected by global constraints (boundary conditions in Newtonian cosmology and propagating spatial constraints in GR 3+1–foliated models) assures, however, that any property of the regional domain is ‘known’ to the rest of the model universe. Thus, any approximate assumption made must be rendered consistent with the whole model universe. Otherwise, an error is introduced that – as in the case of Rácz et al. (2017) – results in a fictitious global contribution of backreaction. To avoid this error, the generic situation has to be implemented, accounting in particular for backreaction on regional domains (i.e. deviations from spherical symmetry). It is of equal importance that the domains of the volume partition are be joined properly to obtain consistency with the global constraint. Furthermore, although an evolution of the background away from a homogeneous Hubble flow is a generic property of GR cosmologies on any scale,9 it lacks a physical basis within the Newtonian framework. The separate universe conjecture (Dai, Pajer & Schmidt 2015, and references therein), i.e. making restrictive assumptions on regional domains and considering the region to evolve separated from the environment, has to be put under careful scrutiny, especially in non-linear modelling attempts. Another example for an uncontrolled implementation of this conjecture is the recent paper by Roukema (2017) that suggests an alternative explanation for structure-emerging global backreaction within the Newtonian framework. While Roukema (2017) uses generic properties on subdomains, in particular including backreaction on regional domains,10 he implements the assumption of ‘silent virialization’11 by stabilizing collapsed domains with the consequence of obtaining an excess global expansion over the assumed background expansion. However, virialization in collisionless Newtonian systems can be modelled exactly by including velocity dispersion, as is explained in point (B) of Section 1. This cannot lead to a global backreaction, as explained in footnote 7. Furthermore, it does not help to appeal to the algebraic similarity between the Newtonian backreaction model of Buchert et al. (2000) and the corresponding relativistic model of Buchert, Nayet & Wiegand (2013) through the claim that ‘virialization’ induces an excess of expansion in a ‘relativistically realistic’ (Roukema 2017) situation, while the actual realization of the model is still confined to the Newtonian setting. Also in relativistic N-body simulations (that include virialization), the spatial GR constraints and their propagation has to be controlled, representing what in Newtonian theory are the boundary conditions. These simulations may also face the constraints discussed here for Newtonian models, if periodic boundary conditions on deviations from a fixed background model are implemented. A promising route to background-free investigations of the average properties of GR universe models is furnished by multiscale volume partitions that have been investigated during the last decade (a comprehensive list of works may be found in recent reviews, e.g. in Buchert & Räsänen 2012, Wiltshire 2014). In particular, the exact volume partitioning formulas for GR models given in Buchert & Carfora (2008) and Wiegand & Buchert (2010) can be profitably used to demonstrate the consequences of any regional assumption made. (The present debate motivated us to explicitly demonstrate these consequences for an exact volume partition of Newtonian models in a forthcoming paper.) It may well be possible to model regional domains of the real Universe using Newtonian equations, but this does not provide justification to model the entire observable Universe by a single Newtonian solution. Kaiser’s remarks appear to assume that small regions combine into a single cosmology in a trivial way (i.e. that the ‘background’ of all small regions should be identical). A kaleidoscopic modelling of the Universe by consecutive Newtonian domains restricts the generality substantially, for example, by calculating the average Ricci curvature of the model universe. Kaiser’s discussion of the relevance of the BE theorem for GR constitutes an opinion. More efforts to justify that the Universe is everywhere close to the same background model have been invested by Green & Wald (2014). However, their assumptions are too restrictive to apply to cosmological backreaction: see (Buchert et al. 2015) and references therein, which provide an extensive discussion of the physics of cosmological backreaction. Finally, quoting statements of early papers that pioneered the subject risks being anachronistic if developments in the literature over the subsequent 15 yr are not taken into account. (The remarks cited by Kaiser are still agreed upon today, but nuanced.) Acknowledgements This work is part of a project that, in the final stage, has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement ERC advanced grant 740021–ARTHUS, PI: T. Buchert). Thanks to Jürgen Ehlers in memoriam with whom I share the results of BE, and thanks to Nick Kaiser, Martin Kerscher, Pierre Mourier, Jan Ostrowski, Boud Roukema, Syksy Räsänen and David Wiltshire for correspondence, discussions and feedback, as well as to an anonymous referee for an insightful and very careful report. 1 The three corrections (1a), (1b) and (1c) in Buchert (2017a) of Kaiser’s re-derivation of part of BE’s results in Kaiser (2017a) have now been included (apart from the erroneous factor 3/2 in K17, equation (9), which should read 2/3) after consulting (Buchert 2017a) (N. Kaiser, private communication), except that Kaiser still claims that his derivation is more general. 2 We shall explain this in point (C). Note that kinematical backreaction also vanishes for spherical regions embedded into a reference background Buchert, Kerscher & Sicka 2000, Buchert 2011, section 7.2. This latter provides a compact proof of the Newton iron sphere theorem. Kaiser quotes instead the work of Einstein & Straus (1945), which is out of context. That backreaction vanishes for a single spherical region is well known. 3 Kaiser confuses the scale-dependent volume expansion rate $$H_{\cal D}$$ with the background expansion rate H. It is trivial that ‘There is no Newtonian backreaction on a(t) from structure’ (K17), but there is Newtonian backreaction on $$a_{\cal D}(t)$$ from structure, except on the global scale where $$a_{\cal D}(t)$$ reduces to a(t). Kaiser discusses how the equation for a(t) can be chosen to recover the usual Newtonian equations for point particles, but this does not imply that the average motion of the point particles follows the same a(t). This dismisses the non-local aspect of averaging contained in $$a_{\cal D}(t)$$. 4 Kaiser now cites this paper without pointing out the limitations explained by Ellis & Gibbons (2014). 5 In general relativity, the averaged equations have been given by including isotropic stresses (Buchert 2001). 6 Recall that in situations of multiple streams, space-dependent variables cease to be functions and we have to move to a description in phase space. 7 Keeping the mean field assumption, it is straightforward to include velocity dispersion: Taking velocity moments of the phase-space density, the Euler–Jeans equation replaces the Euler equation featuring multistream stresses Πij. The additional source, $$-\frac{1}{\varrho }\frac{{\partial} }{{\partial} x_j}\Pi _{ij} =: - \Psi _i$$, results in a divergence term, $$- \langle {\partial} \Psi _k / {\partial} x_k \rangle _{\cal D}$$, adding to the kinematical backreaction term in (1), and leaving the main result of BE on the global vanishing of the kinematical backreaction term unchanged. There is no such term in Kaiser’s derivation. 8 In GR, this holds globally and can be traced back to first principles like the non-conservation of intrinsic curvature (equation (13) in Buchert 2000b, Buchert (2008) and Buchert & Carfora 2008). 9 The quantitative importance of the backreaction effect in GR is a result of the observation that the average model evolves away from a predefined FLRW (Friedmann–Lemaître–Robertson–Walker) background due to the non-conservation of curvature (see footnote 8). A scale-dependent ‘background’ (affected by structure formation) has been investigated in the GR context (Roy & Buchert 2012), together with global instability properties of the FLRW model (Roy et al. 2011). 10 Roukema (2017) employs the analytical backreaction model of Buchert et al. (2000) that is based on exact average properties, and an analytical model for the evolution of inhomogeneities that has earlier been well-tested against N-body simulations. 11 The word ‘silent’ here refers to the separate universe conjecture rather than to so-called ‘silent model universes’ in GR, (cf. Bolejko 2017b, and references therein) that have vanishing magnetic part of the Weyl tensor. Note that even ‘silent model universes’ do not obey the separate universe conjecture, since the spatial GR constraints and their consistent propagation have to be obeyed, which again reflects the non-locality of the gravitational interaction. REFERENCES Bolejko K., 2017a, J. Cosmology Astropart. Phys. , 06, 025 CrossRef Search ADS   Bolejko K., 2017b, preprint (arXiv:1707.01800v3) Buchert T., 2000a, JGRG , 9, 306, (arXiv:gr-qc/0001056) Buchert T., 2000b, Gen. Rel. Grav. , 32, 105 CrossRef Search ADS   Buchert T., 2001, Gen. Rel. Grav. , 33, 1381 CrossRef Search ADS   Buchert T., 2008, Gen. Rel. Grav. , 40, 467 CrossRef Search ADS   Buchert T., 2011, Class. Quant Gravity , 28, 164007 CrossRef Search ADS   Buchert T., 2017a, preprint (arXiv:1704.00703v1) Buchert T., Carfora M., 2008, Class. Quant Gravity , 25, 195001 CrossRef Search ADS   Buchert T., Domínguez A., 2005, A&A , 438, 443 CrossRef Search ADS   Buchert T., Ehlers J., 1997, A&A , 320, 1 (BE) Buchert T., Räsänen S., 2012, Ann. Rev. Nucl. Part. Sci. , 62, 57 CrossRef Search ADS   Buchert T., Kerscher M., Sicka C., 2000, Phys. Rev. D , 62, 043525 CrossRef Search ADS   Buchert T., Nayet C., Wiegand A., 2013, Phys. Rev. D , 87, 123503 CrossRef Search ADS   Buchert T.et al.  , 2015, Class. Quant Gravity , 32, 215021 CrossRef Search ADS   Dai L., Pajer E., Schmidt F., 2015, J. Cosmology Astropart. Phys. , 10, 059 CrossRef Search ADS   Ehlers J., Buchert T., 1997, Gen. Rel. Grav. , 29, 733 CrossRef Search ADS   Einstein A., Straus E. G., 1945, Rev. Mod. Phys , 17, 120 CrossRef Search ADS   Ellis G. F. R., Gibbons G. W., 2014, Class. Quant Gravity , 31, 025003 CrossRef Search ADS   Green S., Wald R. M., 2014, Class. Quant Gravity , 31, 234003 CrossRef Search ADS   Kaiser N., 2017a, preprint (arXiv:1703.08809v1) Kaiser N., 2017b, MNRAS , 469, 744 (K17) Rácz G., Dobos Lszl., Beck Rbert., Szapudi Istvn., Csabai Istv., 2017, MNRAS , 469, L1 CrossRef Search ADS   Roukema B. F., 2017, preprint (arXiv:1706.06179v1) Roy X., Buchert T., 2012, Class. Quant Gravity , 29, 115004 CrossRef Search ADS   Roy X.et al.  , 2011, Class. Quant Gravity , 28, 165004 CrossRef Search ADS   Wiegand A., Buchert T., 2010, Phys. Rev. D , 82, 023523 CrossRef Search ADS   Wiltshire D. L., 2014, in Novello M., Perez Bergliaffa S.E., eds., Cosmology and Gravitation: Proc. XVth Brazilian School , Cambridge Scientific Publishers, p. 203 © 2017 The Author Published by Oxford University Press on behalf of the Royal Astronomical Society

### Journal

Monthly Notices of the Royal Astronomical Society: LettersOxford University Press

Published: Jan 1, 2018

## You’re reading a free preview. Subscribe to read the entire article.

### DeepDyve is your personal research library

It’s your single place to instantly
that matters to you.

over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month ### Explore the DeepDyve Library ### Search Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly ### Organize Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place. ### Access Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals. ### Your journals are on DeepDyve Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more. All the latest content is available, no embargo periods. DeepDyve ### Freelancer DeepDyve ### Pro Price FREE$49/month
\$360/year

Save searches from
PubMed

Create lists to

Export lists, citations