Next Article in Journal
User Behavior on Online Social Networks: Relationships among Social Activities and Satisfaction
Previous Article in Journal
Burrows–Wheeler Transform Based Lossless Text Compression Using Keys and Huffman Coding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrable and Superintegrable Potentials of 2d Autonomous Conservative Dynamical Systems

by
Antonios Mitsopoulos
1,*,
Michael Tsamparlis
1 and
Andronikos Paliathanasis
2,3
1
Department of Astronomy-Astrophysics-Mechanics, Faculty of Physics, University of Athens, Panepistemiopolis, 157 83 Athens, Greece
2
Institute of Systems Science, Durban University of Technology, Durban 4000, South Africa
3
Instituto de Ciencias Físicas y Matemáticas, Universidad Austral de Chile, Valdivia 5090000, Chile
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(10), 1655; https://doi.org/10.3390/sym12101655
Submission received: 31 August 2020 / Revised: 27 September 2020 / Accepted: 29 September 2020 / Published: 10 October 2020
(This article belongs to the Special Issue Geometrization of PDEs and Their Solution by Means of Symmetries)

Abstract

:
We consider the generic quadratic first integral (QFI) of the form I = K a b ( t , q ) q ˙ a q ˙ b + K a ( t , q ) q ˙ a + K ( t , q ) and require the condition d I / d t = 0 . The latter results in a system of partial differential equations which involve the tensors K a b ( t , q ) , K a ( t , q ) , K ( t , q ) and the dynamical quantities of the dynamical equations. These equations divide in two sets. The first set involves only geometric quantities of the configuration space and the second set contains the interaction of these quantities with the dynamical fields. A theorem is presented which provides a systematic solution of the system of equations in terms of the collineations of the kinetic metric in the configuration space. This solution being geometric and covariant, applies to higher dimensions and curved spaces. The results are applied to the simple but interesting case of two-dimensional (2d) autonomous conservative Newtonian potentials. It is found that there are two classes of 2d integrable potentials and that superintegrable potentials exist in both classes. We recover most main previous results, which have been obtained by various methods, in a single and systematic way.

1. Introduction

The precise meaning of the solution of a system of differential equations can be cast in several ways [1]. We say that we have determined a closed-form solution for a dynamical system when we have determined a set of explicit functions describing the variation of the dependent variables in terms of the independent variable(s). On the other hand, when we have proved the existence of a sufficient number of independent explicit first integrals and invariants for the dynamical system, we say that we have found an analytic solution of the dynamical equations. In addition, an algebraic solution is found when one has proved the existence of a sufficient number of explicit transformations which permit the reduction of the system of differential equations to a system of algebraic equations. A feature that is central to each of these three equivalent prescriptions of integrability is the existence of explicit functions that are first integrals/invariants or the coefficient functions of the aforementioned transformations.
A first integral (FI) of a dynamical system is a scalar I defined on the phase space of the system such that d I / d t = 0 . The FIs are classified according to the power of the momenta. The linear FIs (LFIs) are linear in the momenta, the quadratic FIs (QFIs) contain products of two momenta, and so on. A dynamical system of n degrees of freedom is called integrable if it admits n (functionally) independent FIs which are in involution [2], that is, their Poisson brackets are zero, i.e., { I i , I j } = 0 . The maximum number of independent FIs that a dynamical system of n degrees of freedom can have is 2 n 1 , and when this is the case, an integrable system is called superintegrable. The above apply to all dynamical systems which are described by dynamical equations independently if they are Lagrangian or Hamiltonian. If the dynamical system is Hamiltonian, then the FIs are defined equivalently by the requirement { H , I } = 0 , where H is the Hamiltonian function of the system.
FIs are important for the determination of the solution and the study of dynamical systems. In particular, when a dynamical system is integrable, then (in principle) the solution of the dynamical equations can be found by means of quadratures. Such dynamical systems are characterized as Liouville integrable. For this reason, the systematic computation of FIs is a topic of active interest for a long time, perhaps by the time of the early Mechanics. Originally, the FIs were concerned in the field of the geometry of surfaces (see, for example, in [3]) where an attempt was made to compute all integrable and superintegrable 2d surfaces. A new dimension to the topic gave the introduction of the theorem of Noether in 1918 [4] which prevails the topic since then. More recently, one more systematic method, but less general than the Noether one, was presented in which one assumes a general form of the QFI and then uses the condition d I / d t = 0 , or the { H , I } = 0 , to find a system of simultaneous equations involving the coefficients defining I (see, for instance, in [5,6,7,8]). The solution of that system of conditions provides us with all the QFIs admitted by a given dynamical system.
The determination of integrable and superintegrable systems is a topic which is in continuous investigation. Obviously a universal method which computes the FIs for all types of dynamical equations independently of their complexity and degrees of freedom is not available. For this reason, the existing studies restrict their considerations to flat spaces or spaces of constant curvature of low dimension (see, e.g., in [3,9,10,11,12,13,14,15] and references therein). The prevailing cases involve the autonomous conservative dynamical systems with two degrees of freedom and the classification of the potential functions in integrable and superintegrable. A comprehensive review of the known integrable and superintegrable 2d autonomous potentials is given in [16].
Besides the two methods mentioned above, other approaches have appeared. For example, Koenigs [17] used coordinate transformations in order to solve the system of equations resulting from the condition { H , I } = 0 . This solution of that system of equations gives the general functional form of the QFIs and the superintegrable free Hamiltonians, that is, the ones which possess two more QFIs—in addition to the Hamiltonian—which are functionally independent. Koenig’s method has been generalized in several works (see in [18] and references cited therein) for two dimensional autonomous conservative systems.
In the present work, we follow the method which uses the solution of the simultaneous system of equations resulting from the condition d I / d t = 0 . This approach has been used extensively, see, e.g., in [5,6,7]; however, always for special cases only. In this work, we use Theorem 1 proved in [19], which gives the general solution of this system in terms of the collineations and the Killing tensors (KTs) of the kinetic metric in the configuration space. This solution is systematic and covariant therefore can be used in higher dimensions and for curved configuration spaces. Furthermore, it is shown that it is directly related to the Noether approach.
Theorem 1 is applied to the case of 2d autonomous conservative dynamical systems in order to determine the integrable and the superintegrable potentials. It is found that the integrable potentials are classified in Class I and Class II and that superintegrable potentials exist in both classes. All potentials together with their QFIs are listed in tables for easy reference. All the results listed in the review paper of [16] as well as in more recent works (see, e.g., in [14,15]) are recovered while some new ones are found which admit time-dependent LFIs and QFIs.

2. Gauged Noether Symmetries and QFIs

We consider an autonomous conservative dynamical system of n degrees of freedom q a with kinetic energy T = 1 2 γ a b q ˙ a q ˙ b where q ˙ a = d q a d t . We define in the configuration space of the system the kinetic metric γ a b by the requirement γ a b = 2 T q ˙ a q ˙ b . When the dynamical system is regular, that is, det 2 T q ˙ a q ˙ b 0 , it can be shown that the dynamical equations can be written in the form
q ¨ a = Γ b c a q ˙ b q ˙ c V ( q ) , a
where Γ b c a are the Riemann connection coefficients defined by the kinetic metric γ a b , V ( q ) stands for the conservative forces, a comma indicates the partial derivative and the Einstein summation convention is used. Finally, the metric γ a b is used for lowering and raising the indices.
The main methods for the determination of the FIs of Equation (1) are
  • The theorem of Noether which is the standard one and requires a Lagrangian.
  • The direct method (see, e.g., in [5,6,7,18,20]), which is not applied widely, uses only the dynamical equations and involves the solution of the system of equations resulting from the condition d I / d t = 0 .
