Time Fractional Schrodinger Equation Revisited

The time fractional Schrodinger equation (TFSE) for a nonrelativistic particle is derived on the basis of the Feynman path integral method by extending it initially to the case of a “free particle” obeying fractional dynamics, obtained by replacing the integer order derivatives with respect to time by those of fractional order. The equations of motion contain quantities which have “fractional” dimensions, chosen such that the “energy” has the correct dimension [ML/T].The action S is defined as a fractional time integral of the Lagrangian, and a “fractional Planck constant” is introduced. The TFSE corresponds to a “subdiffusion” equation with an imaginary fractional diffusion constant and reproduces the regular Schrodinger equation in the limit of integer order. The present work corrects a number of errors in Naber’s work.The correct continuity equation for the probability density is derived and a Green function solution for the case of a “free particle” is obtained. The total probability for a “free” particle is shown to go to zero in the limit of infinite time, in contrast with Naber’s result of a total probability greater than unity. A generalization to the case of a particle moving in a potential is also given.


Introduction
There has been an explosive research output in recent years in the application of methods of fractional calculus [1][2][3][4][5][6][7][8][9][10][11][12][13] to the study of quantum phenomena .The wellknown Schrodinger equation with a first-order derivative in time and second-order derivatives in space coordinates was given by Schrodinger as an Ansatz.The Schrodinger equation has been generalized to (i) a space fractional Schrodinger equation involving noninteger order space derivatives but retaining first-order time derivative [14][15][16][17][18], (ii) a time fractional Schrodinger equation involving non-integer order time derivative but retaining the second-order space derivatives [19], or (iii) more general fractional Schrodinger equation where both time and space derivatives are of non-integer order [20][21][22][23][24][25][26].The fractional Schrodinger equation has also been obtained by using a fractional generalization of the Laplacian operator [20] and by using a fractional variation principle and a fractional Klein-Gordon equation [36].In all these cases the fractional derivatives employed have been the regular fractional derivatives of the Riemann-Liouville type or the Caputo type (generally used in physical applications with initial conditions) which are both nonlocal in nature.
The fractional derivative which is nonlocal by definition can be made "local" by a limiting process as shown by Kolwankar and Gangal [41].Highly irregular and nowhere differentiable functions can be analyzed locally using these local fractional derivatives.The Heisenberg principle in the fractional context has been investigated using local fractional Fourier analysis [42].
The Schrodinger equation for a free particle has the appearance of a diffusion equation with an imaginary diffusion coefficient.This suggests a method of deriving the Schrodinger equation as has been done using the Feynman path integral technique [43][44][45] based on the Gaussian probability distribution in the space of all possible paths.In other words, the classical Brownian motion leads to the Schrodinger equation in quantum mechanics.As far as deriving the fractional Schrodinger equation is concerned, the path integral approach for the Brownian-like paths for the Levy stable processes which leads to the classical space fractional diffusion equation has been extended to the Levy-like quantum paths leading to the space fractional Schrodinger equation (SFSE) in the seminal papers of Laskin [14][15][16][17][18].It may be noted that in this case, the time derivative is still the integer first-order derivative; only the space part is of Advances in Mathematical Physics fractional order.The SFSE still retains the Markovian character and other fundamental aspects such as the Hermiticity of the Hamilton operator.Parity conservation and the current density have been explored in the space fractional quantum mechanics in terms of the Riesz fractional derivative.Applications of SFSE cover the dynamics of a free particle, particle in an infinite potential well, fractional Bohr atom, and the quantum fractional oscillator.Thus the space fractional Schrodinger equation appears to have been well established [18].This theory has been further generalized recently within the frame work of tempered ultradistributions [39].Thus the theory of SFSE can be considered fully established from the point of view of the Feynman path integral technique.
The fractional time derivative was introduced into the Schrodinger equation by Naber [19] by simply replacing the first-order time derivative by a derivative of non-integer order and retaining the second-order space derivatives intact.The resulting equation is referred to as the time fractional Schrodinger equation (TFSE).He did not derive the TFSE using the path integral or any other method.Naber carried out the time fractional modification to the Schrodinger equation in analogy with time fractional diffusion equation [19] but included the imaginary number  raised to a fractional power (the fractional degree being the same as the fractional order of the time derivative), implying a sort of the Wick rotation.In Naber's opinion [19], the TFSE is equivalent to the usual Schrodinger equation, but with a time-dependent Hamiltonian.He obtained the solutions for a free particle and a particle in a potential well.A lot of subsequent work has been done on the TFSE, mostly based on Naber's work [21][22][23]28], including its generalization into spacetime fractional quantum dynamics by including non-integer order derivatives in both time and space.Yet some basic questions have not been addressed.It has been observed that TFSE describes non-Markovian evolution and that the Green function in the form of the Mittag-Leffler function does not satisfy Stone's theorem on one-parameter unitary groups and the semiclassical approximation in terms of the classical action is not defined [31].There has been no derivation of TFSE on a basis similar to that of SFSE and it is the purpose of the present paper to rectify this lacuna.Since the path integral method of deriving the space fractional part of the Schrodinger equation is well established, the present paper concentrates only on deriving the time fractional Schrodinger equation from the Feynman path integral approach, leading to the time fractional Schrodinger equation as given by Naber.It may be pointed out that some results of Naber, such as the total probability being greater than unity, are difficult to understand physically.Moreover, several major errors in Naber's paper have gone unnoticed and in fact have been repeated by workers who have followed his work.Furthermore, some of these authors have introduced errors of their own.Since many of the conclusions in Naber's paper are based on derivations which include these errors, it calls for a reexamination of Naber's generalization to the time fractional Schrodinger equation.
The present paper derives the time fractional Schrodinger equation using the Feynman path integral technique.It concentrates on the time fractional part only and not on the space fractional part as the theory of the latter has been well established in the works of Laskin [14][15][16][17][18]. Furthermore, the paper considers only the Caputo-type nonlocal fractional derivatives and not the local fractional derivatives discussed earlier.The paper starts from a generalization of the classical dynamics into fractional dynamics of a free particle and then adapting the Feynman technique derives the correct equations for TFSE.It is demonstrated that Naber's result of probability being greater than unity is spurious and is a result of the ad hoc raising of the imaginary number  to a fractional power.The correct continuity equation for the probability density is also derived.The paper concludes with some new results.

