Eisenhart Lift of $2$--Dimensional Mechanics

The Eisenhart lift is a variant of geometrization of classical mechanics with $d$ degrees of freedom in which the equations of motion are embedded into the geodesic equations of a Brinkmann-type metric defined on $(d+2)$-dimensional spacetime of Lorentzian signature. In this work, the Eisenhart lift of $2$-dimensional mechanics on curved background is studied. The corresponding $4$-dimensional metric is governed by two scalar functions which are just the conformal factor and the potential of the original dynamical system. We derive a conformal symmetry and a corresponding quadratic integral, associated with the Eisenhart lift. The energy--momentum tensor is constructed which, along with the metric, provides a solution to the Einstein equations. Uplifts of $2$-dimensional superintegrable models are discussed with a particular emphasis on the issue of hidden symmetries. It is shown that for the $2$-dimensional Darboux--Koenigs metrics, only type I can result in Eisenhart lifts which satisfy the weak energy condition. However, some physically viable metrics with hidden symmetries are presented.


Introduction
It is known since Eisenhart's work on the geometrization of classical mechanics [1] that any dynamical system with d degrees of freedom q i , i = 1, . . . , d, which is governed by the Lagrangian L, can be embedded into the geodesic equations of the Brinkmann-type metric 2Ldt 2 − dtdv, where t is the temporal variable and v is an extra coordinate. When analyzing the geodesic equations, one finds that t is affinely related to the proper time τ , the equations of motion for q i (t) coincide with those following from the Lagrangian L, while the evolution of v(t) is unambiguously fixed provided q i (t) are known. The initial dynamical system is thus recovered by implementing the null reduction along v. a e-mail: a.p.fordy@leeds.ac.uk b e-mail: galajin@tpu.ru Curiously enough, the original publication of [1] did not receive much attention by physicists and had soon fallen into oblivion. Several decades passed before the method was rediscovered in [2], which paved the way for various physical applications (see [3] and references therein).
Particularly interesting geometries result from uplifts of integrable and superintegrable systems. Constants of the motion polynomial in momenta give rise to Killing vectors and Killing tensors, the rank of the latter being equal to the degree of the polynomial. Killing vectors are associated with a clear statement of symmetry. They result from infinitesimal coordinate transformations which leave the form of a metric invariant. On the other hand, higher order Killing tensors (associated with hidden symmetries) have no such simple interpretation. Second order Killing tensors are associated with separation of variables of the Hamilton-Jacobi equation, with Carter's integration of the geodesic equations in the Kerr metric [4] being the prime example in General Relativity. Although quite a few physically meaningful spacetimes have been constructed, which admit one or several second rank Killing tensors (for a recent review see [5]), no solution to the vacuum Einstein equations admitting higher rank Killing tensors is presently known. 1 This empirical barrier of rank-2 seems rather puzzling. It may be a technical issue but perhaps something more fundamental lies behind it. The study of Lorentzian metrics admitting higher rank Killing tensors within Eisenhart's approach generated an extensive recent literature [8][9][10][11][12][13]. While the geometric reformulation of Newtonian mechanics brings mostly aesthetic advantages, the construction of Brinkmann-type metrics with hidden symmetries is a source of new results.
Thus far attention was mostly drawn to integrable mechanics in flat space (for some curved space examples see [8,9]).
Although the interrelationship between geometric characteristics of the Eisenhart metric and those of a Riemannian metric underlying mechanics on a curved background is generally rather complicated, the analysis is greatly simplified for 2D case, because 2-dimensional manifolds are conformally flat.
In recent years there has been a burst of activity in the identification and classification of superintegrable systems, both classical and quantum (see the review [14] and references therein). Most of the interest is in Hamiltonians which are in "natural form" (the sum of kinetic and potential energies), with the kinetic energy being quadratic in momenta and therefore associated with a (pseudo-)Riemannian metric. When an n-dimensional space is either flat or constant curvature, it possesses the maximal group of isometries, which is of dimension 1 2 n(n + 1). In this case, the kinetic energy is actually the second order Casimir function of the symmetry algebra (see [15]). Furthermore all higher order integrals of the geodesic equations are built out of the corresponding Noether constants by just taking polynomial expressions in them. Whilst most of the classification results and examples which occur in applications correspond to flat or constant curvature spaces, there are well known examples of conformally flat spaces (but not constant curvature), possessing quadratic invariants, which are clearly not just quadratic expressions in Noether constants. Specifically, there are the metrics found by Koenigs [16], which are described and analysed in [17,18]. There are other examples of conformally flat spaces (but not constant curvature), possessing one Noether constant and a cubic integral (classified in [19] and further studied and generalised in [20,21]), which again cannot be represented as a cubic expression in the isometry algebra. Being conformally flat, these spaces do have an abundant supply of conformal symmetries and in [22,23] a method was proposed for building quadratic and higher order invariants from appropriate polynomial expressions in conformal invariants.
The goal of this paper is to study the Eisenhart lift of 2-dimensional mechanics in curved space. Specifically, we consider the relationship of the curvature of the 2-, and 4dimensional geometries and the structure of the Einstein tensor. We also consider the conformal symmetries of the Eisenhart lift and use this to build an additional quadratic invariant. We are particularly interested in constructing physically admissible energy-momentum tensors in a purely geometric way and to derive equations which connect geometric characteristics to those of matter.
The work is organized as follows. In Sect. 2, the equations of motion of a generic 2dimensional mechanical system in curved space are embedded into the geodesic equations of a Brinkmann type 4dimensional metric. The latter is determined by two scalar functions which are just the conformal factor and the poten-tial of the original 2D mechanics. We see that first integrals can always be lifted from the 2 to the 4-dimensional domain, with the addition of 2 further involutive integrals, thus preserving Liouville integrability. In Sect. 3 we show that the Eisenhart lift has a conformal symmetry for a large class of 2D Hamiltonians. This conformal symmetry is then used to construct a new first integral for the Eisenhart lift.
In Sect. 4 we discuss the energy-momentum tensor of the 4-dimensional lift. Restrictions on the energy-momentum tensor, which follow from the weak and strong energy conditions, are formulated. In Sect. 5 we discuss some specific Hamiltonian systems, comparing Liouville integrability with superintegrability. It is shown that 2D Darboux-Koenigs metrics result in 4D solutions which violate the weak energy condition in types II, III, and IV, but obey it in the invariant region x > 0 for type I. A physically viable metric, admitting a rank 2 Killing tensor, is constructed by uplifting a superintegrable model on S 2 , as well as other models with the additional functional freedom of being just Liouville integrable.
Some final remarks are gathered in the concluding Sect. 6.

