## Geometric Linearization for Constraint Hamiltonian Systems

Andronikos Paliathanasis 1 2 3 , , ∗

1 Institute of Systems Science, Durban University of Technology

PO Box 1334, Durban 4000, Republic of South Africa

2 School for Data Science and Computational Thinking, Stellenbosch University

44 Banghoek Rd, Stellenbosch 7600, South Africa

- 3 Departamento de Matem´ aticas, Universidad Cat´lica del Norte, o Avda. Angamos 0610, Casilla 1280 Antofagasta, Chile

August 5, 2024

## Abstract

This study investigates the geometric linearization of constraint Hamiltonian systems using the Jacobi metric and the Eisenhart lift. We establish a connection between linearization and maximally symmetric spacetimes, focusing on the Noether symmetries admitted by the constraint Hamiltonian systems. Specifically, for systems derived from the singular Lagrangian

$$L \left ( N, q ^ { k }, \dot { q } ^ { k } \right ) = \frac { 1 } { 2 N } g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } - N V ( q ^ { k } ), \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \Delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta \delta$$

where N and q i are dependent variables and dim g ij = n , the existence of n n ( +1) 2 Noether symmetries is shown to be equivalent to the linearization of the equations of motion. The application of these results is demonstrated through various examples

∗ Email: anpaliat@phys.uoa.gr

of special interest. This approach opens new directions in the study of differential equation linearization.

Keywords: Constraint Hamiltonian systems; Lie symmetries; Linearization;

## 1 Introduction

The theory of invariant infinitesimal transformations, known as Lie symmetry analysis [1-3], is a powerful tool for the analytic treatment of nonlinear differential equations. The concept of symmetry is based on the existence of invariant functions in the solution space of the differential equations [4]. These invariant functions are mainly used to simplify differential equations and to derive solutions [5-7]. Lie symmetry analysis systematically treats dynamical systems and has been widely explored in various areas of applied mathematics, see for instance [8-19] and references therein..

Symmetry in physics is crucial because it is related to conservation laws [20-22]. The conservation laws of momentum for free particles, energy for conservative systems, Kepler's laws of motion, and many others are connected to infinitesimal transformations that keep the differential equations invariant [5]. The existence of a sufficiently large number of Lie symmetries for a given differential equation allows for solving the equation by repeated order reduction using similarity transformations or by determining a sufficient number of first integrals [6]. Additionally, Lie symmetries can be used to classify dynamical systems through their admitted invariant functions [7].

This latter characteristic is important for the concept of global linearization of differential equations via point transformations. Specifically, nonlinear differential equations can be linearized if they admit a specific number of Lie symmetries [4], as will be detailed in the following sections. In this study, we focus on the global linearization of a family of constraint Hamiltonian systems. We derive a set of geometric criteria that allow us, for the first time, to express nonlinear dynamical systems with fewer symmetries in a linear form.

We employ two different methods for the geometrization of non-relativistic dynamical systems with conservative forces, the Einsenhart-Duval lift and the Jacobi metric. In the Jacobi metric approach the conservative force is absorbed in the kinetic terms, such that the dynamical effects of the force to be part of the new geometry which describes the dynamical system. This is achieved by a reparametrization of the independent variables.

There are various applications of Jacobi metric for the study of dynamical systems [23-25], while the method has been generalized for the description of relativistic systems [26,27]. For a more general discussion we refer the reader to [28].

On the other hand, in the Eisenhart-Duval method the conservative forces are embedded into an extended geometry by introducing new independent variables and degrees of freedom. The method was introduced by Eisenhart in early of the twentieth century [29] and it was rediscovered latter by Duval et al. in [30]. The Eisenhart-Duval method is the nonrelativistic case of the Kaluza-Klein framework [31,32]. Due to the simplicity of the method, it has been widely applied for the study of nonlinear dynamical systems. The Eisenhart lift of two-dimensional mechanical systems with or without varying mass term is discussed in [33]. The superintegrability property of three-dimensional Newtonian systems investigated by using the Eisenhart lift in [34]. In [35] the Eisenhart lift applied for the study of the Toda chain for the study of nearest neighbor interacting particles on a line. The geodesic description of the two fixed centers problem by using the Eisenhart lift was studied in [36].

Nevertheless, the Eisenhart lift has been applied not only in classical mechanics but also for the analysis of the Schr¨dinger [37,38] and of the Dirac [39] equations of quantum o mechanics. There is a plethora of applications of the Eisenhart lift in relativistic physics and cosmology, see for instance [40,41,41,42] and references therein. The extended minisuperspace by the Eisenhart lift for the study of quantum cosmology introduced in [43]. Recently, in [44] the Eisenhart lift employed for the derivation of the analytic solutions for the field equations in scalar field cosmology. It was found that the field equations can be linearized in the framework of minimally coupled scalar field cosmological theory for a Friedmann-Lemaˆ ıtre-Robertson-Walker background geometry and for an exponential potential for the description of the mass for the scalar field. Moreover, in [45] it was found that the cosmological constant in Friedmann-Lemaˆ ıtre-Robertson-Walker geometry can be recovered by mean of the application of the hidden symmetries of for the extended minisuperspace in quantum cosmology.

Furthermore, [46] the Eisenhart lift employed to write the equation of motion for the oscillator in terms of the free particle. That is a well known result provided by the symmetry analysis which relates the two dynamical systems by sharing the common symmetry group [47]. We remark that by geometrize a given dynamical systems it is feasible to employ important results from the differential geometry in a systematic way.

This property is considered in the following in order to perform a geometric linearization for a class of constraint Hamiltonian systems. We shall see that there exist a relation between the geometric properties of the extended geometry in the Eisenhart lift with the admitted symmetries of the original system. The structure of the paper is as follows:

In Section 2, we discuss the mathematical tools primarily applied in this study. In particular we discuss the Lie symmetry analysis for differential equations and we focus in the case of second-order differential. We also present Noether's theorems, which play an important role in the main analysis of this study. Furthermore, in Section 3, we provide a literature review on the global linearization of differential equations through Lie symmetry analysis.

The family of dynamical systems considered in this work are introduced in Section 4. Specifically, we examine a family of dynamical systems described by a singular Lagrangian function, leading to equations of motion where the Hamiltonian function is constrained. Additionally, we discuss in detail two geometrization approaches for dynamical systems within this family. Section 5 includes the main results of this study, where we derive a new set of geometric conditions and criteria under which constraint Hamiltonian systems can be transformed into an equivalent system of a free particle in flat space, thereby becoming linearizable.

We demonstrate the application of this geometric approach in Section 6, where we explore the linearization of nonlinear dynamical systems of special interests in gravitational physics. The geometric linearization of the Szekeres system with or without the cosmological constant is discussed. Furthermore we investigate the linearization of the minisuperspace for the gravitational model in a static spherically symmetric spacetime with a charge, that is, the exact solution for the Reissner-Nordstr¨m black hole metric can o be constructed by the solution of the free particle. Last but not least, dynamical systems with interest in Newtonian mechanics are discussed. Finally, in Section 7, we summarize our results.

## 2 Preliminaries

In this Section, we briefly discuss the basic mathematical elements necessary for this study.

## 2.1 Lie symmetries of differential equations

In the following lines, we review the basic definitions for the Lie symmetry analysis of ordinary differential equations.

Let us assume the n -dimensional dynamical system of µ th-order differential equations of the form

$$q ^ { ( \mu ) i } = \omega ^ { i } \left ( t, q ^ { k }, \dot { q } ^ { k }, \ddot { q } ^ { k }, \dots, q ^ { ( \mu - 1 ) i } \right ) & & ( 1 ) \\ \i \text{dependent variables, and } q ^ { i } \, \text{are the dependent variables, } q ^ { i } = ( q ^ { 1 }, q ^ { 2 }, \dots, q ^ { N } ).$$

in which t is the independent variables, and q i are the dependent variables, q i = ( q , q 1 2 , ..., q N ) . Moreover, a dot represents the total derivative with respect to the independent variable t , i.e.,

$$\dot { q } ^ { i } = \frac { d q ^ { i } } { d t }, \, \ddot { q } ^ { i } = \frac { d ^ { 2 } q ^ { i } } { d t ^ { 2 } }, \, \dots, \, q ^ { ( \mu ) } = \frac { d ^ { \mu } q } { d t ^ { \mu } }.$$

Consider the infinitesimal transformation

$$\bar { t } \ = \ t + \varepsilon \xi \left ( t, q ^ { k } \right ) \\ \bar { q } ^ { i } \ = \ q ^ { i } + \varepsilon \eta ^ { i } \left ( t, q ^ { k } \right ) \\ \intertext { r f i d } \, X = \xi \left ( t, q ^ { k } \right ) \partial _ { t } + \eta ^ { i } \left ( t, q ^ { k } \right ) \partial _ { i }.$$

with generator the vector field X = ξ ( t, q k ) ∂ t + η i ( t, q k ) ∂ i . The vector field q [ µ ] is the µ th-prolongation of q in the Jet space J V = { t, q i , q ˙ i , q ¨ i , ..., q ( µ ) } defined as

$$X ^ { [ \mu ] } = X + \eta _ { [ 1 ] } ^ { i } \partial _ { \dot { q } ^ { i } } + \eta _ { [ 2 ] } ^ { i } \partial _ { \ddot { q } ^ { i } } + \dots + \eta _ { [ \mu ] } ^ { i } \partial _ { q ^ { ( \mu ) i } },$$

in which η i [1] , η i [2] , ..., η i [ µ ] are given by the following expressions

$$\eta _ { [ 1 ] } ^ { i } \ = \ \begin{array} { c c } \dot { \eta } ^ { i } - \dot { q } ^ { i } \dot { \xi }, \\ i \quad &. i \quad.. i \end{array} \quad & ( 6 )$$

$$\eta _ { [ 2 ] } ^ { i } \ = \ \dot { \eta } _ { [ 2 ] } ^ { i } - \ddot { q } ^ { i } \dot { \xi }, & & ( 7 ) \\ & & \theta \end{matrix}$$

...

(8)

$$\eta _ { [ \mu ] } ^ { i } \ = \ \dot { \eta } _ { [ \mu - 1 ] } ^ { i } - q ^ { ( \mu ) i } \dot { \xi }.$$

Then we say that the system of differential equations (1) will be invariant under the application of the infinitesimal transformations (3), (4) if and only if there exists a function λ such that the following condition holds [4]

where

$$[ X ^ { [ \mu ] }, A ] = \lambda A,$$

$$\left [ X ^ { [ \mu ] }, A \right ] & = X ^ { [ \mu ] } A - A X ^ { [ \mu ] } \, & ( 1 1 ) \\ & \quad 5$$

and operator A is Hamilton's vector

$$A = \frac { \partial } { \partial t } + \dot { q } ^ { i } \frac { \partial } { \partial q ^ { i } } + \ddot { q } ^ { i } \frac { \partial } { \partial \dot { q } ^ { i } } + \dots + \omega ^ { i } \left ( t, q ^ { k }, \dot { q } ^ { k }, \ddot { q } ^ { k }, \dots, q ^ { ( \mu - 1 ) i } \right ) \frac { \partial } { \partial q ^ { ( \mu ) i } }. \quad \ ( 1 2 ) \\ \text{Condition} \left ( 1 0 \right ) \text{is known as the Lie symmetry condition, and } X \text{ is a Lie point symmetry}$$

Condition (10) is known as the Lie symmetry condition, and X is a Lie point symmetry for the dynamical system (1).

Assume that f i is a solution vector for the dynamical system (1), that is Af i = 0, then, the Lie symmetry condition (10) becomes [ X [ µ ] , A ] f = λAf, that is, [4]

or equivalently

$$X \left ( A f ^ { i } \right ) = 0,$$

$$\eta _ { [ \mu ] } ^ { i } = X ^ { [ \mu - 1 ] } \omega ^ { i } \left ( t, q ^ { k }, \dot { q } ^ { k }, \ddot { q } ^ { k }, \dots, q ^ { ( \mu - 1 ) i } \right ). & & ( 1 4 ) \\ \text{if the latter linear system defines the functional form of the generator vector}$$

The solution of the latter linear system defines the functional form of the generator vector X for the infinitesimal transformation (3), (4).

## 2.2 Second-order differential equations

In this study, we focus on second-order differential equations of the form [5]

$$\ddot { q } ^ { i } & = \omega ^ { i } \left ( t, q ^ { k }, \dot { q } ^ { k } \right ). & ( 1 5 ) \\ \text{ents } \eta ^ { [ 1 ] }, \, \eta ^ { [ 2 ] } & \text{ of the second-prolongation reads } [ 5 ]$$

Therefore, the components η [1] , η [2] of the second-prolongation reads [5]

$$\eta ^ { [ 1 ] } & = \eta ^ { i } _ {, t } + q ^ { i } \left ( q ^ { k } \eta ^ { i } _ {, k } - \xi _ {, t } \right ) - \dot { q } ^ { i } \dot { q } ^ { k } \xi _ {, k }, \\ \ell ^ { i } _ {, t k } - \xi _ {, t t } \right ) \dot { q } ^ { k } + \left ( \eta ^ { i } _ {, k r } - 2 \xi _ {, t q } \right ) \dot { q } ^ { k } \dot { q } ^ { r } - \dot { q } ^ { i } \dot { q } ^ { k } \dot { q } ^ { r } \xi _ {, k r } + \left ( \eta ^ { i } _ {, k } - 2 \xi _ {, t } - 3 \xi _ {, k r } \dot { q } ^ { r } \right ) \ddot { q } ^ { k }. \ \ ( 1 7 )$$

η i [2] = η i ,tt +2 ( η i ,tk -ξ ,tt ) ˙ q k + ( η i ,kr -2 ξ ,tq ) ˙ q k ˙ q r -˙ q q i ˙ k ˙ q ξ r ,kr + ( η i ,k -2 ξ ,t -3 ξ ,kr ˙ q r ) ¨ q k . thus the symmetry conditions (14) becomes [5]

$$( 1 7 )$$

$$\begin{array} { r c l c r c l } \cdot & \cdot & \cdot & \cdot \\ & 0 & = & \eta _ { t t } ^ { i } + 2 \left ( \eta _ { t k } ^ { i } - \xi _ { t t } \right ) \dot { q } ^ { k } + \left ( \eta _ { k r } ^ { i } - 2 \xi _ { t q } \right ) \dot { q } ^ { k } \dot { q } ^ { r } \\ & & - \dot { q } ^ { i } \dot { q } ^ { k } \dot { q } ^ { r } \xi _ { k r } + \left ( \eta _ { k } ^ { i } - 2 \xi _ { t } - 3 \xi _ { k r } \dot { q } ^ { r } \right ) \omega ^ { i } \left ( t, q ^ { k }, \dot { q } ^ { k } \right ) \\ & & - \omega _ {, t } ^ { i } - \omega _ {, k } ^ { i } \eta ^ { k } - \omega _ {, \dot { q } ^ { j } } ^ { i } \left ( \eta _ {, t } ^ { j } + q ^ { j } \left ( q ^ { k } \eta _ {, k } ^ { j } - \xi _ {, t } \right ) - \dot { q } ^ { j } \dot { q } ^ { k } \xi _ {, k } \right ). \\ \text{Lie symmetries have numerous applications.  They are used to derive invariant func-} \\ & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdots & \end{array}$$

Lie symmetries have numerous applications. They are used to derive invariant functions and conservation laws. Furthermore, they are applied to categorize differential equations and establish criteria for when differential equations are equivalent to a linear differential equation under a point transformation. In the following lines, we discuss Noether's theorems for the construction of conservation laws and geometric linearization through symmetries.

## 2.3 Noether's theorems

Let's turn our attention to the case where the second-order dynamical system (15) arises from the variation of the action integral

$$S = \int L \left ( t, q, \dot { q } \right ) d t, \\ \dots \cdot \cdot \cdot \cdot \cdot$$

where L t, q, q ( ˙) is the so-called Lagrangian function.

In 1918 [48], Emmy Noether published two groundbreaking theorems that relate the symmetries of the variation principle to conservation laws in dynamical systems.

The first theorem states that if, under the application of the infinitesimal transformations (3), (4), there exists a function f such that the following condition holds

$$X ^ { [ 1 ] } L + L \frac { d \xi } { d t } = \frac { d f } { d t },$$

then X is called a Noether symmetry. It is evident that Noether symmetries are Lie symmetries for the dynamical system (15), but the converse is not necessarily true; that is, Noether symmetries form a subalgebra of the Lie symmetries for the dynamical system. The function f in (20) represents a boundary term introduced to account for infinitesimal changes in the action integral due to infinitesimal changes in the boundary of the domain, caused by infinitesimal transformations of the variables in the action integral.

Noether's second theorem provides a systematic way to derive conservation laws. In particular, if X = ξ∂ t + η∂ i is a Noether symmetry for the dynamical system (15), then the function is a conservation law, that is, A (Φ) = 0; consequently X (Φ) = 0.

$$\Phi \left ( t, q ^ { k }, \dot { q } ^ { k } \right ) & = \xi \left ( \dot { q } ^ { i } \frac { \partial L } { \partial \dot { q } ^ { i } } - L \right ) - \eta ^ { i } \frac { \partial L } { \partial q ^ { i } } + f, \\ \text{ation law, that is, } A \left ( \Phi \right ) & = 0 ; \text{ consequently } X \left ( \Phi \right ) & = 0.$$

For a recent extended discussion of Noether's theorems we refer the reader to [49].

## 3 Linearization through symmetries

One of the significant findings of Lie symmetry analysis is the linearization criterion for a second-order ordinary differential equation. The linearization process is crucial because it provides a method to express the analytic solution in terms of simple functions, thereby offering a better understanding of the dynamics of the dynamical system.

Sophus Lie's theorem states that if a second-order ordinary differential equation ¨ = q ω t, q, q ( ˙) admits eight Lie symmetries forming the sl (3 , R ) Lie algebra, then there exists a transformation that can bring the equation to the form of the free particle equation ¯ q ′′ = 0. For further discussion, see [50] and references therein.

For third-order ordinary differential equations, various approaches have been developed to address their linearization. Criteria for linearization through point and contact symmetries of third-order differential equations have been established in [53, 54]. Additionally, the Sundman transformation as a method for linearization has been studied in [55,56]. More recently, the Cartan equivalency method was considered in [57], where a straightforward procedure was established for the linearization of third-order differential equations using a four-dimensional Lie algebra [57].

The linearization of higher-order ordinary differential equations through Lie symmetries has been extensively discussed in [47, 51, 52] and related references. However, the linearization of partial differential equations holds special interest and has yielded many important results, as seen in [58].