The two methods are related as follows.
In the Noether approach the Noether symmetries are generated by vector fields (We restrict our considerations to vector fields in the jet space J 1 ( t , q , q ˙ ) .).
X = ξ ( t , q , q ˙ ) t + η a ( t , q , q ˙ ) q a
whose first prolongation X [ 1 ] in the jet space J 1 ( t , q , q ˙ ) is given by
X [ 1 ] = ξ t + η a q a + η ˙ a q ˙ a ξ ˙ q ˙ a .
Let L = T V be the Lagrangian of the dynamical system. The Noether symmetries of Equation (1) are its Lie symmetries which satisfy in addition the Noether condition
X 1 L + L ξ ˙ = f ˙
where f ( t , q , q ˙ ) is the Noether or the gauge function. According to the theorem of Noether, a Noether symmetry produces the FI
I = f L ξ L q ˙ a η a ξ q ˙ a .
The Noether symmetries are classified in a formal way in two classes (The original paper of Noether does not distinguish these classes. For a recent enlightening discussion of Noether theorem see in [21] and references therein.):
(a)
The point Noether symmetries whose generators are vector fields on the augmented configuration space { t , q a } and usually lead to LFIs.
(b)
The generalized Noether symmetries whose generators are vector fields in the jet space J 1 ( t , q , q ˙ ) which produce FIs of higher degree.
In the present work, we restrict our considerations to generalized Noether symmetries in the first jet space J 1 ( t , q , q ˙ ) , which produce QFIs.
The 2d autonomous potentials which admit point Noether symmetries have already been classified in [12] and more recently recovered and extended in [13]. Furthermore, in [13] it has been shown that the generators of the point Noether symmetries are the elements of the homothetic algebra of the kinetic metric. Obviously that firm result is not expected to apply in the case of generalized Noether symmetries which form an infinite dimensional Lie group.
It is well known [22] that the generalized Noether symmetries have one extra degree of freedom (being special generalized Lie symmetries) which is removed if we consider the gauge condition ξ = 0 , which we assume to be the case. Therefore, the (gauged) Noether symmetries we consider are generated by vector fields of the form X = η a ( q , q ˙ , q ¨ , . . ) q a and accordingly the Noether condition and the corresponding FI are simplified as follows,
X [ 1 ] L = f ˙ , I = f L q ˙ a η a .
In the direct method, one assumes for the QFI the generic expression
I = K a b ( t , q ) q ˙ a q ˙ b + K a ( t , q ) q ˙ a + K ( t , q )
where K a b ( t , q ) , K a ( t , q ) , and K ( t , q ) are unknown tensor quantities and demands the condition (Equivalently, if the system is Hamiltonian, one requires { H , I } I t = 0 where { . , . } is the Poisson bracket.) d I d t = 0 . This condition leads to a system of simultaneous equations among the coefficients K a b ( t , q ) , K a ( t , q ) , and K ( t , q ) whose solution provides all the LFIs and the QFIs of the system of this functional form. The involvement of the specific dynamical system is in the replacement of the term q ¨ a whenever it appears from the dynamical Equation (1).
The direct approach is related to the Noether symmetries because once one has determined the QFI the generator of the corresponding gauged Noether symmetry and the Noether function follow immediately. Indeed for a gauged Noether symmetry (in the gauge ξ = 0 ) relation Equation (4) becomes
I = f L q ˙ a η a .
Replacing L = T V ( q ) we find
I = f η a γ a b q ˙ b = f η a q ˙ a
and using Equation (6) it follows
η a = K a b q ˙ b K a , f = K
that is, we obtain directly the Noether generator and the Noether function from the QFI I by reading the coefficients K a b ( t , q ) , K a ( t , q ) and K ( t , q ) , respectively. It can be proved that (a) the set { K a b q ˙ b K a ; K } does satisfy the gauged Noether condition X [ 1 ] L = f ˙ and (b) the QFI I defined in Equation (6) is not in general Noether invariant (as it is the case for the point Noether symmetries—see proposition 2.2 in [23]). Finally, the (gauged) point Noether symmetries, which are defined by the vector K a ( K a b = 0 ) , give the LFIs whereas the (gauged) generalized Noether symmetries with K a b 0 give the QFIs.

3. The QFIs of an Autonomous Conservative Dynamical System

It is known (see, e.g., in [5,7]) that condition d I / d t = 0 leads to the following system of equations (using the dynamical Equations (1) to replace q ¨ a whenever it appears the condition d I / d t = 0 is written)
K ( a b ; c ) q ˙ a q ˙ b q ˙ c + K a b , t + K ( a ; b ) q ˙ a q ˙ b + K a , t + K , a 2 K a b V , b q ˙ a + K , t K a V , a = 0 .
(10) K ( a b ; c ) = 0 (11) K a b , t + K ( a ; b ) = 0 (12) 2 K a b V , b + K a , t + K , a = 0 (13) K , t K a V , a = 0 .
Here, round/square brackets indicate symmetrization/antisymmetrization of the enclosed indices and a semicolon denotes the Riemannian covariant derivative. In the special case of a scalar function, for example, the potential V, it holds that V , a = V ; a .
Condition K ( a b ; c ) = 0 implies that K a b is a Killing tensor (KT) of order 2 (possibly zero) of the kinetic metric γ a b . Because γ a b is autonomous we assume
K a b ( t , q ) = g ( t ) C a b ( q )
where g ( t ) is an arbitrary analytic function and C a b ( q ) ( C a b = C b a ) is a KT of order 2 of the metric γ a b . This choice of K a b and Equation (11) indicates that we set
K a ( t , q ) = f ( t ) L a ( q ) + B a ( q )
where f ( t ) is an arbitrary analytic function and L a ( q ) , B a ( q ) are arbitrary vectors. With these choices the system of Equations (10)–(13) becomes
(14) g ( t ) C ( a b ; c ) = 0 (15) g , t C a b + f ( t ) L ( a ; b ) + B ( a ; b ) = 0 (16) 2 g ( t ) C a b V , b + f , t L a + K , a = 0 (17) K , t ( f L a + B a ) V , a = 0 .
Conditions Equations (14)–(17) must be supplemented with the integrability conditions K , a t = K , t a and K , [ a b ] = 0 for the scalar K . The integrability condition K , a t = K , t a gives—if we make use of Equations (16) and (17)—the PDE
f , t t L a + f , t L b A a b + f L b V , b ; a + B b V , b ; a 2 g , t C a b V , b = 0 .
Condition K , [ a b ] = 0 gives the equation known as the second-order Bertrand–Darboux PDE
2 g C [ a c V , c ; b ] f , t L a ; b = 0
where indices enclosed between vertical lines are overlooked by symmetrization or antisymmetrization symbols.
Finally, the system of equations which we have to solve consists of Equations (14)–(19). The general solution of that system in terms of the collineations of the kinetic metric is given in the following Theorem (see [19]).
Theorem 1. 
The functions g ( t ) , f ( t ) are assumed to be analytic so that they may be represented by polynomial expansion as follows,
g ( t ) = k = 0 n c k t k = c 0 + c 1 t + + c n t n
f ( t ) = k = 0 m d k t k = d 0 + d 1 t + + d m t m
where n , m N , or may be infinite, and c k , d k R . Then, the independent QFIs of an autonomous conservative dynamical system are the following.
Integral 1.
I 1 = t 2 2 L ( a ; b ) q ˙ a q ˙ b + C a b q ˙ a q ˙ b + t L a q ˙ a + t 2 2 L a V , a + G ( q )
where C a b , L ( a ; b ) are KTs, L b V , b , a = 2 L ( a ; b ) V , b and G , a = 2 C a b V , b L a .
Integral 2.
I 2 = t 3 3 L ( a ; b ) q ˙ a q ˙ b + t 2 L a q ˙ a + t 3 3 L a V , a t B ( a ; b ) q ˙ a q ˙ b + B a q ˙ a + t B a V , a
where L a , B a are such that L ( a ; b ) , B ( a ; b ) are KTs, L b V , b , a = 2 L ( a ; b ) V , b , and B b V , b , a = 2 B ( a ; b ) V , b 2 L a .
Integral 3.
I 3 = e λ t L ( a ; b ) q ˙ a q ˙ b + λ e λ t L a q ˙ a + e λ t L a V , a
where λ 0 , L a is such that L ( a ; b ) is a KT and L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a .
It can be checked that the FIs listed above produce all the potentials which admit a LFI or a QFI given in [13] and are due to point Noether symmetries. Since, as shown above, these FIs also follow form a gauged velocity dependent Noether symmetry we conclude that there does not exist a one-to-one correspondence between Noether FIs and the type of Noether symmetry. For example the FI of the total energy (Hamiltonian) E = 1 2 γ a b q ˙ a q ˙ b + V ( q ) (case Integral 1 for L a = 0 and C a b = γ a b 2 ) is generated by the point Noether symmetry ξ = 1 , η a = 0 ; f = 0 and also by the gauged generalized Noether symmetry ξ = 0 , η a = 1 2 γ a b q ˙ b ; f = V ( q ) .
The FI I 2 ( L a = 0 ) for B a be a HV with conformal factor ψ = c o n s t is generated by the point Noether symmetry ξ = 2 ψ t , η a = B a ; f = c t such that B a V , a + 2 ψ V + c = 0 , and also by the gauged generalized Noether symmetry ξ = 0 , η a = t ψ γ a b q ˙ b + B a ; f = t B a V , a .
As a final example, we consider the FI I 3 λ for the gradient HV L a = Φ ( q ) , a where Φ ; a b = ψ γ a b with ψ = c o n s t ). This FI is generated by the point Noether symmetry
ξ = 2 ψ λ e λ t , η a = e λ t Φ ( q ) , a ; f = λ e λ t Φ ( q ) c λ e λ t
where λ , c are non-zero constants and Φ , a V , a = 2 ψ V λ 2 Φ + c , and also by the gauged generalized Noether symmetry
ξ = 0 , η a = e λ t λ ψ γ a b q ˙ b + e λ t Φ , a ; f = e λ t λ Φ , a V , a .

