Arxiv:2012.10213V2 [Gr-Qc] 22 Mar 2021 Tensor, Called the “Law of Existence of Absolute Rotation” Candidate for the Dark Energy

The 1+3-Newton-Cartan system and Newton-Cartan cosmology

Quentin Vigneron∗ Univ Lyon, Ens de Lyon, Univ Lyon1, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, F–69007, Lyon, France (Dated: March 23, 2021) We perform a covariant 1+3 split of the Newton-Cartan equations. The resulting 3-dimensional system of equations, called the 1+3-Newton-Cartan equations, is structurally equivalent to the 1+3- Einstein equations. In particular it features the momentum constraint, and a choice of adapted coordinates corresponds to a choice of shift vector. We show that these equations reduce to the classical Newton equations without the need for special Galilean coordinates. The solutions to the 1+3-Newton-Cartan equations are equivalent to the solutions of the classical Newton equations if space is assumed to be compact or if fall-off conditions at infinity are assumed. We then show that space expansion arises as a fundamental field in Newton-Cartan theory, and not by construction as in the classical formulation of Newtonian cosmology. We recover the Buchert-Ehlers theorem for the general expansion law in Newtonian cosmology. Keywords: Newton-Cartan theory; 3+1 and 1+3 formalisms; general relativity; backreaction

I. INTRODUCTION the fluid 4-velocity, similarly to the 1+3 and 3+1 split in general relativity. The resulting system of equations The Newton-Cartan theory is a formulation of New- will be called the 1+3-Newton-Cartan equations. We will ton’s theory of gravitation in a 4-dimensional Galilei show that there is no need to choose specific coordinates, manifold. The structure on this manifold called a Galilei like Galilean coordinates, to recover the classical New- structure, follows from Newton’s ideas of absolute time ton equations. We will also complete the work of [19] as and absolute space. Similarly to general relativity, the we will study space expansion in Newton-Cartan without goal behind the Newton-Cartan formulation is to describe assuming spatial homogeneity or perturbations. the gravitational force with a spacetime connection. The In Sec.II we recall the definition and properties of physical equations constraining this connection, called Galilei spacetimes. SectionIII presents the construction the Newton-Cartan equations, are equations relating the of the 1+3-Newton-Cartan equations. We solve the sys- Riemann tensor associated to this connection and the tem of equations in Sec.IV and show that space expan- energy content. sion arises as a fundamental field of the theory. After The Newton-Cartan theory, originally introduced as a defining the gravitational field and observers in Sec.V, we spacetime geometrisation of Newton’s theory [6,7], has finally compare the solutions of the 1+3-Newton-Cartan then been developed to study the Newtonian limit (e.g. equations to the solutions of the classical Newton equa- Ref. [15]) and to define post-Newtonian approximations tions with a homogeneous deformation in Sec.VI. to general relativity (e.g. Refs. [10, 22]). Ehlers [11] also The main application we envision for the 1+3-Newton- proposed a unification of Newton-Cartan theory and gen- Cartan equations is to define a non-Euclidean Newtonian eral relativity, within his frame theory (see the overview theory, i.e. a theory which is locally equivalent to New- of literature related to the frame theory in Ref [5]). ton’s theory of gravitation but in a compact 3-manifold The Newton-Cartan equations can be written as 3- with a non-Euclidean Thurston geometry. This would be dimensional equations if adapted coordinates are cho- a powerful tool to study the effects of a non-Euclidean sen. However, to recover the classical Newton equa- global geometry on the structure formation in cosmol- tions, special Galilean coordinates are often assumed (e.g. ogy, as well as on the ‘backreaction’ of these structures [8,9, 11]), and an additional constraint on the Riemann on the global expansion of the Universe (see Ref. [2]): a arXiv:2012.10213v2 [gr-qc] 22 Mar 2021 tensor, called the “law of existence of absolute rotation” candidate for the dark energy. We discuss the possibility by Ehlers [11] needs to be added so that the Coriolis field of defining such a non-Euclidean Newtonian theory with does not depend on space. our new framework in Sec.VII. Newton-Cartan theory was also applied to cosmology in Ref. [19] where the authors introduced a fluid model with a 4-velocity similarly to general relativity. They II. GALILEI SPACETIMES projected the Newton-Cartan equations with respect to this velocity and then assumed homogeneity to derive the A. Notations expansion law in Galilei coordinates. The aim of this paper is to present the covariant 1+3- We denote a tensor of any type, except scalars, in bold split of the Newton-Cartan equations with respect to (example: g). In the case where the type is of impor- tance, a tensor of type (n, m) will feature n over-bars and m under-bars (example: g for a type (0, 2) tensor). ∗ [email protected]; [email protected] ¯ We define the symmetric part T(ab), the antisymmetric 2 part T[ab] and the symmetric traceless part Thabi of a not possible to define a spacetime metric compatible with rank-2 tensor T as the connection (4) (see Chapter 12 of Ref. [16]). Thus raising and lowering indices is not possible a priori. Then 1 1 T(ab) := (Tab + Tba); T[ab] := (Tab − Tba); when defining new tensors for the first time, we will use 2 2 the over and under bars notation introduced in Sec.IIA T Thabi := T(ab) − gab, to avoid confusion. Once they have been defined, we will D however return to the simpler bold notation. where g is the metric of the manifold on which T is de- The exact 1-form τ defines a foliation {Σt}t∈R in M, fined and D the dimension of this manifold. where Σt are spatial hypersurfaces in M defined as the An antisymmetrisation over three indices is defined as level surfaces of the scalar field t, with τ = dt. From the knowledge of τ and h, the connection ∇ 1 is not unique. Its coefficients Γγ are defined up to a T[abc] := Ta[bc] + Tc[ab] + Tb[ca] . (1) αβ 3 timelike vector B¯ and a two form κ as follows: ¯ An index which should not be included in a antisym- ¯ γ B γ µγ metrisation (over two or three indices) or a symmetrisa- Γαβ = Γαβ + 2τ(ακβ)µh , (5) tion is denoted between vertical bars. For instance, in the case of two indices antisymmetrisation: where 1 B γ γµ B 1 B γ T[a|b|c] := (Tabc − Tcba) . (2) Γ := h ∂ b − ∂ b + B ∂ τ , (6) 2 αβ (α β)µ 2 µ αβ (α β) The Lie derivative on a manifold M of a tensor T along B a vector field A¯ is denoted L T . The Lie derivative does and where b is the projector orthonormal to the vector A ¯ not commute with the metric, so for instance, for a rank- B with ¯ a 1 tensor B, g(LAB¯ , ·) 6= LAB. We will then use LAB , ¯ B µ B µβ β β respectively LABa, to denote the coordinate components bαµB := 0 ; bαµh := δα − ταB . (7) of LAB¯ , respectively LAB. Then for a vector A and¯ a tensor T on a manifold M, We have the following properties we have µ B α µ[α B β] B ∇µB = 0 ; h ∇µB = 0. (8) L T a1... := Ac∇ T a1... (3) A b1... c b1... B j where ∇ is the connection associated with the coeffi- ↓ B ¯γ X a1... c X ...c... aj cients Γ . + T ...c...∇bi A − T b1...∇cA , αβ ↑ i i j The connection coefficients (5) naturally define a Rie- σ c mann tensor R αβγ with the formula where ∇ is the Levi-Civita connection of M. The ↑ no- i tation means that c is the ith index. σ σ σ µ R αβγ := 2 ∂[βΓ + 2 Γ Γ . (9) Finally, we denote indices running from 0 to 3 by Greek γ]α µ[β γ]α letters (α, β, γ, ...) and indices running from 1 to 3 by The Ricci and Bianchi identities are still satisfied for this Roman letters (a, b, c, ...). Riemann tensor:

