Preprint
Article

A Theorem on Separated Transformations of Basis Vectors of Polynomial Space and Its Applications in Special Polynomials and Related Lie Algebra

Altmetrics

Downloads

173

Views

94

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

25 February 2023

Posted:

27 February 2023

You are already at the latest version

Alerts
Abstract
The present paper introduces a method of basis transformation of vector fields that is specifically applicable to polynomials space and differential equations with certain polynomials solutions such as Hermite, Laguerre and Legendre polynomials. The method based on separated transformation of vector space basis by a set of operators that are equivalent to the formal basis transformation and connected to it by linear combination with projection operators. Applying the Forbenius covariants yields a general method that incorporates the Rodrigues formula as a special case in polynomial space. Using the Lie algebra modules, specifically , on polynomial space results in isomorphic algebras whose Cartan sub-algebras are Hermite, Laguerre and Legendre differential operators. Commutation relations of these algebras and Baker-Campbell-Hausdorff formula gives new formulas for special polynomials
Keywords: 
Subject: Computer Science and Mathematics  -   Mathematics

1. Introduction

In mathematical physics and specifically quantum mechanics, the solution of many problems requires solving the differential equations and their eigenvalues problem . Hermite, Laguerre and Legendre polynomials and related differential equations are among the most applicable eigenvalues problem in physics and mathematics [1,2]. Schrodinger equation for hydrogen atom reduces to Legendre differential equation and quantum harmonic oscillator requires Hermite polynomials and related differential equation. The well-known Rodrigues formula yields the solutions (eigenfunctions) of many of these differential equations [3]. In this paper we interpret the Rodrigues formula as the transformation of some specific basis in polynomial space to another polynomials (eigenfunctions) of associated differential equations. This approach is feasible, provided that the transformation operator to be considered as a set of operators acting on each basis separately. It is shown that the overall action of these operators is equivalent to a single linear operator. As an example, the change of basis vectors in two dimension can be made by a matrix of rank 2. The action of this matrix could be equivalent with the actions of two different matrices that acts on each basis separately. The relation between these operators achieved by applying projection operators as is proved in the section 2. The connection between separated basis transformation and umbral composition has been revealed by a theorem in section 2. It is proved that by knowing the first two polynomials of Hermite , Laguerre and Legendre polynomials the related differential equations could be retrieved by using the method based on the separated transformation of original basis and Forbenius covariants as projection operators. The examples in section 2 clarifies the details of this method. By using the Rodrigues formula as the separated operators acting on the original basis, we acquire the form of related differential equations. In section 3, we introduce the Lie algebra modules on vector space of polynomials. The s l 2 ,   R  and s l ( 2 ,   c )  has been known as the Lie algebras connected to symmetries in polynomial and monomial space [4,5]. We prove that the conjugation (similarity transformation) of generators of these Lie algebras, yields isomorphic algebras that their Cartan subalgebras are Hermite, Laguerre and Legendre differential operators. The raising and lowering operators has been introduced in section sections 3.1 and 3.2. In section 3.8, applying Baker-Campbell-Hausdorff formula [6] on the basis of these isomorphic algebras, gives new relations of Hermite and Laguerre polynomials and their generating functions. Section 3.6 proves and represents a general form of differential-operator representations of s l ( 2 ,   R ) . In section 3.7, we propose a technique for solutions of differential equations based on raising operators acquired from associated Lie algebras.

2. Separated Operators of Basis Transformation

Let V  be a n-dimensional vector space with basis vectors e 1   ,   e 2   ,     e n   . The linear operator that transforms these basis to another basis e ' 1   , e ' 2   ,     e ' n , normally is defined as a unique linear operator O  in the matrix form. In present theory we define a set of linear operators O 1   ,   O 2   ,   . . .   O n  ; each one acts separately on the corresponding basis as follows:
O i e i = e ' i
The result of the action of operator O  and the set of O i  on the initial basis are the same, but O i   as separated basis transformations allow to choose a wide range of O i   operators whose overall transformations are equivalent to the operation of O . On the other hand, in many problems such as Rodrigues type formulas the separated basis transformation O i   are more accessible than overall operator O . In the context of differential operators, O i  could be regarded as the operators that transform the initial basis (monomials) in polynomial space to another basis as for example we observe in Rodrigues formula for Laguerre polynomials as the solutions (eigenfunctions) for Laguerre differential equation:
L n = 1 n ! ( D 1 ) n   x n
This equation can be interpreted as the transformation of initial basis ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   )  to new basis L n  (Laguerre polynomials) by the action of the operator
O n = 1 n ! ( D 1 ) n
Respect to Rodrigues formula, for all related differential equations such as Legendre, Chebyshev and Bessel Equations, there are separate and independent operators for each basis. Therefor we can apply a set of operators O i  instead of a unique operator O  to transform the initial basis in polynomials space. This method obviates the need to find the unique linear operator O  with the same action on all initial bases. We will show the relation between O i  and O  in proposition 2.1 after defining the projection operators as follows.
Let introduce the projection operators P 1   ,   P 2   ,   . . .   ,   P n   by the definition:
P i V = V i e i   ;   P i e j = δ i j   e i
Where V V  is a vector expanded as:
V = i V i e i
From (2) we have: i P i = I  ; P i P j = P i δ i j   as the main condition for projection operators ( I  is identity operator).
Remark 1
projection operators defined in (2), are linear operators.
We show that by basis transformations according to (1), the projection operators P i  are transformed as:
P i ' = O i P i ( j O j P j ) 1
Theorem 2.1. the generalized form of projection operator under separated basis transformations O i e i = e ' i  is:
P i ' = O i P i ( j O j P j ) 1
Proof . respect to basis transformation O i e i = e ' i  and (2) we have:
e ' i = P i ' e ' i = P i ' O i e i = P i ' ( j O j P j ) e i
Again, with substitution e ' i = O i e i  (5) reads as:
O i e i = P ' i ( j O j P j ) e i
It is valid for all e i . with identity O i e i = O i P i e i   , equation (6) becomes:
O i P i = P ' i ( j O j P j )
Or :   P i ' = O i P i ( j O j P j ) 1
As we expected.
Proposition 2.1. The transformation of projection operator defined in (4) is equivalent to a similarity transformation
P i ' =   (   k O k P k ) P i ( j O j P j ) 1 = O i P i ( j O j P j ) 1
Proof . expansion of the first two terms on left side considering P i P j = P i δ i j   yields:
( k O k P k ) P i = O i P i )
Therefor we have the right side of (8).
Equation (8) implies the similarity transformation of P i  under the basis transformation made by operator O = i O i P i  . Therefor the operator O  is the linear operator for transforming all basis   e i . Operator O  also yields transformations of all linear operators O  in vector space V  under basis transformation   e ' i = O e i  by the similarity transformation O ' = O O O 1 = (   k O k P k ) O ( j O j P j ) 1
Remark 2
Equation (4) meets the projection operator conditions:
a)
i P i ' = i O i P i ( j O j P j ) 1 = ( i O i P i ) ( j O j P j ) 1 = I
P i ' = O P i O 1 = k O k P k ) P i ( j O j P j ) 1 = O i P i ( j O j P j ) 1
The action of O  and O j  on a basis e j  is the same.
O e j = i O i P i e j = O j P j e j = O j e j = e ' i
This implies that the action of O  and O j  on e j  is equivalent.
b) Respect to (4) and (9) we conclude:
P i ' P j ' =   O i P i ( k O k P k ) 1 O j P j ( l O l P l ) 1 = [ ( k O k P k ) P i ( j O j P j ) 1 ] [ ( k O k P k ) P j ( j O j P j ) 1 ]
Middle terms in right side of (7) reduce to identity operator and thus:
P i ' P j ' = ( k O k P k ) P i P j ( j O j P j ) 1
Recalling P i P j = P i δ i j  and (4) we obtain:
P i ' P j ' = O i P i ( j O j P j ) 1 δ i j   = P i ' δ i j  
This proves the idempotency of P i '  i.e., P i ' P i ' = P i '
Where P i '   denoted as posterior probability analogy.
Equation (4) is the unique formula for P i '  and other forms in spite of their validity for satisfaction of projection operator conditions i.e., equation (3), are not the right candidates. As an example, we may propose this formula for P i ' :
P i ' = ( j P j O j ) 1 P i O i
It is straightforward to investigate that this definition is compatible with conditions (3) but if we multiply both sides by   e ' i  we obtain:
( j P j O j ) P i ' e ' i = P i O i e ' i  
Respect to (1) and (2) we get:
( j P j O j ) e ' i = P i O i e ' i
One of the solutions results in a false outcome:
( j P j O j ) = P i O i
Proposition 2.2. The product of projection operators is associative.
Proof. Let projection operators P i '   and P i ' '  correspond to   O i   and O i '  as transformation groups of coordinates i.e.
P i ' = O i P i ( j O j P j ) 1
And P i ' ' = O i ' P i ' ( j O j '   P i ' ) 1  
Substitution of first relation into above equation results in:
P i ' ' = O i ' O i P i ( j O j P j ) 1   [ ( j O j '   O j P j ) ( j O j P j ) 1 ) ] 1   P i ' ' = O i ' O i P i ( j O j P j ) 1 [ ( j O j P j ) ( j O j '   O j P j ) 1 ]
After vanishing of two central terms:
P i ' ' = O i ' O i P i ( j O j '   O j P j ) 1
This is compatible with (4) by replacing   O i  with O i ' ' = O i ' O i . Therefore, the corresponding projection operator for two consecutive transformation O i  and O i '  is equivalent the projection operator of O i ' ' = O i ' O i  transformation.
  P i ' ' = O i ' ' P i ( j O i ' ' P j ) 1
