# Time Fractional Schrodinger Equation Revisited

Authors: B. N. Narahari Achar; Bradley T. Yale; John W. Hanneken
Journal: Advances in Mathematical Physics
 (2013)

---

## Abstract

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[ML2/T2]. 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.


---

## Body

## 1. Introduction

There has been an explosive research output in recent years in the application of methods of fractional calculus [1–13] to the study of quantum phenomena [14–42]. The well-known 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–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–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–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–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 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 i 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–23, 28], including its generalization into space-time 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–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 i to a fractional power. The correct continuity equation for the probability density is also derived. The paper concludes with some new results.


## 2. Feynman Path Integral Method

The starting point for the Feynman method [43–46] is the classical Lagrangian L=L(x,x˙,t) and the action S=∫L(x,x˙,t)dt 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:
(1)mdxdt=p,dpdt=F.
Integrating with respect to time yields
(2)x=x0+1m∫0tp(τ)dτ,(3)p=p0+∫0tF(τ)dτ.
In the usual notation the Lagrangian is given by
(4)L(x,x˙,t)=T-V
and the action is given by
(5)S=∫0tL(x,x˙,τ)dτ.
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) K(xb,tb;xa,ta) for a particle starting from a position xa at time ta to reach a position xb at a later time tb, which arises from the contributions from all trajectories from xa to xb:
(6)K(xb,tb;xa,ta)=∑all pathsϕ[x(t)],
where the contribution from each of the paths has the form
(7)ϕ[x(t)]=const.exp[iSℏ].
Here S 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
(8)K(xb,tb;xa,ta)=∫abexp[iℏ∫tatbLdt]𝔇x(t).
The symbol 𝔇 indicates the fact that the operation of integration is carried over all paths from a to b.The wave functionψ(xb,tb) gives the total probability amplitude to arrive at xb at tb satisfying (9), where the integral is taken over all possible values of xa(9)ψ(xb,tb)=∫-∞∞K(xb,tb;xa,ta)ψ(xa,ta)dxa.
The kernel K can be computed by first carrying out a “time slicing” operation by dividing the time interval from ta to tb into N segments of duration
(10)ε=tb-taN,
where
(11)ta=t0<t1<t2<⋯<tN-1<tN=tb;(12)K(xb,tb;xa,ta)=limε→01A∭⋯∫exp[iS[b,a]ℏ]=limε→01A×dx1Adx2Adx3A⋯dxN-1A,
where
(13)S[b,a]=∫tatbL(x,x˙,t)dt
is a line integral taken over the trajectory passing through the point x(t). The constant A 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 ta to a time tb, when tb differs from ta 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  (i/ℏ) times the classical action for the infinitesimal time interval ε=tb-ta. With an obvious change of notation  xb=x, xa=x0, ta=t, tb=t+ε and using the fact that the particle is free, (2) yields
(14)p0=m(x-x0)ε
and (4) yields
(15)L=T=m(x-x0)22ε2
and (5) yields for the action
(16)S=εL=m(x-x0)22ε.
As a consequence, (9) becomes
(17)ψ(x,t+ε)=∫-∞∞1Aexp[im(x-x0)22ε]ψ(x0,t)dx0.
If x differs appreciably from x0, the exponential in (17) oscillates very rapidly and the integral over x0 contributes a very small value and only those paths which are very close to x give significant contributions. Changing the variable in the integral from x0 to η=x-x0 makes it ψ(x0,t)=ψ(x+η,t). Since both ε and η are small quantities, ψ(x,t+ε) may be expanded in Taylor’s series and only up to terms of order ε are retained. On the right-hand side, ψ(x+η,t) 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
(18)ψ(x,t)+ε∂ψ∂t=∫-∞∞1Aexp[-mη22iℏε]×(ψ(x,t)+η∂ψ∂x+η22∂2ψ∂x2)dη.
On the right-hand side the middle term vanishes on integration. It follows by equating the leading terms on both sides
(19)ψ(x,t)=∫-∞∞1Aexp[-mη22iℏε]ψ(x,t)dη
Hence(20a)A=∫-∞∞exp[-mη22iℏε]dη=2πiℏεm,(20b)∫-∞∞1Aexp[-mη22iℏε](η22∂2ψ∂x2)dη=εiℏ2m∂2ψ∂x2.Equating the remaining terms results in
(21)∂ψ∂t=iℏ2m∂2ψ∂x2.
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 termV(x,t) in the Lagrangian L=T-V in (15). This will necessitate [43] incorporating an additional factor {1-(iε/ℏ)V(x,t)} in (18), which becomes
(22)ψ(x,t)+ε∂ψ∂t=∫-∞∞1Aexp[-mη22iℏε]{1-iεℏV(x,t)}×(ψ(x,t)+η∂ψ∂x+η22∂2ψ∂x2)dη.
Then (21) becomes
(23)∂ψ∂t=iℏ2m∂2ψ∂x2-iℏVψ.
Multiplying both sides by -ℏ/i results in the standard Schrodinger equation of quantum mechanics:
(24)iℏ∂ψ∂t=-ℏ22m∂2ψ∂x2+Vψ.
It is to be noted that the imaginary number i on the left-hand side is not arbitrary; it arises from the coefficient of the potential term V in (23) but ultimately from the coefficient of S/ℏ 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.


