Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Mirror Symmetry for New Physics beyond the Standard Model in 4D Spacetime
Previous Article in Journal
Frequency Shift of Photons in the Spacetime of Deformed RN BH
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Symmetries of Systems with the Same Jacobi Multiplier

by
Gabriel González Contreras
1,† and
Alexander Yakhno
2,*
1
Cátedra CONACyT, Universidad Autónoma de San Luis Potosí, San Luis Potosí 78000, Mexico
2
Departamento de Matemáticas, Universidad de Guadalajara, Guadalajara 44430, Mexico
*
Author to whom correspondence should be addressed.
Current address: NanoScience Technology Center and Department of Physics, University of Central Florida, Orlando, FL 32826, USA.
Symmetry 2023, 15(7), 1416; https://doi.org/10.3390/sym15071416
Submission received: 17 June 2023 / Revised: 13 July 2023 / Accepted: 13 July 2023 / Published: 14 July 2023
(This article belongs to the Section Mathematics)

Abstract

:
The concept of the Jacobi multiplier for ordinary differential equations up to the second order is reviewed and its connection with classical methods of canonical variables and differential invariants is established. We express, for equations of the second order, the Jacobi multiplier in terms of integrating factors for reduced equations of the first order. We also investigate, from a symmetry point of view, how two different systems with the same Jacobi multiplier are interrelated. As a result, we determine the conditions when such systems admit the same two-dimensional Lie algebra of symmetries. Several illustrative examples are given.

1. Introduction

Developing methods for finding solutions of ordinary differential equations (ODEs) has been a problem of interest over many centuries, with significant contributions of great mathematicians such as Euler, Jacobi, Lie, Poincaré and many others. It is very well known that the existence of sufficient number of first integrals (constants of motion) allows one to find the general solution of a given system of ordinary differential equations or reduce their order.
One way to determine a first integral is related with the concept of a group of continuous transformations and their invariants. The existence of such a group, admitted by given ODEs, reveals symmetry properties of differential equations and facilitate the process of integration. The foundations of the symmetry group analysis are contained in the works of Sophus Lie [1,2]. The actual formulation of the methods of modern group analysis of differential equations is described in various classical textbooks such as [3,4,5,6] and many others.
At present, symmetry methods are widely used to construct particular solutions (and sometimes general ones) of both linear and nonlinear equations. Based on the admitted symmetry algebras, one can analyze the integrability, classify the parameters included in the equations, describe conservation laws, search for coordinate systems in which second-order linear equations admit the separation of variables, etc. In the present work, we focus on the well-known but not completely studied concept of the Jacobi multiplier and its application to ODEs.
In Ref. [7], Jacobi introduced the concept of a multiplier M (now called the Jacobi multiplier), which satisfies a first-order partial differential equation, associated with a system of n first-order ODEs for n unknown functions. Moreover, he introduced what is now called the principle of the last multiplier μ , which states the following: if one uses n 1 first integrals to express n 1 functions, then the last first integral can be expressed in quadratures with the help of an integrating factor (multiplier) μ , that satisfies a similar equation for M.
The generalization and the basic properties of the Jacobi multiplier M and its relation to the concept of symmetry of a system of ODEs are described in [8] (see also [9]). Nucci, Leach and co-authors revived the interest of the Jacobi multiplier in the context of group analysis as a useful tool for calculating symmetries and their corresponding Lagrangians. For example, in Ref. [10], the last multiplier is used to calculate nonlocal symmetries for the Euler–Poinsot system; in Ref. [11], the first integrals of the Kepler problem are obtained by the method of Jacobi’s last multiplier; in a series of papers [12,13,14,15], the connection of the Jacobi multiplier with Lagrangians is analyzed for various mechanical systems. Let us mention also two recent works given in Ref. [16] and Ref. [17]. The first one provides the survey of the geometric theory of Jacobi multipliers in terms of jet manifolds in the spirit of Vinogradov’s school [18,19]. The second paper gives and up-to-date review of the Jacobi multiplier from a geometric perspective and shows some of its physical applications, such as finding constants of motion, which generalizes the so-called Hojman symmetry. Because the literature is vast, we recommend to the interested reader Ref. [20] and references therein for a detailed explanation of the use and applications of the Jacobi last multiplier.
The main objective of this work is to study a set of different Euler–Lagrange ODEs of the second order with different Lagrangians but with the same Jacobi multiplier and analyze the relation of their corresponding symmetries for such ODEs. The article is organized as follows: First, in Section 2, we determine the relationship between the symmetries of two first-order ODEs with the same integrating factor. For the case of a first-order ODE, the Jacobi multiplier is just the integrating factor. This section is more methodological in nature, but it is important for a more immediate understanding. In Section 3, we review the definition, main properties and connection of the Jacobi multiplier with the method of differential invariants for integrating a second-order ODE. Moreover, we establish the relation of the Jacobi multiplier for second-order ordinary differential equations, with integrating factors of two first-order ODEs, expressed in terms of invariants of an admitted symmetry. In this sense, the Jacobi multiplier can be considered as an integrating factor of a second-order ODE. Moreover, we consider the relation between the Lagrangian, obtained by means of the Jacobi multiplier, and the Noether theorem. Finally, for systems with the same Jacobi multiplier but with different Lagrangians, we find the necessary conditions for which two different ODEs share the same symmetries. We summarize our conclusions in Section 4.

2. First-Order Ordinary Differential Equation

2.1. Symmetries and Integrating Factor