Remark 3
As is proved, the operators j O j P j  and O  are equivalent operators. Respect to equation (9), the projection operator P i  transforms as a similarity transformation under the action of operator j O j P j , therefor the initial basis should be transformed by this operator and consequently j O j P j  and O  are equivalent.
We show the identity (4) is also valid in function space, where the linear projection operators are defined.
Proposition 2.3. Let V  be a n-dimensional function space over the real field F  with a set of orthogonal basis functions   φ i  and inner product defined on a closed interval a , b .The same definition in (2) can be applied on these basis
P i F = c i φ i = φ i a b φ i F x
Where c i = φ i , F = a b φ i F x  (14)
Are the coefficients in expansion of square integrable function   F  in the basis φ i  calculated by inner product of φ i  and F  over the interval a , b .
Proof . with the identity (4) we conclude:
P i ' ( j O j P j ) = O i P i
Then we have: P i ' ( j O j P j ) F = O i P i F
Respect to (13) and (14) we have:
P i ' ( j O j φ j φ j F x ) = O i φ i a b φ i F x
Regarding (1) we can choose φ ' i = O i φ i  as transformed basis that results in:
P i ' ( j φ ' j a b φ j F x ) = φ ' i a b φ i F x
With the definition (2) of projection operator we have:
P i ' φ ' j = δ i j   φ ' j
Therefor the equation (16) respect to (14) reads as:
φ ' i a b φ i F x = c i φ ' i c i φ ' i = c i φ ' i  
So, the identity (9) is valid for operators in function spaces.
Proposition 2.4. Differential operators with certain eigenvalues and eigenfunctions can be linearly expanded by their projection operators.
Proof. Let the differential operator D  is characterized by eigenfunctions relation:
D φ i = λ i φ i
The eigenfunctions φ i  are linearly independent and are the basis vectors, i.e.: P i φ j = δ i j   φ j  .
Where P i  is the projection on i-th subspace, then by the identity:
D φ i = λ i φ i = ( j λ j P j )   φ i
The validity of this equation for all φ i   yields:
D = j λ j P j
That proves the proposition.
Theorem 2.2. Let the initial basis   e i  correspond to some set of linearly independent non-homogenous polynomials such as the regular bases ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   ) . After transforming the bases by equation O i e i = e ' i  to new bases   e ' i   which correspond the new linearly independent polynomials Ρ n ( x ) , if D  denoted as the differential operator with e i  or equivalently   x n  (n-th exponent of x) as its eigenfunctions (or eigenvector), then the corresponding differential operator D '  with eigenfunctions Ρ n ( x )  can be obtained by the relation:
D ' = ( k λ ' k O k P k ) ( j O j P j ) 1
Where λ ' k  are eigenvalues of D ' .
Proof. Respect to equation (19) the expansion of   D '  in terms of P i '  reads as:
D ' = i λ ' i P i '
Where P j '  are projection operators onto the i-th subspace ( i.e., e ' i ). Substitution of P i '  in equation (21) by equation (4) results in:
D ' = i λ ' i O i P i ( j O j P j ) 1 = ( i λ ' i O i P i ) ( j O j P j ) 1
This proves the theorem.

Projection operators in terms of resolvents

Associated to any differential operator in Hilbert space there are projection operators in terms of their resolvent
i.e.
P i = c ν i d λ ( λ I D ) 1   and   P i ' = c ν ' i d λ ( λ I D ' ) 1
Using (21), (22) we obtain:
c ν ' i d μ μ I D ' = ( O i c ν i d λ λ I D ) ( j O j c ν j d λ λ I D ) 1
For a unique transformation O = O i  for all φ i we get:
c ν ' i d μ μ I D ' = O ( c ν i d λ λ I D ) O 1
with expansion of resolvents d λ λ I D  as a Neumann infinite series (polynomial), it is proved that the corresponding differential operator after action of operator O  on the base functions   φ i  as defined in proposition 2.1, can be presented by a similarity transformation:
D ' = O D O 1
Example 2.1 . Eigenfunctions of the differential operator D = d d x  could be found as φ n = e n x . Transforming by φ ' n = O φ n = x φ n = x e n x  , The resulting corresponding differential operator respect to proposition 2.1 after substituting O = x  reads as:
D ' = x D x 1 = x 1 x 2 + x 1 D = 1 x + D
Action of this operator on x e n x  gives:
( 1 x + D ) x e n x = e n x + e n x + n x e n x = n x e n x
Thus, the eigenfunctions of this operator are x e n x  as expected. Because of the similarity relations of operators D  and D ' , their eigenvalues are identical.
Theorem 2.3. Let the linearly independent monomials p m ( x )   and q m ( x )   of polynomials Ρ n x  and Q n x  of degree n  are connected by the operators O i  as defined in equation (1) i.e.,
q m x = O m p m ( x )
Denote P m  as projection operators that project functions of variable x  on basis p m ( x )  with the definition of equation (2) P m p n x = δ m n   p n x
Then the operator O = m O m P m  acts as umbral composition on polynomial Ρ n x .
Proof. Let expand Ρ n ( x )   in terms of monomial basis  p m ( x )  
Ρ n x = n a n m p m ( x )  
Then action of O  on Ρ n x  gives
O Ρ n x = i O i P i ( m a n m p m ( x ) ) = i O i a n i p i ( x ) = i a n i q i x
This implies that the action of O  on  Ρ n ( x )  replaces the monomials p m ( x )  with q i x  while the coefficients a n m  in the expansion remains unchanged. This means that the operation of O  is equivalent with umbral composition by the definition
Ρ n o Q = m a n m q m x
This definition coincides the action of O  on Ρ n ( x ) . Application of this theorem for finding the generating function of Hermite polynomials has been shown in section 3.1.
Forbenius covariant of operators
For other representation of projection operator in terms of differential operator we apply the Forbenius covariants [7] as projection operators (matrices) which are the coefficient of Sylvester's formula. For a differential operator D  in polynomial space, the projection operator on the one-dimensional eigenfunction subspaces are given by
P l = k = 1 n D λ k λ l λ k k l
These operators act on the functions in function space and yields their projections on basis φ i  which are the eigenfunctions of D  with corresponding eigenvalues λ i .
Similarity transformation
Respect to equation (9) and related proposition, if we substitute O i P i  with ( k O k P k )  respect to the identity, P i P j = δ i j   P i   we have:
P i ' = k O k P k P i ( j O j P j ) 1 = O i P i ( j O j P j ) 1
This equation is a similarity transformation of P i  under the operator k O k P k . This similarity transformation corresponds to the basis transformation O i e i = e ' i . Actually, k O k P k  as an operator O ^  transforms all basis e i  to   e ' i  and corresponds the coordinate transformation matrix. From this equation we can deduce similarity transformation for other operators provided that the operators in similarity transformation have common eigenvalues. Therefor the differential operators with identical eigenvalues could be related by similarity transformation. As an example, differential operator D = x d d x = x D  with basis (eigenfunction) φ n = x n  transforms to another differential operator D '  with eigenfunction   φ ' i  after the basis transformation φ ' i = O i φ i . Therefor we have the similarity transformation:
D ' = k O k P k D ( j O j P j ) 1
If all O j  are the same namely   O , (34) will be reduced to:
D ' = O k P k D ( O j P j ) 1 D ' = O D O 1
In these cases that the single operator transforms all bases, the exact closed form of related differential operator could be derived by this method. However, for cases with separate O i , the validity of the retrieved differential operator relies on the action on the first two polynomials as
we show in next sections. The following example clarifies the method.
Example 2.2 :
Let the vector space V  spanned by the linearly independent basis ( 1 ,   e x ,   e 2 x ,   )  which are the eigenfunctions of operator D = d d x  .If these basis transforms to the new set of basis by multiplying with   e x 2 2  i.e., ( e x 2 2 ,   e x 2 2 + x ,   e x 2 2 + 2 x ,   )  then the corresponding operator with these new basis as its eigenfunctions could be obtained by (34). In this case O k = O = e x 2 2 . Thus, the equation (34) reduces to:
D ' = O D O 1 D ' = e x 2 2 D   e x 2 2
The term D   e x 2 2  is not just the derivative of   e x 2 2  , but an operator that is equal to:
D   e x 2 2 = d d x (   e x 2 2 ) +   e x 2 2 D
Then we have: D ' = e x 2 2 D e x 2 2 = e x 2 2 [ d d x ( e x 2 2 ) + e x 2 2 D ]
D ' = e x 2 2 ( x e x 2 2 + e x 2 2 D ) D ' = ( x + D )
The eigenfunctions of this operator are e x 2 2 + n x  with eigenvalues n as expected. It is noteworthy to note that the expression for probabilist’s Hermite polynomial H e 1  with the definition:
H e 1 = e x 2 2   d d x   e x 2 2
Differs from D ' , because in this definition the term d d x   e x 2 2  is not an operator but merely the derivative of   e x 2 2 .
It is easy to prove that any function of D '  can be expanded in terms of D  as follows:
f D ' = O f ( D ) O 1
Separated basis transformation method based on Forbenius covariants
Another approach to find D '  in terms of D  and   O j  is to apply the Forbenius covariant operators as projection operators as mentioned in (32).
P k = l = 1 N D λ l λ k λ l l k
These operators are projection operators onto the k-th one-dimensional sub-space (basis) [8].
substituting these projectors in equation (34) results in
D ' = k O k l = 1 N D λ l λ k λ l D ( j O j l = 1 N D λ l λ j λ l ) 1 l j
N  denoted as the dimension of function or polynomial space.
This method in comparison with previous methods are more applicable because the calculation of inverse of a product of differential operators is easier than other methods.
In the following sections we introduce an applicable method to find D '  in terms of D . Taking into account the equation (22) we have:
D ' = ( i λ ' i O i P i ) ( j O j P j ) 1
If all O i  are the same i.e., O i = O , then (37) reduces to
D ' = O ( i λ ' i P i ) O 1
The condition of identical eigenvalues for D '  and D  is not necessary in equation (37) and the case of identical eigenvalues are special case of this equation. we apply this equation restricted to the first two polynomials i.e., two-dimensional polynomial space. Substitution of P i  in (37) by Forbenius covariants (35) yields an applicable method as we will show in examples. It is noteworthy to recall that the term j O j P j  stands for a linear operator (equivalent to a matrix) that transforms the basis ( 1 , x , x 2 ,   )  of polynomials space to another basis. For example, it transforms basis ( 1 , x , x 2 ,   )  to Hermite polynomial H e n  as new linearly independent basis by the techniques that is presented in next section.

