Cosmological aspects of the Eisenhart–Duval lift

A cosmological extension of the Eisenhart–Duval metric is constructed by incorporating a cosmic scale factor and the energy-momentum tensor into the scheme. The dynamics of the spacetime is governed by the Ermakov–Milne–Pinney equation. Killing isometries include spatial translations and rotations, Newton–Hooke boosts and translation in the null direction. Geodesic motion in Ermakov–Milne–Pinney cosmoi is analyzed. The derivation of the Ermakov–Lewis invariant, the Friedmann equations and the Dmitriev–Zel’dovich equations within the Eisenhart–Duval framework is presented.

An important tool for solving the EMP equation is a constant of the motion known as the Ermakov-Lewis invariant -which is in fact an exact invariant for any TDHO. It is then natural to inquire about the relationship between the ED lift and the EMP equation. Can one obtain the Ermakov-Lewis invariant using the ED lift?
There is also a considerable body of work in which methods based on the EMP equation have been applied to cosmology [8][9][10][11][12][13][14][15], in particular to finding solutions and the symmetries of the Friedmann equations. This leads to the further question of what insights can be gained into cosmology using the ED lift?
A first attempt to incorporate a cosmic scale factor into the ED scheme has been reported recently [16]. Yet, because these considerations relied upon the EMP equation with a constant frequency, the resulting spacetime turned out to be stationary. It is then natural to wonder whether a dynamical spacetime with cosmological features can also be constructed along the same lines.
The purpose of this paper is to address the questions raised above, focusing on the cosmological aspects of the ED lift.
The paper is organized as follows. In Sect. 2 we provide a thorough review of the EMP equation for readers who may not be familiar with the rather extensive literature on this subject.
In Sect. 3.1 a brief account of the ED lift is given. It is shown that any Newtonian mechanical system can be represented by a metric which solves the Einstein equations provided the energy momentum tensor is chosen in a suitable way. In Sect. 3.2 the Ermakov-Lewis constant is obtained by using a conformal transformation of the associated ED metric. An explicit formula which depends upon a solution of the EMP equation is presented.
Then in Sect. 3.3 the ED lift involving a scale factor analogous to that in the Friedmann metrics is given. It is shown that the dynamics of the spacetime is governed by the EMP equation. Section 3.4 contains a detailed derivation of the conformal symmetries of this modified metric. Geodesic motion in EMP cosmoi is analyzed in Sect. 3.5.
Sect. 4.1 contains some introductory material on Friedmann-Lemaître-Robertson-Walker (FLRW) metrics, establishing the terminology and the conventions. The Einstein equations reduce to the Friedmann and Raychaudhuri equations; in the case of radiation and dark energy, they reduce to the EMP equation.
Sect. 4.2 deals with conformal symmetries and what we refer to as temporal diffeomorphisms, which may be used to provide explicit solutions to the Friedmann equations. These cases include the reduction to a time-independent harmonic oscillator.
Sect. 4.3 contains a derivation of the ED lift of geodesic motion in FLRW spacetimes. This allows describing the system's dynamical symmetries, and relating them to the null geodesics of the EMP spacetime of Sect. 3.5. The result is generalized to obtain the lift of the Dmitriev-Zel'dovich equations.
Finally in Sect. 4.4 we regard the Friedmann equations as a dynamical system with a constraint and provide a generalised ED lift of it.
In Sect. 5 we summarize our results. Throughout the paper summation over repeated indices is understood.

The Ermakov-Milne-Pinney Equation
In this section we shall give a brief account of the EMP equation [5][6][7] and of its uses. The most important of these are: • In dynamics: using its solutions to map problems involving TDHO's to problems involving time-independent harmonic oscillators (TIHO's). • In quantum mechanics: using its solutions to map problems involving the one-dimensional time-independent Schrödinger equation into Schrödinger equations which have explicit solutions.
For historical remarks about the origins of the subject the reader may consult [17]. More ambitious in scope is the Master Thesis [18]. The reader should be aware that in many papers either Milne, or Pinney is omitted from the name of the equation. The order of authors used here merely reflects the order of publication.

Comparison of two time dependent harmonic oscillators
This is perhaps the most direct approach to the EPM equation.
Apparently due originally to Dingle [19], this proceeds as follows [20]. We have two TDHO equations which may be obtained by varying respectively the two integrals We seek a temporal diffeomorphism τ = τ (t) taking one integral into the other up to a boundary term. If and where or This equation may be read in different ways.
• If the solutions Q(τ ) of equation (II.1b) are known and we want the solutions q(t) of Eq. (II.1a), we need to solve for f (t) and hence obtain τ (t) by integration • Conversely, if the solutions q(t) of (II.1a) are known and we want the solutions of Eq. (II.1b), we need to solve Eq. (II.7) for f (t).
In both cases the procedure is complicated by the need to express 2 (τ ) as a function of t. Because this involves the integral (II.8) giving τ (t), this leads in general to very complicated equations, for which the only feasible strategy appears to involve extensive iterations of the type described in [20].

The Ermakov-Milne-Pinney equation and the Ermakov-Lewis invariant
Considerable simplification results if we assume where λ is a constant which need not necessarily be positive. This leads to the original form of the Ermakov-Milne-Pinney equation: which may be thought of as a TDHO with a non-linear term. We remark here that the inverse-square potential is known to have the Schrödinger-type 2:1 scaling symmetry [21].
The EMP equation has the remarkable property [5,7] that its general solution may be expressed in terms of two linearly independent solutions, q 1 (t) and q 2 (t), of (II.1a).
where A, B, C are constants satisfying W is the constant Wronskian of the pair of linearly independent solutions q 1 , q 2 . Given a solution of the EMP Eq. (II.7) one has the identity which equates the conserved Hamiltonian (II.13a) of the TIHO (II.1b) with the conserved Ermakov-Lewis invariant (II.13b) of the TDHO (II.1a) [22,23]. This derivation is essentially the same as the one given in [24]. One may regard the Lewis invariant as an exact form of the adiabatic invariant for a TDHO, and as such, has obvious applications to the quantum theory as will be seen shortly.
Subsequent to the paper by Lewis, the idea was extended to charged particles moving in a spatially uniform magnetic field and a time dependent harmonic potential [25,26]. This is in effect covered by the theorems of Larmor and Kohn which allow one to eliminate the magnetic field by passing to a rotating frame [27][28][29]. Despite its great experimental, astronomical and technological importance, we shall not pursue this issue further.
We mention for completeness that generalizations of the Ermakov-Lewis invariant (II.13b) were obtained by considering modifications of the pair of EMP Eqs. (II.1a) and (II.10) [30,31].
Further insight is provided by the observation by Eliezer and Gray [32] that two TDHO's of the form (II.1a), q 1 (t) and q 2 (t), may be considered as a particle in the Euclidean plane E 2 , moving under the influence of an axisymmetric harmonic potential. Moreover, if one introduces polar coordinates r, θ such that q 1 + iq 2 = r e iθ , the equations of motion of the Lagrangian imply the conservation of the angular momentum (per unit mass) Then substitution back to the radial equation of motion gives r + ω 2 (t)r = h 2 r 3 (II. 16) which, in the case λ = h 2 > 0, coincides with the EMP Eq. (II.10) if r = f . If one substitutes q = r cos θ, f = r, λ = h 2 into (II.13b), one finds that the Ermakov-Lewis quantity is in fact the square of the angular momentum 1 (II.17) A list of cases for which there are explicit solutions of the radial Eq. (II.16) are given in [32,34]. Because of its application to cosmology which we will describe later, we reproduce it here.
• ω = constant = ω 0 , The orbit in the plane is an ellipse.
which yields a logarithmic spiral. • where J n , Y n are Bessel functions. If k > 0 , then r → 0 as t → ∞, and if k < 0 then r → ∞ as t → ∞. 1 The case λ < 0 is obtained by replacing the Euclidean plane E 2 by the Minkowski plane E 1,1 [33].

The one-dimensional Schrödinger equation and the Liouville-Green-Jeffreys-Wentzel-Kramers-Brillouin approximation
If we now think of the temporal coordinates t, τ as spatial coordinates x, X , the complexified positions q(t), Q(τ ) as wave functions ψ(x), (X ) and the frequencies ω(x) and (τ ) as wave vectors k(x) and K (x), we may transcribe all of the fore-going theory into quantum mechanical language. The TDHO Eqs. (II.1a)-(II.1b) become  [35], which contains a systematic procedure for obtaining asymptotic expansions, provides a list of references and justifies the addition of the names of Liouville, Green and Jeffreys to the better known Wentzel, Brillouin, and Kramers.
If we choose K 2 = X 2 we find formal bound state solutions of (II.22a) which take the form where the dependence of r (x) on E has been suppressed. where we have re-instated the dependence of r on E as well as x. 3 The quantum mechanical treatment presented above was based on the formal equivalence of the TDHO Eq. (II.1a) and the time-independent Schrödinger Eq. (II.22a) effected by exchanging the time coordinate for the spatial coordinate of the latter. A deeper connection may be obtained by considering a wave packet solution of the time-dependent Schrödinger equation in one spatial dimension [37][38][39][40], is the real valued expectation value of the position of the "classical" trajectory and p = mq is the classical momentum; the functions N (t) and K (t) not relevant for the following discussion. In [37] it is claimed that substitution of (II.26) into (II.25) leads to the TDHO Eq. (II.1a) and the complex Riccati-type equation then f satisfies the EMP Eq. (II.10); the Ermakov-Lewis quantity I given by (II.13b) is a dynamical invariant. We remark en passant that the EMP equation may be incorporated into the formalism of supersymmetric quantum mechanics [36,37].

The Ermakov-Milne-Pinney equation and Madelung's hydrodynamic transcription
The time-dependent Schrödinger equation with a general potential V (x), can be brought into the hydrodynamical form : Setting = √ ρ e i S h , defining the velocity v = d x dt = 1 m ∇S and taking real and imaginary parts, one gets equations which resemble Euler's equations for a compressible irrotational fluid, called a Madelung fluid [41] ∂ρ ∂t For a stationary state one has ∂ S ∂t = −E, ∂ρ ∂t = 0. One may also require that the probability density be independent of time [37]. In general one typically assumes that ∇S = 0 and hence √ ρ satisfies the linear equation On the other hand in one spatial dimension, taking note of (III.17), one may alternatively assume for some constant C. This leads to the non-linear EMP Eq. (II.10) for f = √ ρ(x) [37], 3 with Similar applications of the EMP equation in this context, including those with magnetic fields, may be found in [42,43].

