Next Article in Journal
Subclass of Analytic Functions Connected with Double Zeta Function
Next Article in Special Issue
Noether Symmetries and Conservation Laws in Non-Static Plane Symmetric Spacetime
Previous Article in Journal
A Dynamic Opposite Learning-Assisted Grey Wolf Optimizer
Previous Article in Special Issue
Noether Symmetries of Lemaitre-Tolman-Bondi Metric via Rif Tree Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quadratic First Integrals of Constrained Autonomous Conservative Dynamical Systems with Fixed Energy

by
Antonios Mitsopoulos
1,* and
Michael Tsamparlis
2,3,*
1
Department of Astronomy-Astrophysics-Mechanics, Faculty of Physics, University of Athens, Panepistemiopolis, 15783 Athens, Greece
2
NITheCS, National Institute for Theoretical and Computational Sciences, Pietermaritzburg 3201, KwaZulu-Natal, South Africa
3
TCCMMP, Theoretical and Computational Condensed Matter and Materials Physics Group, School of Chemistry and Physics, University of KwaZulu-Natal, Pietermaritzburg 3201, KwaZulu-Natal, South Africa
*
Authors to whom correspondence should be addressed.
Symmetry 2022, 14(9), 1870; https://doi.org/10.3390/sym14091870
Submission received: 1 August 2022 / Revised: 23 August 2022 / Accepted: 1 September 2022 / Published: 7 September 2022
(This article belongs to the Special Issue Noether and Space-Time Symmetries in Physics)

Abstract

:
We consider autonomous conservative dynamical systems which are constrained with the condition that the total energy of the system has a specified value. We prove a theorem which provides the quadratic first integrals (QFIs), time-dependent and autonomous, of these systems in terms of the symmetries (conformal Killing vectors and conformal Killing tensors) of the kinetic metric. It is proved that there are three types of QFIs and for each type we give explicit formulae for their computation. It is also shown that when the autonomous QFIs are considered, then we recover the known results of previous works. For a zero potential function, we have the case of constrained geodesics and obtain formulae to compute their QFIs. The theorem is applied in two cases. In the first case, we determine potentials which admit the second of the three types of QFIs. We recover a superintegrable potential of the Ermakov type and a new integrable potential whose trajectories for zero energy and zero QFI are circles. In the second case, we integrate the constrained geodesic equations for a family of two-dimensional conformally flat metrics.

1. Introduction

We consider autonomous conservative dynamical systems of the form
q ¨ a = Γ b c a ( q ) q ˙ b q ˙ c V , a ( q )
which are subjected to the fixed energy constraint
1 2 γ a b ( q ) q ˙ a q ˙ b + V ( q ) = E 0
where q a with a = 1 , 2 , , n are the generalized coordinates of the configuration space of the system, n is the dimension of the configuration space, a dot over a letter indicates derivation with respect to (wrt) the parameter t (time) along the trajectory q a ( t ) , a comma denotes partial derivative, Γ b c a ( q ) are the Riemannian connection coefficients defined by the kinetic metric γ a b ( q ) of the system, V ( q ) is the potential function of the system and E 0 is a fixed constant (i.e., the total energy –Hamiltonian– of the system). Moreover, Einstein’s summation convention is applied and the kinetic metric γ a b ( q ) is used for lowering/raising the tensorial indices.
This type of constrained systems is of particular interest in many areas of physics. For example, in astrophysics and celestial mechanics, autonomous conservative systems with a fixed energy are used in order to model the galactic motion. In such studies (see e.g., [1,2,3,4,5]), it is assumed that the galaxy is described by an axisymmetric autonomous potential and the authors look for a third first integral (FI) in addition to the well-known FIs of energy and angular momentum. The problem is still open; however, it seems that additional FIs may exist in certain regions of the energy domain. Similar considerations have been performed in non-relativistic Quantum Mechanics. Furthermore, in General Relativity, Equations (1) and (2) with V = 0 define constrained geodesic trajectories (timelike, spacelike, or null; depending upon the value of E 0 ) in a Riemannian spacetime which have numerous applications [6,7].
The most important role of the FIs is the assessment of the integrability of a dynamical system, constrained or not. Therefore, it is important that systematic methods are developed which will allow the determination of FIs. Concerning the case of the non-quantum constrained dynamical systems, there have been developed two types of such methods: (a) The method which uses the dynamical equations directly (see Section 1.1), and (b) The mini-superspace Lagrangian description (see Section 1.2). Each method follows a different approach, and its suitability depends on the particular application considered. The general description of each of these approaches is as follows.

1.1. The Dynamical Equations Method

There are three major approaches in this method which are the following:

1.1.1. The Geodesic Approach

Because the mth-order polynomials in velocities FIs of the geodesic equations in a Riemannian space are known (see e.g., [8,9]), one transfers the dynamical equations to the geodesic equations of another Riemannian space. One approach is to introduce two extra degrees of freedom in the configuration space. This is done by treating the time t as a generalized coordinate q n + 1 = t and by introducing an additional dimension q n + 2 via the fixed energy constraint. This approach has been originated by Eisenhart [6] who defined, in the ( n + 2 ) -dimensional expanded configuration space, the so-called Eisenhart metric whose geodesic equations produce both the dynamical equations and the energy constraint. A different approach is the Jacobi geometrization procedure (see e.g., [10,11]). According to this method, instead of introducing new generalized coordinates, one uses an energy-dependent reparameterization of the parameter of the dynamical equations to define a new metric, the Jacobi metric, whose timelike geodesics produce the dynamical equations and the constraint. The new metric is conformally related to the kinetic metric of the system with conformal factor the scalar quantity E 0 V ( q ) .
The approach of the Jacobi metric has been applied in [12,13,14,15] in order to assess the integrability of two-dimensional (2d) autonomous conservative systems constrained at fixed and arbitrary energy. In particular, in [13,14], the integrable and superintegrable 2d Newtonian potentials that admit additional autonomous cubic and quartic FIs have been determined. However, certain drawbacks exist concerning the Jacobi metric when one increases the number of degrees of freedom, or the order of the FIs. This is due to the fact that: (i) The conformal Killing vectors (CKVs) and the conformal Killing tensors (CKTs) of a general metric are not easy to compute, and (ii) The Jacobi metric assumes a reparameterization from the ‘physical time’ of the system to the so-called ‘Jacobi time’ which affects the physical interpretation of the results.

1.1.2. The Lie Symmetry Approach

This is a systematic and interesting approach which has been developed in [16,17,18] and is given as follows:
One considers the unconstrained dynamical Equation (1) and looks for Lie point symmetries generated by the infinitesimal point transformations:
q ¯ a ( t ¯ ) = q a ( t ) + η a ( q ) ε , t ¯ = t + ξ ( t ) ε
where 0 < ε 1 , the generator ξ ( t ) = 2 ϕ ( q ( t ) ) d t ξ ˙ = 2 ϕ ( q ) and the function ϕ is evaluated along trajectories. These Lie point symmetries are referred to as trajectory collineations (TCs). It is proved that the generators η a ( q ) of the TCs are projective collineations (PCs) of the kinetic metric γ a b with projective function ϕ ( q ) which is related to the potential V with the condition:
4 ϕ V , a + γ a b V ; b c η c γ b c V , b η a ; c = 0 L η V , a + 4 ϕ V , a = 0
where L η denotes the Lie derivative with respect to η a (see also [19]). Furthermore, it is shown that the Lie derivative of a known FI (e.g., the total energy) with respect to a TC is also a FI. This last result allows for a systematic evaluation of FIs once the TCs have been determined. These FIs have been called related FIs (RFIs). It is also shown that the quadratic RFIs share the same structure constants with the Lie algebra of TCs.
In [18], the above approach was extended to the constrained autonomous conservative dynamical systems in the sense that one requires that the point transformations (3) are constrained by the further requirement that the energy has a fixed value E 0 . These new transformations map constrained solutions into constrained solutions and, hence, they are Lie point symmetries of the constrained system. These symmetries have been called natural trajectory collineations (NTCs). It is proved that the generators η a ( q ) of the NTCs are CKVs of the kinetic metric γ a b with conformal factor ϕ + c , where c is a constant, which are constrained by the new ‘selection rule’
V , a η a + 2 ( ϕ c ) ( V E 0 ) = 0 .
It has been also shown that the formulation of the RFIs holds the same for the NTCs.

1.1.3. The Direct Approach

In this approach, one assumes quadratic FIs (QFIs) of the functional form
I = K a b ( t , q ) q ˙ a q ˙ b + K a ( t , q ) q ˙ a + K ( t , q )
where the coefficients are symmetric tensors depending on q a and t. The linear FIs (LFIs) are also included for K a b = 0 .
The condition that (6) is a FI of the constrained system (1) and (2) is
d I d t = ψ ( t , q ) + X c ( t , q ) q ˙ c γ a b q ˙ a q ˙ b + 2 ( V E 0 )
where ψ ( t , q ) is an arbitrary scalar and X a ( t , q ) is an arbitrary vector. It is to be noted that condition (7) is due to Hilbert’s zero-theorem (see Chapter IV, par. 8, pp. 166–167 in [20]). According to this theorem, since d I d t is a third order polynomial in the velocities (after replacing q ¨ a from (1)) and the constraint (2) is quadratic, there exists a first order ( 3 2 = 1 ) polynomial ψ ( t , q ) + X c ( t , q ) q ˙ c such that the quantity d I d t is expressed as in equation (7). The term ψ ( t , q ) + X c ( t , q ) q ˙ c plays the role of a ‘Lagrangian multiplier’. Using the dynamical equations (1) to replace the terms q ¨ a whenever they appear, condition (7) leads to a polynomial equation in the velocities q ˙ a .
At this point, there are two approaches. Either one looks for QFIs which hold for all values of the velocities on the constraint surface (we call them exact QFIs), or one looks for QFIs which are non-local, that is, they are valid only on trajectories specified by certain velocities. In the first approach, one sets all the coefficients of the powers of the velocities equal to zero and obtains a system of partial differential equations (PDEs) of the unknown quantities K a b , K a , K and the dynamical quantities V ( q ) , E 0 . The solution of this system of PDEs provides the exact QFIs I. In the second approach, one does the same for all powers of the velocities which are of degree equal or higher than two, that is, the terms of zeroth and first order are excluded. In this case, one finds the conditional QFIs which are non-local QFIs, that is, they contain quantities within integrals. This is the type of FIs given, e.g., in [21,22,23].

1.2. The Mini-Superspace Lagrangian Method