## 3. 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]:
(25)0CDtβf(t)=1Γ(n-β)∫0tfn(τ)dτ(t-τ)β+1-n(n-1<β<n).
In the limit β→n, the Caputo derivative becomes the ordinary nth derivative of the function.The fractional integral of orderβ is defined by
(26)0Itβf(t)=1Γ(β)∫0tf(τ)(t-τ)β-1dτ.
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
(27)x=x0+1mfΓ(α/2)∫0tpf(τ)(t-τ)α/2-1dτ,pf=pf0+1Γ(α/2)∫0tF(τ)(t-τ)α/2-1dτ,
where as usual x0, pf0 refer to the initial position and initial value of the “fractional momentum” pf, respectively. It is to be noted that the variables x, t are still the space and time variables and have the dimensions of length and time [L] and [T], respectively. However, the dynamical quantities obtained by the operation of fractional derivation have different dimensions; for example, “fractional velocity” with the notation 0CDtα/2x=x˙α/2 would have the units [L/Tα/2]. The dimension for the parameter “mf” in the fractional momentum, pf=mfx˙α/2, is no longer just [M] but has to be chosen [47] so that the fractional quantity pf2/2mf has the dimensions of energy [ML2/T2]. Thus the dimension of the parameter “mf” is [MT2-α] and the fractional momentum has the dimensions [ML/T2-α/2]. The Lagrangian in (4) when generalized has the dimensions of energy. Of course, all quantities regain the standard dimensions in the limit α→2.


## 4. Time Fractional Schrodinger Equation for a Free Particle

There are two possible generalizations of the action integral in (5) used in fractional dynamics [11]:
(28)SI=∫0tL(x,x˙α/2,t)dt,(29)SII=1Γ(α/2)∫0tL(x,x˙α/2,τ)(t-τ)α/2-1dτ.
Since the Lagrangian has the dimensions of energy, the dimensions of action defined in (28) and (29), SI and SII, are different.The dimensions ofSI are the same as that of the regular action, namely, [ML2/T], but that of SII is [ML2/T2-α/2]. In the Newtonian limit α→2, SII→ 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” ℏf with appropriate dimensions must be introduced in order to render the argument of the exponential in (7) dimensionless.For a “free particle” (27) yields
(30)x=x0+pf0tα/2mfΓ(1+α/2),(31)pf=pf0.
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
(32)pf0=(x-x0)mfΓ(1+α/2)εα/2
and hence
(33)L=mfΓ(1+α/2)(x-x0)22εα.
But for action, there are two choices:
(34)SI=mfΓ2(1+α/2)(x-x0)22εα-1,(35)SII=mfΓ(1+α/2)(x-x0)22εα/2.Making the appropriate changes, the equations for the evolution of the wave function in the two cases are(36)ψI(x,t+ε)=∫-∞∞1AIexp[imfΓ2(1+α/2)(x-x0)22ℏεα-1]×ψI(x0,t)dx0,(37)ψII(x,t+ε)=∫-∞∞1AIIexp[imfΓ(1+α/2)(x-x0)22ℏfεα/2]×ψII(x0,t)dx0.
Changing the variable in the integrals from x0 to η=x-x0 as before and introducing two constants
(38)aI=mfΓ2(1+α/2)2iℏ,aII=mfΓ(1+α/2)2iℏf.
Equations (36) and (37) can be written as
(39)ψI(x,t+ε)=∫-∞∞1AIexp[-aIη2εα-1]ψI(x+η,t)dη,(40)ψII(x,t+ε)=∫-∞∞1AIIexp[-aIIη2εα/2]ψII(x+η,t)dη.
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
(41)ψI,II(x,t+ε)=ψI,II(x,t)+0CDtγψI,II(x,t)εγΓ(γ+1)+⋯.
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(x,t) becomes
(42)ψI(x,t)+0CDtγψI(x,t)εγΓ(γ+1)=∫-∞∞1AIexp[-aIη2εα-1]×{ψI(x,t)+η∂ψI∂x+η22∂2ψI∂x2}dη.
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) yieldsψ
I
(
x
,
t
)
=
(
1
/
A
I
)
π
ε
α
-
1
/
a
I
ψ
I
(
x
,
t
), requiring that AI=πεα-1/aI. Equation (42) reduces to
(43)ψI(x,t)+0CDtγψI(x,t)εγΓ(γ+1)={ψI(x,t)+14εα-1aI∂2ψI(x,t)∂x2}.
Equating the remaining terms requires that the powers of ε must be the same on both sides; that is, γ=α-1. Inserting the value of aI and simplifying (43) yield
(44)0CDtα-1ψI(x,t)=iℏ2mfΓ(α)Γ2(1+α/2)∂2ψI(x,t)∂x2.
Similarly, by expanding (40) the equation for ψII(x,t) becomes
(45)ψII(x,t)+0CDtγψII(x,t)εγΓ(γ+1)=∫-∞∞1AIIexp[-aIIη2εα/2]×{ψII(x,t)+η∂ψII∂x+η22∂2ψII∂x2}dη.
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 AII=πεα/2/aII. Equation (45) reduces to
(46)ψII(x,t)+0CDtγψII(x,t)εγΓ(γ+1)={ψII(x,t)+14εα/2aII∂2ψII(x,t)∂x2}.
In the remaining terms, the powers of ε must be the same. This requires γ=α/2. Inserting the value of  aII and simplifying (46) yield the equation for the wave function ψII(x,t):
(47)0CDtα/2ψII(x,t)=iℏf2mf∂2ψII(x,t)∂x2.
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
(48)∂ψ(x,t)∂t=iℏ2m∂2ψ(x,t)∂x2
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).The coefficient of the space derivative term in (47) has the dimension [L2/Tα/2], corresponding to the fractional diffusion coefficient. Thus (47) can be considered a time fractional diffusion equation with an imaginary fractional diffusion coefficient, just as (48), the regular Schrodinger equation, can be considered to be a diffusion equation with an imaginary diffusion coefficient. Thus all the mathematical machinery of time fractional diffusion theory [5, 50–59] can be imported advantageously.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, ∂tβψ, with β=α/2 will be used. Naturally, 0<β≤1, and β→1 yields the regular first-order time derivative, denoted by ∂t. After incorporating these changes and defining a new constant Df=ℏf/2mf, (47) becomes
(49)∂tβψ(x,t)=iDf∂2ψ(x,t)∂x2(0<β≤1).Equation (49) can be solved by a combination of the Fourier and Laplace transform methods [52–54].


## 5. Probability Current and the Continuity Equation

The probability densityρ is defined by ρ=ψ*ψ=ψψ*. The complex conjugate wave function satisfies
(50)∂tβψ*(x,t)=-iDf∂2ψ*(x,t)∂x2.
There is an identity satisfied by the Caputo derivative [5]
(51)∂t1-β∂tβf(t)=∂tf(t),
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
(52)∂tρ=∂t(ψψ*)=(∂tψ)ψ*+ψ(∂tψ*).
Inserting from (51) gives
(53)∂tρ=(∂t1-β∂tβψ)ψ*+ψ(∂t1-β∂tβψ*).
Substituting from (49) and (50), (53) yields
(54)∂tρ=(∂t1-β(iDf∂2ψ(x,t)∂x2))ψ*+ψ(∂t1-β(-iDf∂2ψ*(x,t)∂x2)).
Equation (54) can be rewritten after factoring out the constant and interchanging the order of space and time derivatives as
(55)∂tρ=iDf{(∂2∂x2∂t1-βψ)ψ*-ψ(∂2∂x2∂t1-βψ*)}.
Introducing a new function ψ~=∂t1-βψ, (55) can be written as
(56)∂tρ=iDf{(∂2∂x2ψ~)ψ*-ψ(∂2∂x2ψ~*)}.
Defining a probability current density given by
(57)Jx=-iDf[ψ~*∂ψ∂x-ψ∂ψ~*∂x],
equation (56) can be written as
(58)∂tρ+∂Jx∂x=0.
This is the time fractional version of the continuity equation. In the limit β→1, ψ~→ψ  ψ~*→ψ* and (58) reproduce the continuity equation of standard quantum mechanics [49]. It may be noted that (58) differs from Naber’s result ((24) in [19]) and will be discussed later.


## 6. Solution for TFSE for a Free Particle

The solution for the TFSE for a free particle under the conditionsψ
(
x
,
0
)
=
ψ
0
(
x
); ψ(x,t)→0, |x|→∞, t>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(59)ψ~^(k,s)=12π∫-∞∞e-ikx[∫0∞e-stψ(x,t)dt]dx
equation (49) reduces to
(60)sβψ~^(k,s)-sβ-1ψ~(k,0)=-iDfk2ψ~^(k,s)
resulting in
(61)ψ~^(k,s)=ψ~(k,0)sβ-1sβ+iDfk2.Applying the inverse Laplace transform, (61) yields
(62)ψ~(k,t)=ψ~(k,0)Eβ,1(-iDfk2tβ)
in terms of the Mittag-Leffler function defined by a series or by the inverse Lapalce transform [3]:
(63)Eα,β(z)=∑n=0∞znΓ(αn+β)=L-1{sα-βsα-z}.The Green function solution can be written as(64)ψ(x,t)=∫-∞∞ψ(x-ξ)Gβ(ξ,t)dξ,
where the Green function is given by the inverse Fourier transform of the Mittag-Leffler function
(65)Gβ(x,t)=12π∫-∞∞eikxEβ,1(-iDfk2tβ)dk,
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
(66)Mβ(z)=W(-z;-β,1-β)=∑n=0∞(-z)nn!Γ(-βn+1-β)0<β<1,
where the Wright function is defined by [54]
(67)W(z;λ,μ)=∑n=0∞znn!Γ(λn+μ)λ>-1,μ∈C.
The Green function in (65) is then given by
(68)Gβ(x,t)=121iDftβ/2Mβ/2(|x|iDftβ/2).
In the context of fractional diffusion, the function Mβ/2(z) belongs to the Wright type of probability densities characterized by the similarity variable z=|x|/Dftβ/2, where Df is the fractional diffusion coefficient, with ∫0∞Mβ/2(z)dz=1. Furthermore, the probability densities are non-Markovian and exhibit a variance consistent with slow anomalous diffusion [49–51], σβ2(t)=(2/Γ(β+1))Dftβ.In the limit ofβ→1, the probability density function goes over to the Gaussian
(69)M1/2(z)=1πe-z2/4
corresponding to regular diffusion. A plot of the reduced probability density function is given in Figure 1.Figure 1
Reduced probability density function.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 fromxa at ta to xb at tb. 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 t→∞; the case discussed by Naber [19] and will be discussed in the next section.


## 7. 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):
(70)Dt1-νψ(t,x)=1Γ(1-ν)×∫0tdψ(τ,x)dτdτ(t-τ)ν(0<ν<1).
This is incorrect, as can be checked easily by taking the limit ν→1. The left-hand side →ψ(t,x) 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
(71)Dt1-νψ(t,x)=1Γ(ν)∫0tdψ(τ,x)dτdτ(t-τ)1-ν.
As a consequence, the weight factor in (N18) should be (t-τ)1-ν and not (t-τ)-ν, 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).
(72)Dt1-νDtνy(t)=dydt-[Dtνy(t)]t=0t1-νΓ(ν).
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(73)Dt-νDtνy(t)=y(t)-y(0).
This yields  Dt1-νDtνy(t)=dy/dt=∂ty(t)  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(74)Ψ=Ψ0Eν(ω(-it)ν),
where Eν(z) in Naber’s notation corresponds to the Mittag-Leffler function Eν,1(z) 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(75)A(t)=12πi∫γ-i∞γ+i∞estsν-1A0dssν-σiν
along a Hankel contour and considering the contribution from the residue of the pole s0=σ1/νi 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
(76)Ψ=Ψ0{e-iω1/νtν-Fν(ω(-i)ν,t)}.
He argues that in the limit of t→∞ 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 t→∞.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):
(77)ψ~(k,t)=12πi∫γ-i∞γ+i∞estsβ-1ψ~(k,0)dssβ+iDfk2.
As usual, the Bromwich contour is replaced by the Hankel contour. The two contributions arise from (i) the residue at the pole at s0=(Dfk2)1/β(-i)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
(78)Residueψ~(k,0)=e(-i)1/β(Dfk2)1/βtβ
and corresponds to the first term on the right-hand side in (76). However,
(79)(-i)1/β=(e-iπ/2)1/β=(e-iπ/2β)=(cosπ2β-isinπ2β).
Therefore, (78) yields
(80)ψ~(k,t)ψ~(k,0)=Residueψ~(k,0)=et(Dfk2)1/βcos(π/2β)e-it(Dfk2)1/βsin(π/2β).
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 t→∞. Thus the contribution from the residue due to the pole also →0. Thus the entire wave function →0 in the limit of t→∞, 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 ats0=σ1/νi, whereas in the present paper, the pole occurs (using the same notation as Naber) at s0=σ1/ν(-i)1/ν. The simple i 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 t→∞, namely, zero total probability. Naber’s result is a direct consequence of his choice to raise the power of the imaginary number i 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 iS/ℏ in the Feynman propagator. If it is necessary to include a Wick rotation, the power of i 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 t→∞. 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.


## 8. 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 termV(x,t) in the Lagrangian L=T-V. This will necessitate incorporating an additional term -V(x,t)εα/2/Γ(1+α/2) in the right-hand side of (35) and an additional factor -(i/hf)(V(x,t)εα/2/Γ(1+α/2)) in the exponential in (37). This results in changing (45) into
(81)ψII(x,t)+0CDtβψII(x,t)εβΓ(β+1)=∫-∞∞1AIIexp[-aIIη2εα/2]{1-ihfV(x,t)εα/2Γ(1+α/2)}×{ψII(x,t)+η∂ψII∂x+η22∂2ψII∂x2}dη.Equation (47) then becomes
(82)0CDtα/2ψII(x,t)=iℏf2mf∂2ψII(x,t)∂x2-ihfV(x,t)ψ(x,t).
Multiplying both sides by -ℏf/i results in the FTSE for a particle in a potential field as
(83)iℏf∂tβψ(x,t)=-ℏf22mf∂2ψ(x,t)∂x2-V(x,t)ψ(x,t)
which reduces to the standard Schrodinger equation of quantum mechanics:
(84)iℏ∂ψ∂t=-ℏ22m∂2ψ∂x2+Vψ.
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).


## 9. 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 actionS 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 S, 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 t→it 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 m/Mp [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 numberi 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 i is connected with the phase of the action S 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.



---

*Source: 290216-2013-07-28.xml*