/ / Article

Research Article | Open Access

Volume 2020 |Article ID 3170130 | https://doi.org/10.1155/2020/3170130

Z. E. Musielak, N. Davachi, M. Rosario-Franco, "Lagrangians, Gauge Functions, and Lie Groups for Semigroup of Second-Order Differential Equations", Journal of Applied Mathematics, vol. 2020, Article ID 3170130, 11 pages, 2020. https://doi.org/10.1155/2020/3170130

# Lagrangians, Gauge Functions, and Lie Groups for Semigroup of Second-Order Differential Equations

Revised21 May 2020
Accepted08 Jun 2020
Published22 Jun 2020

#### Abstract

A set of linear second-order differential equations is converted into a semigroup, whose algebraic structure is used to generate novel equations. The Lagrangian formalism based on standard, null, and nonstandard Lagrangians is established for all members of the semigroup. For the null Lagrangians, their corresponding gauge functions are derived. The obtained Lagrangians are either new or generalization of those previously known. The previously developed Lie group approach to derive some equations of the semigroup is also described. It is shown that certain equations of the semigroup cannot be factorized, and therefore, their Lie groups cannot be determined. A possible solution of this problem is proposed, and the relationship between the Lagrangian formalism and the Lie group approach is discussed.

#### 1. Introduction

Let be a set of all linear second-order differential equations (ODEs) of the form , where , and and are smooth () and ordinary ( and ) functions to be determined. In general, contains a wide variety of ODEs that can be separated into two subsets, one with known solutions and the other with unknown solutions. According to Murphy [1], there are more than five hundred ODEs in whose solutions are known, and many of those solutions are given by the special functions (SF) of mathematical physics as defined in [1, 2].

An algebraic structure can be added to , and as a result, the set becomes a semigroup. Let be a semigroup and let “+” be the binary operation of addition, such that . Moreover, the associativity axiom requires that . With the above conditions being obeyed by all elements of , the semigroup also preserves the commutative rule , which means that is a commutative semigroup. It must be pointed out that the algebraic structure imposed on the commutative semigroup allows generating both new and known ODEs (see Section 2) that are also elements of .

Let us also point out that an identity element may also be formally introduced into the semigroup by taking and finding the solution with and being constants. Having defined this identity element, the semigroup becomes a monoid , and since the original is commutative, is also commutative. The monoid structure of the considered ODEs is an interesting property that may be studied further in future papers; however, in this paper, only the commutative semigroup of ODEs is investigated.

The first main objective of this paper is to establish the Lagrangian formalism for all ODEs of by using standard and nonstandard Lagrangians. The standard Lagrangians (SLs) are typically expressed as the difference between terms that can be identified as the kinetic and potential energy [3]. On the other hand, for nonstandard Lagrangians (NSLs), originally introduced by Arnold [4] who referred to them as nonnatural Lagrangians, identification of kinetic and potential energy terms may not be obvious. Among different applications, Alekseev and Arbuzov [5] used the NSLs to formulate the Yang-Mills field theory. Extensive discussions of methods to derive NSLs can be found in some previous [68] and more recent work [9, 10]. Other important methods to obtain NSLs were developed by El-Nabulsi [10], Nucci [1114], Cariñena et al. [1517], and Saha and Talukdar [18].

Lagrangians are not unique, which means that some extra terms may be added to them and they give the same original equation. These are null (NLs) or trivial Lagrangians (TLs) as they make the Euler-Lagrange (E-L) equation to vanish identically [19, 20]. In other words, the NLs can be added to any Lagrangian without changing the derivation of the original equation. It is also required that the NLs can be expressed as the total derivative of a scalar function, which is called a gauge function [19, 21, 22]. The null Lagrangians have been extensively studied in mathematics [2328] but only with a couple of applications in physics [29, 30].

In theoretical physics, the Lagrangian formalism is central to any classical or quantum theory of particles, waves, or fields, and the fundamental equations of modern physics are typically derived from given Lagrangians [3]. If the original equations are given but their Lagrangians are not known, then the inverse (or Helmholtz) problem of the calculus of variations must be solved [31]. In this paper, the inverse problem for all ODEs of is solved, and their SLs, NLs, and NSLs are derived; the existence of the obtained Lagrangians is verified by using the Helmholtz conditions [32], and it is demonstrated that the derived NSLs violate the Helmholtz conditions.

As the second main objective of this paper, the Lie group approach is applied to the ODEs of . A Lie group for a given ODE of must be known in advance, so its irreducible representations (irreps) can be determined and used to obtain the corresponding solutions of this ODE. Moreover, operators of the Lie algebra of the selected Lie group can be used to obtain the original ODE. This work has been done so far for the ODEs whose solutions are the SF of mathematical physics [3337] (called here the SF ODEs), and its limitations are explored in this paper. Now, if the original equation is given but its Lie group is not known, then the original equation must be factorized [38, 39], and the resulting lower-order ODE must be related to the Lie algebra, from which the corresponding Lie group is determined. In this paper, we establish the validity of the Lie group approach to ODEs of (beyond the SF ODEs) and demonstrate the existence of limits on the factorization for some ODEs of .

