Low Dimensional Vessiot-Guldberg-Lie Algebras of Second-Order Ordinary Differential Equations

A direct approach to non-linear second-order ordinary differential equations admitting a superposition principle is developed by means of Vessiot-Guldberg-Lie algebras of a dimension not exceeding three. This procedure allows us to describe generic types of second-order ordinary differential equations subjected to some constraints and admitting a given Lie algebra as Vessiot-Guldberg-Lie algebra. In particular, well-known types, such as the Milne-Pinney or Kummer-Schwarz equations, are recovered as special cases of this classification. The analogous problem for systems of second-order differential equations in the real plane is considered for a special case that enlarges the generalized Ermakov systems.


Introduction
The theory of Lie systems, i.e., systems of non-autonomous first-order ordinary differential equations admitting a (generally nonlinear) superposition principle, is an old one and emerges principally from the pioneering work of Lie, Vessiot and Guldberg in the late 19th century [1].For a long time considered as a particular technique for differential equations, the importance of Lie systems in physical applications, particularly in the context of integrable systems, as well as in control theory, has motivated extensive studies on the subject in the last few decades (see, e.g., [2,3]) that have led to natural generalizations of the notion of Lie systems [4][5][6].These generalizations of the classical analytical formulation allow an elegant and effective description of Lie systems in terms of distributions in the sense of Fröbenius and the theory of fiber bundles, hence offering a much wider spectrum of applications, as well as their adaptation to quantum systems (the excellent treatise on the geometrical foundations of Lie systems [2] contains an extensive and updated list of references enumerating these applications).
In this work, we reconsider, mainly for the case of scalar second-order ordinary differential equations (scalar SODE Lie systems), a direct approach to Vessiot-Guldberg-Lie algebras in dimensions not exceeding three.Based on the particularities of the case of two-dimensional Vessiot-Guldberg-Lie algebras, the admissible realizations in terms of vector fields for three dimensions are obtained, enabling a direct construction of differential equations associated with these algebras.
We stress the fact that the multiple results concerning Lie algebras of vector fields could be used (see, e.g., [7]) and that the superposition problem has been studied from this perspective by various authors [8][9][10].The purpose of our direct construction is of a practical nature, in order to illustrate that cleverly-chosen constraints on realizations allow the construction of quite general differential equations (systems) admitting Vessiot-Guldberg-Lie algebras.Definition: A system of ordinary differential equations: is said to admit a fundamental system of solutions if the general solution of the system can be expressed in terms of a finite number m of functionally independent particular solutions {y 1 , • • • , y m } and n significant constants {C 1 , • • • , C n }: where We say that Equation ( 2) is a superposition rule for the system Equation (1), while the set {y 1 , • • • , y m } is called a fundamental system of solutions.The vector field in R n defined by: is further called the t-dependent vector field associated with the system Equation (1).It is immediate to verify that there is a one-to-one correspondence between systems and vector fields of this type.The basic criterion concerning the existence of fundamental systems of solutions is due to Lie himself [1] and can be formulated as follows [2,11]:

Lie-Scheffers Theorem [1]:
The system of ODEs Equation (1) possesses a generic fundamental system of solutions if it can be represented in the form: such that the vector fields: span an r-dimensional Lie algebra L VG .Moreover, the so-called Lie-condition: This criterion implies in particular that the t-dependent vector field Equation (3) can be written as: The Lie algebra L VG is usually called a Vessiot-Guldberg-Lie algebra of the system.The existence of such an algebra hence ensures the existence of a (nonlinear) superposition principle [2,12,13], which can (at least formally) be derived in terms of the differential invariants associated with the generators of a Vessiot-Guldberg-Lie algebra.
It shall however be observed that, in contrast to the Lie algebra of point symmetries of differential equations [11,12], a Vessiot-Guldberg-Lie algebra generally does not provide an invariant of the system.Actually, the same system Equation (1) may admit different superposition rules, hence non-isomorphic Vessiot-Guldberg-Lie algebras that can have different dimensions.The Riccati equation and other ODEs related to it constitute the canonical examples for this pathology [2].In this sense, the problem differs considerably from the construction of differential equations admitting certain types of symmetry algebras [14].

Vessiot-Guldberg-Lie Algebras with r ≤ 3 for Scalar SODE Systems
In this section, realizations of Vessiot-Guldberg-Lie algebras in low dimensions are analyzed directly for the case of a scalar second-order ODE, corresponding formally to the case of first-order systems with n = 2.
A scalar second-order ODE: ẍ − F (t, x, ẋ) = 0 (7) can always be rewritten as a first-order system: with n = 2.In this case, the t-dependent vector field X(t, x, v) has the form: We say that Equation ( 7) is a SODE Lie system if it possesses a Vessiot-Guldberg-Lie algebra L VG Geometrically, the vector field X(t, x, v) takes values in the Lie algebra L VG ; hence, in particular, for two given values t 0 = t 1 , the vector fields: must be elements of L VG .Considering the determinant of the matrix of coefficients of these vector fields, we find that: Now, if dim L VG = 1, then Equation (10) must vanish for all values of t, from which we conclude that: for all t ∈ R, implying that the constraint ∂F ∂t = 0 must be satisfied.The SODE Lie system thus reduces to: ẍ − F (x, ẋ) = 0 (11) and the generator of L VG is given by The substitution ẋ = p further reduces Equation (11) to a first-order equation, recovering the well-known property that autonomous second-order ODEs admit a superposition principle.The same conclusion holds for autonomous SODE Lie systems in higher dimensions [5,13].
The simplest case of a non-trivial Vessiot-Guldberg-Lie algebra is given by a Lie algebra in dimension two, as it is either Abelian or isomorphic to the affine algebra a 2 on the plane [15] (in the following, the reader is led to this reference for the elementary properties of Lie algebras; for classifications in low dimensions, see [15,16]).Suppose thus that the system Equation (8) admits a Vessiot-Guldberg-Lie algebra L VG of dimension two (it is implicitly assumed that it does not admit a smaller dimensional Lie algebra L VG ).By Equation (6), we have the identity: where X 1 , X 2 are two independent elements of L VG .Using again the vector fields Equation ( 9), there exist t 0 = t 1 , such that the determinant Equation ( 10) does not vanish: showing that the vector fields are linearly independent.Without loss of generality, we can define a basis of L VG as: We rewrite the generators of Equation ( 13) as: where, by Equation ( 12), it follows that for some non-constant function F 2 (t) and g (x, v) = 0.The commutator is given by: The only possibility that X 1 , X 2 generate a two-dimensional Lie algebra is that [X 1 , X 2 ] = λX 1 + µX 2 for some λ = 0.This assumption leads to the equations: Solving the first of these equations and inserting the result in the second, we obtain the general solution: A further linear change allows us to suppose that µ = 0.The t-dependent vector field is thus given by: ∂ ∂v and the associated SODE Lie system is easily seen to be an equation of the Liouville type (in various applications, this type of SODE Lie system is also known as the generalized Buchdahl equation [4]): This equation is well known to admit an integrating factor, a first integral being given by: In this sense, a two-dimensional Vessiot-Guldberg-Lie L VG does not provide new results to those already derived from the theory of differential equations.

dim L VG = 3
For the case of three-dimensional Vessiot-Guldberg-Lie algebras, the analysis of the generic form of the realizations of the algebra generators is quite similar to that already developed for the two-dimensional algebras.Suppose that Equation ( 7) does admit a Vessiot-Guldberg-Lie algebra with dim L VG = 3.By the same argument used in two dimensions, the vector fields: for some arbitrary values t 0 = t 1 .Further, the commutator: also belongs to the Vessiot-Guldberg-Lie algebra L VG .The requirement on the dimension of L VG implies that [X 1 , X 2 ] must be linearly independent on X 1 , X 2 , as otherwise L VG reduces to the affine Lie algebra a 2 .Therefore, there exists a third vector field } constitute a basis of L VG and satisfying the identity: for some constants α, β, γ, where necessarily γ = 0.As a consequence, a three-dimensional Vessiot-Guldberg-Lie algebra L VG can always be realized in the following way: In particular, the t-dependent vector field is given by For the latter realization, the remaining commutators are given by: As the component in ∂ ∂x of X 2 is zero, a change of basis always enables us to reduce the first commutator to the form: where µ = 0.By means of a scaling change in X 2 , we can further suppose that µ = 1 (we observe that [X 1 , X 2 ] = αX 2 cannot appear by the previously indicated reasons).The remaining brackets have the form: for some constants and a routine computation shows that the Jacobi identity is fulfilled whenever the relations b 2 = −a 1 and: are satisfied.At this stage, the problem corresponds essentially to classifying three-dimensional Lie algebras, a well-known problem solved in any standard reference [15].We merely indicate that it suffices to consider the following three cases: We observe that in the commutators Equation ( 16), the role of X 1 and X 2 is symmetric, implying that Cases (a) and (c) are equivalent.For each of the remaining possibilities, the commutators Equation ( 16) can be further reduced by successive changes of basis.As this is a routine computation, we skip the details and will only indicate the resulting three-dimensional Lie algebras.
If λ = 0, a similar reasoning with the Jacobi identity and the realization Equation (15) shows that only one isomorphism class can appear, corresponding to the following values of the constants in Equation ( 16): Proposition 1: If the vector fields Equation ( 15) generates a three-dimensional Vessiot-Guldberg-Lie algebra L VG , then it belongs to one of the following types: Concerning the isomorphism class of these Lie algebras, we have that r 1,1 and r 2 are isomorphic to the simple Lie algebra sl(2, R), while r 1,−1 is isomorphic to so(3).The remaining algebras are solvable.Following the list of [16], we have the isomorphisms r 3,λ A 3,8 , r 4 A 3,3 , r 5 A 3,8 and r 6,λ A 3,9 .Once the possible Lie algebras g have been deduced, the precise functions in the realization Equation (15), as well as the differential equations admitting g as Vessiot-Guldberg-Lie algebra L VG are obtained by solving the successive conditions imposed by the commutators of g.As an illustration of the procedure, we give the details for the Lie algebras r 1,ε , the remaining cases being obtained in a completely analogous way.
Let L VG have brackets (18) with ε = ±1.From the two first generators: it follows that the third generator has the form: In principle, four conditions on f (x, v) and g(x, v) corresponding to the components of the remaining commutators [X 1 , X 3 ] and [X 2 , X 3 ] must be solved.These are explicitly given by: As [X 2 , X 3 ] = X 1 holds, from the component in ∂ ∂x , we extract the first-order PDE: with generic solution: for some function G (x). Inserting g (x, v) into the commutators, the component in ∂ ∂x of the commutator [X 1 , X 3 ] further leads to the condition: The latter PDE has general solution:

d h dx
It turns out that with this choice of f (x, v), the commutator [X 2 , X 3 ] = X 1 is identically satisfied.It merely remains to consider the condition specified by the component in ∂ ∂v of [X 1 , X 3 ], given by: This constraint equation shows that the realization Equations ( 24) and (25) depends on an arbitrary function G(x).For any pair of functions {G (x) , k (x)} satisfying Equation (26), the commutators Equation ( 18) is satisfied and defines a three-dimensional Vessiot-Guldberg-Lie algebra for the (nonlinear) second-order ODE: Recall that for ε = 1, the Vessiot-Guldberg-Lie algebra is isomorphic to sl (2, R), while for ε = − 1, it is isomorphic to so (3).

Second-Order ODEs with Three-Dimensional Vessiot-Guldberg-Lie Algebra
For each of the Lie algebras g listed in Proposition 1, we enumerate the realization by vector fields, the possible constraints on the functions in their components and the generic SODE Lie system admitting g as a Vessiot-Guldberg Lie algebra. Realization: Second-order differential equation: Constraints: Realization: Second-order differential equation: Constraints: Realization: Realization: Realization: Realization:

∂ ∂v
Second-order differential equation: We observe that, in most cases, the preceding SODE Lie systems always involve a term quadratic in the velocity, with the exception of Equation ( 36), where more general functions of ẋ are allowed.

Examples
The preceding generic realizations of Vessiot-Guldberg-Lie algebras must in particular cover the already known SODE Lie systems for which the condition dim L VG ≤ 3 holds.
Among the low dimensional Lie algebras, the case r 2 sl (2, R) plays a special role, as it is related to various of the most relevant and best studied cases of SODE Lie systems (see, e.g., [2,4,5] and the references therein).
1.The well-known Milne-Pinney equation: appears making the choice G (x) = − 1 2 x in Equation (29).Here, the constraint Equation ( 30) is given by: The solution of this equation is immediate and equals k (x) = c x −3 .Now, setting F 2 (t) = −2F (t), we recover Equation (38).In particular, for c = 0 and F (t) = ω 2 (t), we get the harmonic oscillator with time-dependent frequency.The Lie algebra L VG is realized by the vector fields: ∂ ∂v We observe that the Milne-Pinney equation is the only SODE Lie system with Vessiot-Guldberg-Lie algebra r 2 , such that the differential equation does not contain a factor in ẋ2 .2. The Kummer-Schwarz equation: is also recovered from Equation (29).In this case, we take G (x) = x, leading to the constraint Equation (30): The solution is given by k (x) = c x 3 .The vector fields realizing the Lie algebra r 2 are: We observe that both equations actually arise as particular cases for the value a = 0 of the more general SODE Lie system: where a, b, m ∈ R and m = 1.

SODE Lie Systems in the Plane
The analogous problem of determining the systems of second-order ordinary differential equations: in the plane admitting a Vessiot-Guldberg-Lie algebra L VG is considerably more complicated, as the algebras are realized in R 4 with the coordinates {x 1 , x 2 , v 1 = ẋ1 , v 2 = ẋ2 }, and the constraints determined by the commutators lead to (nonlinear) partial differential equations that can only seldom be solved in full generality.In this sense, a direct approach as developed for scalar equations is not practical, unless some simplifications in the realizations are introduced.
Using the straightforward generalization of the realization Equation ( 14), it can be easily shown that the most general SODE Lie system admitting a two-dimensional Vessiot-Guldberg-Lie algebra L VG has the generic form: where we specifically have and with commutator [X 1 , X 2 ] = λX 1 (for λ = 0, the system is autonomous and clearly admits a superposition principle [2]).
This result suggests an ansatz based on realizations of the type Equations ( 18)-( 23) and such that the components of the vector field X 2 are linear functions in the coordinates {x 1 , x 2 , v 1 , v 2 }.We will only consider one case, namely that corresponding to the special case: for non-vanishing constants λ, µ.Without loss of generality, we can always suppose that µ = 1.Taking: the vector field X 3 = [X 1 , X 2 ] is given by: As the commutator [X 2 , X 3 ] does not have nonzero components in ∂ ∂x 1 and ∂ ∂x 2 , the only possibility is that: for some constant α.This means that with this realization, the possible Vessiot-Guldberg-Lie algebras are r 2 , r 4 and r 5 .Starting from the assumption Equation ( 42), the components in ∂ ∂v 1 and ∂ ∂v 2 imply the constraints: Introducing the auxiliary variable W 0 = λx 1 v 2 − x 2 v 1 , the solution to the system can be written as: With this choice of functions, it follows from the Jacobi identity that the bracket [X 1 , X 3 ] must have the following form: [X 1 , The term in X 2 is inessential, so that without loss of generality, we can suppose that γ = 0. Analyzing the components in ∂ , we see that the conditions: must be satisfied.This can only happen if ε = −α, from which we conclude that: With these values, it only remains to evaluate the components in ∂ ∂v 1 and ∂ ∂v 2 of [X 1 , X 3 ]: in terms of the variables {x 1 , x 2 , W 0 }, the conditions read: We thus obtain the functions: From the t-dependent vector field X (t, x, v) = X 1 + F 2 (t) X 2 and replacing the v i by the first derivatives ẋi , we obtain the following SODE Lie system admitting the Lie algebra generated by Equations ( 40) and (41) as the Vessiot-Guldberg-Lie algebra: From the commutators, we see at once that the algebra is isomorphic to r 2 = sl (2, R) if α = 0 and to r 4 if α = 0.In particular, for the special values λ = 1, α = 2 and F 2 (t) = −ω 2 (t), the preceding SODE Lie system reduces to: corresponding to the well-known generalizations of Ermakov systems (see, e.g., [5,17] and the references therein).We remark that a slight variation in the preceding argumentation also allows us to obtain SODE Lie systems having the Lie algebra r 5 as the Vessiot-Guldberg-Lie algebra.The system X (t, x, v) = X 1 + F 2 (t) X 2 in particular has the form: x

Conclusions
By means of a direct approach, scalar SODE Lie systems admitting a Vessiot-Guldberg-Lie algebra of a dimension at most three have been constructed.These equations encompass some of the well-known types of SODE Lie systems, like the Liouville-type equations and the Milne-Pinney and Kummer-Schwarz equations.The case of SODE Lie systems in the plane, much more cumbersome due to the complicated partial differential equations arising from the constraints, has been considered for a special case that however contains an important class, namely the generalized Ermakov systems [5].
The problem of determining Vessiot-Guldberg-Lie algebras corresponds formally to studying the realization of Lie algebras in terms of vector fields in R n , and thus, the ansatz could be formulated in terms of the primitive classes of realizations [7,8].Such an approach, combined with the lattice of subalgebras of a Lie algebra and the obstructions obtained for low dimensional Vessiot-Guldberg-Lie algebras, could eventually lead to a complete classification of scalar SODE Lie systems.From a pure practical perspective, however, it can be sometimes more convenient to choose a realization that generalizes some fixed realization, with the purpose of obtaining a class of equations or systems that naturally generalize already known types of SODE Lie systems.Proceeding along these lines also enables us to consider realizations depending on parameters, thus enabling the use of non-isomorphic Lie algebras simultaneously.In this context, an interesting problem concerns the combination of the notion of contractions of (realizations of) Lie algebras with that of Vessiot-Guldberg-Lie algebras, especially in connection with the recently developed notion of strong contractions [18].