In the context of systems of nonlinear differential equations, the existence of transformations that linearize these equations is highly significant, as it offers an alternative approach to studying the integrability of dynamical systems. Due to the complexity of the problem, various criteria have been proposed in the literature for the linearization of higher-dimensional dynamical systems. For example, in [59], it was demonstrated that a system of second-order differential equations admitting four Lie point symmetries forming the A 4 1 , or A 4 2 , Lie algebra [60] can be transformed into a linear form through a point transformation. Furthermore, the introduction of complex Lie symmetries has led to new linearization criteria for systems of second-order ordinary differential equations, as discussed in the series of studies [61-63].

Nevertheless, the linearization process is inherently a geometric approach. Therefore, studies in the literature have shown that a system of differential equations is linearizable based on specific geometric properties [64]. For systems of second-order ordinary differential equations of the form

$$\dot { q } ^ { i } + \alpha ^ { i } _ { s r j } \left ( t, q ^ { k } \right ) \dot { q } ^ { s } \dot { q } ^ { r } \dot { q } ^ { j } + \beta ^ { i } _ { r j } \left ( t, q ^ { k } \right ) \dot { q } ^ { r } \dot { q } ^ { j } + \gamma ^ { i } _ { r } \left ( t, q ^ { k } \right ) \dot { q } ^ { r } + \delta ^ { i } \left ( t, q ^ { k } \right ) = 0, \\ \text{the coefficients } \alpha ^ { i } _ { s r j }, \, \beta ^ { i } _ { r j }, \, \gamma ^ { i } _ { r }, \, \delta ^ { i } \text{ can be related to the connection coefficient of extended}$$

the coefficients α i srj , β i rj , γ i r , δ i can be related to the connection coefficient of extended manifold [65], where if the connection has zero curvature, that is, the geometry is flat, then

there exist a point transformation which linearize the latter system, for some applications of this method we refer the reader in [66].

Recently, in [69], a novel approach to geometric linearization for a family of secondorder differential equations was discovered. It was found that this new family of dynamical systems can be linearized by introducing new dependent variables. The solution of the extended system can then be expressed in terms of solutions to linear equations. In a similar vein, it was discovered that solutions to Einstein's field equations for certain gravitational models can also be represented in terms of linear equations [70].

In the following section, we establish a new geometric criterion for the global linearization of a family of constraint Hamiltonian systems.

## 4 Constraint Hamiltonian systems

We introduce the Lagrangian function

$$L \left ( N, q ^ { k }, \dot { q } ^ { k } \right ) & = \frac { 1 } { 2 N } g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } - N V ( q ^ { k } ), \\ \text{is a second-rank tensor with inverse } g ^ { i j }, \text{ and } N \left ( t \right ), q ^ { i } \left ( t \right ) \text{ are the } \left ( 1 + n \right )$$

For the Lagrangian (22) it follows det ∣ ∣ ∣ ∂ 2 L ∂Q A ∂Q B ∣ ∣ ∣ = 0, where Q A = ( N t ,q ( ) i ( )). t Hence, the dynamical system described by the Lagrangian (22) is singular.

where g ij ( q k ) is a second-rank tensor with inverse g ij , and N t ( ), q i ( ) t are the (1 + n ) degrees of freedom.

Variation with respect to the variable N t ( ), leads to the constraint equation ∂L ∂N = 0, i.e.

$$\frac { 1 } { 2 N ^ { 2 } } g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + V ( q ^ { k } ) = 0.$$

We introduce the momentum p i = 1 N g ij ˙ q j , thus, the latter constraint reads

$$H \left ( q, p \right ) \equiv \frac { 1 } { 2 } g ^ { i j } p _ { i } p _ { j } + V ( q ^ { k } ) = 0.$$

Function H q,p ( ) represents the Hamiltonian of (22), which is constraint due to existence of expression (24). Dynamical systems described by Lagrangians of the form (22) are of special interests in gravitational physics [67] and in other physical theories (see the discussion in [68]). Nevertheless, any regular Hamiltonian system of the form

$$H _ { R } \left ( q, p, h \right ) \equiv \frac { 1 } { 2 } g ^ { i j } p _ { i } p _ { j } + U ( q ^ { k } ) = h,$$

can be written in the singular form by absorbing the integration constant h , inside the potential, i.e. U q k = V ( q k ) -h .

( ) An important characteristic of this family of constrained dynamical systems is their invariance under time-reparametrization. Below, we demonstrate this property.

The action integral related to the Lagrangian function (22) is

$$S = \int \left ( \frac { 1 } { 2 N } g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } - N V ( q ^ { k } ) \right ) d t, \\.. \quad. \quad. \quad.$$

and the corresponding equations of motion are

$$\ddot { q } ^ { i } + \Gamma _ { j k } ^ { i } \dot { q } ^ { j } \dot { q } ^ { k } + V ^ {, i } - 2 ( \ln N ) ^ {, i } \left ( \frac { 1 } { 2 } g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + V ( q ^ { k } ) \right ) \ & = \ 0, \\ \frac { 1 } { 2 N ^ { 2 } } g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + V ( q ^ { k } ) \ & = \ 0,$$

or due to the constraint equation the second-order differential equation reads

$$\ddot { q } ^ { i } + \Gamma ^ { i } _ { j k } \dot { q } ^ { j } \dot { q } ^ { k } + V ^ {, i } = 0,$$

where Γ i jk is the Levi-Civita connection related to the tensor g ij .

Under the change of the independent variable M τ dτ ( ) = dt , it follows

$$S = \int \tilde { L } \left ( \tilde { N }, q ^ { k }, \frac { d q ^ { k } } { d \tau } \right ) d \tau \, \ \tilde { N } = N \left ( \tau \right ) M \left ( \tau \right ),$$

where is the conformally related Lagrangian. The corresponding equations of motion are derived

$$\tilde { L } \left ( \tilde { N }, q ^ { k }, \frac { d q ^ { k } } { d \tau } \right ) = \frac { 1 } { 2 \tilde { N } } g _ { i j } \frac { d q ^ { i } } { d \tau } \frac { d q ^ { j } } { d \tau } - \tilde { N } V ( q ^ { k } ) \, \\ \text{rmallv related I.arranion} \ \text{The corresponding equations of motion are derived}$$

$$\frac { d ^ { 2 } q ^ { i } } { d \tau ^ { 2 } } + \Gamma _ { j k } ^ { i } \frac { d q ^ { j } } { d \tau } \frac { d q ^ { k } } { d \tau } + V ^ {, i } = 0.$$

$$\frac { 1 } { 2 N ^ { 2 } } g _ { i j } \frac { d q ^ { i } } { d \tau } \frac { d q ^ { j } } { d \tau } + V ( q ^ { k } ) = 0$$

It is evident that the equations of motion remain invariant under a time-reparametrization where the parameter N t ( ),does not affect the dynamical behavior of the constrained system. Specifically, the conformally related Lagrangians L N,q ,q ( k ˙ k ) , ˜ L ( ˜ N,q , k dq k dτ ) have common equations of motions.

## 4.1 Symmetry analysis

The Lie and Noether symmetries for the constraint Hamiltonian systems with Lagrangian of the form (22) have been investigated in details before in [71]. Owing to the constraint (23) the symmetry analysis differs from that of the regular systems, see for example [72]. Indeed for regular dynamical systems and for the equations of motion of the form (29) the Lie symmetries are constructed by the elements of the projective algebra of the connection Γ i jk . Additionally, Noether symmetries are linked to homothetic symmetries of the metric tensor g ij [72].

Indeed if Φ ( t, q k , p k ) is a conservation law for the dynamical system with Hamiltonian (24), it holds

$$\dots \quad \frac { d \Phi } { d t } \equiv \frac { \partial \Phi } { \partial t } + \{ \Phi, H \} = 0,$$

where { } , is the Poisson bracket.

Nevertheless in order to make use of the Hamiltonian constraint the latter condition can be relaxed such that as follows

$$\Phi _ { \epsilon } \quad \frac { d \Phi } { d t } \equiv \frac { \partial \Phi } { \partial t } + \{ \Phi, H \} = \chi H \approx 0,$$

in which χ is a conformal factor.

Thus from (35) it follows that the symmetries of constraint Hamiltonian systems are determined by the conformal symmetries of the metric tensor g ij [71] .

## 4.2 Geodesic description

The process of geometrizing dynamical systems of the form (22) involves formulating the dynamical system as a set of geodesic equations, where the potential term can be interpreted as part of the geometry. Two different approaches for this are the Jacobi metric and the Eisenhart lift, both of which incorporate the potential; that is, conservative forces, into the geometry.

In Jacobi metric approach (for a recent discussion see [28]), we introduce ˜ = N 1 V , such that the Lagrangian to be

$$&, \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \\ & \bar { L } \left ( q ^ { k }, \dot { q } ^ { k } \right ) = \frac { 1 } { 2 } \bar { g } _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } - 1 \, \, \bar { g } _ { i j } = V \left ( q \right ) g _ { i j }. \\ & \ n \text{straint equation (23) becomes}$$

Thus, the constraint equation (23) becomes

$$\frac { 1 } { 2 } \bar { g } _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + 1 = 0$$

which means that Lagrangian (36) describes the time/space -like geodesic equations for the conformal metric ¯ g ij . It is well known that the geodesic equations of the metric ¯ g ij which can be linearized are those of the flat space, that is, ¯ g ij should be flat.