In this method, one introduces an extra dimension in the configuration space, the lapse function, and defines the mini-superspace Lagrangian of the system. The Euler–Lagrange (E-L) equations of this Lagrangian produce the dynamical equations and the equation of the constraint for E 0 = 0 . This method is a particular case of a more general procedure applied on degenerate Lagrangian constrained systems. These are constrained systems whose Hessian vanishes, i.e., the quantity det 2 L q ˙ α q ˙ β = 0 where q α denotes the generalized coordinates. In this case, a Hamiltonian cannot be defined because the Legendre transformation has no inverse. Since the Hamiltonian is vital to Quantum Mechanics, the Dirac–Bergmann prescription has been developed [24,25,26,27] to deal with this type of systems. In general terms, this formalism has as follows. The constraints are equations of the form ϕ A ( q , p ) = 0 with A = 1 , 2 , , r , where r measures the degeneracy of the singular Lagrangian (i.e., the degeneracy of the Hessian) and p α = L q ˙ α are the conjugate momenta. Relations ϕ A ( q , p ) = 0 are used to eliminate the r degenerate coordinates. These relations are called primary constraints and their time derivatives, secondary constraints.
Using the primary constraints, one splits the coordinates in two groups. The first group contains n r coordinates which define a non-vanishing Hessian, and the second group contains the remaining r coordinates. This approach leads to conditional FIs which are non-local, that is, they are expressed in terms of integrals and are generated by the CKVs and the CKTs instead of the Killing vectors (KVs) and the Killing tensors (KTs) which generate the exact QFIs. The conditional FIs can be used in the same way as the exact FIs; for example, in order to integrate the geodesic equations. A detailed and clear discussion of this approach with applications can be found in [22].
In the present work, we determine the QFIs (6) for the constrained dynamical systems (1) and (2) by using the direct method and the geometry generated by the kinetic metric of the system.
The structure of the paper is as follows.
In Section 2, we derive the main result of the present paper which is Theorem 1. It is shown that there are at most three types of exact QFIs associated with an autonomous conservative dynamical system constrained on a fixed energy level. Explicit formulae are given which provide these FIs in terms of the geometric symmetries of the kinetic metric. We note that these QFIs can be autonomous or time-dependent. In Section 2.1, we consider the autonomous LFIs and QFIs derived from Theorem 1, and we recover the FIs found in [18]. In Section 3, we recall the basic facts concerning the CKTs of order two, which are necessary in the computation of the QFIs of Theorem 1. In Section 4, we apply Theorem 1 to find 2D Newtonian potentials which admit a QFI of the second type I ( ) 2 . We consider the cases where the vector L a ( q ) is either a HV or a special CKV (SCKV). In the first case, we obtain a superintegrable potential of the Ermakov type and compute the trajectory for specific values of the constants. In the second case, we find a new class of integrable potentials and show that the orbits for E 0 = 0 and I ( 0 ) 2 = 0 are circles. In Section 5, we specialize Theorem 1 to the case of constrained geodesics and we collect our results in Theorem 2. Applying Theorem 2, we recover all previous results in the literature (see e.g., [8,28,29]) for an arbitrary energy level E 0 , and in Propositions 1 and 2 we determine the QFIs of null ( E 0 = 0 ) and non-null ( E 0 0 ) geodesics, respectively. In Section 6, we consider the non-null geodesics of a conformally flat metric and recover the results of [22] which were found using the rather complicated Dirac–Bergmann prescription. Finally, in Section 7, we draw our conclusions.

2. QFIs of Constrained Autonomous Conservative Dynamical Systems with Fixed Energy

We consider autonomous conservative dynamical systems defined by Equations (1) and (2), and we look for QFIs of the general form (6) defined by the condition (7).
Using the dynamical Equation (1) to replace the quantities q ¨ a whenever they appear, condition (7) gives
0 = K ( a b ; c ) X ( a γ b c ) q ˙ a q ˙ b q ˙ c + K a b , t + K ( a ; b ) ψ γ a b q ˙ a q ˙ b + + K a , t + K , a 2 K a b V , b 2 ( V E 0 ) X a q ˙ a + K , t K a V , a 2 ( V E 0 ) ψ .
Assuming that there are no constraints on the velocities q ˙ a , we demand that the polynomial Equation (8) is satisfied for all values of q ˙ a ; therefore, the coefficient of each power of q ˙ a must vanish. This leads to the system of PDEs:
K ( a b ; c ) = X ( a γ b c )
K ( a ; b ) = ψ γ a b K a b , t
K , a = 2 K a b V , b + 2 ( V E 0 ) X a K a , t
K , t = K a V , a + 2 ( V E 0 ) ψ
where a semicolon denotes Riemannian covariant derivative.
Equation (9) implies that K a b is a CKT of order two of the kinetic metric γ a b with associated vector X a . Contracting (9) with γ b c , we find
X a = 1 n + 2 K b b ; a + 2 K b a ; b .
Contracting Equation (10) with γ a b , we obtain
ψ = 1 n K a ; a + K a a , t .
The system of PDEs (9)–(12) must be supplemented with the integrability conditions K , [ a t ] = 0 and K ; [ a b ] = 0 for the scalar K. We have:
K a , t t 2 K a b , t V , b + K b V , b , a + 2 ( V E 0 ) ψ , a X a , t = 0
2 K [ a | c | V , c ; b ] K [ a ; b ] , t + 2 ( V E 0 ) X [ a ; b ] = 0
where round (square) brackets indicate symmetrization (antisymmetrization) of the enclosed indices, and indices enclosed between vertical lines are overlooked by (anti-)symmetrization symbols.
Finally, the system of PDEs which we have to solve consists of Equations (9)–(12) and (15)–(16), where the quantities X a and ψ are given by (13) and (14), respectively.
We state the solution of the above system of PDEs in Theorem 1. The proof is given in the Appendix A.
Theorem 1.
The independent QFIs of the autonomous conservative dynamical system (1) subject to the fixed energy constraint (2) are the following:
Integral 1.
I ( ) 1 = k = 1 t 2 k 2 k L ( 2 k 1 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + k = 1 t 2 k 1 L ( 2 k 1 ) a q ˙ a + k = 1 t 2 k 2 k L ( 2 k 1 ) a V , a + G ( q )
where C ( 0 ) a b ( q ) and L ( 2 k 1 ) ( a ; b ) ( q ) for k = 1 , 2 , , are CKTs with associated vectors X ( 0 ) a ( q ) and Y ( 2 k 1 ) a ( q ) , respectively, while the vectors L ( 2 k 1 ) a ( q ) and the function G ( q ) satisfy the conditions:
L ( 2 1 ) b V , b , a = 2 L ( 2 1 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 2 1 ) a
L ( 2 k 1 ) b V , b , a = 2 L ( 2 k 1 ) ( a ; b ) V , b 2 k ( 2 k + 1 ) L ( 2 k + 1 ) a 2 ( V E 0 ) Y ( 2 k 1 ) a , k = 1 , 2 , , 1
G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a L ( 1 ) a ( > 0 ) .
Integral 2.
I ( ) 2 = k = 0 t 2 k + 1 2 k + 1 L ( 2 k ) ( a ; b ) q ˙ a q ˙ b + t 2 k L ( 2 k ) a q ˙ a + t 2 k + 1 2 k + 1 L ( 2 k ) a V , a
where L ( 2 k ) ( a ; b ) ( q ) for k = 0 , 1 , , are CKTs with associated vectors Y ( 2 k ) a ( q ) and the vectors L ( 2 k ) a ( q ) satisfy the conditions:
L ( 2 ) b V , b , a = 2 L ( 2 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 2 ) a
L ( 2 k ) b V , b , a = 2 L ( 2 k ) ( a ; b ) V , b 2 ( k + 1 ) ( 2 k + 1 ) L ( 2 k + 2 ) a 2 ( V E 0 ) Y ( 2 k ) a , k = 0 , 1 , , 1 .
Integral 3.
I ( e ) = e λ t L ( a ; b ) q ˙ a q ˙ b + λ L a q ˙ a + L a V , a
where λ 0 is an arbitrary constant and L ( a ; b ) ( q ) is a reducible CKT with associated vector Y a ( q ) such that
L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a 2 ( V E 0 ) Y a .
Notation: In I ( ) α with α = 1 , 2 , the index ( ) indicates the degree of the time-dependence, whereas the index α is used in order to distinguish between the two different types of independent QFIs. Moreover, in Equation (20), the quantity L ( 1 ) a ( > 0 ) indicates that the vector L ( 1 ) a exists only when m > 0 .
We note that in the case that the involved second order CKTs of Theorem 1 are KTs, the associated vectors X ( 0 ) a and Y ( N ) a vanish and Theorem 1 reduces to Theorem 3 of [29] for Q a = V , a . Then, the associated independent QFIs hold for an arbitrary energy level E 0 .

2.1. The Autonomous LFIs/QFIs of Theorem 1

Theorem 1 contains only two autonomous FIs. These are the following:
a.
The QFI
J 1 = C a b ( q ) q ˙ a q ˙ b + G ( q )
where C a b is a second order CKT with associated vector X a ( q ) such that G , a = 2 C a b V , b + 2 ( V E 0 ) X a .
The QFI (26) is derived from the FI I ( ) 1 given in (17) for time-dependence = 0 . To simplify the notation, we have set C ( 0 ) a b = C a b and X ( 0 ) a b = X a b .
b.
The LFI
J 2 = L a ( q ) q ˙ a
where L a is a CKV with conformal factor ψ ( q ) such that L a V , a = 2 ( V E 0 ) ψ .
The LFI (27) is derived from the FI I ( ) 2 given in (21) for time-dependence = 0 if we assume that L ( 0 ) a L a ( q ) is a CKV with conformal factor ψ ( q ) . Indeed, we have
L ( a ; b ) = ψ γ a b Y ( 0 ) a = ψ , a
and the condition
L b V , b , a = 2 ψ γ a b V , b 2 ( V E 0 ) ψ , a L a V , a + 2 ( V E 0 ) ψ = c = c o n s t .
Therefore, the associated QFI becomes
I ( 0 ) 2 = t ψ γ a b q ˙ a q ˙ b + L a V , a + L a q ˙ a = t 2 ( V E 0 ) ψ + L a V , a + L a q ˙ a = L a q ˙ a + c t
which for c = 0 gives the autonomous LFI (27).
The autonomous FIs (26) and (27) are those found in Section VI of [18]. In particular: a. The QFI (26) coincides with Equation (6.5) of [18] and the additional conditions are Equations (6.6) and (6.7); b. The LFI (27) coincides with Equation (6.1) of [18] and the additional conditions are Equations (6.3) and (6.4). Moreover, the QFI given in Equation (4.7) of [30] is derived from the QFI (26) for G = E 0 = 0 . We note that the authors in [30] determine this QFI by using the contact symmetries of the constrained system in the mini-superspace Lagrangian formalism.

3. Conformal Killing Tensors (CKTs) of Order Two

In this section, we recall some basic results concerning the CKTs which will be used in applications of Theorem 1 to be considered in the next sections.
A second order CKT [31] in an m-dimensional Riemannian manifold with local coordinates x a and metric g a b ( x ) is a symmetric tensor U a b such that
U ( a b ; c ) = u ( a g b c ) U { a b ; c } = u { a g b c }
where u a is the vector associated to the CKT U a b and curly brackets denote cyclic permutation of the enclosed indices.
By contracting (28) with g a b , we find the vector
u a = 1 m + 2 U ; a + 2 U b a ; b
where U U a a is the trace of U a b and m is the dimension of the manifold. The following terminology and results apply to CKTs:
a.
If u a = 0 , U a b is a second order KT or an improper CKT.
b.
If u a 0 , U a b is called a proper CKT.
c.
If u a is a KV, U a b is called a homothetic KT (HKT) [32].
d.
If the trace U = 0 , then u a = 2 m + 2 U b a ; b and U a b is called a trace-free CKT.
e.
If u a is a gradient (i.e., u a = u , a where u = u ( x ) is a scalar), then U a b is called a CKT of gradient type.
f.
If f is an arbitrary function, then U ¯ a b = U a b + f g a b is a CKT with associated vector u ¯ a = u a + f ; a .
g.
If U a b is a CKT of gradient type, then C a b = U a b u g a b is a second order KT.
We can construct new CKTs as follows:
(i)
If f is an arbitrary function, then f g a b is a gradient CKT with associated vector f , a .
(ii)
If X a and Y a are CKVs with conformal factors, respectively, ψ X and ψ Y , then the symmetrized tensor product X ( a Y b ) is a CKT with associated vector ψ X Y a + ψ Y X a .
(iii)
If T a b and T ¯ a b are CKTs with associated vectors, respectively, V a and W a , then the linear combination λ T a b + μ T ¯ a b , where λ and μ are arbitrary constants, is a CKT with associated vector λ V a + μ W a .
From (i)–(iii), we get the following general result [31,33,34].
If an m-dimensional manifold admits M CKVs X K a with conformal factors ψ K where K = 1 , 2 , , M , then
U a b = f g a b + c K L X K ( a X | L | b )
is a CKT with associated vector
u a = f , a + c K L ψ K X L a + ψ L X K a
where f is an arbitrary function, c K L are arbitrary constants and the summation is over the inequality 1 K L M . In flat spaces, all second order CKTs are of the form (30).

4. Example 1: Constrained Orbits of Potentials V ( x , y ) in E 2 That Admit FIs of the Form I ( 0 ) 2 at Fixed Energy E 0 = 0

In this case, the generalized coordinates q a = ( x , y ) , the kinetic metric γ a b = δ a b = d i a g ( 1 , 1 ) , the fixed energy level E 0 = 0 and the constraint (2) becomes
γ a b q ˙ a q ˙ b + 2 V = 0 x ˙ 2 + y ˙ 2 + 2 V = 0 .
We note that we can always insert the value of the constant E 0 into the potential; therefore, the case E 0 = 0 has no effect on the generality of our discussion.
For time-dependence = 0 and fixed energy E 0 = 0 , the FI (21) becomes (To simplify the notation, we set L ( 0 ) a ( q ) L a ( q ) )
I ( 0 ) 2 = t L ( a ; b ) q ˙ a q ˙ b + L a q ˙ a + t L a V , a
where L ( a ; b ) is a second order CKT of δ a b with associated vector Y a such that
L b V , b , a = 2 L ( a ; b ) V , b 2 V Y a .
We consider various cases concerning the vector L a .

4.1. L a = ( x , y ) Is the Homothetic Vector (HV)

If the vector L a = ( x , y ) is the HV of E 2 , then L ( a ; b ) = δ a b and Y a = 0 .
Using the constraint (32), the QFI (33) reduces to the LFI
I ( 0 ) 2 I 2 = t ( x ˙ 2 + y ˙ 2 ) = 2 V + x x ˙ + y y ˙ + t ( x V , x + y V , y ) = x x ˙ + y y ˙ + c t = r r ˙ + c t
and the condition (34) becomes
x V , x + y V , y + 2 V = c
where we have set I ( 0 ) 2 I 2 , r 2 = x 2 + y 2 and c is an arbitrary constant.
Solving the PDE (36), we find the potential
V ( x , y ) = F y x r 2 + c 2
where F is an arbitrary smooth function of its argument.
The potential (37) admits also the Ermakov QFI (see e.g., [35,36])
I 1 = ( x y ˙ y x ˙ ) 2 + 2 F ;
therefore, it is superintegrable. Using the FIs (35) and (38), the constraint (32) becomes:
x ˙ 2 + y ˙ 2 + 2 F r 2 + c = 0 ( x x ˙ + y y ˙ ) 2 + ( x y ˙ y x ˙ ) 2 + 2 F + c r 2 = 0 ( I 2 c t ) 2 + I 1 + c r 2 = 0
where I 1 and I 2 are arbitrary constants.
Integrating the LFI (35), we find
r 2 · = 2 ( I 2 c t ) r 2 = c t 2 + 2 I 2 t + c 1 r ( t ) = c t 2 + 2 I 2 t + c 1
where c 1 is an arbitrary constant.
Replacing (40) in (39), we get the condition
I 1 = I 2 2 c c 1 .
Using polar coordinates x = r cos θ and y = r sin θ , and conditions (40) and (41), the Ermakov QFI (38) is written
r 2 θ ˙ = I 2 2 c c 1 2 F θ ˙ = I 2 2 c c 1 2 F c t 2 + 2 I 2 t + c 1 .
Because F is a function of y x = tan θ , i.e., F ( tan θ ) , Equation (42) cannot be integrated to give θ ( t ) .
In order to integrate (42), we assume c = 0 and F = k = c o n s t . Then, the potential (37) reduces to the well-known Newton–Cotes potential (see e.g., [37,38])
V = k r 2
and Equations (40) and (42) give the solution:
r ( t ) = 2 I 2 t + c 1 , θ ( t ) = 1 4 k 2 I 2 2 ln 2 I 2 t + c 1 + θ 0
which implies the orbit (with energy E 0 = 0 !)
r = A e B θ
where θ 0 is an integration constant, B 1 1 2 k I 2 2 and A e B θ 0 . It has been checked that the solution (44) satisfies the polar Euler–Lagrange equation r ¨ = r θ ˙ 2 + 2 k r 3 .
The orbit (45) is a logarithmic spiral (or miraculous spiral) which has the following properties: (a) For B > 0 , the size of the spiral increases outward as θ increases. (b) For B < 0 , the size of the spiral decreases inward as θ increases.
We note that for the potential (43), Equation (42) produces the LFI of the angular momentum (as expected for a central potential).

4.2. L a Is a Special CKV (SCKV)

It is well-known that E 2 admits two SCKVs:
(i)
B ( 1 ) a = x 2 y 2 2 x y with conformal factor ψ 1 = x .
(ii)
B ( 2 ) a = x y y 2 x 2 2 with conformal factor ψ 2 = y .
We assume that the vector L a = x 2 y 2 2 x y is the SCKV B ( 1 ) a of E 2 with conformal factor ψ = x . Then, L ( a ; b ) = x δ a b and Y a = ψ , a = 1 , 0 .
Using the constraint (32), we find that the QFI (33) reduces to the LFI
I ( 0 ) 2 = t x x ˙ 2 + y ˙ 2 = 2 V + x 2 y 2 2 x ˙ + x y y ˙ + t x 2 y 2 2 V , x + x y V , y = x 2 y 2 2 x ˙ + x y y ˙ + c t
and the associated condition (34) becomes
x 2 y 2 2 V , x + x y V , y + 2 x V = c
where c is an arbitrary constant.
For c = 0 , the PDE (47) gives the potential
V = M y r 2 r 4
where r 2 = x 2 + y 2 and M is an arbitrary smooth function of its argument, the LFI (46) becomes
I ( 0 ) 2 = x 2 y 2 2 x ˙ + x y y ˙
and the constraint (32) is written
x ˙ 2 + y ˙ 2 + 2 M r 4 = 0 .
We can use the LFI (49) and the zero energy constraint (50) in order to find an orbit with E 0 = 0 for the integrable potential (48).
First, we assume the additional requirement I ( 0 ) 2 = 0 . Then, the LFI (49) becomes
x 2 y 2 2 x ˙ + x y y ˙ = 0 .
Using polar coordinates, Equation (51) gives:
cos θ r ˙ + r sin θ θ ˙ = 0 ln r · = tan θ θ ˙ ln r = tan θ d θ r = c 1 cos θ
where c 1 0 is an arbitrary constant. The orbits (52) are circles with center ( c 1 2 , 0 ) and radius | c 1 | 2 . Indeed, we have
r = c 1 cos θ r 2 = c 1 r cos θ x 2 + y 2 = c 1 x x c 1 2 2 + y 2 = c 1 2 4 .
In order to compute the θ ( t ) , we use the constraint (50) which is written
r ˙ 2 + r 2 θ ˙ 2 + 2 M sin θ r r 4 = 0 .
Replacing the orbit (52) in (53), we obtain
θ ˙ 2 = 2 M tan θ c 1 c 1 6 cos 4 θ t t 0 = c 1 3 cos 2 θ d θ 2 M tan θ c 1
where t 0 is an integration constant.
Therefore, at the fixed energy E 0 = 0 , we have found the new class of integrable potentials (48) and we have shown that for I ( 0 ) 2 = 0 the orbits of these potentials are circles.

5. The LFIs/QFIs of the Constrained Geodesic Equations

Another important area in which we apply Theorem 1 is the determination of the LFIs and QFIs of the constrained geodesic equations in an n-dimensional Riemannian manifold with metric γ a b ( q ) . In this case, V = 0 and Theorem 1 takes the following form.
Theorem 2.
The independent QFIs of the geodesic equations
q ¨ a + Γ b c a ( q ) q ˙ b q ˙ c = 0
subject to the quadratic constraint
γ a b ( q ) q ˙ a q ˙ b = 2 E 0
where Γ b c a are the Riemannian connection coefficients defined by the metric γ a b ( q ) and E 0 is an arbitrary fixed constant, are the following:
Integral 1.
I ( ) 1 = k = 1 t 2 k 2 k L ( 2 k 1 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + k = 1 t 2 k 1 L ( 2 k 1 ) a q ˙ a + G ( q )
where C ( 0 ) a b ( q ) and L ( 2 k 1 ) ( a ; b ) ( q ) for k = 1 , 2 , , are CKTs with associated vectors X ( 0 ) a ( q ) and Y ( 2 k 1 ) a ( q ) , respectively, while the vectors L ( 2 k 1 ) a ( q ) and the function G ( q ) satisfy the conditions:
E 0 Y ( 2 1 ) a = 0
E 0 Y ( 2 k 1 ) a = k ( 2 k + 1 ) L ( 2 k + 1 ) a , k = 1 , 2 , , 1
G , a = 2 E 0 X ( 0 ) a L ( 1 ) a ( > 0 ) .
Integral 2.
I ( ) 2 = k = 0 t 2 k + 1 2 k + 1 L ( 2 k ) ( a ; b ) q ˙ a q ˙ b + t 2 k L ( 2 k ) a q ˙ a
where L ( 2 k ) ( a ; b ) ( q ) for k = 0 , 1 , , are CKTs with associated vectors Y ( 2 k ) a ( q ) and the vectors L ( 2 k ) a ( q ) satisfy the conditions:
E 0 Y ( 2 ) a = 0
E 0 Y ( 2 k ) a = ( k + 1 ) ( 2 k + 1 ) L ( 2 k + 2 ) a , k = 0 , 1 , , 1 .
Integral 3.
I ( e ) = e λ t L ( a ; b ) q ˙ a q ˙ b + λ L a q ˙ a
where λ 0 is an arbitrary constant and L ( a ; b ) ( q ) is a reducible CKT with associated vector Y a ( q ) such that
L a = 2 E 0 λ 2 Y a .
When there is no constraint of the form (56), that is, the value of E 0 is arbitrary, the CKTs of Theorem 2 reduce to KTs (i.e., the associated vectors X ( 0 ) a = Y ( k ) a = 0 ) and, as expected, the QFIs of Theorem 2 produce the well-known ones (see [8], Table IV of [28], and Section X, par. A in [29]).

5.1. The QFIs for Null Constrained Geodesic Equations: Case E 0 = 0

Applying Theorem 2 for E 0 = 0 , we find the following result.
Proposition 1.
The independent QFIs of the null geodesic equations
q ¨ a + Γ b c a ( q ) q ˙ b q ˙ c = 0 ,
that is, subject to the constraint
γ a b ( q ) q ˙ a q ˙ b = 0
where Γ b c a are the Riemannian connection coefficients defined by the metric γ a b ( q ) , are the following:
I 1 = C a b q ˙ a q ˙ b , I 2 = t 2 2 G ; a b q ˙ a q ˙ b t G , a q ˙ a + G ( q ) , I 3 = t L ( a ; b ) q ˙ a q ˙ b + L a q ˙ a
where the quantities C a b , G ; a b and L ( a ; b ) are second order CKTs of γ a b .
Comparing the QFIs (68) for null geodesics with the three QFIs of Table 4 in [28] for unconstrained geodesics, we may be tempted to say that they are the same. However, such a claim is not correct! Observe that in the former case the involved second order symmetric tensors are CKTs; whereas, in the latter case, they are KTs and the corresponding QFIs hold for an arbitrary E 0 . This fundamental difference lies in the fact that the null geodesics are subject to the constraint (67).

5.2. The QFIs of Non-Null (Spacelike or Timelike) Constrained Geodesic Equations: Case E 0 0

Applying Theorem 2 for E 0 0 , we have the following result.
Proposition 2.
The independent QFIs of the geodesic equations
q ¨ a + Γ b c a ( q ) q ˙ b q ˙ c = 0
subject to the quadratic constraint
γ a b ( q ) q ˙ a q ˙ b = 2 E 0 0
where Γ b c a are the Riemannian connection coefficients defined by the metric γ a b ( q ) and E 0 0 is an arbitrary fixed constant, are the following:
Integral 1.
I ( ) 1 = k = 1 t 2 k 2 k L ( 2 k 1 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + k = 1 t 2 k 1 L ( 2 k 1 ) a q ˙ a + G ( q )
where L ( 2 1 ) ( a ; b ) ( q ) is a KT, C ( 0 ) a b ( q ) and L ( 2 k 1 ) ( a ; b ) ( q ) for k = 1 , 2 , , 1 are CKTs with associated vectors X ( 0 ) a ( q ) and Y ( 2 k 1 ) a = k ( 2 k + 1 ) E 0 L ( 2 k + 1 ) a , respectively, and the function G ( q ) is such that
G , a = 2 E 0 X ( 0 ) a L ( 1 ) a ( > 0 ) .
Integral 2.
I ( ) 2 = k = 0 t 2 k + 1 2 k + 1 L ( 2 k ) ( a ; b ) q ˙ a q ˙ b + t 2 k L ( 2 k ) a q ˙ a
where L ( 2 ) ( a ; b ) is a KT and L ( 2 k ) ( a ; b ) ( q ) for k = 0 , 1 , , 1 are CKTs with associated vectors Y ( 2 k ) a = ( k + 1 ) ( 2 k + 1 ) E 0 L ( 2 k + 2 ) a .
Integral 3.
I ( e ) = e λ t L ( a ; b ) q ˙ a q ˙ b + λ L a q ˙ a
where λ 0 is an arbitrary constant and L ( a ; b ) ( q ) is a reducible CKT with associated vector Y a = λ 2 2 E 0 L a .
In the following section, we apply the results of Section 5 to a number of cases considered in [22], and we show that we recover the answers without using the complex arguments of the mini-superspace Lagrangian method and the conditional FIs.

6. Example 2: The Non-Null Constrained Geodesic Equations of the Metric γ ab = f ( x , y ) 0 1 1 0

Consider the constrained dynamical system of two degrees of freedom q a = ( x , y ) , which is described by the dynamical equations:
x ¨ + f , x f x ˙ 2 = 0
y ¨ + f , y f y ˙ 2 = 0
f ( x , y ) x ˙ y ˙ = E 0 .
Equations (75) and (76) are subjected to the quadratic constraint (77), where f ( x , y ) is an arbitrary smooth function and E 0 is a non-zero fixed constant.
From the constraint Equation (77), we read the 2d metric
γ a b = f ( x , y ) 0 1 1 0 .
We calculate that the Riemann tensor is
R 1212 = f , x y f , x f , y f = R 2 γ 1212
where γ a b c d γ a c γ b d γ a d γ b c and γ 1212 = f 2 . However, the metric (78) is not in general a metric of constant curvature because the Ricci scalar R is not constant.
The non-vanishing Riemannian connection coefficients for the metric (78) are Γ 11 1 = f , x f and Γ 22 2 = f , y f . We note then that the dynamical Equations (75) and (76) are the geodesic Equation (69) of the metric (78) provided t is the affine parameter. The constraint Equation (77) is of the required form (70).
The metric (78) has been discussed for various cases of the function f ( x , y ) in [22], using the rather complicated method of Dirac constraints [26,27] and the associated conditional FIs. Our purpose here is to solve these cases, using the FIs of Proposition 2.
We shall need the CKVs of the metric (78) in order to compute from them the KVs and the CKTs of this metric. For the metric (78), we find the two-parameter family of CKVs
B a = f ( x , y ) F 1 ( y ) F 2 ( x ) B = B a a = F 2 ( x ) x + F 1 ( y ) y
whose conformal factor
ψ ( x , y ) = F 2 f , x + F 1 f , y + f d F 1 d y + d F 2 d x 2 f
where F 1 ( y ) and F 2 ( x ) are arbitrary smooth functions.
Using the general formulae (30) and (31), we construct for the metric (78) the second order CKT
C a b = f 2 ( x , y ) A 1 ( y ) 0 0 A 2 ( x )
with associated vector
X a = f , y A 1 ( y ) + f 2 d A 1 d y f , x A 2 ( x ) + f 2 d A 2 d x
where A 1 ( y ) and A 2 ( x ) are arbitrary smooth functions.
To find the KVs, we require ψ = 0 . Then, condition (81) implies that
F 2 f , x + F 1 f , y + f d F 1 d y + d F 2 d x = 0 .

6.1. The Case of a Space of Constant Curvature: LFIs

We require the metric (78) to be a metric of constant curvature. Since the metric (78) is 2d and satisfies the condition (79), it must also satisfy the additional requirement that R = c o n s t .
Without loss of generality, we set R = 4 k where k 0 is an arbitrary non-zero constant. Replacing this R in (79), we get the non-linear PDE:
f , x y f , x f , y f 2 f 2 k = 0
whose solution is f ( x , y ) = k ( x + y ) 2 . Therefore, the metric of constant curvature is
γ a b = k ( x + y ) 2 0 1 1 0 .
The constrained system (75)–(77) becomes:
x ¨ 2 x + y x ˙ 2 = 0
y ¨ 2 x + y y ˙ 2 = 0
k ( x + y ) 2 x ˙ y ˙ = E 0
A space of constant curvature has three KVs. Therefore, the system of the constrained geodesic Equations (87)–(89) admits three independent LFIs and, as a result, it is superintegrable. To find the KVs, we use Equation (84) which in this case becomes
x d F 1 d y + y d F 2 d x + y d F 1 d y 2 F 1 + x d F 2 d x 2 F 2 = 0 .
From the leading terms x d F 1 d y and y d F 2 d x in (90), we deduce that the function F 1 ( y ) must be of the form
F 1 ( y ) = c 1 y 2 + c 2 y + c 3
where c 1 , c 2 , c 3 are arbitrary constants.
Replacing (91) in (90), we find
F 2 ( x ) = c 1 x 2 + c 2 x c 3 .
Substituting the functions (91) and (92) in (80), we find for the metric (86) the following three KVs (one for each constant c i , i = 1 , 2 , 3 ):
L ( 1 ) a = k ( x + y ) 2 y 2 x 2 , L ( 2 ) a = k ( x + y ) 2 y x , L ( 3 ) a = k ( x + y ) 2 1 1 .
Replacing the KVs (93) in the FI (73) for = 0 , we find the following three LFIs:
I 1 = k ( x + y ) 2 y 2 x ˙ x 2 y ˙ , I 2 = k ( x + y ) 2 y x ˙ + x y ˙ , I 3 = k ( x + y ) 2 x ˙ y ˙ .
These LFIs establish the superintegrability of the geodesic Equations (87) and (88).

Integration of the Constrained Geodesics (87)–(89)

To simplify the notation, we introduce the constants a 0 = E 0 k 0 and a D = I D k with D = 1 , 2 , 3 .
Then, the fixed energy constraint (89) and the LFIs (94) are written:
a 0 ( x + y ) 2 = x ˙ y ˙
a 1 ( x + y ) 2 = y 2 x ˙ x 2 y ˙
a 2 ( x + y ) 2 = y x ˙ + x y ˙
a 3 ( x + y ) 2 = x ˙ y ˙ .
Solving Equations (97) and (98) wrt x ˙ and y ˙ , we find that:
x ˙ = ( a 2 + a 3 x ) ( x + y )
y ˙ = ( a 2 a 3 y ) ( x + y ) .
Replacing (99) and (100) in the remaining Equations (95) and (96), we obtain, respectively, the following:
a 0 = ( a 2 a 3 y ) ( a 2 + a 3 x )
0 = ( a 3 x + a 2 ) y 2 + ( a 3 x 2 a 1 ) y x ( a 2 x + a 1 ) .
Solving the quadratic algebraic Equation (102) wrt y, we find the solution
y = a 2 x + a 1 a 3 x + a 2 .
We note that the solution y = x is not acceptable, because of the requirement x + y 0 ; otherwise, the function f = k ( x + y ) 2 is not well-defined.
Replacing (103) into (101), we obtain the relation
a 0 = a 2 2 a 1 a 3 0
where a 0 is a non-zero constant fixed by the energy E 0 of the geodesics. This means that the constants a 1 , a 2 and a 3 take only those values which satisfy the relation (104) as fixed by the ‘energy’ constant a 0 . We note also that the orbit (103) coincides with Equation (6.17) of [22], if we rename the constants as follows: a 1 = κ 2 2 2 κ 1 2 R , a 2 = 2 κ 1 κ 2 and a 3 = 2 κ 1 2 R .
Replacing the orbit (103) in (99), we find:
x ˙ = a 3 x 2 + 2 a 2 x + a 1 t = d x a 3 x 2 + 2 a 2 x + a 1 .
There are two cases to be considered.
(a)
Case a 3 = 0 .
In this case, from Equations (103) and (105), we find the parametric geodesic solutions:
x ( t ) = c 1 e 2 a 2 t + c 2 , y ( t ) = c 1 e 2 a 2 t c 2
where c 1 is an arbitrary constant and the constant c 2 a 1 2 a 2 . Equation (104) implies that the constant a 2 is fixed by the energy constraint a 0 = a 2 2 0 , while the remaining Equation (100) is satisfied identically.
(b)
Case a 3 0 .
In this case, from Equations (103) and (105), we find the parametric geodesic solutions:
x ( t ) = a 0 a 3 tan a 0 t + c 0 a 2 a 3 , y ( t ) = a 0 a 3 1 tan a 0 t + c 0 + a 2 a 3
where c 0 is an arbitrary constant and a 0 = a 2 2 a 1 a 3 < 0 . From Equation (104), the constant a 1 = a 2 2 a 3 a 0 a 3 , while the remaining Equation (100) is satisfied identically.
Remark 1.
The choice f ( x , y ) = x made in [22] (see Section 6, Equation (6.1)) leads to the flat Lorentzian metric and the solution of the geodesic equations is straightforward. Indeed, the coordinate transformation q a = ( x , y ) q ¯ a = ( u , v ) given by the relations:
u = y x 2 4 , v = y + x 2 4
brings the metric γ a b = 0 x x 0 to its canonical form η a b = d i a g ( 1 , + 1 ) .
In the coordinates ( u , v ) , the solution of the constrained geodesic equations is:
u ( t ) = k 1 t + k 2 , v ( t ) = k 3 t + k 4
where k 1 , k 2 , k 3 and k 4 are arbitrary constants, while the associated quadratic constraint (70) implies that
η a b q ¯ ˙ a q ¯ ˙ b = 2 E 0 0 u ˙ 2 + v ˙ 2 = 2 E 0 E 0 = 1 2 k 3 2 k 1 2 .
Therefore, in this case, the application of neither the exact nor the conditional FIs is necessary.

6.2. The Case of a Metric That Does Not Possess KVs: QFIs

In this case, we shall use the QFI I ( ) 1 given in (71) for = 0 (To simplify the notation, we set C ( 0 ) a b = C a b ( q ) and X ( 0 ) a = X a ( q ) ), that is,
I ( 0 ) 1 = C a b ( q ) q ˙ a q ˙ b + G ( q )
where C a b is a CKT of the metric (78) with associated vector X a ( q ) and the function G ( q ) is such that
G , a = 2 E 0 X a .
Replacing C a b from (82) in (108), the QFI is written
I ( 0 ) 1 = f 2 A 1 ( y ) x ˙ 2 + A 2 ( x ) y ˙ 2 + G ( q ) .
Replacing X a from (83) in the condition (109), we obtain the following system of PDEs:
G , x = 2 E 0 f , y A 1 ( y ) + f 2 d A 1 d y , G , y = 2 E 0 f , x A 2 ( x ) + f 2 d A 2 d x .
Taking the integrability condition G , [ x y ] = 0 of the function G ( x , y ) , we find the second order PDE:
f , y y A 1 ( y ) f , x x A 2 ( x ) + 3 2 f , y d A 1 d y f , x d A 2 d x + f 2 d 2 A 1 d y 2 d 2 A 2 d x 2 = 0 .
Equation (112) cannot be solved in full generality because it has three unknown functions f ( x , y ) , A 1 ( y ) and A 2 ( x ) ; therefore, it is an overdetermined PDE. In order to find a solution, we must fix two of the three unknowns.
First, we choose f ( x , y ) = x 3 e y ( x + e y ) .
For this choice, the metric (78) is (see Equation (6.26) of [22])
γ a b = x 3 e y ( x + e y ) 0 1 1 0
with Ricci scalar R = 2 x 3 ( x + e y ) 3 and Riemann tensor R 1212 = x 3 e 2 y x + e y . The metric (113) does not possess KVs (or non-trivial second order KTs) and, therefore, to assess the integrability of its corresponding constrained geodesics, we should focus to the proper CKVs and CKTs.
The constrained system (75)–(77) becomes:
x ¨ + 1 x + e y + 3 x x ˙ 2 = 0
y ¨ + e y x + e y + 1 y ˙ 2 = 0
x 3 e y ( x + e y ) x ˙ y ˙ = E 0 .
In order to solve (112), we have to fix either A 1 ( y ) or A 2 ( x ) . We choose A 2 ( x ) = 0 and the PDE (112) gives the first degree polynomial in x
A 1 + 3 2 d A 1 d y + 1 2 d 2 A 1 d y 2 x + e y 4 A 1 + 3 d A 1 d y + 1 2 d 2 A 1 d y 2 = 0
whose solution is A 1 ( y ) = e 2 y . We note that if instead of A 2 ( x ) = 0 we had chosen A 1 ( y ) = 0 , the PDE (112) would have given the trivial solution A 2 ( x ) = 0
Having fixed the three unknown functions:
f ( x , y ) = x 3 e y ( x + e y ) , A 1 ( y ) = e 2 y , A 2 ( x ) = 0
the system of PDEs (111) becomes:
G , x = 2 E 0 x 3 , G , y = 0
from which we find the function
G = E 0 2 x 4 .
Replacing (118) and (120) in (110), we find for the constrained non-null geodesic Equations (114)–(116) the QFI
I 1 = x 6 ( x + e y ) 2 x ˙ 2 + E 0 2 x 4 .
Using the QFI (121) and the quadratic constraint (116), the constrained non-null geodesics (114)–(116) can be integrated as follows.
First, we assume that I 1 = 0 . Then, the QFI (121) gives
E 0 = 2 x 2 ( x + e y ) 2 x ˙ 2 .
Replacing (122) in the constraint (116), we find the trajectory
y ˙ x ˙ = 2 ( x + e y ) x e y d y d x = 2 ( x + e y ) x e y y = ln c 1 x 2 2 x .
where c 1 is a constant. This result coincides with Equation (6.27c) of [22].
Moreover, replacing (123) in (122), we find
c 1 x 3 x 2 2 x ˙ 2 = E 0 2 > 0
since E 0 0 . Therefore, E 0 < 0 and Equation (124) implies that
c 1 x 3 x 2 d x = ± E 0 2 d t t = ± 2 E 0 c 1 4 x 4 x 3 3 + t 0
where t 0 is an integration constant.
Remark 2.
As we have seen, t is the affine parameter of the geodesics. If we assume E 0 = 1 2 and a reparameterization t = t ( τ ) such that in the new parameter x ( τ ) = τ , Equation (125) determines the parameter transformation
d t d τ = ± 2 τ 2 c 1 τ 1 N ( τ ) .
The function N ( τ ) is the lapse function of the corresponding mini-superspace Lagrangian (see Equation (6.27a) of [22]).

6.3. The Case of a Class of Integrable Lorentzian Toda Systems: QFIs

Another case where QFIs are required is a class of 2d integrable Lorentzian Toda systems [39] which is equivalent to the constrained dynamical system (75)–(77) for (see Equation (6.35) of [22])
f ( x , y ) = k 1 e 2 b 1 x + ( b 2 b 1 ) y + k 2 e 1 2 b 1 x + b 3 y
where k 1 , k 2 , b 1 , b 2 , b 3 are arbitrary non-zero constants and b 1 b 2 .
In this case, the metric (78) is
γ a b = k 1 e 2 b 1 x + ( b 2 b 1 ) y + k 2 e 1 2 b 1 x + b 3 y 0 1 1 0 .
We note that this metric is in fact the Jacobi metric of the system, and it is flat only when b 3 = 2 ( b 2 b 1 ) .
In order to solve (112), we choose A 1 ( y ) = 0 . Then, the PDE (112) gives A 2 ( x ) = e 2 b 1 x .
Having fixed the three unknown functions:
f ( x , y ) = k 1 e 2 b 1 x + ( b 2 b 1 ) y + k 2 e 1 2 b 1 x + b 3 y , A 1 ( y ) = 0 , A 2 ( x ) = e 2 b 1 x
the system of PDEs (111) becomes:
G , x = 0 , G , y = 2 k 1 b 1 E 0 e 2 ( b 2 b 1 ) y
from which we find the function
G = k 1 b 1 E 0 b 1 b 2 e 2 ( b 2 b 1 ) y .
Replacing (128) and (130) in (110), we find for the constrained non-null geodesic Equations (75)–(77) with f ( x , y ) given by (126) the QFI
I 1 = k 1 e 2 b 1 x + ( b 2 b 1 ) y + k 2 e 1 2 b 1 x + b 3 y 2 e 2 b 1 x y ˙ 2 + k 1 b 1 E 0 b 1 b 2 e 2 ( b 2 b 1 ) y .
Using the QFI (131) and the corresponding quadratic constraint, the integration of the constrained non-null geodesics for I 1 = 0 is straightforward (see previous examples). We conclude that the conditional LFI used in [22] is again not necessary because we can solve the problem by using instead the exact QFI (131).

7. Conclusions

The assessment of the integrability of a dynamical system requires the knowledge of ‘enough’ in number (functionally) independent FIs in involution. In this work, we have considered constrained autonomous conservative dynamical systems where the constraint is a specified value of the total energy. This constraint defines in the configuration space a surface on which the trajectories of the dynamical system evolve. We proved a theorem which allows the systematic determination of autonomous and time-dependent exact QFIs of these constrained dynamical systems in terms of the symmetries of the kinetic metric; the latter being defined by the dynamical equations. Specifically, it is found that these QFIs are generated by CKVs and CKTs of the kinetic metric. Furthermore, there are three types of QFIs and specific formulae were given for their determination. It is to be noted that the time-dependent FIs appear for the first time in the literature of constrained systems.
In order to test the consequences and the validity of the theorem, we have considered various applications from previous works in the topic. As a first application, we required the potential to admit a QFI of the second type, and we found a superintegrable potential of the Ermakov type and an integrable potential which for zero energy and zero QFI gives orbits which are circles. Next, in order to discuss the case of geodesics, we set the potential function equal to zero and determined the corresponding QFIs for the cases of null and timelike/spacelike constrained geodesics. The well-known QFIs of the unconstrained geodesics (see e.g., [8,28,29]) are recovered as special cases. As an application in this area, we considered a general example which provides the various cases discussed in [22] using the Dirac–Bergmann formalism. We recovered all results of [22] and showed that one can integrate the constrained systems using exact FIs instead of the more complicated non-local FIs.
Theorem 1 opens new directions in the study of the integrability of constrained dynamical systems. It would be interesting to generalize Theorem 1 for other types of constraints, less trivial than the fixed energy constraint; and to reexamine well-known problems in celestial mechanics using the constraint of specified energy with the possibility to obtain new QFIs.

Author Contributions

Conceptualization, A.M. and M.T.; methodology, A.M. and M.T.; formal analysis, A.M.; writing—original draft preparation, A.M.; writing—review and editing, A.M. and M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Since the considered constrained dynamical system is autonomous, we should use the polynomial method described in [38] in order to solve the system of PDEs. According to this method, one assumes a general polynomial expression in the variable t for both the quantities K a b ( t , q ) and K a ( t , q ) , and replaces these expressions in the system of PDEs.
In particular, we assume that the CKT
K a b ( t , q ) = C ( 0 ) a b ( q ) + P = 1 m 1 C ( P ) a b ( q ) t P P
where C ( P ) a b ( q ) , P = 0 , 1 , , m 1 , is a sequence of second rank symmetric tensors, and the vector
K a ( t , q ) = M = 0 m 2 L ( M ) a ( q ) t M
where L ( M ) a ( q ) , M = 0 , 1 , , m 2 , are arbitrary vectors.
We note that both powers m 1 and m 2 in the above polynomial expressions may be infinite.
Substituting (A1) and (A2) in Equations (13) and (14), we obtain:
X a = X ( 0 ) a + P = 1 m 1 X ( P ) a t P P
ψ = 1 n M = 0 m 2 L ( M ) a ; a t M + P = 1 m 1 C ( P ) a a t P 1
where
X ( P ) a = 1 n + 2 C ( P ) b b ; a + 2 C ( P ) b a ; b , P = 0 , 1 , , m 1 .
Then, the PDEs (9)–(12) and (15)–(16) become:
0 = C ( 0 ) ( a b ; c ) X ( 0 ) ( a γ b c ) + P = 1 m 1 C ( P ) ( a b ; c ) X ( P ) ( a γ b c ) t P P
0 = P = 1 m 1 C ( P ) a b 1 n C ( P ) c c γ a b t P 1 + M = 0 m 2 L ( M ) ( a ; b ) 1 n L ( M ) c ; c γ a b t M
0 = K , a 2 C ( 0 ) a b V , b 2 ( V E 0 ) X ( 0 ) a + M = 1 m 2 M L ( M ) a t M 1 2 P = 1 m 1 C ( P ) a b V , b + ( V E 0 ) X ( P ) a t P P
0 = K , t 2 n ( V E 0 ) P = 1 m 1 C ( P ) a a t P 1 M = 0 m 2 L ( M ) a V , a + 2 n ( V E 0 ) L ( M ) a ; a t M
0 = M = 2 m 2 M ( M 1 ) L ( M ) a t M 2 + 2 P = 1 m 1 1 n ( V E 0 ) C ( P ) b b , a ( V E 0 ) X ( P ) a C ( P ) a b V , b t P 1 + M = 0 m 2 L ( M ) b V , b , a + 2 n ( V E 0 ) L ( M ) b ; b a t M
0 = 2 C ( 0 ) [ a c V , c ; b ] + 2 ( V E 0 ) X ( 0 ) [ a ; b ] M = 1 m 2 M L ( M ) a ; b t M 1 + 2 P = 1 m 1 C ( P ) [ a c V , c ; b ] + ( V E 0 ) X ( P ) [ a ; b ] t P P .
Integrating Equation (A9), we find the scalar
K = 2 n ( V E 0 ) P = 1 m 1 C ( P ) a a t P P + M = 0 m 2 L ( M ) a V , a + 2 n ( V E 0 ) L ( M ) a ; a t M + 1 M + 1 + G ( q )
where G ( q ) is an arbitrary function.
Replacing (A12) in Equation (A8), we get the condition
0 = G , a 2 C ( 0 ) a b V , b 2 ( V E 0 ) X ( 0 ) a + M = 1 m 2 M L ( M ) a t M 1 + M = 0 m 2 L ( M ) b V , b , a + 2 n ( V E 0 ) L ( M ) b ; b a + 2 n V , a L ( M ) b ; b t M + 1 M + 1 + 2 P = 1 m 1 1 n ( V E 0 ) C ( P ) b b , a + 1 n V , a C ( P ) b b C ( P ) a b V , b ( V E 0 ) X ( P ) a t P P .
We note that the integrability conditions (A10) and (A11) of the scalar K are trivially satisfied, because the assumptions (A1) and (A2) bring the PDE (12) into the integrable form (A9) whose integration gives directly the scalar K ( t , q ) (see Equation (A12)) in terms of the arbitrary function G ( q ) . Hence, the integrability conditions of K are substituted by the integrability conditions G , [ a b ] = 0 of G, which in the case of the Euclidean plane E 2 produce the well-known second order Bertrand–Darboux equation [36,40]. However, it can been checked that both Equations (A10) and (A11) are always satisfied identically from the solutions of the other equations of the system.
Moreover, Equation (A6) implies that the second rank symmetric tensors C ( P ) a b ( q ) with P = 0 , 1 , , m 1 are second order CKTs of the kinetic metric γ a b with associated vectors the quantities X ( P ) a given by (A5). We conclude that the system of equations we have to solve consists of Equations (A7) and (A13).
To solve Equations (A7) and (A13), we consider the following cases
I. Case m 1 = m 2 m (both m 1 and m 2 are finite)
In this case, Equations (A7) and (A13) become:
0 = k = 0 m 1 C ( k + 1 ) a b 1 n C ( k + 1 ) c c γ a b + L ( k ) ( a ; b ) 1 n L ( k ) c ; c γ a b t k + L ( m ) ( a ; b ) 1 n L ( m ) c ; c γ a b t m
0 = G , a 2 C ( 0 ) a b V , b 2 ( V E 0 ) X ( 0 ) a + L ( 1 ) a + k = 1 m 1 L ( k 1 ) b V , b , a 2 C ( k ) a b V , b + k ( k + 1 ) L ( k + 1 ) a + 2 n ( V E 0 ) C ( k ) b b , a + L ( k 1 ) b ; b a n X ( k ) a + 2 n V , a C ( k ) b b + L ( k 1 ) b ; b t k k + L ( m 1 ) b V , b , a 2 C ( m ) a b V , b + 2 n ( V E 0 ) C ( m ) b b , a + L ( m 1 ) b ; b a n X ( m ) a + 2 n V , a C ( m ) b b + L ( m 1 ) b ; b t m m + L ( m ) b V , b , a + 2 n ( V E 0 ) L ( m ) b ; b a + 2 n V , a L ( m ) b ; b t m + 1 m + 1 .
Equation (A14) implies that L ( m ) a is a CKV of γ a b with conformal factor ψ ( m ) = 1 n L ( m ) a ; a and the CKTs
C ( k ) a b = f ( k ) γ a b L ( k 1 ) ( a ; b ) , k = 1 , 2 , , m
where f ( k ) 1 n C ( k ) c c + L ( k 1 ) c ; c .
From (A16), we have:
L ( k 1 ) ( a ; b ) = f ( k ) γ a b C ( k ) a b L ( k 1 ) ( ( a ; b ) ; c ) = f ( k ) , ( a γ b c ) C ( k ) ( a b ; c ) = X ( k ) ( a γ b c ) L ( k 1 ) ( ( a ; b ) ; c ) = Y ( k 1 ) ( a γ b c )
where the vector Y ( k 1 ) a f ( k ) , a X ( k ) a . Therefore, the vectors L ( k 1 ) a , k = 1 , 2 , , m , produce the second order reducible CKTs L ( k 1 ) ( a ; b ) of γ a b with associated vectors
Y ( k 1 ) a = f ( k ) , a X ( k ) a = f ( k ) , a 1 n + 2 C ( k ) b b ; a + 2 C ( k ) b a ; b = f ( k ) , a 1 n + 2 n f ( k ) , a L ( k 1 ) b ; b a + 2 f ( k ) δ a b γ b c L ( k 1 ) ( c ; a ) ; b = 1 n + 2 L ( k 1 ) b ; b a + 2 γ b c L ( k 1 ) ( c ; a ) ; b
which, as expected, are in accordance with the defining relation (29).
From Equation (A15), we find the following conditions (We have replaced the CKTs C ( k ) a b with k = 1 , 2 , , m from the relations (A16)):
L ( m ) b V , b , a = 2 ( V E 0 ) ψ ( m ) , a 2 ψ ( m ) V , a L ( m ) a V , a + 2 ( V E 0 ) ψ ( m ) = s = c o n s t
L ( m 1 ) b V , b , a = 2 L ( m 1 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( m 1 ) a
L ( k 1 ) b V , b , a = 2 L ( k 1 ) ( a ; b ) V , b k ( k + 1 ) L ( k + 1 ) a 2 ( V E 0 ) Y ( k 1 ) a , k = 1 , 2 , , m 1
G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a L ( 1 ) a ( m > 0 ) .
The notation L ( 1 ) a ( m > 0 ) indicates that the vector L ( 1 ) a exists only when m > 0 . If m = 0 , then L ( 1 ) a = 0 .
From (A12), we find that the scalar
K = k = 0 m 1 2 n ( V E 0 ) C ( k + 1 ) a a + L ( k ) a V , a + 2 n ( V E 0 ) L ( k ) a ; a t k + 1 k + 1 + L ( m ) a V , a + 2 n ( V E 0 ) L ( m ) a ; a t m + 1 m + 1 + G ( q ) .
Then, the QFI (6) is
I ( m ) = k = 1 m f ( k ) γ a b t k k k = 1 m L ( k 1 ) ( a ; b ) t k k + C ( 0 ) a b q ˙ a q ˙ b + k = 0 m t k L ( k ) a q ˙ a + L ( m ) a V , a + 2 ( V E 0 ) ψ ( m ) t m + 1 m + 1 + k = 1 m L ( k 1 ) a V , a + 2 ( V E 0 ) f ( k ) t k k + G ( q ) = k = 1 m L ( k 1 ) ( a ; b ) t k k + C ( 0 ) a b q ˙ a q ˙ b + k = 0 m t k L ( k ) a q ˙ a + L ( m ) a V , a + 2 ( V E 0 ) ψ ( m ) t m + 1 m + 1 + k = 1 m L ( k 1 ) a V , a t k k + G ( q ) + k = 1 m f ( k ) γ a b q ˙ a q ˙ b + 2 ( V E 0 ) t k k = 0
I ( m ) = k = 1 m L ( k 1 ) ( a ; b ) t k k + C ( 0 ) a b q ˙ a q ˙ b + k = 0 m t k L ( k ) a q ˙ a + s t m + 1 m + 1 + k = 1 m L ( k 1 ) a V , a t k k + G ( q ) .
In (A22), the index ( m ) indicates the degree of the time-dependence of the coefficients, the quantities C ( 0 ) a b and L ( k ) ( a ; b ) for k = 0 , 1 , , m 1 are second order CKTs with associated vectors X ( 0 ) a and Y ( k ) a , respectively, L ( m ) a is a CKV with conformal factor ψ ( m ) , while the constant s, the vectors L ( k ) a and the function G ( q ) satisfy the conditions (A17)–(A20).
In what follows, we show that the QFI I ( m ) given in (A22) consists of two independent QFIs.
For small values of m, we have:
- For m = 0 .
The QFI is
I ( 0 ) = C ( 0 ) a b q ˙ a q ˙ b + L ( 0 ) a q ˙ a + s t + G ( q )
where C ( 0 ) a b is a CKT with associated vector X ( 0 ) a , L ( 0 ) a is a CKV with conformal factor ψ ( 0 ) such that L ( 0 ) a V , a + 2 ( V E 0 ) ψ ( 0 ) = s and the function G ( q ) satisfies the condition G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a .
The QFI (A23) consists of the independent FIs:
I ( 0 , 1 ) = C ( 0 ) a b q ˙ a q ˙ b + G ( q ) I ( 0 , 2 ) = L ( 0 ) a q ˙ a + s t .
- For m = 1 .
The QFI is
I ( 1 ) = t L ( 0 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + t L ( 1 ) a q ˙ a + L ( 0 ) a q ˙ a + s t 2 2 + t L ( 0 ) a V , a + G ( q )
where C ( 0 ) a b and L ( 0 ) ( a ; b ) are CKTs with associated vectors X ( 0 ) a and Y ( 0 ) a , respectively, L ( 1 ) a is a CKV with conformal factor ψ ( 1 ) such that L ( 1 ) a V , a + 2 ( V E 0 ) ψ ( 1 ) = s , while the vector L ( 0 ) a and the function G ( q ) satisfy the conditions:
L ( 0 ) b V , b , a = 2 L ( 0 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 0 ) a G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a L ( 1 ) a .
The QFI (A24) consists of the independent FIs:
I ( 1 , 1 ) = C ( 0 ) a b q ˙ a q ˙ b + t L ( 1 ) a q ˙ a + s t 2 2 + G ( q ) I ( 1 , 2 ) = t L ( 0 ) ( a ; b ) q ˙ a q ˙ b + L ( 0 ) a q ˙ a + t L ( 0 ) a V , a .
- For m = 2 .
The QFI is
I ( 2 ) = t 2 2 L ( 1 ) ( a ; b ) t L ( 0 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + t 2 L ( 2 ) a q ˙ a + t L ( 1 ) a q ˙ a + L ( 0 ) a q ˙ a + s t 3 3 + t 2 2 L ( 1 ) a V , a + t L ( 0 ) a V , a + G ( q )
where C ( 0 ) a b and L ( k ) ( a ; b ) for k = 0 , 1 are CKTs with associated vectors X ( 0 ) a and Y ( k ) a , respectively, L ( 2 ) a is a CKV with conformal factor ψ ( 2 ) such that L ( 2 ) a V , a + 2 ( V E 0 ) ψ ( 2 ) = s , while the vectors L ( k ) a and the function G ( q ) satisfy the conditions:
L ( 1 ) b V , b , a = 2 L ( 1 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 1 ) a L ( 0 ) b V , b , a = 2 L ( 0 ) ( a ; b ) V , b 2 L ( 2 ) a 2 ( V E 0 ) Y ( 0 ) a G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a L ( 1 ) a .
The QFI (A25) consists of the independent FIs:
I ( 2 , 1 ) = t 2 2 L ( 1 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + t L ( 1 ) a q ˙ a + t 2 2 L ( 1 ) a V , a + G ( q ) I ( 2 , 2 ) = t L ( 0 ) ( a ; b ) q ˙ a q ˙ b + t 2 L ( 2 ) a q ˙ a + L ( 0 ) a q ˙ a + s t 3 3 + t L ( 0 ) a V , a .
- For m = 3 .
I ( 3 ) = t 3 3 L ( 2 ) ( a ; b ) t 2 2 L ( 1 ) ( a ; b ) t L ( 0 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + t 3 L ( 3 ) a q ˙ a + t 2 L ( 2 ) a q ˙ a + t L ( 1 ) a q ˙ a + L ( 0 ) a q ˙ a + s t 4 4 + t 3 3 L ( 2 ) a V , a + t 2 2 L ( 1 ) a V , a + t L ( 0 ) a V , a + G ( q )
where C ( 0 ) a b and L ( k ) ( a ; b ) for k = 0 , 1 , 2 are CKTs with associated vectors X ( 0 ) a and Y ( k ) a , respectively, L ( 3 ) a is a CKV with conformal factor ψ ( 3 ) such that L ( 3 ) a V , a + 2 ( V E 0 ) ψ ( 3 ) = s , while the vectors L ( k ) a and the function G ( q ) satisfy the conditions:
L ( 2 ) b V , b , a = 2 L ( 2 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 2 ) a L ( 1 ) b V , b , a = 2 L ( 1 ) ( a ; b ) V , b 6 L ( 3 ) a 2 ( V E 0 ) Y ( 1 ) a L ( 0 ) b V , b , a = 2 L ( 0 ) ( a ; b ) V , b 2 L ( 2 ) a 2 ( V E 0 ) Y ( 0 ) a G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a L ( 1 ) a .
The QFI (A26) consists of the independent FIs:
I ( 3 , 1 ) = t 2 2 L ( 1 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + t 3 L ( 3 ) a q ˙ a + t L ( 1 ) a q ˙ a + s t 4 4 + t 2 2 L ( 1 ) a V , a + G ( q ) I ( 3 , 2 ) = t 3 3 L ( 2 ) ( a ; b ) t L ( 0 ) ( a ; b ) q ˙ a q ˙ b + t 2 L ( 2 ) a q ˙ a + L ( 0 ) a q ˙ a + t 3 3 L ( 2 ) a V , a + t L ( 0 ) a V , a .
By mathematical induction, it is proved that the QFI I ( m ) consists of the following two independent QFIs:
a.
I ( ) 1 = k = 1 t 2 k 2 k L ( 2 k 1 ) ( a ; b ) + C ( 0 ) a b q ˙ a q ˙ b + k = 1 t 2 k 1 L ( 2 k 1 ) a q ˙ a + k = 1 t 2 k 2 k L ( 2 k 1 ) a V , a + G ( q )
where C ( 0 ) a b and L ( 2 k 1 ) ( a ; b ) for k = 1 , 2 , , are CKTs with associated vectors X ( 0 ) a and Y ( 2 k 1 ) a , respectively, and the vectors L ( 2 k 1 ) a and the function G ( q ) satisfy the conditions:
L ( 2 1 ) b V , b , a = 2 L ( 2 1 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 2 1 ) a L ( 2 k 1 ) b V , b , a = 2 L ( 2 k 1 ) ( a ; b ) V , b 2 k ( 2 k + 1 ) L ( 2 k + 1 ) a 2 ( V E 0 ) Y ( 2 k 1 ) a , k = 1 , 2 , , 1 G , a = 2 C ( 0 ) a b V , b + 2 ( V E 0 ) X ( 0 ) a L ( 1 ) a ( > 0 ) .
b.
I ( ) 2 = k = 0 t 2 k + 1 2 k + 1 L ( 2 k ) ( a ; b ) q ˙ a q ˙ b + t 2 k L ( 2 k ) a q ˙ a + t 2 k + 1 2 k + 1 L ( 2 k ) a V , a
where L ( 2 k ) ( a ; b ) for k = 0 , 1 , , are CKTs with associated vectors Y ( 2 k ) a and the vectors L ( 2 k ) a satisfy the conditions:
L ( 2 ) b V , b , a = 2 L ( 2 ) ( a ; b ) V , b 2 ( V E 0 ) Y ( 2 ) a L ( 2 k ) b V , b , a = 2 L ( 2 k ) ( a ; b ) V , b 2 ( k + 1 ) ( 2 k + 1 ) L ( 2 k + 2 ) a 2 ( V E 0 ) Y ( 2 k ) a , k = 0 , 1 , , 1 .
II. Case(MDPI: Please confirm if keep bold and underline.) m 1 m 2 . ( m 1 or m 2 may be infinite)
We find QFIs that are subcases of those found in Case I and Case III below.
III. Both m 1 and m 2 are infinite.
In this case, the only non-trivial solution is for K a b = e λ t C a b ( q ) and K a = e μ t L a ( q ) , where C a b is a second order CKT with associated vector B a ( q ) and the constants λ μ 0 .
Then, the scalar (A12) is
K = e μ t μ L a V , a + 2 n ( V E 0 ) L a ; a + 2 n ( V E 0 ) C a a e λ t + G ( q )
and the PDEs (A7) and (A13) become:
0 = λ e λ t C a b 1 n C c c γ a b + e μ t L ( a ; b ) 1 n L c ; c γ a b
0 = G , a + e μ t μ L b V , b , a + 2 n ( V E 0 ) L b ; b a + 2 n V , a L b ; b + μ 2 L a + 2 e λ t 1 n ( V E 0 ) C b b , a + 1 n V , a C b b C a b V , b ( V E 0 ) B a .
We consider the following subcases:
(1)
For λ μ .
Equation (A28) implies that C a b = f γ a b is a second order gradient CKT with associated vector B a = f , a , where f ( q ) is a smooth function, and L a is a CKV with conformal factor ρ ( q ) .
From the condition (A29), we get:
G = c o n s t 0 , L a = 1 μ 2 L b V , b + 2 ( V E 0 ) ρ , a .
Recall that by contracting with γ a b the relations C a b = f γ a b and L ( a ; b ) = ρ γ a b , we find f = 1 n C a a and ρ = 1 n L a ; a , respectively.
The scalar (A27) is written
K = e μ t μ L a V , a + 2 ( V E 0 ) ρ + 2 ( V E 0 ) f e λ t + G ( q ) .
The QFI reduces to the LFI
I e ( λ μ ) = e λ t f γ a b q ˙ a q ˙ b + 2 ( V E 0 ) f e λ t = 0 + e μ t L a q ˙ a + e μ t μ L a V , a + 2 ( V E 0 ) ρ = e μ t L a q ˙ a + e μ t μ L a V , a + 2 ( V E 0 ) ρ
where L a = 1 μ 2 L b V , b + 2 ( V E 0 ) ρ , a is a gradient CKV with conformal factor ρ ( q ) .
(2)
For λ = μ .
In this case, Equation (A28) implies that the CKT λ C a b = g γ a b L ( a ; b ) , where g ( q ) is a smooth function. Then, L ( a ; b ) is a reducible CKT with associated vector Y a g , a λ B a .
From the remaining condition (A29), we get:
G = c o n s t 0 , L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a 2 ( V E 0 ) Y a .
Recall that by contracting the relation λ C a b = g γ a b L ( a ; b ) , we find g = 1 n λ C a a + L a ; a .
The scalar (A27) is written
K = e λ t λ L a V , a + 2 ( V E 0 ) g + G ( q ) .
The QFI is (We multiply by the non-zero constant λ )
I e ( λ = μ ) = e λ t L ( a ; b ) q ˙ a q ˙ b + λ L a q ˙ a + L a V , a + g e λ t γ a b q ˙ a q ˙ b + 2 ( V E 0 ) = 0 = e λ t L ( a ; b ) q ˙ a q ˙ b + λ L a q ˙ a + L a V , a
where L ( a ; b ) is a reducible CKT with associated vector Y a such that L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a 2 ( V E 0 ) Y a .
We note that the LFI I e ( λ μ ) is derived from I e ( λ = μ ) in the case that L a is a CKV.
The above completes the proof of Theorem 1.

References

  1. Contopoulos, G. A third integral of motion in a galaxy. Z. Astrophys. 1960, 49, 273. [Google Scholar]
  2. Lynden-Bell, D. Stellar Dynamics: Only isolating integrals should be used in Jeans’ Theorem. Mon. Not. R. Astron. Soc. 1962, 124, 1–9. [Google Scholar] [CrossRef]
  3. Contopoulos, G. On the Existence of a Third Integral of Motion. Astron. J. 1963, 68, 1. [Google Scholar] [CrossRef]
  4. Hénon, M.; Heiles, C. The Applicability of the Third Integral of Motion: Some Numerical Experiments. Astron. J. 1964, 69, 73. [Google Scholar] [CrossRef]
  5. Contopoulos, G. A Review of the “Third” Integral. Math. Eng. 2020, 2, 472. [Google Scholar] [CrossRef]
  6. Eisenhart, L.P. Dynamical trajectories and geodesics. Ann. Math. 1928, 30, 591. [Google Scholar] [CrossRef]
  7. Benn, I.M. Geodesics and Killing tensors in mechanics. J. Math. Phys. 2006, 47, 022903. [Google Scholar] [CrossRef]
  8. Katzin, G.H.; Levine, J. Geodesic first integrals with explicit path-parameter dependence in Riemannian space-times. J. Math. Phys. 1981, 22, 1878. [Google Scholar] [CrossRef]
  9. Mitsopoulos, A.; Tsamparlis, M. Higher order first integrals of autonomous dynamical systems. J. Geom. Phys. 2021, 170, 104383. [Google Scholar] [CrossRef]
  10. Pin, O.C. Curvature and Mechanics. Adv. Math. 1975, 15, 269. [Google Scholar] [CrossRef]
  11. Abraham, R.; Marsden, J.E. Foundations of Mechanics; Addison-Wesley Publishing Company, Inc.: Boston, MA, USA, 1978. [Google Scholar]
  12. Rosquist, K.; Pucacco, G. Invariants at fixed and arbitrary energy. A unified geometric approach. J. Phys. A Math. Gen. 1995, 28, 3235. [Google Scholar] [CrossRef] [Green Version]
  13. Karlovini, M.; Rosquist, K. A unified treatment of cubic invariants at fixed and arbitrary energy. J. Math. Phys. 2000, 41, 370. [Google Scholar] [CrossRef]
  14. Karlovini, M.; Pucacco, G.; Rosquist, K.; Samuelson, L. A unified treatment of quartic invariants at fixed and arbitrary energy. J. Math. Phys. 2002, 43, 4041. [Google Scholar] [CrossRef]
  15. Pucacco, G.; Rosquist, K. Configurational invariants of Hamiltonian systems. J. Math. Phys. 2005, 46, 052902. [Google Scholar] [CrossRef]
  16. Katzin, G.H.; Levine, J. Related First Integral Theorem: A Method for Obtaining Conservation Laws of Dynamical Systems with Geodesic Trajectories in Riemannian Spaces Admitting Symmetries. J. Math. Phys. 1968, 9, 8. [Google Scholar] [CrossRef]
  17. Katzin, G.H. Related integral theorem II. A method for obtaining quadratic constants of the motion for conservative dynamical systems admitting symmetries. J. Math. Phys. 1973, 14, 1213. [Google Scholar] [CrossRef]
  18. Levine, J.; Katzin, G.H. Symmetry mappings of constrained dynamical systems and an associated realted integral theorem. J. Math. Phys. 1973, 14, 1886. [Google Scholar] [CrossRef]
  19. Tsamparlis, M.; Paliathanasis, A. Two-dimensional dynamical systems which admit Lie and Noether symmetries. J. Phys. A Math. Theor. 2011, 44, 175202. [Google Scholar] [CrossRef]
  20. Hodge, W.V.D.; Pedoe, D. Methods of Algebraic Geometry; Cambridge University Press: Cambridge, UK, 1994; Volume I. [Google Scholar]
  21. Dimakis, N.; Terzis, P.A.; Christodoulakis, T. Contact symmetries of constrained quadratic Lagrangians. J. Phys. Conf. Ser. 2016, 670, 012021. [Google Scholar] [CrossRef]
  22. Dimakis, N.; Terzis, P.A.; Christodoulakis, T. Integrability of geodesic motions in curved manifolds through nonlocal conserved charges. Phys. Rev. D 2019, 99, 104061. [Google Scholar] [CrossRef]
  23. Dimakis, N. Hidden symmetries from distortions of the conformal structure. Phys. Rev. D 2022, 106, 024043. [Google Scholar] [CrossRef]
  24. Dirac, P.A.M. Generalized Hamiltonian Dynamics. Can. J. Math. 1950, 2, 129. [Google Scholar] [CrossRef]
  25. Anderson, J.L.; Bergmann, P.G. Constraints in Covariant Field Theories. Phys. Rev. 1951, 83, 1018. [Google Scholar] [CrossRef]
  26. Dirac, P.A.M. Generalized Hamiltonian Dynamics. Proc. R. Soc. Lond. A 1958, 246, 326. [Google Scholar] [CrossRef]
  27. Dirac, P.A.M. Lectrures on Quantum Mechanics; Yeshiva University Press: New York, NY, USA, 1964. [Google Scholar]
  28. Tsamparlis, M.; Mitsopoulos, A. Quadratic first integrals of autonomous conservative dynamical systems. J. Math. Phys. 2020, 61, 072703. [Google Scholar] [CrossRef]
  29. Tsamparlis, M.; Mitsopoulos, A. First integrals of holonomic systems without Noether symmetries. J. Math. Phys. 2020, 61, 122701. [Google Scholar] [CrossRef]
  30. Terzis, P.A.; Dimakis, N.; Christodoulakis, T.; Paliathanasis, A.; Tsamparlis, M. Variational contact symmetries of constrained Lagrangians. J. Geom. Phys. 2016, 101, 52. [Google Scholar] [CrossRef]
  31. Rani, R.; Edgar, S.B.; Barnes, A. Killing tensors and conformal Killing tensors from conformal Killing vectors. Class. Quantum Gravity 2003, 20, 1929. [Google Scholar] [CrossRef]
  32. Prince, G. Homothetic Killing tensors. Phys. Lett. A 1983, 97, 133. [Google Scholar] [CrossRef]
  33. Walker, M.; Penrose, R. On Quadratic First Integrals of the Geodesic Equations for Type {22} Spacetimes. Commun. Math. Phys. 1970, 18, 265. [Google Scholar] [CrossRef]
  34. Weir, G.J. Conformal Killing tensors in reducible spaces. J. Math. Phys. 1977, 18, 1782. [Google Scholar] [CrossRef]
  35. Hietarinta, J. Direct methods for the search of the second invariant. Phys. Rep. 1987, 147, 87. [Google Scholar] [CrossRef]
  36. Mitsopoulos, A.; Tsamparlis, M.; Paliathanasis, A. Integrable and superintegrable potentials of 2d autonomous conservative dynamical systems. Symmetry 2020, 12, 1655. [Google Scholar] [CrossRef]
  37. Ibragimov, N.H.; Kara, A.H.; Mahomed, F.H. Lie-Bäcklund and Noether Symmetries with Applications. Nonlinear Dyn. 1998, 15, 115. [Google Scholar] [CrossRef]
  38. Mitsopoulos, A.; Tsamparlis, M. Quadratic first integrals of time-dependent dynamical systems of the form q ¨ a = Γ a bc q ˙ b q ˙ c ω ( t ) Q a ( q ) . Mathematics 2021, 9, 1503. [Google Scholar] [CrossRef]
  39. Gavrilov, V.R.; Melnikov, V.N. Integration of D-dimensional cosmological models with two factor spaces by reduction to the generalized Emden-Fowler equation. Theor. Math. Phys. 1998, 114, 335. [Google Scholar] [CrossRef] [Green Version]
  40. Darboux, G. Sur un probléme de mècanique. Arch. Neerl. Sci. 1901, 6, 371. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mitsopoulos, A.; Tsamparlis, M. Quadratic First Integrals of Constrained Autonomous Conservative Dynamical Systems with Fixed Energy. Symmetry 2022, 14, 1870. https://doi.org/10.3390/sym14091870

AMA Style

Mitsopoulos A, Tsamparlis M. Quadratic First Integrals of Constrained Autonomous Conservative Dynamical Systems with Fixed Energy. Symmetry. 2022; 14(9):1870. https://doi.org/10.3390/sym14091870

Chicago/Turabian Style

Mitsopoulos, Antonios, and Michael Tsamparlis. 2022. "Quadratic First Integrals of Constrained Autonomous Conservative Dynamical Systems with Fixed Energy" Symmetry 14, no. 9: 1870. https://doi.org/10.3390/sym14091870

APA Style

Mitsopoulos, A., & Tsamparlis, M. (2022). Quadratic First Integrals of Constrained Autonomous Conservative Dynamical Systems with Fixed Energy. Symmetry, 14(9), 1870. https://doi.org/10.3390/sym14091870

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop