- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Advances in Mathematical Physics
Volume 2013 (2013), Article ID 290216, 11 pages
Time Fractional Schrodinger Equation Revisited
Physics Department, University of Memphis, Memphis, TN 38152, USA
Received 29 April 2013; Revised 1 July 2013; Accepted 2 July 2013
Academic Editor: Dumitru Baleanu
Copyright © 2013 B. N. Narahari Achar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
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 . The action 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.
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 , 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  and by using a fractional variation principle and a fractional Klein-Gordon equation . 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 . 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 .
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 . This theory has been further generalized recently within the frame work of tempered ultradistributions . 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  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  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 , 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 . 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 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 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 . 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 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 but also in the limit of small time interval . The Kernel is proportional to the exponential of times the classical action for the infinitesimal time interval . With an obvious change of notation , , , 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 , the exponential in (17) oscillates very rapidly and the integral over contributes a very small value and only those paths which are very close to give significant contributions. Changing the variable in the integral from to makes it . 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 HenceEquating 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). This will necessitate  incorporating an additional factor 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 left-hand 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.
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 . The Caputo derivative of order is defined by : 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 . Then (2) and (3) become where as usual , 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 would have the units . The dimension for the parameter “” in the fractional momentum, , is no longer just but has to be chosen  so that the fractional quantity has the dimensions of energy . Thus the dimension of the parameter “” is and the fractional momentum has the dimensions . The Lagrangian in (4) when generalized has the dimensions of energy. Of course, all quantities regain the standard dimensions in the limit .
4. Time Fractional Schrodinger Equation for a Free Particle
There are two possible generalizations of the action integral in (5) used in fractional dynamics : Since the Lagrangian has the dimensions of energy, the dimensions of action defined in (28) and (29), and , are different.
The dimensions of are the same as that of the regular action, namely, , but that of is . In the Newtonian limit , 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 to 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  in time, with fractional derivative of order , and keeping only the lowest-order term in yields 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 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) yields
, requiring that . Equation (42) reduces to Equating the remaining terms requires that the powers of must be the same on both sides; that is, . Inserting the value of and simplifying (43) yield Similarly, by expanding (40) the equation for 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 . Equation (45) reduces to In the remaining terms, the powers of must be the same. This requires . Inserting the value of and simplifying (46) yield the equation for the wave function : 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 , 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 .
It should be noted that since , with being the limiting case, the order of the time fractional derivative is 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 . Furthermore, the order 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 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 , 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, , with will be used. Naturally, , and yields the regular first-order time derivative, denoted by . After incorporating these changes and defining a new constant , (47) becomes
5. 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  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 , (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. In the limit , and (58) reproduce the continuity equation of standard quantum mechanics . It may be noted that (58) differs from Naber’s result (() in ) 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
; , , 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
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  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  The Green function in (65) is then given by In the context of fractional diffusion, the function belongs to the Wright type of probability densities characterized by the similarity variable , where is the fractional diffusion coefficient, with . Furthermore, the probability densities are non-Markovian and exhibit a variance consistent with slow anomalous diffusion [49–51], .
In the limit of , 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  and will be discussed in the next section.
7. Comments about Some Results in Naber’s Work 
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 , 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 (), reproduced here for convenience in (70): This is incorrect, as can be checked easily by taking the limit . 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 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 .
(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  and reproduced later in the current notation for This yields 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 .
(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 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 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 , assuming that the wave function was initially normalized. Since , 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 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 and corresponds to the first term on the right-hand side in (76). However, Therefore, (78) yields The right-hand side of (80) has an amplitude term and an oscillatory factor. Because , 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 , whereas in the present paper, the pole occurs (using the same notation as Naber) at . 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 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.
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 term in the Lagrangian . This will necessitate incorporating an additional term in the right-hand side of (35) and an additional factor in the exponential in (37). This results in changing (45) into
Equation (47) then becomes Multiplying both sides by results in the FTSE for a particle in a potential field as which reduces to the standard Schrodinger equation of quantum mechanics: 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 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.
- K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.
- K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 1993.
- I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
- S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives: Theory and Applications, Gordon and Breach Science Publishers, Yverdon, Switzerland, 1993.
- A. Carpinteri and F. Mainardi, Fractals and Fractional Calculus in Continuum Mechanics, vol. 378 of CISM Courses and Lectures, Springer, Vienna, Austria, 1997.
- R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific Publishing, River Edge, NJ, USA, 2000.
- B. J. West, M. Bologna, and P. Grigolini, Physics of Fractal Operators, Institute for Nonlinear Science, Springer, New York, NY, USA, 2003.
- R. L. Magin, Fractional Calculus in Bioengineering, Begell House, Redding, Conn, USA, 2006.
- L. Debnath and D. Bhatta, Integral Transforms and Their Applications, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2nd edition, 2007.
- R. Hermann, Fractional Calculus, World Scientific, Singapore, 2011.
- A. B. Malinowska and D. F. M. Torres, Introduction to the Fractional Calculus of Variations, Imperial College Press, London, UK, 2012.
- D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo, Fractional Calculus: Models and Numerical Methods, vol. 3 of Series on Complexity, Nonlinearity and Chaos, World Scientific Publishing, Hackensack, NJ, USA, 2012.
- A. K. GolmanKhaneh, Investigations in Dynamics:With Focus on Fractal Dynamics, Lap Lambert Academic Publishing, 2012.
- N. Laskin, “Fractional quantum mechanics,” Physical Review E, vol. 62, no. 3, pp. 3135–3145, 2000.
- N. Laskin, “Fractional quantum mechanics and Lévy path integrals,” Physics Letters A, vol. 268, no. 4–6, pp. 298–305, 2000.
- N. Laskin, “Fractional Schrödinger equation,” Physical Review E, vol. 66, no. 5, Article ID 056108, 7 pages, 2002.
- N. Laskin, “Fractals and quantum mechanics,” Chaos, vol. 10, no. 4, pp. 780–790, 2000.
- N. Laskin, “Principles of fractional quantum mechanics,” http://arxiv.org/abs/1009.5533.
- M. Naber, “Time fractional Schrödinger equation,” Journal of Mathematical Physics, vol. 45, no. 8, pp. 3339–3352, 2004.
- R. K. Saxena, R. Saxena, and S. L. Kalla, “Solution of space-time fractional Schrödinger equation occurring in quantum mechanics,” Fractional Calculus & Applied Analysis, vol. 13, no. 2, pp. 177–190, 2010.
- S. Wang and M. Xu, “Generalized fractional Schrödinger equation with space-time fractional derivatives,” Journal of Mathematical Physics, vol. 48, no. 4, Article ID 043502, 10 pages, 2007.
- J. Dong and M. Xu, “Space-time fractional Schrödinger equation with time-independent potentials,” Journal of Mathematical Analysis and Applications, vol. 344, no. 2, pp. 1005–1017, 2008.
- X. Guo and M. Xu, “Some physical applications of fractional Schrödinger equation,” Journal of Mathematical Physics, vol. 47, no. 8, Article ID 082104, 9 pages, 2006.
- X. Y. Jiang, “Time-space fractional Schrödinger like equation with a nonlocal term,” European Physical Journal: Special Topics, vol. 193, no. 1, pp. 61–70, 2011.
- J. Dong, “Fractional Green's function for the time-dependent scattering problem in the space-time-fractional quantum mechanics,” http://arxiv.org/abs/1301.3206.
- J. Dong and M. Xu, “Applications of continuity and discontinuity of a fractional derivative of the wave functions to fractional quantum mechanics,” Journal of Mathematical Physics, vol. 49, no. 5, Article ID 052105, 16 pages, 2008.
- E. Scalas, D. Baleanu, F. Mainardi, and A. Mura, “Fractional calculus and the Schrodinger equation,” in Proceedings of the 2nd IFAC Workshop on Fractional Differentiation and Its Applications, Porto, Portugal, July 2006.
- M. Bhatti, “Fractional Schrödinger wave equation and fractional uncertainty principle,” International Journal of Contemporary Mathematical Sciences, vol. 2, no. 17–20, pp. 943–950, 2007.
- S. Wang, M. Xu, and X. Li, “Green's function of time fractional diffusion equation and its applications in fractional quantum mechanics,” Nonlinear Analysis: Real World Applications, vol. 10, no. 2, pp. 1081–1086, 2009.
- A. Tofighi, “Probability structure of time fractional schrodinger equation,” Acta Physica Polonica A, vol. 116, no. 2, pp. 114–118, 2009.
- A. Iomin, “On fractional time quantum dynamics,” Physical Review E, vol. 80, no. 2, Article ID 022103, 4 pages, 2009.
- A. Iomin, “Fractional-time Schrödinger equation: fractional dynamics on a comb,” Chaos, Solitons & Fractals, vol. 44, no. 4-5, pp. 348–352, 2011.
- A. M. Shahin, E. Ahmed, and Y. A. Omar, “On fractional order quantum mechanics,” International Journal of Nonlinear Science, vol. 8, no. 4, pp. 469–472, 2009.
- P. Rozmej and B. Bandrowski, “On fractional Schrödinger equation,” Computational Methods in Science and Technology, vol. 16, pp. 191–194, 2010.
- H. Ertik, D. Demirhan, H. Şirin, and F. Büyükkılıç, “Time fractional development of quantum systems,” Journal of Mathematical Physics, vol. 51, no. 8, Article ID 082102, 15 pages, 2010.
- S. I. Muslih, O. P. Agrawal, and D. Baleanu, “A fractional Schrödinger equation and its solution,” International Journal of Theoretical Physics, vol. 49, no. 8, pp. 1746–1752, 2010.
- V. B. L. Chaurasia and D. Kumar, “Solutions of Unified fractional Schrödinger equation,” ISRN Mathematical Physics, vol. 2012, Article ID 935365, 7 pages, 2012.
- A. K. Mahata, “On fractional Schrödinger equation and its application,” International Journal of Scientific Engineering Research, vol. 4, pp. 1–5, 2013.
- A. L. de Paoli and M. C. Rocca, “The fractionary Schrödinger equation, Green functions and ultradistributions,” Physica A, vol. 392, no. 1, pp. 111–122, 2013.
- E. C. Oliveira, F. S. Costa, and J. Vaz Jr., “The fractional Schrödinger equation for delta potentials,” Journal of Mathematical Physics, vol. 51, no. 12, Article ID 123517, 16 pages, 2010.
- K. M. Kolwankar and A. D. Gangal, “Fractional differentiability of nowhere differentiable functions and dimensions,” Chaos, vol. 6, no. 4, pp. 505–513, 1996.
- X. J. Yang, D. Baleanu, and J. A. T. Machado, “Mathematical aspects of the Heisenberg uncertainty principle within local fractional Fourier analysis,” Boundary Value Problems, vol. 2013, article 131, 2013.
- R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, New York, NY, USA, 1965.
- H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets, World Scientific Publishing, Singapore, 2006.
- L. S. Schulman, Techniques and Applications of Path Integration, John Wiley & Sons, New York, NY, USA, 1981.
- M. Kac, Probability and Related Topics in Physical Sciences, Interscience, London, UK, 1957.
- B. N. Narahari Achar, J. W. Hanneken, N. Becker, and M. Morgan, “Intrinsic dissipation in fractional order dynamics,” in Proceedings of the International Workshop on Fractional Differentiation and Its Applications (FDA '10), Badajoz, Spain, October 2010.
- B. N. Narahari Achar and J. W. Hanneken, “Fractional order dynamics under gravity,” in Proceedings of the International Workshop on Fractional Differentiation and Its Applications (FDA '10), Badajoz, Spain, October 2010.
- L. I. Schiff, Quantum Mechanics, McGraw-Hill, New York, NY, USA, 1968.
- F. Mainardi, P. Paradisi, and R. Gorenflo, “Probability distributions generated by fractional order diffusion equations,” in Econophysics: An Emerging Science, J. Kertesz and I. Kondov, Eds., Kluwer Academic Publishers, Dordrecht, The Netherlands, 1998.
- R. Gorenflo, A. Iskenderov, and Y. Luchko, “Mapping between solutions of fractional diffusion-wave equations,” Fractional Calculus & Applied Analysis, vol. 3, no. 1, pp. 75–86, 2000.
- F. Mainardi, Y. Luchko, and G. Pagnini, “The fundamental solution of the space-time fractional diffusion equation,” Fractional Calculus & Applied Analysis, vol. 4, no. 2, pp. 153–192, 2001.
- F. Mainardi, “The time fractional diffusion-wave equation,” Radiofisica, vol. 38, no. 1-2, pp. 20–36, 1996.
- F. Mainardi, “The fundamental solutions for the fractional diffusion-wave equation,” Applied Mathematics Letters, vol. 9, no. 6, pp. 23–28, 1996.
- F. Mainardi and G. Pagnini, “The role of the Fox-Wright functions in fractional sub-diffusion of distributed order,” Journal of Computational and Applied Mathematics, vol. 207, no. 2, pp. 245–257, 2007.
- F. Mainardi, A. Mura, and G. Pagnini, “The M-Wright function in time-fractional diffusion processes: a tutorial survey,” International Journal of Differential Equations, vol. 2010, Article ID 104505, 29 pages, 2010.
- R. Gorenflo, Y. Luchko, and F. Mainardi, “Analytical properties and applications of the Wright function,” Fractional Calculus & Applied Analysis, vol. 2, no. 4, pp. 383–414, 1999.
- R. Gorenflo, F. Mainardi, D. Moretti, G. Pagnini, and P. Paradisi, “Fractional diffusion: probability distributions and random walk models,” Physica A, vol. 305, no. 1-2, pp. 106–112, 2002.
- F. Mainardi, Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models, Imperial College Press, London, UK, 2010.