The Jacobi metric is obtained through a conformal transformation where ¯ gij is conformally related to the metric tensor gij . Conformally related spacetimes share the same conformal symmetries, implying that the symmetries and conservation laws of the original dynamical system also exist for the Jacobi metric. Notably, conformal symmetries in this context transform into isometries. On the other hand, the Eisenhart lift approach involves augmenting the dimensionality of the system [29,73].

Geometrization is achieved by introducing additional dimensions through new dependent variables. The potential term becomes integrated into the newly extended metric tensor. These extended spaces possess isometries that correspond to Noetherian conservation laws for the geodesic equations. When these conservation laws are applied within the extended system, the original dynamics are recovered.

In the Eisenhart approach, we write the extended Lagrangian function

$$L _ { l i f t } \left ( N, q ^ { k }, \dot { q } ^ { k }, z, \dot { z } \right ) & = \frac { 1 } { 2 N } \left ( g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + \frac { 1 } { V ( q ^ { k } ) } \dot { z } ^ { 2 } \right ), \\ \text{asraint equation and conservation law}$$

with constraint equation and conservation law

$$g _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + \frac { 1 } { V ( q ^ { k } ) } \dot { z } ^ { 2 } = 0 \, \ \frac { 1 } { V \left ( q ^ { k } \right ) } \dot { z } = I _ { 0 }.$$

Thus by replacing the ˙ = z I V 0 ( q k ) in the equations of motion and the constraint equation, the original system is recovered when ( I 0 ) 2 = 1.

Lagrangian function (38) describe a sets of null geodesic equations in the extended manifold with metric with A, B = i, j.., z , are the indices of the extended geometry.

$$\text{ metric} \hat { g } _ { A B } \left ( q ^ { k }, z \right ) = g _ { i j } \left ( q ^ { k } \right ) + \frac { 1 } { V \left ( q ^ { k } \right ) } \delta _ { A } ^ { z } \delta _ { B } ^ { z }, & & ( 4 0 ) \\ i.., z, & \text{ are the indices of the extended geometry.}$$

## 5 Geometric Linearization

In this section, we will establish new geometric criteria for the global linearization of the equations of motion described by the Lagrangian (22).

We will utilize the two geometrization approaches described previously. Specifically, for the equivalent Lagrangian (36) obtained through the Jacobi metric approach, we will

employ the Eisenhart lift. This allows us to describe the dynamical system using the equivalent/extended Lagrangian function

$$\bar { L } _ { l i f t } = \frac { 1 } { 2 } \bar { g } _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + \frac { 1 } { 2 } \dot { z } ^ { 2 }$$

with constraints

$$\bar { g } _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } + \dot { z } ^ { 2 } = 0 \, \ \dot { z } = \pm 1.$$

Lagrangian (41) describes the equations of motion for the null geodesics of the extended metric tensor ˜ g AB q k , z = ¯ g ij q k + δ z A B δ z .

( ) ( ) Consequently, if the metric tensor ˜ g AB is conformally flat then there exist a point transformation ( q A , z ) = Q A such that ˜ g AB = e 2 U η AB , where η AB is the flat space in diagonal coordinates. Thus, Lagrangian (41) and the constraint (42) becomes

$$\bar { L } _ { l i f t } = e ^ { 2 U } \eta _ { A B } \dot { Q } ^ { A } \dot { Q } ^ { B },$$

$$\eta _ { A B } \dot { Q } ^ { A } \dot { Q } ^ { B } = 0.$$

Therefore after the change of the independent variable dτ = e -2 U dt , we end with the equivalent geodesic equations

$$\ddot { Q } ^ { A } = 0$$

which are the equations of motion for the free particle in the flat space.

Thus, regarding the linearization of a dynamical system described by the singular Lagrangian function (22) the following theorem states.

Theorem 1 The n -dimensional, with n ≥ 2 , constraint Hamiltonian system described by the Lagrangian function (22) is globally linearizable if one of the following equivalent statements are true

- A) The admitted (non trivial) Noether symmetries of the constraint Hamiltonian system are n n ( +1) 2 .
- B) The Jacobi metric ¯ g ij is maximally symmetric.
- C) The extended 1 + n decomposable space with metric ˜ g AB is conformally flat.

Proof. All statements of Theorem 1 are equivalent. As we discussed before if the extended 1+ n decomposable space with metric ˜ g AB is conformally flat then the equations of motions are linearizable through point transformations in the space { q A , z } .

Nevertheless, ˜ g AB is conformally flat if and only if the n -dimensional space with metric ¯ g ij is maximally symmetric (see Proposition two in [74]).

However, when ¯ g ij is maximally symmetric it means ¯ g ij admits n n ( +1) 2 isometries, as many as the Noether symmetries for the original Lagrangian (22).

The inverse proof is straightforward.

It is important to note that two-dimensional spaces are maximally symmetric when they have constant curvature. Conversely, three-dimensional spaces are conformally flat when the Cotton-York tensor vanishes. Therefore, for the case of two-dimensional systems, from Theorem 1, we derive the following corollary.

Corollary 2 The two-dimensional constraint Hamiltonian systems described by the Lagrangian function (22) are globally linearizable if one of the following equivalent states are true

- A) The admitted Noether symmetries of the constraint Hamiltonian system are three.
- B) There exist a point transformation where the two-dimensional Jacobi metric can be written as ¯ g ij = diag ( e 2 U x,y ( ) , e 2 U x,y ( ) ) , and U x,y ( ) is a solution of the equation U ,xx + U ,yy +2 κe 2 U = 0 , in which κ is the curvature of the two-dimensional space and V ( x, y ) = e -2 U x,y ( )
- C) The Cotton-York tensor C µνκ = R µν κ ; -R κν µ ; + 1 4 ( R g ; ν µκ -R g ; κ µν ) for the threedimensional extended metric ˜ g AB has zero components.

Corollary 3 The one-dimensional constraint Hamiltonian system described by the Lagrangian function (22) is always globally linearizable for arbitrary potential function.

Proof. Consider the one-dimensional constraint Hamiltonian system with Lagrangian function L N,q,q ( ˙) = 1 2 N ˙ q 2 -NV ( q ), then after the application of Jacobi's and Eisenhart's methods we find the equivalent dynamical system of null geodesics for the two-dimensional space. All two-dimensional spaces are conformally flat which means that the null geodesics can be globally linearizable. In particular the extended two-dimensional space is described by the decomposable line element ds 2 = 1 V ( q ) dq 2 + dz 2 . Thus, under the change of variable √ 1 V ( q ) dq = dY , it follows ds 2 = dY 2 + dz 2 , which is the flat space.

On the other hand, the Jacobi metric approach leads to the equivalent system of time/space -like geodesic equations for the one-dimensional line element ds ¯ 2 = 1 V ( q ) dq 2 . Because the one-dimensional space is the flat space the point transformation √ 1 V ( q ) dq = dY leads to the linearization of the dynamical system.

## 6 Applications

In the following lines, we demonstrate the application of Theorem 1 for the geometric linearization of Hamiltonian systems of special interests.

## 6.1 Exponential interaction

Consider two particles which interact with an exponential law, similar to that of the Toda lattice. The Lagrangian which describes the dynamical system is

$$L \left ( q _ { 1 }, \dot { q } _ { 1 }, q _ { 2 }, \dot { q } _ { 2 } \right ) & = \frac { 1 } { 2 } \left ( \dot { q } _ { 1 } ^ { 2 } + \dot { q } _ { 2 } ^ { 2 } \right ) - V _ { 0 } e ^ { 2 \gamma ( q _ { 1 } - q _ { 2 } ) }. \\ \text{Iamiltonian function is } H & = h.$$

where the Hamiltonian function is H = h .

We introduce the equivalent singular Lagrangian

$$\bar { L } \left ( N, q _ { 1 }, \dot { q } _ { 1 }, q _ { 2 }, \dot { q } _ { 2 } \right ) & = \frac { 1 } { 2 N } \left ( \dot { q } _ { 1 } ^ { 2 } + \dot { q } _ { 2 } ^ { 2 } \right ) - N \left ( V _ { 0 } e ^ { q _ { 1 } - q _ { 2 } } - h \right ), \\ \text{constraint equation}$$

with constraint equation in which p 1 = ˙ , q 1 p 2 = ˙ . q 2

$$\bar { H } \equiv \frac { 1 } { 2 } \left ( p _ { 1 } ^ { 2 } + p _ { 2 } ^ { 2 } \right ) + V _ { 0 } e ^ { q _ { 1 } - q _ { 2 } } - h = 0, \\ i _ { 1 }, \, p _ { 2 } = \dot { q } _ { 2 }.$$

We employ the Eisenhart lift and we write the equivalent Lagrangian of the form (38), that is,

$$\text{na is}, & \bar { L } _ { l i f t } \left ( N, q _ { 1 }, \dot { q } _ { 1 }, q _ { 2 }, \dot { q } _ { 2 } \right ) = \frac { 1 } { 2 N } \left ( \left ( \dot { q } _ { 1 } ^ { 2 } + \dot { q } _ { 2 } ^ { 2 } \right ) + \frac { 1 } { \left ( V _ { 0 } e ^ { q _ { 1 } - q _ { 2 } } - h \right ) } \dot { z } ^ { 2 } \right ), \\ \text{where the extended three-dimensional metric reads}$$

where the extended three-dimensional metric reads