Bose-Einstein condensates
Among the many applications of the EMP equation is one to Bose-Einstein condensates (BEC's) [44,45]. The authors start with the spherically symmetric but timedependent Gross-Pitaevski equation in d spatial dimensions with radial coordinate r , which may be derived from the Hamiltonian Then, using standard "moment methods" in the literature, an equation is obtained for a quantity y(t) [44,45] with K a constant. The case d = 2 is evidently special and in that case, setting f = √ y(t) and ω 2 (t) = 2λ(t), one obtains (II.7), with (II.38)

Eisenhart-Duval lift of Newtonian mechanics
As originally formulated in [1][2][3], the ED lift provides a geometric description of a mechanical system with d degrees of freedom x i , i = 1, . . . , d, and the potential energy U (t, x) in terms of geodesics of the Lorentzian metric on a (d + 2)dimensional spacetime Computing the Christoffel symbols and analyzing the geodesic equations one finds that t is affinely related to the proper time τ where κ is a constant, x i obeys Newton's Eq. (passing from τ to t) while the dynamics of v is fixed from the condition that the geodesic is null or time-like where = 0 is for null geodesics and = −1 for timelike geodesics. Newtonian mechanics is thus recovered by implementing the null reduction along v [1]. Note that the original construction [1][2][3] dealt with null geodesics only. The spacetime (III.1) belongs to the Kundt class as it admits the covariantly constant null Killing vector field The latter can be used to construct the trace-free energymomentum tensor in a geometrically rather appealing way 4 where (y) 2 is an arbitrary function (the energy density). The only non-vanishing component of ξ μ is ξ t = − 1 2 , which gives T tt = d 8π 2 while the rest vanishes. Because 4 The factor d 2π , d being the dimension of the x-subspace, is chosen for further convenience. ξ μ ∂ μ = 0 holds true if does not depend on v, the energymomentum tensor is conserved, The Ricci tensor which derives from (III.1) has only one nonzero component while the scalar curvature vanishes, Given the ED metric (III.1) and the energy-momentum tensor (III.6), the Einstein equations imply that the contribution of the cosmological term necessarily vanishes, = 0, thus reducing (III.9) to Only the (tt)-component is non-trivial; it links to U , as A particularly interesting example of the ED geometry occurs if U and are t-independent. Setting and interpreting G as Newton's constant and ρ(x) as the mass density, one recovers the Newton equation for the gravitational potential.
In conclusion, Newtonian mechanics can be represented in terms of a metric which solves the Einstein Eqs. (III.10) provided the energy-momentum tensor is chosen in the form (III.6).

Lifting the Ermakov-Lewis invariant
The aim of this subsection is to show how applying the ED lift to a TDHO allows one to obtain the Ermakov-Lewis invariant by performing a conformal transformation of the metric [1][2][3].
The equations of motion of a TDHO follow from the Hamiltonian Let H = −p t and let The ED metric on the extended ("Bargmann" [2,3]) spacetime M = {R 3 , g} is given by, The null geodesics, considered as unparametrised curves, of a Lorentzian metric g μν are the same as for any conformally related Lorentzian metric. They are given by the Hamiltonian flow on T M with Hamiltonian 16) subject to the constraint Hamilton's equations following from the Hamiltonian (III. 16) with constraint (III.17) for two conformally related metrics g μν and 2 (x)g μν are identical. The vector field V = ∂ ∂v is a null Killing vector field and thus To obtain the motion in 2 spacetime dimensions we perform a Marsden-Weinstein reduction (referred to as "ignoration of cyclic coordinates" in old fashioned books) bearing in mind the constraint (III.17) and setting p v = 1 to recover the Hamiltonian system (III.13). We now set for some function f (t) to be determined. In (Q, t, v) coordinates the metric may be cast in the form where A is the one-form and g(t) is an arbitrary function of time. If A is closed, d A = 0, then we may define a new coordinateṽ by in terms of which the metric becomes Now introducing a new time coordinate τ by we have The condition for A to be closed is where an arbitrary constant, then the metric becomes As noted above, two conformally related metrics have the same null geodesics, although their affine parameters differ. From this we deduce that the null geodesics of the static metric inside the braces in (III.28) are the same as those for our time dependent metric (III.15).
If we now perform a null reduction on the Killing vector ∂ ∂ṽ , we obtain the Hamiltonian of a TIHÕ

Ermakov-Milne-Pinney cosmology
Now we show that the EMP equation naturally arises if one incorporates a cosmic scale factor into a suitably chosen Bargmann metric, when the energy-momentum tensor is chosen in a proper way. Consider a (d + 2)-dimensional spacetime parametrized by the coordinates y μ = (t, v, x i ), i = 1, . . . , d, and endowed with the Lorentzian metric [16] where a(t) is an arbitrary function and γ is a constant. For a fixed value of t the line element in the d-dimensional slice parametrized by x i is given by a(t) 2 dx i dx i . Therefore a(t) 2 may be interpreted as a cosmic scale factor. The metric (III.30) possesses, as does its conventional counterpart (III.1), a covariantly constant null Killing vector field, namely (III.5). The choice of (III.30) will be justified a posteriori in Sect. 4. Constructing the energy-momentum tensor following the prescription (III.6) and specifying to the class of energy densities which depend on t only, = (t), from Eq. (III.10) one finds Thus the dynamics of (III.30) is governed by the EMP equation. The instance of = const has been discussed recently in [16], in which case the EMP equation reduces to conformal mechanics in one dimension [21].
One can learn more about the geometry of (III.30) by analyzing the geodesic equations. Computing the Christoffel symbols, one concludes that t is affinely related to the proper time, (III.32) The coordinate x i obeys in turn the oscillator-like equation in which we passed from τ to t, while the evolution of v is fixed from the condition that the geodesic be null or timelike. Equation (III.33) prompts one to introduce the conformal time which brings the metric (III.30) to the form where a(η) = a(t (η)), and t (η) is the inverse of η(t) in (III.34). Equation (III.35) is an analog of the flat (k=0) FLRW cosmological model in which the Minkowski metric has been changed into the simplest PP-wave; the Friedmann equation which determines the evolution of the cosmic scale factor is replaced, in this case, by the EMP equation.
To conclude this section we note that the coordinate transformation brings (III.30) to the form where we omitted the primes, or, in view of (III.31), 38) which is the d + 2 dimensional Bargmann metric associated with an isotropic oscillator in (d, 1) dimensions with time-dependent frequency (t), the motions of which are the projections of the null geodesics of (III.38).