4. The Determination of the QFIs

From Theorem 1 follows that for the determination of the QFIs the following problems have to be solved.
  • Determine the KTs of order 2 of the kinetic metric γ a b .
  • Determine the special subspace of KTs of order 2 of the form C a b = L ( a ; b ) where L a is a vector.
  • Determine the KTs satisfying the constraint G , a = 2 C a b V , b .
  • Find all KVs L a of the kinetic metric which satisfy the constraint L a V , a = s where s is a constant, possibly zero.
We note that constraints a. and b. depend only on the kinetic metric. Because the kinetic energy is a positive definite non-singular quadratic 2-form, we can always choose coordinates in which this form reduces either to δ a b or to A ( q ) δ a b . As we know the KTs and all the collineations of a conformally flat metric (of Euclidean or Lorentzian character) [24], we already have the results for all autonomous (Newtonian or special relativistic) conservative dynamical systems.
The involvement of the potential function is only in the constraints c. and d. which also depend on the geometric characteristics of the kinetic metric. There are two different ways to proceed.

4.1. The Potential V q Is Known

In this case the following procedure is used.
(a)
Substitute V in the constraints L a V , a = s and G , a = 2 C a b V , b and find conditions for the defining parameters of L a and C a b .
(b)
From these conditions determine L a , C a b .
(c)
Substitute C a b in the constraint G , a = 2 C a b V , b and find the function G ( q ) .
(d)
Using the above results write the FI I in each case and determine directly the gauged Noether generator and the Noether function.
(e)
Examine if I can be reduced to simpler independent FIs or if it is new.

4.2. The Potential V q Is Unknown

In this case the following algorithm is used.
(a)
Compute the KTs and the KVs of the kinetic metric.
(b)
Solve the PDE L a V , a = s or the ( The integrability conditions for the scalar G are very general PDEs from which one can find only special solutions by making additional simplifying assumptions (e.g., symmetries) involving L a , C a b , and V ( q ) itself. Therefore, one does not find the most general solution. For example, in [25] it is required that the QFI I is axisymmetric, that is, ϕ [ 1 ] I = 0 , where ϕ i [ 1 ] = y x + x y y ˙ x ˙ + x ˙ y ˙ is the first prolongation of the rotation ϕ i = y x + x y . It is proved easily that in this case we have also the constraints L ϕ K a = 0 and L ϕ K a b = 0 .) G , [ a b ] = 0 and find the possible potentials V ( q ) .
(c)
Substitute the potentials and the KTs found in the constraint G , a = 2 C a b V , b and compute the function G ( q ) .
(d)
Write the FI I for each potential and determine the gauged Noether generator and the Noether function.
(e)
Examine if I can be reduced further to simpler independent FIs or if it is a new FI.
In the following sections we assume the potential is not given and apply the second procedure. For that we need first the geometric quantities of the 2d Euclidean plane E 2 .

5. The Geometric Quantities of E 2

Using well-known results (see also in [11,20]), we state the following.
- E 2 admits two gradient Killing vectors (KVs) x , y whose generating functions are x , y respectively and one non-gradient KV (the rotation) y x x y . These vectors can be written collectively
L a = b 1 + b 3 y b 2 b 3 x
where b 1 , b 2 , b 3 are arbitrary constants, possibly zero.
- The general KT of order 2 in E 2 is
C a b = γ y 2 + 2 α y + A γ x y α x β y + C γ x y α x β y + C γ x 2 + 2 β x + B
from which follows
C a b ( q ) q ˙ a q ˙ b = γ y 2 + 2 α y + A x ˙ 2 + 2 γ x y α x β y + C x ˙ y ˙ + γ x 2 + 2 β x + B y ˙ 2
where α , β , γ , A , B , C are arbitrary constants.
- The vectors L a generating KTs of E 2 of the form C a b = L ( a ; b ) are
L a = 2 β y 2 + 2 α x y + A x + a 1 y + a 4 2 α x 2 + 2 β x y + a 3 x + B y + a 2
where a 1 , a 2 , a 3 , a 4 are arbitrary constants.
- The KTs C a b = L ( a ; b ) in E 2 generated from the vector Equation (25) are
C a b = L ( a ; b ) = L x , x 1 2 ( L x , y + L y , x ) 1 2 ( L x , y + L y , x ) L y , y = 2 α y + A α x β y + C α x β y + C 2 β x + B
where (Note that L a in Equation (25) is the sum of the non-proper ACs of E 2 and not of its KVs which give C a b = 0 ). 2 C = a 1 + a 3 . Observe that these KTs are special cases of the general KTs Equation (23) for γ = 0 .
According to Theorem 1 the above are common to all 2d Newtonian systems and what changes in each particular case are the constraints G , a = 2 C a b V , b and L a V , a = s which determine the potential V ( q ) .

6. Computing the Potentials and the FIs

The application of Theorem 1 in the case of E 2 indicates that there are three different ways to find potentials that admit QFIs (other than the Hamiltonian):
(1) The constraint L a V , a = s , which leads to the PDE
( b 1 + b 3 y ) V , x + ( b 2 b 3 x ) V , y s = 0 .
(2) The constraint G , a = 2 C a b V , b , which leads to the second order Bertrand–Darboux PDE ( G , x y = G , y x )
0 = ( γ x y + α x + β y C ) ( V , x x V , y y ) + γ ( y 2 x 2 ) 2 β x + 2 α y + A B V , x y 3 ( γ x + β ) V , y + 3 ( γ y + α ) V , x .
(3) The constraint L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a with (For λ = 0 this constraint is a subcase of G , a = 2 C a b V , b , hence only the case λ 0 must be considered.). λ 0 and the integrability condition L b V , b , x y = L b V , b , y x , which lead to the PDEs
0 = ( 2 β y 2 + 2 α x y + A x + a 1 y + a 4 ) V , x x + ( 2 α x 2 + 2 β x y + a 3 x + B y + a 2 ) V , x y + + ( 6 α x + 2 a 3 + a 1 ) V , y + 3 ( 2 α y + A ) V , x + λ 2 ( 2 β y 2 + 2 α x y + A x + a 1 y + a 4 )
0 = ( 2 α x 2 + 2 β x y + a 3 x + B y + a 2 ) V , y y + ( 2 β y 2 + 2 α x y + A x + a 1 y + a 4 ) V , x y + + 3 ( 2 β x + B ) V , y + ( 6 β y + 2 a 1 + a 3 ) V , x + λ 2 ( 2 α x 2 + 2 β x y + a 3 x + B y + a 2 )
0 = ( α x + β y C ) ( V , x x V , y y ) + 2 β x + 2 α y + A B V , x y 3 β V , y + 3 α V , x + + λ 2 2 ( 6 α x 6 β y + a 1 a 3 ) , 2 C = a 1 + a 3 .
For α = β = 0 and a 1 = a 3 Equation (31) reduces to Equation (28). Therefore, in order to find new potentials one of these conditions must be relaxed. This case of finding potentials is the most difficult because the problem is overdetermined, i.e., we have a system of three PDEs Equations (29)–(31) and only one unknown function, the V ( x , y ) .
In the following sections, we solve these constraints and find the admitted potentials which, as a rule, are integrable. Subsequently, we apply Theorem 1 to each of these potentials in order to compute the admitted FIs and determine which of those are integrable and in particular superintegrable.

7. The Constraint L a V , a = s

The constraint L a V , a = s gives Equation (27) which can be solved using the method of the characteristic equation.
To cover all possible occurrences we have to consider the following cases. (a) b 3 = 0 and b 1 0 (KVs x and x , y ) ; (b) b 3 = b 1 = 0 and b 2 0 (KV y ) ; and (c) b 3 0 (KVs y x x y ; x , y x x y ; and y , y x x y ) . For each case the solution is shown in the following table.
CaseKV V ( x , y )
a b 3 = 0 , b 1 0 s b 1 x + F ( b 1 y b 2 x )
b b 3 = b 1 = 0 , b 2 0 s b 2 y + F ( x )
c b 3 0 s b 3 tan 1 y + b 1 b 3 x + b 2 b 3 + F ( b 1 y + b 3 2 y 2 b 2 x + b 3 2 x 2 )
We shall refer to the above solutions as Class I potentials. In order to determine if these potentials admit QFIs, we apply Theorem 1 to the following potentials resulting from the table above.
V 1 = c x + F ( y b x ) V 2 = c y + F ( x )
V 3 = c tan 1 y + b 1 x + b 2 + F x 2 + y 2 2 + b 1 y b 2 x .
Before we continue we recall that if I 1 , I 2 , , I k are FIs of a given dynamical system then any function f ( I 1 , , I k ) is also a FI of the dynamical system.

7.1. The Potential V 1 = c x + F ( y b x )

Case a. b = 0 and F = λ y .
The potential reduces to V 1 a = c x + λ y .
The irreducible FIs are
L 11 = x ˙ + c t , L 12 = y ˙ + λ t , Q 11 = 1 2 x ˙ 2 + c x , Q 12 = 1 2 y ˙ 2 + λ y .
We note that Q 11 + Q 12 = 1 2 ( x ˙ 2 + y ˙ 2 ) + V = H the Hamiltonian. We compute { Q 11 , Q 12 } = 0 , { L 11 , Q 11 } = c .
The FIs I 1 = Q 11 + Q 12 , I 2 = λ L 11 c L 12 = λ x ˙ c y ˙ and I 3 = Q 11 are functionally independent and satisfy the relations
{ I 1 , I 2 } = { I 1 , I 3 } = 0 , { I 2 , I 3 } = c λ .
Therefore, the potential V 1 a is superintegrable.
We note that the FIs I 2 , Q 11 are, respectively, the FIs (3.1.4) and (3.2.20) of [16].
Case b. d 2 F d w 2 0 and w y b x .
The irreducible FIs are
L 21 = x ˙ + b y ˙ + c t , L 22 = t ( x ˙ + b y ˙ ) ( x + b y ) + c 2 t 2 , Q 21 = ( x ˙ + b y ˙ ) 2 + 2 c ( x + b y ) .
For F ( y b x ) = 1 2 λ 2 y 2 and b = 0 , we have the potential V 1 b = c x 1 2 λ 2 y 2 , λ 0 , which admits the additional time-dependent FI L 23 = e λ t ( y ˙ λ y ) . Observe also that in this case Q 21 reduces to Q e 1 = 1 2 x ˙ 2 + c x which using the Hamiltonian generates the QFI
Q e 2 H Q e 1 = 1 2 y ˙ 2 1 2 λ 2 y 2 .
The LFI L 21 ( c = 0 ) is the (3.1.4) of [16].
We compute { H , L 21 } = L 21 t = c because L 21 is a time-dependent FI.
The potential of the case b is integrable because { H , Q 21 } = 0 .
Moreover,
{ H , L 22 } = L 21 = L 22 t , { L 21 , L 22 } = 1 + b 2 ,
{ Q 21 , L 21 } = 2 c ( 1 + b 2 ) = 2 c { L 21 , L 22 } , { Q 21 , L 22 } = 2 ( 1 + b 2 ) L 21 = 2 { L 21 , L 22 } L 21 .
For the special case V = c x 1 2 λ 2 y 2 we have
{ H , L 23 } = { Q e 2 , L 23 } = λ L 23 = L 23 t , { Q e 1 , L 23 } = 0 .
The triplet Q e 1 , Q e 2 , L 23 proves that this potential is superintegrable.
We note that in [16] only the Class II potentials (to be considered in the next section) are examined for superintegrability (see in [16] p. 108 (3.2.34)–(3.2.36)).
For c 0 , the potential V 1 = c x + F ( y b x ) is not included in [16] because the author seeks for autonomous LFIs of the form (3.1.1) and in that case s = 0 .

7.2. The Potential V 2 = c y + F ( x )

We consider the case F = d 2 F d x 2 0 , because otherwise we retrieve the potential V 1 a discussed above.
The irreducible FIs are
L 31 = y ˙ + c t , Q 31 = 1 2 x ˙ 2 + F ( x ) , Q 32 = 1 2 y ˙ 2 + c y .
Therefore, the potential V 2 is integrable. This potential is also of the form V = F 1 ( x ) + F 2 ( y ) , which is the (3.2.20) of [16].
For F ( x ) = 1 2 λ 2 x 2 , we obtain the potential V 2 a = c y 1 2 λ 2 x 2 , λ 0 , which admits the additional FI L 32 = e λ t ( x ˙ λ x ) . This potential is superintegrable because of the functionally independent triplet Q 31 , Q 32 , and L 32 .

7.3. The Potential V 3 = c tan 1 y + b 1 x + b 2 + F x 2 + y 2 2 + b 1 y b 2 x

We find the time-dependent LFI
L 51 = y x ˙ x y ˙ + b 1 x ˙ + b 2 y ˙ + c t .
For c = 0 this potential is integrable. For c 0 we do not know.
- For c = 0 and F = λ x 2 + y 2 2 + b 1 y b 2 x , λ 0 , the independent FIs are
L 41 = y x ˙ x y ˙ + b 1 x ˙ + b 2 y ˙ , Q 41 = 1 2 x ˙ 2 + 1 2 λ x 2 λ b 2 x , Q 42 = 1 2 y ˙ 2 + 1 2 λ y 2 + λ b 1 y ,
Q 43 = x ˙ y ˙ + λ ( x y + b 1 x b 2 y ) .
Observe that Q 41 + Q 42 = H is the energy of the system. The LFI L 41 is the (3.1.6) of [16]. The functionally independent triplet H, L 41 , Q 41 proves that this potential is superintegrable. We have
{ H , L 41 } = { H , Q 41 } = 0 , { L 41 , Q 41 } = Q 43 + λ b 1 b 2 .
If b 1 = b 2 = 0 and λ = k 2 0 we obtain the superintegrable (A subcase of the above superintegrable potential is the potential V 3 a = λ x 2 + y 2 2 + b 1 y b 2 x .) potential V 3 b = 1 2 k 2 ( x 2 + y 2 ) which admits the additional time-dependent LFIs
L 42 ± = e ± k t ( x ˙ k x ) , L 43 ± = e ± k t ( y ˙ k y ) .
We also compute
{ L 41 , Q 42 } = Q 43 λ b 1 b 2 , { L 41 , Q 43 } = 2 Q 41 2 Q 42 + λ ( b 2 2 b 1 2 )
{ Q 41 , Q 42 } = 0 , { Q 41 , Q 43 } = { Q 43 , Q 42 } = λ L 41 .
In Section 4 of the work in [26], the author has found the superintegrable ClassI potentials V 1 a and V 3 a .
We note that in the review in [16], the time-dependent LFIs of the potentials V 1 a , V 2 are not discussed. In general, in [16] all the time-dependent FIs are ignored, although they can be used to decide the superintegrability of the system.

7.4. Summary

We collect the results for the Class I potentials in the following tables.
PotentialRef. [16]LFIs and QFIs
V 3 ( c 0 ) = c tan 1 y + b 1 x + b 2 + F x 2 + y 2 2 + b 1 y b 2 x - L 51 = y x ˙ x y ˙ + b 1 x ˙ + b 2 y ˙ + c t
Integrable potentials
V 1 = c x + F ( y b x ) , d 2 F d w 2 0 , w y b x - L 21 = x ˙ + b y ˙ + c t ,
L 22 = t ( x ˙ + b y ˙ ) ( x + b y ) + c 2 t 2 ,
Q 21 = ( x ˙ + b y ˙ ) 2 + 2 c ( x + b y )
V 2 = c y + F ( x ) , F 0 (3.2.20) L 31 = y ˙ + c t , Q 31 = 1 2 x ˙ 2 + F ( x ) ,
Q 32 = 1 2 y ˙ 2 + c y
V 3 ( c = 0 ) (3.1.6) L 51 ( c = 0 )
Superintegrable Potentials
PotentialRef. [16]LFIs and QFIs
V 1 a = c x + λ y (3.1.4),
(3.2.20)
L 11 = x ˙ + c t , L 12 = y ˙ + λ t ,
Q 11 = 1 2 x ˙ 2 + c x , Q 12 = 1 2 y ˙ 2 + λ y
V 1 b = c x 1 2 λ 2 y 2 , λ 0 (3.2.20) L 11 , L 22 ( b = 0 ) = t x ˙ x + c 2 t 2 , L 23 = e λ t ( y ˙ λ y ) ,
Q 2 e 1 = 1 2 x ˙ 2 + c x , Q 2 e 2 = 1 2 y ˙ 2 1 2 λ 2 y 2
V 2 a = c y 1 2 λ 2 x 2 , λ 0 (3.2.20) L 31 , Q 31 a = 1 2 x ˙ 2 1 2 λ 2 x 2 , Q 32 , L 32 = e λ t ( x ˙ λ x )
V 3 a = λ x 2 + y 2 2 + b 1 y b 2 x , λ 0 (3.1.6) L 41 = y x ˙ x y ˙ + b 1 x ˙ + b 2 y ˙ , Q 41 = 1 2 x ˙ 2 + 1 2 λ x 2 λ b 2 x ,
Q 42 = 1 2 y ˙ 2 + 1 2 λ y 2 + λ b 1 y , Q 43 = x ˙ y ˙ + λ ( x y + b 1 x b 2 y )
V 3 b = 1 2 k 2 ( x 2 + y 2 ) , k 0 (3.1.5) L 41 b = y x ˙ x y ˙ , Q 41 b = 1 2 x ˙ 2 1 2 k 2 x 2 ,
Q 42 b = 1 2 y ˙ 2 1 2 k 2 y 2 , Q 43 b = x ˙ y ˙ k 2 x y ,
L 42 ± = e ± k t ( x ˙ k x ) , L 43 ± = e ± k t ( y ˙ k y )

8. The Constraint G , a = 2 C ab V , b

In this case, we have the PDE Equation (28)
0 = ( γ x y + α x + β y C ) ( V , x x V , y y ) + γ ( y 2 x 2 ) 2 β x + 2 α y + A B V , x y 3 ( γ x + β ) V , y + 3 ( γ y + α ) V , x .
The potentials which follow from this equation we call Class II potentials. This equation cannot be solved in full generality (see also in [16]), therefore we consider various cases which produce the known FIs. We emphasize that the potentials we find in this section are only a subset of the possible potentials which will follow from the general solution of Equation (32). However the important point here is that we recover the known results with a direct and unified approach which can be used in the future by other authors to discover new integrable and superintegrable potentials in E 2 and in other spaces.
(1) γ 0 , A = B and α = β = C = 0 . Then C a b = γ y 2 + A γ x y γ x y γ x 2 + A and Equation (32) becomes
x y ( V , x x V , y y ) + ( y 2 x 2 ) V , x y 3 x V , y + 3 y V , x = 0
whose solution gives
V 21 = F 1 y x d 1 x 2 + d 2 y 2 + F 2 ( x 2 + y 2 )
where d 1 , d 2 are arbitrary constants.
For the subcase d 1 = d 2 = 1 with A = 0 , we find the QFI
I 11 = ( y x ˙ x y ˙ ) 2 + 2 F 1 y x = ( r 2 θ ˙ ) 2 Φ ( θ )
where r 2 = x 2 + y 2 and θ = tan 1 y x . This is the well-known Ermakov–Lewis invariant; see also (3.2.11) of the work in [16].
For d 1 0 , the potential Equation (34) is written equivalently
V 21 = F 1 y x x 2 + c y 2 + F 2 ( x 2 + y 2 )
where c is an arbitrary constant.
This potential admits QFIs for F 1 = c k y 2 + k x 2 x 2 + ( 2 c ) y 2 . Therefore,
V 21 a = k x 2 + ( 2 c ) y 2 + F 2 ( x 2 + y 2 ) = k x 2 + y 2 + F 2 ( x 2 + y 2 )
with the QFI
I 11 a = ( y x ˙ x y ˙ ) 2 + 2 k ( c 1 ) y 2 x 2 + ( 2 c ) y 2 = ( y x ˙ x y ˙ ) 2 + 2 k ( 1 ) y 2 x 2 + y 2
where 2 c .
For d 1 = 0 , d 2 0 , the potential V 21 becomes
V 21 = F 1 y x y 2 + F 2 ( x 2 + y 2 ) .
This potential admits QFIs for F 1 = k y 2 2 x 2 + y 2 . Then,
V 21 b = k 2 x 2 + y 2 + F ( x 2 + y 2 )
with the QFI
I 11 b = ( y x ˙ x y ˙ ) 2 + k y 2 2 x 2 + y 2 .
Observe that V 21 b is of the form V 21 a ( c = 3 / 2 ) or V 21 a ( = 1 / 2 ) with k ¯ 2 k . Therefore, V 21 b is included in case V 21 a .
(2) γ = 1 and α = β = B = C = 0 . Then C a b = y 2 + A x y x y x 2 and Equation (32) becomes
x y ( V , x x V , y y ) + ( y 2 x 2 + A ) V , x y 3 x V , y + 3 y V , x = 0 .
For A = 0 Equation (38) reduces to Equation (33).
For A 0 the PDE Equation (38) gives the Darboux solution
V 22 = F 1 ( u ) F 2 ( v ) u 2 v 2
where r 2 = x 2 + y 2 , u 2 = r 2 + A + ( r 2 + A ) 2 4 A x 2 1 / 2 and v 2 = r 2 + A ( r 2 + A ) 2 4 A x 2 1 / 2 .
We find the QFI (see (3.2.9) of the work in [16]).
I 21 = ( y x ˙ x y ˙ ) 2 + A x ˙ 2 + v 2 F 1 ( u ) u 2 F 2 ( v ) u 2 v 2 .
(3) γ = 1 , B = A , C = ± i A 0 and α = β = 0 . Then,
C a b = y 2 + A x y ± i A x y ± i A x 2 A
and Equation (32) gives again a potential of the form Equation (39), but with u 2 = r 2 + r 4 4 A ( x ± i y ) 2 1 / 2 and v 2 = r 2 r 4 4 A ( x ± i y ) 2 1 / 2 .
We find the QFI (see (3.2.13) in [16])
I 31 = ( y x ˙ x y ˙ ) 2 + A ( x ˙ ± i y ˙ ) 2 + v 2 F 1 ( u ) u 2 F 2 ( v ) u 2 v 2 .
(4a) α = 1 and β = γ = A = B = C = 0 . Then,
C a b = 2 y x x 0
and Equation (32) becomes
x ( V , x x V , y y ) + 2 y V , x y + 3 V , x = 0
which gives the potential
V 24 = F 1 ( r + y ) + F 2 ( r y ) r
where r 2 = x 2 + y 2 .
We find the QFI (see (3.2.15) in [16])
I 41 = x ˙ ( y x ˙ x y ˙ ) + ( r + y ) F 2 ( r y ) ( r y ) F 1 ( r + y ) r .
(4b) β = 1 and α = γ = A = B = C = 0 . Then,
C a b = 0 y y 2 x
and Equation (32) becomes
y ( V , x x V , y y ) 2 x V , x y 3 V , y = 0
which gives the potential
V 24 b = F 1 ( r + x ) + F 2 ( r x ) r
where r 2 = x 2 + y 2 .
We find the QFI
I 41 b = y ˙ ( x y ˙ y x ˙ ) + ( r + x ) F 2 ( r x ) ( r x ) F 1 ( r + x ) r
Observe that the potential Equation (46) is just the Equation (43) after the rotation x y . All the results of the case 4b can be derived from the case 4a if we apply the transformation x y . For this reason the case 4b is ignored when we search for integrable systems, but in superintegrability the PDE Equation (45) shall be proved useful (see superintegrable potential Equation (74) in Section 8.1).
(5) α = 1 , β = i , A = B = i 4 , C = 1 4 and γ = 0 . Then,
C a b = 2 y + i 4 x + i y + 1 4 x + i y + 1 4 2 i x i 4
and Equation (32) becomes
( x i y 1 4 ) V , x x V , y y + 2 y + i x + i 4 V , x y + 3 i V , y + 3 V , x = 0 .
This is written equivalently
( x i y ) x + i y 2 V 1 4 ( x i y ) 2 V + 3 ( x + i y ) V = 0
and gives the potential
V 25 = w 1 / 2 F 1 ( z + w ) + F 2 ( z w )
where z = x + i y and w = x i y .
We find the QFI (see (3.2.17) in [16])
I 51 = ( y x ˙ x y ˙ ) ( x ˙ + i y ˙ ) + i 8 ( x ˙ i y ˙ ) 2 + i 1 z w F 1 ( z + w ) + + i 1 z w F 2 ( z w ) .
(6) α = 1 , β = i and γ = A = B = C = 0 . Then,
C a b = 2 y x ± i y x ± i y 2 i x
and Equation (32) becomes
( x i y ) V , x x V , y y + 2 y ± i x V , x y ± 3 i V , y + 3 V , x = 0
from which follows
V 26 = F 1 ( z ) r + F 2 ( z )
where F 2 = d F 2 d z and z = x ± i y .
We find the QFI (see (3.2.18) in [16])
I 61 = ( y x ˙ x y ˙ ) ( x ˙ ± i y ˙ ) i z V + i F 2 ( z ) .
(7) A B 0 , A B and α = β = γ = C = 0 . Then, C a b = A 0 0 B .
Equation (32) becomes
( A B ) V , x y = 0 V , x y = 0
which gives the separable potential
V 27 = F 1 ( x ) + F 2 ( y ) .
We find the irreducible QFIs (see (3.2.20) of [16])
I 71 a = x ˙ 2 + 2 F 1 ( x ) , I 71 b = y ˙ 2 + 2 F 2 ( y ) .
It can be shown that there are four special potentials of the potential Equation (55) which admit additional time-dependent QFIs and are superintegrable. These are
7a. The potential
V 271 = k 1 x + c 1 2 + k 2 y + c 2 2
admits the independent FIs
I 72 a = t 2 2 y ˙ 2 + t ( y + c 2 ) y ˙ t 2 k 2 ( y + c 2 ) 2 1 2 y 2 c 2 y I 72 b = t 2 2 x ˙ 2 + t ( x + c 1 ) x ˙ t 2 k 1 ( x + c 1 ) 2 1 2 x 2 c 1 x .
7b. The potential
V 272 = F 1 ( x ) + k 2 y + c 2 2
admits the FI I 72 a .
7c. The potential
V 273 = F 2 ( y ) + k 1 x + c 1 2
admits the FI I 72 b .
7d. The potential (see [8])
V 274 = λ 2 8 ( x 2 + y 2 ) λ 2 4 c 1 x + c 2 y k 1 ( x + c 1 ) 2 k 2 ( y + c 2 ) 2
admits the independent FIs.
I 73 a = e λ t x ˙ 2 + λ ( x + c 1 ) x ˙ λ 2 4 ( x + c 1 ) 2 + 2 k 1 ( x + c 1 ) 2 I 73 b = e λ t y ˙ 2 + λ ( y + c 2 ) y ˙ λ 2 4 ( y + c 2 ) 2 + 2 k 2 ( y + c 2 ) 2 .
In all the above relations, λ , c 1 , c 2 , k 1 , k 2 are arbitrary constants.
(8) C 0 and α = β = γ = 0 .
Then, C a b = A C C B and Equation (32) becomes
C ( V , y y V , x x ) + ( A B ) V , x y = 0 .
Solving Equation (56), we find the potential
V 28 = F 1 y + b 0 x + b 0 2 + 1 x + F 2 y + b 0 x b 0 2 + 1 x
where b 0 A B 2 C .
This potential admits the QFI
I 81 = A x ˙ 2 + B y ˙ 2 + 2 C x ˙ y ˙ + ( A + B ) V + 2 C b 0 2 + 1 ( F 1 F 2 ) .
We note that b 0 ( A , B , C ) . Here, A , B , C are parameters of the potential and therefore cannot be taken as independent parameters of the QFI.
For b 0 = 0 , we have A = B , V , y y V , x x = 0 and the potential reduces to
V 28 ( b 0 = 0 ) = F 1 ( y + x ) + F 2 ( y x )
which is the solution of the 1d-wave equation.
For the potential Equation (59) we find the QFI
I 82 = x ˙ y ˙ + F 1 ( y + x ) F 2 ( y x ) .
(9) A = 2 , C = ± i and α = β = γ = B = 0 . Then,
C a b = 2 ± i ± i 0
and Equation (32) becomes
i ( V , x x V , y y ) + 2 V , x y = 0 .
Solving Equation (60), we find the potential
V 29 = r 2 F 1 ( z ) + F 2 ( z )
where F 1 = d 2 F 1 d z 2 and z = x ± i y .
This potential admits the QFI (see (3.2.21) in [16])
I 91 = x ˙ ( x ˙ ± i y ˙ ) + V 29 + 2 z F 1 ( z ) 2 F 1 ( z ) .
Observe that for the trivial KT C a b = A δ a b the condition G , a = 2 C a b V , b gives
G , a = 2 A V , a G = 2 A V
for all potentials V ( x , y ) . Therefore, we recover the trivial result that all 2d-potentials V ( x , y ) admit the QFI
I = A ( x ˙ 2 + y ˙ 2 + 2 V ) = 2 A H .
Comparing with previous works, we see that the potentials V 21 a and V 28 are new. The potential V 274 ( c 1 = c 2 = 0 ) is mentioned in [8].

