Int. J. Contemp. Math. Sciences, Vol., 007, no. 9, 943-950 Fractional Schrödinger Wave Equation and Fractional Uncertainty Principle Muhammad Bhatti Department of Physics and Geology University of Texas Pan American Edinburg, TX 7854-999 USA Abstract Free particle wavefunction of the fractional Schrödinger wave equation is obtained. The wavefunction of the equation is represented in terms of generalized three-dimensional Green s function that involves fractional powers of time as variable t α. It is shown that the wavefunction corresponding to the integral order Schrödinger wave equation follows as special case of that of the corresponding Schrödinger equation with fractional derivatives with respect to time. The wavefunction is obtained using Laplace and Fourier transforms methods and eventually the wavefunction is expressed in terms of Mittag-Leffler function. Heisenberg Uncertainty principle is deduced from the wavefunction of the fractional Schrödinger equation using the integral value of fractional parameter α=. Introduction Recent applications of fractional equations to number of systems such as those exhibiting enormously slow diffusion or sub-diffusion have given opportunity for physicists to study even more complicated systems. Those systems include charge transport in amorphous semiconductors, the relaxation in polymer systems, fluid mechanics and viscoelasticity and Hall Effect. The generalized diffusion equation allows describing complex systems with anomalous behavior in much the same way as simpler systems. Anyway, fractional calculus is now considered as a practical technique in many branches of science including physics []. Several authors including [-9] have discussed many examples of homogeneous fractional ordinary differential equations, homogenous fractional
944 Muhammad Bhatti diffusion equation and homogenous wave equations. Debnath [7-9] considered solutions of fractional order homogeneous and nonhomogenous partial differential equations and integral equations in fluid mechanics. The practical use of fractional calculus is underlined by the fact that Laplace transform of the operator 0 D α t f(t) has the simple form, L { 0D α t f(t) } = s α L{f(t)}, where L represents Laplace transform. In this paper, we are concerned with fractional time derivatives with 0< α<representing as the fractional parameter. The physical significance of the parameter α is unknown. The fractional Schrödinger wave equation may represent wave equation for quasi-particles when the parameter takes fractional values between 0 <α<, otherwise it will represent particle wave equation for α =. We set β = 0 in the α-th fractional derivative β D α t = dn dx n β D α n t, specifies t = 0 as the starting of the system s time evolution with 0 α [3]. If this definition is applied to the diffusion equation for a particle at the origin of coordinates in n-dimensional space, it s mean squared displacement is proportional to time t. In a variety of physical systems, however, the simple scaling law is violated [0, ]. The fractional diffusion equation provided the scaling law as t α for the mean square distance []. In the the following, we use the fractional calculus to solve the Schrödinger wave equation for a free moving particle. The solution to this equation is presented in a closed form and the uncertainty principle is deduced for the special case α= for a free particle. The Fractional Schrödinger Wave Equation We will consider the solution of the fractional Schrödinger equation in the form given by α ψ(x, t) t α = i m ψ(x, t), (.) where = + +, is the Planck s constant divided by π, m x y z is the mass and ψ(x, t) is a wave function of the particle. Also we set a = i~ m as a complex constant and solve this equation (.) with the following initial and boundary conditions ψ(x,0) =ψ 0 (x) ψ(x, t) 0 as x, t>0. (.) We apply the joint Laplace transform with respect to t and the Fourier transform with respect to x (see reference [3]) defined by
Fractional Schrödinger wave equation 945 ψ(k, s) = e i k x d 3 x e st ψ(x, t) dt. (.3) π 0 Where and are used to denote the Laplace and the Fourier transforms respectively, k and s are the Fourier and the Laplace transform variables respectively. Application of the joint transform to equations (.) and (.) give s α ψ(k, s) s α ψ(k, 0) = a (i k) ψ(k, s). (.4) We combine the terms and take the inverse Laplace transform of equation (.4) to yield { } s ψ(k, s) = L α ψ(k, 0) = G s α + a k α (k,s) ψ(k, 0). (.5) Next, we use the formula for the inverse Laplace transform to express the solution in Mittag-Leffler function { } m!s L α β = t mα+β E (m) (s α + a) m+ α,β ( atα ), (.6) where E α,β (z) is the Mittag-Leffler function (see Erdèlyi, 955) defined by the series z m E α,β (z) =,α>0, β > 0 (.7) Γ(mα+ β) m=0 where Γ(x) is the Gamma function and E (m) dm α,β (z) = dx E α,β(z) (.8) m Application of the inverse Laplace transform combined with the formula in equation (.5) and using equation (.6) yields solution ψ(k,t)= ψ(k, 0) E α, ( a k t α )=G α (k,t) ψ(k, 0), (.9) where k = k x +k y +k z and the Green function is defined as G α (k,t) = E α, ( ak t α ) = ( ak t α ) m m=0. Alternatively, the Green function in Eq. Γ(mα+) (.5) in momentum space may be simplified and the complex integration in variable s can be carried out to obtain s α e st G α (k, t) = 0 s α +a k ds = { e ak t for α = sin(πα) πα 0 e (a k x) /α t x +x cos(πα)+ dx for 0 <α< (.0)
946 Muhammad Bhatti Both forms of the Green function in Eqs. (.9) and (.0) agree for α = and α=. For α =, the expression for Green function is G /(k,t) = E /, ( a k t / )=e a k 4 t erfc(ak t / ). Finally, we take the Fourier transform of equation (.9) to obtain the solution of the wave equation (.). ψ(x,t) = e i k x (π) 3/ ψ0 (k) E α, ( a k t α )d 3 k. (.) We assume that the Fourier transform of the initial wave function at t=0 to be as ψ(k) = e i k ξ ψ (π) 3/ 0 (ξ) d 3 ξ. (.) By the convolution theorem of the Fourier transform (see Debnath [3]), the solution in equation (.) may be expressed in the form ψ(x,t) = G α (x ξ,t) ψ 0 (ξ) d 3 ξ, (.3) where the Green function is given by G α (x, t) = e i k x E (π) 3 α, ( a k t α )d 3 k. (.4) For the case α =, the fractional wave equation (.) reduces to the Schrödinger wave equation. In this special case, the solution (.3) after integration over k reduces to the familiar form ψ(x, t) = e (x ξ) 4 at ψ 0 (ξ)d 3 ξ, (.5) (4πat) 3 where the 3-D Green function for α = in spatial coordinates is given by G(x, t) = x e 4 at (4πat) 3/. (.6) In equation (.5), the result E, (z) =e z is used. The solution is in perfect agreement with the standard solution of the integer order α = of the wave equation (.). We also calculated [] the mean square distance <r > in spherical coordinates using the Green function given in equation (.3). The Green function is expressed in spherical coordinates given by
Fractional Schrödinger wave equation 947 G α (r, t) = 4π π (π) 3 0 kr J /(kr) E α, ( ak t α ) k dk. (.7) By taking the Laplace transform of equation (.7) and integrating over k, we get Ḡ α (r, s) = 4π ar sα e r a s α/. (.8) Taking the inverse Laplace transform of the above equation (.8) we get G α (r, t) = 4π (at α ) 3 z W ( z, α, α). (.9) Where W (z, α, β) = z n n=0 is called Wright Function [see Erdèlyi, n!γ(αn+ β) 955] and z = r at α/. We define a new function for odd values of n = m + for β = 0 as: B(z, α) = m=0 z m+ (m + )!Γ(α(m + )). (.0) For the special case α =, the B-function B( z, α ) reduces to simple exponential form of type B( z, )= z π ze 4 (.) and the 3-D Green function in equation (.9) which reduces to simple form given by G (r, t) = (4πat) 3 e r 4 at. (.) It is also convenient to express the Green function in equation (.8) after taking inverse Laplace transform and carrying out complex integration, we obtain G α (r, t)= ( )α+ π ar 0 e tx r a x cos(πα/) x α sin[ r a x sin(πα/)] dx. (.3) For α =, the above equation also gives the same answer as in equation (.).
948 Muhammad Bhatti 3 Uncertainty Principle The Green function in Eq. (.) G (x, t) = e x (4πat) / 4 at for α = is used to verify the uncertainty principle in one dimension (-D). The solution of the Schrödinger equation is then written in the form ψ(r, t)= G (ξ,t) ψ 0 (r ξ) dξ. (3.) We choose normalized initial Gaussian wave function ψ 0 (r) = e (πb ) 4, where constant b represents an initial Gaussian width of the initial wave function. Substituting the expression for Green function in Eq. (3.) and performing the integration, we obtain wave function r ψ(r, t)= (πb ) 4 ( + i t/b ) e b (+i ~ t/b ). (3.) Using this wave function in Eq.(3.), we can obtain the quantum mechanical probability density by taking the complex conjugate and multiplying with the wave function itself given as P (r, t) = ψ(r, t) = π bc e r b r b c, (3.3) where the constant c = ( + ~ t ). Mean square distance is calculated m b 4 using the quantum mechanical probability given by This expression Eq. given by r = π 0 x r e b c dx = bc (bc). (3.4) (3.4) provides uncertainty in the position which is Δr = b ( + t ). m b4 (3.5) It is straightforward to calculate expectation value of the momentum square with integral value of α= [4]. Δp = k = b. (3.6) The uncertainty principle is just the product of both expressions given in equations (3.5 and 3.6)
Fractional Schrödinger wave equation 949 Δr Δp = ( + t m b ) (3.7) 4 The above expression with t=0 reduces to the known uncertainty principle given by Δr Δp =. (3.8) These calculations show that fractional calculus is a powerful tool that may be used to solve the problems in general and the reduction leads to the known results when integral values of α are used. References [] R. Hilfer, ed., Applications of Fractional Calculus in Physics, World Scientific, River Edge, N. J. (000). [] A. Carpinteri, and Mainardi (eds.), Fractals and Fractional Calculus in Continuum Mechanics, Springer Verlag, New York (997). [3] I. Podlubny, Fractional Differential Equations, Academic Press, Boston, 999. [4] H. Beyer and S. Kempfle, Definition of Physically Consistent Damping Laws with Fractional Derivatives, Z. Angew Mech., 75 (995) 63-635. [5] S. Kempfle, and L. Gaul, Global Solutions of Fractional Linear Differential Eqautions, Z. Angew. Math. And Mech., 76(996) 57-57. [6] W. R. Schneider and W. Wyss, Fractional Difussion and Wave Equations, J. Math. Phys., 30 (989) 34-44. [7] L. Debnath, Fractional Integral and Fractional Differential Eqautions in Fluid Mechanics, Fractional Calculus and Applied Analysis, (to appear 003 a). [8] L. Debnath, Recent Developments in Fractional Calculus and Its Applications to Science and Engineering, Internat. Jour. Math and Math. Sci.
950 Muhammad Bhatti 003 (003 b). [9] L. Debnath, Nonlinear Partial Differential Equations for Scientists and Engineers, Birkhauser Verlag, (997). [0] M. F. Shlesinger, G. M. Zaslavsky, J. Klafter, Nature 3, 363 (993). [] J.-P. Bouchaud, A. Georges, Phys. Rep. 95, 7 (990). [] On Fractional Schrödinger and Dirac Equations, Muhammad Bhatti and Lokenath Debnath, International Journal of Pure and Applied Mathematics, 5, - (004). [3] L. Debnath, Integral Transform and Their Applications, CRC Press CityplaceBoca Raton (995). [4] J. J. Sakurai, Advanced Quantum Mechanics, Addison-Wesley (978). Received: December 7, 006