First, we review some well-known facts from the theory of point symmetries for ordinary differential equations (ODEs) of the first order given by
β ( x , y ) d x α ( x , y ) d y = 0 d x α = d y β
or in normal form
F : = y f ( x , y ) = 0 , f ( x , y ) : = β / α .
Let us remember that the problem of solving the ODE (1) is equivalent to determining an invariant of the operator A : = α x + β y , i.e., a function ω = ω ( x , y ) such that
A ( ω ) = 0 ,
then ω = const . is the integral of (1). We use t to denote t for the rest of the paper.
The standard way to find the point symmetry of (1) [3]
S = ξ ( x , y ) x + η ( x , y ) y
is to apply the first prolongation of S
S 1 = S + ζ y , ζ = D x ( η ) y D x ( ξ ) , D x : = x + y y
to equation F = 0
S 1 ( F ) F = 0 = 0 ,
and to solve the first-order linear partial equations for the unknown functions ξ and η .
Taking into account that D x F = 0 = A / α , criterion (2) takes the form
S ( α ) A ( ξ ) α = S ( β ) A ( η ) β .
If one considers the relation for the commutator [5,9]
[ S , A ] = S A A S = λ ( x , y ) A ,
then we obtain Equation (3) exactly, excluding the function λ . Obviously, there is a trivial symmetry of the form S 0 = δ ( x , y ) A for any smooth function δ , because the above equality (4) holds for λ = A ( δ ) .
If one knows a nontrivial symmetry S, then the application of (4) to integral ω gives
S ( A ( ω ) ) A ( S ( ω ) ) = A ( S ( ω ) ) = λ A ( ω ) = 0 ;
therefore, S ( ω ) = h ( ω ) for some function h, which represents the well-known property of an admitted symmetry to transform a solution of a given equation into another solution of the same equation. Thus, the integral ω satisfies the following algebraic system of equations for ω , x : = x ( ω ) , ω , y : = y ( ω ) , respectively, given by:
A ( ω ) = α ω , x + β ω , y = 0 , S ( ω ) = ξ ω , x + η ω , y = h ( ω ) .
The determinant of this system Δ : = β ξ α η is nonzero; therefore, we can write
d ω = h ( ω ) β Δ d x α Δ d y = 0 .
If we demand μ : = 1 / Δ to be the integrating factor for (1)
β Δ , y = α Δ , x
we arrive at expression (3). Thus, to obtain a symmetry S, one can use the equivalent relations: (2), (4) or (5). This connection reflects the well-known fact that in general, the problem of finding a point symmetry of (1) has the same degree of difficulty as the problem of finding its integrating factor.
If an integrating factor μ for (1) is known, then expressing ξ = Δ / β + α η / β , we obtain
S = Δ / β + α η / β x + η y = Δ β x + η β A = : S x + S 0 ,
and for η = β ξ / α Δ / α , by analogy, we have
S = ξ x + β ξ / α Δ / α y = Δ α y + ξ α A = : S y + S 0 ,
where we have taken into account that η β A and ξ α A are trivial symmetries S 0 . Thus, any equation of (1) with integrating factor μ = μ ( x , y ) admits the symmetry (which we call the canonical symmetry)
S x = ( β μ ) 1 x or S y = ( α μ ) 1 y .

2.2. Symmetries of Two ODEs with the Same Integrating Factor

Let us now consider another ODE different from (1), but with the same integrating factor μ ( x , y )
β ˜ ( x , y ) d x α ˜ ( x , y ) d y = 0 ,
where α β ˜ β α ˜ 0 . Then, the relation between the canonical symmetry S x for (1) and the symmetry S ˜ x = ( β ˜ μ ) 1 x for (6) is very simple
S ˜ x = β β ˜ S x .
Due to the presence of the trivial symmetry S 0 = δ ( x , y ) A , it is easy to show that any two equations with the same integrating factor actually have the same symmetry. Indeed, from the following equality
S = S x + δ A = S ˜ x + δ ˜ A ˜ ,
we obtain
δ ˜ β ˜ = δ β , δ = μ 1 1 β ˜ / β α β ˜ β α ˜ .
This result, to the best of our knowledge, has not been reported in the literature.
For example, let us consider two equations. The first one is with homogeneous coefficients
2 d x x 2 y 2 d y = 0 ,
and the second one is the linear equation
cos x 2 y x d x + d y = 0 .
Both equations have the same integrating factor μ = x 2 ; therefore, their canonical symmetries are given by:
S x = x 2 2 x , S ˜ x = x 3 x cos x 2 y x .
It is easy to see, using (7) and (8), that Equations (9) and (10) admit the same following symmetry
S = x 2 ( x 2 + y 2 ) x + x y 2 ( 2 x + 2 y x cos x ) y 2 y 2 2 x y + x 2 cos x ,
which gives the same integrating factor μ for both equations.

3. Second-Order Ordinary Differential Equation

3.1. Symmetries and the Jacobi Multiplier