8.1. The Superintegrable Potentials

When a potential belongs to two of the above nine Class II cases simultaneously it is superintegrable (e.g., potentials (3.2.34)–(3.2.36) of [16]), because in that case the potential admits two more autonomous FIs other than the Hamiltonian. From the above results, we find the following Class II superintegrable potentials (see also in [14,15]).
(S1) The potential (see (3.2.34) in [16], case (b) in [8,14])
V s 1 = k 2 ( x 2 + y 2 ) + b x 2 + c y 2
where k , b , c are arbitrary constants.
This is of the form Equation (34) for d 1 = d 2 = 1 ,
F 1 y x = b y x 2 + c x y 2 , F 2 ( x 2 + y 2 ) = k 2 ( x 2 + y 2 ) + b + c x 2 + y 2
and also of the separable form Equation (55). Therefore, V s 1 admits the additional QFIs
(64) I s 1 a = ( y x ˙ x y ˙ ) 2 + 2 b y 2 x 2 + 2 c x 2 y 2 (65) I s 1 b = 1 2 x ˙ 2 + k 2 x 2 + b x 2 (66) I s 1 c = 1 2 y ˙ 2 + k 2 y 2 + c y 2 .
We note that V s 1 k = λ 2 4 , b = k 1 , c = k 2 , λ 0 , is the V 274 for c 1 = c 2 = 0 , and therefore admits also the time-dependent FIs I 73 a , I 73 b .
(S2) Potentials of the form Equations (43) and (55). Then we have to solve the systems of PDEs Equation (42) and V , x y = 0 . We find
V s 2 = k 1 2 ( x 2 + 4 y 2 ) + k 2 x 2 + k 3 y
and the QFIs
(68) I s 2 a = x ˙ ( y x ˙ x y ˙ ) k 1 y x 2 + 2 k 2 y x 2 k 3 2 x 2 (69) I s 2 b = 1 2 x ˙ 2 + k 1 2 x 2 + k 2 x 2 (70) I s 2 c = 1 2 y ˙ 2 + 2 k 1 y 2 + k 3 y .
where k 1 , k 2 , k 3 are arbitrary constants.
This is the superintegrable potential of case (a) in [14]. Note that the QFI I 3 a given in [14] is not correct. The correct is the I s 2 a of Equation (68) above.
We remark that the potential (3.2.35) in [16] is superintegrable only for b = 4 a in which case the potential becomes V s 2 for k 1 = 2 a , k 2 = c and k 3 = 0 .
(S3) Potentials of the form Equations (34) and (43). We solve the system of PDEs Equations (33) and (42). We find
V s 3 = k 1 x 2 + k 2 r + k 3 y r x 2
and the QFIs
(72) I s 3 a = ( y x ˙ x y ˙ ) 2 + 2 k 1 y 2 x 2 + 2 k 3 r y x 2 (73) I s 3 b = x ˙ ( y x ˙ x y ˙ ) + 2 k 1 y x 2 + k 2 y r + k 3 x 2 + 2 y 2 r x 2
where r 2 = x 2 + y 2 .
The superintegrable potential Equation (71) is symmetric ( x y ) to the superintegrable potential of case (c) of [14]. Indeed in order to find the superintegrable potential of [14] we simply consider the case leading to the potential of the form V 24 of Equation (43) for β = 1 instead of α = 1 .
We note that if we rename the constants in Equation (71) as k 1 = b + c , k 2 = a , k 3 = c b we recover the superintegrable potential (3.2.36) of [16]. Indeed, we have
V s 3 = a r + b r + y + c r y r .
(S4) If we substitute the solution Equation (43) of the PDE Equation (42) in the PDE Equation (45), we find that for
F 1 ( r + y ) = k 1 + k 2 r + y , F 2 ( r y ) = k 3 r y
both PDEs Equations (42) and (45) are satisfied simultaneously. Therefore, the potential (see case (d) in [14])
V s 4 = k 1 r + k 2 r + y r + k 3 r y r
is superintegrable with additional QFIs
I s 4 a = x ˙ ( y x ˙ x y ˙ ) + k 1 y r + k 3 ( r + y ) r y k 2 ( r y ) r + y r I s 4 b = y ˙ ( x y ˙ y x ˙ ) + G ( x , y )
where G , x + y V s 4 , y = 0 and G , y + y V s 4 , x 2 x V s 4 , y = 0 .
We note that in case (d) in [14] the corresponding QFIs I 2 d and I 3 d are not correct, because { H , I 2 d } 0 and { H , I 3 d } 0 . Moreover, this superintegrable potential is the case (E20) in [15], and it is not mentioned in the review [16].
In the following tables, we collect the results on Class II potentials with the corresponding reference to the review paper [16].
Integrable Potentials
PotentialRef. [16]LFIs and QFIs
V 21 = F 1 y x x 2 + y 2 + F 2 ( x 2 + y 2 ) (3.2.10) I 11 = ( y x ˙ x y ˙ ) 2 + 2 F 1 y x
V 21 a = k x 2 + y 2 + F 2 ( x 2 + y 2 ) - I 11 a = ( y x ˙ x y ˙ ) 2 + 2 k ( 1 ) y 2 x 2 + y 2
V 22 = F 1 ( u ) F 2 ( v ) u 2 v 2 ,
u 2 = r 2 + A + ( r 2 + A ) 2 4 A x 2 1 / 2 ,
v 2 = r 2 + A ( r 2 + A ) 2 4 A x 2 1 / 2
(3.2.7,8) I 21 = ( y x ˙ x y ˙ ) 2 + A x ˙ 2 + v 2 F 1 ( u ) u 2 F 2 ( v ) u 2 v 2
V 23 = F 1 ( u ) F 2 ( v ) u 2 v 2 ,
u 2 = r 2 + r 4 4 A ( x ± i y ) 2 1 / 2 ,
v 2 = r 2 r 4 4 A ( x ± i y ) 2 1 / 2
(3.2.7,12) I 31 = ( y x ˙ x y ˙ ) 2 + A ( x ˙ ± i y ˙ ) 2 + v 2 F 1 ( u ) u 2 F 2 ( v ) u 2 v 2
V 24 = F 1 ( r + y ) + F 2 ( r y ) r (3.2.15) I 41 = x ˙ ( y x ˙ x y ˙ ) + ( r + y ) F 2 ( r y ) ( r y ) F 1 ( r + y ) r
V 25 = w 1 / 2 F 1 ( z + w ) + F 2 ( z w ) ,
z = x + i y , w = x i y
(3.2.17) I 51 = ( y x ˙ x y ˙ ) ( x ˙ + i y ˙ ) + i 8 ( x ˙ i y ˙ ) 2 +
   + i 1 z w F 1 ( z + w ) +
   + i 1 z w F 2 ( z w )