Rσ = 0, (10) B. Galilei structure [αβγ] σ ∇[µR |α|βγ] = 0. (11) This section is largely inspired on the presentation of Galilei spacetimes by Künzle (1972) [14]. The definitions (4) leads to the following additional rela- A Galilei spacetime is a 4-dimensional differentiable tions: manifold M with a Galilei structure (τ , h¯, ∇), where τ ¯ ¯ ¯ τ Rµ = 0, (12) is an exact 1-form, h¯ is a symmetric (2,0)-tensor of rank µ αβγ αµ µ(α β) 3, with h τµ = 0, and ∇ is a connection compatible h R µγσ = 0. (13) ¯ with τ and h¯, called a Galilei connection: µ ¯ The Ricci tensor Rαβ is defined as Rαβ := R αµβ. αβ ∇ατβ = 0 ; ∇γ h = 0. (4)

µ A vector u¯ is called timelike if u τµ = 1, and an (n,0)- α C. Coordinates in a Galilei spacetime ↓ ...µ... tensor T is called spatial if τµT = 0 for all α ∈ 1, n . No spacetime metric, i.e. a symmetric (0,2)-tensorJ ofK In this section, we deﬁne objects which will be used in rank 4, is part of the Galilei structure. Furthermore it is the construction of the 1+3-Newton-Cartan equations. 3

1. Adapted coordinate systems. 4. Spatial covariant derivative.

α A coordinate system {x }α=0,1,2,3, associated to the In an adapted coordinate system, the spatial projection ¯ γ coordinate basis vectors {∂α}α=0,1,2,3, is said to be Γab of the connection coeﬃcients (5) are adapted to the foliation {Σt}t∈R if the three vectors ¯ ¯ 1 1 {∂a}a=1,2,3 are spatial and ∂0 is timelike. Such coordi- Γγ = δγ hcd ∂ h − ∂ h , (17) nates are not unique and are determined up to the spatial ab c (a b)d 2 d ab vector freedom in the deﬁnition of the timelike vector ∂0. for any timelike vector B chosen in the relation (5). This In any adapted coordinate system B comes from the fact that bab = hab for any timelike

0 α1...αn a1...an α1 αn vector B. τα = δ ; T = T δ ...δ , (14) α a1 an The coeﬃcients (17) correspond to the coeﬃcients of the Levi-Civita connection D of the metric h on Σ. Then where T is a spatial tensor. βγ α1...αn the pull-back of h ∇γ T on Σ with T a spatial tensor gives

βγ α1...αn b a1...an 2. Pull-back. h ∇γ T → D T , (18)

α1...γ...αn and the pull-back of the divergence ∇γ T gives The relation (14) shows that any spatial tensor T is to- a1...an α1...γ...αn a1...c...an tally determined by its components T in an adapted ∇γ T → DcT . (19) coordinate system. We can then consider that T a1...an are the components of a tensor living in a 3-dimensional (hereafter 3D) Riemannian manifold Σ whose metric con- III. THE 1+3-NEWTON-CARTAN EQUATIONS travariant components are hab, thus defining a pull-back α ...α a ...a T 1 n → T 1 n . A. The Newton-Cartan equations As Σ is a Riemannian manifold, indices of tensor components on this manifold can be raised and lowered with The Newton-Cartan (NC) equations are: the metric h on Σ. We then define Ta1...an as µα ∇µT = 0, (20) c1...cn Ta ...a := T hc a ... hc a (15) µν 1 n 1 1 n n Rαβ = τατβ (4πGτµτν T − Λ) , (21) µ[α β] ab h R = 0, (22) where hab is the inverse matrix of h and corresponds to (γσ)µ the covariant components of the Riemannian metric on ¯ Σ. where Λ is the cosmological constant and T is symmetric and corresponds to the stress-energy tensor of the matter. Equation (20) is the energy and momentum conservation; Eq. (21) is the equivalent to the Einstein equation 3. Classes of adapted coordinates. and links the geometry of M to its energy content; Eq. (22) is the Trautman-Künzle condition. Given a timelike vector u, one can characterise with ¯ respect to this vector any adapted system {∂α}α=0,1,2,3 Remark. In most of the literature, the condition (22) is by introducing a vector β¯ as called the Trautman condition, citing Trautman (1963) [24].2 However Trautman originally gave the condition µ[α β] β := ∂0 − u. (16) h R γσµ = 0, i.e. without the symmetrisation. This original condition is stronger than (22). In particular The vector β is spatial and is called the shift vector of it already implies the proportionality Rαβ ∝ τατβ, i.e. ¯ Eq. (21). Based on a count of the remaining degrees of the system {∂α}α=0,1,2,3 with respect to u. The shift vector defines a class of adapted coordinate systems, de- freedom in the Riemann tensor, Künzle (1972) [14] pro- u posed instead the condition (22). It has the advantage noted Xβ . This class is the set of all adapted coordinate systems whose shift vector with respect to u is β. The of still implying the irrotationality of the gravitational systems inside a class are related by time-independent field and the closeness of κ, i.e. the original reasons for spatial changes of coordinates. the introduction of the condition by Trautman in 1963,