Let us now consider a second-order differential equation given in the following form
x ¨ = f ( t , x , x ˙ )
for the unknown function x = x ( t ) . Introducing a new function v : = x ˙ = d x / d t , we obtain the system F : = ( F 1 , F 2 ) = 0 of two first-order equations:
F 1 : = x ˙ v = 0 , F 2 : = v ˙ f ( t , x , v ) = 0 ,
or in its symmetric form
d t = d x v = d v f .
We can consider the following operator A : = t + v x + f v , similarly to what we had for the first-order ODEs. The general integral of (12) is given by a set of two functionally independent first integrals ω i = ω i ( t , x , v ) = c i (constants of motion), such that
d ω i F = 0 = 0 , | J x v | : = ( ω 1 , ω 2 ) ( x , v ) 0 .
From (13) it follows that A ( ω i ) = 0 , i.e., functions ω i are invariants of the operator A. Since any other solution Ω of (12) is an invariant of A and is a function of ω 1 and ω 2 , the Jacobian reads
( Ω , ω 1 , ω 2 ) ( t , x , v ) = 0 ,
i.e.,
J x v Ω , t J t v Ω , x + J t x Ω , v = 0 , | J t v | : = ( ω 1 , ω 2 ) ( t , v ) , | J t x | : = ( ω 1 , ω 2 ) ( t , x ) .
Comparing the above relation with A ( Ω ) = 0 , and denoting M : = J x v , we obtain
M v = J t v , M f = J t x .
The function M is called the Jacobi multiplier [9]. Due to the relation of M with Jacobians, M is a solution of the equation
M , t + ( M v ) , x + ( M f ) , v = 0 .
The form of M depends on the selection of ω 1 , 2 . If one has ω ˜ 1 , ω ˜ 2 as an independent first integrals, then function M ˜ = M ( ω ˜ 1 , ω ˜ 2 ) ( ω 1 , ω 2 ) is a Jacobi multiplier, and in general, any function of the form M Ω is a Jacobi multiplier. In particular, any multiplier M is defined up to an arbitrary nonzero constant factor.
Any point symmetry of system (12) has the form
S = ξ ( t , x , v ) t + η 1 ( t , x , v ) x + η 2 ( t , x , v ) v
and represents the so-called contact symmetry for Equation (11) if ξ or η 1 depends on v (with the nonlinear dependence of η 1 on v if ξ , v = 0 ). The application of the first prolongation of S
S 1 = S + ζ 1 x ˙ + ζ 2 v ˙ , ζ 1 = D t ( η 1 ) x ˙ D t ( ξ ) , ζ 2 = D t ( η 2 ) v ˙ D t ( ξ ) , D t = t + x ˙ x + v ˙ v
to system F = 0
S 1 F 1 F = 0 = 0 , S 1 F 2 F = 0 = 0
gives two linear equations for determining ξ , η 1 and η 2 :
η 2 = A ( η 1 ) v A ( ξ ) , S ( f ) = A ( η 2 ) f A ( ξ ) ,
where we have taken into account that D t F = 0 = t + v x + f v = A .
Excluding the function λ from the commutator relation
[ S , A ] = λ ( t , x , v ) A
leads to (16). Once again, there is a trivial symmetry S 0 = δ ( t , x , v ) A . For simplicity, one can introduce the following operator (evolutionary form of symmetry [6]):
S ^ : = S ξ A = : η ^ 1 x + η ^ 2 v .
The function η ^ 1 is called in Ref. [18] a generating function of symmetry. For the operator S ^ , we have
[ S ^ , A ] = 0 η ^ 2 = A ( η ^ 1 ) , S ^ ( f ) = A ( η ^ 2 ) , η ^ 1 0 .
Thus, to determine the unknown function η ^ 1 , we have to solve one second-order linear partial differential equation (defining equation) given by
S ^ ( f ) = A ( A ( η ^ 1 ) ) .
As in the case of the first-order ODE, there is an infinite number of point symmetries S ^ for system (12) and in general, the problem of finding them has the same degree of difficulty as solving the given system. Of course, if we make an additional supposition about the form of the symmetries or about the form of f ( t , x , v ) , then the problem of determining the admitted symmetries is drastically simplified. For example, in order to find all admitted point symmetries of Equation (11), we can assume that η ^ 1 = a ( t , x ) v + b ( t , x ) , which enables us to split Equation (18) with respect to v. In the case of the so-called Hojman symmetries [17,21], when f , v = 0 , Equation (18) is simplified to A ( A ( η ^ 1 ) ) η ^ 1 f , x = 0 . The derivation of the above equation with respect to v leads to the following expression A ( η ^ , x 1 + η ^ , v 2 ) = 0 , which gives the first integral div S ^ : = η ^ , x 1 + η ^ , v 2 = const . The simplest case f = 0 leads to the general solution η ^ 1 = a ( v , x t v ) t + b ( v , x t v ) .
Now, let us consider the general case. Let ω = ω ( t , x , v ) be a nontrivial invariant of the symmetry operator S ^ , i.e.,
S ^ ( ω ) = η ^ 1 ω , x + η ^ 2 ω , v = 0 , ω t , ω const . ;
then, the general invariant of S ^ has the form I = I ( t , ω ) . Operator A is a symmetry of the above equation, because [ A , S ^ ] = 0 , therefore A ( ω ) = : h ( t , ω ) . If h 0 , then ω 1 = ω . If h 0 , then taking into account that A ( ω ) = d ω / d t , we have the first-order equation
d ω d t = h ( t , ω ) .
If one is able to determine the general integral of the above equation, then one obtains the first integral ω 1 = ω 1 ( t , ω ) = c 1 of (12).
Let us note that S ^ ( ω 1 ) = ω , ω 1 S ^ ( ω ) = 0 , that is, ω 1 obtained from (20) is an invariant of symmetry S ^ . Moreover, ω , v 1 0 (if ω , v 1 = 0 , then ω , v = 0 and from (19), ω , x = 0 , which gives a trivial invariant, or η ^ 1 = 0 and so from (17), η ^ 2 = 0 ); therefore, we can relate v with a new variable z : = ω 1 ( t , x , v ) = c 1 . Using z, we can rewrite system (13) with only one first-order equation
v ( t , x , z ) d t d x = 0 ,
and the operator A takes the form A z = t + v ( t , x , z ) x . Thus, the nontrivial symmetry S ^ allows us to reduce the order of the system (12). Let us note that for Equation (11) the described procedure is known as a method of differential invariants, because the invariant ω (differential invariant of the first order) depends on the derivative v = x ˙ , and A ( ω ) is the differential invariant of the second order.
The integrating factor μ = μ ( t , x , z ) of (21) satisfies the equation μ , t + μ v ( t , x , z ) , x = 0 similar to (14); therefore, μ is a Jacobi multiplier, called Jacobi’s last multiplier, because it permits us to determine the last integral ω 2 = c 2 . Finally, the knowledge of an invariant of symmetry S ^ and Jacobi’s last multiplier μ allows us to solve the system (12) in quadratures.
On the other hand, if it is possible to determine a function u ( t , x , v ) such that S ^ ( u ) = 1 , then the Jacobian of transformation u = u ( t , x , v ) , ω = ω ( t , x , v ) , reads as
| J u ω | : = ( u , ω ) ( x , v ) = u , x ω , v u , v ω , x = ω , v η ^ 1 0 ,
where the functions ( u , ω ) can be taken as new variables. Moreover, S ^ ( A ( u ) ) = A ( S ^ ( u ) ) = 0 ; therefore, A ( u ) = : a ( t , ω ) and A u ω = t + a ( t , ω ) u + h ( t , ω ) ω , i.e., in new variables (called canonical coordinates), the system F = 0 admits the symmetry S ^ u ω = u , which represents a translation in u. If one knows the first integral ω 1 ( t , ω ) = c 1 obtained from (20), then u a ( t , ω ( t , c 1 ) ) d t = ω 2 = c 2 is the second first integral. As an example, let us consider [22] (it is a particular case of the class of ODEs from [23]). A = t + v x + v t 2 + t + v x 2 t 1 v , which admits the symmetry S ^ = exp t ( t 2 + t + v ) / x x + t 2 + t + v x v . The canonical variables are u = x e t ω , ω = ( t 2 + t + v ) / x , and A u ω = t t ( t + 1 ) e t ω u ; therefore,
ω 1 = ω = c 1 , ω 2 = x e c 1 t + t ( t + 1 ) e c 1 t d t = c 2 .
It is interesting to note, that in this particular example, there is no symmetry for A with a linear dependence on v, that is, Equation (11) corresponding to A is completely integrable but does not admit any point symmetry, only the contact one. The importance of contact (dynamical) symmetries for the integration of ODEs of the second order is discussed, for example, in [23].