Finally, the Lagrangian formalism and the Lie group approach, the two independent methods of obtaining the ODEs of , are compared. The advantages and disadvantages of each method are discussed, and it is shown that the Lagrangian formalism can be established for all considered ODEs; however, the Lie group approach is only limited to some ODEs that form a subsemigroup of . Applications of the obtained results to selected problems in mathematical physics are also discussed.

The outline of the paper is as follows: a procedure of generating new ODEs in is described and applied to some ODEs in Section 2; the Lagrangian formulation for all ODEs of is established by using standard and nonstandard Lagrangians and applied to selected ODEs in Section 3; the Lie group approach and factorization methods are presented in Section 4; comparison between the two methods is discussed in Section 5; and conclusions are given in Section 6.

#### 2. Commutative Semigroup of Differential Equations

An interesting property of is that its algebraic operation can be used to generate novel ODEs by simply adding two different elements of the semigroup. Since any two elements (two chosen ODEs) of can be added up, the results may be rather surprising as shown by a couple of examples below. Let and , where , , and are smooth and arbitrary functions, be two operators, and let and be two ODEs of . Then, the binary operation of allows writing where and .

Here are some interesting examples. All Bessel (regular, modified, spherical, and modified spherical) equations are elements of . To add the regular and modified Bessel equations, we have , and , where is either real or integer constant, which gives

Since the eigenvalue of this equation is (see Section 4), the resulting ODE is neither regular nor modified Bessel equation, but instead, it is known as the Euler equation [1]. Addition of other Bessel equations gives novel Bessel and Euler equations.

To show these novel equations, let a general form of Bessel equation be where and . In addition, and are constants, whose values correspond to four Bessel equations and the Euler equation (see Table 1). In addition, is either integer or real. Note that the form of Eq. (3) considered here is more general than that used in [10] as it allows reducing the general Bessel equation to the Euler equations (see Table 1).

 Known and derived equations Regular Bessel 1 1 Modified Bessel 1 -1 Spherical Bessel 2 1 Modified spherical Bessel 2 -1 Semispherical Bessel 3/2 1 Modified semispherical Bessel 3/2 -1 Regular Euler 1 0 Spherical Euler 2 0 Semispherical Euler 3/2 0
Novel equations derived by using the binary operation of S.

Note that for , a new name “semispherical” was introduced for both Bessel and Euler equations. Additions of Bessel equations required to obtain the Euler equation and the novel equations (see Table 1) are presented in Table 2. It must be pointed out that further additions of the Bessel equations and the derived equations are also possible but not considered here.

 Derived equations Addition of Bessel equations Semispherical Bessel Regular and spherical Modified semispherical Bessel Modified and modified spherical Regular Euler Regular and modified Spherical Euler Spherical and modified spherical Semispherical Euler Regular and modified spherical Semispherical Euler Modified and spherical
Novel equations derived by using the binary operation of S.

Another interesting result is obtained when the regular Legendre equation with and is added up together with the associated Legendre equation with and , where and are constants. The resulting ODE is again the associated Legendre equation given by where . This shows that only the constant is affected by the addition, but the form of the resulting ODE does not change.

The considered simple examples show that numerous novel ODEs can be generated by using the algebraic operation of the semigroup . Since the process of adding the elements of to each other can be carried on ad infinitum, there are infinite different resulting ODEs, and all of them are the elements of .

#### 3. Lagrangian Formalism and Its Lagrangians

Let be a functional that depends on an ordinary and smooth function , and let be defined by an integral over a smooth function , which is called Lagrangian, and depends on , and on , so . The principle of least action or Hamilton’s principle [3] requires that , where is the variation defined as the Fréchet derivative of with respect to . Using the condition , the E-L equations are obtained and their solutions give that makes the action stationary. The described procedure is the basis of the classical calculus of variations, and it works well when the Lagrangian is either already given or it can be determined for a given physical system.

The basis for the Lagrange formalism is the jet-bundle theory [25, 26]. Let and be differentiable manifolds of dimensions and , respectively, and let be a fibred bundle structure, with being the canonical projection of the fibration. Let be the th jet bundle, with , , and . Then, an ODE of order is called locally variational (or the E-L type) if, and only if, there exists a local real function constrainted by the condition . For the ODEs of , and is a local Lagrangian. Such local Lagrangians are not unique as other Lagrangians may also exist and they would give the same original equations when substituted into the E-L equations. Now, if Lagrangians are not known, then the problem of finding them is the inverse (or Helmholtz) variational problem [31, 32, 40]. One of the main goals of this paper is to solve this problem for all ODEs of .

##### 3.1. Standard Lagrangians

In the previous work [6, 7], special forms of standard Lagrangians were derived for . These results are generalized in Proposition 1 by allowing the Lagrangians to depend on smooth and ordinary functions that are to be determined.