Feynman Path Integral Method
The starting point for the Feynman method [43][44][45][46] is the classical Lagrangian  = (, ẋ , ) and the action  = ∫ (, ẋ , ) constructed from it.However, in view of the generalization to fractional calculus methods to be carried out later, the equations of motion of a classical particle in one dimension in the usual notation are considered first: Integrating with respect to time yields In the usual notation the Lagrangian is given by and the action is given by An outline of the Feynman path integral method is presented following very closely the account given by Feynman and Hibbs [43].The essence of the Feynman path integral approach to quantum mechanics is in the probability amplitude (also known as the propagator or the Green function) (  ,   ;   ,   ) for a particle starting from a position   at time   to reach a position   at a later time   , which arises from the contributions from all trajectories from   to   : where the contribution from each of the paths has the form Here  is the action defined in (5) and ℎ is Planck's constant, the quantum of action.The time integral of the Lagrangian is to be taken along the path in question.Restricting to one dimension, the probability amplitude can be written as The symbol D indicates the fact that the operation of integration is carried over all paths from  to .
The wave function (  ,   ) gives the total probability amplitude to arrive at   at   satisfying (9), where the integral is taken over all possible values of The kernel  can be computed by first carrying out a "time slicing" operation by dividing the time interval from   to   into  segments of duration where where is a line integral taken over the trajectory passing through the point ().The constant  is a normalizing factor.The Schrodinger equation for a free particle in one dimension is derived by considering a special case of (9), which describes the evolution of the wave function from a time   to a time   , when   differs from   by an infinitesimal amount  and applied to the case of a free particle.This step is based on the fact that the semiclassical approximation is valid not only in the limit of ℎ → 0 but also in the limit of small time interval [45].The Kernel is proportional to the exponential of (/ℎ) times the classical action for the infinitesimal time interval  =   −  .With an obvious change of notation   = ,   =  0 ,   = ,   =  +  and using the fact that the particle is free, (2) yields and ( 4) yields and ( 5) yields for the action As a consequence, (9) becomes If  differs appreciably from  0 , the exponential in (17) oscillates very rapidly and the integral over  0 contributes a very small value and only those paths which are very close to  give significant contributions.Changing the variable in the integral from  0 to  =  −  0 makes it ( 0 , ) = ( + , ).
Since both  and  are small quantities, (,  + ) may be expanded in Taylor's series and only up to terms of order  are retained.On the right-hand side, ( + , ) may be expanded in Taylor's series in powers of , retaining terms up to second order in  (the integral involving the first-order term vanishes).Then (17) becomes On the right-hand side the middle term vanishes on integration.It follows by equating the leading terms on both sides Hence Equating the remaining terms results in This can be recognized as the diffusion equation with an imaginary diffusion coefficient or the Schrodinger equation for a free particle in quantum mechanics.These considerations can be easily extended to the case of a particle moving in a potential field by incorporating a potential term (, ) in the Lagrangian  =  −  in (15).