Applications of separated basis transformation method

In the Sturm Liouville problem and related differential equations and their specific solutions such as Hermite, Laguerre, Legendre and Jacobi polynomials, the transformation of basis in function space seems to be an interesting subject. For example, transformation of basis   ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   )  under the multiple differentiation which is compatible with Rodrigues’ formula to derive Hermite polynomials, presented as follows:
H e n = e D 2 2 x n
Where D = d d x  .
In the case of Laguerre polynomials, we have the transformation:
L n = 1 n ! ( D 1 ) n x n
Respect to our theory, these transformations are compatible with the operator action of O n  separately on basis   ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   ) , for Hermite polynomial we have:
O n = O = e D 2 2
And for Laguerre polynomials:
O n = 1 n ! ( D 1 ) n
We introduce the operator x D  as the unique operator with basis ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   )  as its associated eigenfunctions with eigenvalues ( 0,1 , 2 ,   . . .   ) :
x D   x n = n x n
Therefor we can use the equation (38) to find the differential operator that its polynomials are determined by applying related O n  on basis ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   )  as in (39) and (40). By substitution of P k  and O n  in equations (36) and (37) and the Forbenius covariants (35) and D = x D  in equation (36) we recover the corresponding differential equations of eigenfunctions such as H e n  and L n  as presented in net examples. The presented technique uses the first two polynomials i.e., the two-dimensional space of polynomials with monomials of order 1 and 0. This facilitates the calculation of desired differential equations and shows that if an infinite set of polynomials present the eigenfunctions of a unique differential operator, then applying this technique for the first two polynomials gives the exact form of related differential equation. We clarify this method by the following proposition
Proposition 2.5. Let the set of linearly independent polynomials   P n  are the eigenfunctions of a differential operator D '  and the set of original basis [ 1   ,   B x ,   B 2 x   ,     ,   B n x ]   are the eigenfunctions of differential operator D . Then applying the Forbenius covariant operator defined in (35) and equations (36) and (37) for the first two polynomials (eigenfunctions) P 0  and P 1 , yields the corresponding differential operator D '  from D .
First, we prove this proposition for Rodrigues formula as the action of operators O n   on the initial basis   B n x  in polynomial space to transform them to new basis P n  that correspond to the desired differential operator D ' (i.e., differential equation) as its eigenfunctions.

Rodrigues’ formula as a special case of separated basis transformation