Proposition 1. Let be the standard Lagrangian given by and let , and be smooth and ordinary functions to be determined. The Lagrangian can be used to obtain the original equation , if and only if, , and .

Proof. Substitution of into the E-L equation gives or if . This E-L equation can be converted into the original equation if and , which concludes the proof.
Other implications of the results of Proposition 1 are now presented by the following corollaries.

Corollary 2. The function is given as , where is an integration constant and .

Corollary 3. The Euler-Lagrange equations do not supply enough constraints to determine both and , thus, one of these functions must be specified. Different choices of or lead to different standard Lagrangians.

For a given physical problem with specified intitial conditions, more constraints on the functions and can be imposed. However, in this paper, we keep our analysis general and do not consider any specific physical problem with given initial conditions.

We may use the result of Corollary 3 to define three Lagrangians that play special roles in this paper. In general, many different Lagrangians can be obtained but considering all of them is out of the scope of this paper. Let us begin with the choice of , which sets . The result is a Lagrangian that is called here the minimal Lagrangian , and it is defined as

It is easy to verify that this Lagrangian is the simplest standard Lagrangian that allows obtaining the original equation directly from the E-L equations (6). Let us also point out that the same Lagrangian is obtained when , which requires that depends on and that their relationship is , where is another integration constant.

We may also choose = const, which gives and the so-called middle Lagrangian is obtained where is a null Lagrangian [10]. The middle Lagrangian is a new standard Lagrangian.

Finally, we take , which gives , and the maximal Lagrangian can be written as where is another null Lagrangian. The Lagrangian was also previously found [7], and it was demonstrated that this Lagrangian gives the original equation upon substitution into the E-L equation. It is easy to verify that even if the obtained Lagrangians are different they lead to the same original equation .

##### 3.2. Null Lagrangians

A well-known result is that the difference between two Lagrangians that give the same equation can be written as the total derivative of scalar function [20]. This means that we may write where the functions , , and exist on the configuration manifold (see Section 2); however, it may not always be possible to find a single-value function on the entire [20]. Since the null Lagrangians can be added to any Lagrangian without modifying the resulting original equations, we may call them “gauge transformations” in reference to the gauge transformation of electromagnetism. Moreover, the above results demonstrate that the functions , , and are the gauge functions and that they exist [20].

It is easy to find , where is an integration constant; however, the other two functions require to be specified, which will not be done here to keep our approach as general as possible. Important consequences of the fact that the null Lagrangians allow finding the gauge functions are now summarized in the following corollaries.

Corollary 4. The null Lagrangians and are the solutions that satisfy the E-L equations but cannot be used to derive

Corollary 5. The null Lagrangians and can be added to any known Lagrangian without making any changes in the resulting original equation.

Finally, let us point out that both and depend on and but are independent of . Since the first-derivative term with its coefficient can be removed from by using the integral transformation given by Proposition 10 (see Section 3.2), the null Lagrangians for such transformed ODEs would be zero. Our results show that with and , the minimal Lagrangian is the simplest and most fundamental standard Lagrangian for all ODEs of .

##### 3.3. Nonstandard Lagrangians

We consider nonstandard Lagrangians to be given in the following general form [9, 10]. where , and are smooth functions to be determined. Comparing to given by Eq. (5), it is seen that both Lagrangians have three arbitrary functions; however, the forms of these Lagrangians are significantly different.

Studies of the Lagrangian previously performed [8] were limited to only a few special cases that allowed uniquely determine the functions , , and . These functions were determined by substituting into the E-L equations. First, it was shown that separates from the other two functions, and as a result, it can be evaluated independently. Second, the evaluation of and is connected, and they can only be obtained if a Riccati equation is solved but neither [7] nor [8] gives the required solutions. Therefore, here, we solve the resulting Riccati equation (9) in the following proposition.

Proposition 6. Let be the Riccati equation, and let , where is an arbitrary function of the nonstandard Lagrangian given by Eq. (12). The solution of this Riccati equation is if, and only if, the function satisfies .

Proof. We transform the Riccati equation by introducing a new variable , which is related to by with , and obtain where .
Now, we transform Eq. (15) by using which gives [10], and the solution is where must satisfy . This concludes the proof.

Having obtained the solution of the Riccati equation, we may use it to calculate the functions and and find the following nonstandard Lagrangian where . It is seen that the derived Lagrangian depends explicitly on , which is a solution to that becomes an auxiliary condition for the developed Lagrangian formalism based on non-standard Lagrangians.

The obtained nonstandard Lagrangians were previously derived but only for some ODEs whose solutions are special functions of mathematical physics [6, 8, 10]. Those previous results are here generalized to all ODEs of . Both standard and nonstandard Lagrangians allow obtaining the original equations; however, the nonstandard Lagrangians lead to a new phenomenon in calculus of variations, which is the requirement of the auxiliary condition in order to obtain

##### 3.4. Applications to Bessel and Euler Equations

We now derive Lagrangians for the general Bessel equation given by Eq. (3). More specifically, we find the minimal, middle, and maximal standard Lagrangians using Eqs. (8), (9), and (10), respectively, as well as the nonstandard Lagrangian given by Eq. (18).