V 26 = F 1 ( z ) r + F 2 ( z ) ,
F 2 = d F 2 d z , z = x ± i y
(3.2.18) I 61 = ( y x ˙ x y ˙ ) ( x ˙ ± i y ˙ ) i z V + i F 2 ( z )
V 27 = F 1 ( x ) + F 2 ( y ) (3.2.20) I 71 a = 1 2 x ˙ 2 + F 1 ( x ) , I 71 b = 1 2 y ˙ 2 + F 2 ( y )
V 28 = F 1 y + b 0 x + b 0 2 + 1 x +
   + F 2 y + b 0 x b 0 2 + 1 x , b 0 A B 2 C
- I 81 = A x ˙ 2 + B y ˙ 2 +
   + 2 C x ˙ y ˙ + ( A + B ) V + 2 C b 0 2 + 1 ( F 1 F 2 )
V 28 ( b 0 = 0 ) = F 1 ( y + x ) + F 2 ( y x ) - I 82 = x ˙ y ˙ + F 1 ( y + x ) F 2 ( y x )
V 29 = r 2 F 1 ( z ) + F 2 ( z ) , F 1 = d 2 F 1 d z 2 ,
z = x ± i y
(3.2.21) I 91 = x ˙ ( x ˙ ± i y ˙ ) + V 29 + 2 z F 1 ( z ) 2 F 1 ( z )
Superintegrable Potentials
PotentialRef. [16]LFIs and QFIs
V s 1 = k 2 ( x 2 + y 2 ) + b x 2 + c y 2 (3.2.34) I s 1 a = ( y x ˙ x y ˙ ) 2 + 2 b y 2 x 2 + 2 c x 2 y 2 ,
I s 1 b = 1 2 x ˙ 2 + k 2 x 2 + b x 2 , I s 1 c = 1 2 y ˙ 2 + k 2 y 2 + c y 2
- For k = 0 : I 72 a , I 72 b
where c 1 = c 2 = 0 , k 1 = b , k 2 = c
- For k = λ 2 4 0 : I 73 a , I 73 b
where c 1 = c 2 = 0 , k 1 = b , k 2 = c
V s 2 = k 1 2 ( x 2 + 4 y 2 ) + k 2 x 2 + k 3 y (3.2.35)
for k 3 = 0
I s 2 a = x ˙ ( y x ˙ x y ˙ ) k 1 y x 2 + 2 k 2 y x 2 k 3 2 x 2 ,
I s 2 b = 1 2 x ˙ 2 + k 1 2 x 2 + k 2 x 2 , I s 2 c = 1 2 y ˙ 2 + 2 k 1 y 2 + k 3 y
V s 3 = k 1 x 2 + k 2 r + k 3 y r x 2 (3.2.36) I s 3 a = ( y x ˙ x y ˙ ) 2 + 2 k 1 y 2 x 2 + 2 k 3 r y x 2 ,
I s 3 b = x ˙ ( y x ˙ x y ˙ ) + 2 k 1 y x 2 + k 2 y r + k 3 x 2 + 2 y 2 r x 2
V s 4 = k 1 r + k 2 r + y r + k 3 r y r - I s 4 a = x ˙ ( y x ˙ x y ˙ ) + k 1 y r + k 3 ( r + y ) r y k 2 ( r y ) r + y r ,
I s 4 b = y ˙ ( x y ˙ y x ˙ ) + G ( x , y )
V 271 = k 1 x + c 1 2 + k 2 y + c 2 2 (3.2.20) I 71 a , I 71 b ,
I 72 a = t 2 2 y ˙ 2 + t ( y + c 2 ) y ˙ t 2 k 2 ( y + c 2 ) 2 1 2 y 2 c 2 y ,
I 72 b = t 2 2 x ˙ 2 + t ( x + c 1 ) x ˙ t 2 k 1 ( x + c 1 ) 2 1 2 x 2 c 1 x
V 272 = F 1 ( x ) + k 2 y + c 2 2 (3.2.20) I 71 a , I 71 b , I 72 a
V 273 = F 2 ( y ) + k 1 x + c 1 2 (3.2.20) I 71 a , I 71 b , I 72 b
V 274 = λ 2 8 ( x 2 + y 2 ) λ 2 4 c 1 x + c 2 y
k 1 ( x + c 1 ) 2 k 2 ( y + c 2 ) 2 , λ 0
(3.2.20) I 71 a , I 71 b ,
I 73 a = e λ t x ˙ 2 + λ ( x + c 1 ) x ˙ λ 2 4 ( x + c 1 ) 2 + 2 k 1 ( x + c 1 ) 2 ,
I 73 b = e λ t y ˙ 2 + λ ( y + c 2 ) y ˙ λ 2 4 ( y + c 2 ) 2 + 2 k 2 ( y + c 2 ) 2

9. The Constraint L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a

The integrability condition of the constraint L b V , b , a = 2 L ( a ; b ) V , b λ 2 L a gives the PDE Equation (31).
As mentioned above in Section 6, in order to find new potentials from the PDE Equation (31) one (or both) of the conditions α = β = 0 and a 1 = a 3 must be relaxed. However, if we do find a new potential, this solution should satisfy also the remaining PDEs Equations (29) and (30) in order to admit the time-dependent QFI I 3 given in case Integral 3 of Theorem 1. New potentials which admit the QFI I 3 shall be referred as Class III potentials.
We note that the PB { H , I 3 } = I 3 t 0 . Therefore, to find a new integrable potential we should find a Class III potential admitting two independent FIs of the form I 3 , say I 3 a and I 3 b , such that { I 3 a , I 3 b } = 0 .
After relaxing one, or both, of the conditions α = β = 0 and a 1 = a 3 we found that the only non-trivial Class III potential is the superintegrable potential V 3 b = λ 2 2 ( x 2 + y 2 ) (see Section 7.3) found for α 0 or β 0 above. Therefore, there are no new Class III potentials.

10. Using FIs to Find the Solution of 2d Integrable Dynamical Systems

In this section we consider examples which show how one uses the 2d (super-)integrable potentials to find the solution of the dynamical equations.
(1) The superintegrable potential V 3 b = 1 2 k 2 ( x 2 + y 2 ) where k 0 .
We find the solution by using the time-dependent LFIs L 42 ± = e ± k t ( x ˙ k x ) and L 43 ± = e ± k t ( y ˙ k y ) . Specifically, we have
e k t ( x ˙ k x ) = c 1 + e k t ( x ˙ + k x ) = c 1 x ˙ k x = c 1 + e k t x ˙ + k x = c 1 e k t x ( t ) = c 1 2 k e k t c 1 + 2 k e k t .
Similarly for the LFIs L 43 ± we find
y ( t ) = c 2 2 k e k t c 2 + 2 k e k t .
Here, c 1 ± , c 2 ± are arbitrary constants.
(2) The integrable potential V 2 = c y + F ( x ) where F 0 .
Using the LFI L 31 = y ˙ + c t = c 1 , we find directly y ( t ) = c 2 t 2 + c 1 t + c 2 where c , c 1 , c 2 are arbitrary constants.
Using the QFI 2 Q 31 = x ˙ 2 + 2 F ( x ) = c o n s t = c 3 , we have
d x d t = ± 2 F ( x ) + c 3 1 / 2 d t = ± 2 F ( x ) + c 3 1 / 2 d x t = ± 2 F ( x ) + c 3 1 / 2 d x + c 0
where c 0 is an arbitrary constant. The inverse function of t = t ( x ) is the solution of the system. If the function F ( x ) is given, the solution can be explicitly determined.
(3) For the integrable potential V 27 = F 1 ( x ) + F 2 ( y ) by using the QFIs
I 71 a = 1 2 x ˙ 2 + F 1 ( x ) and I 71 b = 1 2 y ˙ 2 + F 2 ( y )
we find
t = c 1 2 F 1 ( x ) 1 / 2 d x + c 0 , t = c 2 2 F 2 ( y ) 1 / 2 d y + c 3
where c 0 , c 1 = 2 I 71 a , c 2 = 2 I 71 b , c 3 are constants.

11. Conclusions

Using Theorem 1 we have reproduced in a systematic way most known integrable and superintegrable 2d potentials of autonomous conservative dynamical systems. The method used being covariant it is directly applicable to spaces of higher dimensions and to metrics with any signature and curvature.
We have found two classes of potentials, and in each class we have determined the integrable and the superintegrable potentials together with their QFIs. As the general solution of the PDE Equation (28) is not possible, we have found the potentials due to certain solutions only. New solutions of this equation will lead to new integrable and possibly superintegrable 2d potentials.
It appears that the most difficult part in the application of Theorem 1 to higher dimensions and curved configuration spaces is the determination of the KTs. The use of algebraic computing is limited once one considers higher dimensions since then the number of the components of the KT increases dramatically. Fortunately, today new techniques in Differential Geometry have been developed (see, e.g., in [24,27,28,29,30]), especially in the case of spaces of constant curvature and decomposable spaces, which can help to deal with this problem.

Author Contributions

Performed the calculations and wrote a first version of the manuscript, A.M.; Provided the methodology, supervised and revised the paper, M.T.; Contributed into the discussion and interpretation of the results, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

A.P. acknowledges financial support of Agencia Nacional de Investigación y Desarrollo—ANID through the program FONDECYT Iniciación grant no. 11180126. Additionally, by Vicerrectoría de Investigación y Desarrollo Tecnológico at Universidad Catolica del Norte.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paliathanasis, A.; Leach, P.G.L. Nonlinear ordinary differential equations: A discussion on symmetries and singularities. Int. J. Geom. Methods Mod. Phys. 2016, 13, 1630009. [Google Scholar] [CrossRef] [Green Version]
  2. Arnold, V.I. Mathematical Methods of Classical Mechanics; Springer: Berlin/Heidelberg, Germany, 1989; pp. 272–284. [Google Scholar]
  3. Darboux, G. Sur un problème de mécanique. Arch. Neerl. Sci. 1901, VI, 371. [Google Scholar]
  4. Noether, E. Invariant Variation Problems. Transp. Theory Statist. Phys. 1971, 1, 186. [Google Scholar] [CrossRef] [Green Version]
  5. 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]
  6. Katzin, G.H.; Levine, J. Dynamical symmetries and constants of the motion for classical particle systems. J. Math. Phys. 1974, 25, 1460. [Google Scholar] [CrossRef]
  7. Kalotas, T.M.; Wybourne, B.G. Dynamical noether symmetries. J. Phys. A Math. Gen. 1982, 15, 2077. [Google Scholar] [CrossRef]
  8. Fris, J.; Mandrosov, V.; Smorodinsky, Y.A.; Uhlir, M.; Winternitz, P. On higher symmetries in quantum mechanics. Phys. Lett. 1965, 16, 354. [Google Scholar] [CrossRef]
  9. Whittaker, E.T. A Treatise on the Analytical Dynamics of Particles and Rigid Bodies, 4th ed.; Cambridge University Press: Cambridge, UK, 1959; pp. 334–335, (1st ed. 1904). [Google Scholar]
  10. Dorizzi, B.; Grammaticos, B.; Ramani, A. A new class of integrable systems. J. Math. Phys. 1983, 24, 2282. [Google Scholar] [CrossRef]
  11. Thompson, G. Darboux’s problem of quadratic integrals. J. Phys. A Math. Gen. 1984, 17, 985. [Google Scholar] [CrossRef]
  12. Sen, T. Lie symmetries and integrability. Phys. Lett. A 1987, 122, 327. [Google Scholar] [CrossRef]
  13. Tsamparlis, M.; Paliathanasis, A. Two-dimensional dynamical systems which admit Lie and Noether symmetries. J. Phys. A Math. Theor. 2011, 44, 175. [Google Scholar] [CrossRef] [Green Version]
  14. Rañada, M.F. Superintegrable n = 2 systems, quadratic constants of motion, and potentials of Drach. J. Math. Phys. 1997, 38, 4165. [Google Scholar] [CrossRef]
  15. Kalnins, E.G.; Kress, J.M.; Pogosyan, G.S.; Miller, W., Jr. Completeness of superintegrability in two-dimensional constant-curvature spaces. J. Phys. A Math. Gen. 2001, 34, 4705. [Google Scholar] [CrossRef]
  16. Hietarinta, J. Direct methods for the search of the second invariant. Phys. Rep. 1987, 147, 87. [Google Scholar] [CrossRef]
  17. Koenigs, G. Sur les géodésiques a intégrales quadratiques. A note. In Lecons sur la Théorie Générale des Surfaces; Darboux, G., Ed.; Chelsea Publishing Company: New York, NY, USA, 1972; Volume IV, pp. 368–404. [Google Scholar]
  18. Daskaloyannis, C.; Ypsolantis, K. Unified treatment and classification of superintegrable systems with integrals quadratic in momenta on a two-dimensional manifold. J. Math. Phys. 2006, 47, 042904. [Google Scholar] [CrossRef] [Green Version]
  19. Tsamparlis, M.; Mitsopoulos, A. Quadratic first integrals of autonomous conservative dynamical systems. J. Math. Phys. 2020, 61, 072703. [Google Scholar] [CrossRef]
  20. Karpathopoulos, L.; Tsamparlis, M.; Paliathanasis, A. Quadratic conservation laws and collineations: A discussion. J. Geom. Phys. 2018, 133, 279. [Google Scholar] [CrossRef] [Green Version]
  21. Hadler, A.K.; Paliathanasis, A.; Leach, G.L. Noether’s Theorem and Symmetry. Symmetry 2018, 10, 744. [Google Scholar]
  22. Stephani, H. Differential Equations: Their Solutions Using Symmetry; Cambridge University Press: New York, NY, USA, 1989. [Google Scholar]
  23. Sarlet, W.; Cantrijn, F. Generalizations of Noether’s theorem in classical mechanics. SIAM Rev. 1981, 23, 467. [Google Scholar] [CrossRef] [Green Version]
  24. Rani, R.; Edgar, S.B.; Barnes, A. Killing tensors and conformal Killing tensors from conformal Killing vectors. Class. Quant. Grav. 2003, 20, 1929. [Google Scholar] [CrossRef] [Green Version]
  25. Markakis, C. Constants of motion in stationary axisymmetric gravitational fields. MNRAS 2014, 441, 2974. [Google Scholar] [CrossRef] [Green Version]
  26. Adlam, C.M.; McLenaghan, R.G.; Smirnov, R.G. An orbit analysis approach to the study of superintegrable systems in the Euclidean plane. Phys. Atom. Nucl. 2007, 70, 486. [Google Scholar] [CrossRef]
  27. Kalnins, E.G.; Miller, W., Jr. Killing tensors and variable separation for Hamilton-Jacobi and Helmholtz equations. SIAM J. Math. Anal. 1980, 11, 1011. [Google Scholar] [CrossRef]
  28. Coll, B.; Ferrando, J.J.; Sáez, J.A. On the geometry of Killing and conformal tensors. J. Math. Phys. 2006, 47, 062503. [Google Scholar] [CrossRef] [Green Version]
  29. Crampin, M. Structural equations for a special class of conformal Killing tensors of arbitrary valence. Rep. Math. Phys. 2008, 62, 241. [Google Scholar] [CrossRef]
  30. Garfinkle, D.; Glass, E.N. Killing tensors and symmetries. Class. Quant. Grav. 2010, 27, 095004. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Mitsopoulos, A.; Tsamparlis, M.; Paliathanasis, A. Integrable and Superintegrable Potentials of 2d Autonomous Conservative Dynamical Systems. Symmetry 2020, 12, 1655. https://doi.org/10.3390/sym12101655

AMA Style

Mitsopoulos A, Tsamparlis M, Paliathanasis A. Integrable and Superintegrable Potentials of 2d Autonomous Conservative Dynamical Systems. Symmetry. 2020; 12(10):1655. https://doi.org/10.3390/sym12101655

Chicago/Turabian Style

Mitsopoulos, Antonios, Michael Tsamparlis, and Andronikos Paliathanasis. 2020. "Integrable and Superintegrable Potentials of 2d Autonomous Conservative Dynamical Systems" Symmetry 12, no. 10: 1655. https://doi.org/10.3390/sym12101655

APA Style

Mitsopoulos, A., Tsamparlis, M., & Paliathanasis, A. (2020). Integrable and Superintegrable Potentials of 2d Autonomous Conservative Dynamical Systems. Symmetry, 12(10), 1655. https://doi.org/10.3390/sym12101655

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