3.2. Relations between Jacobi Multiplier, Last Jacobi Multiplier and Symmetries

Incidentally, there is a simple relation between the Jacobi multiplier and the coefficients of two nonproportional symmetries S ^ 1 : = η ^ 11 x + η ^ 12 v and S ^ 2 : = η ^ 21 x + η ^ 22 v admitted by system (12).
Let us consider Δ = η ^ 11 η ^ 22 η ^ 12 η ^ 21 0 . Then, [9] (par. 152)
Δ M = η ^ 11 η ^ 12 η ^ 21 η ^ 22 J x v = η ^ 11 η ^ 12 η ^ 21 η ^ 22 ω , x 1 ω , x 2 ω , v 1 ω , v 2   = S ^ 1 ( ω 1 ) S ^ 1 ( ω 2 ) S ^ 2 ( ω 1 ) S ^ 2 ( ω 2 ) = S ^ 1 ( ω 1 ) S ^ 2 ( ω 2 ) S ^ 1 ( ω 2 ) S ^ 2 ( ω 1 ) .
Any admitted symmetry transforms a solution of the given system into another solution of the same system, that is S ^ i ( ω j ) = h i j ( ω 1 , ω 2 ) , and we obtain
Δ M = H ( ω 1 , ω 2 ) , A ( H ) = 0 .
Therefore, from the property that M Ω is a Jacobi multiplier if A ( Ω ) = 0 , function Δ 1 is a Jacobi multiplier.
Using the variables t, x, z for Equation (21), the symmetries S ^ 1 , S ^ 2 take the following form
Z ^ 1 : = S ^ 1 ( x ) x , Z ^ 2 : = S ^ 2 ( x ) x + S ^ 2 ( z ) z ,
and the integrating factor (last Jacobi multiplier) reads as follows:
μ 1 ( t , x , z ) = S ^ 1 ( x ) 0 S ^ 2 ( x ) S ^ 2 ( z ) = S ^ 1 ( x ) S ^ 2 ( z ) = : Δ z , Δ z = ω , v 1 Δ .
Let S ^ 1 and S ^ 2 be two nonproportional symmetries of system (12), forming the basis of two-dimensional Lie algebra. If the Lie algebra is not Abelian, then it is possible to determine S ¯ i = a i S ^ 1 + b i S ^ 2 , a i , b i R such that [ S ¯ 1 , S ¯ 2 ] = S ¯ 1 . Here, and after, we assume that the Lie algebra S ^ 1 , S ^ 2 is Abelian or [ S ^ 1 , S ^ 2 ] = S ^ 1 .
Let us start to reduce the order of system (12) using the symmetry S ^ 1 , i.e., selecting its invariant S ^ 1 ( ω ) = 0 , then expressing operator A in terms of the variables t, x and ω , we obtain
A ω : = t + v ( t , x , ω ) x + h ( t , ω ) ω ,
and the given system is reduced to the following one
d t = d x v ( t , x , ω ) = d ω h ( t , ω ) .
Note that Equation (20) is included in the above relations. The second symmetry S ^ 2 transforms to S ^ 2 ω : = S ^ 2 ( x ) x + S ^ 2 ( ω ) ω and obviously commutes with A ω . Moreover, S ^ 2 ( ω ) 0 (else Δ = 0 ). Considering operator S ^ 2 ω in variables t and ω and applying all the facts from Section 2 for ODE h ( t , ω ) d t d ω = 0 , we obtain its integrating factor μ ω : = 1 / S ^ 2 ( ω ) .
If the symmetries S ^ i commute, [ S ^ 1 , S ^ 2 ] = 0 , then S ^ 2 is the symmetry of equation S ^ 1 ( ω ) = 0 , so S ^ 2 ( ω ) = h 2 ( t , ω ) . In the case when [ S ^ 1 , S ^ 2 ] = S ^ 1 , we have [ S ^ 1 , S ^ 2 ] ( ω ) = S ^ 1 ( S ^ 2 ( ω ) ) = 0 , that is, function S ^ 2 ( ω ) is an invariant of operator S ^ 1 ; therefore, once again, S ^ 2 ( ω ) = h 2 ( t , ω ) . In any case, μ ω does not depend on x explicitly. Therefore, the first integral does not either: ω 1 = ω 1 ( t , ω ) = z . Using the variables t, x and v ( t , x , z ) the symmetry S ^ 1 takes the form S ^ 1 ( x ) x and produces the integrating factor μ x : = 1 / S ^ 1 ( x ) = 1 / η ^ 11 for Equation v ( t , x , z ) d t d x = 0 . Collecting all integrating factors, we have the following relations between them:
Δ 1 = ω , v μ x μ ω , μ = μ x μ ω / ω , ω 1 .
To the best of our knowledge, the above relation has not been reported in the literature. Relations (22) allow us to write down Equation (20) as a differential one-form [23]
d ω 1 = Δ 1 d t d x d v 1 v f 0 η 1 η 2 = Δ 1 ω , ω 1 A ( ω ) d t + d ω = 0 .
Taking into account (22), the above equation is equivalent to μ ω ( h ( t , ω ) d t d ω ) = 0 and is therefore integrable.
Let us note that in canonical coordinates ( u , ω ) , such that S ^ 1 ( u ) = 1 , S ^ 1 ( ω ) = 0 , we have
Δ u ω = 1 0 S ^ 2 ( u ) S ^ 2 ( ω ) = S ^ 2 ( ω ) = 1 / μ ω = | J u ω | Δ .
The quantity S 1 S 2 : = ξ 1 η 2 η 1 ξ 2 for two point symmetries S 1 = ξ 1 ( t , x ) t + η 1 ( t , x ) x and S 2 = ξ 2 ( t , x ) t + η 2 ( t , x ) x of Equation (11) is used in the general scheme for the integration in quadratures of such ODEs [1,9] with Lie algebra S 1 , S 2 . In [24], expression S 1 S 2 is called the pseudo scalar product. It is easy to see that Δ = ( S 1 S 2 ) f + η ^ 11 ζ 12 η ^ 21 ζ 11 , with Δ formed by two nonproportional canonical symmetries S ^ 1 , 2 corresponding to S 1 , 2 , and ζ 1 i are coefficients of the first prolongation of S i .
Another interesting relation between the Jacobi multiplier M and symmetry S ^ was pointed out by Bianchi [9] (par. 152) and is given by
A ( S ^ ( ln M ) + div S ^ ) = 0 .
In the case of a divergence-free symmetry S ^ , one obtains S ^ ( ln M ) = Ω ( ω 1 , ω 2 ) . Using the well-known property of a symmetry to transform any solution of a given system to a solution of the same system, i.e., S ^ ( H ( ω 1 , ω 2 ) ) = Ω ( ω 1 , ω 2 ) we conclude that A ( M ) = 0 .
As an example, let us consider A = t + v x + 2 v ( v + 1 ) / ( x t ) v with two commuting symmetries S ^ 1 = v + 1 t x x + v + 1 x t v , S ^ 2 = ( 1 v ) x + 2 v + v 2 t x v and Δ = ( v + 1 ) 3 / ( t x ) 2 . It is easy to see that ω = ( v + 1 ) / ( x t ) is an invariant of S ^ 1 and S ^ 2 ( ω ) = ω 2 . Therefore, μ ω = ω 2 is the integrating factor of Equation (20) for h ( t , ω ) = A ( ω ) = ω 2 . The integration gives the first integral ω 1 = t + ω 1 . For Equation (21) with v ( t , x , z ) = x z z t , the integrating factor is μ x = 1 / S ^ 1 ( x ) = t z . Let us note that the last Jacobi multiplier μ = 1 / Δ z in this case coincides with μ x , because S ^ 2 ( z ) = 1 .
If one knows one symmetry S ^ 1 and Δ , then, to obtain the second admitted symmetry S ^ 2 , we have the following three equations:
η ^ 22 = A ( η ^ 21 ) , Δ = η ^ 11 η ^ 22 η ^ 12 η ^ 21 ,
S ^ 2 ( f ) = A ( η ^ 22 ) ,
but it is easy to verify that Equation (24) holds if (23) is valid. Therefore, from (23), we obtain one equation of the first order to determine η ^ 21
A η ^ 21 η ^ 11 = Δ η ^ 11 2
instead of an equation of the second order (18).

3.3. Jacobi Multiplier and Noether Theorem

It is known [18] that the Euler–Lagrange equation gives a system of differential equations for some Lagrangian if and only if the so-called operator of universal linearization of the system is autoconjugated (see [19] Chapter 5, par. 4). Noether’s first theorem establishes a connection between the conservation laws of equations arising from the variational principle with Lagrangian symmetries. Until now, the use of this theorem has provided a practical tool for calculating the conservation laws of the Euler–Lagrange equations. In the case when a given equation is not related with a Lagrangian, one can use the concept of the generating function of conservation laws [18]. Since any second-order ODE can be seen as the Euler–Lagrange equation for some Lagrangian, we now study the relation between the Jacobi multiplier and Noether’s theorem.
Suppose we have the following Lagrangian L = L ( t , x , v ) ; then, Equation (11) is obtained by means of the Euler–Lagrange equation
L , x D t ( L , v ) = 0 ,
where operator D t is the operator given in (15). For the system (13), we have
E ( L ) = 0 , E : = x A v ;
in the case when L , v v 0 , we can obtain an expression for f ( t , x , v ) in terms of the Lagrangian
f ( t , x , v ) = L , x L , t v v L , x v L , v v .
The derivation of the above expression with respect to v leads to Equation (14) [25] (par. 123) with M = L , v v . That is, in the case of a nonzero Jacobi multiplier, one can obtain a Lagrangian by quadratures
L M : = M d v d v + Q ( t , x ) v + R ( t , x ) ,
with
Q , t R , x = M , x d v d v A M d v ,
where the right side of the above relation does not depend on v due to (14).
For two nonproportional symmetries S ^ 1 , 2 , the function Δ 1 is a Jacobi multiplier, so one can take M = Δ 1 0 and obtain Lagrangian L M at least in quadratures. The relaxed condition for a point symmetry S ^ 1 to be a variational symmetry, i.e., S ^ 1 ( L M ) = A ( B ) for some function B ( t , x , v ) , gives the relation
η ^ 11 L , x M + η ^ 12 L , v M = A ( B ) A η ^ 11 L , v M B = 0 .
That is, we get the first integral
Ω ( ω 1 , ω 2 ) = η ^ 11 M d v + Q ( t , x ) B ( t , x , v ) = const .
This is the main idea of Noether’s theorem, which relates a variational symmetry with a constant of motion. Let us note that function B ( t , x , v ) is more general than gauge function G ( t , x ) , which can be added to any Lagrangian in the form of a total derivative with respect to t. If we consider a Lagrangian of the form L = L M + A ( G ( t , x ) ) , then
S ^ 1 ( L ) = S ^ 1 ( L M ) + S ^ 1 ( A ( G ) ) = S ^ 1 ( L M ) + η ^ 11 A ( G , x ) + A ( η ^ 11 ) G , x = S ^ 1 ( L M ) + A ( η ^ 11 G , x ) .
Comparing the above relation with the relaxed condition, we obtain
B ( t , x , v ) = η ^ 11 G , x ( t , x ) ,
which is a very particular form, and actually, G , x enters into the function Q in (27).
Let us consider as an example the classical harmonic oscillator x ¨ = k 2 x , which admits the scale symmetry S ^ 1 = x x + v v and the time-translation symmetry S ^ 2 = t A = v x + k 2 x v with Δ = k 2 x 2 + v 2 = M 1 . Introducing “polar” coordinates v = r sin ϕ , k x = r cos ϕ , we obtain A r ϕ = t k ϕ , S ^ 1 r ϕ = r r , S ^ 2 r ϕ = k ϕ , M = r 2 . Therefore,
L M = v k x arctan v k x 1 2 ln ( v 2 + k 2 x 2 ) = ϕ tan ϕ ln r .
This Lagrangian was mentioned in [26,27]. Condition S ^ 1 r ϕ ( L ) = 1 = A r ϕ ( B ) is fulfilled for B = t . From (27), we obtain the first integral
Ω 1 = x ϕ r cos ϕ + t = const . k t + ϕ = c 1 .
Taking into account that in “polar” coordinates, r ˙ = 0 r = c 2 = const , we have the well-known solution v = c 2 sin ( c 1 k t ) and k x = c 2 cos ( c 1 k t ) .