In this section we prove the compatibility of Rodrigues’ formula with our presented techniques and show that substitution of O n  in equation (43) by Rodrigues’ formula transformation, yields the corresponding differential operators and equations.
Proof
Due to the presented theory, we showed that if the bases   e n  of a vector space V  which are the eigenfunctions of differential operator D , are transformed separately by operators   O n , the transformed differential operator obeys the equation (43) i.e.,
D ' = ( i λ ' i O i P i ) ( j O j P j ) 1
We check the basis transformation by Rodrigues formula [3]:
P n = 1 ω D n [ ω B n x ]
Where ω  defined by the relation ω ω ' = A B ' B  with A  as a polynomial of first degree.
If we choose monomial   B n ( x )   as the original basis of vector space:
[ 1   ,   B x ,   B 2 x   ,     ,   B n x ]
The Rodrigues formula could be chosen as the action of operator:
O n = 1 ω D n [ ω   . ]
On these basis. Therefor it is a special case of separated basis transformation.
The suitable operator with eigenfunctions   B n x  can be presented as
B x B ' x D   B n x =   n B n x
Where B ' x  denoted as the derivative of B x  The B x B ' x D  should replace D  in (35):
P k = l = 1 N B x B ' x D λ l λ k λ l
Respect to (43) by replacing   O n  by Rodrigues formula and P n  by equation (46) we get:
D ' = ( i λ ' i 1 ω D n [ ω l = 1 N B x B ' x D λ l λ k λ l   ] ) ( j O j P j ) 1
Let O 1 = ( j O j P j ) 1  then we obtain:
D ' = ( i λ ' i 1 ω D n [ ω l = 1 N B x B ' x D λ l λ k λ l   ] ) O 1
Taking into account the 2-dimensional space, and using (44) and (47) we have:
λ ' 0 = 0   ,   O 1 = 1 ω D [ ω   . ]   ,   P 1 = B x B ' x D λ ' 1 = 1 λ ' 1 B x B ' x D
Thus the (47) reads as:
D ' = λ ' 1 1 ω D ω 1 λ ' 1 B x B ' x D O 1 = 1 ω D ω B x B ' x D O 1
D ' = ( ω ω ' B B ' D + B ' 2 B ' ' B B ' 2 D + B B ' D 2 ) O 1
If we assume ω ω ' = A B ' B  ( as a crucial assumption in Rodrigues formula) this equation reduces to:
D ' = A B ' 1 D + D B ' ' B B ' 2 D + B B ' D 2 O 1
D ' = ( A B ' D B ' ' B B ' 2 D + B B ' D 2 ) O 1
Acting both side on P 1  as the second eigenfunction of D ' , we have:
D ' P 1 = ( A B ' D B ' ' B B ' 2 D + B B ' D 2 ) O 1 P 1
The term O 1 P 1  equals B x , thus:
D ' P 1 = ( A B ' D B ' ' B B ' 2 D + B B ' D 2 ) B x
D ' P 1 = A B ' B ' ' B B ' 2 + B B ' D B ' = A B ' ' B B ' + B B ' ( B ' ' + B ' D ) D ' P 1 = B D + A D ' D 1 D P 1 = B D + A
The term   D P 1 will be a constant α , thus:
α D ' D 1 = B D + A Or :   α D ' = B D 2 + A D
This implies that Rodrigues formula gives the solutions (or eigenfunctions) of the differential operator B D 2 + A D  and related differential equation up to a constant coefficient   α . i.e.,
B D 2 + A D y = β y
The following examples clarify this technique for some polynomials.
Example 2.3: Laguerre differential equation
Let we intend to find the differential equation which corresponds to a set of linearly independent polynomials in variable x . For example, we are given a few first Laguerre polynomials i.e., ( 1 ,   1 x ,   )  and we know the operator that maps the standard basis ( 1 , x , x 2 ,   )  to Laguerre basis i.e., operator presented in (42).
We can recover the corresponding Laguerre differential equation (operator) via the formula:
D ' = ( i λ ' i O i P i ) ( j O j P j ) 1  
Proof in 2 dimension (first 2 polynomials)
We restrict calculation in 2-dimensional polynomial space with basis ( 1 , x ) .These polynomials are transformed by O i  to Laguerre polynomials in the same dimension i.e., ( 1 , x + 1 ) . Thus, the corresponding operator x D  will be transformed to Laguerre differential operator by the equation (43).
Substitution of O i  by equation (42) and taking λ ' i  as the eigenvalues of Laguerre differential equation in 2-dimensional space of polynomials and replacing projection operators   P i  for basis ( 1 , x )  by equation (35) into equation (37) results in:
D ' = ( i = 0 1 λ ' i O i P i ) ( j = 0 1 O j P j ) 1
We have λ ' i = λ i   and λ 0 = 0   ,   λ 1 = 1 . Then equation (44) reduces to:
D ' = O 1 P 1 ( O 0 P 0 + O 1 P 1 ) 1
By equations (35) and (42) we obtain:
O 0 = 1   ,   O 1 = D 1   ,   P 0 = l = 0 1 D λ 1 λ 0 λ 1 = x D 1 1   ,   P 1 = l = 0 1 D λ 0 λ 1 λ 0 = x D 1
Therefor we have:
O 1 P 1 = D 1 x D = D + x D 2 x D
And: D ' = O 1 P 1 ( O 0 P 0 + O 1 P 1 ) 1 =  ( D + x D 2 x D  )( D + x D 2 2 x D + 1 ) 1  (57)
If we denote the ( D + x D 2 x D )  as D , we can reduce the equation (46) as follows:
D ' = O 1 P 1 ( O 0 P 0 + O 1 P 1 ) 1 = D   (   D D + 1 ) 1
The term
O 0 P 0 + O 1 P 1 = D D + 1 = O ^
Is the linear operator which transforms the basis ( 1 , x )  to Laguerre basis ( 1 , x + 1 )  and vice versa If we restrict the action of operators to 2-dimensional polynomial space. Therefor we have:
O ^ 1 ( x + 1 ) = x
If we act both sides of (47) on another basis x + 1  we get as well:
D ' x + 1 = D O ^ 1 x + 1 = D x = D ( x + 1 )
Or briefly: D ' x + 1 = D ( x + 1 )  (61)
(Note that D . 1 = 0 )
The equation (61) implies that the action of both operators D '  and D  on basis ( 1 , x + 1 )  are identical and therefor the simplest form of operator D '  which its eigenfunctions are Laguerre polynomials and its related transformation operators are O i , reads as:
D ' = D = ( x D 2 x D + D )
This is the exactly the Laguerre differential equation with positive eigenvalues, i.e.:
x D 2 x D + D y = n y
Action of this operator on the first basis i.e., “1” gives 0 as the first eigenvalue and therefor the required conditions for validity of this differential operator are met.
Proof in 3 dimension (first 3 polynomials)
In 3-dimension with basis ( 1 , x + 1 ,   1 2 ( x 2 4 x + 2 )  of Laguerre polynomial and ( 1 , x , x 2 )  of original basis, considering eigenvalues   λ 0 = 0   ,   λ 1 = 1   ,   λ 2 = 2 ,   the D '  reads as:
D ' = ( i λ ' i O i P i ) ( j O j P j ) 1 = ( O 1 P 1 + 2 O 2 P 2 ) ( O 0 P 0 + O 1 P 1 + O 2 P 2 ) 1 D ' = ( O 1 P 1 + 2 O 2 P 2 ) O ^ 1
Here O ^ 1  denotes the last term in (64). Acting both side on basis ( x + 1 )  results in:
D ' ( x + 1 ) = ( O 1 P 1 + 2 O 2 P 2 ) O ^ 1 ( x + 1 )
Respect to (60) and the identity P 2 x = 0  we have:
D ' x + 1 = O 1 P 1 x
In this dimension P 1 can be find as:
P 1 = l = 0 2 D λ l λ 1 λ l = x D λ 0 1 λ 0 x D λ 2 1 λ 2 = x D x D 2 1 2 = x D ( x D 2 )
Then (65) reads as:
D ' x + 1 = ( D 1 ) x D ( x D 2 ) x D ' x + 1 = D 1 x D x D ' x + 1 = [ D 1 x D ] ( x + 1 ) D ' x + 1 = ( x D 2 x D + D ) ( x + 1 )
This proves:
D ' = ( x D 2 x D + D )
As the Laguerre differential operator.
Example 2.4: Hermite differential equation
The same technique could be applied to derive Hermite differential equation by the formula (37). Because all O n  that transforms basis ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   )  to Hermite polynomials are equal to   O  as is shown in (41), after getting P k  by (35) and substitute them in (37) we have:
D ' = ( i λ ' i O P i ) ( j O P j ) 1 = O ( i λ ' i P i ) ( O j P j ) 1 = O ( i λ ' i P i ) ( j P j ) 1 O 1
Respect to j P j = 1  we get:
D ' = O ( i λ ' i P i ) O 1
Expanding the sum for eigenvalues λ ' i = λ i = 0,1   and substitution of O  by (41)we have:
D ' = e D 2 2 ( P 1 ) e D 2 2
From (35) we calculate P 1  as:
P 1 = l = 0 1 D λ 0 λ 1 λ 0 = D 0 1 0 = D
We know ( 1 ,   x ,   x 2 ,   x 3 ,   . . .   )  are the eigenfunctions of x D , thus by D = x D  we have:
D ' = e D 2 2 ( x D ) e D 2 2
This equation can be interpreted as a similarity transformation that maps x D  into D '  after basis changes. This will be hold just for the cases that eigenvalues are common between x D and D '  as we see in Hermite and Laguerre differential equations.
Expansion of e D 2 2  and   e D 2 2  results in:
D ' = ( 1 D 2 2 + D 4 8 + ) ( x D ) ( 1 + D 2 2 + D 4 8 + )
D ' = ( 1 D 2 2 + D 4 8 + ) ( x D + x D 3 2 + x D 5 8 + ) D ' = x D D 2 2 x D + D 4 8 x D + + ( 1 D 2 2 + D 4 8 + ) ( x D 3 2 + x D 5 8 + )
In the 2-dimensional space of polynomials the orders higher than 2 for D n  will be omitted, as it could be verified by action of both side on basis x  . By omitting the higher orders, we obtain:
D ' = ( 1 D 2 2 ) x D D ' = x D D 2 2 x D = x D 1 2 D ( D + x D 2 ) D ' = x D 1 2 ( D 2 + D 2 + x D 3 )
Omitting x D 3  results in:
D ' = x D D 2 = ( D 2 x D )
This is the well-known Hermit probabilist's Hermite differential operator with Hermite polynomial as its eigenfunctions and positive eigenvalues 0 , 1 , 2 , …. as its eigenvalues.
Example 2.5: Legendre differential equation
For Legendre polynomials we have:
P n = 1 2 n n ! D n ( x 2 1 ) n
That transforms the basis set S = 1 , ( x 2 1 )   ,   ( x 2 1 ) 2 ,    to Legendre polynomials. We can choose the appropriate operator D  whose eigenfunctions are these basis. Simply we write:
D = x 2 1 2 x D
Eigenfunctions of this operator are members of the set S .
The transforming operator is
O n = 1 2 n n ! D n
In this case the eigenvalues of D  and D '  (Legendre differential operator) are not identical and therefor the similarity transformation is not valid. However, we can apply the equation (37) after determining the P i  from (35).
For calculating P i  in 2-dimension, we have:
P 0 = l = 0 1 D λ 1 λ 0 λ 1 = x 2 1 2 x D 1 1 = 1 x 2 1 2 x D
P 1 = l = 0 1 D λ 0 λ 1 λ 0 = x 2 1 2 x D 0 1 0 = x 2 1 2 x D
Noe we get:
i λ ' i O i P i = 2 O 1 x 2 1 2 x D = 2 1 2 D x 2 1 2 x D = D x 2 1 2 x D
From equation (37) and (63) we get:
D ' = ( i λ ' i O i P i ) ( j O j P j ) 1 D ' = D x 2 1 2 x D ( j O j P j ) 1 = D x 2 1 2 x D O 1
Where O 1 = j O j P j 1 .
By action of both sides of (75) on basis x  as the second basis of Legendre polynomials in two dimension, we have:
D ' x = D x 2 1 2 x D O 1 x
Respect to O 1 x = x 2 1 , (65) reads as:
D ' x = D x 2 1 2 x D ( x 2 1 ) D ' ( D 1 D ) x = D x 2 1 2 x ( 2 x ) D ' D 1 ( D x ) = D ( x 2 1 ) D ' D 1 = D ( x 2 1 ) D ' = D ( x 2 1 ) D = D ( 1 x 2 ) D
Expansion of (77) reads as:
D ' = ( 1 x 2 D 2 2 x D ]
Which is the Legendre differential operator with positive eigenvalues n n + 1 .

3. Hermit, Laguerre, and Legendre differential operator as Cartan subalgebra of s l 2 , R  and s u ( 2 )

Let g l V  denote the linear transformation that maps vector space V  onto itself. In this section we present isomorphic Lie algebras to s l ( 2 , R )  defined by s l ( 2 , R )  module on vector space V  which is a linear map ϕ  defined by ϕ   :   s l 2 , R g l ( V )  that preserves the commutator relations of s l ( 2 , R )  algebra [4,8].
ϕ a , b = ϕ a , ϕ ( b ) a , b s l ( 2 , R )
This representation is s l ( 2 , R )  module on vector space V .
First, we review the structure of irreducible vector field representation of s l ( 2 , R ) .The generators of this algebra in matrix representation are as follows:
H = 1 2 1 0 0 1   ,   X = 0 0 1 0   ,   and   Y = 0 1 0 0
The commutation relations for this representation of s l ( 2 , R ) are:
X , Y   = 2 H   ,   H , X   = X   ,   H , Y   = Y
Let define a representation of s l ( 2 , R )  as its module on V  that preserves commutation relations by differential operators as its generators:
h = x D n 2   ,   e = D = x   ,   f = x 2 D n x
With the similar commutation relations
e ,   f = 2 h   ,   h ,   e = e   ,   h ,   f = f
The Cartan sub-algebra H = h  produces a decomposition of representation space:
V = V j
V j  are the eigenspace (eigenfunction) of generator h  as Cartan sub-algebra of s l 2 , R  and provide the solutions to the related differential equation.
h V j = j V j
In present paper the eigenspaces V j  are one dimensional and coincide the basis of polynomial space. These basis are called weight vectors. For a finite dimensional representation there is a highest weight j = n  that determines the dimension of representation space by d i m V = n + 1 . As an example, the Cartan subalgebra of s l 2 , R  can be represented by h = x D  with x n  as its weight vectors (eigenfunctions) and integer n  as eigenvalues. Due to the properties of s l ( 2 , R ) , the operator e   acts as lowering operator A  and f   as raising operator A + . The action of these operator on representation basis (eigenfunction) of h  lowers or raise the power of x n .
e V j = α V j 1   ,   f V j = β V j + 1
In the following sections we will construct a set of isomorphic Lie algebras to s l 2 , R  based on differential operators of Hermite, Laguerre and Legendre equations whose Cartan sub-algebras are Hermit and Laguerre differential operators. These algebras could be derived by similarity transformations (conjugation) of generators of s l ( 2 , R )  defined in equation (79). The similarity transformation is achieved by the transforming operator by which the original polynomial space basis transforms to the deemed polynomial i.e., Hermite, Laguerre and Legendre polynomials as transformed basis. These operators could be derived from Rodrigues’ formula as has been shown in previous examples. For each algebra there exist a set of lowering and raising operators that derives the recursion equations for related polynomials.

3.1. Associated Lie Algebra of Hermite Differential Operator