Symmetries as conformal Killing isometries
Finding the symmetries of time dependent harmonic oscillators generated an extensive literature in the early 1980's, including [46][47][48][49]. A discussion in terms of canonical transformations is in [51]. An alternative approach is presented below in terms of the ED lift. For simplicity, we stick to Eq. (III.38). Following [2,3], the symmetries of a non-relativistic system in (d, 1) dimensions can be obtained as a subgroup of the conformal symmetries of the d + 2 dimensional Bargmann manifold: one selects those conformal transformations that leave invariant the covariantly constant null vector ∂ v .
In our case, the symmetries given in Ref. [49] are seen to be consistent with the Schrödinger group in d dimensions -which is a subgroup of the conformal group of the extended (Bargmann) spacetime. Now we re-derive the above-mentioned symmetries in the specific case of our EMP spacetime of Sect.3.3.
Consider indeed a generic infinitesimal transformation one gets a coupled set of partial differential equations to fix δt, δv, δx i and . Omitting details, we present the (conformal) isometries of (III.38), where , ν, κ, ω i j = −ω ji are constant infinitesimal parameters, x 2 = x i x i , the functions μ i (t), ρ i (t), and λ(t) obey the ordinary differential equations while the conformal factor is Taking into account the order of the differential Eq. (III.42), one finds that the transformations (III.41) involve 6 + 4d + independent infinitesimal parametersthe same as that of the flat-space conformal group in d + 2 dimensions.
From Eq. (III.43) one concludes that the κ,λ(t), , and ρ i (t)-transformations give rise to conformal Killing vectors, 5 while the ν, μ i (t), ω i j and constant λ-transformations generate Killing vectors. In view of (III.42), the isometry with constant parameter λ is only possible for constant frequency which corresponds to a stationary spacetime. Time translations are broken in general.
We notice that, while distinct conformal isometries act on the coordinates as in (III.41), the functions λ(t), ρ i (t) and μ i (t) are not independent. Obviously, ρ i (t) differs from μ i (t) by an inessential constant. If μ satisfies the first equation in (III.42), then λ(t) = μ 2 (t) will satisfy the second. Moreover, if μ 1 and μ 2 are independent solutions of the first equation, then the three independent solutions of the second are given by μ 2 1 , μ 2 2 , μ 1 μ 2 . The deeper reason why this happens is that the conformal Killing vectors form a Lie algebra, which we identify below with so(2, 2 + d).
The interpretation of the (conformal) Killing vectors above and their algebra becomes more transparent if one switches to conformal time, (III.35). It is evident that all (conformal) isometries of the PP-wave metric will be automatically transmitted to become (conformal) isometries of (III.35), because the two expressions only differ by the cosmic scale factor a(η). In particular, all the symmetry transformations of (III.44), which do not involve η explicitly, will be transformed into the Killing vectors of (III.35), while those affecting η will be transmitted into conformal Killing vectors. Our clue is that the metric above is conformally flat. This follows from the vanishing of the Weyl tensor. Skipping details we merely mention that it can also be seen, explicitly, by applying the Arnold transformation, see [4,38,39].
The group of conformal transformations of any (1, 1+d)dimensional conformally flat spacetime is isomorphic to that of Minkowski space, explaining the "coincidence" we noted earlier.
From the first equation in (III.42) one finds, where the infinitesimal parameters α i and β i are associated with spatial translations and Newton-Hooke boosts [52]. The second equations generate where the infinitesimal parameters θ , σ , ρ, linked to time translations, special conformal transformations, and dilatations form an so(2, 1) subalgebra. Along with spatial rotations described by ω i j , the μ i (η) and λ(η)-transformations form the conformal Newton-Hooke algebra. For a detailed discussion of the Schrödinger and conformal Newton-Hooke algebras and their realizations in spacetime see e.g. [53][54][55].
As it follows from (III.41), the isometries of (III.44) also involve the translation in the v-direction, while the set of conformal isometries contains the κ, and -transformations.
Note that the κ-transformation in (III.41) and (III.43) is realized in a way analogous to conventional dilatation in the Schrödinger algebra. It derives from the latter by replacing the temporal variable t by the "null" coordinate v; it has appeared before in the context of gravitational waves [56]. The -transformation is in turn an analog of special conformal transformation, again t replaced by v. It is straightforward to verify that, along with the translations in the vdirection, v = v+ν, they form an so(2, 1) subalgebra. Interestingly enough and extending the Galilei-Carroll "duality" [57], the latter acts upon the null coordinate v in very much the same way as so(2, 1) entering the conformal Newton-Hooke algebra affects the temporal coordinate t.
To summarize, the algebra of vector fields which involve both Killing and conformal Killing vectors can be identi-fied with so(2, 2+d), the conformal Newton-Hooke algebra being its subalgebra.
Having identified the conformal isometries of the metric (III.38), the symmetries of the underlying classical system in one fewer dimension (i.e., the time-dependent oscillator) could now be derived. Skipping details, we just mention that implementing the null reduction along v, the SO(2, 1) conformal subgroup with parameters ν, κ and in (III.41) will be broken, allowing us to recover the Schrödinger symmetry found in [49]. The generators are conveniently identified using the formulae in Sect. 3 of [50].

Geodesic motion in Ermakov-Milne-Pinney cosmoi
Having established a link between the EMP equation and the ED lift, let us study the geodesic motion in EMP cosmoi. The analysis is facilitated by switching to the conformal time (III.35) which allows one to solve the geodesic equations by quadrature, where τ 0 , κ, α i , β i , v 0 are constants of integration and α 2 = α i α i , αβ = α i β i . It is assumed that = 0 for null geodesics and = −1 for time-like geodesics. Before proceeding, it is worth mentioning that, because for null geodesics Eq. (III.48a) admit a particular solution the metric (III.35) is formulated in a reference frame comoving with a light signal which travels along the v-axis. This correlates with the fact that (III.35) differs from a PP-wave by a scale factor only. The first line in (III.48a) defines η as an implicit function of the proper time τ . Although in most cases of interest the integral η η 0 a 2 (η)dη cannot be evaluated exactly, Eq. (III.48a) prove to be sufficient to comprehend a qualitative behaviour of geodesics. Indeed, Eq. (III.48b) defines an ellipse. By making use of rotational invariance, one can set the ellipse to lie, say, in the x 1 x 2 -plane. As v evolves with time, geodesics in the EMP cosmology wrap around the elliptic cylinder, v being its axis.
For null geodesics ( = 0), the trajectory is a closed loop and the motion is periodic (see Fig. 1). For time-like geodesics ( = −1), the orbit wraps around the cylinder remaining in a compact region of space for some time, then   (x 1 , x 2 , v) parametrizing the two horizontal and one vertical axes, respectively the -term in Eq. (III.48c) starts dominating and the particle escapes (see Fig. 2). Given the coordinate system in which the metric (III.35) is formulated, this happens because massive particles travel slower than a reference frame comoving with a light signal which propagates along the v-axis.
Little is known about analytic solutions to the EMP equation. Assuming that the energy density 2 (t) decreases with time (which happens in an expanding universe), one has (see, where n is a rational number such that ν is positive, h is a positive constant, and J ν , Y ν are Bessel functions. Note that in these cases a(t) is a monotonically increasing convex function which tends to a fixed nonzero value as t → 0.
Other interesting examples are provided by negative integer ν; then ν should be replaced by |ν| in the expression for a(t) in Eq. (III.50b). These models are represented by monotonically increasing convex functions starting at a(0) = 0. The instance t = 0 can be interpreted as the Big Bang. A typical example is shown in Fig. 3.
A notable simplification takes place forä = 0, when the solution is expressed in terms of elementary functions This case corresponds to a linearly expanding universe in which the energy density of matter decreases consistently with the inverse square law. The general solution to the geodesic equations is, where t 0 , κ, α i , β i , v 0 are constants of integration. 6 From now on we switch back to cosmic time t.
Although the graphs representing the cosmic scale factors in the EMP cosmology look quite reasonable, the geodesic motion is apparently unrealistic. This happens because Eq. (III.48b) defines an ellipse. An obvious cure is to generalize the construction to the case of time-dependent γ which will alter the qualitative behaviour of geodesics. It proves sufficient to replace γ in (III.30) by an arbitrary function γ (t) which will then show up on the right hand side of Eq. (III.31), viewed as an algebraic equation to fix γ (t) in terms of a(t) and (t). In this way one can model a reasonable geodesic behavior in the generalized EMP cosmoi by properly choosing the cosmic scale factor and the energy density. Note, however, that, as the associated geodesic equations involve a time-dependent oscillator, with frequency γ (t) 2 , finding an analytic solution may be complicated.

Friedmann-Lemaître-Robertson-Walker spacetimes
There has been considerable interest in the past few years in applying the ideas circling around the EMP equation to problems in cosmology, both classical and quantum, [8][9][10][11][12][13][14][15]. In this section we shall establish some connections. Before doing so we shall begin by establishing our notations and conventions.
A Friedmann-Lemaître-Robertson-Walker (FLRW) spacetime is one carrying the metric where d 2 K is a 3-dimensional metric of constant curvature K . The scale factor a(t) is governed by three equations of which, ifȧ = 0, only two are independent. These are the Raychaudhuri equation the Friedmann equation and the conservation equatioṅ where ρ is the energy density, P is the pressure and is the cosmological constant.

Matter models
The simplest models for the matter content are those in which it consists of a perfect fluid with constant barotropic index. This means that where γ , or equivalently w, is constant. It follows from (IV.4) that One has the following special cases • γ = 1 corresponds to pressure-free matter.
• γ = 4 3 to radiation. • γ = 0 to a cosmological constant which is equivalent to a constant energy density ρ = 8π G . • γ = 2 is the largest value of γ consistent with the dominant energy condition. If the energy density is positive it is sometimes known as "stiff matter". If γ > 2 there exist solutions for which the scale factor blows up in finite time. • γ = 2 3 corresponds to a gas of cosmic strings. • γ = 1 3 corresponds to a gas of membranes. • If γ < 2 3 , there exist solutions exhibiting a "big rip", that is for which the scale factor blows up in finite time [58] .
Another simple model is to suppose that the matter consists of a scalar field φ with potential V (φ) It follows from (IV.4) thaẗ 2 this behaves like stiff matter, γ = 2. If |V (φ)| >> 1 2φ 2 this behaves like a cosmological constant.

Cosmic clocks and temporal diffeomorphisms
The coordinate t is called cosmic time and may be measured by atomic clocks. To good approximation it also coincides with local astronomical time. For some purposes it is convenient to introduce another coordinate η called conformal time, defined by η = dt a(t) , (IV.9) in terms of which the FLRW metric is manifestly conformally static If χ is radial distance, conformal time may be measured by light clocks. In other words, the ticks of a light clock correspond to light rays being bounced between two co-moving mirrors. This procedure is often attributed to Marzke and Wheeler [59], but in fact it goes back to Einstein, Robb, Milne, Whitrow and Walker. A survey of explicit solutions of the Friedmann equations for matter with a constant barotropic index both as a function of cosmic time and of conformal time may be found in [60]. Concerns about the notion of a Big Bang and the beginning of the universe have led to consideration of other cosmic clocks or choices of what have variously been called time scales or time graduations [61][62][63]. One obvious choice is the temperature or equivalently the redshift where T 0 ≈ 3K and the subscript 0 denotes the present epoch. Now if a(t) runs from 0 at the Big Bang to ∞ in the infinite future then T runs from ∞ to 0. In fact given three successive instants of time t 1 ≺ t 2 ≺ t 3 the temperature changes, T 12 = T (t 1 ) T (t 2 ) etc behave multiplicatively Of course merely introducing a new time coordinate will not turn a geodesically incomplete spacetime into a geodesically complete one, but it does underscore the point that rejecting a geodesically incomplete spacetime as being singular requires a physical justification, not merely a convenient mathematical definition in order to prove theorems. If K > 0 this is a simple harmonic oscillator. If K < 0 this is an "inverted" simple harmonic oscillator. Both may be construed as a free particle moving in a 2-dimensional Newton-Hooke spacetime with coordinates (τ, x) [52]. The systems are invariant under the conformal extensions of Newton-Hooke-groups N H ± [52,55].
A related construction works if K = 0 but we include a cosmological constant term or dark energy term.
If > 0 multiplication of (IV.3) by a 3γ leads to an upside down simple harmonic oscillator equation for a One may calculate the jerk [64] j = a 2 d 3 a dt 3 da dt . (IV.20) If γ = 1, j = 1 for all times. Otherwise it merely starts from a value greater than one and monotonically decreases to 1 at late times. For example for radiation, γ = 4 3 , the jerk starts from 3.
If < 0 the motion is converted to a standard harmonic oscillator. For radiation γ = 4 3 and (IV.2) becomes which is a relationship between dynamical variables of order 2 in the momenta. Then an important step is to transform this relation into a homogeneous one by introducing a conserved momentum, p v -which is the physical equivalent of constructing a projective conic from a standard one, as described in [68]. In our case the process transforms (IV.27) into This dynamical relation can be recovered by studying the null geodesics of the Hamiltonian which is associated with the metric After a translation, v =ṽ − m η a(η )dη , we obtain In particular, if the spatial slices are flat then the metric is conformally flat. Its conformal symmetries contain the Schrödinger group as a subgroup, providing us with symmetries also of the original FLRW Hamiltonian H F L RW in (IV.25). Another consequence is that the null geodesics of the ED lift just constructed are the same as the null geodesics of the EMP spacetime, described in Sect.3.3.
We conclude this subsection by generalizing our findings to N particles in a FLRW spacetime, subject to their mutual gravitational interaction. The Hamiltonian we use is where the p A i are the spatial momenta of the individual particles and r A their coordinates.
Our Hamiltonian relies on the assumption that the proper times of all particles can be synchronized and described by the same parameter. This is possible because the spacetime is spatially homogeneous. As before, we set the Hamiltonian equal to a constant value −M 2 /2, calculate the lowenergy limit and make the relationship homogeneous in the momenta. We obtain the null Hamiltonian where again we have redefined the v variable in order to remove the term Ma(η) p 2 v . Going back to cosmic time, with p η = ap t , we obtain which belongs to the same conformal class as the ED lift of the Hamiltonian presented in [70]. This gives rise to the Dmitriev-Zel'dovich equations [52,69,70]. If we consider, as is done in almost all numerical studies of large scale structure, the motion of N nonuniformities about a FLRW background metric, they are governed by the Dmitriev-Zeldovich equations [69,70] d dt where a, b = 1, 2, . . . , N label the non-uniformities and m a are their masses. It is simple to verify that if one replaces every x a in (IV.35) by (IV.37) This system Eisenhart-Duval lifts to Up to signs and factors this coincides with the metric (III.30) as long as γ = 0 and G = 0. In other words a free nonrelativistic particle in a background Friedmann universe lifts to a null geodesic of (III.30) provided γ = 0. In fact the timedependent harmonic oscillator in (III.30) looks very much like a cosmological constant term.

The Eisenhart-Duval lift of the Friedmann equations
The action density of Einstein gravity coupled to a scalar field is To understand better what is happening we note that in setting g tt = −1 in the Friedmann metric we have used up some gauge freedom and as a consequence the L does not depend upon sufficiently many independent variables in order to obtain the full set of equations of motion. To rectify this problem, on considers the more general metric (IV. 46) We see that the lapse function N enters (IV.46) as a Lagrange multiplier and varying the LagrangianL with respect to N and then setting N = 1 one gets the Friedmann equation, i.e. the Hamiltonian constraint (IV.44). In fact in terms of the canonical momenta we have 7 We are now in the situation to perform an ED lift or Marsden-Weinstein oxidation, we set p t = −H (IV. 48) and introduce an ignorable momentum p v and rewrite (IV.47) as Note that if we set N = 1 and compute the 4-metric then the coefficients of da 2 and dφ 2 have the opposite sign. If we regard N as a further coordinate, then we have a Galilean or Newton-Cartan metric.

Conclusion
Let us recapitulate the results obtained in this work: • Any Newtonian mechanical system can be described in terms of the Eisenhart-Duval metric which solves the Einstein equations (III.10). The key ingredient involved in the construction is the energy-momentum tensor (III.6) built out of the covariantly constant null Killing vector field (III.5) and a proper energy density function.