3.4. Systems with the Same Jacobi Multiplier

Let us now consider a modified Lagrangian given by
L ˜ : = L M g ( t , x ) d x ,
where g ( t , x ) is an arbitrary function. Note that L ˜ and L M have the same Jacobi multiplier M by construction. Using Euler–Lagrange Equation (26), we obtain
f ˜ = f g ( t , x ) M .
Thus, we have two operators A and A ˜ : = A Δ g ( t , x ) v , which represent two different systems with different Lagrangians but with the same Jacobi multiplier. A natural question arises: if Δ is formed by two nonproportional symmetries S ^ 1 , 2 of A, is it possible to get two symmetries S ˜ 1 , 2 of A ˜ as functions of η ^ i j ?
Let us determine at least when symmetry S ^ 1 of A is a symmetry for A ˜ , i.e., when S ˜ 1 = S ^ 1 . From the commutator [ S ^ 1 , A ˜ ] = 0 , we obtain η ^ , v 11 = 0 and η ^ , x 11 = S ^ 1 ln g Δ . Condition g , v = 0 leads to the following relation
g ( t , x ) = I ( t , ω ) Δ , v ,
where I ( t , ω ) is a general invariant of operator S ^ 1 , i.e., S ^ 1 ( I ) = 0 (see also (19)). The specific form of function I should be selected to avoid the dependence of g on v. If there exists a symmetry S ˜ 2 = η ˜ 21 x + η ˜ 22 v which commutes with S ^ 1 , then we obtain S ^ 1 ( η ˜ 21 ) = S ˜ 2 ( η 11 ) . Using the relation Δ = η ^ 11 η ˜ 22 η ˜ 21 η ^ 12 , we have S ^ 1 η ˜ 21 / η ^ 11 = 0 , i.e.,
η ˜ 21 = η ^ 11 I ( t , ω ) .
Thus, there is a symmetry S ˜ 2 if function I ( t , ω ) in (29) satisfies the following equation (similar to (23))
A ˜ ( I ) = I , t + I , ω ( v ω , x + f ˜ ω , v ) = Δ η ^ 11 2 .
Let us consider the same example with a harmonic oscillator (see previous subsection). The general invariant of S ^ 1 is I ( t , ω ) = I ( t , v / x ) . To determine the modified operator A ˜ which commutes with S ^ 1 and has the same Δ , we start from (28) and obtain the necessary condition for g
g ( t , x ) = I ( t , v / x ) 2 v .
Taking into account that g , v = 0 , we obtain g = τ ( t ) / x . Thus, any equation x ¨ = f ˜ = k 2 x ( 1 + τ ( t ) ) τ ( t ) v 2 / x admits the scale symmetry S ^ 1 .
To satisfy Equation (30), which takes the form
I , t ( 1 + τ ) I , ω ( k 2 + ω 2 ) = k 2 + ω 2 ,
we assume I , t = 0 , so τ = c = const . and I = ω / ( c + 1 ) . Thus, η ˜ 21 , which defines the second symmetry of the modified equation, satisfies the relation (29), i.e., η ˜ 21 = v / ( c + 1 ) . Finally, we obtain the equation of the “hidden” oscillator x ¨ = k 2 x ( 1 + c ) c x ˙ 2 / x , which admits the symmetry
S ˜ 2 = v c + 1 x + k 2 x + c c + 1 v 2 x v ,
which is nothing more than t = ( c + 1 ) S ˜ 2 + A ˜ .

