Hamiltonian Formulation of Classical Fields with Fractional Derivatives
Total Page:16
File Type:pdf, Size:1020Kb
Hamiltonian Formulation of Classical Fields with Fractional Derivatives A. A. Diab 1, R. S. Hijjawi 2, J. H. Asad 3, and J. M. Khalifeh 4 1General Studies Department- Yanbu Industrial College- P. O Box 30436- Yanbu Industrial City- KSA 2Dep. of Physics- Mutah University- Karak- Jordan 3Dep. of Physics- Tabuk University- P. O Box 741- Tabuk 71491- Saudi Arabia 4Dep. of Physics-Jordan University- Amman 11942 - Jordan Abstract An investigation of classical fields with fractional derivatives is presented using the fractional Hamiltonian formulation. The fractional Hamilton's equations are obtained for two classical field examples. The formulation presented and the resulting equations are very similar to those appearing in classical field theory. Keywords: Fractional derivatives, Lagrangian and Hamiltonian formulation. 1 1. Introduction Fractional calculus is an extension of classical calculus. In this branch of mathematics, definitions are established for integrals and derivatives of arbitrary non-integer (even complex) order, such /1 2 as d f (t) . It has started in 1695 when Leibniz presented his analysis of dt /1 2 the derivative of order 1/2 . Then it is developed mainly as a theoretical aspect of mathematics, being considered by some of the great names in mathematics such as Euler, Lagrange, and Fourier, among many others. This branch of mathematics has experienced a revival of interest and has been used in very diverse topics such as fractal theory 1viscoelasticy 2, electrodynamics 3,4 , optics 5,6 , and thermodynamics 7. The literature of fractional calculus started with Leibniz and today is growing rapidly 8-14 . Recently fractional derivatives have played a significant role in physics, mathematics, engineering, and pure and applied mathematics in recent years 14-20 . Several attempts have been made to include non conservative forces in the Lagrangian and the Hamiltonian mechanics. Riewe 20,21 presented a new approach to mechanics that allows one to obtain the equations for non conservative systems using fractional derivatives. Eqab et al 22-24 developed a general formula for determining the potentials of arbitrary forces, conservative and non conservative, using the Laplace transform of fractional integrals and fractional derivatives. Recently, they 24 also constructed the Hamiltonian formulation of discrete and continuous fields in terms of fractional derivatives. The fractional Hamiltonian is not uniquely defined; it seems that there are several choices of fractional Hamiltonian giving the same classical limit, i.e., the same classical Hamiltonian. This present paper is a generalization of the above work on Hamilton's equations for classical fields with Riemann- Liouville fractional derivatives. This paper is organized as follows: In Sec.2, the definitions of fractional derivatives are discussed briefly. In Sec.3, the fractional form of Euler- Lagrangian equation in terms of functional derivative of the full Lagrangian is presented. In Sec.4, the fractional form of Euler-Lagrange equation in terms of momentum density is investigated. Section 5 is devoted to the equations of motion in terms of Hamiltonian density in fractional form. In section 6, we introduce two examples of classical fields leading to Schrödinger equation and Dirac equation in fractional 2 derivative forms. The work closes with some concluding remarks (section 7). 2. Definitions of Fractional Derivatives In this section two different definitions of the fractional derivatives (left and right Riemann- Liouville fractional derivatives) are discussed. These definitions are used in the Hamiltonian formulation and the solution of examples leading to the equations of motion of the fractional order. The left Riemann-Liouville fractional derivative (LRLFD) reads as 25: 1 d x f (τ ) Dα f (x) = ( ) n dτ . (1) a x Γ(n −α) dx ∫ (x −τ )α −n+1 a The right Riemann-Liouville fractional derivative (RRLFD) reads as 25: b α 1 d n f (τ ) x Db f (x) = (− ) ∫ α −n+1 dτ . (2) Γ(n −α) dx x (τ − x) Here α is the order of the derivative such that n −1≤ α ≤ n and is not equal to zero. If α is an integer, these derivatives are defined in the usual sense, i.e., d Dα f (x) = ( )α f (x) ; a x dx d Dα f (x) = (− )α f (x) ; α = 2,1 ,... (3) x b dx 3. Fractional form of Euler-Lagrangian Equation in terms of Functional derivative of the full Lagrangian L We start our formalism by taking the Lagrangian density to be a function of field amplitude ψ and its fractional derivatives with respect to space and time as: L L α β α β = (,ψaxD ψ , xb D ψ , at D ψ , tb Dt ψ ,) . (4) Euler-Lagrange equation for such Lagrangian density in fractional form can be given as 26 : ∂∂LLα β ∂ L α ∂ L β ∂ L +xbDα + ax D β + tb D α + at D β = 0 . ∂∂ψaxD ψ ∂ xb D ψ ∂ at D ψ ∂ tb D ψ (5) 3 Now we can write the full Lagrangian L as: 3 L= ∫ L d r . (6) Using the variational principle, we can write: δ∫Ldt= δ ∫∫L d3 rdt = ∫∫ ( δ L ) dtd 3 r = 0 . (7) Using Eq. (4), the variation of is: L L L l L ∂∂α ∂ β ∂ α δ=+ δψα δψ()axD + β δψ () xb D + α δψ () at D + ∂∂ψψaxD ∂ xb D ψ ∂ at D ψ ∂L β β δ(tD b ψ )= 0 . (8) ∂ tD b ψ Substituting Eq. (8) into Eq. (7), and using the following commutation relation, ∂Lα ∂ L α αδ()axD ψ= α ax D () δψ . (9) ∂axDψ ∂ ax D ψ we get, ∂∂L Lα ∂l β ∂ L α ∫∫ [()δψ+αaxD () δψ + β xb D () δψ + α δ () at D ψ + ∂∂ψψaxD ∂ xb D ψ ∂ at D ψ ∂L β 3 β δ(tD b ψ )] d r dt = 0 . (10) ∂ tD b ψ Integrating by parts the second and the third terms in Eq. (10) with respect to space, we get: ∂∂Lα L β ∂l ∂ L α 0[=++∫∫dtx D bα a D x β ][() dτδψ + ∫∫ dt α δa D t ψ + ∂∂ψψψaxD ∂ xb D ∂ at D ψ ∂L β β δ(tD b ψ )] d τ . (11) ∂ tD b ψ Now, we can take the integration over space in the previous equation and convert it into summation, thus: ∂∂LLα β ∂ L ∂ L α ∑++xbD ax D δψ ii δτ + ∑ δ( atii D ψ ) δτ + ∂∂ψψψDα ∂ D β ∂ D α ψ i ax xb ii at i ∂L β +∑ δ(tbiD ψ ) δτ i = 0 . (12) ∂ Dβψ i t b i 4 We can write Eq. (12) in terms of Lagrangian density as: L ∑[δ ]i δτ i = 0 . (13) i Where the left hand side in Eqs.(12 and 13) represents the variation of L (i.e. δL ) which is now produced by independent variations inδψ i , α β α δ ( a Dt ψ i ) and δ (t Db ψ i ) . Suppose now that all δψ i ,δ (a Dt ψ i ) and β δ (t Db ψ i ) are zeros except for a particular δψ j . It is natural to define the functional derivative of the full Lagrangian ( ∂/L) with respect to ψ , ( Dαψ ) and ( D βψ ) for a point in the j-th cell to the ratio of δL to δψ 27. a t t b j ∂/L δL = lim . (14) ∂/ψ δτ j →0 δψ j δτ j Using Eq. (12), and note that the left hand side represents we get: ∂/L ∂Lα ∂ L β ∂ L = +xbDα + ax D β . (15a) ∂/ψψ ∂ ∂axD ψ ∂ xb D ψ ∂/Lδ L ∂ L = = . (15b) αlim α α ∂/atDψδτ j →0 δτ jatj δ( D ψ ) ∂ at D ψ ∂/Lδ L ∂ L = = . (15c) βlim β β ∂/tbDψδτ j →0 δτ jtbj δ( D ψ ) ∂ tb D ψ Now, using Eq. (15) we can rewrite Eq. (5) Euler- Lagrange equation in terms of the full Lagrangian L using functional derivative in fractional form: ∂/L α ∂/L β ∂/L +t Db α + a Dt β = 0. (16) ∂/ψ ∂/ a Dt ψ ∂/ t Db ψ It is worth mentioning that for α, β → 1, Eq. (16) reduces to the usual Euler- Lagrange equation for the classical fields 27 . With the help of Eqs. (12 and 16), we can write the variation of full Lagrangian in terms of α β functional derivatives and variations of ψ , a Dt ψ and t Db ψ as: ∂/L ∂/L α ∂/L β 3 δL = ∫ (δψ ) + α δ (a Dt ψ ) + β δ (t Db ψ )d r . (17) ∂/ψ ∂/ a Dt ψ ∂/ t Db ψ 5 4. Fractional form of Euler-Lagrange Equation in terms of Momentum Density The right side fractional form of momentum can be written as 27 : a δL Pj = α . (18a) δ aDt ψ j Using Esq. (12 and 15b) we get: a ∂L ∂/ L Pj=αδτ j = α δτ j . (18b) ∂atjDψ ∂/ atj D ψ From Eq. (18b), we can define the right side form of momentum density πα as: ∂/L ∂ L (πα ) j =α = α . (19) ∂/atjDψ ∂ atj D ψ Now, taking the left fractional derivative for Eq. (19), one gets: ∂/L Dα (π ) = Dα t b α j t b α . (20) ∂/ a Dt ψ j Repeating same steps above for left side fractional form of momentum density π β , we get: ∂/L ∂ L (π β ) j =β = β . (21) ∂/tbjDψ ∂ tbj D ψ Now, taking the right fractional derivative for Eq. (21), one gets: ∂/L Dβ (π ) = Dβ a t β j a t β . (22) ∂/ t Db ψ j Now, substituting Eqs. (20 and 22) into Eq. (16), we get: ∂/L = −[ D β π + Dαπ ]. (23) ∂/ψ a t β t b α 6 The above equation represents the fractional form of Euler- Lagrange equation in terms of momentum density and the functional derivative of full Lagrangian.