Open Access
Issue
Wuhan Univ. J. Nat. Sci.
Volume 27, Number 3, June 2022
Page(s) 201 - 210
DOI https://doi.org/10.1051/wujns/2022273201
Published online 24 August 2022

© Wuhan University 2022

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

0 Introduction

As we all know, there are two ways to describe a mechanical system, the Lagrangian Method and the Hamiltonian Method. The two methods are transformed through Legendre transformation. However, if the system studied is a singular one, there will exist inherent constraints when it is transformed from the Lagrangian to the Hamiltonian, and the system is called a constrained Hamiltonian system (also called constrained canonical system). The constrained Hamiltonian system has great relationship with gauge field theory, which dominates particle physics. Besides, the constrained Hamiltonian systems can also be applied to the theory of condensed matter and some other research areas. Therefore, the study of the constrained Hamiltonian system is a hot topic, and much progress has been made, especially on the aspect of symmetry. For example, the standard constrained canonical systems and their quantization, standard canonical symmetries, symmetries of the path integral, canonical symmetries for higher-order singular systems, etc. Readers who are interested, please see Refs. [1-3] for details.

Recently, fractional calculus has gained much attention and has been applied to many fields [4-10]. Because fractional derivatives are elegant tools to deal with the dissipative forces for the nonconservative systems [11,12], fractional calculus of variations was studied by many scholars, such as Klimek [13], Agrawal [14,15], Muslih et al [16], Rabei et al [17], Malinowska and Torres [18], and so on [19-21]. Especially, Frederico and Torres [22] studied the fractional variational problems within mixed integer and Riemann-Liouville fractional derivatives and established the following fractional Euler-Lagrange equation

Lqi-ddtLq˙i+DtDt2αLDt1Dtαqi=0Mathematical equation(1)

where L=L(t,qj,q˙j,Dt1Dtαqj)Mathematical equation, 0<α<1Mathematical equation, qjMathematical equation is the generalized coordinate, q˙j=dqj/dtMathematical equation is the generalized velocity, t1DtαqjMathematical equation means the Riemann-Liouville fractional derivative of qjMathematical equation, qjC2([t1,t2];R)Mathematical equation, L(,,,)C2([t1,t2]×Rn×Rn×Rn;R)Mathematical equation, i,j=1,2,,nMathematical equation.

After the differential equations of motion are established, one needs to find the solutions to them. Symmetry methods play the key role in solving differential equations [23-52]. Particularly, Lie symmetry, which means the invariance of the differential equations of motion, is one of the most important methods to help reduce the freedoms of differential equations (see Refs.[33-39] for a review). Therefore, in this paper, we intend to establish the Lagrange and Hamilton equations for the fractional singular Lagrangian L=L(t,qj,q˙j,Dt1Dtαqj)Mathematical equation, and then study Lie symmetry and conserved quantity for the corresponding dynamic system.

It is noted that Eq. (1) in Ref. [22] is not correct for the Lagrangian L=L(t,qj,q˙j,Dt1Dtαqj)Mathematical equation. Therefore, we begin with the fractional calculus of variations within mixed integer and Riemann-Liouville fractional derivatives. In this paper, we set 0<α<1Mathematical equation.

1 Preliminaries

Several definitions and properties of fractional calculus are listed [10, 14, 21].

Let f(t)Mathematical equation be a function, t(t1,t2)Mathematical equation, then the left and the right Riemann-Liouville fractional integrals are defined as

t1Dt-αf(t)=1Γ(α)t1t(t-τ)α-1f(τ)dτ, α>0Mathematical equation(2)

tDt2-αf(t)=1Γ(α)tt2(τ-t)α-1f(τ)dτ, α>0Mathematical equation(3)

the left and the right Riemann-Liouville fractional derivatives are defined as

t1Dtαf(t)=1Γ(n-α)(ddt)nt1t(t-τ)n-α-1f(τ)dτ, n-1α<nMathematical equation(4)

tDt2αf(t)=1Γ(n-α)(-ddt)ntt2(τ-t)n-α-1f(τ)dτ, n-1α<nMathematical equation(5)

the left and the right Caputo fractional derivatives are defined as

t1CDtαf(t)=1Γ(n-α)t1t(t-τ)n-α-1(ddτ)nf(τ)dτ, n-1α<nMathematical equation(6)

tCDt2αf(t)=1Γ(n-α)tt2(τ-t)n-α-1(-ddτ)nf(τ)dτ, n-1α<nMathematical equation(7)

The left and the right Riemann-Liouville fractional derivatives satisfy the following formulae for integration by parts

t1t2g(t)Dt1Dtαf(t)dt=t1t2f(t)DtCDt2αg(t)dt+j=0n-1t1Dtj+α-nf(t)(-D)n-1-jg(t)|t1t2Mathematical equation(8)

t1t2g(t)DtDt2αf(t)dt=t1t2f(t)Dt1CDtαg(t)dt-j=0n-1tDt2j+α-nf(t)(D)n-1-jg(t)|t1t2Mathematical equation(9)

where D=d/dtMathematical equation.

Specially, when α1Mathematical equation, we have

t1Dt1f(t)=Dt1CDt1f(t)=df(t)/dt=f˙(t) , DtDt21f(t)=DtCDt21f(t)=-df(t)/dt=-f˙(t)Mathematical equation(10)

2 Fractional Constrained Hamiltonian System

2.1 Fractional Euler-Lagrange Equation

Hamilton action within mixed integer and Riemann-Liouville fractional derivatives is

I[qj()]=t1t2L(t,qj,q˙j,Dt1Dtαqj)dt, j=1,2,,nMathematical equation(11)

The isochronous variation

δI=0Mathematical equation(12)

with qj(t1)=qj1Mathematical equation, qj(t2)=qj2Mathematical equation and δDt1Dtαqj=Dt1DtαδqjMathematical equation is called Hamilton principle within mixed integer and Riemann-Liouville fractional derivatives, where δMathematical equation means the isochronous variation.

From Eq. (12), using the integration by parts formula (Eq. (8)), we have

δI=t1t2(Lqiδqi+Lq˙iδq˙i+LDt1DtαqiδDt1Dtαqi)dt=0Mathematical equation(13)

where

t1t2(Lq˙iδq˙i)dt=t1t2(Lq˙iddtδqi)dt=t1t2[ddt(Lq˙iδqi)-(ddtLq˙i)δqi]dtMathematical equation

=-t1t2[(ddtLq˙i)δqi]dt+(Lq˙iδqi)|t1t2=-t1t2[(ddtLq˙i)δqi]dtMathematical equation(14)

t1t2(LDt1DtαqiδDt1Dtαqi)dt=t1t2(LDt1DtαqiDt1Dtαδqi)dtMathematical equation

=t1t2(δqiDtCDt2αLDt1Dtαqi)dt+(t1Dtα-1δqiLDt1Dtαqi)|t1t2Mathematical equation

=t1t2(δqiDtCDt2αLDt1Dtαqi)dt+L(t2)Dt1Dtαqit1t2[(t2-t)-αΓ(1-α)δqi(t)]dtMathematical equation(15)

Substituting Eqs. (14) and (15) into Eq. (13), considering the independence of δqiMathematical equation, i=1,2,,nMathematical equation and the arbitrariness of the interval [t1,t2]Mathematical equation, we get

Lqi-ddtLq˙i+DtCDt2αLDt1Dtαqi+L(t2)Dt1Dtαqi(t2-t)-αΓ(1-α)=0, i=1,2,,nMathematical equation(16)

where L(t2)Dt1Dtαqi=LDt1Dtαqi(t2,qj(t2),q˙j(t2),Dt1Dtαqj(t2))Mathematical equation. Eq. (16) is called the fractional Euler-Lagrange equation within mixed integer and Riemann-Liouville fractional derivatives.

Remark 1   Equation (16) is different from Eq. (1), which discards the second terms of Eq. (15).

For the fractional Lagrangian system (Eq. (16)), define the element HijMathematical equation of the Hessian matrix [Hij]Mathematical equation as

Hij=2Lq˙iq˙j, i,j=1,2,,nMathematical equation(17)

If det[Hij]=0Mathematical equation, the Hessian matrix is said to be degenerated and the corresponding Lagrangian is called singular. Furthermore, if the fractional Lagrangian system (Eq. (16)) is a singular one, then there exist inherent constraints when it is transformed into a fractional constrained Hamiltonian system.

2.2 Fractional Inherent Constraint

Define the generalized momenta and the Hamiltonian as

pi=L(t,qj,q˙j,Dt1Dtαqj)q˙i, piα=L(t,qj,q˙j,Dt1Dtαqj)Dt1DtαqiMathematical equation(18)

HC=piq˙i+piα(t1Dtαqi)-L(t,qj,q˙j,Dt1Dtαqj), i,j=1,2,,nMathematical equation(19)

In this paper, we assume that t1Dtαqi=xi(t,qj,q˙j,pjα)Mathematical equation, i,j=1,2,,nMathematical equation and rank[Hij]=RMathematical equation, 0R<nMathematical equation.

When 1R<nMathematical equation, we have

q˙A=fA(t,q,pα,pA,q˙B)Mathematical equation(20)

where q˙A=(q˙1,q˙2,,q˙R)Mathematical equation, fA=(f1,f2,,fR)Mathematical equation, q=(q1,q2,,qn)Mathematical equation, pA=(p1,p2,,pR)Mathematical equation, q˙B=(q˙R+1,q˙R+2,,q˙n)Mathematical equation and pα=(p1α,p2α,,pnα)Mathematical equation. From Eq. (18), we get

pi=L(t,qj,q˙j,Dt1Dtαqj)q˙i=gi(t,q,q˙A,q˙B,pα)=gi(t,q,pA,q˙B,pα), i=1,2,,nMathematical equation(21)

If i=1,2,,RMathematical equation, then Eq. (21) is reasonable. When i=R+1,R+2,,nMathematical equation, giMathematical equation will not depend on q˙BMathematical equation, i.e.,

pB=gB(t,q,pA,pα)Mathematical equation(22)

where pB=(pR+1,pR+2,,pn)Mathematical equation and gB=(gR+1,gR+2,,gn)Mathematical equation. We write Eq. (22) as

ϕB(t,q,p,pα)=pB-gB(t,q,pA,pα)=0Mathematical equation(23)

where ϕB=(ϕR+1,ϕR+2,,ϕn)Mathematical equation, p=(p1,p2,,pn)Mathematical equation. Renumber Eq. (23), we can obtain

ϕa(t,q,p,pα)=0, a=1,2,,n-R, 1R<nMathematical equation(24)

Equation (24) is the fractional inherent constraint, and we call it the fractional primary constraint. It is noted that the fractional primary constraint comes from the definition of the generalized momenta (Eq. (18)) rather than the fractional Euler-Lagrange equation (Eq. (16)).

Similarly, when R=0Mathematical equation, we can get

ϕa(t,q,p,pα)=0, a=1,2,,nMathematical equation(25)

Therefore, from Eqs. (24) and (25), we have

ϕa(t,q,p,pα)=0, a=1,2,,n-R, 0R<nMathematical equation(26)

Equation (26) is the inherent constraint and is called fractional primary constraint.

2.3 Fractional Constrained Hamilton Equation

From Eqs. (18) and (19), we obtain

δHC=q˙iδpi+δpiαDt1Dtαqi-LqiδqiMathematical equation(27)

where

Lqi=ddtLq˙i-DtCDt2αLDt1Dtαqi-L(t2)Dt1Dtαqi(t2-t)-αΓ(1-α)=p˙i-DtCDt2αpiα-piα(t2)(t2-t)-αΓ(1-α)Mathematical equation(28)

Because HC=HC(t,q,p,pα)Mathematical equation, we have

δHC=HCpiδpi+HCqiδqi+HCpiαδpiαMathematical equation(29)

From Eq. (26), we have

δϕa(t,q,p,pα)=ϕaqiδqi+ϕapiδpi+ϕapiαδpiα=0Mathematical equation(30)

Making use of the Lagrange multipliers λa(t)Mathematical equation, a=1,2,,n-RMathematical equation, 0R<nMathematical equation, from Eqs. (27) - (30), we get

q˙i=HCpi+λaϕapi, Dt1Dtαqi=HCpiα+λaϕapiα, p˙i=DtCDt2αpiα+piα(t2)(t2-t)-αΓ(1-α)-HCqi-λaϕaqiMathematical equation(31)

Equation (31) is called the fractional constrained Hamilton equation (also called fractional constrained canonical equation) within mixed integer and Riemann-Liouville fractional derivatives.

Remark 2   When the element t1DtαqjMathematical equation in L(t,qj,q˙j,Dt1Dtαqj)Mathematical equation is omitted, the fractional constrained Hamilton equation (Eq. (31)) will reduce to the classical constrained Hamilton equation in Ref. [2].

2.4 Consistency Condition

In this section, the main purpose is to solve the Lagrange multipliers λa(t)Mathematical equation, a=1,2,,n-RMathematical equation, 0R<nMathematical equation.

Let F=F(t,q,p,pα)Mathematical equation, G=G(t,q,p,pα)Mathematical equation, then Poisson bracket is defined as

{F,G}=FqiGpi-FpiGqiMathematical equation(32)

Let HT=HC+λaϕaMathematical equation, a=1,2,,n-RMathematical equation, 0R<nMathematical equation, then we have

{ϕa,HT}=ϕaqiq˙i-ϕapi(tCDt2αpiα+piα(t2)(t2-t)-αΓ(1-α)-p˙i)Mathematical equation(33)

It follows from ϕ˙a=0Mathematical equation, a=1,2,,n-RMathematical equation, 0R<nMathematical equation, that

{ϕa,HC}+λb{ϕa,ϕb}+ϕat+ϕapiαp˙iα+ϕapi(tCDt2αpiα+piα(t2)(t2-t)-αΓ(1-α))=0Mathematical equation(34)

where a,b=1,2,,n-RMathematical equation, 0R<nMathematical equation, i=1,2,,nMathematical equation. Eq. (34) is called the consistency condition of the fractional primary constraint.

3 Lie Symmetry and Conserved Quantity

3.1 Lie Symmetry

Lie symmetry means the invariance of the differential equations of motion (Eq. (31)) under the infinitesimal transformations

t¯=t+Δt, q¯i(t¯)=qi(t)+Δqi, p¯i(t¯)=pi(t)+Δpi, p¯iα(t¯)=piα(t)+ΔpiαMathematical equation(35)

whose expanding expressions are

t¯=t+θξ0(t,q,p,pα)+ο(θ), q¯i(t¯)=qi(t)+θξi(t,q,p,pα)+ο(θ)Mathematical equation

p¯i(t¯)=pi(t)+θηi(t,q,p,pα)+ο(θ), p¯iα(t¯)=piα(t)+θηiα(t,q,p,pα)+ο(θ)Mathematical equation(36)

where θMathematical equation is a small parameter, ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation are called infinitesimal generators of the infinitesimal transformations. In this paper, we only consider the linear part of θMathematical equation and neglect its higher order.

Expanding Eq. (31) and denoting

q˙i=si(t,q,p,pα)Mathematical equation(37)

t1Dtαqi=hi(t,q,p,pα)Mathematical equation(38)

p˙i=DtCDt2αpiα+fi(t,q,p,pα)Mathematical equation(39)

For Eq. (37), we have

q¯˙i-si(t¯,q¯,p¯,p¯α)=dq¯idt¯-si(t+Δt,q+Δq,p+Δp,pα+Δpα)Mathematical equation

=q˙i-si(t,q,p,pα)+θ[ξ˙i-q˙iξ˙0-X(0)(si)]Mathematical equation(40)

where X(0)=ξ0t+ξkqk+ηkpk+ηkαpkαMathematical equation, k=1,2,,nMathematical equation.

For Eq. (38), we have

t¯1Dt¯αq¯i-hi(t¯,q¯,p¯,p¯α)=Dt1Dtαqi-hi(t,q,p,pα)Mathematical equation

+θ[t1Dtα(ξi-q˙iξ0)-X(0)(hi)+ξ0ddtDt1Dtαqi-1Γ(1-α)ddt((t-t1)-αqi(t1)ξ0(t1))]Mathematical equation(41)

where ξ0(t1)=ξ0(t1,q(t1),p(t1),pα(t1))Mathematical equation.

For Eq. (39), we have

p¯˙i-Dt¯CDt¯2αp¯iα-fi(t¯,q¯,p¯,p¯α)=p˙i-DtCDt2αpiα-fi(t,q,p,pα)+θ[η˙i-ξ0ddtDtCDt2αpiαMathematical equation

-p˙iξ˙0-DtCDt2α(ηiα-p˙iαξ0)+1Γ(1-α)(t2-t)-αp˙iα(t2)ξ0(t2)-X(0)(fi)]Mathematical equation(42)

where ξ0(t2)=ξ0(t2,q(t2),p(t2),pα(t2))Mathematical equation.

For the fractional primary constraint Eq. (26), we have

ϕa(t¯,q¯,p¯,p¯α)=ϕa(t,q,p,pα)+θX(0)(ϕa)Mathematical equation(43)

Lie symmetry requires that the coefficients of θMathematical equation are zero in Eqs. (40) - (43), then we obtain

ξ˙i-q˙iξ˙0=X(0)(si), Dt1Dtα(ξi-q˙iξ0)+ξ0ddtDt1Dtαqi-qi(t1)ξ0(t1)Γ(1-α)ddt(t-t1)-α=X(0)(hi),Mathematical equation

η˙i-p˙iξ˙0-ξ0ddtDtCDt2αpiα-DtCDt2α(ηiα-p˙iαξ0)+1Γ(1-α)(t2-t)-αp˙iα(t2)ξ0(t2)=X(0)(fi)Mathematical equation(44)

and

X(0)(ϕa)=0Mathematical equation(45)

Equation (44) is called determined equation and Eq. (45) is called limited equation.

Definition 1   For the fractional constrained Hamiltonian system within mixed integer and Riemann-Liouville fractional derivatives (31), if the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation satisfy the determined equation (44), then the corresponding symmetry is called Lie symmetry.

Definition 2   For the fractional constrained Hamiltonian system within mixed integer and Riemann-Liouville fractional derivatives (31), if the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation satisfy both the determined equation (44) and the limited equation (45), then the corresponding symmetry is called weak Lie symmetry.

However, if we consider the deduction process of the fractional constrained Hamilton equation (31), then another condition will be exposed on the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation, and a new Lie symmetry will be defined.

Substituting Eq. (36) into Eq. (30), we get

ϕaqi(ξi-q˙iξ0)+ϕapi(ηi-p˙iξ0)+ϕapiα(ηiα-p˙iαξ0)=0Mathematical equation(46)

Equation (46) is called additional limited equation.

Definition 3   For the fractional constrained Hamiltonian system within mixed integer and Riemann-Liouville fractional derivatives (31), if the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation satisfy the determined equation (44), the limited equation (45) and the additional limited equation (46), then the corresponding symmetry is called strong Lie symmetry.

3.2 Conserved Quantity

Lie symmetry may not lead to a conserved quantity. However, under some conditions, a conserved quantity can be deduced from Lie symmetry. Of course, different conditions lead to different results.

Definition 4   A quantity IMathematical equation is called a conserved quantity if and only if dI/dt=0Mathematical equation.

Theorem 1   For the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation, which meet the determined equation (44), if there exists a gauge function G(t,q,p,pα)Mathematical equation satisfying the following structural equation

piξ˙i+piαDt1Dtα(ξi-q˙iξ0)+(piαddtDt1Dtαqi-HCt)ξ0-HCqiξi+λaϕapiαηiαMathematical equation

+λaϕapiηi-qi(t1)ξ0(t1)Γ(1-α)piαddt(t-t1)-α+(piαDt1Dtαqi-HC)ξ˙0+G˙=0Mathematical equation(47)

then the following Lie symmetry conserved quantity exists for the fractional constrained Hamiltonian system (31)

I=piξi+(piαDt1Dtαqi-HC)ξ0+{piαDt1Dτα(ξi-q˙iξ0)-(ξi-q˙iξ0)[τCDt2αpiα+piα(t2)(t2-τ)-αΓ(1-α)]}dτMathematical equation

-qi(t1)ξ0(t1)Γ(1-α)t1tpiαddτ(τ-t1)-αdτ+G=constMathematical equation(48)

Proof   Using Eqs. (26), (31) and (47), we have

ddtI=p˙iξi+piξ˙i+(p˙iαDt1Dtαqi+piαddtDt1Dtαqi-HCt-HCqiq˙i-HCpip˙i-HCpiαp˙iα)ξ0Mathematical equation

-(ξi-q˙iξ0)[tCDt2αpiα+piα(t2)(t2-t)-αΓ(1-α)]-qi(t1)ξ0(t1)Γ(1-α)piαddt(t-t1)-α+G˙Mathematical equation

+(piαDt1Dtαqi-HC)ξ˙0+piαDt1Dtα(ξi-q˙iξ0)Mathematical equation

=HCqiξi-λaϕapiαηiα-λaϕapiηi+(p˙iαDt1Dtαqi-HCqiq˙i-HCpip˙i-HCpiαp˙iα)ξ0Mathematical equation

+p˙iξi-(ξi-q˙iξ0)[tCDt2αpiα+piα(t2)(t2-t)-αΓ(1-α)]Mathematical equation

=(ξi-q˙iξ0)(HCqi-DtCDt2αpiα-piα(t2)(t2-t)-αΓ(1-α)+p˙i)-λaϕapiαδpiα-λaϕapiδpi=0Mathematical equation

This proof is completed.

Theorem 2   For the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation, which meet the determined equation (44) and the limited equation (45), if there exists a gauge function G(t,q,p,pα)Mathematical equation satisfying the structural equation (47), then a weak Lie symmetry conserved quantity (Eq. (48)) exists for the fractional constrained Hamiltonian system (31).

Theorem 3   For the infinitesimal generators ξ0Mathematical equation, ξiMathematical equation, ηiMathematical equation and ηiαMathematical equation, which meet the determined equation (44), the limited equation (45) and the additional limited equation (46), if there exists a gauge function G(t,q,p,pα)Mathematical equation satisfying the structural equation (47), then a strong Lie symmetry conserved quantity (Eq. (48)) exists for the fractional constrained Hamiltonian system (31).

Remark 3   When the element t1DtαqjMathematical equation in L(t,qj,q˙j,Dt1Dtαqj)Mathematical equation is omitted, Theorem 1 - Theorem 3 will reduce to the integer order cases, which are consistent with the results in Ref. [39].

4 An Example

We study Lie symmetry for the following fractional singular Lagrangian system

L=q˙1q2-q1q˙2+(q1)2+(q2)2+12[(t1Dtαq1)2+(t1Dtαq2)2]Mathematical equation(49)

From Eq. (49) we have

p1=Lq˙1=q2, p2=Lq˙2=-q1, p1α=LDt1Dtαq1=Dt1Dtαq1, p2α=LDt1Dtαq2=Dt1Dtαq2Mathematical equation(50)

HC=p1q˙1+p2q˙2+p1αDt1Dtαq1+p2αDt1Dtαq2-L=12[(p1α)2+(p2α)2]-(q1)2-(q2)2Mathematical equation(51)

Therefore, there exist two fractional primary constraints

ϕ1=p1-q2=0, ϕ2=p2+q1=0Mathematical equation(52)

From Eq. (34), we get

λ1=-q2-12[tCDt2αp2α+p2α(t2)(t2-t)-αΓ(1-α)] , λ2=q1+12[tCDt2αp1α+p1α(t2)(t2-t)-αΓ(1-α)]Mathematical equation(53)

Then Eq. (31) gives the following fractional constrained Hamilton equation

q˙1=-q2-12[tCDt2αp2α+p2α(t2)(t2-t)-αΓ(1-α)] , q˙2=q1+12[tCDt2αp1α+p1α(t2)(t2-t)-αΓ(1-α)] ,Mathematical equation

p˙1=q1+12[tCDt2αp1α+p1α(t2)(t2-t)-αΓ(1-α)] , p˙2=q2+12[tCDt2αp2α+p2α(t2)(t2-t)-αΓ(1-α)] ,Mathematical equation

t1Dtαq1=p1α , Dt1Dtαq2=p2αMathematical equation(54)

The determined equation (44) gives

ξ˙1-q˙1ξ˙0=-ξ0t[12p2α(t2)(t2-t)-αΓ(1-α)]-ξ2Mathematical equation

ξ˙2-q˙2ξ˙0=ξ0t[12p1α(t2)(t2-t)-αΓ(1-α)]+ξ1Mathematical equation

t1Dtα(ξ1-q˙1ξ0)+ξ0ddtDt1Dtαq1-1Γ(1-α)q1(t1)ξ0(t1)ddt(t-t1)-α=η1αMathematical equation

t1Dtα(ξ2-q˙2ξ0)+ξ0ddtDt1Dtαq2-1Γ(1-α)q2(t1)ξ0(t1)ddt(t-t1)-α=η2αMathematical equation

η˙1-p˙1ξ˙0-ξ0ddtDtCDt2αp1α-DtCDt2α(η1α-p˙1αξ0)+1Γ(1-α)(t2-t)-αp˙1α(t2)ξ0(t2)Mathematical equation

=ξ0t[12p1α(t2)(t2-t)-αΓ(1-α)]+ξ1Mathematical equation

η˙2-p˙2ξ˙0-ξ0ddtDtCDt2αp2α-DtCDt2α(η2α-p˙2αξ0)+1Γ(1-α)(t2-t)-αp˙2α(t2)ξ0(t2)Mathematical equation

=ξ0t[12p2α(t2)(t2-t)-αΓ(1-α)]+ξ2Mathematical equation(55)

The limited equation (45) gives

-ξ2+η1=0, ξ1+η2=0Mathematical equation(56)

The additional limited equation (46) gives

-ξ2+η1+(q˙2-p˙1)ξ0=0, ξ1+η2-(p˙2+q˙1)ξ0=0Mathematical equation(57)

Taking calculation, we find that

ξ0=1, ξ1=ξ2=0, η1=η2=0, η1α=η2α=0, G=0Mathematical equation(58)

satisfy the structural equation (47) and they also meet the determined equation (55) under the condition p1α(t2)=p2α(t2)=0Mathematical equation. Therefore, Lie symmetry conserved quantity can be obtained from Theorem 1 as

I=t1t(p1αddτDt1Dταq1+p2αddτDt1Dταq2-q˙1DτCDt2αp1α-q˙2DτCDt2αp2α)dτMathematical equation

-[12(p1α)2+12(p2α)2+(q1)2+(q2)2]Mathematical equation(59)

Because the solution (58) also meets the limited equation (56) and the additional limited equation (57), Eq. (59) is also a strong Lie symmetry conserved quantity.

Specially, if the element t1DtαqjMathematical equation in L(t,qj,q˙j,Dt1Dtαqj)Mathematical equation is omitted, then this example will reduce to the integer order case, in which the strong Lie symmetry conserved quantity has the form

I=(q1)2+(q2)2Mathematical equation(60)

5 Conclusion

Lie symmetry and conserved quantity are studied for the fractional constrained Hamiltonian system within mixed integer and Riemann-Liouville fractional derivatives here. The study of the fractional singular system just begins, and a lot of work is deserving to be done, such as, other symmetry methods, singular systems with only fractional derivatives, singular systems on time scales, and so on.

References

  1. Li Z P. Classical and Quantal Dynamics of Contrained Systems and Their Symmetrical Properties [M]. Beijing: Beijing Polytechnic University Press, 1993(Ch). [Google Scholar]
  2. Li Z P. Contrained Hamiltonian Systems and Their Symmetrical Properties [M]. Beijing: Beijing Polytechnic University Press, 1999(Ch). [Google Scholar]
  3. Li Z P, Jiang J H. Symmetries in Constrained Canonical Systems [M]. Beijing: Science Press, 2002. [Google Scholar]
  4. Kusnezov D, Bulgac A, Dang G D. Quantum Lévy processes and fractional kinetics [J]. Phys Rev Lett, 1999, 82: 1136-1139. [NASA ADS] [CrossRef] [Google Scholar]
  5. Naber M. Time fractional Schrdinger equation [J]. J Math Phys, 2004, 45: 3339-3352. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  6. Muslih S I, Agrawal O P, Baleanu D. A fractional Dirac equation and its solution [J]. J Phys A Math Theor, 2010, 43(5): 055203. [NASA ADS] [CrossRef] [Google Scholar]
  7. Metzler R, Klafter J. The random walk's guide to anomalous diffusion: A fractional dynamics approach [J]. Phys Rep, 2000, 339(1): 1-77. [CrossRef] [Google Scholar]
  8. Herrmann R. Gauge invariance in fractional field theories [J]. Phys Lett A, 2008, 372(34): 5515-5522. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  9. Lazo M J. Gauge invariant fractional electromagnetic fields [J]. Phys Lett A, 2011, 375(41): 3541-3546. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  10. Wu Q, Huang J H. Fractional Calculus [M]. Beijing: Tsinghua University Press, 2016(Ch). [Google Scholar]
  11. Riewe F. Nonconservative Lagrangian and Hamiltonian mechanics [J]. Phys Rev E, 1996, 53(2): 1890-1899. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
  12. Riewe F. Mechanics with fractional derivatives [J]. Phys Rev E, 1997, 55(3): 3581-3592. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  13. Klimek M. Lagrangian and Hamiltonian fractional sequential mechanics [J]. Czech J Phys, 2002, 52(11): 1247-1253. [NASA ADS] [CrossRef] [Google Scholar]
  14. Agrawal O P. Formulation of Euler-Lagrange equations for fractional variational problems [J]. J Math Anal Appl, 2002, 272(1): 368-379. [CrossRef] [MathSciNet] [Google Scholar]
  15. Agrawal O P. Fractional variational calculus in terms of Riesz fractional derivatives [J]. J Phys A Math Theor, 2007, 40(24): 6287-6303. [CrossRef] [Google Scholar]
  16. Muslih S I, Baleanu D. Hamiltonian formulation of systems with linear velocities within Riemann-Liouville fractional derivatives [J]. J Math Anal Appl, 2005, 304(2): 599-606. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  17. Rabei E M, Nawafleh K I, Hijjawi R S, et al. The Hamilton formalism with fractional derivatives [J]. J Math Anal Appl, 2007, 327(2): 891-897. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  18. Malinowska A B, Torres D F M. Introduction to the Fractional Calculus of Variations [M]. London: Imp Coll Press, 2012. [CrossRef] [Google Scholar]
  19. Zhang H B, Chen H B. Generalized variational problems and Birkhoff equations [J]. Nonlinear Dyn, 2016, 83(1-2): 347-354. [CrossRef] [Google Scholar]
  20. Song C J, Agrawal O P. Hamiltonian formulation of systems described using fractional singular Lagrangian [J]. Acta Appl Math, 2021, 172: 9. [CrossRef] [Google Scholar]
  21. Luo S K, Xu Y L. Fractional Birkhoffian mechanics [J]. Acta Mech, 2015, 226(3): 829-844. [CrossRef] [MathSciNet] [Google Scholar]
  22. Frederico G S F, Torres D F M. Fractional Noether's theorem with classical and Riemann-Liouville derivatives [C]// Proceedings of the 51st IEEE Conference on Decision and Control . Washington D C: IEEE Press, 2012: 6885-6890. [Google Scholar]
  23. Song C J, Zhang Y. Conserved quantities for Hamiltonian systems on time scales [J]. Appl Math Comput, 2017, 313: 24-36. [MathSciNet] [Google Scholar]
  24. Song C J, Zhang Y. Noether symmetry and conserved quantity for fractional Birkhoffian mechanics and its applications [J]. Fract Calc Appl Anal, 2018, 21(2): 509-526. [CrossRef] [MathSciNet] [Google Scholar]
  25. Zhang Y, Tian X. Conservation laws for Birkhoffian systems of Herglotz type [J]. Chin Phys B, 2018, 27(9): 090502. [NASA ADS] [CrossRef] [Google Scholar]
  26. Zhang Y, Tian X. Conservation laws of nonholonomic nonconservative system based on Herglotz variational problems [J]. Phys Lett A, 2019, 383(8): 691-696. [CrossRef] [MathSciNet] [Google Scholar]
  27. Zhang Y. Lie symmetry and invariants for a generalized Birkhoffian system on time scales [J]. Chaos, Solitons and Fractals, 2019, 128: 306-312. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  28. Ding J J, Zhang Y. Conserved quantities and adiabatic invariants of fractional Birkhoffian system of Herglotz type [J]. Chin Phys B, 2020, 29(4): 044501. [NASA ADS] [CrossRef] [Google Scholar]
  29. Zhou Y, Zhang Y. Noether symmetries for fractional generalized Birkhoffian systems in terms of classical and combined Caputo derivatives [J]. Acta Mech, 2020, 231(7): 3017-3029. [CrossRef] [MathSciNet] [Google Scholar]
  30. Zhai X H, Zhang Y. Mei symmetry and new conserved quantities of time-scale Birkhoff's equations [J]. Complexity, 2020, 2020: 1691760. [Google Scholar]
  31. Mei F X. Analytical Mechanics (II) [M]. Beijing: Beijing Institute of Technology Press, 2013(Ch). [Google Scholar]
  32. Mei F X. Symmetry and Conserved Quantity of Constrained Mechanical Systems [M]. Beijing: Beijing Institute of Technology Press, 2004(Ch). [Google Scholar]
  33. Mei F X, Wu H B. Dynamics of Constrained Mechanical Systems [M]. Beijing: Beijing Institute of Technology Press, 2009. [Google Scholar]
  34. Mei F X, Wu H B, Zhang Y F. Symmetries and conserved quantities of constrained mechanical systems [J]. Int J Dynam Control, 2014, 2(3): 285-303. [CrossRef] [Google Scholar]
  35. Cai P P, Fu J L, Guo Y X. Lie symmetries and conserved quantities of the constraint mechanical systems on time scales [J]. Reports on Mathematical Physics, 2017, 79(3): 279-298. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  36. Han Y L, Wang X X, Zhang M L, et al. Lie symmetry and approximate Hojman conserved quantity of Appell equations for a weakly nonholonomic system [J]. Nonlinear Dyn, 2013, 71(3): 401-408. [CrossRef] [Google Scholar]
  37. Chen X W, Li Y M, Zhao Y H. Lie symmetries, perturbation to symmetries and adiabatic invariants of Lagrange system [J]. Phys Lett A, 2005, 337(4-6): 274-278. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  38. Ding N, Fang J H. Lie symmetry and conserved quantities for nonholonomic Vacco dynamical systems [J]. Commun Theor Phys, 2006, 46(2): 265-268. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  39. Zhang Y, Xun Y. Lie symmetries of constrained Hamiltonian system with the second type of constraints [J]. Acta Phys Sin, 2001, 50(5): 816-819. [CrossRef] [Google Scholar]
  40. Song C J, Shen S L. Noether symmetry method for Birkhoffian systems in terms of generalized fractional operators [J]. Theor Appl Mech Lett, 2021, 11(6): 100298. [CrossRef] [Google Scholar]
  41. Jia Q L, Wu H B, Mei F X. Noether symmetries and conserved quantities for fractional forced Birkhoffian systems [J]. J Math Anal Appl, 2016, 442(2): 782-795. [CrossRef] [MathSciNet] [Google Scholar]
  42. Xia L L, Chen L Q. Conformal invariance of Mei symmetry for discrete Lagrangian systems [J]. Acta Mech, 2013, 224(9): 2037-2043. [CrossRef] [MathSciNet] [Google Scholar]
  43. Wang P, Xue Y. Conformal invariance of Mei symmetry and conserved quantities of Lagrange equation of thin elastic rod [J]. Nonlinear Dyn, 2016, 83(4): 1815-1822. [CrossRef] [Google Scholar]
  44. Jiang W A, Luo S K. Mei symmetry leading to Mei conserved quantity of generalized Hamiltonian system [J]. Acta Phys Sin, 2011, 60(6): 060201. [Google Scholar]
  45. Zhou S, Fu H, Fu J L. Symmetry theories of Hamiltonian systems with fractional derivatives [J]. Sci China: Phys Mech Astro, 2011, 54(10): 1847-1853. [NASA ADS] [Google Scholar]
  46. Song C J. Mei symmetry for constrained mechanical systems on time scales [J]. Transactions of Nanjing University of Aeronautics and Astronautics, 2021, 38(2): 353-360. [Google Scholar]
  47. Song C J, Cheng Y. Noether symmetry method for Hamiltonian mechanics involving generalized operators [J]. Adv Math Phys, 2021, 2021: 1959643. [Google Scholar]
  48. Song C J, Cheng Y. Noether's theorems for nonshifted dynamic systems on time scales [J]. Appl Math Comput, 2020, 374: 125086. [MathSciNet] [Google Scholar]
  49. Zhang Y. Mei's symmetry theorems for non-migrated Birkhoffian systems on a time scale [J]. Acta Phys Sin, 2021, 70(24): 244501. [CrossRef] [Google Scholar]
  50. Zhang Y. Mei symmetry and conservation laws for time-scale nonshifted Hamilton equations [J]. Adv Math Phys, 2021, 2021: 7329399. [Google Scholar]
  51. Zhang Y, Tian X, Zhai X H, et al. Hojman conserved quantity for time scales Lagrange systems [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(10): 2829-2837. [Google Scholar]
  52. Zhang Y, Cai J X. Noether theorem of Herglotz-type for nonconservative Hamilton systems in event space [J]. Wuhan University Journal of Natural Sciences, 2021, 26(5): 376-382. [Google Scholar]

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.