4. Conclusions and Discussion

Firstly, we analyzed how the symmetries of two first-order ODEs with the same integrating factor were related. We found out, as a general result, that two such equations admitted the same point symmetry.
Secondly, we established the relation between the Jacobi multiplier for second-order ordinary differential equations and the integrating factors of two first-order ODEs, expressed in terms of the invariant of an admitted symmetry. We pointed out that the knowledge of one admitted symmetry and the Jacobi multiplier allowed us to reduce the order of the equation for the second admitted symmetry.
Moreover, we found out that there existed a set of essentially different Euler–Lagrange ODEs of the second order with different Lagrangians, but with the same Jacobi multiplier. For such systems, we obtained the conditions for which the admitted symmetries were the same. The general problem of relating the symmetries of such systems is still open. Of course, for every second-order ODE, there are infinitely many contact transformations transforming it into any other given second-order equation. However, how to determine such transformations? Probably, the indicated connection between the pseudo scalar product and a nonzero Jacobi multiplier and the use of canonical variables may shed some light on this issue, and it will be the subject of our further investigation.
The obtained results are important for exploring the quantization of Hamiltonian systems in which one can use the Lie symmetries to obtain the quantum-correct Schrödinger equation by preserving the Noether symmetries [28,29,30,31]. Moreover, one can use these results to introduce new nonlinear systems to explore their classical and quantum dynamics, as was recently done in Ref. [32].

Author Contributions

Conceptualization, formal analysis, investigation, G.G.C. and A.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Consejo Nacional de Ciencia y Tecnología: 1757, A1-S-43579.

Data Availability Statement

Not applicable.

Acknowledgments