2 Paper written in French and available on the website of the 1 Actually ∂¯ 0 only needs not to be spatial. But by convention we Bibliothèque Nationale de France at the following web page: take it to be timelike. https://gallica.bnf.fr/ark:/12148/bpt6k4007z/f639.image. 4 but without the stronger proportionality implication. 2. 1+3-Newton-Cartan equations That is why we propose to call this condition the Trautman-Künzle condition. In this section we project the NC equations (20)-(22) with respect to τ , u, ub and h. For the remainder of this paper, we will only consider The conservation equation (20) projected along ub and that the matter is a fluid described by a timelike vector τ gives u¯ and which stress-energy tensor T is3 µ Luρ = −ρθ − ∇µq , (27) αβ α β αβ (α β) αβ T := ρu u + ph + 2q u + π , (23) u α µα µα ρ a = −h ∇µP − ∇µπ (28) α α u µ να να where ρ is the mass density, p the pressure, q¯ the heat flux − [Luq + q θ + 2 bµν q (Θ + Ω )] . and π¯ the anisotropic stress of the fluid. By definition µ µα u µν The NC equation (21) projected respectively twice q τµ = 0, π τµ = 0 and bµν π = 0. along u, along u and h, and twice along h gives:

u µ Luθ = −4πGρ + Λ + ∇µ a (29) B. 1+3 split in Newton-Cartan u u αβ µν u u αβ µν − bαµ bβν Θ Θ + bαµ bβν Ω Ω , µα αµ αµ The basics behind the 1+3 split of the NC equations h ∇µθ − ∇µ (Θ + Ω ) = 0, (30) µα νβ is to decompose the Ricci tensor along and normal to h h Rµν = 0. (31) the fluid velocity u. To do so, we first introduce the kinematical variables of the fluid in Sec. IIIB1 and then In the Trautman-Künzle condition (22), the indices α realise the split in Sec. IIIB2. and β are spatial. Then only the indices γ and σ need to This section also completes Sec. 4 in Ref. [19] by taking be split along u and h. The condition projected respec- a general fluid, and using all the NC equations. tively twice along u, along u and h, and twice along h gives:

αβ u µ[α β]ν µ[α u β] 1. The kinematical variables LuΩ = 4 bµν Θ Ω + h ∇µ a , (32) µ[α βγ] Similar to general relativity (GR), we introduce the h ∇µΩ = 0, (33) ¯ ¯ expansion tensor Θ and the vorticity tensor Ω of the hµγ hνσhζ[αRβ] = 0. (34) fluid as the projection orthogonal to the fluid of the 4- (µν)ζ velocity gradient ∇u: For the system to be closed, the relation (26) needs to be added. αβ u (β α)µ ν αβ u [β α]µ ν Θ := P ν h ∇µu ;Ω := P ν h ∇µu , (24) u β u βµ β β with P α := bµαh = δα − ταu . We denote θ := C. The equations u µν bµν Θ . u We also introduce the acceleration a¯ of the fluid 4- The 1+3-NC equations (26)-(34) are all scalar or spa- velocity u as tial equations on M. By pulling them back they become

u α µ α 3D-equations living on a 3D-Riemannian manifold Σ.A a := u ∇µu . (25) pull-back is defined for each Σt, i.e. each time t. This u implies that the geometrical properties (Riemann tensor The tensors a, Θ and Ω are all spatial. and metric) of the manifold Σ, as long as all the other We have the following additional relation tensors defined on it, are parametrised by the time. In Sec IIC4 we detailed the pull-back of the spatial αβ 1 αβ Θ = − Luh . (26) derivative, but to fully write the 1+3-NC equations as 2 3D-equations, there remains to pull-back the operator Lu present in the evolution equations. This is done by Remark. Relation (26) was originally introduced by u Ref. [23] as the definition for Θ. Note that the relation introducing a class Xβ of coordinates. Then Lu applied u u γσ 1 u on a spatial tensor T becomes, under the pull-back, “ bαγ bβσΘ = 2 Lu bαβ” given in Ref. [15] is incor- rect. α1...αn a1...an LuT → (∂t − Lβ) T , (35)

a1...an where the Lie derivative LβT applied on the spatial components of a spatial tensor corresponds to the Lie derivative on Σ. 3 For an electromagnetic ﬂuid in NC theory, see Ref. [15]. 5