We search for a Lie algebra L H  isomorphic to s l ( 2 , R )  algebra with generators to be defined based on Hermite differential operators . Here we apply the transformation operator e D 2 2  as described in (41) for Hermite polynomials to derive similarity transformations (conjugation) of s l ( 2 , R )  bases as follows:
X 1 = e D 2 2 h   e D 2 2  , X 2 = e D 2 2 e e D 2 2  , X 3 = e D 2 2 f e D 2 2  (80) Equations (69) are the similarity transformations of Lie algebra L H , that results in an algebra with basis X i  isomorphic to L H .
Then for X 1  we have:
X 1 = e D 2 2 f e D 2 2
X 1 = e D 2 2 ( x D n 2 ) e D 2 2
Respect to (55) this equation reduces to:
X 1 = D H ' n 2
Where D H ' = x D D 2  as proved in (57), denoted as Hermite differential operator.
For   X 2  we get:
X 2 = e D 2 2 D e D 2 2
Since the operator D  is commutable with both e D 2 2  and e D 2 2 , we have:
X 2 = D e D 2 2 e D 2 2 = e D 2 2 e D 2 2 D = D
Similarly, for X 3 :
X 3 = e D 2 2 ( x 2 D n x ) e D 2 2 X 3 = e D 2 2 ( x 2 D ) e D 2 2 e D 2 2 ( n x ) e D 2 2 X 3 = e D 2 2 x 2 e D 2 2 D n e D 2 2 x e D 2 2
To calculate this generator, first we know from (57) that:
D H ' = e D 2 2 ( x D ) e D 2 2
Because D  commutes with   e D 2 2  we obtain:
D H ' = e D 2 2 x e D 2 2 D
Or : D H ' D 1 = e D 2 2 x e D 2 2
With: D H ' D 1 = x D D 2 D 1 = x D  (86)
Therefor we have: x D = e D 2 2 x e D 2 2
Multiplying this with itself results in:
( x D ) 2 = e D 2 2 x e D 2 2 e D 2 2 x e D 2 2 = e D 2 2 x 2 e D 2 2
With substitutions, equation (84) reads as:
X 3 = ( x D ) 2 D n ( x D )
Then the list for generators of this representation of s l 2 , R  is:
X 1 = D H ' n 2   ,   X 2 = D   ,   X 3 = ( x D ) 2 D n x D x D = ( x D ) ( D H ' n )
The Cartan subalgebra of this algebra is X 1 = D H ' n 2 .
Clearly these generators span the Lie algebra L H  isomorphic to s l ( 2 , R ) , which is a representation for an isomorphism of s l ( 2 , R ) . The commutation relations can be checked as:
X 1   ,   X 2 = D H ' n 2 D D D H ' n 2 = D = X 2
X 2   ,   X 3 = D [ x D ) 2 D n x D [ x D ) 2 D n x D D = 2   ( x D D 2 n 2 ) = 2 X 1
For X 1   ,   X 3  , first we note: X 3 = ( x D ) ( D H ' n ) , and we use D H '  instead X 1 without any change in commutator result. Thus, we have:
X 1   ,   X 3 = D H ' x D D H ' n x D D H ' n D H '
Due to the identity: D H ' n D H ' = D H ' D H ' n
The equation (80) becomes:
X 1   ,   X 3 = D H ' x D D H ' n x D D H ' D H ' n X 1   ,   X 3 = [ D H ' x D x D D H ' ] D H ' n X 1   ,   X 3 = [ D H ' x D x D D H ' ] D H ' n
Substitution of D H '  by x D D 2  gives:
X 1   ,   X 3 = [ D H ' x D x D ( x D D 2 ) ] D H ' n
Replacing operator x D  with its equivalence D x 1  results in:
X 1   ,   X 3 = [ D H ' x D x D ( D x 1 D 2 ) ] D H ' n X 1   ,   X 3 = [ D H ' x D + x D x D D ( x D ) ] D H ' n X 1   ,   X 3 = [ D H ' + 1 x D D ] x D D H ' n X 1   ,   X 3 = [ D H ' + 1 x D D ] x D D H ' n X 1   ,   X 3 = x D D H ' n = X 3
This proves the isomorphism of the Lie algebra L H  with basis X 1   ,   X 2 , X 3  with s l 2 , R .

Lowering and Raising operators of Hermite Polynomials and its Generating function

In this section we introduce the raising and lowering operators of Hermite polynomials which act on vector space representation of s l ( 2 , R ) . We denote raising and lowering operators as A + and A  respectively. These operators act on the weight vectors which are eigenfunctions of X 1 or D H '  i.e., the Hermite polynomials H e n . As an example, for Lie algebra L H  the following relations could be considered.
(1)
Due to the properties of s l ( 2 , R ) algebra the generator X 2 acts as a lowering operator A . This implies that:
D H e n = n H e n 1
(2)
Consecutive action of the X 1 and X 2 generators on the eigenfunction H e n of X 1 (i.e., the Hermite polynomial of degree n ) results in lowering of polynomial degree. Respect to (81):
X 1 X 2 H e n = ( D H ' n 2 ) D H e n = ( D H ' n 2 ) H e n 1 = n ( n 2 1 ) H e n 1
This means that the operator X 1 X 2  acts as a lowering (ladder) operator A   in the subspaces spanned by the Cartan subalgebra X 1  of L H .
(3)
The raising operator can be derived from equation (85) and (86):
D H ' = e D 2 2 x D e D 2 2 = e D 2 2 x e D 2 2 D D H ' D 1 = e D 2 2 x e D 2 2
If we act the right side of (94) on a Hermite polynomial of degree n ,  respect to equation (41)
we get:
e D 2 2 x e D 2 2 H e n = e D 2 2 x O 1 H e n = e D 2 2 x . x n = e D 2 2 x n + 1 = O x n + 1 = H e n + 1
Thus (94) and (95) yields:
D H ' D 1 H e n = x D D 2 D 1 H e n = ( x D ) H e n = H e n + 1
Therefor the operator x D  acts as raising operator   A +  in the associated vector space spanned by H e n .
(4)
If this method be repeated for X 1 X 3 operator, we have:
X 1 X 3 H e n = ( D H ' n 2 ) [ ( x D ) 2 D n x D ] H e n = ( D H ' n 2 ) ( x D ) 2 D H e n n ( x D ) H e n
Taking into account (95) and (96) we deduce:
X 1 X 3 H e n = ( D H ' n 2 ) ( x D ) 2 H e n 1 n H e n + 1 X 1 X 3 H e n = ( D H ' n 2 ) H e n + 1 n H e n + 1 X 1 X 3 H e n = ( n 2 + 1 ) H e n + 1 n H e n + 1 = ( 1 n 2 ) H e n + 1
Clearly the operator X 1 X 3  acts as a raising operator A + .
The results of this section can be used to derive recursive relations for Hermits polynomials as follows:
Any combination of operators involved in (92),(93),(95),(96) and (97) results in a recursive relation for Hermite polynomials.
(5)
The generating function of Hermit polynomial can be derived by a method based on theorem 3.2 as follows.
By expansion of e t x  and acting the operator O  defined in (41) on it and taking into account the umbral property of O  proved in theorem 2.3. We have
g x , t = O e t x = e D 2 2 m = 0 n t m x m m ! = m = 0 n t m m ! H e m
Recall that the e t x  are eigenfunctions of the operator e D 2 2
e D 2 2 e t x = 1 D 2 2 + .   .   . e t x = e t x t 2 2 e t x + .   .   . = e t 2 2 e t x
Therefor we can replace D  by t  in equation (*)
g x , t = e t 2 2 e t x = m = 0 n t m m ! H e m g x , t = e t 2 2 + t x = m = 0 n t m m ! H e m
This yields the Hermite polynomial generating function.

3.2. Associated Lie Algebra of Laguerre Differential Operator

For Laguerre polynomials the similarity transformation of the original basis of s l 2 , R will be obtained by operators in equation (42):
O n = 1 n ! ( D 1 ) n
For global transformation respect to the definition, we have:
O = n O n P n = n 1 n ! ( D 1 ) n P n Y 1 = O h O 1   ,   Y 2 = O e O 1   ,   Y 3 = O f O 1
These generators construct a Lie algebra L L isomorphic to both s l 2 , R and L H . Replacing h , e , f respect to (68) we get:
Y 1 = O ( x D n 2 )   O 1   ,   Y 2 = O D O 1   ,   Y 3 = O ( x 2 D n x ) O 1
For Y 1 due to equation (36), substituting D by ( x D n 2 ) and D ' by D L ' n 2   simply we obtain:
Y 1 = O x D n 2 O 1 = O x D O 1 n 2 = D L ' n 2
Because of complex structures of O and P n , we calculate the raising operator by recursive relation:
( n x D ) L n = n L n 1
We use the operator n x D as a lowering operator A . Substitution of L n by O x n gives rise to:
n x D O x n = n L n 1
Multiplying both side from the left by O 1 :
O 1 n x D O x n = n O 1 L n 1 O 1 n x D O x n = n x n 1
Action of left side on   x n equals the derivative of x n , then we have:
O 1 n x D O D Or   n x D = O D O 1
The left side should be replaced by its operator equivalent i.e.,
n x D L n = ( D L ' x D ) L n Thus :   Y 2 = O D O 1 = D L ' x D
For Y 3 we need O x O 1 :
O x O 1 O D O 1 = O x D O 1 = D L ' O x O 1 D L ' x D = D L ' Or   O x O 1 = D L ' D L ' x D 1
And for O x 2 D O 1 :
O x 2 D O 1 = D L ' D L ' x D 1 D L ' D L ' x D 1 ( D L ' x D ) O x 2 D O 1 = D L ' D L ' x D 1 D L '
Then Y 3 reads as:
Y 3 = D L ' D L ' x D 1 D L ' n D L ' D L ' x D 1 Y 3 = D L ' D L ' x D 1 ( D L ' n )
This operator acts as raising operator. Eventually for representation of s l 2 , R in basis of Laguerre polynomial and related differential is an algebra L L with generators:
Y 1 = D L ' n 2   ,   Y 2 = D L ' x D   ,   Y 3 = D L ' D L ' x D 1 ( D L ' n )
To prove the isomorphism of L L and s l 2 , R , first, we calculate the commutation relation Y 1   ,   Y 2   :
Y 1   ,   Y 2   = D L ' n 2 D L ' x D D L ' x D D L ' n 2 Y 1   ,   Y 2   = D L ' D L ' x D D L ' x D D L ' Y 1   ,   Y 2   = D L ' x D + x D D L '
We know: x D 2 + D , x D = x D 2 + D (102)
Because D L ' = x D 2 + D x D , after substitution in(100) we have:
D L ' + x D   , x D = D L   ' ,   x D   = x D 2 + D = D L   ' + x D = Y 2
Or: Y 1 , Y 2 = Y 2
This is compatible with s l 2 , R algebra.
For Y 2 , Y 3 we have:
Y 2 Y 3 = D L ' x D D L ' D L ' x D 1 ( D L ' n ) Y 2 Y 3 = D L ' D L ' D L ' x D 1 ( D L ' n ) x D D L ' D L ' x D 1 ( D L ' n )
Respect to (101) and (103), in second term, substitution of x D D L ' by Y 2 + D L ' x D yields:
x D D L ' D L ' x D 1 ( D L ' n ) = ( D L ' x D Y 2 ) D L ' x D 1 ( D L ' n ) = D L ' x D D L ' x D D L ' x D 1 ( D L ' n ) = D L ' x D D L ' x D D L ' x D 1 ( D L ' n ) = D L ' x D D L ' x D 1 ( D L ' n ) ( D L ' n )
Replacing second term of (***) by this, yields:
Y 2 Y 3 = D L ' D L ' D L ' x D 1 ( D L ' n ) D L ' x D D L ' x D 1 D L ' n + ( D L ' n ) = D L ' D L ' D L ' x D 1 ( D L ' n ) x D D L ' x D 1 ( D L ' n ) + ( D L ' n ) = D L ' D L ' D L ' x D 1 ( D L ' n ) x D D L ' x D 1 ( D L ' n ) + ( D L ' n ) = D L ' ( D L ' x D ) D L ' x D 1 ( D L ' n ) + ( D L ' n ) = ( D L ' + 1 ) ( D L ' n ) Y 2 Y 3 = ( D L ' + 1 ) D L ' n = D L ' 2 n 1 D L ' n
For Y 3 Y 2 we have:
Y 3 Y 2 = D L ' D L ' x D 1 ( D L ' n ) D L ' x D = D L ' D L ' x D 1 D L ' D L ' x D n D L ' D L ' x D 1 D L ' x D Y 3 Y 2 = D L ' D L ' x D 1 D L ' D L ' x D n D L '
Replacement of D L ' D L ' x D with relations of Y 1 , Y 2 gives:
Y 3 Y 2 = D L ' D L ' x D 1 ( Y 2 + D L ' x D D L ' ) n D L ' Y 3 Y 2 = D L ' D L ' x D 1 ( D L ' x D ( D L ' 1 ) n D L ' Y 3 Y 2 = D L ' ( D L ' 1 ) n D L ' = D L ' D L ' 1 n = D L ' 2 ( n + 1 ) D L '
Thus :
Y 2   ,   Y 3   = Y 2 Y 3 Y 3 Y 2 = 2 D L ' n = 2 D L ' n 2 = 2 Y 1
This proves isomorphism of L L and s l 2 , R as expected.

Lowering and Raising operators of Laguerre Polynomials and its Generating function

Applying the method used to derive lowering and raising operators for Hermite polynomial could be repeated for Laguerre polynomials too. Respect to the properties of Lie algebra L L , the generator Y 2 acts as lowering operator A and Y 3 acts as raising operator A + on the weight vectors L n which are the eigenfunctions of Y 1 or D L ' :
Y 2 L n = ( D L ' x D ) L n = n L n x D L n Y 2 L n = n x D L n = n L n 1
The action of Y 1 Y 2 on L n is also a lowering operator:
Y 1 Y 2 L n = ( D L ' n 2 ) ( D L ' x D ) L n Y 1 Y 2 L n = ( D L ' n 2 ) ( n x D ) L n = n ( D L ' n 2 ) L n 1 = n ( n 2 1 ) L n 1
To derive raising operator due to the equation (100) we have:
D L ' D L ' x D 1 = O x O 1
Action of both side on L n gives:
D L ' D L ' x D 1 L n = O x x n = O x n + 1 = L n + 1
Thus, the operator D L ' D L ' x D 1 acts as the raising operator A + in weight vector space of Laguerre polynomials.
Proposition 3.1
The generating function of Laguerre polynomial is derived by projection operator method.
proof
due to umbral properties of operator O = n O n P n , as we proved in theorem 2.3, we have:
g x , t = O 1 + x t + x 2 t 2 + = n O n P n 1 + x t + x 2 t 2 + = n O n x n t n = n L n t n
Substitution the series in x t powers with 1 1 x t and the identity e x D n e x = ( D 1 ) n gives
g x , t = n O n P n 1 1 x t = n 1 n ! e x D n e x ( P n 1 1 x t ) = n 1 n ! e x D n e x ( P n t n u n 1 x u )
If u 0 denoted as extractor coefficient operator for u 0 = 1 , Then the term P n t n u n 1 x u is equivalent to u 0 t n u n 1 x u = P n t n u n 1 x u
This yields
g x , t = n 1 n ! e x D n e x u 0 t n u n 1 x u = e x n 1 n ! ( t u ) n D n e x u 0 1   1 x u
Respect to Tylor series
f x + α = f x + α f ' x + α 2 2 ! f ' ' x +
We get
g x , t = e x n 1 n ! ( t u ) n D n e x u 0 1 1 x u = e x   e ( x + t u   ) u 0 1 1 ( x + t u   ) u g x , t = u 0 e ( t u   ) 1 ( x + t u   ) u = 1 1 t u 0 e ( t u   ) 1 x u 1 t
Expansion of the right side in terms of u with some algebra results in Laguerre generating function
g x , t = n L n t n = 1 1 t   e ( x t 1 t )

3.3. Associated Lie Algebra of Legendre Differential Operator

The main difference between Legendre differential operator and Hermite or Laguerre differential operator is its eigenvalues. For Hermite and Laguerre differential operators the eigenvalue are the same as the eigenvalues of original differential operator x D . The eigenvalues of x D are integers n
Correspond to eigenfunctions x n . The Hermite and Laguerre differential operators have the same eigenvalues and therefore we can apply the similarity transformation O x D O 1 to derive both operators from x D . Note that operator   O is defined specific for each differential operator. For Legendre differential operator the eigenvalues are n n + 1 which differs from eigenvalues of operator D =   x 2 1 2 x D whose eigenvalues are integers n and eigenfunctions are ( x 2 1 ) n . In this case we alter the original operator D to turn the same eigenvalues n n + 1 . This allows us to use similarity transformation O D O 1 to construct Legendre associated Lie algebra isomorphic to s l 2 , R . Let to add n 2 to D and act the result on the original basis ( x 2 1 ) n .
D + n 2 ( x 2 1 ) n = x 2 1 2 x D + n 2 ( x 2 1 ) n = n n + 1 ( x 2 1 ) n
Therefor we choose D + n 2 for similarity transformation of the form O D + n 2 O 1 . Now we search for a Lie algebra L P isomorphic to s l ( 2 , R ) algebra with generators to be defined based on Legendre differential operators. We define the following generators for Lie Algebra of Legendre Differential Operator.
Z 1 = O h ' O 1   ,   Z 2 = O e ' O 1   ,   Z 3 = O f ' O 1
The generators h ' , e ' , f ' are different from h , e , f defined for s l ( 2 , R ) in previous sections. These operators are defined to be compatible for original basis ( x 2 1 ) n . An isomorphic algebra to s l 2 , R with generators h ' , e ' , f ' represented as
h ' = x 2 1 2 x D + n 2   ,   e ' = D 2 x   ,   f ' = ( x 2 1 ) 2 2 x D n ( ( x 2 1 )
The commutation relations of these basis are:
h ' , e ' = x 2 1 2 x D + n 2 D 2 x D 2 x x 2 1 2 x D + n 2 = 1 4 x D 1 x D 1 x D x D = D 2 x = e '
For h ' , f ' we use the identity
h ' = 1 x 2 1 f ' + n h ' , f ' = 1 x 2 1 f ' + n   ,   f ' = 1 x 2 1 f '   ,   f ' + n   ,   f ' = 1 x 2 1 f '   ,   f ' h ' , f ' = 1 x 2 1 f ' f ' 1 x 2 1 f '
Some algebra shows
h ' , f ' = f '
With these commutation relations, respect to Jacobi identity we have
e ' , f ' = h '
This proves that generators h ' , e ' , f ' gives an isomorphic algebra to s l ( 2 , R ) . Based on these basis and conjugation them with operator O which is defined for Legendre polynomials in equation(71), we could derive its adjoint algebra with basis that are formed by Legendre differential operator. Due to (34) and common eigenvalues of and h ' and D L ' (not be confused with D L ' for Laguerre differential operator) we have
Z 1 = O h ' O 1 = O x 2 1 2 x D + n 2 O 1 = D L '
For another basis it is required to calculate O ( x 2 1 ) O 1 . The action of this operator on Legendre polynomial P n gives
O ( x 2 1 ) O 1 P n = O ( x 2 1 ) ( x 2 1 ) n = O ( x 2 1 ) n + 1 = P n + 1
This implies that O ( x 2 1 ) O 1   acts as raising operator and is equivalent to f '
f ' = O ( x 2 1 ) O 1
This equation and (113) gives
D L ' = O x 2 1 2 x D + n 2 O 1 = O x 2 1 2 x D O 1 + n 2 = O ( x 2 1 ) O 1 O 1 2 x D O 1 + n 2 D L ' = f ' ( O e ' O 1 )
Or Z 2 = O e ' O 1 = f ' 1 D L ' (114)
For Z 3 respect to (**) and (*) we have
Z 3 = O f ' O 1 = O ( x 2 1 ) 2 2 x D n ( ( x 2 1 ) O 1 Z 3 = O ( x 2 1 ) O 1 O x 2 1 2 x D O 1 n O ( x 2 1 ) O 1 Z 3 = f ' D L ' n 2 n f ' = f ' [ D L ' n n + 1 ]
Thus, the set of generators for Lie algebra of Legendre differential operator are as follows
Z 1 = D L '   ,   Z 2 = f ' 1 D L '   ,   Z 3 = f ' [ D L ' n n + 1 ]

3.4. Adjoint representation of s l 2 , c based on Hermite differential operator

An appropriate representation of s l ( 2 , c ) algebra presented as 9]:
h = 1 2 x D + 1 2   ,   e = i 2 D 2   ,   f = i 2 x 2
The commutation relations of these generators will be unchanged after omitting the imaginary i from e and f yields a representation of s l ( 2 , R ) with commutation relations of equation (79):
h = 1 2 x D + 1 2   ,   e = 1 2 D 2   ,   f = 1 2 x 2
The adjoint representation of elements of this Lie algebra, can be derived by conjugation with any element of the group S L ( 2 , R ) :
A d g X = g X g 1   ,   g S L ( 2 , R )
The element g could be derived by exponential map of generators of s l ( 2 , c ) :
g = e t X
assume X = 1 2 D 2 and t = 1 , then the adjoint representation elements will read as:
A d ( h ) = e D 2 2 h e D 2 2   ,   A d e = e D 2 2 e e D 2 2   ,   A d ( f ) = e D 2 2 f e D 2 2
Respect to equations (81) to (88):
A d h = e D 2 2 x D + 1 2 e D 2 2 = 1 2 D H ' + 1 2
A d e = e D 2 2 1 2 D 2 e D 2 2 = 1 2 D 2 A d f = e D 2 2 1 2 x 2 e D 2 2 = 1 2 ( x D ) 2
The eigenfunctions of h as Cartan subalgebra of s l ( 2 , R ) are x n . After conjugation with e D 2 2 , the adjoint representation’s Cartan subalgebra will be 1 2 D H ' + 1 2 with eigenfunctions or weight vectors 1 2 H e n . The transformation of x n to H e n , respect to (39) is given by the relation:
H e n = e D 2 2 x n = g x n g S L ( 2 , R )
Therefore, the conjugation of generators of algebra s l ( 2 , R ) by an element group g , results in an isomorphic adjoint algebra that its Cartan subalgebra’s weight vectors (eigenfunctions) could be derived by action of the same group element on the eigenfunctions of the original Lie algebra i.e., x n .
If we choose the exponent of generator f as group element g = e x 2 2   S L 2 , R   we have:
A d h = e x 2 2 1 2 x D + 1 2 e x 2 2 = e x 2 2 x D e x 2 2 + 1 2
Due to Example 2: 1 2 e x 2 2 x D e x 2 2 + 1 2 = 1 2 e x 2 2 x e x 2 2 e x 2 2 D e x 2 2 + 1 2 = 1 2 x ( D x ) + 1 2 This implies that the weight vectors of adjoint algebra should be ν n = x n e x 2 2 . And can be verified by the action of x ( D x ) on ν n .
For g = e t x 2 2 we get:
e t x 2 2 x D e t x 2 2 = e t x 2 2 x e t x 2 2 e t x 2 2 x t e t x 2 2 + e t x 2 2 D = x ( x t + D )
And e t x 2 2 x D e t x 2 2 = x ( D x t )
Thus, the eigenfunctions of this operator would be ν n = x n e t x 2 2 .

