Cosmological aspects of the Eisenhart-Duval lift

M. Cariglia¹*, A. Galajinsky²†, G.W. Gibbons³‡, P.A. Horvathy⁴§

¹DEFIS, Universidade Federal de Ouro Preto, MG-Brasil

²Tomsk Polytechnic University, Russia

³D.A.M.T.P., Cambridge University, U.K.

⁴Laboratoire de Mathématiques et de Physique Théorique,
Université de Tours, France

(Dated: April 11, 2018)

Abstract

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 cosmology 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.

**Keywords:** Cosmology, Eisenhart-Duval lift, Ermakov-Milne-Pinney equation, time-dependent oscillator, Newton-Hooke group.

<footnote>* e-mail: marco@iceb.ufop.br</footnote>

<footnote>† e-mail: galajin@tpu.ru</footnote>

<footnote>‡ e-mail: G.W.Gibbons@damtp.cam.ac.uk</footnote>

<footnote>§ e-mail: horvathy@lmpt.univ-tours.fr</footnote>
# Contents

I. Introduction 3

II. The Ermakov-Milne-Pinney Equation 4

A. Comparison of two time-dependent harmonic oscillators 4

B. The Ermakov-Milne-Pinney equation and the Ermakov-Lewis invariant 5

C. The one-dimensional Schrödinger equation and the Liouville-Green-Jeffreys-Wentzel-Kramers-Brillouin approximation 7

D. The Ermakov-Milne-Pinney equation and Madelung's hydrodynamic transcription 9

E. Bose-Einstein condensates 10

III. The Eisenhart-Duval Lift and the Ermakov-Milne-Pinney Equation 10

A. Eisenhart-Duval lift of Newtonian mechanics 10

B. Lifting the Ermakov-Lewis invariant 12

C. Ermakov-Milne-Pinney cosmology 14

D. Symmetries as conformal Killing isometries 15

E. Geodesic motion in Ermakov-Milne-Pinney cosmoi 18

IV. Friedmann-Lemaître-Robertson-Walker Spacetimes 21

A. Matter models 21

B. Cosmic clocks and temporal diffeomorphisms 22

C. Lifting the Friedmann-Lemaître-Robertson-Walker geodesics 25

D. The Eisenhart-Duval lift of the Friedmann equations 27

V. Conclusion 28

Acknowledgments 29

References 29
# I. INTRODUCTION

The Eisenhart-Duval lift [1, 2] provides a convenient framework for treating time dependent dynamical systems and their symmetries [3]. Such systems include time-dependent harmonic oscillators. time-dependent harmonic oscillators and their symmetries have a large literature. In particular, they are conveniently studied using the Ermakov-Milne-Pinney equation [4–6].

An important tool for solving the Ermakov-Milne-Pinney equation is a constant of the motion known as the *Ermakov-Lewis invariant* – which is in fact an exact invariant for any time-dependent harmonic oscillator. It is then natural to inquire about the relationship between the Eisenhart-Duval lift and the Ermakov-Milne-Pinney equation. Can one obtain the Ermakov-Lewis invariant using the Eisenhart-Duval lift?

There is also a considerable body of work in which methods based on the Ermakov-Milne-Pinney equation have been applied to cosmology [7–14], 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 Eisenhart-Duval lift?

A first attempt to incorporate a cosmic scale factor into the Eisenhart-Duval scheme has been reported recently [15]. Yet, since these considerations relied upon the Ermakov-Milne-Pinney 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 Eisenhart-Duval lift.

The paper is organized as follows. In section II we provide a thorough review of the Ermakov-Milne-Pinney equation for readers who may not be familiar with the rather extensive literature on this subject.

In subsection III A a brief account of the Eisenhart-Duval 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 subsection III B the Ermakov-Lewis constant is obtained by using a conformal transformation of the associated Eisenhart-Duval metric. An explicit formula which depends upon a solution of the Ermakov-Milne-Pinney equation is presented.

Then in subsection III C the Eisenhart-Duval 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 Ermakov-Milne-Pinney equation. Subsection IIID contains a detailed derivation of the conformal symmetries of this modified metric. Geodesic motion in Ermakov-Milne-Pinney cosmoi is analyzed in subsection IIIE.

Subsection IVA 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 Ermakov-Milne-Pinney equation.

Subsection IVB 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.

Subsection IVC contains a derivation of the Eisenhart-Duval lift of geodesic motion in FLRW spacetimes. This allows describing the system's dynamical symmetries, and relating them to the null geodesics of the Ermakov-Milne-Pinney spacetime of subsection III E. The result is generalized to obtain the lift of the Dmitriev-Zel'dovich equations.

Finally in subsection IVD we regard the Friedmann equations as a dynamical system with a constraint and provide a generalised Eisenhart-Duval lift of it.

In section V we summarize our results.

Throughout the paper summation over repeated indices is understood.

## II. THE ERMAKOV-MILNE-PINNEY EQUATION

In this section we shall give a brief account of the Ermakov-Milne-Pinney equation [4–6] and of its uses. The most important of these are:

* In dynamics: using its solutions to map problems involving time-dependent harmonic oscillators to problems involving time-independent harmonic oscillators.

* 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 [16]. More ambitious in scope is the Master Thesis [17]. 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.

### A. Comparison of two time-dependent harmonic oscillators

This is perhaps the most direct approach to the EPM equation. Apparently due originally to Dingle [18], this proceeds as follows [19]. We have two time-dependent harmonic oscillator equations

$$ \frac{d^2 q}{dt^2} + \omega^2(t) q = 0, \qquad (II.1a) $$

$$ \frac{d^2 Q}{d\tau^2} + \Omega^2(\tau) Q = 0. \qquad (II.1b) $$

which may be obtained by varying respectively the two integrals

$$ \int \frac{1}{2} \left( \left( \frac{dq}{dt} \right)^2 - \omega^2(t) q^2 \right) dt, \qquad (II.2a) $$

$$ \int \frac{1}{2} \left( \left( \frac{dQ}{d\tau} \right)^2 - \Omega^2(\tau) Q^2 \right) d\tau. \qquad (II.2b) $$
We seek a temporal diffeomorphism $\tau = \tau(t)$ taking one integral into the other up to a boundary term. If

$$f(t) = \left(\frac{d\tau}{dt}\right)^{-\frac{1}{2}}, \qquad (II.3)$$

and

$$q(t) = \left(\frac{d\tau}{dt}\right)^{-\frac{1}{2}} Q(\tau) = f(t)Q(\tau) \qquad (II.4)$$

we find that

$$\int \left( \left( \frac{dQ}{d\tau} \right)^2 - \Omega^2(\tau) Q^2 \right) d\tau = \int \frac{1}{2} \left( \left( \frac{dq}{dt} \right)^2 - \omega^2(t) q^2 \right) dt - \int \frac{d}{dt} \left( q^2 \frac{1}{f^2} \frac{df}{dt} \right) dt, \qquad (II.5)$$

where

$$\omega^2(t) = \frac{\Omega^2(\tau)}{f^4} - \frac{1}{f} \frac{d^2 f}{dt^2} \qquad (II.6)$$

or

$$\frac{d^2 f}{dt^2} + \omega^2 f = \frac{\Omega^2(\tau)}{f^3}. \qquad (II.7)$$

This equation may be read in different ways.

* If the solutions $Q(\tau)$ of equation (II.1b) are known and we want the solutions $q(t)$ of equation (II.1a), we need to solve for $f(t)$ and hence obtain $\tau(t)$ by integration

$$\tau = \int \frac{dt}{f^2}. \qquad (II.8)$$

* Conversely, if the solutions $q(t)$ of (II.1a) are known and we want the solutions of equation (II.1b) we need to solve equation (II.7) for $f(t)$.