The minimal (basic) standard Lagrangian is given by where in Eq. (8).

The middle Lagrangian can be written as where with and being a null Lagrangian that does not depend on either or .

The maximal Lagrangian is where is another null Lagrangian that depends on but is independent of , which is only present in the basic Lagrangian.

According to Eq. (18), the nonstandard Lagrangian becomes with the auxiliary condition . It is seen that depends on but does not depend explicitly on ; the -dependence is only through the solution for .

The obtained Lagrangians are expressed in terms of two parameters and that uniquely determine whether the considered equation is the Bessel or Euler equation (see Table 1). In other words, the standard and nonstandard Lagrangians for each equation can be obtained from Eqs. (19) through (24) by using the values of and given in Table 1.

##### 3.5. Helmholtz Conditions

Our results show that the Lagrangian formalism may be established for all ODEs of by using either the standard or nonstandard Lagrangians. The existence of the obtained Lagrangians must be validated by the Helmholtz conditions [32]. There are three original Helmholtz conditions and two of them are trivially satisfied. The third condition requires a lot of caution, because its original version shows that no Lagrangian can be constructed for the general form of the ODEs used in this paper, namely, . It is easy to demonstrate that is inconsistent with the third Helmholtz condition [19].

To fix the problem of validity of the third Helmholtz condition, we substitute either the minimal, middle, or maximal Lagrangian into the E-L equation, and as expected, all three Lagrangians give the same original equation. The resulting ODE is , and this equation is consistent with the third Helmholtz condition. The problem arises when one more step is performed, namely, dividing the equation by , which gives the original ODE, but the resulting ODE violates the third Helmholtz condition; note that in principle, the division can be performed because in general . In other words, the third Helmholtz condition tells us that Lagrangian exists for , but it does not exist for ; however, it is the latter that we want to derive as it is our original equation. Thus, the origin of the inconsistency is clearly demonstrated.

The problem with the nonstandard Lagrangians is more profound because they do allow obtaining the original equation, but the resulting ODE does not obey the third Helmholtz condition. There are a few special cases that may lead to a result that is consistent with this condition. However, for the original equation , which is derived from the nonstandard Lagrangians, the results remain inconsistent with the third Helmholtz condition. This violation of the third Helmholtz condition is a known phenomenon in the calculus of variations [10].

#### 4. Lie Group Approach

##### 4.1. Special Function Equations

Linear second-order ODEs whose solutions are special functions (SFs) of mathematical physics are called here the special function equations (SFEs). All SFEs are elements of and they form a subsemigroup denoted , with . There are many applications of SFEs and SFs in applied mathematics, physics, and engineering [2, 34]. The Lie group approach provides a unifying perspective of all SFs and their SFEs. Therefore, it is surprising that the approach is not commonly known and that its descriptions can only be found in a very few advanced textbooks of mathematical methods (e.g., [34]) or monographs [35, 36] and reviews [37] devoted to this topic.

Using the previously obtained results [2124], we find that Bessel and spherical Bessel functions are obtained, respectively, from the Euclidean groups and , also known as the ISO(2) and ISO(3) groups. In addition, it was shown that Legendre polynomials and functions are derived, respectively, from unitary groups and unimodular quasi-unitary groups . As long as Hermite polynomials and functions are concerned, they can be obtained from Heisenberg and groups, respectively.

Let us now demonstrate how the procedure of finding special functions works by selecting the group, which consists a subgroup of all translations in , and a subgroup of all rotations in ; thus, the group structure is , where is an invariant subgroup. The procedure requires finding a function that transform as one of irreps of a given group; in the previous work [41], such function was found for the Galilean group of the metric; however, here, it is obtained for the Euclidean group in the following proposition and corollary.

Proposition 7. Let be a translation operator of and let be a smooth and ordinary function. The function transforms as the irreps of if, and only if, satisfies the eigenvalue equations where labels the irreps.

Proof. The action of on is given by where is a translation and . Making the Taylor series expansion of , we find Comparing Eq. (27) to Eq. (26), we obtain the following eigenvalue equations which represent the necessary conditions that transforms as the irreps of [40]. This concludes the proof.

Corollary 8. The solutions to the eigenvalue equations of Proposition 9 are the plane waves , where is an integration constant.

By performing the Fourier-sine and Fourier-cosine expansions of the plane waves of Corollary 8, one finds that each expansion has only two nonzero coefficients, and that one of these coefficients is proportional to , and the other to , with being the label of the coefficients and the resulting SFs being Bessel functions [2, 34]. Moreover, the Lie algebra of the group allows finding the raising and lowering operators. Using these operators, the power series expansions for and , their recurrence relationships, generating functions, and other properties may be obtained, without making any reference to the Bessel equation [3537]. To derive the Bessel equation, may be used; however, it must be noted that the resulting ODE is the Bessel equation written in its second canonical form (see Proposition 9).