The 1+3-Newton-Cartan equations (26)-(34) on Σ can in the former. It is replaced by the ﬂatness of space, i.e. then be written as: evolution equations Eq. (42). Ref. [25] showed that this results from the fact that this equation becomes a relation for the second order, in 1/c2, of the spatial curvature in the Newtonian c (∂t − Lβ) ρ = −ρθ − Dcq , (36) limit. ab ab The 1+3-NC system does not suﬀer from the lack (∂t − Lβ) h = −2Θ , (37) u c of the Hamilton constraint of GR. This is because this (∂t − Lβ) θ = −4πGρ + Λ + Dc a (38) constraint is only needed when only the 1+3-Ricci cd cd cd − Θ Θcd + Ω Ω , equation is considered, without the Raychaudhuri (∂ − L )Ω = D ua , (39) equation. As the latter equation is present, only the t β ab [a b] 1+3-Ricci equation is missing. In other words, in GR from the knowledge of the Raychaudhuri equation and and constraint equations the 1+3-Ricci equation, one can derive the Hamilton

a ac ac constraint. D θ = Dc (Θ + Ω ) , (40)

D[aΩbc] = 0, (41) Remark. The 1+3-NC system does not feature a 3Rab = 0, (42) dependence on the choice of Galilei connection. This is discussed in AppendixA. 3Rd[abc] = 0, (43) with IV. SPACE EXPANSION IN NEWTON-CARTAN u a a ca ρ a = −D P − Dcπ (44) a a ca ca − [(∂t − Lβ) q + q θ + 2qc (Θ + Ω )] . In order to derive the expansion law (in Sec.IVB), we first need solve the constraint equations (40) and (41) (in Equation (43), coming from Eq. (34), is the Ricci iden- Sec.IVA). tity for the spatial Riemann tensor, and thus is not a constraint. A. Solving the constraint equations In the system (37)-(43), the shift vector β is not physical and corresponds to a choice of coordinates. In the Sec.V we will see what choice of β leads to a Equation (41) implies that the 2-form Ω is closed, Galilean coordinate system. which translates into Ω = D w + ω , (45) Remark. Ref. [25] showed that even in the classical ab [a b] ab c formulation of Newton’s theory, the choice of time- where ω is a harmonic 2-form on Σ, i.e. DcD ωab = 0. parametrised coordinate systems is characterised by a The manifold Σ being flat with Eq. (42), we can use vector, similar to a shift vector. the decomposition theorem showed by Straumann [20] to uniquely decompose the expansion tensor into scalar, vector and tensor parts (hereafter SVT decomposition) as D. Discussions Θab = χhab + D(avb) + Ξab, (46) In the 1+3-NC system, the expansion and vorticity with χ a scalar field and Ξ is a transverse-traceless (TT) c ca tensors are not defined as the symmetric and antisym- tensor , i.e. Ξc := 0 and DcΞ := 0. The theorem is metric gradient of a spatial vector. Instead, they are de- valid for a Riemanian metric of constant scalar curvature, fined via the constraints (40) and (41). This is a major with zero traceless Ricci curvature. Fall-off conditions at difference with the classical definition of Newton’s equa- infinity or compactness of Σ also have to be added. tions (see Sec.VIA). In Sec.IVA we will see what the The link between the vectors in the expansion and vor- consequences for those two tensors are. ticity tensors is made by the momentum constraint (40), Because of these constraint equations the 1+3-NC sys- which becomes, with the decomposition (46), tem is nearly formally equivalent to the 1+3-Einstein sys- c tem of equations (see Ref. [18]). This was also spotted by Daχ = D D[cva] − D[cwa] , (47) the seminal paper on the comparison between Newton’s c using the fact that D ωca = 0 for a harmonic 2-form. theory and GR [12]. However in that paper, Θ and Ω The right-hand side (hereafter rhs) is divergence free, were not defined via the constraints (40) and (41), but whereas the left-hand side (hereafter lhs) is vorticity free. directly as the gradients of a spatial vector. Then the Hodge decomposition implies that The main difference between the 1+3-NC and 1+3- c Einstein systems is the missing of the 1+3-Ricci equation Daχ = 0 ; D D[cva] − D[cwa] = 0. (48) 6

The first equation implies that χ is only a function theorem. Equation (52) is the Friedmann equation for of time. The second equation implies that the 2-form the acceleration rate of HΣ with two additional terms cd cd D[cva] − D[cwa] is co-closed, but as it is also exact, we being ΞcdΞ Σ and ωcdω Σ. have D[cva] − D[cwa] = 0. We finally have While we have an evolution equation for χ, this is not the case for the TT tensor Ξ and the harmonic 2-form Θab = χ(t)hab +D(avb) +Ξab ;Ωab = D[cva] +ωab. (49) ω which are totally free, in their spatial and time dependence. We call Ξ the global shear, and the Coriolis We see that in general, the tensor Θ is not the ω field. The reason for this name will be given in Sec.VD. gradient of a vector, but also features non-zero scalar and tensor parts. The same applies for the vorticity tensor which features a non-zero harmonic part ω. The V. GRAVITATIONAL FIELD AND OBSERVERS physical role of these terms will be discussed in the next IN NEWTON-CARTAN THEORY section. A. The gravitational field Remark. From the NC equations (20)-(22), there are no more constraints on the harmonic 2-form ω. However, if one derives these equations from a limit of general In the literature concerning the Newton-Cartan theory, the gravitational field is often defined using the coeffi- relativity, an additional constraint appears on the γ Galilei connection (see Eq. (B4) in AppendixB), which cients Γ αβ of the Galilei connection as (e.g. Refs [8, 11]) eventually implies ω = 0. Then, the only physical choice ga := Γa . (53) on ω that is compatible with general relativity is ω = 0. 00 This is derived in AppendixB. As we did not consider This definition is however valid only in a specific adapted the NC theory as a limit of general relativity in this g coordinate system (defined by a time vector ∂t ) such paper, we will keep ω 6= 0. that

g ∂t := B, (54)