In both cases the procedure is complicated by the need to express $\Omega^2(\tau)$ as a function of $t$. Since this involves the integral (II.8) giving $\tau(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 [19].

## B. The Ermakov-Milne-Pinney equation and the Ermakov-Lewis invariant

Considerable simplification results if we assume

$$\Omega^2(\tau) = \lambda, \qquad (II.9)$$

where $\lambda$ is a constant which need not necessarily be positive. This leads to the original form of the Ermakov-Milne-Pinney equation:

$$\frac{d^2 f}{dt^2} + \omega^2(t) f = \frac{\lambda}{f^3}, \qquad (II.10)$$

which may be thought of as a time-dependent harmonic oscillator 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 [20].
The Ermakov-Milne-Pinney equation has the remarkable property [4, 6] that its general solution may be expressed in terms of two linearly independent solutions, $q_1(t)$ and $q_2(t)$, of (II.1a).

$$f(t) = (Aq_1^2(t) + Bq_2^2(t) + 2Cq_1(t)q_2(t))^{\frac{1}{2}} \qquad (II.11)$$

where A, B, C are constants satisfying

$$AB - C^2 = \frac{\lambda}{W^2}, \qquad W = q_1 \frac{dq_2}{dt} - q_2 \frac{dq_1}{dt}. \qquad (II.12)$$

W is the constant Wronskian of the pair of linearly independent solutions $q_1, q_2$.

Given a solution of the Ermakov-Milne-Pinney equation (II.7) one has the identity

$$H = \frac{1}{2} \left( \left( \frac{dQ}{d\tau} \right)^2 - \lambda Q^2 \right) \qquad (II.13a)$$

$$= \frac{1}{2} \left( \frac{q^2}{f^2} + \left( \frac{dq}{dt} - \frac{df}{dt}q \right)^2 \right) = I, \qquad (II.13b)$$

which equates the conserved Hamiltonian (II.13a) of the time-independent harmonic oscillator (II.1b) with the conserved Ermakov-Lewis invariant (II.13b) of the time-dependent harmonic oscillator (II.1a) [21, 22]. This derivation is essentially the same as the one given in [23].

One may regard the Lewis invariant as an exact form of the adiabatic invariant for a time-dependent harmonic oscillator, 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 [24, 25]. 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 [26, 27]. 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 Ermakov-Milne-Pinney equations (II.1a) and (II.10) [28, 29].

Further insight is provided by the observation by Eliezer and Gray [30] that two time-dependent harmonic oscillators of the form (II.1a), $q_i(t)$ and $q_2(t)$, may be considered as a particle in the Euclidean plane $\mathbb{E}^2$, moving under the influence of an axisymmetric harmonic potential. Moreover, if one introduces polar coordinates $r, \theta$ such that $q_1 + iq_2 = re^{i\theta}$, the equations of motion of the Lagrangian

$$L = \frac{1}{2}(\dot{r}^2 + r^2\dot{\theta}^2 - \omega^2(t)r^2), \qquad (II.14)$$

imply the conservation of the angular momentum (per unit mass)

$$r^2 \dot{\theta} = h. \qquad (II.15)$$

Then substitution back to the radial equation of motion gives

$$\ddot{r} + \omega^2(t)r = \frac{h^2}{r^3} \qquad (II.16)$$
which, in the case $\lambda = h^2 > 0$, coincides with the Ermakov-Milne-Pinney equation (II.10) if $r = f$.

If one substitutes $q = r \cos \theta$, $f = r$, $\lambda = h^2$ into (II.13b), one finds that the Ermakov-Lewis quantity is in fact the square of the angular momentum<fnref n="1" />

$$I = h^2. \qquad (II.17)$$

A list of cases for which there are explicit solutions of the radial equation (II.16) are given in [30, 32]. Because of its application to cosmology which we will describe later, we reproduce it here.

$$\bullet \qquad \omega = \text{constant} = \omega_0, \qquad r^2 = (A \cos \omega_0 t + B \sin \omega_0 t)^2 + \frac{h^2}{A^2 \omega_0^2} \sin^2 \omega_0 t. \qquad (II.18)$$

The orbit in the plane is an ellipse.

$$\bullet \qquad \omega = \frac{b}{t^2}, \qquad \theta - \theta_0 = -\frac{b}{t}, \qquad r = \sqrt{\frac{h}{b}} t = \frac{\sqrt{hb}}{\theta_0 - \theta}. \qquad (II.19)$$

$$\bullet \qquad \omega = \frac{b}{t}, \qquad r^2 = bt, \qquad b = \frac{h}{\sqrt{b^2 - \frac{1}{4}}}, \qquad b \neq \frac{1}{2}, \qquad \theta = \theta_0 + \frac{h}{b} \ln t, \qquad r = b \exp \frac{b}{h}(\theta - \theta_0), \qquad (II.20)$$

which yields a logarithmic spiral.

$$\bullet \qquad \omega = bt^k, \qquad r^2 = \pi h A t \left( J_n^2(2Ab^2t^{k+1}) + Y_n^2(2Ab^2t^{k+1}) \right), \qquad n = \frac{1}{2}(k+1), \qquad (II.21)$$

where $J_n, Y_n$ are Bessel functions. If $k > 0$, then $r \to 0$ as $t \to \infty$, and if $k < 0$ then $r \to \infty$ as $t \to \infty$.

### C. The one-dimensional Schrödinger equation and the Liouville-Green-Jeffreys-Wentzel-Kramers-Brillouin approximation

If we now think of the temporal coordinates $t, \tau$ as spatial coordinates $x, X$, the complexified positions $q(t), Q(\tau)$ as wave functions $\psi(x), \Psi(X)$ and the frequencies $\omega(x)$ and $\Omega(\tau)$ as wave vectors $k(x)$ and $K(x)$, we may transcribe all of the foregoing theory

<footnote n="1">¹ The case $\lambda < 0$ is obtained by replacing the Euclidean plane $\mathbb{E}^2$ by the Minkowski plane $\mathbb{E}^{1,1}$ [31].</footnote>
into quantum mechanical language. The time-dependent harmonic oscillator equations (II.1a)-(II.1b) become

$$ \frac{d^2\psi(x)}{dx^2} + k^2(x)\psi = 0, \qquad (II.22a) $$

$$ \frac{d^2\Psi(x)}{dX^2} + K^2(X)\Psi = 0, \qquad (II.22b) $$

where $k^2 = 2(E - V(x))$ and we have used units in which $\hbar^2 = m$, $E$ is the energy, and $V(x)$ the potential energy.

The relation to the Liouville-Green-Jeffreys-Wentzel-Kramers-Brillouin method may be seen by noting that (II.3) and (II.4), (II.8), and (II.7), are identical to (2), (6), (3) in section 17.122 of [33], 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

$$ \psi(x) \propto r(x) \sin \left( \int_{x_o}^{x} \frac{d\tilde{x}}{r^2(\tilde{x})} \right), \qquad (II.23) $$

where the dependence of $r(x)$ on $E$ has been suppressed. Eqn (II.23) implies the *Milne quantization condition*<fnref n="2" /> for the allowed values of the energy $E$ [5],

$$ \frac{1}{\pi} \int_{-\infty}^{\infty} \frac{dx}{r^2(x, E)} = n + 1, \quad n = 0, 1, 2, 3, \dots, \qquad (II.24) $$

where we have re-instated the dependence of $r$ on $E$ as well as $x$<fnref n="3" />.

The quantum mechanical treatment presented above was based on the formal equivalence of the time-dependent harmonic oscillator equation (II.1a) and the time-independent Schrödinger equation (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 [35–38],

$$ i\hbar \frac{\partial \Psi(x, t)}{\partial t} = \left[ -\frac{\hbar^2}{2m} \frac{\partial^2}{\partial x^2} + \frac{1}{2} m \omega^2(t) x^2 \right] \Psi(x, t) \qquad (II.25) $$

for

$$ \Psi(x, t) = N(t) \exp \left[ i \left( y(t) \tilde{x}^2 + K(t) \right) \right], \qquad (II.26) $$

where $\tilde{x} = x - q(t)$, and

$$ q(t) \langle x \rangle = \int_{-\infty}^{\infty} \bar{\Psi} x \Psi dx \qquad (II.27) $$

is the real valued expectation value of the position of the “classical” trajectory and $\langle p \rangle = m\dot{q}$ is the classical momentum; the functions $N(t)$ and $K(t)$ not relevant for the following

<footnote n="2">² This entitles Milne to his place in the name of (II.10) [34].</footnote>

<footnote n="3">³ Note that $\hbar^2/2m$ in (3) of [34] should be replaced by its inverse.</footnote>
discussion. In [35] it is claimed that substitution of (II.26) into (II.25) leads to the time-
dependent harmonic oscillator equation (II.1a) and the complex Riccati-type equation

$$
\frac{2\hbar}{m}\dot{y} + \left(\frac{2\hbar}{m}y\right)^2 + \omega^2 t = 0. \qquad (II.28)
$$

Moreover, if

$$
\mathrm{Re} y = \frac{m}{2\hbar} = \frac{1}{f^2}, \qquad (\mathrm{II}.29)
$$

then *f* satisfies the Ermakov-Milne-Pinney equation (II.10); the Ermakov-Lewis quantity *I* given by (II.13b) is a dynamical invariant.

We remark *en passant* that the Ermakov-Milne-Pinney equation may be incorporated
into the formalism of supersymmetric quantum mechanics [34, 35].

D. The Ermakov-Milne-Pinney equation and Madelung’s hydrodynamic tran-
scription

The time-dependent Schrödinger equation with a general potential $V(\boldsymbol{x})$,

$$
i\hbar \frac{\partial \Psi}{\partial t} = -\frac{\hbar^2}{2m} \nabla^2 \Psi + V\Psi, \qquad (II.30)
$$

can be brought into the hydrodynamical form : Setting $\Psi = \sqrt{\rho} e^{i \frac{S}{\hbar}}$, defining the velocity
$v = \frac{dx}{dt} = \frac{1}{m} \nabla S$ and taking real and imaginary parts, one gets equations which resemble
Euler's equations for a compressible irrotational fluid, called a *Madelung fluid* [39]

$$
\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \qquad (II.31a)
$$

$$
\frac{\partial \boldsymbol{v}}{\partial t} + \boldsymbol{v} \cdot \nabla \boldsymbol{v} = -\frac{1}{m} \nabla \tilde{V}. \qquad (II.31b)
$$

For a stationary state one has $\frac{\partial S}{\partial t} = -E$, $\frac{\partial \rho}{\partial t} = 0$.

One may also require that the probability density be independent of time [35]. In
general one typically assumes that $\nabla S = 0$ and hence $\sqrt{\rho}$ satisfies the linear equation

$$
-\frac{\hbar^2}{2m} \nabla^2 \sqrt{\rho} + V \sqrt{\rho} = E \sqrt{\rho}. \qquad (II.32)
$$

On the other hand in one spatial dimension, taking note of (III.17), one may alternatively
assume

$$
\frac{dS}{dx} = \frac{C}{\rho} \qquad (II.33)
$$

for some constant *C*. This leads to the *non-linear* Ermakov-Milne-Pinney equation (II.10)
for $f = \sqrt{\rho(x)}$ [35],

$$
\frac{d^2 \sqrt{\rho}}{dx^2} + \omega^2(x) \sqrt{\rho} = \frac{C^2}{\hbar^2} \frac{1}{(\sqrt{\rho})^3} \quad \text{with} \quad \omega^2(x) = \frac{2m}{\hbar} (E - V(x)). \qquad (II.34)
$$

Similar applications of the Ermakov-Milne-Pinney equation in this context, including
those with magnetic fields, may be found in [40, 41].
## E. Bose-Einstein condensates

Among the many applications of the Ermakov-Milne-Pinney equation is one to Bose-Einstein condensates (BEC's) [42].

The authors start with the spherically symmetric but time-dependent Gross-Pitaevski equation in $d$ spatial dimensions with radial coordinate $r$,

$$i \frac{\partial \Psi(r, t)}{\partial t} = -\frac{1}{2} \nabla^2 \Psi(r, t) + \left( \lambda(t) r^2 + \nu(t) |\Psi(r, t)|^2 \right) \Psi(r, t), \quad (II.35)$$

which may be derived from the Hamiltonian

$$H = \frac{1}{2} \int \left( |\nabla \Psi|^2 + \nu(t) |\Psi(r, t)|^4 + 2\lambda(t) r^2 |\Psi(r, t)|^2 \right) r^{d-1} dr. \quad (II.36)$$

Then, using standard “moment methods” in the literature, an equation is obtained for a quantity $y(t)$ [42]

$$\frac{d}{dt} \left( y \ddot{y} - \frac{1}{2} \dot{y}^2 - 2K\nu(t) + 4\lambda(t)y^2 \right) = 2(d-2)Ky \frac{d}{dt} \left( \frac{\nu(t)}{y} \right) \quad (II.37)$$

with $K$ a constant. The case $d=2$ is evidently special and in that case, setting $f = \sqrt{y(t)}$ and $\omega^2(t) = 2\lambda(t)$, one obtains (II.7), with

$$\Omega^2(\tau) = \frac{1}{2}C + K\nu(t). \quad (II.38)$$

# III. THE EISENHART-DUVAL LIFT AND THE ERMAKOV-MILNE-PINNEY EQUATION

## A. Eisenhart-Duval lift of Newtonian mechanics

As originally formulated in [1, 2], the Eisenhart-Duval lift provides a geometric description of a mechanical system with $d$ degrees of freedom $x_i$, $i = 1, \dots, d$, and the potential energy $U(t, x)$ in terms of geodesics of the Lorentzian metric on a $(d+2)$-dimensional spacetime

$$g_{\mu\nu}(y)dy^\mu dy^\nu = -2U(t, x)dt^2 - dtdv + dx_i dx_i, \quad (III.1)$$

where $y^\mu = (t, v, x_i)$. Computing the Christoffel symbols and analyzing the geodesic equations one finds that $t$ is affinely related to the proper time $\tau$

$$\frac{d^2 t}{d\tau^2} = 0 \quad \Rightarrow \quad \frac{dt}{d\tau} = \kappa, \quad (III.2)$$

where $\kappa$ is a constant, $x_i$ obeys Newton's equation (passing from $\tau$ to $t$)

$$\frac{d^2 x_i}{dt^2} + \partial_i U(t, x) = 0, \quad (III.3)$$
while the dynamics of v is fixed from the condition that the geodesic is null or time-like

$$ \frac{dv}{dt} = \frac{dx_i}{dt} \frac{dx_i}{dt} - 2U - \frac{\epsilon}{\kappa^2}, \qquad (III.4) $$

where $\epsilon = 0$ is for null geodesics and $\epsilon = -1$ for time-like geodesics. Newtonian mechanics is thus recovered by implementing the null reduction along v [1]. Note that the original construction [1, 2] dealt with null geodesics only.

The spacetime (III.1) belongs to the Kundt class as it admits the covariantly constant null Killing vector field

$$ \xi^{\mu} \partial_{\mu} = \partial_{v}, \qquad \nabla_{\mu} \xi_{\nu} = 0, \qquad \xi^{2} = 0. \qquad (III.5) $$

The latter can be used to construct the trace-free energy-momentum tensor in a geometrically rather appealing way<fnref n="4" />

$$ T_{\mu\nu} = \frac{d}{2\pi}\Omega(y)^2\xi_{\mu}\xi_{\nu}, \qquad T^{\mu}_{\mu} = 0, \qquad \nabla_{\rho}T_{\mu\nu} = \frac{d}{\pi}\Omega\partial_{\rho}\Omega\xi_{\mu}\xi_{\nu}, \qquad (III.6) $$

where $\Omega(y)^2$ is an arbitrary function (the energy density). The only non-vanishing component of $\xi_\mu$ is $\xi_t = -\frac{1}{2}$, which gives $T_{tt} = \frac{d}{8\pi}\Omega^2$ while the rest vanishes. Since $\xi^\mu \partial_\mu \Omega = 0$ holds true if $\Omega$ does not depend on v, the energy-momentum tensor is conserved,

$$ \nabla^{\mu} T_{\mu\nu} = 0, \qquad (III.7) $$

provided $\Omega = \Omega(t, x)$.

The Ricci tensor which derives from (III.1) has only one nonzero component while the scalar curvature vanishes,

$$ R_{tt} = \partial_i \partial_i U, \qquad R = 0. \qquad (III.8) $$

Given the Eisenhart-Duval metric (III.1) and the energy-momentum tensor (III.6), the Einstein equations

$$ R_{\mu\nu} - \frac{1}{2}g_{\mu\nu}(R + 2\Lambda) = 8\pi T_{\mu\nu} \qquad (III.9) $$

imply that the contribution of the cosmological term necessarily vanishes, $\Lambda = 0$, thus reducing (III.9) to

$$ R_{\mu\nu} = 8\pi T_{\mu\nu}. \qquad (III.10) $$

Only the (tt)-component is non-trivial; it links $\Omega$ to $U$, as

$$ \Omega^2 = \frac{1}{d} \partial_i \partial_i U. \qquad (III.11) $$

A particularly interesting example of the Eisenhart-Duval geometry occurs if $U$ and $\Omega$ are $t$-independent. Setting

$$ \Omega^2(x) = \frac{4\pi G}{d} \rho(x) \qquad (III.12) $$

<footnote n="4">⁴ The factor $\frac{d}{2\pi}$, $d$ being the dimension of the $x$-subspace, is chosen for further convenience.</footnote>
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 equations (III.10) provided the energy–momentum tensor is chosen in the form (III.6).

## B. Lifting the Ermakov-Lewis invariant

The aim of this subsection is to show how applying the Eisenhart-Duval lift to a time-dependent harmonic oscillator allows one to obtain the Ermakov-Lewis invariant by performing a conformal transformation of the metric [1, 2].

The equations of motion of a time-dependent harmonic oscillator follow from the Hamiltonian

$$H = \frac{1}{2} (m^{-1}(t)p_q^2 + m(t)\omega^2(t)q^2). \qquad (III.13)$$

Let $H = -p_t$ and let

$$\mathcal{H}(x^{\mu}, p_{\mu}) = 2p_{t}p_{v} + m^{-1}p_{q}^{2} + m\omega^{2}q^{2}p_{v}^{2} = g^{\mu\nu}p_{\mu}p_{\nu} \qquad (III.14)$$

with $x^{\mu} = (t, v, q)$, $p_{\mu} = (p_t, p_v, p_q)$. The Eisenhart-Duval metric on the extended (“Bargmann” [2]) spacetime $M = \{\mathbb{R}^3, g\}$ is given by,

$$g_{\mu\nu}dx^{\mu}dx^{\nu} = 2dtdv + mdq^2 - m\omega^2q^2dt^2. \qquad (III.15)$$

The null geodesics, considered as unparametrised curves, of a Lorentzian metric $g_{\mu\nu}$ are the same as for any conformally related Lorentzian metric. They are given by the Hamiltonian flow on $T^*M$ with Hamiltonian

$$\mathcal{H} = g^{\mu\nu}p_{\mu}p_{\nu} \qquad (III.16)$$

subject to the constraint

$$\mathcal{H} = 0. \qquad (III.17)$$

Hamilton’s equations following from the Hamiltonian (III.16) with constraint (III.17) for two conformally related metrics $g_{\mu\nu}$ and $\Omega^2(x)g_{\mu\nu}$ are identical.

The vector field $V = \frac{\partial}{\partial v}$ is a null Killing vector field and thus

$$p_v = \text{constant}. \qquad (III.18)$$

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

$$Q = \frac{q}{f} \qquad (III.19)$$
for some function $f(t)$ to be determined. In $(Q, t, v)$ coordinates the metric may be cast in the form

$$ds^2 = 2dt(v + A) + mf^2dQ^2 - Q^2g(t)dt^2, \qquad (III.20)$$

where $A$ is the one-form

$$A = mf\dot{f}QdQ + \frac{1}{2}mQ^2(\dot{f}^2 - f^2\omega^2)dt + \frac{1}{2}Q^2g(t)dt \qquad (III.21)$$

and $g(t)$ is an arbitrary function of time. If $A$ is closed, $dA = 0$, then we may define a new coordinate $\tilde{v}$ by

$$d\tilde{v} = dv + A, \qquad (III.22)$$

in terms of which the metric becomes

$$ds^2 = 2dtd\tilde{v} + mf^2dQ^2 - Q^2g^2dt^2. \qquad (III.23)$$

Now introducing a new time coordinate $\tau$ by

$$dt = mf^2d\tau, \qquad (III.24)$$

we have

$$ds^2 = mf^2\{2d\tau d\tilde{v} + dQ^2 - mf^2Q^2gd\tau^2\}. \qquad (III.25)$$

The condition for $A$ to be closed is

$$\frac{d}{dt}(mf\dot{f}) + m\omega^2 f = \frac{g}{f}. \qquad (III.26)$$

If

$$g = \frac{\Omega^2}{mf^2}, \qquad (III.27)$$

where $\Omega$ an arbitrary constant, then the metric becomes

$$ds^2 = mf^2\{2d\tau d\tilde{v} + dQ^2 - Q^2d\tau^2\}. \qquad (III.28)$$

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 $\frac{\partial}{\partial\tilde{v}}$, we obtain the Hamiltonian of a time-independent harmonic oscillator

$$\tilde{H} = \frac{1}{2}(p_Q^2 + \Omega^2 Q^2). \qquad (III.29)$$
### C. Ermakov-Milne-Pinney cosmology

Now we show that the Ermakov-Milne-Pinney 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^\mu = (t, v, x_i)$, $i = 1, \dots, d$, and endowed with the Lorentzian metric [15]

$$ds^2 = -\frac{\gamma^2 x_i x_i}{a(t)^2} dt^2 - dtdv + a(t)^2 dx_i dx_i, \qquad (III.30)$$

where $a(t)$ is an arbitrary function and $\gamma$ 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. IV.

Constructing the energy-momentum tensor following the prescription (III.6) and specifying to the class of energy densities which depend on $t$ only, $\Omega = \Omega(t)$, from Eq. (III.10) one finds

$$\ddot{a} + \Omega^2(t)a = \frac{\gamma^2}{a^3}. \qquad (III.31)$$

Thus the dynamics of (III.30) is governed by the Ermakov-Milne-Pinney equation. The instance of $\Omega = \text{const}$ has been discussed recently in [15], in which case the Ermakov-Milne-Pinney equation reduces to conformal mechanics in one dimension [20].

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,

$$\frac{d^2 t}{d\tau^2} = 0. \qquad (III.32)$$

The coordinate $x_i$ obeys in turn the oscillator-like equation

$$a^2 \frac{d}{dt} \left( a^2 \frac{d}{dt} x_i \right) + \gamma^2 x_i = 0, \qquad (III.33)$$

in which we passed from $\tau$ to $t$, while the evolution of $v$ is fixed from the condition that the geodesic be null or time-like. Eq. (III.33) prompts one to introduce the conformal time

$$a^2(t) \frac{d}{dt} = \frac{d}{d\eta}, \qquad \eta(t) = \int_{t_0}^{t} \frac{d\tilde{t}}{a^2(\tilde{t})}, \qquad (III.34)$$

which brings the metric (III.30) to the form

$$ds^2 = a^2(\eta) \left( -\gamma^2 x_i x_i d\eta^2 - d\eta dv + dx_i dx_i \right), \qquad (III.35)$$

where $a(\eta) = a(t(\eta))$, and $t(\eta)$ is the inverse of $\eta(t)$ in (III.34). Eq. (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 Ermakov-Milne-Pinney equation.
To conclude this section we note that the coordinate transformation

$$t' = t, \quad x'_i = a(t)x_i, \quad v' = v + a(t)\dot{a}(t)x_i x_i \quad (III.36)$$

brings (III.30) to the form

$$ds^2 = \left(\frac{\ddot{a}}{a} - \frac{\gamma^2}{a^4}\right) x_i x_i dt^2 - dtdv + dx_i dx_i, \quad (III.37)$$

where we omitted the primes, or, in view of (III.31),

$$ds^2 = -\Omega(t)^2 x_i x_i dt^2 - dtdv + dx_i dx_i \quad (III.38)$$

which is the $d+2$ dimensional Bargmann metric associated with an isotropic oscillator in $(d, 1)$ dimensions with time-dependent frequency $\Omega(t)$, the motions of which are the projections of the null geodesics of (III.38).

## D. Symmetries as conformal Killing isometries

Finding the symmetries of time-dependent harmonic oscillators generated an extensive literature in the early 1980's, including [43–46]. A discussion in terms of canonical transformations is in [48]. An alternative approach is presented below in terms of the Eisenhart-Duval lift. For simplicity, we stick to Eq. (III.38).

Following [2], 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 $\partial_v$.

In our case, the symmetries given in Ref. [46] 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 Ermakov-Milne-Pinney spacetime of sec.III C.

Consider indeed a generic infinitesimal transformation

$$t' = t + \delta t, \quad v' = v + \delta v, \quad x'_i = x_i + \delta x_i, \quad (III.39)$$

where $\delta t$, $\delta v$, $\delta x_i$ are arbitrary functions of $t, v, x$. Demanding (III.38) to be invariant under (III.39) up to a conformal factor,

$$ds'^2 = (1 + \Lambda)ds^2, \quad (III.40)$$

one gets a coupled set of partial differential equations to fix $\delta t$, $\delta v$, $\delta x_i$ and $\Lambda$. Omitting
details, we present the (conformal) isometries of (III.38),

$$t' = t + \lambda(t) + \frac{1}{2}\epsilon x^2 - x_i \rho_i(t), \qquad (\text{III.41a})$$

$$x'_i = x_i + \mu_i(t) + \omega_{ij}x_j + \frac{1}{2}\kappa x_i + \frac{1}{2}\dot{\lambda}(t)x_i + \frac{1}{2}v\epsilon x_i \qquad (\text{III.41b})$$

$$+ \frac{1}{2}x^2 \dot{\rho}_i(t) - x_i x_j \dot{\rho}_j(t) - \frac{1}{2}v \rho_i(t), \qquad (\text{III.41c})$$

$$v' = v + \nu + 2\dot{\mu}_i(t)x_i + \frac{1}{2}\ddot{\lambda}(t)x^2 - \frac{1}{2}\epsilon\Omega^2 x^2 x^2 + \kappa v + \frac{1}{2}\epsilon v^2 \qquad (\text{III.41d})$$

$$+ \Omega^2 x^2 x_i \rho_i(t) - v x_i \dot{\rho}_i(t), \qquad (\text{III.41e})$$

where $\epsilon, \nu, \kappa, \omega_{ij} = -\omega_{ji}$ are constant infinitesimal parameters, $x^2 = x_i x_i$, the functions $\mu_i(t), \rho_i(t)$, and $\lambda(t)$ obey the ordinary differential equations

$$\ddot{\mu}_i + \Omega^2 \mu_i = 0, \qquad \ddot{\rho}_i + \Omega^2 \rho_i = 0, \qquad \lambda^{(3)} + 4\Omega^2 \dot{\lambda} + 4\Omega \dot{\Omega} \lambda = 0, \qquad (\text{III.42})$$

while the conformal factor is

$$\Lambda = \kappa + \dot{\lambda}(t) + v\epsilon - 2x_i \dot{\rho}_i(t). \qquad (\text{III.43})$$

Taking into account the order of the differential equations (III.42), one finds that the transformations (III.41) involve $6 + 4d + \frac{d(d-1)}{2} = \frac{(4+d)(3+d)}{2}$ independent infinitesimal parameters — the same as that of the flat-space conformal group in $d+2$ dimensions.

From Eq. (III.43) one concludes that the $\kappa, \dot{\lambda}(t), \epsilon$, and $\rho_i(t)$–transformations give rise to conformal Killing vectors$^5$, while the $\nu, \mu_i(t), \omega_{ij}$ and constant $\lambda$–transformations generate Killing vectors. In view of (III.42), the isometry with constant parameter $\lambda$ is only possible for constant frequency $\Omega$ 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 $\lambda(t), \rho_i(t)$ and $\mu_i(t)$ are not independent. Obviously, $\rho_i(t)$ differs from $\mu_i(t)$ by an inessential constant. If $\mu$ satisfies the first equation in (III.42), then $\lambda(t) = \mu^2(t)$ will satisfy the second. Moreover, if $\mu_1$ and $\mu_2$ are independent solutions of the first equation, then the three independent solutions of the second are given by $\mu_1^2, \mu_2^2, \mu_1\mu_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

$$ds^2 = -\gamma^2 x_i x_i d\eta^2 - d\eta dv + dx_i dx_i, \qquad (\text{III.44})$$

<footnote>$^5$ Note that a combination of the $\kappa$ and $\lambda$-transformations may result in a Killing vector field, provided $\dot{\lambda} = -\kappa$. In view of (III.42), this only happens if $\Omega(t) = g/t$, where $g$ is a constant.</footnote>
will be automatically transmitted to become (conformal) isometries of (III.35), since the two expressions only differ by the cosmic scale factor $a(\eta)$. In particular, all the symmetry transformations of (III.44), which do not involve $\eta$ explicitly, will be transformed into the Killing vectors of (III.35), while those affecting $\eta$ 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 [3, 36, 37].

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.

The explicit form of the (conformal) isometries of (III.44) follows from Eqs. (III.41), (III.42) at $\Omega = 1$ after the substitution $t \to \eta$. In particular, conformal time translations, $\eta \to \eta + \theta$, are now symmetries, because $\gamma = \text{const.}$ in (III.44).

From the first equation in (III.42) one finds,

$$ \mu_i(\eta) = \alpha_i \cos(\eta) + \beta_i \sin(\eta), \qquad (III.45) $$

where the infinitesimal parameters $\alpha_i$ and $\beta_i$ are associated with spatial translations and Newton-Hooke boosts [49]. The second equations generate

$$ \rho_i(\eta) = \tilde{\alpha}_i \cos(\eta) + \tilde{\beta}_i \sin(\eta), \qquad (III.46) $$

which involve the infinitesimal parameters $\tilde{\alpha}_i$, $\tilde{\beta}_i$ and provide contributions to Eqs. (III.41) which are nonlinear in $x_i$. The third equation yields

$$ \lambda(\eta) = \theta + \sigma \cos(2\eta) + \rho \sin(2\eta), \qquad (III.47) $$

where the infinitesimal parameters $\theta$, $\sigma$, $\rho$, linked to time translations, special conformal transformations, and dilatations form an $so(2, 1)$ subalgebra. Along with spatial rotations described by $\omega_{ij}$, the $\mu_i(\eta)$ and $\lambda(\eta)$-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. [50, 51].

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 $\kappa$, and $\epsilon$-transformations.

Note that the $\kappa$-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 [52]. The $\epsilon$-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 $v$-direction, $v' = v + \nu$, they form an $so(2, 1)$ subalgebra. Interestingly enough and extending the Galilei-Carroll “duality” [53], 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 identified 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 less 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 $\nu$, $\kappa$ and $\epsilon$ in (III.41) will be broken, allowing us to recover the Schrödinger symmetry found in [46]. The generators are conveniently identified using the formulae in sec. 3 of [47].

## E. Geodesic motion in Ermakov-Milne-Pinney cosmoi

Having established a link between the Ermakov-Milne-Pinney equation and the Eisenhart-Duval lift, let us study the geodesic motion in Ermakov-Milne-Pinney cosmoi. The analysis is facilitated by switching to the conformal time (III.35) which allows one to solve the geodesic equations by quadrature,

$$ \tau = \kappa \int_{\eta_0}^{\eta} a^2(\tilde{\eta}) d\tilde{\eta} + \tau_0, \qquad (III.48a) $$

$$ x_i(\eta) = \alpha_i \cos(\gamma\eta) + \beta_i \sin(\gamma\eta), \qquad (III.48b) $$

$$ v(\eta) = v_0 - \frac{1}{2}(\alpha^2 - \beta^2)\gamma \sin(2\gamma\eta) + \alpha\beta\gamma \cos(2\gamma\eta) - \epsilon\kappa^2 \int_{\eta_0}^{\eta} a^2(\tilde{\eta}) d\tilde{\eta}, \qquad (III.48c) $$

where $\tau_0$, $\kappa$, $\alpha_i$, $\beta_i$, $v_0$ are constants of integration and $\alpha^2 = \alpha_i\alpha_i$, $\alpha\beta = \alpha_i\beta_i$. It is

<figure><img src="image_1.png" /><figcaption>FIG. 1: Null geodesics in Ermakov-Milne-Pinney cosmoi are closed loops which wrap around elliptic cylinders. The graph lies in three-dimensional space with $(x_1, x_2, v)$ parametrizing the two horizontal and the vertical axes, respectively.</figcaption></figure>

assumed that $\epsilon = 0$ for null geodesics and $\epsilon = -1$ for time-like geodesics.
Before proceeding, it is worth mentioning that, since for null geodesics Eqs. (III.48) admit a particular solution

$$x_i(\eta) = 0, \quad v(\eta) = v_0, \qquad (III.49)$$

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.48) defines η as an implicit function of the proper time τ. Although in most cases of interest the integral $\int_{\eta_0}^{\eta} a^2(\tilde{\eta})d\tilde{\eta}$ cannot be evaluated exactly, Eqs. (III.48) prove to be sufficient to comprehend a qualitative behaviour of geodesics. Indeed, the second line in (III.48) defines an ellipse. By making use of rotational invariance, one can set the ellipse to lie, say, in the $x_1x_2$-plane. As v evolves with time, geodesics in the Ermakov-Milne-Pinney cosmology wrap around the elliptic cylinder, v being its axis.

For null geodesics ($\epsilon = 0$), the trajectory is a closed loop and the motion is periodic (see FIG. 1). For time-like geodesics ($\epsilon = -1$), the orbit wraps around the cylinder remaining in a compact region of space for some time, then the $\epsilon$-term in the last line in (III.48) 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.

<figure><img src="image_2.png" /><figcaption>FIG. 2: Time-like geodesics in Ermakov-Milne-Pinney cosmology wrap around elliptic cylinders. They remain in a compact region of space for some time before escaping. The graph is given in three-dimensional space, ($x_1, x_2, v$) parametrizing the two horizontal and one vertical axes, respectively.</figcaption></figure>

Little is known about analytic solutions to the Ermakov-Milne-Pinney equation. Assuming that the energy density $\Omega^2(t)$ decreases with time (which happens in an expanding
universe), one has (see, e.g., [22] and Sect. 2.4.2 in [32])<fnref n="6" />

$$ \Omega^2(t) = \frac{h^{n-2}}{4\nu^2 t^n}, \quad \nu = \frac{1}{2-n}, \quad n \neq 2, \qquad (\text{III.50a}) $$

$$ a(t) = \sqrt{\pi\gamma\nu t(J_{\nu}^{2}(\lambda) + Y_{\nu}^{2}(\lambda))}, \quad \lambda = \left(\frac{t}{h}\right)^{\frac{1}{2\nu}}, \qquad (\text{III.50b}) $$

where $n$ is a rational number such that $\nu$ is positive, $h$ is a positive constant, and $J_\nu, Y_\nu$ 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 \to 0$.

Other interesting examples are provided by negative integer $\nu$; then $\nu$ should be replaced by $|\nu|$ in the expression for $a(t)$ in Eq. (III.50). 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 $\ddot{a} = 0$, when the solution is expressed in terms of elementary functions

$$ \Omega(t) = \frac{h}{t^2}, \quad a(t) = \sqrt{\frac{\gamma}{h}} t \quad \Rightarrow \quad \Omega(t) = \frac{\gamma}{a(t)^2}. \qquad (\text{III.51}) $$

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,

<figure><img src="image_3.png" /><figcaption>FIG. 3: A graph of the cosmic scale factor in the Ermakov-Milne-Pinney cosmological model at $n = 3$.</figcaption></figure>

$$ t = t_0 + \tau/\kappa, \qquad (\text{III.52a}) $$

$$ x_i(\tau) = \alpha_i \cos(h/t) + \beta_i \sin(h/t), \qquad (\text{III.52b}) $$

$$ v(\tau) = v_0 - \epsilon\kappa^2 t - \alpha\beta\gamma \cos(2h/t) + \frac{1}{2}(\alpha^2 - \beta^2)\gamma \sin(2h/t), \qquad (\text{III.52c}) $$

where $t_0, \kappa, \alpha_i, \beta_i, v_0$ are constants of integration.

<footnote n="6">⁶ From now on we switch back to cosmic time $t$.</footnote>
Although the graphs representing the cosmic scale factors in the Ermakov-Milne-Pinney cosmology look quite reasonable, the geodesic motion is apparently unrealistic. This happens because the second line in (III.48) defines an ellipse. An obvious cure is to generalize the construction to the case of time dependent $\gamma$ which will alter the qualitative behaviour of geodesics. It proves sufficient to replace $\gamma$ in (III.30) by an arbitrary function $\gamma(t)$ which will then show up on the right hand side of Eq. (III.31), viewed as an *algebraic equation* to fix $\gamma(t)$ in terms of $a(t)$ and $\Omega(t)$. In this way one can model a reasonable geodesic behavior in the generalized Ermakov-Milne-Pinney cosmology 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 $\gamma(t)^2$, finding an analytic solution may be complicated.

# IV. FRIEDMANN-LEMAÎTRE-ROBERTSON-WALKER SPACETIMES

There has been considerable interest in the past few years in applying the ideas circling around the Ermakov-Milne-Pinney equation to problems in cosmology, both classical and quantum, [7–14]. 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

$$g_{\mu\nu}dx^{\mu}dx^{\nu} = -dt^2 + a^2(t)d\Sigma_K^2, \qquad (IV.1)$$

where $d\Sigma_K^2$ is a 3-dimensional metric of constant curvature $K$. The scale factor $a(t)$ is governed by three equations of which, if $\dot{a} \neq 0$, only two are independent. These are the *Raychaudhuri equation*

$$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3}(\rho + 3P) + \frac{\Lambda}{3}, \qquad (IV.2)$$

the *Friedmann equation*

$$\left(\frac{\dot{a}}{a}\right)^2 + \frac{K}{a^2} = \frac{8\pi G}{3}\rho + \frac{\Lambda}{3}, \qquad (IV.3)$$

and the *conservation equation*

$$\dot{\rho} + 3\frac{\dot{a}}{a}(\rho + P) = 0, \qquad (IV.4)$$

where $\rho$ is the energy density, $P$ the pressure and $\Lambda$ the cosmological constant.

## A. 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

$$P = (\gamma - 1)\rho = w\rho, \qquad (IV.5)$$
where $\gamma$, or equivalently $w$, is constant. It follows from (IV.4) that

$$ \rho = \frac{\rho_0}{a^{3\gamma}}. \qquad (IV.6) $$

One has the following special cases

* $\gamma = 1$ corresponds to pressure-free matter.

* $\gamma = \frac{4}{3}$ to radiation

* $\gamma = 0$ to a cosmological constant which is equivalent to a constant energy density
  $$ \rho = \frac{\Lambda}{8\pi G}. $$

* $\gamma = 2$ is the largest value of $\gamma$ consistent with the dominant energy condition. If the energy density is positive it is sometimes known as “stiff matter”. If $\gamma > 2$ there exist solutions for which the scale factor blows up in finite time.

* $\gamma = \frac{2}{3}$ corresponds to a gas of cosmic strings

* $\gamma = \frac{1}{3}$ corresponds to a gas of membranes

* If $\gamma < \frac{2}{3}$, there exist solutions exhibiting a “big rip”, that is for which the scale factor blows up in finite time [54].

Another simple model is to suppose that the matter consists of a scalar field $\phi$ with potential $V(\phi)$

$$ \rho = \frac{1}{2}\dot{\phi}^2 + V(\phi), \quad P = \frac{1}{2}\dot{\phi}^2 - V(\phi). \qquad (IV.7) $$

It follows from (IV.4) that

$$ \ddot{\phi} + 3\frac{\dot{a}}{a}\dot{\phi} + V'(\phi) = 0. \qquad (IV.8) $$

If $|V(\phi)| \ll \frac{1}{2}\dot{\phi}^2$ this behaves like stiff matter, $\gamma = 2$. If $|V(\phi)| \gg \frac{1}{2}\dot{\phi}^2$ this behaves like a cosmological constant.

## B. 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 $\eta$ called conformal time, defined by

$$ \eta = \int \frac{dt}{a(t)}, \qquad (IV.9) $$

in terms of which the FLRW metric is manifestly conformally static

$$ ds^2 = a^2(t) \{-d\eta^2 + d\Sigma_K^2\}. \qquad (IV.10) $$
If $\chi$ 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 [55], 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 [56].

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 [57–59]. One obvious choice is the temperature or equivalently the redshift

$$T = \frac{a(t_0)}{a(t)} T_0 = T_0(1+z), \qquad (IV.11)$$

where $T_0 \approx 3K$ and the subscript 0 denotes the present epoch. Now if $a(t)$ runs from 0 at the Big Bang to $\infty$ in the infinite future then $T$ runs from $\infty$ to 0. In fact given three successive instants of time $t_1 \prec t_2 \prec t_3$ the temperature changes, $T_{12} = \frac{T(t_1)}{T(t_2)}$ etc behave multiplicatively

$$T_{13} = T_{12}T_{23}, \qquad (IV.12)$$

but taking logarithms $T = e^{-\Omega}$ gives an additive scale [58, 59]:

$$\Omega_{13} = \Omega_{12} + \Omega_{23}, \qquad (IV.13)$$

which runs from $-\infty$ to $+\infty$. Since from (IV.3) we have

$$d\Omega = \frac{da}{a} = \frac{\dot{a}}{a}dt \qquad (IV.14)$$

the FLRW metric takes the Eternal Peripatetic Form

$$g_{\mu\nu}dx^{\mu}dx^{\nu} = -\frac{d\Omega^2}{\frac{8\pi G\rho}{3} + \frac{\Lambda}{3} - K^2e^{-2\Omega}} + e^{2\Omega}d\Sigma_K^2. \qquad (IV.15)$$

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.

For a perfect fluid with constant barotropic index and vanishing cosmological constant, another temporal diffeomorphism defined by

$$d\tau = \frac{2}{3\gamma - 2}2adt \qquad (IV.16)$$

combined with the definition

$$x = \sqrt{\frac{3}{8\pi G\rho_0}} a^{\frac{3\gamma-2}{2}} \iff a = \left(\frac{8\pi G\rho_0}{3}x^2\right)^{\frac{1}{3\gamma-2}} \qquad (IV.17)$$

reduces the Friedmann equation (IV.3) to that of a time-independent harmonic oscillator

$$(\frac{dx}{d\tau})^2 + Kx^2 = 1. \qquad (IV.18)$$
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 $(\tau, x)$ [49]. The systems are invariant under the conformal extensions of Newton-Hooke-groups $NH_{\pm}$ [49, 51].

A related construction works if $K = 0$ but we include a cosmological constant term or dark energy term.

If $\Lambda > 0$ multiplication of (IV.3) by $a^{3\gamma}$ leads to an upside down simple harmonic oscillator equation for $a^{\frac{3\gamma}{2}}$ with solutions of the form

$$a = \left[ \sqrt{\frac{8\pi G\rho_0}{\Lambda}} \sinh\left(\frac{3\gamma}{2}\sqrt{\frac{\Lambda}{3}}t\right) \right]^{\frac{2}{3\gamma}}. \qquad (\text{IV.19})$$

One may calculate the jerk [60]

$$j = \frac{a^2 \frac{d^3 a}{dt^3}}{\left(\frac{da}{dt}\right)^3} = 1 + \frac{9\gamma(\gamma - 1)}{2 \cosh^2\left(\frac{3\gamma}{2} \sqrt{\frac{\Lambda}{3}} t\right)}. \qquad (\text{IV.20})$$

If $\gamma = 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, $\gamma = \frac{4}{3}$, the jerk starts from 3.

If $\Lambda < 0$ the motion is converted to a standard harmonic oscillator. For radiation $\gamma = \frac{4}{3}$ and (IV.2) becomes

$$\ddot{a} + \frac{|\Lambda|}{3}a = -\frac{8\pi G\rho_0}{a^3}, \qquad (\text{IV.21})$$

which is a variant of the Ermakov-Milne-Pinney equation [11]. Equation (IV.19) gives the particular solution:

$$a = \left[ \sqrt{\frac{8\pi G\rho_0}{|\Lambda|}} \sin\left(2\sqrt{\frac{|\Lambda|}{3}}t\right) \right]^{\frac{1}{2}} = \left[ \sqrt{\frac{8\pi G\rho_0}{-\Lambda}} 2 \sin\left(\sqrt{\frac{|\Lambda|}{3}}t\right) \cos\left(\sqrt{\frac{|\Lambda|}{3}}t\right) \right]^{\frac{1}{2}}. \quad (\text{IV.22})$$

One may check that the solution (IV.22) is a special case of the general solution given by (II.11) and (II.12).

The instance $\gamma = \frac{2}{3}$ or $w = \gamma - 1 = -\frac{1}{3}$ and $\Lambda < 0$ gives the special $k = 0$ case of the so-called simple harmonic universe [14]. It takes the form

$$a = \frac{4\pi G\rho_0}{|\Lambda|} \left(1 - \cos\left(\sqrt{\frac{|\Lambda|}{3}}t\right)\right). \qquad (\text{IV.23})$$

The general solution of (IV.2) subject to the constraint (IV.3) and an origin of time is

$$a = \frac{4\pi G\rho_0}{|\Lambda|} - \sqrt{\left(\frac{4\pi G\rho_0}{|\Lambda|}\right)^2 - \frac{3K}{|\Lambda|}} \cos\left(\sqrt{\frac{|\Lambda|}{3}}t\right). \qquad (\text{IV.24})$$

If $K = 0$ this coincides with our previous solution; if $K > 0$ it is the non-singular cyclic universe of [14].

For more examples of solutions of the Friedmann and Raychaudhuri equations the reader is referred to [61–63].
### C. Lifting the Friedmann-Lemaître-Robertson-Walker geodesics

Having presented the main features of FLRW cosmologies, now we show how the Eisenhart-Duval lift can be applied to geodesic motion in FLRW spacetimes and how to derive its dynamical symmetries. We then extend the result to the Dmitriev-Zel'dovich equations.

Our starting point is the geodesic Hamiltonian

$$H_{FLRW} = \frac{1}{2a^2} (-p_{\eta}^2 + h^{ij} p_i p_j), \qquad (IV.25)$$

where $\eta$ is the conformal time introduced in the previous section, and $h_{ij}$, $i, j = 1, \dots, 3$ is the metric on the 3-dimensional space. Setting $H_{FLRW} = -m^2/2$ and solving for $p_{\eta}$, we get:

$$p_{\eta} = \sqrt{m^2 a^2(\eta) + h^{ij} p_i p_j}. \qquad (IV.26)$$

We identify the low-energy regime as $\frac{p_i p_j}{m^2 a^2(\eta)} \ll 1$. In an expanding universe $p_i$ is constant. Therefore the dynamics will, unlike in the usual non-relativistic case, always enter this regime for large enough $\eta$. In this limit

$$p_{\eta} = ma(\eta) + \frac{h^{ij} p_i p_j}{2ma(\eta)}, \qquad (IV.27)$$

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 [64].

In our case the process transforms (IV.27) into

$$p_{\eta} p_{v} = ma(\eta) p_{v}^{2} + \frac{h^{ij} p_{i} p_{j}}{2ma(\eta)}. \qquad (IV.28)$$

This dynamical relation can be recovered by studying the null geodesics of the Hamiltonian

$$\mathcal{H} = \frac{h^{ij} p_i p_j}{2ma(\eta)} - p_{\eta} p_v + ma(\eta) p_v^2, \qquad (IV.29)$$

which is associated with the metric

$$ds^2 = ma(\eta)h^{ij}dx_i dx_j - 2ma(\eta)d\eta^2 - 2d\eta dv. \qquad (IV.30)$$

After a translation, $v = \tilde{v} - m \int^{\eta} a(\eta')d\eta'$, we obtain

$$ds^2 = ma(\eta)h^{ij}dx_i dx_j - d\eta d\tilde{v}. \qquad (IV.31)$$

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_{FLRW}$ in (IV.25).
Another consequence is that the null geodesics of the Eisenhart-Duval lift just constructed are the same as the null geodesics of the Ermakov-Milne-Pinney spacetime, described in sec.III C.

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

$$ \frac{1}{2a^2} \left[ -p_{\eta}^2 + \sum_{A=1}^{N} h^{ij} p_i^A p_j^A \right] - \sum_{1 \le A < B \le N} \frac{G m_i m_j}{a |\vec{r}^A - \vec{r}^B|}, \qquad (IV.32) $$

where the $p_i^A$ are the spatial momenta of the individual particles and $\vec{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 since the spacetime is spatially homogeneous. As before, we set the Hamiltonian equal to a constant value $-M^2/2$, calculate the low-energy limit and make the relationship homogeneous in the momenta. We obtain the null Hamiltonian

$$ \mathcal{H} = -p_{\eta}p_{\bar{\nu}} + \frac{1}{2Ma} \sum_{A=1}^{N} h^{ij} p_i^A p_j^A - \frac{1}{M} \sum_{1 \le A < B \le N} \frac{G m_i m_j}{|\vec{r}^A - \vec{r}^B|} p_{\bar{\nu}}^2, \qquad (IV.33) $$

where again we have redefined the $v$ variable in order to remove the term $Ma(\eta)p_{\bar{\nu}}^2$. Going back to cosmic time, with $p_{\eta} = ap_t$, we obtain

$$ \mathcal{H} = a(t) \left[ -p_t p_{\bar{\nu}} + \frac{1}{2Ma^2} \sum_{A=1}^{N} h^{ij} p_i^A p_j^A - \frac{1}{Ma} \sum_{1 \le A < B \le N} \frac{G m_i m_j}{|\vec{r}^A - \vec{r}^B|} p_{\bar{\nu}}^2 \right], \qquad (IV.34) $$

which belongs to the same conformal class as the Eisenhart-Duval lift of the Hamiltonian presented in [66].

This gives rise to the Dmitriev-Zel'dovich equations [49, 65, 66]. If we consider, as is done in almost all numerical studies of large scale structure, the motion of N non-uniformities about a FLRW background metric, they are governed by the Dmitriev-Zel'dovich equations [65, 66]

$$ \frac{d}{dt} (a^2(t) \dot{\mathbf{x}}_a) = \frac{1}{a(t)} \sum_{b \neq a} \frac{G m_b (\mathbf{x}_b - \mathbf{x}_a)}{|\mathbf{x}_b - \mathbf{x}_a|} \qquad (IV.35) $$

where $a, b = 1, 2, \dots, N$ label the non-uniformities and $m_a$ are their masses. It is simple to verify that if one replaces every $\mathbf{x}_a$ in (IV.35) by

$$ \mathbf{x}_a + \mathbf{a} + \mathbf{b} \int_0^t \frac{dt'}{a(t')} \qquad (IV.36) $$

then the Dmitriev-Zel'dovich equations are left invariant.

The equations (IV.35) may be derived from a Lagrangian [66, 67]

$$ L = \frac{1}{2} \sum_{1 \le a \le N} a^2(t) m_a \dot{\mathbf{x}}_a^2 + \frac{1}{a(t)} \sum_{1 \le a < b \le N} \frac{G m_a m_b}{|\mathbf{x}_a - \mathbf{x}_b|}. \qquad (IV.37) $$
This system Eisenhart-Duval lifts to

$$ds^2 = \sum_{1 \le a \le N} a^2(t) m_a d\mathbf{x}_a^2 + 2 dt dv + 2 \frac{1}{a(t)} \sum_{1 \le a < b \le N} \frac{G m_a m_b}{|\mathbf{x}_a - \mathbf{x}_b|} dt^2. \quad (IV.38)$$

Up to signs and factors this coincides with the metric (III.30) as long as $\gamma = 0$ and $G = 0$. In other words a free non-relativistic particle in a background Friedmann universe lifts to a null geodesic of (III.30) provided $\gamma = 0$. In fact the time-dependent harmonic oscillator in (III.30) looks very much like a cosmological constant term.

### D. The Eisenhart-Duval lift of the Friedmann equations

The action density of Einstein gravity coupled to a scalar field is

$$\mathcal{L} = \left\{ \frac{1}{16\pi G} R - \frac{1}{2} g^{\mu\nu} \partial_{\nu} \phi \partial_{\mu} \phi - V(\phi) \right\}. \quad (IV.39)$$

Substituting in the FLWR metric and assuming $\phi$ depends on $t$ we find that *up to a total time derivative*

$$\mathcal{L} = \frac{6}{16\pi G} (-a\dot{a}^2 + Ka) + a^3(\frac{1}{2}\dot{\phi}^2 + V(\phi)). \quad (IV.40)$$

Varying $\mathcal{L}$ with respect to $\phi$ and $a$ yields the field equation for $\phi$ and the following second order ODE for $a$

$$\frac{\ddot{a}}{a} + \frac{\dot{a}^2}{a^2} = -\frac{4\pi G}{3}(\dot{\phi}^2 - 4V(\phi)). \quad (IV.41)$$

One may check that (IV.41) is a linear combination of the second order Raychaudhuri equation and the first order Friedmann equation. If one calculates the Hamiltonian density $\mathcal{H}$ of $\mathcal{L}$ one finds that the two canonical momenta are

$$p_{\phi} = a^3 \dot{\phi}, \qquad p_a = -\frac{3}{4\pi G} a \dot{a}, \quad (IV.42)$$

and considered as a function of the velocities $(\dot{a}, \dot{\phi})$, we find that

$$\mathcal{H} = \frac{6}{16\pi G} (-a\dot{a}^2 - Ka) + a^3(\frac{1}{2}\dot{\phi}^2 + V(\phi)). \quad (IV.43)$$

One may check that the Friedmann ODE is the same as the constraint

$$\mathcal{H} = 0. \quad (IV.44)$$

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 $\mathcal{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

$$ds^2 = -N^2(t)dt^2 + a^2(t)d\Sigma_K^2 \quad (IV.45)$$
which however differs from our previous form by no more than a temporal diffeomorphism or clock regraduation. The function $N(t)$ is called the lapse function. It measures the lapse of time relative to cosmic time. Now

$$ \tilde{\mathcal{L}} = \frac{6}{16\pi G}(-N^{-1}a\dot{a}^2 + NKa) + a^3(N^{-1}\frac{1}{2}\dot{\phi}^2 + NV(\phi)). \quad (IV.46) $$

We see that the lapse function $N$ enters (IV.46) as a Lagrange multiplier and varying the Lagrangian $\tilde{\mathcal{L}}$ with respect to $N$ and then setting $N = 1$ yields the Friedmann equation, i.e. the Hamiltonian constraint (IV.44). In fact in terms of the canonical momenta we have$^7$

$$ \mathcal{H} = N\left\{-\frac{2\pi G}{3}\frac{p_a^2}{a} - \frac{3K}{8\pi G} + \frac{1}{2}\frac{p_\phi^2}{a^3} + a^3V(\phi)\right\}. \quad (IV.47) $$

We are now in the situation to perform an Eisenhart-Duval lift or Marsden-Weinstein oxidation, we set

$$ p_t = -\mathcal{H} \quad (IV.48) $$

and introduce an ignorable momentum $p_v$ and rewrite (IV.47) as

$$ 2H_{\text{Eisenhart-Duval}} = 2p_t p_v + 2N \left\{ -\frac{2\pi G}{3} \frac{p_a^2}{a} - \frac{3K}{8\pi G} p_v^2 + \frac{1}{2} \frac{p_\phi^2}{a^3} + a^3 V(\phi) p_v^2 \right\} = 0. \quad (IV.49) $$

Note that if we set $N = 1$ and compute the 4-metric then the coefficients of $da^2$ and $d\phi^2$ have the opposite sign. If we regard $N$ as a further coordinate, then we have a Galilean or Newton-Cartan metric.

# V. 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.

* The celebrated Ermakov-Lewis invariant of a time-dependent harmonic oscillator can be obtained in purely geometric way by applying the Eisenhart-Duval lift.

* A cosmological extension of the Eisenhart-Duval metric is constructed by properly incorporating into the scheme a cosmic scale factor and the energy-momentum tensor. The evolution of 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.

<footnote>$^7$ From now on we drop the tildes.</footnote>
• The algebra of vector fields which involve both Killing and conformal Killing vec-
tors is identified with so(2, 2 + d), the conformal Newton-Hooke algebra being its
subalgebra.

• Geodesic motion in Ermakov-Milne-Pinney cosmoi is described.

* The Eisenhart-Duval lift of geodesics in the Friedmann-Lemaître-Robertson-Walker spacetimes is found and then generalized to the Dmitriev-Zel'dovich equations.

* The derivation of the Friedmann equations within the framework of the Eisenhart-Duval lift is presented.

Acknowledgments

We would like to thank an anonymous referee for informing us that Ermakov’s paper
was preceeded by one of Adolph Steen [68, 69] and therefore it would be more appropriate
to speak of the Steen-Ermakov-Milne-Pinney equation. MC was funded by the CNPq
under project 303923/2015-6, and by a *Pesquisador Mineiro* project n. PPM-00630-17.
AG was supported by the Tomsk Polytechnic University competitiveness enhancement
program.

[1] L. P. Eisenhart, “Dynamical trajectories and geodesics,” *Annals Math. 30* (1929) 591.

[2] C. Duval, G. Burdet, H. Künzle, M. Perrin, “Bargmann structures and Newton–Cartan theory,” *Phys. Rev. D* **31** (1985) 1841; C. Duval, G.W. Gibbons, P.A. Horvathy, “Celestial mechanics, conformal structures and gravitational waves,” *Phys. Rev. D* **43** (1991) 3907, [arXiv:hep-th/0512188].

[3] M. Cariglia, C. Duval, G. W. Gibbons, P. A. Horvathy, “Eisenhart lifts and symmetries of time-dependent systems,” *Annals Phys. 373* (2016) 631, doi:10.1016/j.aop.2016.07.033, [arXiv:1605.01932 [hep-th]].

[4] V. P. Ermakov, “Second order differential equations. Conditions of complete integrability,” *Univ. Izv. Kiev, Series III 9* (1880) 1 (English translation: A. O. Harin, under redaction by P. G. L. Leach, *Appl. Anal. Discrete Math. 2* (2008) 123, doi:10.2298/AADM0802123E).

[5] W. E. Milne, “The numerical determination of characteristic numbers,” *Phys. Rev. 35* (1930) 863.

[6] E. Pinney, “The nonlinear differential equation $y'' + p(x)y + \frac{c}{y^3} = 0,$” *Proc. Amer. Math. Soc.* **1** (1959) 681.
[7] H. C. Rosu, “Ermakov approach for empty FRW minisuperspace oscillators,” [arXiv:gr-qc/9806075].

[8] H. C. Rosu, J. Socorro, “Ermakov approach for minisuperspace oscillators,” *Int. J. Theor. Phys.* **41** (2002) 39, doi:10.1023/A:1013265331607, [arXiv:gr-qc/9908028].

[9] H. C. Rosu, P. Espinoza, M. Reyes, “Ermakov approach for $Q = 0$ empty FRW minisuperspace oscillators,” *Nuovo Cim. B* **114** (1999) 1439, [arXiv:gr-qc/9910070].

[10] H. C. Rosu, P. Espinoza, “An Ermakov study of $Q$ does not equal zero EFRW minisuperspace oscillators,” [arXiv:gr-qc/9912033].

[11] R. M. Hawkins, J. E. Lidsey, “The Ermakov-Pinney equation in scalar field cosmologies,” *Phys. Rev. D* **66** (2002) 023523, doi:10.1103/PhysRevD.66.023523, [arXiv:astro-ph/0112139].

[12] F. L. Williams, P. G. Kevrekidis, “On (2+1)-dimensional Friedmann-Robertson-Walker universes: An Ermakov-Pinney equation approach,” *Class. Quant. Grav.* **20** (2003) L177, doi:10.1088/0264-9381/20/13/101.

[13] H. C. Rosu, P. B. Espinoza, “Ermakov-Lewis angles for one-parameter supersymmetric families of Newtonian free damping modes,” [arXiv:physics/0004014v4 [physics.class-ph]].

[14] P. W. Graham, B. Horn, S. Kachru, S. Rajendran, G. Torroba, “A simple harmonic Universe,” *JHEP* **1402** (2014) 029, doi:10.1007/JHEP02(2014)029, [arXiv:1109.0282v2 [hep-th]].

[15] A. Galajinsky, “Geometry of the isotropic oscillator driven by the conformal mode,” *Eur. Phys. J. C* **78** (2018) 72, doi:10.1140/epjc/s10052-018-5568-8, [arXiv:1712.00742v3 [hep-th]].

[16] R. M. Morris, P. G. L. Leach, “The Ermakov-Pinney Equation: its varied origins and the effects of the introduction of symmetry-breaking functions,” [arXiv:1510.08992v1 [math.CA]].

[17] P. B. Espinoza, *Ermakov-Lewis dynamic invariants*, [arXiv:math-ph/0002005v3].

[18] R. B. Dingle, “The method of comparison equations in the solution of linear second-order differential equations (generalized w.k.b. method),” *Appl. Sci. Res. B* **5** (1956) 345.

[19] A. Kamenshchik, G. Venturi, “Remarks on the method of comparison equations (generalized WKB method) and the generalized Ermakov-Pinney equation,” *Russ. Phys. J.* **52** (2009) 1339, doi:10.1007/s11182-010-9375-4, [arXiv:math-ph/0506017v2].

[20] V. de Alfaro, S. Fubini, G. Furlan, “Conformal invariance in quantum mechanics,” *Nuovo*
*Cim. A.* **34** (1976) 569.

[21] H. R. Lewis, “Classical and quantum systems with time-dependent harmonic oscillator-type Hamiltonians,” *Phys. Rev. Lett.* **18** (1967) 510.

[22] H. R. Lewis, “Class of exact invariants for classical and quantum time-dependent harmonic oscillators,” *J. Math. Phys.* **9** (1968) 1976, https://doi.org/10.1063/1.1664532.

[23] T. Padmanabhan, “Demystifying the constancy of the Ermakov-Lewis invariant for a time dependent oscillator,” [arXiv:1712.07328 [physics.class-ph]].

[24] H. R. Lewis “Motion of a time-dependent harmonic oscillator and of a charged particle in a time-dependent, axially symmetric electromagnetic field,” *Phys. Rev.* **172** (1968) 1313.

[25] H. R. Lewis, W. B. Riesenfeld, “An exact quantum theory of the time-dependent harmonic oscillator and of a charged particle in a time-dependent electromagnetic field,” *J. Math. Phys.* **10** (1969) 1458, https://doi.org/10.1063/1.1664991.

[26] G. W. Gibbons, C. N. Pope, “Kohn’s theorem, Larmor’s equivalence principle and the Newton-Hooke group,” *Annals Phys.* **326** (2011) 1760, doi:10.1016/j.aop.2011.03.003, [arXiv:1010.2455 [hep-th]].

[27] P. M. Zhang, G. W. Gibbons, P. A. Horvathy, “Kohn’s theorem and Newton-Hooke symmetry for Hill’s equations,” *Phys. Rev. D* **85** (2012) 045031, doi:10.1103/PhysRevD.85.045031, [arXiv:1112.4793 [hep-th]]; P. M. Zhang, P. A. Horvathy, K. Andrzejewski, J. Gonera, P. Kosinski, “Newton-Hooke type symmetry of anisotropic oscillators,” *Annals Phys.* **333** (2013) 335, doi:10.1016/j.aop.2012.11.018, [arXiv:1207.2875 [hep-th]].

[28] R. Ray, J. L. Reid, “More exact invariants for the time-dependent harmonic oscillator,” *Phys. Lett. A* **71** (1979) 317.

[29] R. Ray, J. L. Reid, M. Lutsky, “New nonlinear dynamical systems possessing invariants,” *Phys. Lett. A* **84** (1981) 42.

[30] C. J. Eliezer, A. Gray, “A note on the time-dependent harmonic oscillator,” *SIAM J. Appl. Math.* **30** (1976) 463.

[31] J. W. Barrett, G. W. Gibbons, M. J. Perry, C. N. Pope, P. Ruback, “Kleinian geometry and the N=2 superstring,” *Int. J. Mod. Phys. A* **9** (1994) 1457, doi:10.1142/S0217751X94000650, [arXiv:hep-th/9302073].

[32] A. D. Polyanin, V. F. Zaitsev, *Handbook of exact solutions for ordinary differential equations*, 2nd ed., CRC Press, 2003.
[33] Sir Harold Jeffreys, Bertha Swirles (Lady Jeffreys), “Methods of mathematical physics,”
Cambridge University Press, 3rd Edition, (1972).

[34] M. V. Ioffe, H. J. Korsch, “Nonlinear supersymmetric (Darboux) covariance of the Ermakov-Milne-Pinney equation,” *Phys. Lett. A* **311** (2003) 200, doi:10.1016/S0375-9601(03)00495-X, [arXiv:math-ph/0209013].

[35] D. Schuch, “Riccati and Ermakov equations in time-dependent and time-independent quantum systems,” *SIGMA* **4** (2008) 043, doi:10.3842/SIGMA.2008.043, [arXiv:0805.1687 [quant-ph]].

[36] J. Guerrero, F. F. Lopez-Ruiz, “The quantum Arnold transformation and the Ermakov-Pinney equation,” *Phys. Scr.* **87** (2013) 038105, doi:10.1088/0031-8949/87/03/038105, [arXiv:1302.1316 [math-ph]].

[37] F. F. Lopez-Ruiz, J. Guerrero, “Generalizations of the Ermakov system through the quantum Arnold transformation,” *J. Phys. Conf. Ser.* **538** (2014) 012015, doi:10.1088/1742-6596/538/1/012015.

[38] S. P. Kim and W. Kim, “Construction of exact Ermakov-Pinney solutions and time-dependent quantum oscillators,” *J. Korean Phys. Soc.* **69** (2016) 1513, doi:10.3938/jkps.69.1513, [arXiv:1609.00248 [quant-ph]].

[39] E. Madelung, “Quantentheorie in hydrodynamischer Form,” *Z. Phys.* **40** (1926) 322.

[40] R. A. Lee, “Quantum ray equations,” *J. Phys. A: Math. Gen.* **15** (1982) 2761.

[41] R. S. Kaushal, “Quantum analogues of the Ermakov systems and the phase of the quantum wave function,” *Int. J. Theor. Phys.* **40** (2001) 835.

[42] G. Herring, P. G. Kevrekidis, F. Williams, T. Christodoulakis, D. J. Frantzeskakis, “From Feshbach-resonance managed Bose-Einstein condensates to anisotropic universes: some applications of the Ermakov-Pinney equation with time-dependent nonlinearity,” *Phys. Lett. A* **367** (2007) 140, doi:10.1016/j.physleta.2007.01.087, [arXiv:cond-mat/0701756 [cond-mat.other]]; *Phys. Lett. A* **372** (2008) 277, doi.org/10.1016/j.physleta.2007.01.078.

[43] R. Ray, J. L. Reid, “New Ermakov systems from Lie symmetry theory,” *Phys. Lett. A* **89** (1982) 317.

[44] J. R. Ray, J. L. Reid, “Ermakov systems, Noether’s theorem and the Sarlet-Bahar method,” *Lett. Math. Phys.* **4** (1980) 235, doi:10.1007/BF00316679.

[45] J. R. Ray, J. L. Reid, “Noether’s theorem and Ermakov systems for nonlinear equations of
motion,” *Nuovo Cim. A* **59** (1980) 134, doi:10.1007/BF02902329.

[46] G. E. Prince, C. J. Eliezer, “Symmetries of the time-dependent N-dimensional oscillator,” *J. Phys. A: Math. Gen.* **13** (1980) 815.

[47] M. Hassaïne and P. A. Horváthy, “Field-dependent symmetries of a non-relativistic fluid model,” *Ann. Phys. (N. Y.)* **282**, 218 (2000) [math-ph/9904022].

[48] I. A. Pedrosa, “Canonical transformations and exact invariants for dissipative systems,” *J. Math. Phys.* **28** (1987) 2662, https://doi.org/10.1063/1.527707.

[49] G. W. Gibbons, C. E. Patricot, “Newton-Hooke space-times, Hpp waves and the cosmological constant,” *Class. Quant. Grav.* **20** (2003) 5225, doi:10.1088/0264-9381/20/23/016, [arXiv:hep-th/0308200].

[50] C. Duval and P. A. Horvathy, “Non-relativistic conformal symmetries and Newton-Cartan structures,” *J. Phys. A* **42**, 465206 (2009), doi:10.1088/1751-8113/42/46/465206, [arXiv:0904.0531 [math-ph]]; “Conformal Galilei groups, Veronese curves, and Newton-Hooke spacetimes,” *J. Phys. A* **44** (2011) 335203, doi:10.1088/1751-8113/44/33/335203, [arXiv:1104.1502].

[51] A. Galajinsky, “Conformal mechanics in Newton-Hooke spacetime,” *Nucl. Phys. B* **832** (2010) 586, doi:10.1016/j.nuclphysb.2010.02.023, [arXiv:1002.2290v2 [hep-th]].

[52] C. G. Torre, “Gravitational waves: Just plane symmetry,” *Gen. Rel. Grav.* **38** (2006) 653, doi:10.1007/s10714-006-0255-8, [arXiv:gr-qc/9907089].

[53] C. Duval, G. W. Gibbons, P. A. Horvathy and P. M. Zhang, “Carroll versus Newton and Galilei: two dual non-Einsteinian concepts of time,” *Class. Quant. Grav.* **31** (2014) 085016, doi:10.1088/0264-9381/31/8/085016, [arXiv:1402.0657 [gr-qc]].

[54] R. R. Caldwell, M. Kamionkowski, N. N. Weinberg, “Phantom energy and cosmic doomsday,” *Phys. Rev. Lett.* **91** (2003) 071301, doi:10.1103/PhysRevLett.91.071301, [arXiv:astro-ph/0302506].

[55] R. F. Marzke, J. A. Wheeler, “Gravitation as geometry I: the geometry of space-time and the geometrodynamical standard meter,” in *Gravitation and Relativity*, H.-Y. Chiu and W. F. Hoffmann, eds. (W. A. Benjamin, New York, Amsterdam, 1964), p. 40.

[56] T. Harada, B. J. Carr, T. Igata, “Complete conformal classification of the Friedmann-Lemaître-Robertson-Walker solutions with a linear equation of state,” [arXiv:1801.01966 [gr-qc]].
[57] C. W. Misner, “Absolute zero of time,” *Phys. Rev.* **186** (1969) 1328, doi:10.1103/PhysRev.186.1328.

[58] J-M. Lévy-Leblond, “The unbegun Big Bang,” *Nature* **342** (1980) 23.

[59] J-M. Lévy-Leblond, “Did the Big Bang begin ?” *Am. J. Phys.* **58** (1990) 156.

[60] M. Dunajski, G. W. Gibbons, “Cosmic jerk, snap and beyond,” *Class. Quant. Grav.* **25** (2008) 235012, doi:10.1088/0264-9381/25/23/235012, [arXiv:0807.0207 [gr-qc]].

[61] S. Chen, G. W. Gibbons, Y. Li, Y. Yang, “Friedmann’s equations in all dimensions and Chebyshev’s theorem,” *JCAP* **1412** (2014) 035, doi:10.1088/1475-7516/2014/12/035, [arXiv:1409.3352 [astro-ph.CO]].

[62] S. Chen, G. W. Gibbons, Y. Yang, “Explicit integration of Friedmann’s equation with nonlinear equations of state,” *JCAP* **1505** (2015) 020, doi:10.1088/1475-7516/2015/05/020, [arXiv:1502.05042 [gr-qc]].

[63] S. Chen, G. W. Gibbons, Y. Yang, “Friedmann-Lemaître cosmologies via roulettes and other analytic methods,” *JCAP* **1510** (2015) 056, doi:10.1088/1475-7516/2015/10/056, [arXiv:1508.06750 [gr-qc]].

[64] M. Cariglia, “Null lifts and projective dynamics”, *Ann. Phys.* **362** (2015) 642-58, doi:10.1016/j.aop.2015.09.002, [arXiv:1506.00714 [math-ph]].

[65] N A Dmitriev and Ya B Zel'dovich, “The Energy of Accidental Motions in an Expanding Universe” *Soviet Physics JETP* **18** (1964) 793

[66] G. F. R. Ellis and G. W. Gibbons, “Discrete Newtonian Cosmology: Perturbations”, *Class. Quant. Grav.* **32** (2015) 055001, doi:10.1088/0264-9381/32/5/055001, [arXiv:1409.0395 [gr-qc]].

[67] P. J. E. Peebles, “Tracing galaxy orbits back in time”, Astrophys. J. **344** (1989) 53-56

[68] A. Steen, “Om Formen for Integralet af den lineaere Differential aligning af an den Orden”, Overs. d. K. Danske Vidensk. Selsk. Forh. (1874), 1-12.

[69] R. Redheffer, I. Redheffer, “Steen’s 1874 paper: historical survey and translation”. Aequationes Math. **61** (2001) 131-15