As described above, the remarkable relationships between Lie groups and the SFs of the ODEs of have been found. The relationships are important as they unified the well-known and commonly used SFs with simple Lie groups used in physics and applied mathematics. To the best of our knowledge, the Lie group approach has only been applied to the ODEs of , and Lie groups for these equations have been identified. However, the semigroup contains also many other ODEs, which we denote here as , with the subscript standing for the remaining ODEs. In general, the solutions of are either known, but not given in terms of the SFs, or unknown. In the following, we determine whether the Lie group approach can also be applied to all ODEs of , or not.

##### 4.2. Remaining Equations

Now, if a Lie group for a given ODE is not known, then one may factorize the ODE, obtain the eigenvalue equations and the raising and lowering operators [20] that are related to a Lie algebra of infinitesimal group generators from which the corresponding Lie group can be determined [3537]. Using the obtained Lie algebra, the SFs naturally appear as the basis functions of differential operator representations of the algebra. Moreover, once the Lie group is known, its irreps can be found, and if is the matrix of such irreps, then their matrix element , where is a group element, becomes proportional to one of the SFs, or to a product of SFs. This shows that the Lie group approach is established only for those ODEs that are factorized. There are many factorization methods (e.g., [38] and references therein); however, only the Infeld and Hull (IH) factorization [39] shows clear relationships to Lie algebra [3537]; therefore, only the IH method will be used in this paper.

We now want to determine whether there are any ODEs of that cannot be factorized by using the IH method. Our main results are presented in the following proposition.

Proposition 9. Not all ODEs of can be factorized by using the Infeld and Hull factorization method [39].

Proof. The semigroup contains the ODEs of the form . However, in general, the functions , , and may depend on a parameter that can be either integer or real; thus, we may write the above ODEs as . Let be an eigenvalue of the operator , so that the corresponding eigenvalue equation becomes , or in its explicit form is Following the IH factorization method [39], Eq. (29) is converted into a Sturm-Liouville equation and then cast into its first canonical form, which can be written as where the Sturm-Liouville functions and , with , and with and being constants given by the boundary conditions.
The procedure of converting the first canonical form into the second canonical form involves two transformations [39]. Then, the second canonical form of Eq. (30) becomes where is the transformed dependent variable , is given in terms of , , and (see [42] for details). Expressing and by and in terms of and , we get Having obtained , it is shown that Eq. (31) can be factorized if it can be replaced by each of the following equations [42] and where is a function to be determined from [42], and are the raising, , and lowering, , operators that give and if is a solution of Eq. (31).
According to IH [39], there are six types of factorization denoted as , , , , , and . For each type, the functions , , and are given in Table 17 of [39], which shows that the -dependence of the functions and is different for each type of factorization; nevertheless, it remains fixed within the same type. Similarly for whose -dependence is also fixed. However, all three functions depend on some parameters that may change their values from one ODE to another. There are examples of ODEs in Table 17 of [39] that correspond to many known ODEs of and also other examples of ODEs that we would classify as elements of .
Now, using Eq. (32), the function can be obtained for any ODEs of . However, since and are arbitrary functions of , the function may have in principle any form, not necessary one of the forms of established by IH for their factorization types [39]; note that for each factorization type has its -dependence fixed. This shows that there are some ODEs of that cannot be factorized by the IH method. This concludes the proof.

The detailed calculations of casting the ODEs in their first and second canonical forms are lengthly and involve transformations that lead to a complicated form of [42]. We show in the following proposition how these calculations can be significantly simplified.

Proposition 10. Let , and let be an eigenvalue of the operator , so that the corresponding eigenvalue equation is

The following integral transformation converts Eq. (38) directly into its second canonical form (see Eq. (31)).

Proof. Applying the transformation given by Eq. (39) to Eq. (38), we obtain which is the second canonical form of Eq. (38), with Note that the obtained result neither requires casting Eq. (38) into its first canonical form nor using the original IH transformations [39]. This concludes the proof.

The results of Proposition 10 can be used to formulate the following corollaries.

Corollary 11. All ODEs of can be cast into their second canonical form.

Corollary 12. If the ODEs of are cast into their second canonical form, then all their null Lagrangians are zero.

It must also be noted that the null Lagrangians are zero for some ODEs of , such as Airy, Hill, Mathieu, Weber, Wittaker, and other similar ODEs.

Let us apply the results of Proposition 10 to the Bessel equation by taking , , , , and either or as the linearly independent solutions. By performing the transformation given by Eq. (39), the Bessel equation is cast in its second canonical form, and we find , which gives , in agreement with the results given in Section 3.1.

We also want to point out that the IH factorization of the considered ODEs and the Lie group approach are closely related as shown by the following corollary.

Corollary 13. The Lie group approach can only be established for those ODEs that can be factorized using the IH factorization method [39]. The results of Proposition 9 show that there are some ODEs of that cannot be factorized, which means that the Lie group approach cannot be established for these ODEs.

Our description of the Lie group approach showed that only ODEs of can be factorized, and therefore, the Lie group approach can be established for them. We also established that all ODEs of can be converted into their second canonical forms, and that the general form of for these ODEs can be obtained. Despite the fact that can be derived, our results showed that there are some ODEs of that cannot be factorized using the IH method [39].