Eisenhart lift of 2-dimensional mechanics in curved space
In this section we describe the Eisenhart lift, which is similar to the Kaluza-Klein extension, which allows us to consider the motion of a particle in a curved background, under the influence of a potential, as geodesic motion on a larger curved space. First we give a brief review of ideas from geometric mechanics and the relationship between first integrals and Killing tensors.

First integrals in classical mechanics
A particle moving in a curved background, with metric coefficients g i j and scalar potential U (q), can be written in Lagrange form, with L = 1 2 n i, j=1 g i jq iq j − U (q) or, equivalently, in Hamiltonian form with H = 1 2 n i, j=1 g i j p i p j +U (q), where g i j are inverse metric coefficients and p i = n j=1 g i jq j . In either case, we obtain where k i j are the components of the Levi-Civita connection. The simplest way to compute the connection coefficients is to use the Hamiltonian to calculate {{q k , H }, H } and to replace p i by the above formula and then to just equate coefficients to those in (1).
In the Lagrangian framework we have Noether's theorem, which relates symmetries of the Lagrangian to first integrals (Noether constants). This connection is more transparent in the Hamiltonian formulation, being related to first integrals which are of first degree in momenta. In the case of geodesic equations, there is a direct relation to Killing vectors. Given then where ∇ i denotes the covariant derivative. Furthermore, Hamilton's equations for K generate a vector field on configuration space: since the first n components are written entirely in terms of the position variables.K is just the Killing vector corresponding to the Noether constant K . Similarly, if F = n i, j=1 f i j p i p j , f i j the components of a symmetric matrix, then thus defining a rank 2 Killing tensor. In general, if F is an integral which is homogeneously polynomial of degree m in momenta, then the coefficients define a Killing tensor of rank m. Such integrals no longer generate a vector field, such as (2c), on configuration space, so have a different geometric significance. A Killing vector defines a "motion" on configuration space, defined by the dynamical system under which the metric is invariant. For simple cases, such as rotations, this can be explicitly solved to give a transformation of coordinates. For each Killing vector K it is possible (in principle) to find new coordinates such that K = ∂ ∂ Q 1 , in which case the metric coefficients are independent of the variable Q 1 . If we have n commuting Killing vectors K i (the flat case), then it is possible to find coordinates Q 1 , . . . , Q n , such that K i = ∂ ∂ Q i , in which case metric coefficients are constant.
Whilst quadratic integrals have no such simple geometric meaning, they arise in the theory of separation of variables of the Hamilton-Jacobi equation. In fact, a complete solution of the Hamilton-Jacobi equation depends upon n parameters, which can then be written in terms of the dynamical variables to give n mutually commuting quadratic integrals (some of which could be squares of linear integrals).
Whilst n is the maximal number of independent functions which can be in involution, it is possible to have further integrals of the Hamiltonian H , which necessarily generate a non-Abelian algebra of integrals of H . The maximum number of additional independent integrals is n − 1, since the "level surface" of 2n−1 integrals (meaning the intersection of individual level surfaces) is just the (unparameterised) integral curve. Such systems are called superintegrable (maximal when there are 2n − 1 independent integrals). Well known elementary examples are the isotropic harmonic oscillator, the Kepler system and the Calogero-Moser system. The role of superintegrability in both the classical and quantum context is described in [24]. It is much stronger than just Liouville integrability. In this paper we discuss both Liouville and superintegrable examples, particularly in the context of satisfying the weak energy condition.