$$d s _ { l i f t } ^ { 2 } = d q _ { 1 } ^ { 2 } + d q _ { 2 } ^ { 2 } + \frac { 1 } { ( V _ { 0 } e ^ { q _ { 1 } - q _ { 2 } } - h ) } d z ^ { 2 }. \\ \intertext { r e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e f e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e }$$

We calculate the Cotton-York tensor and we found that it is zero when h = 0. Thus, from Corollary 2 it follows that the dynamical system can be written in the equivalent form of the three-dimensional free particle of the flat geometry.

Under the he conformal transformation ds ¯ 2 lift = 1 XY ds 2 lift and the change of variables

$$i x \ = \ - \frac { 1 } { 4 } \left ( 1 + i \right ) \ln Y + \frac { 1 } { 4 } \left ( 1 - i \right ) \ln X,$$

$$i y \ = \ \frac { 1 } { 4 } \, ( 1 - i ) \, ( i \ln X - \ln Y ) \,,$$

the extended space ds ¯ 2 lift becomes

$$d s _ { l i f t } ^ { 2 } = \frac { 1 } { 2 } d X d Y + d z ^ { 2 }.$$

where the null geodesics read

$$\ddot { X } = 0 \, \ddot { Y } = 0 \, \ \ddot { z } = 0,$$

and constraint equation

$$\frac { 1 } { 2 } \dot { X } \dot { Y } + \dot { z } ^ { 2 } = 0.$$

## 6.2 Two-dimensional oscillator with corrections

We assume the following singular Lagrangian, which describes a two-dimensional oscillator with correction terms, that is,

$$L \left ( N, x, \dot { x }, y, \dot { y } \right ) = \frac { 1 } { 2 N } \left ( \dot { x } ^ { 2 } + \dot { y } ^ { 2 } \right ) + N \left ( 1 + \frac { \kappa } { 4 } \left ( x ^ { 2 } + y ^ { 2 } \right ) \right ) ^ { 2 }. \\ \text{We calculate the Jacobi metric, that is, the corresponding line element is of the form}$$

We calculate the Jacobi metric, that is, the corresponding line element is of the form

$$d s _ { J a c o b i } ^ { 2 } = \frac { 1 } { \left ( 1 + \frac { \kappa } { 4 } \left ( x ^ { 2 } + y ^ { 2 } \right ) \right ) ^ { 2 } } \left ( d x ^ { 2 } + d y ^ { 2 } \right ) \\ \text{we observe that is a space of constant $\kappa$ nonzero curvature, that it, it is a}$$

from where we observe that is a space of constant κ nonzero curvature, that it, it is a maximally symmetric space.

We employ Eisenhart's lift and we write the equivalent extended Lagrangian with constraint

$$L _ { l i f t } \left ( N, x, \dot { x }, y, \dot { y } \right ) & = \frac { 1 } { 2 } \frac { \dot { x } ^ { 2 } + \dot { y } ^ { 2 } } { \left ( 1 + \frac { \tilde { x } } { 4 } \left ( x ^ { 2 } + y ^ { 2 } \right ) \right ) ^ { 2 } } + \frac { 1 } { 2 } \dot { z } ^ { 2 }, \\ \i \text{straint} & 1 \quad \dot { x } ^ { 2 } + \dot { u } ^ { 2 } \quad 1 \ \. \quad \.$$

For the three-dimensional extended space with line element

$$\frac { 1 } { 2 } \frac { \dot { x } ^ { 2 } + \dot { y } ^ { 2 } } { \left ( 1 + \frac { \kappa } { 4 } \left ( x ^ { 2 } + y ^ { 2 } \right ) \right ) ^ { 2 } } + \frac { 1 } { 2 } \dot { z } ^ { 2 } = 0. \text{ \quad \ } \text{$\text{$\text{$line element}$} }$$

$$d s _ { L i f t } ^ { 2 } = d s _ { J a c o b i } ^ { 2 } + d z ^ { 2 }$$

we calculate the Cotton-York tensor components which are zero. Hence, the line element (60) is conformally flat.

Under the change of variables

$$x = \frac { X \left ( Z + \sqrt { Z ^ { 2 } + 1 6 \kappa \left ( X ^ { 2 } + Y ^ { 2 } \right ) } \right ) } { 2 \kappa \left ( X ^ { 2 } + Y ^ { 2 } \right ) } \, \ y = \frac { Y \left ( Z + \sqrt { Z ^ { 2 } + 1 6 \kappa \left ( X ^ { 2 } + Y ^ { 2 } \right ) } \right ) } { 2 \kappa \left ( X ^ { 2 } + Y ^ { 2 } \right ) },$$

$$e ^ { \sqrt { K } z } & = 2 \kappa \frac { ( X ^ { 2 } + Y ^ { 2 } ) } { Z + \sqrt { Z ^ { 2 } + 1 6 \kappa \, ( X ^ { 2 } + Y ^ { 2 } ) } } \left ( \frac { 8 \, ( X ^ { 2 } + Y ^ { 2 } ) } { \left ( Z + \sqrt { Z ^ { 2 } + 1 6 \kappa \, ( X ^ { 2 } + Y ^ { 2 } ) } \right ) } + Z \right ), \\ & \text{the three-dimensional space } ( 6 0 ) \, \text{ reads}$$

the three-dimensional space (60) reads

$$d s _ { L i f t } ^ { 2 } = e ^ { \sqrt { K } z ( X, Y, Z ) } \left ( d X ^ { 2 } + d Y ^ { 2 } + d Z ^ { 2 } \right ), \\ \text{desc equations}$$

with null geodesic equations

$$\ddot { X } = 0 \, \ \ddot { Y } = 0 \, \ \ddot { Z } = 0,$$

and constraint ˙ X 2 + ˙ Y 2 + ˙ Z 2 = 0.

This two-dimensional linearizable example can be generalized to n -dimensional system in a similar way.

Corollary 4 The n -dimensional dynamical system with Lagrangian

$$L \left ( N, q ^ { k }, \dot { q } ^ { k } \right ) = \frac { 1 } { 2 N } \left ( \eta _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } \right ) + N \left ( 1 + \frac { \kappa } { 4 } \left ( \eta _ { i j } q ^ { i } q ^ { j } \right ) \right ) ^ { 2 }, \\ e \text{ linearized though the Jacobi metric and the Eisenhart lift.}$$

can be linearized though the Jacobi metric and the Eisenhart lift.

## 6.3 The Szekeres system

Consider the two-dimensional dynamical system described by the regular Lagrangian function [75, 76]

$$L _ { R } \left ( u, \dot { u }, v, \dot { v } \right ) = \dot { u } \dot { v } - \frac { v } { u ^ { 2 } },$$

and equations of motion

$$\ddot { u } + u ^ { - 2 } = 0 \, \ \ddot { v } - 2 v u ^ { - 3 } = 0.$$

Furthermore, the Hamiltonian function is

$$H \equiv \dot { u } \dot { v } + \frac { v } { u ^ { 2 } } = h$$

where h is a constant.

The system described above corresponds to the dynamics governed by the Einstein field equations for Szekeres spacetimes, commonly referred to as the Szekeres system [77]. Szekeres spacetimes represent exact inhomogeneous solutions with a dust fluid and find various applications in gravitational physics. For more detailed information, the interested reader is referred to [78].

We introduce the equivalent singular Lagrangian function

$$L \left ( N, u, \dot { u }, v, \dot { v } \right ) = \frac { \dot { u } \dot { v } } { N } - N \left ( \frac { v } { u ^ { 2 } } - h \right ), \\ \text{rite the original system in the form of a constraint Hamiltonian dynamical}$$

such that to write the original system in the form of a constraint Hamiltonian dynamical system.

The latter singular Lagrangian for arbitrary parameter h has not any Noetherian symmetry. In particular it admits a hidden symmetry related to a quadratic conservation law [75]. However, for h = 0, Lagrangian (69) admits three Noether symmetries

$$X _ { 1 } = \frac { 1 } { v } \partial _ { v } \, \ X _ { 2 } = u ^ { 2 } \partial _ { u } \, \ X _ { 3 } = v \partial _ { v } + 2 u \partial _ { u }.$$

Thus, according to the first statement of Theorem 1, on the surface with h = 0, the equations of motions can be linearized. Moreover, the solution space with initial conditions h = 0, is equivalent to the equations of motion for the two-dimensional flat space.

Indeed, we employ Jacobi's approach and from (69) we define the equivalent geodesic Lagrangian

$$\bar { L } \left ( u, \dot { u }, v, \dot { v } \right ) = \frac { v } { u ^ { 2 } } \dot { u } \dot { v } - 1.$$

Hence, under the change of variables

$$\frac { d u } { u ^ { 2 } } = d U, \ v d v = d V,$$

it follow with equations of motion

and constraint equation

$$L \left ( U, \dot { U }, V, \dot { V } \right ) = \dot { U } \dot { V } - 1, \\ \intertext { t i o n }$$

$$\ddot { U } = 0, \, \ddot { V } = 0$$

$$\dot { U } \dot { V } + 1 = 0.$$

## 6.4 The Cosmological Constant in Szekeres model

The introduction of the cosmological constant in the Szekeres model [79] leads to the modification of the regular Lagrangian (66).

Specifically, the new dynamics follows from the Lagrangian function [80]