B. Space expansion in the Newton-Cartan theory where B is the vector freedom in the deﬁnition of the Galilei connection of Eq. (5). The deﬁnition (53) im-

The expansion rate HΣ(t) of Σ is defined as plies that the gravitational field is the opposite of the µ α α 4-acceleration of B, with B ∇µB = −g , and that it is 1 ∂tVΣ 1 solution of the cosmological Poisson equation {Eq. (16) HΣ(t) := = hθiΣ (t), (50) 3 VΣ 3 in Ref. [8]}. In this paper we are interested in giving a purely co- 1 R p 3 where h•i := • det(hab)d x is the spatial Σ VΣ Σ ordinate independent definition of g. We propose the average over the whole manifold Σ and VΣ := following definition: R p 3 Σ det(hab)d x is the volume of Σ. a a c a a c a u a The definition (50) of the expansion rate shows that g := (∂t − Lβ) v + 2v (Θc + Ωc ) − v Dcv − a . (55) HΣ(t) = χ(t) (51) In the following and in Sec.VC we will see that it is coherent with the standard definition (53). and thus that χ is the space expansion rate of Σ. As we explained in the previous section, the scalar χ Our definition (55) implies that the vorticity equa- enters in the SVT decomposition of the expansion tensor tion (39) becomes Θ. This tensor, along with the vorticity tensor, char- acterises the fluid. Then any part of the decomposition D[agb] = − (∂t − Lβ+v) ωab, (56) of Θ corresponds to a physical fundamental field charac- terising the fluid. In particular, Eq. (51) shows that the and the Raychaudhuri equation (38) becomes expansion, through χ, is a fundamental physical field of c 2 the 1+3-NC system (36)-(44). Dcg = − 4πGρ + Λ − 3 (∂t − Lβ) χ + χ (57) This system does not explicitly feature an evolution cd cd − ΞcdΞ + ωcdω , equation for the scalar χ. Such an equation can be obtained by taking the spatial average of the Raychaudhuri which can be rewritten, using the expansion law (52), as equation (38), which gives follows

3 (∂ − L ) H + H2 = − 4πG hρi + Λ (52) c cd cd t β Σ Σ Σ Dcg = − 4πGρb− Ξ\cdΞ + ω\cdω , (58) cd cd − ΞcdΞ Σ + ωcdω Σ . where fb := f − hfiΣ with f a scalar field. To get this equation we used the fact that divergences Equations (56) and (57) are equivalent to Eqs. (12) and averaged over a compact domain are zero due to Stokes’ (16) in Ref. [8] implying that the vector field g we defined 7 in Eq. (55) can indeed be interpreted as the gravitational this is true only if v is unique from the knowledge of field. This is further confirmed in SecVC. Θ. This is not the case as v is defined up to a spatial The definition (55) could be used as a covariant vector A whose spatial gradient DaAb is zero. For a definition of the gravitational field in GR, as it uses flat space, the solution to the equation DaAb = 0 is a the same kinematical variables Θ and Ω. However it “constant” spatial vector, i.e. corresponding to a global requires the knowledge of the vector v defined from Θ translation. Thus v is defined up to a global translation, using the SVT decomposition. This decomposition is a and therefore a Galilean observer is also defined up to a priori not possible in general in GR as the curvature global translation. orthogonal to the fluid is not necessarily zero. Thus it However once we choose the vector v this fixes the seems non-trivial to adapt the definition (55) in this Galilean observer. We denote this observer with the vec- theory. tor G, where G := u − v. Then the vector V = o − G of a general observer o corresponds to its spatial velocity Remark. Equations (56) and (58) were also obtained with respect to the Galilean observer G. Therefore, the from NC by [9, 11]. But these studies neither assume vector v = u − G is the spatial velocity of the fluid with expansion, nor global shear and therefore do not have respect to the Galilean observer. cd General and Galilean observers are represented with the term Ξ\cdΞ and the averages given by the operator respect to the fluid in Fig.1. b• in the gravitational field source equation (58). Remark. If the Galilei connection is chosen such that u B := G, then in the coordinate class X−v, i.e. ∂t := G, B. General observers we retrieve the usual definition (53) of Refs. [8, 11] for the gravitational field. An observer in the NC theory is described by a timelike vector o. A choice of coordinates can be associated to a choice of observer whose timelike vector is the time basis vector, i.e. o = ∂t. In this sense, solving the NC equations (36)-(43) in such coordinates corresponds to solving the dynamics of the fluid u with respect to the u o observer o. G For a general observer o 6= u, and we define the spatial vector as V β

V := o − u + v. (59) V Σt v We have V = β+v. The acceleration oa of this observer can then be written as Figure 1. Representation of the different vectors involved in o a a c a a = −g + (∂t − Lβ+v) V + V DcV (60) the definition of the fluid u, the Galilean observer G and a c a a a general observer o. These vectors are represented with re- + 2V (Θc + Ωc − Dcv ) , spect to a slice Σt of the foliation {Σt}t∈R. In black are the which simplifies into vectors defining the fluid; in green are the vectors defining a general observer. Note that because there is no global metric o a a c a c a a a in the structure of the embedding Galilei manifold M, the a = −g + ∂tV + V DcV + 2V (χδc + Ξc + ωc ) . (61) orthogonality on the figure between the timelike vector u and the slice Σt has no signification and is just a representation convention. C. Galilean observers

When the observer is chosen such that V = 0, then D. Test observers

oa = −ga. (62) We define a test observer with timelike vector T as a T This corresponds to an observer whose acceleration is the geodesic observer, i.e. a = 0. The equation of motion opposite of the gravitational field created by the fluid u. of these observers with respect to the Galilean observer, Such an observer is called a Galilean observer. A class of i.e. the evolution equation for V , is coordinates associated to a Galilean observer is called a a c a a c a a a ∂tV + V DcV = g − 2V (χδc + Ξc + ωc ) . (63) Galilean class. With the above definition, there seems to be a unique This corresponds to the second law of Newton with the Galilean observer, i.e. the observer o = u − v. However, velocity acceleration on the lhs and three non-inertial 8

a c a c a terms on the rhs: −2 χV , −2V Ξc and −2V ωc . The The first case is used for astrophysical flows, and the first one corresponds to an expansion force; the second second one for cosmological flows. In the second choice, c one to an anisotropic force resulting from the global the trace Hc =: 3H corresponds to the volume expan- shear; the third one is a Coriolis force created by the sion of the periodic boundary conditions imposed on P ; a Coriolis field ω. the symmetric traceless part of Hc corresponds to the c a The term −2V ωc corresponds to the Coriolis force anisotropic expansion of these conditions. Once the ho- created by a global rotation only if the Coriolis field com- mogeneous deformation is introduced, the spatial velocity ponents ωab are constants in Cartesian coordinates. As which describes the fluid is considered to be the peculiar- ω is harmonic along with compactness or fall-off condi- velocity P . 3 tions at infinity, this is the case in the present paper. The vector vH, and consequently v, is defined on R , This confirms the result of Ref. [9]. Note that Ref. [11] implying that the periodic boundary conditions imposed does not suppose closing conditions, but instead adds an on P only effectively define a compact space. Therefore additional constraint to the NC system in order for ω the homogeneous deformation is only a construction to be a global rotation: the “law of existence of absolute allowing for the description of expansion, in an effective rotation” {Eq. (23) in Ref. [11]}. However, this equa- compact space, in the classical formulation of Newton’s tion does not have a relativistic equivalent, and therefore theory. In other words, the expansion rate H quantifies cannot be obtained from a Newtonian limit. the expansion of the boundary conditions imposed by construction on the peculiar-velocity, and not the expansion of a compact space. VI. COMPARISON WITH THE CLASSICAL FORMULATION Remark. Newtonian cosmological simulations use a trace homogeneous deformation, called a Hubble flow, to “sim- We recall the classical Newton equations with a ho- ulate” expansion while using Newton’s theory. In this mogeneous deformation in Sec.VIA, and compare them case the velocity v is with the 1+3-NC equations in Sec.VIB. va = H(t)xa + P a, (68)

in Cartesian coordinates. A. Homogeneous deformation in the classical Newton theory

The classical Newton equations, in their kinematical B. 1+3-NC VS Classical Newton forms, are The expansion and vorticity tensors in the classical (∂t + Lv) ρ = −ρθ, (64) c theory with homogeneous deformation can be written as (∂t + Lv) θ = −4πGρ + Λ + Dca6=grav cd cd Θ = H(t)h + D P + H ;Ω = D P + H . − ΘcdΘ + ΩcdΩ , (65) ab ab (a b) habi ab [c a] [ab] (69) (∂ + L ) ω = D (a ) , (66) t v cd [a 6=grav b] These expressions are equivalent to those in NC [Eq. (49)], if we make the following associations: where a6=grav corresponds to the non-gravitational acceleration, and with P → v ; H → χ ; Hhabi → Ξab ; H[ab] → ωab. (70) Θab := D(avb) ;Ωab := D[avb], (67) As Hab is only a function of time in Cartesian coordi- where D is a flat connection not depending on time. nates, this implies: firstly, that Hhabi is divergence-free For the above system to be well defined, closing con- and falls into the class of TT-tensors, and secondly, that ditions need to be added. There are two possibilities: H[ab] is harmonic. This justifies the last two associations. Furthermore, as the Raychaudhuri and the vorticity (i) the system is isolated, and fall-off conditions are conservation equations are formally the same in both for- r→∞ r→∞ taken at infinity, i.e. ρ −−−→ 0 and v −−−→ 0. mulations [Eqs. (38) and (39) in NC; Eqs. (65) and (66) in classical Newton], then the equations governing the (ii) the velocity v is decomposed into a homogeneous peculiar-velocity P and the homogeneous deformation a a c deformation velocity vH with vH := Hc x (in are exactly the 1+3-NC equations in a compact space c Cartesian coordinates) where H a are functions of with the association (70). This shows that the solutions t only, and a peculiar-velocity P : we have v = to the classical Newton equations with the addition of a vH + P . The peculiar-velocity is periodically de- homogeneous deformation are equivalent to the solutions 3 R a fined on R , with the constraint P dx = 0 in of the 1+3-NC equations. Cartesian coordinates, where the integral is taken As an example, the expansion law we derived in over the volume defined by the periodic conditions. NC [Eq. (52)] is the same as the one in the classical 9 theory, given by the Buchert-Ehlers theorem {Eq. (B.6) written as a 3D-system which is formally equivalent to in Ref. [4]}. This theorem which states, in particular, the 1+3-Einstein system (apart for the 1+3-Ricci equa- that for a Hubble flow the expansion is given by the tion), and has zero Ricci curvature. Then if one wants to Friedmann equation, has been retrieved in the present keep a formal approach in the definition of a NEN theory, paper from the NC theory. and have a possible justification from GR, we think that this should be done with the 1+3-NC system. 4 Remark. From the association Hhabi → Ξab we can This was indirectly proposed by Künzle [15], who interpret Ξ as an anisotropic expansion if its compo- modified the NC equation (21) as follows nents are spatially constant in Cartesian coordinates. However, as the TT condition on Ξ is more general R(t) R − ub = τ τ (4πGτ τ T µν − Λ) , (71) (see Ref. [21] for the study of the TT tensors in flat αβ 3 αβ α β µ ν spaces) and a priori allows for non-constant components in Cartesian coordinates, we expect the physical effects with R(t) a spatial constant. This modification implies of this tensor to be more than an anisotropic expan- the same 1+3-NC equations of the present paper, with sion. If this is indeed the case, the NC theory would Eq. (42) replaced by be slightly more general than the classical Newton theory. R 3Rab = h . (72) 3 ab This corresponds to our proposed formal approach to de- VII. ON NON-EUCLIDEAN NEWTONIAN THEORIES fine a NEN theory. However solving this new system, es- pecially deriving the expansion law, is beyond the scope of this paper. One of the candidate for the dark energy, i.e. the re- This approach being only formal, it lacks of clear justi- cent acceleration of the expansion, is the effect, called fications from GR. These could be found by adapting the ‘backreaction’, of inhomogeneities on this expansion (see Newtonian limit from GR (e.g. [15]) to allow for non-zero Ref. [2]). The Buchert-Ehlers theorem we retrieved shows spatial curvature at the limit. that, in Newton’s theory, an isotropic expansion is necessarily given by the Friedman equation [Eq. (52) with Ξ = 0 and ω = 0], and therefore needs at least a cosmological constant to feature an acceleration. VIII. CONCLUSION This theorem is therefore a major result of cosmology as it tells us that if our Universe has an Euclidean ge- This paper aimed at presenting the equations re- ometry and is well described on small scales by Newton’s sulting from a covariant 1+3-split of the Newton- gravitation, dark energy could not be explained by in- Cartan equations, called the 1+3-Newton-Cartan equa- homogeneities (see Refs. [3, 13]). This result is however tions [Eqs. (36)-(44)], and the solutions to these equa- limited to Euclidean geometries, and therefore needs to tions. The main results are: be generalised for non-Euclidean ones. This requires the definition of a non-Euclidean Newtonian theory (NEN), (i) The 1+3-Newton-Cartan equations have the same i.e. a theory locally equivalent to Newton’s theory but algebraic structure as the 1+3-Einstein evolution defined in a non-Euclidean 3-manifold. In addition to and constraint equations (apart the 1+3-Ricci generalising the Buchert-Ehlers theorem, such a theory equation). In particular, as in the relativistic the- would be a powerful tool to study the effects of topology ory, a choice of spatial coordinates in Newton- on the structure formation. Cartan corresponds to a choice of shift vector. Several attempts have been made to define a NEN theory [1, 17]. Both rely on keeping formally the Poisson (ii) We give a covariant definition of the gravitational equation of the classical Newton theory and assuming field [Eq. (55) in Sec.VA], that the metric is not flat anymore, but has a curvature corresponding to the desired topology. A first issue with (iii) When solving the constraint equations, the space this approach is that it is purely formal, and is not justi- expansion arises as a fundamental physical field in fied from GR. Furthermore, it depends on the equations the theory. This contrasts with the classical theory we consider: in Newton, one can derive the Poisson equa- of Newton where the expansion is only a construction from the Raychaudhuri equation and vice versa, but tion, called homogeneous deformation. if we assume a non-flat metric, this is not the case anymore. So the gravitational field obtained with a non-flat metric, i.e. in a non-Euclidean geometry, would be different depending on the equation we kept to define the 4 Künzle claims to give references for the modification (71), but NEN theory. these are unrelated to this equation. We can therefore consider In this paper, we showed that Newton’s theory can be Ref. [15] to be the first occurrence of this modification. 10

(iv) The solutions to the 1+3-Newton-Cartan equations In this choice of connection, we see that the spatial pro- are equivalent to the solutions of the classical New- jection of κ is the vorticity of the fluid. Then this projec- ton equations with a homogeneous deformation, as- tion cannot be taken to zero as this would be a physical suming fall-off conditions at infinity or spatial com- restriction to the fluid. µα βν pactness. Only if B 6= u, one is allowed to take h κµν h = 0 without loss of generality. But in any case, the tensors B (v) We retrieve the Buchert-Ehlers theorem (Ref. [4]) and κ do not appear in the 1+3-NC equations on Σ, and in the Newton-Cartan theory, with the expansion thus their choice, in addition to having no physical impli- law (52). cations, is not relevant to the solving of these equations. (vi) We show that the Coriolis field should be zero if one Only the choice of adapted coordinates via β defining the derives the Newton-Cartan theory as a limit (e.g. partial time derivative, i.e. the choice of observer, plays Ref. [15]) of general relativity (see AppendixB). a role in Eqs. (36)-(44). We also discussed the possibility of defining a non- Euclidean Newtonian theory from the 1+3-Newton- Cartan equations. Appendix B: ω = 0 from general relativity.

1. The Newtonian limit ACKNOWLEDGEMENTS We consider a manifold M and a Lorentzian structure g This work is part of a project that has received fund- (g, ∇) on this manifold, where g is a Lorentzian metric ing from the European Research Council (ERC) under g the European Union’s Horizon 2020 research and in- and ∇ the Levi-Civita connection associated to g. novation programme (Grant agreement ERC advanced The limit allowing for the recovering of the NC equa- Grant 740021–ARTHUS, PI: Thomas Buchert). I was tions from general relativity is based on a Taylor expan- 2 supported by a ‘spécifique Normalien’ PhD Grant from sion of g in powers of λ := 1/c so that (e.g. Ref. [15]) the École Normale Supérieure de Lyon. I thank Éti- αβ αβ enne Jaupart, Pierre Mourier and Léo Brunswic for a lot g = h + O(λ), (B1) of valuable discussions on technical parts. I would also 1 g = − τ τ + O(1), (B2) like to thank Thomas Buchert and Boudewijn Roukema αβ λ α β for discussions and comments on the manuscript, Guil- laume Laibe for comments on the final manuscript, where τ is an exact 1-form and h is a (2-0)-tensor of rank 3. and Madeleine Pham-Thanh for corrections on the final g manuscript. From Eqs. (B1) and (B2), the connection ∇ can be developed in powers of λ (see the first equation on page 452 in Ref. [15]), giving Appendix A: Choice of Galilei connection g Γγ =B Γγ + τ τ hγµ∂ φ + O(λ), (B3) Throughout this paper we made no assumption on the αβ αβ α β µ choice of Galilei connection, i.e. choice of B and κ in B γ Eq. (5). However once we wrote the 1+3-NC equations where φ is an arbitrary scalar field and Γαβ is defined µ and pulled them back on Σ, the Galilei connection dis- by Eq. (6) with B a vector satisfying B τµ = 1. appears. Then the intrinsic freedom ont the definition of We see that the leading order of the Lorentzian con- this connection disappears too. As in most of the litera- nection corresponds to a Galilei connection. Then from g ture on NC, reasoning is often made with B and κ, we Eqs. (B1)-(B3), the Lorentzian structure (g, ∇) is a detail in the present section the relation between these Galilei structure (τ , h, ∇) at leading order. It is then two tensors and the kinematical variables. possible to develop the Einstein equation, the momentum We have conservation and the Bianchi identity at leading order to αβ µ(α B β) obtain the Newton-Cartan system (20)-(22). Θ = h ∇µu , (A1) αβ µ[α B β] µα βν Ω = h ∇µu + h κµν h , (A2) u α µ B α µ αν a = u ∇µu + 2u κµν h . (A3) 2. Constraint on ω If we choose B = u, then using Eq. (8) we have Contrary to the Galilei structure considered in this αβ µ(α u β) Θ = h ∇µu , (A4) paper, the one obtained from general relativity is con- αβ µα βν strained: Eq. (B3) implies the contraint Ω = h κµν h , (A5) u α µ αν a = 2u κµν h . (A6) µγ γµ 2τ(ακβ)µh = τατβh ∂µφ. (B4) 11

µα βν ab [a b] Then h κµν h = 0, which from Eq. (A2), implies As p := u − B is spatial, then Ω = D p and Ω is exact, which shows that ω = 0 from the unicity of the αβ µ[α B β] Ω = h ∇µu . (B5) Hodge decomposition. Then using Eq. (8), we have

αβ µ[α B β] Ω = h ∇µ (u − B) . (B6)

REFERENCES

[1] Barrow J. D., 2020, Non-Euclidean Newtonian cosmology, Classical and Quantum Gravity, 37, 125007 [2] Buchert T., 2008, Dark Energy from structure: a status report, General Relativity and Gravitation, 40, 467 [3] Buchert T., 2018, On Backreaction in Newtonian cosmology, MNRAS, 473, L46 [4] Buchert T., Ehlers J., 1997, Averaging inhomogeneous Newtonian cosmologies., A&A, 320,1 [5] Buchert T., Mädler T., 2019, Editorial note to: On the Newtonian limit of Einstein’s theory of gravitation (by Jürgen Ehlers), General Relativity and Gravitation, 51, 162 [6] Cartan E., 1923, Sur les variétés à connexion affine et la théorie de la relativité généralisée (première partie), Annales scientifiques de l’École Normale Supérieure, 3e série, 40, 325 [7] Cartan E., 1924, Sur les variétés à connexion affine, et la théorie de la relativité généralisée (première partie) (Suite), Annales scientifiques de l’École Normale Supérieure, 3e série, 41, 1 [8] Dautcourt G., 1990a, On the Newtonian limit of general relativity, Acta Phys. Pol. B, 21, 755 [9] Dautcourt G., 1990b, Cosmological Coriolis fields in the Newton-Cartan theory, General Relativity and Gravitation, 22, 765 [10] Dautcourt G., 1997, Post-Newtonian extension of the Newton - Cartan theory, Classical and Quantum Gravity, 14, A109 [11] Ehlers J., 2019, Republication of: On the Newtonian limit of Einstein’s theory of gravitation, General Relativity and Gravitation, 51 [12] Ellis G. F. R., 1967, Dynamics of Pressure-Free Matter in General Relativity, Journal of Mathematical Physics,8, 1171 [13] Kaiser N., 2017, Why there is no Newtonian backreaction, MNRAS, 469, 744 [14] Künzle H. P., 1972, Galilei and Lorentz structures on space-time : comparison of the corresponding geometry and physics, Annales de l’I.H.P. Physique théorique, 17, 337 [15] Künzle H. P., 1976, Covariant Newtonian limit of Lorentz space-times, General Relativity and Gravitation,7, 445 [16] Misner C. W., Thorne K. S., Wheeler J. A., 1973, Gravitation, San Francisco: W.H. Freeman and Co. [17] Roukema B. F., Różański P. T., 2009, The residual gravity acceleration effect in the Poincaré dodecahedral space, A&A, 502, 27 [18] Roy X., 2014, On the 1+3 Formalism in General Relativity, arXiv e-prints, ADS link, arXiv:1405.6319 [19] Ruede C., Straumann N., 1997, On Newton-Cartan cosmology., Helvetica Physica Acta, 70, 318 [20] Straumann N., 2008, Proof of a decomposition theorem for symmetric tensors on spaces with constant curvature, Annalen der Physik, 520, 609 [21] Tafel J., 2018, All transverse and TT tensors in flat spaces of any dimension, General Relativity and Gravitation, 50, 31 [22] Tichy W., Flanagan É. É., 2011, Covariant formulation of the post-1-Newtonian approximation to general relativity, Phys. Rev. D, 84, 044038 [23] Toupin R. A., 1957, World invariant kinematics, Archive for Rational Mechanics and Analysis,1, 181 [24] Trautman A., 1963, Sur la théorie newtonienne de la gravitation, C. R. Acad. Sci. (Paris), 257 [25] Vigneron Q., 2020, 1 +3 formulation of Newton’s equations, Phys. Rev. D, 102, 124005

When available this bibliography style features three diﬀerent links associated with three diﬀerent colors: links to the journal or editor website are in red, links to the ADS website are in blue and links to the arXiv website are in green.