#### 5. Lagrange Formalism versus Lie Group Approach

In this paper, we have established the Lagrangian formalism based on standard and nonstandard Lagrangians for all ODEs of and demonstrated that the previously formulated Lie group approach is applicable to only some ODEs of . Let be a subsemigroup of all ODEs for which both the Lagrange formalism and the Lie group approach are established, and let be another subsemigroup of all remaining ODEs for which only the Lagrange formalism is established, with . The ODEs of seem to be special as they all belong to , so it would be interesting to find deeper connections between the Lagrange formalism and the Lie group approach for these equations. Similarly, for other ODEs of , which are also elements of . Now, for the ODEs of , the underlying Lie groups are not yet known, so it would be interesting to find them.

Possible solutions to the above problems involve developing the Lagrange formalism jointly with the Lie group approach, which would require replacing the configuration manifold by a manifold associated with a given Lie group, so that . The approach guarantees that the resulting Lagrangian is -invariant; however, it may require deriving “new” variational principles and “new” Euler-Lagrange equations [43, 44]. The resulting invariance of is its important property that is strongly related to the Lie group . Actually, the problem may be reversed, and for every known Lagrangian, its invariance with respect to rotations, translations, and boots may indicate the presence of the underlying Lie groups [45, 46]. In other words, the groups may be identified by investigating the Lagrangian invariance. Moreover, the invariance of guarantees that the original equation derived from this Lagrangian preserved the same invariance [47]. It may be also interesting to explore the relevance of the recently discovered particle-like structure of Lie algebra [48] to both the Lagrange formalism and the Lie group approach.

Let us now point out that the null Lagrangians that were identified in this paper as gauge transformations allow finding the gauge functions, and also likely gauge groups underlying these transformations. The latter is not included in our paper as it will be a subject of future explorations. Specifically, it would be of great interest to discovered relationships between the gauge groups and the Lie groups, whose irreps may be used to obtain the ODEs of mathematical physics with the special function solutions.

The above topics are of great importance in mathematical physics, and their exploration may give new connections between the Lagrange formalism and the gauge groups, and the Lie groups for the ODEs of . It may also help finding Lie groups that underlie the ODEs of . However, these topics are out of the scope of the present paper, but they will be investigated in future papers.

#### 6. Conclusions

A set of general second-order linear ordinary differential equations with nonconstant coefficients was considered, and an algebraic structure (binary addition) was added to this set to form a semigroup. An interesting result is that the operation of the semigroup can be used to produce many new differential equations, which are still elements of the semigroup. Among a few presented examples, novel Bessel equations with their zero eigenvalues were obtained and discussed.

The Lagrangian formalism was established for all ODEs of the semigroup. We solved the inverse variational problem for these equations and derived standard and nonstandard Lagrangians. Among the derived standard Lagrangians, we obtained the minimal, middle, and maximal Lagrangians and demonstrated that they are equivalent Lagrangians. Moreover, we showed that the minimal Lagrangian is the simplest and also the most fundamental, and that the other Lagrangians differ by the so-called null Lagrangians from the minimal one. In our analysis of the null Lagrangians, we found out that these Lagrangians correspond to addition of the total derivative of a scalar function, which is a well-known result in calculus of variations [3, 20], and in this paper, it was identified as gauge transformations. We also used the null Lagrangians to determine the gauge functions for all considered ODEs. As long as nonstandard Lagrangians are concerned, our results showed that these Lagrangians require that calculus of variations is amended by auxiliary conditions that are different for different ODEs. This is a novel phenomenon in calculus of variations that is responsible for violation of the third Helmholtz conditions for these Lagrangians.

A Lie group approach was briefly described and compared to the Lagrange formalism established in this paper. The comparison shows that the Lagrange formalism is more robust as it can be applied to all considered equations, and that the Lie group approach is limited to the equations whose solutions are the special functions of mathematical physics and some other similar equations. We identified a subsemigroup of ODEs that cannot be obtained by the Lie group approach. Moreover, we also suggested alternative approaches to establish deeper connections between the Lagrangian formalism based on the standard and nonstandard Lagrangians, the gauge (still to be determined), and the Lie group approach.

#### Data Availability

All results presented in our paper are easily reproducible by any mathematician or physicist. All what it takes is to follow the steps outlined in the paper, and it is straightforward to derive the presented results. There is no other data in the paper required to be shared.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

We thank two anonymous referees for their valuable comments and suggestions that allow us to significantly improve the original version of this paper. This work was supported by the Alexander von Humboldt Foundation (Z.E.M.) and by the LSAMP Program at the University of Texas at Arlington (M.R.).

#### References