$$L _ { R } ^ { \Lambda } \left ( u, \dot { u }, v, \dot { v } \right ) = \dot { u } \dot { v } - \left ( \frac { v } { u ^ { 2 } } - \Lambda u v \right ), \\ \intertext { e \cos \text{mological constant term} } \text{ }$$

in which Λ is the cosmological constant term.

In a similar approach with before we introduce the singular Lagrangian function

$$L ^ { \Lambda } \left ( N, u, \dot { u }, v, \dot { v } \right ) = \dot { u } \dot { v } - \left ( \frac { v } { u ^ { 2 } } - \Lambda u v - h \right ) \\ \intertext { e Jacobi metric is defined as }$$

Thus, the Jacobi metric is defined as

$$d s _ { J a c o b i } ^ { 2 } & = \frac { 1 } { \left ( \frac { v } { u ^ { 2 } } - \Lambda u v - h \right ) } d u d v. \\ \text{Ricciscalor for the two-dimensional space, it is,}$$

We calculate the Ricciscalar for the two-dimensional space, it is,

$$R _ { J a c o b i } = - \frac { 4 \left ( 2 + \Lambda u ^ { 3 } \right ) h } { u \left ( v \left ( \Lambda u ^ { 3 } - 1 \right ) + h u ^ { 2 } \right ) }. \\ \intertext { a c o b i } \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{array}$$

Hence, for h = 0, the two-dimensional Jacobi metric (78) describes the flat space, that is, the equations of motion can be linearized though a point transformation.

Indeed, under the point transformation the Jacobi metric reads

$$\frac { d u } { ( \frac { 1 } { u ^ { 2 } } - \Lambda u ) } = d U \, \ \frac { d v } { v } = d V, \\ \intertext { a d s } \cdot \theta \quad \dots \dots \quad \theta \dots$$

$$d s _ { J a c o b i } ^ { 2 } = d U d V,$$

which leads to the geodesic equations of the flat space, that is,

$$\ddot { U } = 0 \, \ \ddot { V } = 0.$$

## 6.5 Static spherical symmetric spacetime with charge

Einstein's gravitational field equations for the a static spherical symmetric spacetime with a charge are described by the variation of the singular Lagrangian [81]

$$L ^ { R N } \left ( N, a, a ^ { \prime }, b, b ^ { \prime }, \zeta, \zeta ^ { \prime } \right ) = \frac { 1 } { 2 N } \left ( 8 b a ^ { \prime } b ^ { \prime } + 4 a b ^ { \prime 2 } + 4 \frac { b ^ { 2 } } { a } \zeta ^ { \prime 2 } \right ) + 2 N a,$$

where N,a, b are the scale factors of the background geometry with line element

$$d s ^ { 2 } = - a \left ( r \right ) ^ { 2 } d t ^ { 2 } + N \left ( r \right ) ^ { 2 } d r ^ { 2 } + b ^ { 2 } \left ( r \right ) \left ( d \theta ^ { 2 } + \sin ^ { 2 } \theta d \varphi ^ { 2 } \right ), \\ \intertext { a prime means total derivative with respect the independent parameter $r$, i.e. $a^{\prime}=\frac{d}{d}$}.$$

and a prime means total derivative with respect the independent parameter r , i.e. a ′ = da dr . Function ζ is related to the charge. The solution of the field equations is known as the Reissner-Nordstr¨m black hole [82,83] o

The kinetic term of the singular Lagrangian (83) is defined by a three-dimensional space. The admitted Noether symmetries of Lagrangian (83) are calculated to be six [81], they are [70]

$$\infty \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array} \, \begin{array} { c } \infty \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \col. \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \code. \end{array} \, \begin{array} { c } \infty \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \code. \end{array} \, \infty \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \code. \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \\ X ^ { 1 } = \frac { 2 \zeta } { a b } \partial _ { a } \partial _ { b } \partial _ { \zeta } \,, \ X ^ { 2 } = - a \partial _ { \zeta } \partial _ { \zeta } \,, \ X ^ { 2 } = \partial _ { \zeta } \partial _ { \zeta } \,. \end{array} \, \begin{array} { c } \infty \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \ colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \\ \colon A \text{ of Theorem $1$ states the field equations can be wr} \quad. \quad.$$

Therefore, case A of Theorem 1 states the field equations can be written in the equivalent form of the free particle. This transformation derived before in [70] where the common solution space for a large families of gravitational models investigated.

We employ the Eisenhart lift and we write the extended line element

$$d s _ { l i f t } ^ { 2 \ m } = \frac { 1 } { N } \left ( 8 b \ d a \ d b + 4 a \ d b ^ { 2 } + 4 \frac { b ^ { 2 } } { a } d \zeta ^ { 2 } - \frac { d \psi ^ { 2 } } { 2 a } \right ) \\ \intertext { a n e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e m e s e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e }$$

and under the point transformation

$$a = \sqrt { \frac { X + Y } { X - Y } + \frac { z ^ { 2 } } { \left ( X - Y \right ) ^ { 2 } } }, \ \zeta = \frac { z } { \left ( X - Y \right ) ^ { 2 } },$$

the latter line element becomes

$$d s ^ { R N \ 2 } = \frac { 1 } { n } \frac { X - Y } { \sqrt { X ^ { 2 } - Y ^ { 2 } + z ^ { 2 } } } \left ( 4 d X ^ { 2 } - 4 d Y ^ { 2 } + 4 d z ^ { 2 } - d \psi ^ { 2 } \right ) \\ \text{ere the corresponding null geodesics are written in the linear form}$$

where the corresponding null geodesics are written in the linear form

$$X ^ { \prime \prime } = 0 \, \ Y ^ { \prime \prime } = 0, \ z ^ { \prime \prime } = 0 \, \ \psi ^ { \prime \prime } = 0.$$

and Hamiltonian constraint

$$\dot { X } ^ { 2 } - \dot { Y } ^ { 2 } + \dot { z } ^ { 2 } - \frac { 1 } { 4 } \dot { \psi } ^ { 2 } = 0.$$

## 7 Conclusions

In this study, by using the Jacobi metric and the Eisenhart lift, we established a new criterion for the global linearization of constrained Hamiltonian systems. The requirements for the dynamical system are to be in the form of (22) with the constraint expression (23). The n -dimensional dynamical system must admit n n ( +1) 2 Noether point symmetries, which correspond to a number of n n ( +1) 2 independent conservation laws. This property indicates that the given dynamical system, constrained by equation (23), posses the property of superintegrability.

The linearization of the dynamical systems is achieved through geometry and is based on the linearization of the equivalent system that describes geodesic equations for an extended geometry. The main result of this analysis is summarized in Theorem 1, where three equivalent statements related the number of admitted symmetries, and the geometric characteristics of the Jacobi metric and of the extended Eisenhart metric are given.

Two immediate results are given in Corollaries 2 and 3. Corollary 2 specialize the statements of Theorem 1 for the case of two-dimensional dynamical systems, while Corollary 3 states that all one-dimensional constraint Hamiltonian systems can be globally linearized. This geometric linearization for the one-dimensional systems is possible, because the two-dimensional extended Eisenhart space is always conformally flat.

To demonstrate the application of this new geometric approach, we presented a series of applications of special interest. We focused on two dynamical systems of analytic mechanics and on some gravitational models. Specifically we consider the Szekeres system, which describes inhomogeneous cosmological models with or without the cosmological constant term, and the Reissner-Nordstr¨m black hole. o

A criticism that can be made is that this specific algorithm fails in the case of the simplest maximally symmetric system, which is that of the n -dimensional harmonic oscillator with Lagrangian

We should make it clear that the geometric linearization discussed in this study is based on the existence of symmetries generated by conformal symmetries of the metric tensor g ij , which defines the kinetic term in (22). As far as the symmetries of the oscillator (90) are concerned, they are related to the projective algebra [84] and not to the conformal algebra. Nevertheless, only the isometries are the common subalgebra for the projective

$$L \left ( N, q ^ { k }, \dot { q } ^ { k } \right ) = \frac { 1 } { 2 } \eta _ { i j } \dot { q } ^ { i } \dot { q } ^ { j } - \frac { \omega ^ { 2 } } { 2 } \eta _ { i j } q ^ { i } q ^ { j }. & & ( 9 0 ) \\ \text{make it clear that the geometric linearization discussed in this study is}$$

and conformal algebras. An alternative Eisenhart lift has been proposed in [69], where the oscillator is reduced to the free particle. This is possible should the introduction of an extended space which belongs to the pp-wave geometries. In this case the dimension of the extended space is increased by two and the resulting Eisenhart metric is conformally flat.

Thus, any dynamical system of the form (22) where an extended Eisenhart metric can be constructed to be conformally flat, can be linearized. Which means that Theorem 1 can be generalized for other family of Eisenhart lifts.

This work opens new directions in the study of nonlinear differential equations and demonstrates that geometry is a powerful tool for the study of dynamical systems.

Anatural extension of this geometric consideration is for the study of the Klein-Gordon equation which describes the quantum limit of the constraint Hamiltonian system with Lagrangian (22).

Specifically, the quantization of the constraint Hamiltonian (24) leads to the Yamabe equation [85]

where

$$\hat { \Delta } \Psi + V \left ( q ^ { k } \right ) \Psi = 0, \quad \ \ ( 9 1 ) \\ \hat { \Delta } \Psi + \frac { 0 } { 2 }$$

$$\hat { \Delta } & = \Delta + \frac { n - 2 } { 4 \left ( n - 1 \right ) } R \\ \underline { \partial } \right \rangle \text{ is the manual language operator and } R \text{ is the Riscal scalar for}$$

and ∆ = 1 √ | g | ∂ ∂q i ( √ gg ij ∂ ∂x j ) is the usual Laplace operator and R is the Ricci scalar for the metric tensor g ij which defines the kinetic term and n = dim g ij . The introduction of the second term in (92) is necessary in order equation the equation to be invariant under conformal transformations.

We demonstrate the application of the geometric approach in equation (91) let us study the Yamabe equation for the Szekeres system (69) with h = 0. The equation which describes the quantum system is

Assume now equation

$$\frac { u ^ { 2 } } { v } \left ( \Psi _ {, u v } - \frac { v } { u ^ { 2 } } \Psi \right ) = 0, \ \Psi = \Phi \left ( u, v \right ) \\ \text{$\text{emation} }$$

$$\frac { u ^ { 2 } } { v } \Phi _ {, u v } + \Phi _ {, z z } = 0 \, \ \Phi = \Phi \left ( u, v, z \right ).$$

The vector field ∂ z -i Φ ∂ Φ , is a Lie symmetry for equation (94). The corresponding invariants are { u, v, Φ = Ψ e iz } . Thus by replacing in (94) we end with equation (93).

Under the change of variables (72), and U = X + Y , V = X -Y , equation (94) reads

$$\Phi _ {, X X } - \Phi _ {, Y Y } + \Phi _ {, z z } = 0,$$

which is the wave equation for the three-dimensional flat space. Equation (95) is maximally symmetric and admits ten Lie point symmetries plus the infinity symmetries related to the infinity number of solutions of the linear equation (for more details we refer to [86]). On the other hand, equation (93) admits only three Lie point symmetries (plus the infinity symmetries).Therefore, the symmetries of the maximally symmetric equation (94) which does not survive under the reduction with the invariants { u, v, Φ = Ψ e iz } , becomes nonlocal symmetries, which can be used for the construction of new solutions for the inhomogeneous equation (93) related to solutions for the homogeneous equation (95)

In a future work, we plan to investigate further this geometric consideration for the analysis of other dynamical systems and of partial differential equations.

Data Availability Statements: Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Acknowledgements: The author thanks the support of Vicerrector´ ıa de Investigaci´n y Desarrollo Tecnol´gico (Vridt) at Universidad Cat´lica del Norte through N´cleo o o o u de Investigaci´n Geometr´ ıa Diferencial y Aplicaciones, o Resoluci´n Vridt No - 096/2022 o and Resoluci´n Vridt No - 098/2022. o

## References

- [1] S. Lie,Theorie der Transformationsgrupprn: Vol I, Chelsea, New York (1970)
- [2] S. Lie, Theorie der Transformationsgrupprn: Vol II, Chelsea, New York (1970)
- [3] S. Lie,Theorie der Transformationsgrupprn: Vol III, Chelsea, New York (1970)
- [4] N.H. Ibragimov, CRC Handbook of Lie Group Analysis of Differential Equations, Volume I: Symmetries, Exact Solutions, and Conservation Laws, CRS Press LLC, Florida (2000)
- [5] H. Stephani, Differential Equations: Their Solutions Using Symmetry, Cambridge University Press, New York, (1989)

- [6] P.J. Olver, Applications of Lie Groups to Differential Equations, Springer-Verlag, New York, (1993)
- [7] G.W. Bluman and S. Kumei, Symmetries of Differential Equations, Springer-Verlag, New York, (1989)
- [8] S. Meleshko and N. Samatova, Group classification of the two-dimensional shallow water equations with the beta-plane approximation of coriolis parameter in Lagrangian coordinates, Comm. Nonl. Sci. Num. Sim. 90, 105337 (2020)
- [9] A.A. Chesnokov, Symmetries and exact solutions of the rotating shallow-water equations, J. Appl. Mech. Techn. Phys. 49, 737 (2008)
- [10] M. Torrisi and R. Tracina, Symmetries and Conservation Laws for a Class of FourthOrder Reaction-Diffusion-Advection Equations, Symmetry 15, 1936 (2023)
- [11] L. Blanco Diaz, C. Sardon, F.J. Alburquerque and J. de Lucas, Geometric Numerical Methods for Lie Systems and Their Application in Optimal Control, Symmetry 15, 1285 (2023)
- [12] V.A.Dorodnitsyn and E.I. Kaptsov, Shallow water equations in Lagrangian coordinates: Symmetries, conservation laws and its preservation in difference models, Comm. Nonl. Sci. Num. Sim. 89, 105343 (2020)
- [13] A.V. Aminova, Projective transformations and symmetries of differential equations, Sbornik Mathematics 186 (12), 1711 (1995)
- [14] H. Liu and Y. Yun, Lie Symmetry Analysis and Conservation Laws of Fractional Benjamin-Ono Equation, Symmetry 16, 473 (2024)
- [15] A.H. Bokhari, A. H. Kara, A. R. Kashif, and F.D. Zaman, Noether Symmetries Versus Killing Vectors and Isometries of Spacetimes, International Journal of Theoretical Physics, 45, 1029 (2006)
- [16] A.H. Bokhari, Johnpillai, A.G., Kara, A.H., Mahomed, F.M., and F.D. Zaman, Classification of Static Spherically Symmetric Spacetimes by Noether Symmetries, International Journal of Theoretical Physics 52, 3534, (2013)

- [17] F. Oliveri, Lie Symmetries of Differential Equations: Classical Results and Recent Contributions, Symmetry 2, 658 (2010)
- [18] J. Belmonte-Beitia, V.M. P´ erez-Garc´ ıa, Vadym Vekslerchik, and Pedro J. Torres, Lie Symmetries and Solitons in Nonlinear Systems with Spatially Inhomogeneous Nonlinearities, Phys. Rev. Lett. 98, 064102 (2007)
- [19] S. Jamal, Quadratic integrals of a multi-scalar cosmological model, Mod. Phys. Lett. A. 35, 2050068 (2020)
- [20] S.A. Hojman, A new conservation law constructed without using either Lagrangians or Hamiltonians, J. Phys. A.: Math. Gen. 25, L291 (1992)
- [21] T. Pillay and P.G.L. Leach, Comment on a theorem of Hojman and its generalizations, J. Phys. A.: Math. Gen. 29, 6999 (1996)
- [22] H.-B. Zhang and L.-Q. Chen, The Unified Form of Hojman's Conservation Law and Lutzky's Conservation Law, J. Phys. Soc. Jpn. 74, 905 (2005)
- [23] S. Chanda, G.W. Gibbons and P. Guha, Jacobi-Maupertuis metric and Kepler equation, Int. J. Geom. Meth. Mod. Phys. 14, 1730002 (2017)
- [24] M. Szyd/suppressowski, l Desingularization of Jacobi metrics and chaos in general relativity, J. Math. Phys. 40, 3519 (1999)
- [25] A. Duenas-Vidal and O.L. Andino, The Jacobi metric approach for dynamical wormholes, Gen. Relat. Grav. 55, 9 (2023)
- [26] M. Szydlowski, M. Heller and W. Sasin, Geometry of spaces with the Jacobi metric, J. Math. Phys. 37, 346 (1996)
- [27] G.W. Gibbons, The Jacobi metric for timelike geodesics in static spacetimes, Cl.ass. Quantum Grav. 33, 025004 (2016)
- [28] P. Maraner, On the Jacobi metric for a general Lagrangian system, J. Math. Phys. 60, 112901 (2019)
- [29] L.P. Eisenhart, Dynamical Trajectories and Geodesics, Annals. Math. 30, 591 (1928)

- [30] C. Duval, G. Burdet, H.P. Kunzle and M. Perrin, Bargmann structures and NewtonCartan theory, Phys. Rev. D 1841, 31 (1985)
- [31] T. Kaluza, Zum unit¨tsproblem der physik Itzungsber, Preuss. Akad. Wiss. Berlin a 966 (1921)
- [32] O. Klein, The atomicity of electricity as a quantum theory law, Nature 118, 516 (1926)
- [33] A.P. Fordy and A. Galajinsky, Eisenhart lift of 2-dimensional mechanics, Eur. Phys. J. C 79, 301 (2019)
- [34] M. Cari˜ena, F. Jose, n F.J. Herranz, M. Fern´ andez-Ra˜ada Men´ endez n De Luarca, Superintegrable systems on 3-dimensional curved spaces: Eisenhart formalism and separability, J. Math. Phys. 58, 022701 (2017)
- [35] M. Cariglia and G.W. Gibbons, Generalised Eisenhart lift of the Toda chain, J. Math. Phys. 55, 022701 (2014)
- [36] H. Zhang and Q.-Y. Hao, Eisenhart lift for Euler's problem of two fixed centers, Appl. Math. Compt. 350, 305 (2019)
- [37] K. Fin, S. Karamitsos and A. Pilaftsis, Quantizing the Eisenhart lift, Phys. Rev. D 103, 065004 (2021)
- [38] K. Fin, The Eisenhart Lift. In: Geometric Approaches to Quantum Field Theory, Springer Theses, Springer, Cham (2021)
- [39] M. Cariglia, Hidden symmetries of Eisenhart-Duval lift metrics and the Dirac equation with flux, Phys. Rev. D 86, 084050 (2012)
- [40] M. Cariglia, A. Galajinsky, G.W. Gibbons and P.A. Horvathy, Cosmological aspects of the Eisenhart-Duval lift, Eur. Phys. J. C 78, 314 (2018)
- [41] A. Balcerzak and M. Lisaj, Spinor wave function of the Universe in non-minimally coupled varying constants cosmologies, Eur. Phys. J. C 83, 401 (2023)
- [42] C. Duval, G.W. Gibbons and P.A. Horvathy, Conformal and projective symmetries in Newtonian cosmology, J. Geom. Phys. 112, 197 (2017)

| [43]   | N. Kan, T. Aoyama, T. Hasegawa and K. Shiraishi, Eisenhart-Duval lift for minisu- perspace quantum cosmology, Phys. Rev. D 104, 086001 (2021)                             |
|--------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [44]   | A. Paliathanasis, Classical and quantum solutions in scalar field cosmology via the Eisenhart lift and linearization, Phys. Dark Univ. 44, 101466 (2024)                  |
| [45]   | A. Paliathanasis, Cosmological Constant from Equivalent Transformation in Quan- tum Cosmology, (2024) [arXiv:2405.20683]                                                  |
| [46]   | S. Dhasmana, A. Sen and Z.K. Silagadze, Equivalence of a harmonic oscillator to a free particle and Eisenhart lift, Annals of Physics 434, 168623 (2021)                  |
| [47]   | F.M. Mahomed and P.G.L. Leach, Symmetry Lie Algebras of nth Order Ordinary Differential Equations, J. Math. Anal. Appl. 151, 80 (1990)                                    |
| [48]   | E. Noether, Invariante Variationsprobleme Koniglich Gesellschaft der Wissenschaften Gottingen Nachrichten Mathematik-physik Klasse 2, 235-267 (1918)                      |
| [49]   | A. Halder, A. Paliathanasis and P.G.L. Leach, Noether's Theorem and Symmetry, Symmetry 10, 744 (2018)                                                                     |
| [50]   | W. Sarlet, F.M. Mahomed and P.G.L Leach, Symmetries of non-linear differential equations and linearisation, J. Phys A.: Math. Gen. 20, 277 (1987)                         |
| [51]   | L.M. Berkovich, Factorization and Transformations of Ordinary Differential Equa- tions, Saratov Univ. Publ., Saratov, (1989)                                              |
| [52]   | R. Campoamor-Stursberg and J. Guer´n, o Linearizing Systems of Second-Order ODEs via Symmetry Generators Spanning a Simple Subalgebra, Acta Appl. Math. 127, 105 (2013)   |
| [53]   | N.H. Ibragimov and S.V. Meleshko, Linearization of third-order ordinary differential equations by point and contact transformations, J. Math. Anal. Appl. 308, 266 (2005) |
| [54]   | S.V. Meleshko, On linearization of third-order ordinary differential equations, J. Phys. A.: Math. Gen. 39, 15135 (2006)                                                  |
| [55]   | L.M. Berkovich, Method of Factorization of ordinary differential operators and some of its applications, Appl. Anal. Disc. Math. 1, 122 (2007)                            |

- [56] W. Nakpim and S.V. Meleshko, Linearization of third-order ordinary differential equations by generalized Sundman transformations: X ′′′ + αX = 0, Commun. Nonlinear Sci. Numer. Simulat. 15, 1717 (2010)
- [57] A.Y. Al-Dweik, M.T. Mustafa, F.M. Mahomed and R.S. Alassar, Math. Meth. Appl. Sci. 41, 6955 (2018)
- [58] G.W. Bluman, Potential Symmetries and Linearization. In: Clarkson, P.A. (eds) Applications of Analytic and Geometric Methods to Nonlinear Differential Equations. NATO ASI Series, vol 413. Springer, Dordrecht (1993)
- [59] C. Wafo Soh, F.M. Mahomed, Linearization criteria for a system of second-order ordinary differential equations, Int. J. Non-Linear Mechanics 36, 671 (2001)
- [60] J. Patera, R. T. Sharp, P. Winternitz, and H. Zassenhaus, Invariants of real low dimension Lie algebras, J. Math. Phys. 17, 986 (1976)
- [61] S. Ali, F.M. Mahomed and A. Qadir, Linearizability criteria for systems of two second-order differential equations by complex methods, Nonlinear Dynamics 66, 77 (2011)
- [62] S. Ali, M. Safdar and A. Qadir, Linearization from complex Lie point transformations, J. Appl. Math. 2014, 793247 (2014)
- [63] M. Safdar, S. Ali and F.M. Mahomed, Linearization of systems of four second-order ordinary differential equations, Pranama - J. Phys. 77, 581 (2011)
- [64] A. Qadir, Geometric Linearization of Ordinary Differential Equations, Sigma 3, 1 (2007)
- [65] A.V. Aminova, Projective transformations and symmetries of differential equations, Sbornik Math. 185, 1711 (1995)
- [66] M. Tsamparlis, Linearization of Second-Order Non-Linear Ordinary Differential Equations: A Geometric Approach, Symmetry 15, 2082 (2023)
- [67] M.P. Ryan and L.C. Shepley, Homogeneous Relativistic Cosmologies, Princeton University Press, Princeton, (1975)

- [68] P.A. Terzis, N. Dimakis, T. Christodoulakis, A. Paliathanasis and M. Tsamparlis, Variational contact symmetries of constrained Lagrangians, J. Geom. Phys. 101, 52 (2016)
- [69] A. Paliathanasis, Solving Nonlinear Second-Order ODEs via the Eisenhart Lift and Linearization, Axioms 13, 331 (2024)
- [70] A. Paliathanasis, The Common Solution Space of General Relativity, (2024) [arXiv:2406.01998]
- [71] T. Christodoulakis, N. Dimakis and P.A. Terzis, Lie point and variational symmetries in minisuperspace Einstein gravity, J. Phys. A: Math. Theor. 47, 095202 (2014)
- [72] M. Tsamparlis and A. Paliathanasis, Two-dimensional dynamical systems which admit Lie and Noether symmetries, J. Phys. A: Math. Theor. 44, 175202 (2011)
- [73] A. Galajinsky and I. Masterov, Eisenhart lift for higher derivative systems, Phys. Lett. B 765, 86 (2017)
- [74] M. Tsamparlis, A. Paliathanasis and L. Karpathopoulos, Exact solutions of Bianchi I spacetimes which admit Conformal Killing vectors, Gen. Rel. Grav. 47, 15 (2015)
- [75] A. Paliathanasis and P.G.L. Leach, Symmetries and singularities of the Szekeres system, Phys. Lett. A 381, 1277 (2017)
- [76] A. Paliathanasis, A. Zampeli, T. Christodoulakis and M.T. Mustafa, Quantization of the Szekeres system, Class. Quantum Grav. 35, 125005 (2018)
- [77] P. Szekeres, A class of inhomogeneous cosmological models, Commun. Math. Phys. 41, 55 (1975)
- [78] A. Krasinski, Inhomogeneous cosmological models, Cambridge Univ. Press, Cambridge, UK, (2011)
- [79] J.D. Barrow and J. Stein-Schabes, Inhomogeneous cosmologies with cosmological constant, Phys. Lett. A 103, 315 (1984)
- [80] A. Zampeli and A. Paliathanasis, Quantization of inhomogeneous spacetimes with cosmological constant term, Class. Quantum Grav. 38, 165012 (2021)

- [81] T. Christodoulakis, N. Dimakis, P.A. Terzis, B. Vakili, E. Melas and Th. Grammenos, Minisuperspace canonical quantization of the Reissner-Nordstr¨m black hole o via conditional symmetries, Phys. Rev. D 89, 044031 (2014)
- [82] H. Reissner, Uber die Eigengravitation des elektrischen Feldes nach der Einsteinschen ¨ Theorie, Annals of Physics 50, 106 (1916)
- [83] G. Nordstr¨m, On the Energy of the Gravitational Field in Einstein's Theory, Proo ceedings Series B Physical Sciences 26, 1201 (1918)
- [84] A. Paliathanasis, Projective Collineations of Decomposable Spacetimes Generated by the Lie Point Symmetries of Geodesic Equations, Symmetry 13, 1018 (2021)
- [85] J.M. Lee and T.H. Parker, The Yamabe problem, Bull (NS) Am Math Soc 17, 37 (1987)
- [86] B. Abraham-Shrauner, K.S. Govinder and D.J. Arrigo, Type-II hidden symmetries of the linear 2D and 3D wave equations, J. Phys. A. Math. Gen. 39, 5739 (2006)
Andronikos Paliathanasis
Geometric Linearization for Constraint Hamiltonian Systems
This study investigates the geometric linearization of constraint Hamiltonian
systems using the Jacobi metric and the Eisenhart lift. We establish a
connection between linearization and maximally symmetric spacetimes, focusing
on the Noether symmetries admitted by the constraint Hamiltonian systems.
Specifically, for systems derived from the singular Lagrangian $$ L\left(
N,q^{k},\dot{q}^{k}\right) =\frac{1}{2N}g_{ij}\dot{q}^{i}\dot{q}^{j}-NV(q^{k}),
$$ where $N$ and $q^{i}$ are dependent variables and $\dim g_{ij}=n$, the
existence of $\frac{n\left( n+1\right) }{2}$ Noether symmetries is shown to be
equivalent to the linearization of the equations of motion. The application of
these results is demonstrated through various examples of special interest.
This approach opens new directions in the study of differential equation
linearization.