Advances in Mathematical Physics
This will necessitate [43] incorporating an additional factor {1 − (/ℎ)(, )} in (18), which becomes Then (21) becomes Multiplying both sides by −ℎ/ results in the standard Schrodinger equation of quantum mechanics: It is to be noted that the imaginary number  on the lefthand side is not arbitrary; it arises from the coefficient of the potential term  in ( 23) but ultimately from the coefficient of /ℎ in (7).This minor detail becomes important as will be discussed later in Section 8.These considerations will be generalized for a particle obeying fractional dynamics.It can then be extended to the case of a particle in a potential field by including the potential term and making appropriate changes as will be described later.

Fractional Dynamics of a Free Particle
The first step is to generalize the equations of motion, (1)-(3), by replacing the integrals and derivatives by appropriate fractional integrals and derivatives.For physical problems with well-definable initial conditions the accepted practice is to employ the Caputo fractional derivatives [5].The Caputo derivative of order  is defined by [3]: In the limit  → , the Caputo derivative becomes the ordinary th derivative of the function.
The fractional integral of order  is defined by In generalizing the equations of motion, the second-order time derivative in Newton's law is replaced by a Caputo derivative of order , and the first-order derivative is replaced by a Caputo derivative of order (/2) [47].Then (2) and ( 3) become where as usual  0 ,   0 refer to the initial position and initial value of the "fractional momentum"   , respectively.It is to be noted that the variables ,  are still the space and time variables and have the dimensions of length and time [] and [], respectively.However, the dynamical quantities obtained by the operation of fractional derivation have different dimensions; for example, "fractional velocity" with the notation  0   /2  = ẋ /2 would have the units [/ /2 ].The dimension for the parameter "  " in the fractional momentum,   =   ẋ /2 , is no longer just [] but has to be chosen [47]

Time Fractional Schrodinger Equation for a Free Particle
There are two possible generalizations of the action integral in (5) used in fractional dynamics [11]: Since the Lagrangian has the dimensions of energy, the dimensions of action defined in ( 28) and ( 29),  I and  II , are different.The dimensions of  I are the same as that of the regular action, namely, [ 2 /], but that of  II is [ 2 / 2−/2 ].In the Newtonian limit  → 2,  II → regular action.These dimensional considerations have to be kept in mind in generalizing the Feynman method.In particular, if the choice from ( 29) is made, then a "fractional Planck constant" ℎ  with appropriate dimensions must be introduced in order to render the argument of the exponential in (7) dimensionless.
For a "free particle" (27) yields After carrying out the time slicing operation as in (11) and making the same approximation the evolution of the wave function in an infinitesimal interval of time  can now be obtained.Equation (30) yields and hence But for action, there are two choices: Making the appropriate changes, the equations for the evolution of the wave function in the two cases are Changing the variable in the integrals from  0 to  =  −  0 as before and introducing two constants Equations ( 36) and ( 37) can be written as The left-hand sides of ( 39) and ( 40) can be expanded in fractional Taylor's series [48] in time, with fractional derivative of order , and keeping only the lowest-order term in  yields  I,II (,  + ) =  I,II (, ) In the right-hand side a Taylor expansion with terms up to second order in  with respect to space can be carried out and the two cases will be considered separately.Thus, the equation for  I (, ) becomes On evaluating the integrals on the right-hand side of ( 42), the middle term with the first power of  vanishes.Equating the leading terms on both sides of ( 42 Equating the remaining terms requires that the powers of  must be the same on both sides; that is,  =  − 1. Inserting the value of  I and simplifying (43) yield Similarly, by expanding (40) the equation for  II (, ) becomes On evaluating the integrals on the right-hand side of (45) the middle term vanishes.Equating the leading terms on both sides of (45) as before yields for the normalizing factor In the remaining terms, the powers of  must be the same.This requires  = /2.Inserting the value of  II and simplifying (46) yield the equation for the wave function  II (, ): This completes the derivation based on the Feynman path integral method and ( 44) and ( 47) constitute the time fractional Schrodinger equations for a free particle corresponding to two ways of defining the action integral for a fractional dynamical system.In the limit  → 2, the fractional dynamical system goes over to the regular Newtonian system and in this case both (44) and ( 47) reduce to given earlier, thus recovering the standard Schrodinger equation for a free quantum particle [49].
It should be noted that since  ≤ 2, with  = 2 being the limiting case, the order of the time fractional derivative is ≤ 1 in both ( 44) and ( 47) and cannot exceed 1.This means the time fractional Schrodinger equation as derived from the path integral method always corresponds to the "subdiffusion" case in contrast to the case where the TFSE is obtained by a simple replacement of the first-order time derivative by a fractional order derivative [19].Furthermore, the order /2 of the fractional derivative corresponds to the first-order regular derivative as has been used above in Section 3. Thus it appears that the second method of defining the action leading to (47) is the natural way to generalize to TFSE.In appearance also it is as if the equation has been obtained by replacing all quantities in the Schrodinger equation by an equivalent fractional quantity, except the space derivative.Furthermore, (44) becomes a fractional order integro-differential equation when  < 1 and no longer just a fractional order differential equation.Because of these reasons, it is considered the method of choice to use (47) as the TFSE derived from the path integral method and no further reference will be made to (44).
Although it is possible to introduce additional parameters and cast (47) in a dimensionless form, it has not been done here.However, for convenience, the subscript II is dropped from the wave function; a simplified notation for the Caputo derivative,    , with  = /2 will be used.Naturally, 0 <  ≤ 1, and  → 1 yields the regular first-order time derivative, denoted by   .After incorporating these changes and defining a new constant   = ℎ  /2  , (47) becomes Equation ( 49) can be solved by a combination of the Fourier and Laplace transform methods [52][53][54].

Probability Current and the Continuity Equation
The probability density  is defined by  =  *  =  * .The complex conjugate wave function satisfies There is an identity satisfied by the Caputo derivative [5] where the right-hand side represents the regular first-order derivative.This identity can be used in studying the time derivative of the probability density, given by Inserting from (51) gives Substituting from ( 49) and ( 50), ( 53) yields Equation ( 54) can be rewritten after factoring out the constant and interchanging the order of space and time derivatives as Introducing a new function ψ =  1−  , (55) can be written as Defining a probability current density given by equation ( 56) can be written as This is the time fractional version of the continuity equation.
It may be noted that (58) differs from Naber's result (( 24) in [19]) and will be discussed later.

Solution for TFSE for a Free Particle
The solution for the TFSE for a free particle under the conditions (, 0) =  0 (); (, ) → 0, || → ∞,  > 0 is available in the literature but considered here for purposes of obtaining the Green function for the TFSE.
By applying the combined Fourier and Laplace transforms defined by equation ( 49) reduces to resulting in Applying the inverse Laplace transform, (61) yields in terms of the Mittag-Leffler function defined by a series or by the inverse Lapalce transform [3]: The Green function solution can be written as where the Green function is given by the inverse Fourier transform of the Mittag-Leffler function where the Mittag-Leffler function has been defined [3] before in (63).The Fourier inversion in (65) can be carried out [52-54] using the property that the Mittag-Leffler function is related through the Laplace integral to another special function of the Wright type denoted by where the Wright function is defined by [54]  (; , ) =  The Green function in (65) is then given by ) .
In the context of fractional diffusion, the function  /2 () belongs to the Wright type of probability densities characterized by the similarity variable  = ||/√   /2 , where   is the fractional diffusion coefficient, with ∫ ∞ 0  /2 () = 1.Furthermore, the probability densities are non-Markovian and exhibit a variance consistent with slow anomalous diffusion [49][50][51],  2  () = (2/Γ( + 1))    .In the limit of  → 1, the probability density function goes over to the Gaussian corresponding to regular diffusion.A plot of the reduced probability density function is given in Figure 1.
The Green function for regular diffusion describes a probability density, whereas the corresponding Green function for the Schrodinger equation is the propagator, which describes the probability amplitude for the particle to propagate from   at   to   at   .In exactly the same way, the Green function for time fractional diffusion describes a probability density, whereas the Green function in (68) is the fractional propagator and gives the probability amplitude.Of particular interest is the Fourier component of the wave function in (62) in connection with the total probability as  → ∞; the case discussed by Naber [19] and will be discussed in the next section.

Comments about Some Results in
Naber's Work [19] This section draws attention to some errors in Naber's paper which have gone unnoticed and have been reproduced repeatedly.Naber's equations will be referred to by their number and a prefix N. Naber explicitly states that he uses the Caputo derivative, which has been defined earlier in (25) in this paper.Although Naber does not give the explicit definition of the Caputo fractional derivative in his paper [19], it can be inferred from (NA.3) given in the appendix to his paper.
(a) However, Naber writes in (N16) for a Caputo derivative of order (1 − ]), reproduced here for convenience in (70): This is incorrect, as can be checked easily by taking the limit ] → 1.The left-hand side → (, ) as the derivative of zero order, but the right-hand side → 0 because of the Γ function in the denominator.The correct form of equation is As a consequence, the weight factor in (N18) should be ( − ) 1−] and not ( − ) −] , which Naber uses to give a physical significance to the entity he has introduced.The correct weight factor in (N18) would → 1 in the limit ] → 1.
(b) In (N11), Naber gives an identity satisfied supposedly by fractional Caputo derivatives of order less than 1, reproduced here for convenience in (72).
This identity is not correct and cannot be found anywhere.The identity satisfied by the Caputo derivatives is given in [5] and reproduced later in the current notation for ] < 1 This yields  1−]   ]  () = / =   () and has been used earlier in (52) in this paper.
The incorrect identity has been used by Naber to derive an equation for the probability current, which is obviously incorrect.Unfortunately, the incorrect identity, as given by Naber, in (72) has been repeatedly used in the literature.The correct equation for the probability current has been given in this paper in (59), which reduces to the standard continuity equation for the probability current in regular quantum mechanics [49].
(c) This point concerns the separation of the Mittag-Leffler function with an imaginary argument into an oscillatory part and a part which decays exponentially with time.There is nothing wrong with the derivation itself as given by Naber, and the function under discussion is the Fourier component of the free particle wave function, in his notation where  ] () in Naber's notation corresponds to the Mittag-Leffler function  ],1 () defined in (63).
The Mittag-Leffler function with the complex argument has been separated into an oscillatory part and a part based on the evaluation of the inverse Laplace transform along a Hankel contour and considering the contribution from the residue of the pole  0 =  1/]  together with the contribution from the integral along the two strips on either side of the branch cut, which is a standard procedure.Naber finally gives the solution as He argues that in the limit of  → ∞ the total probability arises basically from the first term and is equal to 1/] 2 , assuming that the wave function was initially normalized.Since ] < 1, the total probability is >1, a result difficult to understand physically.However, it is shown later that the solution derived in this paper yields a probability that → 0 in the limit  → ∞.
The solution corresponding to (74) is given by (62) in this paper.The separation into two parts corresponding to (76) can be carried out by considering the inverse Laplace transform of (61): As usual, the Bromwich contour is replaced by the Hankel contour.The two contributions arise from (i) the residue at the pole at  0 = (   2 ) 1/ (−) 1/ and (ii) the integral along the two strips from 0 to −∞ introduced along the branch cut.The latter yields a contribution which decays in time just as the second term in Naber's equation ( 76).The contribution from the residue is given by Residue and corresponds to the first term on the right-hand side in (76).However, ) . (79) Therefore, (78) yields ψ (, ) ψ (, 0) = Residue ψ (, 0) =  (   2 ) 1/ cos(/2)  −(   2 ) 1/ sin(/2). (80) The right-hand side of (80) has an amplitude term and an oscillatory factor.Because  < 1, cos(/2) is negative, the amplitude factor decays in time exponentially and → 0 in the limit  → ∞.Thus the contribution from the residue due to the pole also → 0. Thus the entire wave function → 0 in the limit of  → ∞, so does the probability density.Hence the total probability also → 0, in contrast to Naber's result.
The question naturally arises why there is this difference in the two results, both of which are concerned with the Mittag-Leffler functions with complex arguments.The reason appears to be that in Naber's case, the pole occurs at  0 =  1/] , whereas in the present paper, the pole occurs (using the same notation as Naber) at  0 =  1/] (−) 1/] .The simple  in Naber's case leads to the purely oscillatory solution and hence to the result of the probability being greater than unity.In our paper the imaginary number raised to the fractional power leads to the exponentially decaying solution and hence leads to the correct limit when  → ∞, namely, zero total probability.Naber's result is a direct consequence of his choice to raise the power of the imaginary number  to the fractional power, so as to incorporate a Wick rotation.However, as had been indicated earlier, this imaginary number cannot be arbitrarily altered as it is connected with the phase factor /ℎ in the Feynman propagator.If it is necessary to include a Wick rotation, the power of  should be changed to ] + 1 in Naber's equation (N9) instead of just ].If this is done, the total probability would properly → 0 in the limit  → ∞.If this is done, then the TFSE can be interpreted as the analytic continuation of the fractional diffusion equation, just as the regular Schrodinger equation can be considered as the analytic continuation of the regular diffusion equation.

TFSE for a Particle in a Potential Field
So far attention has been focused on a free particle.These considerations can be easily extended to the case of a particle moving in a potential field by incorporating a potential term (, ) in the Lagrangian  =  − .This will necessitate incorporating an additional term −(, ) /2 /Γ(1 + /2) in the right-hand side of (35) and an additional factor −(/ℎ  )((, ) /2 /Γ(1+/2)) in the exponential in (37).This results in changing (45) into One final remark needs to be made concerning the use of the Planck units for casting the TFSE in terms of dimensionless quantities.After the derivation, the FTSE in (83) can be cast in dimensionless quantities; however, the appropriate fractional Planck units must be defined to take care of the fractional quantities involved in (83).

Discussion and Conclusions
The TFSE has been derived using the Feynman path integral technique for a nonrelativistic particle.As expected the TFSE looks like a time fractional diffusion equation with an imaginary fractional diffusion constant but pertains to the realm of subdiffusion only, in contrast to Naber's generalization which includes superdiffusion as well.This is understandable because the case considered in this paper pertains to the nonrelativistic case.Relativistic considerations would have to be included for the superdiffusion case, which would lead to the Klein-Gordon equation in the integer order limit.In adapting the Feynman method, it is shown that it is preferable to introduce the action  as a fractional time integral of the Lagrangian and that it is necessary to introduce a "fractional Planck constant." In the limit of integer order, the regular action , the regular Schrodinger equation, and the regular Planck constant are all recovered.It may be of interest to note that there is a fractional Planck constant implied in Naber's work also, although it is not explicitly so stated.His equations are rendered nondimensional by using the Planck units of mass, length, and time and then generalized to fractional derivatives after including a Wick rotation.This implies a change of variable of time  →  so that the imaginary number is raised to the same power as the order of the fractional time derivative involved.However, the Planck units may not be the appropriate quantities as the equations involve quantities of fractional dimension and the equations must be made nondimensional after the generalization to fractional derivatives and not before.A fractional Planck constant does show up in Naber's treatment as well, as the ratio of masses /  [31,32].
A number of errors in Naber's work have been corrected.The correct continuity equation for the probability current is derived and the Green function solution for a free particle is given.The Green function is given in terms of a special type of function, the M-Wright function, which is used extensively in studies of time fractional subdiffusion studies.In the context of time fractional diffusion, the M-Wright function is a probability density function in time, which is non-Markovian and goes over to the Gaussian in the nonfractional limiting case.In the context of TFSE, the M-Wright function gives the propagator, which is a probability amplitude.Probability considerations are accounted for by the usual process of squaring of the amplitude.In particular, it is shown that Naber's result that the total probability is greater than unity in the long-time limit is a spurious result arising out of the operation of arbitrarily raising of the imaginary number  to the power of the same degree as the fractional time derivative invoking a Wick rotation.It is shown that such arbitrary change of the imaginary number cannot be carried out as the imaginary number  is connected with the phase of the action  in the path integral contribution.However, if we desired to consider a Wick rotation, it should be included as an additional increase of the index of the power of the imaginary number.The TFSE for a particle moving in a potential field is also derived.Furthermore, it is suggested that even in studies of fractional classical mechanics, such as those using variational methods, the action integral be expressed as a time fractional integral of the Lagrangian.Further extensions including the solutions to particle subject to different potentials are underway.
so that the fractional quantity   2 /2  has the dimensions of energy [ 2 / 2 ].Thus the dimension of the parameter "  " is [ 2− ] and the fractional momentum has the dimensions [/ 2−/2 ].The Lagrangian in (4) when generalized has the dimensions of energy.Of course, all quantities regain the standard dimensions in the limit  → 2.