1. G. M. Murphy, Ordinary Differential Equations and Their Solutions, Dover Publication, Inc., New York, 2011.
2. A. M. Mathai and H. J. Haubold, Special Functions for Applied Scientists, Springer, New York, 2008. View at: Publisher Site
3. N. A. Doughty, Lagrangian Interaction, Addison-Wesley, New York, 1990.
4. V. I. Arnold, Mathematical Methods of Classical Mechanics, Springer, New York, 1978. View at: Publisher Site
5. A. I. Alekseev and B. A. Arbuzov, “Classical Yang-Mills field theory with nonstandard Lagrangians,” Theoretical and Mathematical Physics, vol. 59, no. 1, pp. 372–378, 1984. View at: Publisher Site | Google Scholar
6. J. L. Cieśliński and T. Nikiciuk, “A direct approach to the construction of standard and non-standard Lagrangians for dissipative-like dynamical systems with variable coefficients,” Journal of Physics A: Mathematical and Theoretical, vol. 43, no. 17, article 175205, 2010. View at: Publisher Site | Google Scholar
7. Z. E. Musielak, “Standard and non-standard Lagrangians for dissipative dynamical systems with variable coefficients,” Journal of Physics A: Mathematical and Theoretical, vol. 41, no. 5, article 055205, 2008. View at: Publisher Site | Google Scholar
8. Z. E. Musielak, “General conditions for the existence of non-standard Lagrangians for dissipative dynamical systems,” Chaos, Solitons & Fractals, vol. 42, no. 5, pp. 2645–2652, 2009. View at: Publisher Site | Google Scholar
9. N. Davachi and Z. E. Musielak, “Generalized Non-Standard Lagrangians,” Journal of Undergraduate Reports in Physics, vol. 29, no. 1, article 100004, 2019. View at: Publisher Site | Google Scholar
10. Z. E. Musielak, N. Davachi, and M. Rosario-Franco, “Special Functions of Mathematical Physics: A Unified Lagrangian Formalism,” Mathematics, vol. 8, no. 3, p. 379, 2020. View at: Publisher Site | Google Scholar
11. R. A. El-Nabusi, “Fractional field theories from multidimensional fractional variational problems,” International Journal of Geometric Methods in Modern Physics, vol. 5, p. 2614, 2008. View at: Google Scholar
12. R. A. El-Nabulsi, “Fractional Dirac operators and deformed field theory on Clifford algebra,” Chaos, Solitons & Fractals, vol. 42, no. 5, pp. 2614–2622, 2009. View at: Publisher Site | Google Scholar
13. R. A. El-Nabulsi, “Fractional action-like variational approach, perturbed Einstein’s gravity and new cosmology,” Fizika B-Journal of Experimental and Theoretical Physics-Zagreb, vol. 19, pp. 103–210, 2010. View at: Google Scholar
14. R. A. El-Nabulsi, “Gravitons in Fractional Action Cosmology,” International Journal of Theoretical Physics, vol. 51, no. 12, pp. 3978–3992, 2012. View at: Publisher Site | Google Scholar
15. M. C. Nucci and P. G. L. Leach, “Lagrangians galore,” Journal of Mathematical Physics, vol. 48, no. 12, article 123510, 2007. View at: Publisher Site | Google Scholar
16. M. C. Nucci and P. G. L. Leach, “Jacobi’s last multiplier and Lagrangians for multidimensional systems,” Journal of Mathematical Physics, vol. 49, no. 7, article 073517, 2008. View at: Publisher Site | Google Scholar
17. M. C. Nucci and P. G. L. Leach, “The Jacobi Last Multiplier and its applications in mechanics,” Physica Scripta, vol. 78, no. 6, article 065011, 2008. View at: Publisher Site | Google Scholar
18. J. F. Carinena, M. F. Ranada, and M. Santander, “Lagrangian formalism for nonlinear second-order Riccati systems: One-dimensional integrability and two-dimensional superintegrability,” Journal of Mathematical Physics, vol. 46, no. 6, article 062703, 2005. View at: Publisher Site | Google Scholar
19. A. Saha and B. Talukdar, “Inverse Variational Problem for Nonstandard Lagrangians,” Reports on Mathematical Physics, vol. 73, no. 3, pp. 299–309, 2014. View at: Publisher Site | Google Scholar
20. P. J. Olver, Applications of Lie Groups to Differential Equations, Springer-Verlag, New York, 1993. View at: Publisher Site
21. M. Crampin and D. J. Saunders, “On null Lagrangians,” Differential Geometry and its Applications, vol. 22, no. 2, pp. 131–146, 2005. View at: Publisher Site | Google Scholar
22. J.-M. Levy-Leblond, “Group-theoretical foundations of classical mechanics: The Lagrangian gauge problem,” Communications in Mathematical Physics, vol. 12, no. 1, pp. 64–79, 1969. View at: Publisher Site | Google Scholar
23. N. Wilson and A. H. Kara, “Equivalent Lagrangians: Generalization, Transformation Maps, and Applications,” Journal of Applied Mathematics, vol. 2012, Article ID 860482, 19 pages, 2012. View at: Publisher Site | Google Scholar
24. P. J. Olver and J. Sivaloganathan, Nonlinearity, vol. 1, p. 389, 1989.
25. D. R. Grigore, “Variational equations and symmetries in the Lagrangian formalism,” Journal of Physics A: Mathematical and General, vol. 28, no. 10, pp. 2921–2937, 1995. View at: Publisher Site | Google Scholar
26. D. Krupka and J. Musilová, “Trivial lagrangians in field theory,” Differential Geometry and its Applications, vol. 9, no. 3, pp. 293–305, 1998. View at: Publisher Site | Google Scholar
27. D. Krupka, O. Krupková, and D. Saunders, “The cartan form and its generalizations in the calculus of variations,” International Journal of Geometric Methods in Modern Physics, vol. 7, no. 4, pp. 631–654, 2011. View at: Publisher Site | Google Scholar
28. R. Vitolo, “On different geometric formulations of Lagrangian formalism,” Differential Geometry and its Applications, vol. 10, no. 3, pp. 225–255, 1999. View at: Publisher Site | Google Scholar
29. D. R. Anderson, D. E. Carlson, and E. Fried, “A Continuum-Mechanical Theory for Nematic Elastomers,” Journal of Elasticity, vol. 56, no. 1, pp. 33–58, 1999. View at: Publisher Site | Google Scholar
30. Z. E. Musielak and T. B. Watson, “Gauge functions and Galilean invariance of Lagrangians,” Physics Letters A, vol. 384, no. 26, article 126642, 2020. View at: Publisher Site | Google Scholar
31. J. Lopuszanski, The Inverse Variational Problem in Classical Mechanics, World Scientific, 1999. View at: Publisher Site
32. H. von Helmholtz, “Ueber die physikalische Bedeutung des Prinicips der kleinsten Wirkung. (Fortsetzung),” Journal für die reine und angewandte Mathematik (Crelles Journal), vol. 1887, no. 100, pp. 213–222, 1887. View at: Publisher Site | Google Scholar
33. J. Douglas, “Solution of the inverse problem of the calculus of variations,” Transactions of the American Mathematical Society, vol. 50, no. 1, p. 71, 1941. View at: Publisher Site | Google Scholar
34. C. D. Cantrell, Modern Mathematical Methods for Physicists and Engineers, Cambridge University Press, Cambridge, 2019.
35. W. Miller, Lie Theory and Special Functions, Academic Press Inc., New York, 1968.
36. N. Vilenkin, Special Functions and the Theory of Group Representations, American Mathematical Society, Providence, RI, 1968. View at: Publisher Site
37. E. G. Kalnins, Revista Colombiana de Matematicas, vol. 21, p. 1, 1997.
38. F. Schwarz, “Decomposition of ordinary differential equations,” Bulletin of Mathematical Sciences, vol. 7, no. 3, pp. 575–613, 2017. View at: Publisher Site | Google Scholar
39. L. Infeld and T. E. Hull, “The Factorization Method,” Reviews of Modern Physics, vol. 23, no. 1, pp. 21–68, 1951. View at: Publisher Site | Google Scholar
40. P. S. Bauer, “Dissipative Dynamical Systems: I,” Proceedings of the National Academy of Sciences of the United States of America, vol. 17, no. 5, pp. 311–314, 1931. View at: Publisher Site | Google Scholar
41. Z. E. Musielak and J. L. Fry, “Physical theories in Galilean space-time and the origin of Schrödinger-like equations,” Annals of Physics, vol. 324, no. 2, pp. 296–308, 2009. View at: Publisher Site | Google Scholar
42. Ş. S. Bayin, Mathematical Methods in Science and Engineering: Bayin/Mathematical, John Wiley & Sons, Hoboken, NJ, 2006. View at: Publisher Site
43. A. Lucas, “Lagrangian mechanics on Lie groups: a pedagogical approach,” November 2011, https://arxiv.org/abs/1111.1275v1. View at: Google Scholar
44. J. E. Marsden, S. Pekarsky, and S. Shkoller, “Symmetry reduction of discrete Lagrangian mechanics on Lie groups,” April 2000, https://arxiv.org/abs/math/0004018v1. View at: Google Scholar
45. Z. Muzsnay and G. Thompson, “Inverse problem of the calculus of variations on Lie groups,” Differential Geometry and its Applications, vol. 23, no. 3, pp. 257–281, 2005. View at: Publisher Site | Google Scholar
46. C. Muriel, J. L. Romero, and P. J. Olver, “Variational -symmetries and Euler-Lagrange equations,” Journal of Differential Equations, vol. 222, no. 1, pp. 164–184, 2006. View at: Publisher Site | Google Scholar
47. G. F. Torres del Castillo and A. Moreno-Ruiz, “Symmetries of the equations of motion that are not shared by the Lagrangian,” May 2017, https://arxiv.org/abs/1705.08446v1. View at: Google Scholar
48. A. M. Vinogradov, “Particle-like structure of Lie algebras,” Journal of Mathematical Physics, vol. 58, no. 7, article 071703, 2017. View at: Publisher Site | Google Scholar

Copyright © 2020 Z. E. Musielak et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.