We obtain the solution of a unified fractional Schrödinger equation. The solution is derived by the application of the Laplace and Fourier transforms in closed form in terms of the Mittag-Leffler function. The result obtained here is quite general in nature and capable of yielding a very large number of results (new and known) hitherto scattered in the literature. Most of results obtained are in a form suitable for numerical computation.
1. Introduction
Recent applications of fractional differential equations to number of systems such as those exhibiting enormously slow diffusion or subdiffusion have given opportunity for physicists to study even more complicated systems. These systems include charge transport in amorphous semiconductors: the relaxation in polymer systems, fluid mechanics, viscoelasticity, and Hall effect. The generalized diffusion equation is studied to describe complex systems with anomalous behavior in much the same way as simpler systems. Fractional calculus is now considered as practical techniques in many branches of applied sciences and engineering. Several authors notably Hilfer [1], Beyer and Kempfle [2], Kempfle and Gaul [3], Schneider and Wyss [4], and Debnath [5–7] have discussed many examples of homogeneous fractional ordinary differential equations and homogeneous fractional diffusion and wave equations.
Laskin [8–10] constructed space fractional quantum mechanics by using Feynman path integrals, the only difference being that Lévy distributions are employed instead of Gaussians distributions for the set of possible paths. The Schrödinger equation thus obtained contains fractional derivatives. Naber [11] has investigated certain properties of time fractional Schrödinger equation by writing the Schrödinger equation in terms of fractional derivatives as dimensionless objects. In recent work, solutions of fractional Schrödinger equations are investigated by Bhatti [12], Chaurasia and Singh [13], Saxena et al. [14], among others.
In the present paper, we obtain the solution of a unified fractional Schrödinger equation. The result obtained here provides an elegant extension of the results given earlier by Bhatti [12], Chaurasia and Singh [13], Saxena et al. [14], and Debnath [5].
The Riemann-Liouville fractional integral of order ν is defined by (Miller and Ross [15, page 45]; Kilbas et al. [16]),Dt-ν0f(x,t)=1Γ(ν)∫0t(t-u)ν-1f(x,u)du,
where Re(ν)>0.
The following fractional derivative of order α>0 is introduced by Caputo [17]; see also Kilbas et al. [16] in the formDtα0f(x,t)=1Γ(m-α)∫0tf(m)(x,τ)dτ(t-τ)α+1-m,m-1<α≤m,Re(α)>0,m∈N,=∂mf(x,t)∂tm,ifα=m,
where ∂mf(x,t)/∂tm is the mth partial derivative of f(x,t) with respect to t.
The Laplace transform of the Caputo derivative is given by Caputo [17]; see also Kilbas et al. [16] in the formL{Dtα0f(x,t)}=sαf̅(x,s)-∑r=0m-1sα-r-1f(r)(x,0+),(m-1<α≤m).
We also need the Weyl fractional operator, defined byDxμ-∞f(x,t)=1Γ(n-μ)(∂∂x)n∫-∞xf(τ,t)(x-τ)μ-n+1dτ,
where n=[μ] is an integral part of μ>0.
Its Fourier transform is (Metzler and Klafter [18, page 59, A.11])F{Dxμ-∞f(x,t)}=(ik)μf*(k,t),
where f*(k,t) is the Fourier transform of f(x,t) with respect to the variable x of f(x,t).
Following the convention initiated by Compte [19], we suppress the imaginary unit in Fourier space by adopting the slightly modified form of above result in our investigations (Metzler and Klafter, [18, p.59, A.12])F{Dxμ-∞f(x,t)}=-|k|μf*(k,t)
instead of (1.5).
2. Unified Fractional Schrödinger Equation
In this section, we will derive the solution of the unified fractional Schrödinger equation (2.1). The result is as follows.
Theorem 2.1.
Consider the following unified fractional Schrödinger equation:
Dtα0ψ(x,t)+aDtβ0ψ(x,t)+bDtγ0ψ(x,t)=iℏ2mDxμ-∞ψ(x,t),-∞<x<∞,t>0,0<α≤1,0<β≤1,0<γ≤1,μ>0,
with the initial conditions:
ψ(x,0)=ψ0(x),-∞<x<∞,ψ(x,t)⟶0as|x|⟶∞,t>0,
where ℏ is Planck’s constant divided by 2π, m is the mass, and ψ(x,t) is a wave function of the particle. Then, for the solution of (2.1), subject to the initial conditions (2.2), there holds the formula
ψ(x,t)=12π∫-∞∞ψ0*(k)∑r=0∞(-1)r∑l=0r(rl)albr-lt(α-γ)r+(γ-β)l×Eα,(α-γ)r+(γ-β)l+1r+1(-η|k|μtα)e-ikxdk+12π∫-∞∞aψ0*(k)∑r=0∞(-1)r∑l=0r(rl)albr-lt(α-γ)r+(γ-β)l+α-β×Eα,(α-γ)r+(γ-β)l+α-β+1r+1(-η|k|μtα)e-ikxdk+12π∫-∞∞bψ0*(k)∑r=0∞(-1)r∑l=0r(rl)albr-lt(α-γ)r+(γ-β)l+α-γ×Eα,(α-γ)r+(γ-β)l+α-γ+1r+1(-η|k|μtα)e-ikxdk,
where η=iℏ/2m and Eα,βγ(z) is the generalized Mittag-Leffler function [20].
Proof.
Applying the Laplace transform with respect to the time variable t on both the sides of (2.1) and using the initial conditions (2.2), we get
sαψ̅(x,s)-sα-1ψ0(x)+asβψ̅(x,s)-asβ-1ψ0(x)+bsγψ̅(x,s)-bsγ-1ψ0(x)=η-∞Dxμψ̅(x,s),
where η=iℏ/2m.
If we apply the Fourier transform with respect to variable x and use the formula (1.6), it yields
sαψ̅*(k,s)-sα-1ψ0*(k)+asβψ̅*(k,s)-asβ-1ψ0*(k)+bsγψ̅*(k,s)-bsγ-1ψ0*(k)=-η|k|μψ̅*(k,s).
Solving for ψ̅*(k,s), it gives
ψ̅*(k,s)=sα-1ψ0*(k)sα+asβ+bsγ+η|k|μ+asβ-1ψ0*(k)sα+asβ+bsγ+η|k|μ+bsγ-1ψ0*(k)sα+asβ+bsγ+η|k|μ.
On taking the inverse Laplace transform of (2.6) and applying the formula (Saxena et al. [21]), it is seen that
ψ*(k,t)=ψ0*(k)∑r=0∞(-1)r∑l=0r(rl)albr-lt(α-γ)r+(γ-β)l×Eα,(α-γ)r+(γ-β)l+1r+1(-η|k|μtα)+aψ0*(k)∑r=0∞(-1)r∑l=0r(rl)albr-lt(α-γ)r+(γ-β)l+α-β×Eα,(α-γ)r+(γ-β)l+α-β+1r+1(-η|k|μtα)+bψ0*(k)∑r=0∞(-1)r∑l=0r(rl)albr-lt(α-β)r+(γ-β)l+α-γ×Eα,(α-γ)r+(γ-β)l+α-γ+1r+1(-η|k|μtα).
Finally, the required solution (2.3) is obtained by taking inverse Fourier transform of (2.7).
3. Special Cases
If we take b=0 in (2.1), then we obtain the following result.
Corollary 3.1.
Consider the following fractional Schrödinger equation
Dtα0ψ(x,t)+aDtβ0ψ(x,t)=iℏ2mDxμ-∞ψ(x,t),-∞<x<∞,t>0,0<α≤1,0<β≤1,μ>0,
with the initial conditions:
ψ(x,0)=ψ0(x),ψ(x,t)⟶0,as|x|⟶∞,t>0,
where ℏ is Planck’s constant divided by 2π, m is the mass, and ψ(x,t) is the wave function of the particle. Then, for the solution of (3.1), subject to the initial conditions (3.2), there holds the formula
ψ(x,t)=12π∫-∞∞ψ0*(k)∑r=0∞(-b)rEα,(α-β)r+1r+1(-η|k|μtα)e-ikxdk+12π∫-∞∞aψ0*(k)∑r=0∞(-b)rEα,(α-β)r+α-β+1r+1(-η|K|μtα)e-ikxdk.
If we set a=0=b in (2.1), then we arrive at the following result recently obtained by Saxena et al. [14].
Corollary 3.2.
The solution of the following fractional Schrödinger equation:
Dtα0ψ(x,t)=iℏ2mDxμ-∞ψ(x,t),-∞<x<∞,t>0,0<α≤1,μ>0,
with the initial conditions
ψ(x,0)=ψ0(x),ψ(x,t)⟶0,as|x|⟶∞,t>0,
where ℏ is the Planck’s constant divided by 2π, m is the mass, and ψ(x,t) is a wave function of the particle, is given by
ψ(x,t)=12π∫-∞∞ψ0*(k)tαEα,1(-η|k|μtα)e-ikxdk.
Finally, on taking a=0=b and μ = 2 in (2.1), then we arrive at the following result given by Bhatti [12].
Corollary 3.3.
Consider the following fractional Schrödinger equation:
Dtα0ψ(x,t)=iℏ2m∇2ψ(x,t),-∞<x<∞,t>0,0<α≤1,
with the initial conditions:
ψ(x,0)=ψ0(x),ψ(x,t)⟶0,as|x|⟶∞,t>0,
where ℏ is Planck’s constant divided by 2π, m is the mass, and ψ(x,t) is a wave function of the particle. Then, for the solution of (3.7), under the initial conditions (3.8), there holds the relation:
ψ(x,t)=12π∫-∞∞ψ0*(k)tαEα,1(-η|k|2tα)e-ikxdk.
4. Conclusion
In this paper, we have introduced a unified fractional Schrödinger equation and established solution for the same. The solution has been developed in terms of the generalized Mittag-Leffler function in a compact and elegant form with the help of Laplace and Fourier transforms and their inverses. All the results derived in this paper are in a form suitable for numerical computation. The fractional Schrödinger equation discussed in the present article contains a number of known (may be new also) fractional Schrödinger equations. The results obtained in the present paper provide an extension of the results given by Bhatti [12], Chaurasia and Singh [13], Debnath [5], and Saxena et al. [14].
Acknowledgment
The authors are grateful to Professor H. M. Srivastava, University of Victoria, Canada, for his kind help and valuable suggestions in the preparation of this paper.
HilferR.2000River Edge, NJ, USAWorld ScientificBeyerH.KempfleS.Definition of physical consistent damping laws with fractional derivatives199575623635KempfleS.GaulL.Global solutions of fractional linear differential equations19967625715722-s2.0-33749480385SchneiderW. R.WyssW.Fractional diffusion and wave equations19893011341442-s2.0-0001553919DebnathL.Fractional integrals and fractional differential equations in fluid mechanics20036119155DebnathL.20052ndBasel, SwitzerlandBirkhäuserDebnathL.Recent applications of fractional calculus to science and engineering2003200354341334422-s2.0-1784437219310.1155/S0161171203301486LaskinN.Fractals and quantum mechanics20001047807902-s2.0-0034363921LaskinN.Fractional quantum mechanics and Levy path integrals20002684–62983052-s2.0-000041530910.1016/S0375-9601(00)00201-2LaskinN.Fractional Schrödinger equation200266572-s2.0-4134908476110.1103/PhysRevE.66.056108056108NaberM.Distributed order fractional sub-diffusion200412123322-s2.0-454429916010.1142/S0218348X04002410BhattiM.Fractional Schrödinger wave equation and fractional uncertainty principle20072943952ChaurasiaV. B. L.SinghJ.Application of sumudu transform in schödinger equation occurring in quantum mechanics2010457284328502-s2.0-78650298757SaxenaR. K.SaxenaR.KallaS. L.Computational solution of a fractional generalization of the Schrödinger equation occurring in quantum mechanics20102165141214172-s2.0-7795086762510.1016/j.amc.2010.02.041MillerK. S.RossB.1993New York, NY, USAJohn Wiley & SonsKilbasA. A.SrivastavaH. M.TrujilloJ. J.2006Amsterdam, The NetherlandsElsevierCaputoM.1969Bologna, ItalyZanichelliMetzlerR.KlafterJ.The random walk's guide to anomalous diffusion: a fractional dynamics approach200033911772-s2.0-0002641421CompteA.Stochastic foundations of fractional dynamics1996534419141932-s2.0-0000582309PrabhakarT. R.A singular integral equation with a generalized Mittag-Leffler function in the kernel197119715SaxenaR. K.MathaiA. M.HauboldH. J.Solultions of fractional reaction-diffusion equations in terms of Mittag-Leffler functions200615117