Conformally flat spaces in 2 dimensions
The Lagrangian kinetic energy is directly related to the metric ds 2 = n i, j=1 g i j dq i dq j . In 2 dimensions, all metrics are conformally flat, so there exist coordinates (q 1 , q 2 ) = (x, y), with respect to which the metric takes the form with only one independent curvature component, given by the scalar curvature where f = f x x + f yy . The Eisenhart lift of this metric incorporates the potential function U (x, y) in the definition of L: We require both ϕ(x, y) and U (x, y) to be strictly positive, at least in some region of the (x, y) plane which is invariant under the geodesic flow. We also require that the signature of the metric is (+, +, +, −).
Since the 2D metric satisfies the Einstein vacuum equations, the Einstein tensor for the Eisenhart lift has at most only 3 non-zero components, which are easily computed: with R taking the same value (4b). These can be equated with the energy momentum tensor, which naturally lies in the tensor spaces

The Hamiltonian formulation of the Eisenhart lift
The Hamiltonian in 2 dimensions takes the form The Eisenhart lift of this Hamiltonian is In the 2-dimensional context the symbol "t" represents the ∂ p x , etc, whilst for the Eisenhart lift, t is a coordinate in spacetime, whose evolution in the new "time parameter" τ can be found by using Hamilton's equations from (5b): Equations (6c) immediately give that p v and p t are first integrals, so we can set where κ and E are constants. The evolution of (x, y, p x , p y ) just gives the original equations associated with the 2dimensional Hamiltonian (5a). If, in 2 dimensions, there is a quadratic first integral of (5a): then, in 4 dimensions, this is "homogenised" to give and { p v , F} = {p t , F} = 0, so H, F, p v , p t are in involution. Thus Liouville integrability is preserved in this construction (see also [8]).
In a similar way we can extend higher order integrals from 2 dimensions to the 4-dimensional Eisenhart lift.

A conformal invariant in 4 D with an additional quadratic invariant
Whilst the 2-dimensional metric is conformally flat (with an infinite number of conformal Killing vectors), the 4dimensional extension (4c) is not (unless the 2D curvature R = 0 and U (x, y) takes a very simple form). For the Hamiltonian (5b) we ask that there exists a first degree (in momenta) conformal invariant where a 1 , a 2 , a 3 are functions of (x, y, v, t), whilst a 4 (t) is a function of only t. This equation is homogeneously quadratic in momenta, so gives us 10 equations for the coefficients a i , which can be partially solved. We quickly find that a 1 and a 2 can be written in terms of a potential w(x, y, v, t), with a 1 = w x , a 2 = −w y and with w satisfying the 2-dimensional Laplace equation The equations are generally fairly complicated, but simplified by assuming that w is independent of v and t.
First the coefficients of p x p v and p y p v give us formulae for w x x x and w xxy , which are compatible ((w x x ) xy = (w x x ) yx is identically satisfied) and can be integrated to give after which we have explicit formulae for a 3 and a 4 , giving 2 and U (x, y) satisfying Since we have {t, H } = −2 p v and (9d), it follows that is a first integral. This is in addition to p v , p t and any first integrals inherited from the 2-dimensional system. These additional first integrals satisfy 3.1 Some solutions of the system (9) We can read Eq. (9) in two ways. We can choose a solution of Laplace's equation (9a) to insert into the remaining equations, to be solved for ϕ(x, y) and U (x, y), and then analyse the resulting system. This analysis could be of the resulting Hamiltonian system and/or the geometry of the resulting 4-dimensional metric. Is it an Einstein vacuum metric? If not, what is its energy-momentum tensor?
Alternatively, we could start with a 2-dimensional metric (i.e., the function ϕ(x, y)) and then solve Eq. (9) for w(x, y), giving the conformal symmetry K (and the value of b), and then solve (9e) for the compatible family of potential functions U (x, y).
With this choice of w(x, y), we easily find that and K takes the form K = 2(x p x + yp y ) + (c 3 t + (2b − c 4 )v) p v + c 4 t p t . Thus, with this particular w(x, y), we have two arbitrary functions of a single variable in the definition of the metric (4c), each with this conformal symmetry K . We look at a particular example in Sect. 5.3.

Starting with ϕ(x, y) = 1
With this choice of ϕ(x, y), we easily find that with U (x, y) satisfying A simple form of solution corresponds to c 10 = c 01 = c 11 = This is the general Hénon-Heiles potential (with the harmonic terms suppressed), which is not integrable for arbitrary α, β. However, the choice α = 1, β = 2 corresponds to an integrable case (associated with the KdV hierarchy [25]), having a second quadratic integral, with Eisenhart Lift The full Poisson algebra of integrals With H, p v , p t , F in involution, the system has retained its Liouville integrability, and has gained one extra integral.

Starting with ϕ(x, y) = x 2
With this choice of ϕ(x, y), the 2D metric has constant curvature. We easily find that w(x, y) = c 01 y + This is an integrable potential (in 2 dimensions), having a second quadratic integral, with Eisenhart Lift Since b = 0, the above K is a first integral. The full Poisson algebra of integrals p v , With H, p v , p t , F in involution, the system has retained its Liouville integrability, and has gained one extra integral. This metric belongs to the classification by Darboux and Koenigs of metrics (see [17]) with only one Killing vector and two second order Killing tensors (necessarily not constant curvature). We easily find that A simple form of solution corresponds to c 01 = c 3 = 0, This is a superintegrable potential (in 2 dimensions) (see [17]), having two quadratic integrals, with Eisenhart Lifts The Poisson relations of these integrals are with the additional relations With H, p v , p t , F 1 in involution, the system has retained its Liouville integrability. In 4 degrees of freedom, we can have at most 7 functionally independent integrals. We have 7 integrals, but their Jacobian has rank 6 as a result of the algebraic relation which is a deformation of a similar relation in the 2dimensional context.

Energy-momentum tensor and Einstein equations
We saw in Eq. (6c) that the Hamiltonian (5b) always has first integrals p v and p t , corresponding to a pair of commuting Killing vectors where ∂ μ = ∂ ∂z μ and z μ = (x, y, v, t), the first of which is time-like while the second is null and covariantly constant, which implies (4c) belongs to the class of Brinkmann or Ppwave spacetimes.
The Eisenhart metric (4c) does not solve the vacuum Einstein equations unless U and log(ϕ) are harmonic functions. In this case the original 2D metric (4a) is actually flat (see (4b)). We saw that the general form of the Einstein tensor of the metric (4c) is given by (4d). With Einstein's equations the formula (4d) gives the form of T μν to be where and are given by with (11a) implying that Equation (11c) can be regarded in two different ways. Given the pair ( , ), (11c) provide the partial differential equations to fix the metric (4c) in a way compatible with the Einstein equations (11a). Vice versa, assuming the Lorentzian metric (4c) is given, then (11c) algebraically determine the matter characteristics ( , ), which fix the energy-momentum tensor (11b).
Let us discuss physical conditions to be imposed on the energy-momentum tensor (11b). The energy density measured by an observer moving along a time-like geodesic where we have used the form of the matrix (4d), with Equations (6) Since for time-like geodesics E > 0, 3 it suffices to require It turns out that (12c) also ensures the strong energy condition T μν − 1 2 g μν T λ λ dz μ dτ dz ν dτ ≥ 0, which reduces to The latter is automatically satisfied, as for time-like geodesics E > 1 2κ 2 . Recall that for a geodesic congruence which has a vanishing rotation tensor (hypersurface orthogonal), the strong energy condition implies that the geometry exerts a focussing effect on time-like geodesics. This is a consequence of the Raychaudhuri equation.
We conclude this section with a remark on the physical meaning of and . Raising the indices in (11b), one finds two non-vanishing components In order to correctly interpret and , it proves instructive to make recourse to the flat space which occurs at U = 0 and ϕ = 1 One sees that t and v are actually the double null coordinates. Implementing the coordinate transformation one brings (14) to the standard Minkowski metric, while the energy-momentum tensor acquires the form (T μν ∂ μ ∂ ν = Tμν∂μ∂ν)

Hidden symmetries and integrable models
We have already seen in Sect. 2.1 how each first integral of homogeneous degree m corresponds to a Killing tensor of 3 It follows from g μν dz μ dτ dz ν dτ = −1 and Eqs. (5a), (7a), and (7b) that E = H 2d + 1 2κ 2 > 0. Note that E can be arbitrarily large.
rank m of the metric corresponding to the kinetic energy. This can then be homogenised to sit within the Eisenhart lift, as seen in Sect. 2.3. In this section we consider some further examples of integrable and superintegrable systems and investigate the physical properties of the corresponding energy-momentum tensors.

The Darboux-Koenigs metrics
The Darboux-Koenigs metrics, classified in 1898 by Koenigs (see [16]), possess exactly one Killing vector (hence are not constant curvature) and a pair of second order Killing tensors (one of which is necessarily functionally dependent). There are 4 such metrics, characterised by the Killing vector K = ∂ y and the function ϕ(x, y) of (5a): where a > 1. The existence of the Killing vector and tensors means that the geodesic equations are superintegrable.
Computing the scalar curvature (11c) 3 , 4(a + cos x) 3 , one concludes that the Eisenhart uplifts of the Darboux-Koenigs metrics of types II, III, and IV violate the weak energy condition (12c), while type-I models are viable provided the domain x > 0 is chosen. In [17,18] there is a classification of potential functions U (x, y) for which superintegrability is maintained. In general, U (x, y) breaks the Killing isometry associated with K = ∂ y and makes two second rank Killing tensors functionally independent. Focusing on type-I models, one reveals two options [17]. The first is described by the (separable) Hamiltonian where b i ≥ 0 are constants, and two quadratic integrals of motion which obey the non-linear algebra jointly with F 3 = {F 1 , F 2 }: Demanding to be non-negative one is led to set b 1 = 0. The region x > 0 is invariant under the Hamiltonian flow: supposing x > 0 at t = 0, then H = E > 0 at t = 0 and therefore, for all t ≥ 0. Hence, noting that all terms within the parentheses of (18) are positive definite, we have x > 0 for all t ≥ 0.
Another type-I model, which is compatible with the weak energy conditions (12c), is governed by the (separable) Hamiltonian where b 1 , b 2 ≥ 0 are constants, and integrals of motion Along with F 3 = {F 1 , F 2 }, they form the Poisson algebra The corresponding is positive-definite Using the same argument as before, the domain x > 0 is invariant under the Hamiltonian flow, so x > 0 for all t ≥ 0. We will have more to say about separable systems in Sect. 5.3.1.

A superintegrable system on the two-dimensional sphere
Consider the 2-dimensional superintegrable systems on a constant curvature space, admitting two quadratic constants of the motion. To have positive we focus on models on the 2-dimensional sphere (classified in [26]), which are invariant under the rotation group, whose infinitesimal generators form a simple Poisson algebra: satisfying {L i , L j } = i jk L k , with Casimir L 2 1 + L 2 2 + L 3 3 = p 2 θ + sin −2 θ p 2 φ . As an example, consider the system (S9) in [26], written in spherical coordinates where (θ, p θ ) and (φ, p φ ) are canonical pairs and α, β, γ are free parameters (coupling constants). The model is characterized by two quadratic integrals of motion [26,27] which along with H form a functionally independent set. We define F 3 = 1 4 {F 1 , F 2 } = L 1 L 2 L 3 + "first order terms", where the leading order term is determined by the Poisson relations of L i . This is not functionally independent of H, F 1 , F 2 and satisfies the polynomial constraint (21c) below The leading order term is easily found by considering L 2 1 L 2 2 L 2 3 . We can use this constraint to determine the final two relations of our Poisson algebra: Whilst (21a) is not explicitly written in the conformally flat form of (5a), we can clearly replace the original kinetic energy by that on a 2-dimensional sphere, with metric dθ 2 + sin 2 θ dφ 2 , replace U (x, y) by U (θ, φ), and impose the Einstein equations to find = 2 and = α 2 (5 + cos 2φ) 2 sin 4 θ sin 4 φ + β 2 (5 − cos 2φ) 2 sin 4 θ cos 4 φ + γ 2 (2 − cos 2θ) cos 4 θ + 2(α 2 + β 2 + (α 2 − β 2 ) cos 2φ) cos 2θ sin 2 2φ sin 4 θ .
Combining the terms involving α 2 , β 2 and γ 2 separately, one can verify that is positive definite. Killing tensors associated with constants of the motion (21b) arise in the usual way.

Some Liouville integrable systems in cartesian coordinates
Superintegrable systems have very rigid choices of both ϕ(x, y) and U (x, y), so leave no room for manoeuvering and to be non-negative. On the other hand, Liouville integrable (including separable) systems can have potentials with arbitrary functions, which can be judiciously chosen as in the examples below.

A simple separable system
Consider the Poisson commuting pair where ϕ(x) and u(y) are positive definite functions. The Eisenhart metric (4c) associated with (22) admits the second rank Killing tensor corresponding to F (of (7d)), whose upper index form has non-zero components f yy = 1, f vv = 4u(y), which can be lowered with (4c) to give non-vanishing components If we choose ϕ(x) = x −2 and u(y) = y 2 , then from (11c) we find which clearly obey the weak energy conditions (12c). With these choices of ϕ(x) and U (x, y) = 1 2 ϕ(x)u(y), Eq. (9) give the conformal symmetry and the additional quadratic integral F K (see (10a)), also defining a second rank Killing tensor, together with the following (non-zero) Poisson relations:

The flat metric with quartic potential
Consider the Poisson commuting pair which is a well known integrable system, separable in parabolic coordinates. In fact, there are infinitely many homogeneous polynomial potentials, including the Hénon-Heiles potential of Sect. 3.1.2, which are separable in parabolic coordinates. The even degree polynomials give an energy-momentum tensor (of the Eisenhart lift) satisfying the weak energy condition. The above quartic potential gives (using (11c)): = 36(6x 2 + y 2 ), = 0, which clearly obey the weak energy conditions (12c). With these choices of ϕ(x) and U (x, y), equations (9) give the conformal symmetry and the additional quadratic integral F K (see (10a)), with both F and F K defining second rank Killing tensors. These functions satisfy the following (non-zero) Poisson relations:

Specific case from Sect. 3.1.1
In Sect. 3.1.1, we specified w(x, y) = x 2 − y 2 (in the definition of the conformal symmetry K ) and derived the corresponding general form of ϕ(x, y) and U (x, y) in terms of 2 arbitrary functions of a single variable: A y x and B y x . Equation (11c) give and as differential expressions in A and B. We don't propose here to analyse these expressions to determine all specific forms of A and B for which the weak energy conditions (12c) is satisfied, but clearly there would be many viable cases. If we choose A y x ≡ 1, b = 2, then we have the flat metric with ϕ(x, y) = 1, so = R = 0. If we now specify the two parameters c 3 = 0, c 4 = −2, then so just need to choose the form of B y x so that this expression is non-negative for all (x, y). For example U (x, y) = d 1 x 4 + d 2 x 2 y 2 + d 3 y 4 , = 2(6d 1 + d 2 )x 2 + 2(d 2 + 6d 3 )y 2 , so is positive definite whenever the two coefficients are positive. In all cases, this has conformal symmetry K = 2(x p x + yp y + 3vp v − t p t ) with {K , H } = 4H, regardless of whether or not we choose an integrable case, such as (d 1 , d 2 , d 3 ) = (16, 12, 1), as above.

Conclusion
In this paper we have considered the Eisenhart lift of a fairly general 2-dimensional metric, with particular interest in comparing the Hamiltonian properties of the 2D system and its 4D lift. Whilst the original 2D metric (4a) was conformally flat, its Lorentzian counterpart (4c) fails to be so. Computing the Weyl tensor, one can verify that its vanishing requires log(ϕ) = 0, which in turn implies that the metric (4a) is flat. This makes the cosmological applications in the spirit of a recent work [28] problematic. Furthermore, it is worth recalling that reductions of the Goryachev-Chaplygin and Kovalevskaya tops, which are obtained by discarding a cyclic variable, result in 2D integrable systems in curved space possessing cubic and quartic integrals of motion, respectively. The corresponding Eisenhart metrics and higher rank Killing tensors were constructed in [8]. One can verify that, while > 0, fails to be positive definite in the whole domain thus ruling out these examples from the physically acceptable list.
In our derivation of the conformal symmetry (9c), we made several restrictions on the coefficients a i (x, y, v, t), so an obvious question is whether a more general solution can be found, or even some algebra of independent conformal symmetries can be found. Even with these restrictions, we found that a given function w(x, y) led to ϕ(x, y) and U (x, y), depending upon arbitrary functions. Particular choices can lead to energy-momentum tensors obeying the weak energy condition. There is therefore the question of how to choose these functions in a more systematic way and, if possible, to classify these.
It would be interesting to generalise the analysis in this work to d > 2 mechanics on conformally flat backgrounds. The first problem is that the Einstein tensor would not have such a simple decomposition as (4d), so the energymomentum tensor would be more difficult to analyse. Secondly, the calculation of the conformal symmetry, as in Sect. 3, would be considerably more complex.
Another open problem is whether Lorentzian metrics admitting third or higher rank Killing tensors, linked to the work in [20,21] give rise to physically meaningful solutions.