3.5. Representation of s u ( 2 ) and Hermite differential operator

Let introduce the basis a 1 , a 2 , a 3 of s u ( 2 ) given by
a 1 = 1 0 0 1   ,   a 2 = 0 1 1 0   ,   a 3 = 0 1 1 0
With commutation relations
a 2 , a 1 = a 3   ,   a 3 , a 2 = a 1   ,   a 1 , a 3 = a 2
These commutation relations coincide the complexified algebra of s u ( 2 ) that is the same as complexified s l ( 2 , R ) .
Comparing these basis with the generators of s l ( 2 , R ) presented in (78) reveals the relations
a 1 = 2 H   ,   a 2 = ( X + Y )   ,   a 3 = ( X Y )
Conjugation of these basis with an element of the group S L ( 2 , R ) gives the adjoint representation of s l ( 2 , R ) . Let use the operator introduced in (41) to derive Hermite polynomials from monomials x n . The similarity transformations
X 1 ' = O a 1 O 1 X 2 ' = O a 2 O 1 X 3 ' = O a 3 O 1
X 1 ' = 2 e D 2 2 H e D 2 2   ,   X 2 ' = e D 2 2 ( X + Y ) e D 2 2   ,   X 3 ' = e D 2 2 ( X Y ) e D 2 2
Substituting the basis H , X , Y by (78) and (79) gives
X 1 ' = 2 e D 2 2 ( x D ) e D 2 2   ,   X 2 ' = e D 2 2 ( D + x 2 D n x ) e D 2 2   ,   X 3 ' = e D 2 2 ( D x 2 D + n x ) e D 2 2
Thus, by equations (82) to (87) we get
X 1 ' = 2 ( D H ' n 2 )   ,   X 2 ' = D + ( x D ) ( D H ' n )   ,   X 3 ' = D ( x D ) ( D H ' n )
The commutation relations of these basis coincide the complexified algebra of s u ( 2 ) and as well s o ( 3 ) , the algebra of rotation group in 3-dimensional space.

3.6. General form of differential-operator representations of s l ( 2 , R )

Theorem 3.1
Denote by   B x any function of x and choose a set of its ordered integer exponents as linearly independent basis [ 1   ,   B x ,   B 2 x   ,     ,   B n x ] , then the set of generators
h = B B ' D   ,   e = D B '   ,   f = B 2 B ' D n B
Satisfy the commutation relations of s l ( 2 , R ) and yields an isomorphic algebra to it.
Proof
h , e = B B ' D D B ' D B ' B B ' D B B ' D B ' B ' '   B ' 2 D + D 2 B ' 1 B ' B ' 2 B ' '   B B ' 2 D + B B ' D 2 = D B ' = e h , f = B B ' D B 2 B ' D n B B 2 B ' D n B B B ' D = B 2 B ' D n B = f
By Jacobi identity, these two commutation relations imply the third commutation relation
e , f = 2 h
Thus, the above generators are representation of the algebra s l ( 2 , R ) based on an arbitrary linearly independent basis [ 1   ,   B x ,   B 2 x   ,     ,   B n x ] of polynomial space.
Assume these basis be transformed to new linearly independent basis P n by the equation
P n = O B n x
Where, O denoted as an operator that introduced in proposition (2.1 ) and equation (8) i.e., O = j O j P j acts on B n x as the n-th power of B x . Associated algebra of polynomials P n can be derived as the similarity transformation or adjoint representation of s l ( 2 , R ) as defined in examples. Note that the corresponding differential operator D P is derived by D P = O h O 1 . The generators of related associated algebra are
X 1 = O h O 1   ,   X 2 = O e O 1   ,   X 3 = O f O 1
In this setting O B O 1 can be acts as a raising operator for P n basis
O B O 1 P n = O B B n = O B n + 1 = P n + 1
Therefor we could apply this operator as raising operator A + .
O B O 1 = A +
By this substitution, The general form of generators could be derived
X 1 = O h O 1 = O B B ' D O 1 = D P D P = O B O 1 O D B ' O 1 = A + O D B ' O 1
( A + ) 1 D P = O D B ' O 1
Consequently, respect to (123) and (125) for X 2 we get
X 2 = O e O 1 = O D B ' O 1 = ( A + ) 1 D P
And for X 3
X 3 = O A + O 1 = O B 2 B ' D n B O 1 = O B O 1 O B B ' D O 1 n O B O 1 X 3 = A + D P n = A + D P n
Thus, the generators
X 1 = D P   ,   X 2 = ( A + ) 1 D P   ,   X 3 = A + D P n
Form an algebra L P as a representation of s l ( 2 , R ) .
The polynomials P n are the eigenfunctions of X 1 as weight vectors of Cartan subalgebra of L P .
As an example, the generators of Hermite algebra can be derived by this formula regarding the raising operator Of Hermite polynomials i.e., A + = x D
X 1 = D H ' X 2 = ( A + ) 1 D H ' = x D 1 D H ' = x D 1 x D D 2 = x D 1 x D D = D X 3 = ( x D ) ( D H ' n )
As it is expected.
The Lie algebra L P is the general form of representation of s l ( 2 , R ) whose weight vectors are eigenfunctions of arbitrary differential operator D P . This implies that for any differential equation with eigenfunction problem, we can apply the corresponding algebra L P and its raising operator to derive its solutions as described below.

3.7. Solutions to Differential equations by Raising operator method

In this section we apply the raising operators of the Lie algebra associated with differential operators defining the related differential equations to derive its solutions. We start with a known differential equation and first two solutions i.e., the first two eigenfunctions with the lowest eigenvalues. Then by the definition of raising operator   A + defined by (126), we derive this operator by restriction to 2 dimension of polynomial space and using the first two terms of j O j P j and Forbenius covariant operator, the entire eigenfunction (solutions) of the differential equation could be derived.
Example 3.1
As an example, for Laguerre differential equation, if we know the first two monomial i.e., L 0 = 1 and L 1 = x + 1 as the trivial eigenfunctions, respect to equation (126) the raising operator is
A + = O B O 1
Where operator O transforms the basis x n to Laguerre polynomials L n .
For Laguerre differential equation by B = x , the raising operator appears as
A + = O x O 1 = ( j O j P j ) x O 1
A + = O 1 P 1 x O 1
And acting both side on 1 as the first monomial we get
A + . 1 = O 1 P 1 x O 1 . 1
By O 1 . 1 = 1 and P 1 x = x and by the   O 1 = D 1 , this equation yields
A + . 1 = ( D 1 ) x 1 = ( A + ) 1 ( D 1 ) x
The action of operator ( A + ) 1 ( D 1 ) on x is the same as D , then we have the identity
( A + ) 1 D 1 = D ( A + ) 1 D 1 D 1 = 1 ( A + ) 1 1 D 1 = 1
And this gives
A + = 1 D 1
Applying this operator on the first two Laguerre polynomials gives the nth solution
L n = ( A + ) n . 1 = ( 1 D 1 ) n . 1
This method can be applied for any differential operator to find its eigenfunctions or ordered solutions.
Example 3.2
For Hermite differential equation to derive O 1 due to equation (37) for 2 dimension we have
D H ' = ( λ ' 0 O 0 P 0 + λ ' 1 O 1 P 1 ) ( j O P j ) 1
We assume λ ' 0 = 0 , λ ' 1 = 1 and O = ( j O P j ) 1
D H ' = O 1 P 1 O 1
Acting both side on first basis x definition for projection operator P 1   , gives
D H '   x = O 1 P 1 O 1 x D H '   x = O 1 P 1 x D H '   x = O 1 ( x D ) x
This equation shows both operators in the equation are equivalent
D H ' = O 1 ( x D )
Substitution for D H ' and action of D 1   on both sides, yields
( x D D 2 ) D 1 = O 1 ( x D ) D 1 x D = O 1 x
Or x D x 1 = O 1 (136)
Respect to A + = O 1 B O 1 we get
A + = x D x 1 x O 1
Acting both side on 1 as the first 1 basis
A + . 1 = x D x 1 x O 1 . 1
A + . 1 = x D x 1 x . 1
Thus, we have
A + = x D
With raising operator, we derive all Hermits eigenfunctions as solutions to its differential equation
H e n = ( A + ) n . 1 = ( x D ) n . 1

3.8. Baker-Campbell-Hausdorff formula application for Lie Algebras of Differential Operators

A specific version of Baker-Campbell-Hausdorff formula implies that if the commutator relation of a Lie algebra generators X 1   ,   X 2 meets the equation [6]:
X 1   ,   X 2 = s   X 2
With s R , then the BCH formula reduces to
e X 1 e X 2 = exp X 1 + s X 2 1 e s
Adjoint representation of s l 2 , c as defined in equations (116) and (118) represented by generators
A d h = 1 2 D H ' + 1 2   ,   A d e = 1 2 D 2   ,   A d f = 1 2 ( x D ) 2
That obey the commutation relations in equations (79)
1 2 D H '   , 1 2 D 2   = 1 2 D 2
Multiplying by 1   yields
1 2 D H '   , 1 2 D 2   = 1 2 D 2
Due to (139) and (140) we have
e D H '   e 1 2 D 2 = exp D H ' + D 2 2 ( 1 e )
Acting both sides on x n by equation (39) yields
e D H '   e 1 2 D 2 x n = exp D H ' + D 2 2 ( 1 e ) x n e D H ' H e n = exp D H ' + D 2 2 ( 1 e ) x n e n H e n = exp D H ' + D 2 2 ( 1 e ) x n H e n = exp D H ' n + D 2 2 ( 1 e ) x n
This is a new relation the converts x n to H e n and alternative to the classic relation:
H e n = e 1 2 D 2 x n
This technique is also applicable to other differential operators such as Laguerre and Legendre differential operators.

References

  1. Ismail, Mourad, Mourad EH Ismail, and Walter van Assche. Classical and quantum orthogonal polynomials in one variable. Vol. 13. Cambridge university press, 2005.
  2. Arfken, George B., and Hans-Jurgen Weber. "Mathematical methods for physicists." (1972).
  3. Rasala, Richard. "The Rodrigues formula and polynomial differential operators." Journal of Mathematical Analysis and Applications 84.2 (1981): 443-482.
  4. Post, Gerhard, and Nico van den Hijligenberg. "gl (λ) and differential operators preserving polynomials." Acta Applicandae Mathematica 44 (1996): 257-268.
  5. Turbiner, A. "Lie algebras and polynomials in one variable." Journal of Physics A: Mathematical and General 25.18 (1992): L1087. [CrossRef]
  6. Matone, Marco. "Closed form of the Baker–Campbell–Hausdorff formula for the generators of semisimple complex Lie algebras." The European Physical Journal C 76 (2016): 1-7. [CrossRef]
  7. Horn, Roger A., and Charles R. Johnson. "Topics in matrix analysis Cambridge university press." Cambridge, UK (1991).
  8. Howe, Roger E., and Eng Chye Tan. Non-Abelian Harmonic Analysis: Applications of SL (2,c). Springer Science & Business Media, 2012. [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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated