Eisenhart Lift of $2$--Dimensional Mechanics

TL;DR

This paper explores the Eisenhart lift of 2D mechanics on curved backgrounds, deriving symmetries, constructing solutions to Einstein equations, and analyzing superintegrable models with hidden symmetries.

Contribution

It provides a detailed study of the Eisenhart lift for 2D systems, including conformal symmetries, energy-momentum tensors, and conditions for physically viable metrics.

Findings

Derived conformal symmetry and quadratic integral.

Constructed energy--momentum tensor satisfying Einstein equations.

Identified conditions under which Darboux--Koenigs metrics meet energy conditions.

Abstract

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.

Equations

$$\frac{d^2{q}^k}{d{t}^2}+\sum _{i,j=1}^n{\Gamma }_{ij}^k\frac{d{q}^i}{dt}\frac{d{q}^j}{dt}=-\sum _{i=1}^n{g}^{ki}{\partial }_{i}U,k=1,\dots ,n,$$

$$H_0=\frac{1}{2}\sum _{i,j=1}^n{g}^{ij}{p}_i{p}_j\text{mbox}andK=\sum _{i=1}^n{a}^i(\mathbf{q})p_i,$$

$$\{K,H_0\}=0\Rightarrow {\nabla }_{(i}{a}_{j)}=0,$$

$$\hat{K}=\sum _{i=1}^n{a}^i(\mathbf{q})\frac{\partial }{\partial q^i},$$

$$\{F,H_0\}=0\Rightarrow {\nabla }_{(i}{f}_{jk)}=0,$$

$$\frac{d{q}^i}{ds}=a^i(\mathbf{q}),$$

$$d{s}^2=\frac{d{x}^2+d{y}^2}{\phi (x,y)},$$

$$R=\phi \Delta \log (\phi ),$$

$$d{s}^2=\frac{d{x}^2+d{y}^2}{\phi (x,y)}-dvdt-2U(x,y)d{t}^2.$$

$$R_{ij}-\frac{1}{2}Rg_{ij}=\left(\begin{array}[]{cccc}0&0&0&0\\ 0&0&0&0\\ 0&0&0&\frac{1}{4}R\\ 0&0&\frac{1}{4}R&\varphi\Delta U+RU\end{array}\right),$$

$$H^{2d}=\frac{1}{2}\phi (x,y)(p_x^2+p_y^2)+U(x,y).$$

$$H=\frac{1}{2}\phi (x,y)(p_x^2+p_y^2)+4U(x,y)p_v^2-2p_v{p}_t.$$

$$\frac{dx}{d\tau }=\phi p_x,\frac{dy}{d\tau }=\phi p_y,\frac{dv}{d\tau }=8U{p}_v-2p_t,\frac{dt}{d\tau }=-2p_v,$$

$$\frac{d{p}_x}{d\tau }=-\frac{1}{2}{\phi }_x(p_x^2+p_y^2)-4U_x{p}_v^2,\frac{d{p}_y}{d\tau }=-\frac{1}{2}{\phi }_y(p_x^2+p_y^2)-4U_y{p}_v^2,$$

$$\frac{d{p}_v}{d\tau }=0,\frac{d{p}_t}{d\tau }=0.$$

$$p_v=-\frac{1}{2}\kappa ,$$

$$p_t=-\kappa \mathcal{E}$$

$$F^{2d}=\sum _{i,j=1}^2{f}^{ij}(x,y)p_i{p}_j+w(x,y),\text{mbox}with\{H^{2d},F^{2d}\}=0,$$

$$F=\sum _{i,j=1}^2{f}^{ij}(x,y)p_i{p}_j+4w(x,y)p_v^2,\text{mbox}satisfying\{H,F\}=0,$$

$$K=a_1{p}_x+a_2{p}_y+a_3{p}_v+a_4{p}_t,\text{mbox}satisfying\{K,H\}=\sigma (x,y,v,t)H,$$

$$w_{xx}+w_{yy}=0.$$

$$w_{xx}=\frac{1}{2}(w_x{\partial }_x-w_y{\partial }_y)\log (\phi )+b=\frac{1}{2}\{\log (\phi ),K\}+b,$$

$$K=w_x{p}_x-w_y{p}_y+((2b-c_4)v+c_{3}t)p_v+c_{4}t{p}_t,$$

$$\{K,H\}=2bH,$$

$$\{U,K\}=(w_x{\partial }_x-w_y{\partial }_y)U=2(b-c_4)U-\frac{1}{2}{c}_3.$$

$$F_K=K{p}_v+btH$$

$$\{p_v,p_t\}=0,\{p_v,F_K\}=(c_4-2b)p_v^2,\{p_t,F_K\}=-(c_3{p}_v+c_4{p}_t)p_v-bH.$$

$$\phi (x,y)=x^{2-b}A\left(\frac{y}{x}\right),U(x,y)=\frac{c_3}{4(b-c_4)}+x^{b-c\text{-}4}B\left(\frac{y}{x}\right),$$

$$w(x,y)$$

$$K$$

$$(c_{10}+bx+c_{11}y)U_x-(c_{01}+c_{11}x-by)U_y=2(b-c_4)U-\frac{1}{2}{c}_3.$$

$$U(x,y)=\alpha x{y}^2+\beta x^3,K=\frac{1}{2}b(2x{p}_x+2y{p}_y+5v{p}_v-t{p}_t).$$

$$F=p_y(y{p}_x-x{p}_y)+y^2(4x^2+y^2)p_v^2.$$

$$\{p_v,p_t\}=\{p_v,F\}=\{p_t,F\}=0,\{p_v,F_K\}=-\frac{5b}{2}{p}_v^2,$$

$$\{p_t,F_K\}=\frac{1}{2}b(p_v{p}_t-2H),\{F,F_K\}=-b{p}_{v}F.$$

$$w(x,y)=c_{01}y+\frac{1}{2}{c}_{10}(x^2-y^2)+\frac{1}{6}{c}_{11}(3x^{2}y-y^3),b=0.$$

Full text

Keywords: Eisenhart lift, conformally flat metrics, Darboux–Koenigs metrics, Hamiltonian system, superintegrability.

MSC: 17B63,37J35,70G45,70H06.

PACS: 02.40.Ky, 04.20.Jb, 11.10.Ef, 11.30.-j, 11.30.Na.

1 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,\dots,d$, which is governed by the Lagrangian $\mathcal{L}$, can be embedded into the geodesic equations of the Brinkmann–-type metric $2\mathcal{L}dt^{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 $\tau$, the equations of motion for $q^{i}(t)$ coincide with those following from the Lagrangian $\mathcal{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$.

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.111Ricci–flat metrics of ultrahyperbolic signature, which admit higher rank Killing tensors, were constructed in [6]. 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 [7]–[12]. 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 [7, 8]). 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 [13] 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 $\frac{1}{2}n(n+1)$. In this case, the kinetic energy is actually the second order Casimir function of the symmetry algebra (see [14]). 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 [15], which are described and analysed in [16, 17]. There are other examples of conformally flat spaces (but not constant curvature), possessing one Noether constant and a cubic integral (classified in [18] and further studied and generalised in [19, 20]), 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 [21, 22] 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 $4$–dimensional 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 Section 2, the equations of motion of a generic $2$–dimensional mechanical system in curved space are embedded into the geodesic equations of a Brinkmann type $4$–dimensional metric. The latter is determined by two scalar functions which are just the conformal factor and the potential 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 Section 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 Section 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 Section 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 Section 6.

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

2.1 First Integrals in Classical Mechanics

A particle moving in a curved background, with metric coefficients $g_{ij}$ and scalar potential $U({\bf q})$, can be written in Lagrange form, with ${\mathcal{L}}=\frac{1}{2}\sum_{i,j=1}^{n}g_{ij}\dot{q}^{i}\dot{q}^{j}-U({\bf q})$ or, equivalently, in Hamiltonian form with $H=\frac{1}{2}\sum_{i,j=1}^{n}g^{ij}p_{i}p_{j}+U({\bf q})$, where $g^{ij}$ are inverse metric coefficients and $p_{i}=\sum_{j=1}^{n}g_{ij}\dot{q}^{j}$. In either case, we obtain

[TABLE]

where $\Gamma^{k}_{ij}$ 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

[TABLE]

since the first $n$ components are written entirely in terms of the position variables. $\hat{K}$ is just the Killing vector corresponding to the Noether constant $K$.

Similarly, if $F=\sum_{i,j=1}^{n}f^{ij}p_{i}p_{j}$, $f^{ij}$ the components of a symmetric matrix, then

[TABLE]

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

[TABLE]

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=\frac{{\partial}}{{\partial}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},\dots,Q_{n}$, such that $K_{i}=\frac{{\partial}}{{\partial}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 [23]. 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.

2.2 Conformally Flat Spaces in $2$ Dimensions

The Lagrangian kinetic energy is directly related to the metric $ds^{2}=\sum_{i,j=1}^{n}g_{ij}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

[TABLE]

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:

[TABLE]

with $R$ taking the save value (4b). These can be equated with the energy momentum tensor, which naturally lies in the tensor spaces ${\partial}_{t}\otimes{\partial}_{v}+{\partial}_{v}\otimes{\partial}_{t}$ and ${\partial}_{t}\otimes{\partial}_{t}$.

2.3 The Hamiltonian Formulation of the Eisenhart Lift

The Hamiltonian in 2 dimensions takes the form

[TABLE]

In the $2$–dimensional context the symbol ”$t$” represents the Hamiltonian ”time-parameter”, with $\frac{dx}{dt}=\frac{{\partial}H^{2d}}{{\partial}p_{x}}$, etc, whilst for the Eisenhart lift, $t$ is a coordinate in spacetime, whose evolution in the new ”time parameter” $\tau$ can be found by using Hamilton’s equations from (5b):

[TABLE]

Equations (6c) immediately give that $p_{v}$ and $p_{t}$ are first integrals, so we can set

[TABLE]

If, in 2 dimensions, there is a quadratic first integral of (5a):

[TABLE]

then, in 4 dimensions, this is ”homogenised” to give

[TABLE]

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

In a similar way we can extend higher order integrals from 2 dimensions to the $4$–dimensional Eisenhart lift.

3 A Conformal Invariant in $4D$ with an Additional Quadratic Invariant

Whilst the $2$–dimensional metric is conformally flat (with an infinite number of conformal Killing vectors), the $4$–dimensional 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

[TABLE]

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

[TABLE]

First the coefficients of $p_{x}p_{v}$ and $p_{y}p_{v}$ give us formulae for $w_{xxx}$ and $w_{xxy}$, which are compatible ($(w_{xx})_{xy}=(w_{xx})_{yx}$ is identically satisfied) and can be integrated to give

[TABLE]

after which we have explicit formulae for $a_{3}$ and $a_{4}$, giving222Here and below $b$ and $c_{i}$ designate constants.

[TABLE]

with

[TABLE]

and $U(x,y)$ satisfying

[TABLE]

Since we have $\{t,H\}=-2p_{v}$ and (9d), it follows that

[TABLE]

3.1 Some Solutions of the System (9)

We can read Equations (9) in two ways.

We can choose a solution of Laplace’s equation (9a) to insert into the remaining equations, to be solved for $\varphi(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 (ie the function $\varphi(x,y)$) and then solve Equations (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)$.

3.1.1 Starting with $w(x,y)=x^{2}-y^{2}$

With this choice of $w(x,y)$, we easily find that

[TABLE]

and $K$ takes the form $K=2(xp_{x}+yp_{y})+(c_{3}t+(2b-c_{4})v)p_{v}+c_{4}tp_{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 Section 5.3.

3.1.2 Starting with $\varphi(x,y)=1$

With this choice of $\varphi(x,y)$, we easily find that

[TABLE]

with $U(x,y)$ satisfying

[TABLE]

A simple form of solution corresponds to $c_{10}=c_{01}=c_{11}=c_{3}=0$, giving $U(x,y)=x^{2\left(1-\frac{c_{4}}{b}\right)}\psi\left(\frac{y}{x}\right)$. If we choose $c_{4}=-\frac{1}{2}b$ and $\psi(z)=\alpha z^{2}+\beta$, we obtain

[TABLE]

This is the general Hénon-Heiles potential (with the harmonic terms suppressed), which is not integrable for arbitrary $\alpha,\,\beta$. However, the choice $\alpha=1,\,\beta=2$ corresponds to an integrable case (associated with the KdV hierarchy [24]), having a second quadratic integral, with Eisenhart Lift

[TABLE]

The full Poisson algebra of integrals $p_{v},p_{t},F,F_{K}$ is

[TABLE]

With $H,p_{v},p_{t},F$ in involution, the system has retained its Liouville integrability, and has gained one extra integral.

3.1.3 Starting with $\varphi(x,y)=x^{2}$

With this choice of $\varphi(x,y)$, the 2D metric has constant curvature. We easily find that

[TABLE]

A simple form of solution corresponds to $c_{01}=c_{11}=c_{3}=0$, giving $U(x,y)=x^{-\frac{2c_{4}}{c_{10}}}\psi\left(\frac{y}{x}\right)$. If we choose $c_{10}=-\frac{1}{2}c_{4}$ and $\psi(z)=1+4z^{2}$, we obtain

[TABLE]

This is an integrable potential (in 2 dimensions), having a second quadratic integral, with Eisenhart Lift

[TABLE]

Since $b=0$, the above $K$ is a first integral. The full Poisson algebra of integrals $p_{v},p_{t},K,F$ is

[TABLE]

With $H,p_{v},p_{t},F$ in involution, the system has retained its Liouville integrability, and has gained one extra integral.

3.1.4 The Darboux–Koenigs Metric $D_{1}$ with $\varphi(x,y)=x^{-1}$

This metric belongs to the classification by Darboux and Koenigs of metrics (see [16]) with only one Killing vector and two second order Killing tensors (necessarily not constant curvature). We easily find that

[TABLE]

A simple form of solution corresponds to $c_{01}=c_{3}=0,c_{4}=2b$, giving $U(x,y)=x^{-3}\psi\left(\frac{y}{x}\right)$. If we choose $\psi(z)=\mu z^{-2}$, we obtain

[TABLE]

This is a superintegrable potential (in 2 dimensions) (see [16]), having two quadratic integrals, with Eisenhart Lifts

[TABLE]

together with

[TABLE]

The Poisson relations of these integrals are

[TABLE]

with the additional relations

[TABLE]

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

[TABLE]

which is a deformation of a similar relation in the $2$–dimensional context.

4 Energy–Momentum Tensor and Einstein Equations

We saw in Equation (6c) that the Hamiltonian (5b) always has first integrals $p_{v}$ and $p_{t}$, corresponding to a pair of commuting Killing vectors

[TABLE]

where $\partial_{\mu}=\frac{\partial}{\partial z^{\mu}}$, 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 Pp-wave spacetimes.

The Eisenhart metric (4c) does not solve the vacuum Einstein equations unless $U$ and $\log(\varphi)$ 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

[TABLE]

Equations (11c) can be regarded in two different ways. Given the pair $(\Omega,\Sigma)$, (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 $(\Omega,\Sigma)$, 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 $z^{\mu}=z^{\mu}(\tau)$ is

[TABLE]

The latter is automatically satisfied, as for time–like geodesics $\mathcal{E}>\frac{1}{2\kappa^{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 $\Omega$ and $\Sigma$. Raising the indices in (11b), one finds two non–vanishing components

[TABLE]

In order to correctly interpret $\Omega$ and $\Sigma$, it proves instructive to make recourse to the flat space which occurs at $U=0$ and $\varphi=1$

[TABLE]

One sees that $t$ and $v$ are actually the double null coordinates. Implementing the coordinate transformation

[TABLE]

one brings (14) to the standard Minkowski metric, while the energy–momentum tensor acquires the form ($T^{\mu\nu}\partial_{\mu}\partial_{\nu}={\tilde{T}}^{\tilde{\mu}\tilde{\nu}}\partial_{\tilde{\mu}}\partial_{\tilde{\nu}}$)

[TABLE]

At this point $\frac{1}{8\pi}\left(\Omega+\frac{1}{2}\Sigma\right)$ can be identified with the energy density, with $-\frac{1}{8\pi}\Omega$ being the energy flux density in the direction orthogonal to the $(x,y)$–plane, while $\frac{1}{8\pi}\left(\Omega-\frac{1}{2}\Sigma\right)$ is the only non–vanishing component of the stress tensor.

5 Hidden Symmetries and Integrable Models

We have already seen in Section 2.1 how each first integral of homogeneous degree $m$ corresponds to a Killing tensor of rank $m$ of the metric corresponding to the kinetic energy. This can then be homogenised to sit within the Eisenhart lift, as seen in Section 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.

5.1 The Darboux–Koenigs Metrics

The Darboux–Koenigs metrics, classified in 1898 by Koenigs (see [15]), 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={\partial}_{y}$ and the function $\varphi(x,y)$ of (5a):

[TABLE]

where $a>1$. The existence of the Killing vector and tensors means that the geodesic equations are superintegrable.

Computing the scalar curvature (11c)

[TABLE]

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 [16, 17] 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={\partial}_{y}$ and makes two second rank Killing tensors functionally independent. Focusing on type–I models, one reveals two options [16]. The first is described by the (separable) Hamiltonian

[TABLE]

where $b_{i}\geq 0$ are constants, and two quadratic integrals of motion

[TABLE]

which obey the non–linear algebra jointly with $F_{3}=\{F_{1},F_{2}\}$:

[TABLE]

Demanding $\Omega$ to be non–negative

[TABLE]

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\geq 0$. Hence, noting that all terms within the parentheses of (18) are positive definite, we have $x>0$ for all $t\geq 0$.

Another type–I model, which is compatible with the weak energy conditions (12c), is governed by the (separable) Hamiltonian

[TABLE]

where $b_{1},b_{2}\geq 0$ are constants, and integrals of motion

[TABLE]

Along with $F_{3}=\{F_{1},F_{2}\}$, they form the Poisson algebra

[TABLE]

The corresponding $\Omega$ is positive–definite

[TABLE]

Using the same argument as before, the domain $x>0$ is invariant under the Hamiltonian flow, so $x>0$ for all $t\geq 0$.

We will have more to say about separable systems in Section 5.3.1.

5.2 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 $\Sigma$ we focus on models on the two–dimensional sphere (classified in [25]), which are invariant under the rotation group, whose infinitesimal generators form a simple Poisson algebra:

[TABLE]

satisfying $\{L_{i},L_{j}\}=\epsilon_{ijk}L_{k}$, with Casimir $L_{1}^{2}+L_{2}^{2}+L_{3}^{3}=p_{\theta}^{2}+\sin^{-2}\theta p_{\phi}^{2}$.

As an example, consider the system $(S9)$ in [25], written in spherical coordinates

[TABLE]

Whilst (LABEL:sp) is not explicitly written in the conformally flat form of (5a), we can clearly replace the original kinetic energy by that on a two–dimensional sphere, with metric $d\theta^{2}+\sin^{2}{\theta}d\phi^{2}$, replace $U(x,y)$ by $U(\theta,\phi)$, and impose the Einstein equations to find $\Sigma=2$ and

[TABLE]

Combining the terms involving $\alpha^{2}$, $\beta^{2}$ and $\gamma^{2}$ separately, one can verify that $\Omega$ is positive definite. Killing tensors associated with constants of the motion (21b) arise in the usual way.

5.3 Some Liouville Integrable Systems in Cartesian Coordinates

Superintegrable systems have very rigid choices of both $\varphi(x,y)$ and $U(x,y)$, so leave no room for manoeuvering $\Sigma$ and $\Omega$ 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.

5.3.1 A Simple Separable System

Consider the Poisson commuting pair

[TABLE]

where $\varphi(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

[TABLE]

If we choose $\varphi(x)=x^{-2}$ and $u(y)=y^{2}$, then from (11c) we find

[TABLE]

which clearly obey the weak energy conditions (12c).

With these choices of $\varphi(x)$ and $U(x,y)=\frac{1}{2}\varphi(x)u(y)$, equations (9) give the conformal symmetry

[TABLE]

and the additional quadratic integral $F_{K}$ (see (10a)), also defining a second rank Killing tensor, together with the following (non-zero) Poisson relations:

[TABLE]

5.3.2 The Flat Metric with Quartic Potential

Consider the Poisson commuting pair

[TABLE]

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 Section 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)):

[TABLE]

which clearly obey the weak energy conditions (12c).

With these choices of $\varphi(x)$ and $U(x,y)$, equations (9) give the conformal symmetry

[TABLE]

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:

[TABLE]

5.3.3 Specific Case from Section 3.1.1

In Section 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 $\varphi(x,y)$ and $U(x,y)$ in terms of 2 arbitrary functions of a single variable: $A\left(\frac{y}{x}\right)$ and $B\left(\frac{y}{x}\right)$. Equations (11c) give $\Omega$ and $\Sigma$ 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\left(\frac{y}{x}\right)\equiv 1,\,b=2$, then we have the flat metric with $\varphi(x,y)=1$, so $\Sigma=R=0$. If we now specify the two parameters $c_{3}=0,c_{4}=-2$, then

[TABLE]

so just need to choose the form of $B\left(\frac{y}{x}\right)$ so that this expression is non–negative for all $(x,y)$. For example

[TABLE]

so $\Omega$ is positive definite whenever the two coefficients are positive. In all cases, this has conformal symmetry

[TABLE]

regardless of whether or not we choose an integrable case, such as $(d_{1},d_{2},d_{3})=(16,12,1)$, as above.

6 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 $\Delta\log(\varphi)=0$, which in turn implies that the metric (4a) is flat. This makes the cosmological applications in the spirit of a recent work [27] 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 [7]. One can verify that, while $\Sigma>0$, $\Omega$ 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 $\varphi(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 energy–momentum tensor would be more difficult to analyse. Secondly, the calculation of the conformal symmetry, as in Section 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 [19, 20] give rise to physically meaningful solutions.

Acknowledgements

This work was supported by the Tomsk Polytechnic University competitiveness enhancement program.