A.Y. thanks the program PRO-SNI 2023 of University of Guadalajara (México). G.G.C. would like to acknowledge support from the program Cátedras CONACyT through project 1757 and from project A1-S-43579 of SEP-CONACyT Ciencia Básica and Laboratorio Nacional de Ciencia y Tecnología de Terahertz. G.G.C. is also grateful to Saiful I. Khondaker for his kind hospitality at the NanoScience Technology Center at the University of Central Florida.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lie, S. Vorlesungen Über Continuierliche Gruppen mit Geometrischen und Anderen Anwendungen; Druck and Verlag von B. G. Teubner: Leipzig, Germany, 1893. [Google Scholar]
  2. Lie, S. Classification und Integration von gewöhnlichen Differentialgleichungen zwischen xy, die eine Gruppe von Transformationen gestatten. Math. Ann. 1888, 32, 213–281. [Google Scholar] [CrossRef]
  3. Ovsyannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  4. Ibragimov, N.H. Transformation Groups Applied to Mathematical Physics; Springer: Dordrecht, The Netherlands, 1985. [Google Scholar]
  5. Stephany, H. Differential Equations: Their Solution Using Symmetries; Cambridge University Press: Cambridge, UK, 1989. [Google Scholar]
  6. Olver, P.J. Application of Lie Groups to Differential Equations, 2nd ed.; Springer: New York, NY, USA, 1993. [Google Scholar]
  7. Jacobi, C.G.J. Sul principio dell’ultimo moltiplicatore e suo uso come nuovo principio generale di meccanica. G. Arcadico Sci. Lett. Arti 1844, 99, 129–146. [Google Scholar]
  8. Lie, S. Verallgemeinerung und neue Verwertung der Jacobischen Multiplikatortheorie. Christiania Forh. 1874, XIV, 255–274, reprinted in Gesammelte Abh. 1973, III, 188–206. [Google Scholar]
  9. Bianchi, L. Lezioni sulla Teoria dei Gruppi Continui Finiti di Trasformazioni; Enrico Spoerri: Pisa, Italy, 1903. [Google Scholar]
  10. Nucci, M.C.; Leach, P.G.L. Jacobi’s last multiplier and the complete symmetry group of the Euler–Poinsot system. J. Nonlin. Math. Phys. 2002, 9, 110–121. [Google Scholar] [CrossRef] [Green Version]
  11. Nucci, M.C.; Leach, P.G.L. Jacobi’s last multiplier and symmetries for the Kepler problem plus a lineal story. J. Phys. A Math. Gen. 2004, 37, 7743–7753. [Google Scholar] [CrossRef]
  12. Nucci, M.C.; Leach, P.G. Lagrangians galore. J. Math. Phys. 2007, 48, 123510. [Google Scholar] [CrossRef] [Green Version]
  13. Nucci, M.; Leach, P. The Jacobi Last Multiplier and its applications in mechanics. Phys. Scr. 2008, 78, 065011. [Google Scholar]
  14. Nucci, M.; Leach, P. An old method of Jacobi to find Lagrangians. J. Nonlin. Math. Phys. 2009, 16, 431–441. [Google Scholar] [CrossRef] [Green Version]
  15. Nucci, M.; Tamizhmani, K. Lagrangians for dissipative nonlinear oscillators: The method of Jacobi last multiplier. J. Nonlin. Math. Phys. 2010, 17, 167–178. [Google Scholar] [CrossRef] [Green Version]
  16. Cariñena, J.; De Lucas, J.; Rañada, M.F. Jacobi multipliers, non-local symmetries, and nonlinear oscillators. J. Math. Phys. 2015, 56, 063505. [Google Scholar]
  17. Cariñena, J.F.; Fernández-Núñez, J. Jacobi Multipliers in integrability and the inverse problem of mechanics. Symmetry 2021, 13, 1413. [Google Scholar] [CrossRef]
  18. Vinogradov, A.M. Local symmetries and conservation laws. Acta Appl. Math. 1984, 2, 21–78. [Google Scholar] [CrossRef]
  19. Krasil’shchik, I.S.; Vinogradov, A.M. Symmetries and Conservation Laws for Differential Equations of Mathematical Physics; American Mathematical Society: Providence, RI, USA, 1999. [Google Scholar]
  20. Nucci, M.C. Jacobi last multiplier and Lie symmetries: A novel application of an old relationship. J. Nonlin. Math. Phys. 2005, 12, 284–304. [Google Scholar] [CrossRef] [Green Version]
  21. Hojman, S.A. A new conservation law constructed without using either Lagrangians or Hamiltonians. J. Phys. A Math. Gen. 1992, 25, L291–L295. [Google Scholar] [CrossRef]
  22. Ibragimov, N. Opit Gruppovogo Analiza Obiknovennih Differencialnih Uravneniy; Znanie: Moscow, Russia, 1991. (In Russian) [Google Scholar]
  23. González-López, A. Symmetry and integrability by quadratures of ordinary differential equations. Phys. Lett. A 1988, 133, 190–194. [Google Scholar] [CrossRef]
  24. Ibragimov, N.K. Group analysis of ordinary differential equations and the invariance principle in mathematical physics. Russ. Math. Surv. 1992, 47, 89–156. [Google Scholar] [CrossRef] [Green Version]
  25. Whittaker, E.T. A Treatise on the Analytical Dynamics of Particles and Rigid Bodies, 2nd ed.; Cambridge University Press: Cambridge, UK, 1917. [Google Scholar]
  26. Havas, P. The range of application of the Lagrange formalism. Nuovo C 1957, 5, 363–388. [Google Scholar] [CrossRef]
  27. Bersani, A.M.; Caressa, P. Lagrangian descriptions of dissipative systems: A review. Math. Mech. Solids 2021, 26, 785–803. [Google Scholar] [CrossRef]
  28. Gubbiotti, G.; Nucci, M.C. Quantization of quadratic Liénard-type equations by preserving Noether symmetries. J. Math. Anal. Appl. 2015, 422, 1235–1246. [Google Scholar] [CrossRef]
  29. Nucci, M.C. Quantizing preserving Noether symmetries. J. Nonlin. Math. Phys. 2013, 20.3, 451–463. [Google Scholar] [CrossRef] [Green Version]
  30. Nucci, M.C. Quantization of classical mechanics: Shall we Lie? Theor. Math. Phys. 2011, 168, 994–1001. [Google Scholar] [CrossRef]
  31. Nucci, M.C. From Lagrangian to quantum mechanics with symmetries. J. Phys. Conf. Ser. 2012, 2012 380, 012008. [Google Scholar] [CrossRef]
  32. Contreras Gonzalez, G. Classical and quantum dynamics of over-damped non linear systems. Rev. Mex. Física 2023, 69, 010701-1. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

González Contreras, G.; Yakhno, A. Symmetries of Systems with the Same Jacobi Multiplier. Symmetry 2023, 15, 1416. https://doi.org/10.3390/sym15071416

AMA Style

González Contreras G, Yakhno A. Symmetries of Systems with the Same Jacobi Multiplier. Symmetry. 2023; 15(7):1416. https://doi.org/10.3390/sym15071416

Chicago/Turabian Style

González Contreras, Gabriel, and Alexander Yakhno. 2023. "Symmetries of Systems with the Same Jacobi Multiplier" Symmetry 15, no. 7: 1416. https://doi.org/10.3390/sym15071416

APA Style

González Contreras, G., & Yakhno, A. (2023). Symmetries of Systems with the Same Jacobi Multiplier. Symmetry, 15(7), 1416. https://doi.org/10.3390/sym15071416

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