SciELO - Scientific Electronic Library Online

 
vol.68 issue4A dynamical dark energy solution to hubble tension in the light of the multimessenger eraA non-renormalizable neutrino mass model with S3⊗Z2 symmetry author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista mexicana de física

Print version ISSN 0035-001X

Rev. mex. fis. vol.68 n.4 México Jul./Aug. 2022  Epub May 19, 2023

https://doi.org/10.31349/revmexfis.68.040703 

Research

Initial value problem for a Caputo space-time fractional Schrödinger equation for the delta potential

S. Saberhaghparvar1 

H. Panahi1  * 

1Department of Physics, Faculty of Science, University of Guilan, Rasht 51335-1914, Iran.


Abstract

In this paper, we consider a Caputo space-time fractional Schrödinger equation for the delta potential. To solve the equation, we use the joint Laplace and Fourier transforms on the spatial and time coordinates, respectively. After applying the integral transformations, we use the special initial and boundary physical conditions obtained by trial and error; these special initial conditions involve considering the initial spatial wave function in terms of the Mittag-Leffler function. Consequently, using the fractional calculus, we obtain the wave functions and corresponding eigenvalues. Finally, to verify the solution, we recover the standard case corresponding to α →1 and β → 1.

Keywords: The fractional Schrödinger equation; Caputo space-time fractional derivative; Mittag-Leffler function

1. Introduction

The fractional calculus is a generalization of the usual calculus, where derivatives and integrals are defined for arbitrary real numbers. In some phenomena, the fractional operators simulate the phenomena better than ordinary derivatives and normal integrals. Fractional calculus has been widely used in science and engineering1-3; in particular, the study of fractional differential equations such as fractional wave equations have received great deal of attention due to the importance of propagation of electromagnetic waves and vibrating strings4,5. Moreover, the fractional Schrödinger equation has been recently studied in many fields, such as the obstacle problem, phase transitions and anomalous diffusion6-12, etc. Fractional calculus began with Leibniz (1695-1697) and Euler’s speculations (1730); later on, Riemann, Liouville, Grünwald and Letnikov13-16 provided the general definition of fractional derivatives. In 2000, Laskin provided the first application of fractional calculus in quantum mechanics by employing the path integral formulation over Lévy paths and showed that the corresponding equation of motion is indeed the space fractional Schrödinger equation17-20. In Refs.21-24 the authors used the method developed by Laskin to study some new physical applications, its non-locality and the consistency of its solutions. On the other hand, Naber introduced in 2004 the time fractional Schrödinger equation25, where he investigated the wave function of the time fractional Schrödinger equation for free particles and potential wells, and obtained the corresponding wave functions in terms of the Mittag-Leffler function. An overview of quantum mechanics and quantum dynamics with time fractional derivatives can be found in Refs.26,27. In Refs.28,29 the authors solved the space-time fractional Schrödinger equation for free particle and for square potential well using the integral transform approach; moreover, in Ref.30, the Schrödinger equation with Riesz fractional derivative was solved using the scaling transformations. Also, the numerical solution of fractional Schrödinger differential equation with the Dirichlet condition has been investigated in Ref.31. The space-time fractional Schrödinger equation has been obtained by replacing the space and time derivatives of integer order by the derivatives of non-integer order. There are several ways for generalizing an integer-order derivative to the arbitrary case including the Riemann-Liouville and Caputo definitions. In Ref.32, the definitions of the fractional Caputo and Riemann-Liouville derivatives and p-Laplacian operator 𝜙p have been used and the existence of solutions for nonlinear fractional differential has been investigated. Another application of fractional calculus can be found in Ref.33, where the fractional order model HIV/AIDS with the Liouville-Caputo and Atangana-Baleanu-Caputo derivatives has been considered. In this paper, we use the Caputo fractional derivative, defined as13

0cDtαf(t)=1Γ(n-α)0tf(n)(τ)dτ(t-τ)α+1-n,(n-1)<α<n, (1)

where for the case α → n, the Caputo derivative becomes an ordinary n-th order derivative of the function ƒ(t).

We consider the space-time fractional Schrödinger equation for the delta potential, then, by imposing a special initial condition we obtain the corresponding wave function by using the joint Laplace and Fourier transforms with respect to time and coordinates, respectively. The Laplace transform of the Caputo derivative is given by13

L0cDtαft=sαFs-k=0n-1sα-k-1fk0, (2)

where

F(s)=0e-stf(t)dt

The paper is organized as follows. In 2, we introduce the space-time fractional Schrödinger equation and a special initial condition. In 3, we obtain the wave function corresponding to the imposed initial condition in terms of Fox’s H-function by applying the joint Laplace and Fourier transforms. In 4, the time dependent energy eigenvalues are determined by an asymptotic expansion of Fox’s H-function. The paper ends with the conclusions in 5 and two appendices.

2. The space-time fractional Schrödinger equation for delta potential

The one-dimensional time-dependent Schrödinger equation is given by

iφx,tt=-22m2x2φx,t+Vxφx,t, (3)

where m is the mass of the particle and ℏ is the Planck’s constant; in an analogous way, the one-dimensional space-time fractional Schrödinger equation is written as25

iαη0cDtαφ(x,t)=-12mc2(mc)2β0cDx2βφ(x,t)+V(x)φ(x,t), (4)

where 0 < ( ≤ 1, 0.5 < β ≤ 1, η = mc2(ℏ/mc2)(. The scaling factors iαη and -(1/2)mc2(ℏ/mc) have been added to equalize units on both sides of the equation. It is easy to see that for the special case ( = 1 and β = 1, the space-time fractional Schrödinger equation (4) reduces to the ordinary Schrödinger equation (3). In Ref.22, Eq. (4) was studied for a free particle and for a potential well assuming only the space fractional state; the wave functions and the corresponding eigenvalues were calculated. It was observed that in the case of the free particle, the wave function could be obtained in terms of the Fox’s H-function while for the potential well, the method of separating the variables was used and the solutions of the system were obtained in terms of the space fractional parameter; the some problems were investigated in Ref.25 assuming that only the time fractional case was taken into account. The wave function in the case of the free particle was obtained as the time dependent Mittag-Leffler function multiplied by the spatial function of the free particle, while in the case of the potential well, energy dependent on the time fractional parameter was obtained. Further, in Ref.34, the space fractional Schrödinger equation for single and double Dirac delta potential was considered; the wave functions were determined in terms of the Fox H-function, and the energy eigenvalues were obtained by using the Fourier transform and the Riesz fractional derivative; while the authors assume only the spatial derivative of Schrödinger equation as the Riesz fractional derivative, here, we intend to investigate the Caputo space-time fractional Schrödinger equation for the Dirac delta potential V(x) = -v0 δ(x), that is

iαη0cDtαφ(x,t)=-12mc2(mc)2β0cDx2βφ(x,t)-v0δ(x)φ(x,t). (5)

We solve the problem analytically by using the joint Laplace and Fourier transforms; to this aim, first we consider the following physical boundary and initial conditions

φ(x,0)=g(x),φ(x,t)0as|x|,φ(0,t)=0n1,φx(n)(x,0)=0, (6)

Choosing an arbitrary function for g(x) one can solve the problem numerically, but in this article we try to choose a specific function for the special case for g(x) so that the problem can be solved analytically. As a result, we shall show that for the special case α = β = 1, the solution of the ordinary Schrödinger equation for Dirac delta potential is obtained. Now to solve the Eq. (5), we apply the joint Laplace transform on the spatial coordinate and the Fourier transform on the time coordinate defined by35

φ~¯(k,s)=12π-+e-ikxdx0e-stφ(x,t)dt, (7)

where the symbols for (-) and (∼) are used to denote the Laplace and Fourier transforms, respectively, and k and s are the Fourier and Laplace transform variables.

3. Wave function of one-dimensional space-time fractional Schrödinger equation for the imposed initial condition

By applying the joint Laplace and Fourier transforms to Eq. (5) and using the boundary conditions of Eq. (6), we obtain the following equation

φ̃¯k,s=sα-1g̃ksα+Fik2βiαη, (8)

where F=(1/2)mc2(ℏ/mc). The inverse Laplace transform of Eq. (8) gives

φ̃k,t=g̃kEα-Fik2βtαiαη. (9)

In the above calculation, we have used the following formula36

L-1m!sα-βsα±a2m+1=tαm+β-1Eα,βma2tα, (10)

where Eα,β (z) is the Mittag-Leffler function defined as37

Eα,βz=n=0znΓnα+β,α>0,βC. (11)

Now, we assume the function g(x) as

gx=AEα'-λxα',α'C. (12)

This choice, as mentioned above, can help us to solve the problem analytically and obtain the analytical solution of the ordinary Schrödinger equation, α = 1 and β = 1. Needless to say, for other choice of the g(x), the analysis is difficult, and perhaps impossible; accordingly, the problem must be solved numerically. Applying the Fourier transform to the relation (12), we get

g~(k)=A2π-+e-ikxEα'(-λ|x|α')dx=2A2π0os(kx)H1,21,1λxα'(0,1)(0,1),(0,α')dx, (13)

where Hp,qm,n&#091;x|-&#093; denotes the Fox’s H-function, which, in terms of the Mittag-Leffler function is given by38

H1,21,1-z0,10,1,0,α'=Eα'z. (14)

In Appendices A and B, we have reproduced some useful relations about on the Fox’s H-function and its integral according to38; using (B3), (13) is calculated as

g̃k=2A2λ1α'H2,32,11λk2α'1-1α',1,0,α'0,α'2,1-1α',1,12,α'2. (15)

Hence, by substituting (21) into (9), we get

φ̃k,t=2A2λ1α'H2,32,11λk2α'1-1α',1,0,α'0,α'2,1-1α',1,12,α'2Eα-Fik2βtαiαη. (16)

Now by calculating the inverse Fourier transform of (16), using the integral containing of two Fox’s H-functions, an taking into account the properties of Fox’s H-function38,39, we get the wave function φ(x,t) in the following form

φx,t=2Ax-1α'βλ1α'H2,0:1,2;2,30,1:1,1;2,1-4Fiαη1βtαβx2λ-2α'x212;1,1,0;1,1:0,1β;1-1α',2α',0,2-:0,1β,0,αβ;0,1,1-1α',2α',12,1, (17)

where β > - (1/2), (α/β) > 0, α + β < 2, ά < (4/3); the definition of the Fox’s H-function of two variables

Hxy

is given in Eqs. (A.4-A.7). It should be noted here that the acceptable values of α and β above are consistent with the values of α and β, i.e., 0 < α ≤ 1 and 0.5 < β ≤ 1. Before calculating the energy eigenvalues, let us show that the obtained wave function coincides with the wave function of the ordinary Schrödinger equation for the delta function potential; to this aim, we inspect the wave function 𝜙(x, t) in limits of α ( 1, β ( 1, ά ( 1. Using the integral relations of the Fox’s H-function, Eqs. (A.4-A.7), and after some calculations, we arrive at

φ(x,t)α1β1α'1=2A|x|-1λH2,0:0,1;1,10,1:1,1;1,1-4Ftiη|x|2λ-2|x|2|12;1,1,(0;1,1):-;(0,1)-:(0,1);(0,1)=2A|x|-1λπ12πin=0(-1)nn!-4Ftiη|x|2L1Γ(-η')Γ(1+η')Γ(12+n+η')Γ(-n-η')4λ-2|x|2n'dη'. (18)

Then, using the properties of the Gamma function and Fox’s H-function given in Ref.40, the function in (18) is simplified to

(x,t)α1β1α'1=2A|x|-1λπn=01n!-4Ftiη|x|2nH2,00,24|x|2λ2(12-n,1),(-n,1)-=AeiFtλ2ηe-λ|x|=Aeitλ22me-λ|x|(19)

Therefore, up to a normalization coefficient A, it is easy to see that the above equation is in agreement with the wave function of the ordinary Schrödinger equation for the delta function potential; choosing λ=v×m/2 and A=v×m/, the normalized wave function is obtained.

In Fig. 1 we have plotted the probability of the particle presence for some values of β by considering α = ά = 1. Also, the time-dependent effect for different choices of α and β is illustrated in Fig. 2. From Fig. 1 it can be seen that by increasing the parameter β, the probability of the particle presence decreases; the result coincides with Fig.1 in Ref.36 for the non-fractional time derivative of the Schrödinger equation. It can be seen from Fig. 2 that the probability of the particle presence decreases over time; the probability decreases in the central points (x = 0) and increases at the points far from the center. In Fig. 2 the dependence of the probability on time at those points close to the well (x = 0) illustrated for both fractional modes, space and time. Also, it is seen that by decreasing the fractional state (α ≈ 1 and β ≈ 1), the dependence of the probability on time disappears.

Figure 1 Plot of Eq. (17) as a function of x for α = 1, ά = 1, t = 0.1 and different values of β; for simplicity, we have chosen m = ℏ = c = λ = 1. 

Figure 2 The probability of presence of the particle in dependence on x for different values of time. 

In the next section, we shall obtain the energy eigenvalues of the system.

4. The energy eigenvalues of the delta potential for the imposed initial condition

According to Ref.41, for the space-time Schrödinger equation, the energy eigenvalues can be obtained from the equation

Eα=iαη-+φ*(x,t)0cDtαφ(x,t)dx. (20)

First, let us rewrite the wave function (18) in another form of the Fox’s H-function as

φ(x,t)=Aiα'β(λ)1α'(Fiαη)12βtα2βH2,0:1,2;2,30,1:1,1;2,1×-4(Fiαη)1βtαβ|x|2λ-2α'|x|2|(1;1,1),(12;1,1):(12β,1β);(1-1α',2α'),(0,2)-:(12β,1β),(α2β,αβ);(0,1),(1-1α',2α'),(12,1). (21)

Then, using the Caputo fractional derivative of the Fox’s H-function, given in Ref.40, one can calculate the derivative of the first order of the Fox’s H-function of two variables; substituting the obtained result into Eq. (20), we get

Eα=2iαηA2t-α-αβ(α')2β2λ2α'(F2η2)12β0+H2,0:1,2;2,30,1:1,1;2,1-4F(-i)αη1βtαβx2λ-2α'x2(1;1,1),(12;1,1):(12β,1β);(1-1α',2α'),(0,2)-:(12β,1β),(α2β,αβ);(0,1),(1-1α',2α'),(12,1)×H2,0:1,2;2,30,1:1,1;2,1-4Fiαη1βtαβx2λ-2α'x2(1;1,1),(12;1,1):(12β,1β);(1-1α',2α'),(0,2)-:(12β,1β),(α+α2β,αβ);(0,1),(1-1α',2α'),(12,1)dx.

To compute the integral, we first express both Fox’s H-functions of two variables in terms of the Mellin-Barnes type integral given in Eqs. (A.4-A.7); then, following the procedure in Ref.35, we obtain

Eα=2iαηA2t-α-αβ(α')2β2λ2α'Fη1β-14π2L1L2L3Γ(1α'-2α't')Γ(t')Γ(1-1α'+2α't')Γ(12-t')Γ(2t')(λ2α')t'Γ(1-12β-1βξ')Γ(12β+1βξ')Γ(1-α-α2β-αβξ')-14F(i)αη1βtαβξ'×(λ2α')η'Γ(-ξ'-η')Γ(12+ξ'+η')0-14π2L4Γ(-s-t')Γ(12+s+t')Γ1-12β-1βsΓ12β+1βsΓ1-α2β-αβs×-x24F(-i)αη1βtαβsdsx2η'+2ξ'+2t'dxdξ'dt'dη'. (23)

Now, if we use the Mellin transform theorem and substitute s → -η’ -t’ -ξ’ - (½) in Eq. (23), the energy eigenvalue is obtained in terms of the Fox’s H-function of three variables as

Eα=2iαηA2t-α-αβ(α')2β2λ2α'Fη1β-4F(-i)αη1βtαβ12L1L22L3Γ1α'-2α't'Γ(t')Γ1-1α'+2α't'Γ12-t'Γ(2t')Γ1-12β-1βξ'Γ12β+1βξ'Γ1-α-α2β-αβξ'Γ1+ξ'β+t'β+η'βΓ-ξ'β-t'β-η'βΓ1+αξ'β+αt'β+ηα'β-4λ2α'F(-i)αη1βtαβt'×-4λ2α'F(-i)αη1βtαβη'(-1)αξ'βdt'dξ'dη'=2iαηA2t-α-αβ(α')2β2λ2α'Fη1β-4F(-i)αη1βtαβ12H1,2:2,1;3,2;3,21,1:1,1;1,2;1,2(-1)αβ-4λ2α'F(-i)αη1βtαβ-4λ2α'F(-i)αη1βtαβ(0;1β,1β,1β):(1-12β,1β),(1-α-α2β,αβ);(1,1),(1α',2α'),(12,1);(1,1),(1α',2α'),(12,1)(0;1β,1β,1β),(0;αβ,αβ,αβ):(1-12β,1β);(1α',2α'),(1,2);(1α',2α'),(1,2), (24)

which is valid under the following conditions

β>-14,12<α'<43,α+β<2.

Thus, it is seen that the energy eigenvalue (24) is time dependent. Also for complex values of the parameters α and β, the energy eigenvalue E may be complex for an unstable system.

Now, let us investigate the energy eigenvalue (24) for the special cases α → 1, β → 1 and ά → 1, that is

Eαα1β1α'1=-2iη2A2t-2λ2Fπ4Ftiη12 L1 L2  L3Γ(-ξ'-η'-t')Γξ'+32× Γt'Γ1-t'Γη'Γ1-η'-1ξ'λ2Ftiηt'λ2Ftiηη'dξ'dt'dη'==-4η32A2(t)-32(i)12λ2πF12n=0in!H1,0:1,1;1,10,1:1,1;1,1iηλ2Ftiηλ2Ft-12-n:0,1;0,1-:0,1;0,1=-4(η)32A2(t)-32(i)12λ2π(F)12n=0in!H1,22,1λ2Ftiη(1,1)(32+n,1),(2,1) (25)

In the above relation, we have used identities of the Fox’s H-function of two variables given in Refs.38-40. We also need to use the expansion of the Fox’s H-function at small and large times. For small values of t → 0, using the Fox’s H-function40, the energy is obtained as E=-A2Fλ=-mv02/22, which is in accordance with the energy eigenvalue of the ordinary Schrödinger equation. On the other hand, the energy eigenvalue tends to zero as t → ∞.

5. Conclusions

In this paper, we have studied the Caputo space-time fractional Schrödinger equation for the delta potential. To solve this equation, we have used the joint Laplace transform on the spatial coordinate and the Fourier transform on the time coordinate; then we have obtained an assumption by trial and error, so that the problem can be solved analytically. In other words, assuming that the initial wave function is the Mittag-Leffler function, then the wave function and the energy eigenvalues can be obtained in terms of the Fox’s H-function of two and three variables, respectively. On the other hand, assuming other functions as the initial wave function, solving the corresponding equation is difficult and it must be solved numerically. We also investigated the wave function for particular values of the fractional parameters α → 1 and β → 1 and obtained the results of the standard Schrödinger equation with the delta potential. Moreover, we studied the behavior of the particle presence probability. In Fig. 1, for the case of non-fractional time derivative (α = 1), it was observed that by increasing the spatial fractional parameter (, the probability of the presence decreases which is consistent with the result given in Ref.36. In Fig. 2, we have demonstrated the dependence of probability on the time while it was not the case in the standard Schrödinger equation. We have also shown that the probability of presence is decreased as time increases such that the probability is decreased in the central points (x = 0) and at points farther from the center, the probability of presence is increased. It was also observed that by decreasing the fractional state (α ≈ 1 and β ≈ 1), the dependence of the probability of presence on time disappears. In the end, we investigated the energy eigenvalue for the special cases α → 1, β → 1 and ά = 1; then, we checked the approximate behavior of energy eigenvalue at small and large times. It was observed that while for large times, the energy of the system dissipates, for small times, the energy is in accordance with that of the standard Schrödinger equation for the delta potential.

References

1. A. Iomin, G. M. Zaslavsky, Quantum manifestation of Lévy-type flights in a chaotic system. Chem. Phys. 284 (2002) 3, htt://doi.org/10.1016/S0301-0104(02)00532-3. [ Links ]

2. R. Metzler and J. Klafter, The random walk’s guide to anomalous diffusion: a fractional dynamics approach. Phys. Rep. 339 (2009) 1, htt://doi.org/10-1016/S0370-1573(00)00070-3. [ Links ]

3. R. Metzler and J. Klafter, the restaurant at the random walk: recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A 37 (2004) 161, htt://doi.org/10.1088/0305-4470/37/31/R01. [ Links ]

4. A. Khan, T. S. Khan, M. I. Syam, and H. Khan, Analytical solutions of time fractional wave equation by double Laplace transform method. Eur. Phys. J. Plus. 134 (2019) 163, htt://doi.org/10.1140/epjp/i2019-12499-y. [ Links ]

5. Z. M. Odibt, S. Momani, Approximate solutions for boundary value problems of time fractional wave equation. Appl. Math. Comput. 181 (2006) 767, htt://doi.org/10.1016/j.amc.2006.02.004. [ Links ]

6. M. Cheng, Bound state for the fractional Schrödinger equation with unbounded potential. J. Math. Phys. 53 (2012) 043507, htt://doi.org/10.1063/1.3701574. [ Links ]

7. P. Felmer, A. Quaas, and J. Tan, Positive solutions of the non-linear Schrödinger equation with the fractional Laplacian. Proc. Roy. Soc. Edinb. A 142 (2012) 1237, htt://doi.org/10.1017/S0308210511000746. [ Links ]

8. B. Guo and D. Huang, Existence and stability of standing waves for nonlinear fractional Schrödinger equations. J. Math. Phys. 53 (2012) 15, htt://doi.org/10.1063/1.4746806. [ Links ]

9. N. Laskin, Fractional Schrödinger equation. Phys. Rev. E. 66 (2002) 056108, htts://doi.org/10.1103/PhysRevE.66.056108. [ Links ]

10. S. Secchi, Ground state solutions for nonlinear fractional Schrödinger equations in ℝN. J. Math. Phys. 54 (2013) 031501 htt://doi.org/10.1063/1.4793990. [ Links ]

11. L. Silvestre, Regularity of the obstacle problem for a fractional power of the Laplace operator. Comm. Pure Appl. Math. 60 (2007) 67, htt://doi.org/10.1002/cpa.20153. [ Links ]

12. W. Zhang, X. H. Tang and J. Zhang, Infinitely many radial and non-radial solutions for a fractional Schrödinger equation. Comput. Math. Appl. 71 (2016) 737, htt://doi.org/10.1016/j.camwa.2015.12.036. [ Links ]

13. I. Podlubny, Fractional Differential Equations, 1st ed. (Academic Press, San Diego, 1998). [ Links ]

K. S. Millar, and B. Ross, An introduction to the fractional calculus and fractional differential equations, (John Wiley and Sons, New York, 1993). [ Links ]

15. K. Oldham, and J. Spanier, The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order (Elsevier Science, Amsterdam, 1974) [ Links ]

16. J. Tenreiro Machado, V. Kiryakova, and F. Mainardi, Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 16 (2011) 1140, htt://doi.org/10.1016/j.cnsns.2010.05.027. [ Links ]

17. N. Laskin, Fractional quantum mechanics and Lévy path integrals. Phys. Lett. A 268 (2000) 298, htt://doi.org/10.1016/S0375-9601(00)00201-2. [ Links ]

18. N. Laskin, Fractional quantum mechanics. Phys. Rev. E. 62 (2000) 3135, htt://doi.org/10.1103/PhysRevE.62.3135. [ Links ]

19. N. Laskin, Fractals and quantum mechanics. Chaos. 10 (2000) 780, htt://doi.org/10.1063/1.1050284. [ Links ]

20. N. Laskin, Fractional Schrödinger equation. Phys. Rev. E. 66 (2002) 056108, htts://doi.org/10.1103/PhysRevE.66.056108. [ Links ]

21. S. Ş. Bayin, On the consistency of the solutions of the space fractional Schrödinger equation. J. Math. Phys. 53 (2012) 042105, htts://doi.org/10.1063/1.4705268. [ Links ]

22. X. Guo and, M. Xu, Some physical applications of fractional Schrödinger equation. J. Math. Phys. 47 (2006) 082104, htts://doi.org/10.1063/1.2235026. [ Links ]

23. M. Jeng, S.-L.-Y. Xu, E. Hawkins, and J. M. Schwarz, On the locality of the fractional Schrödinger equation. J. Math. Phys. 51 (2010) 062102, htts://doi.org/10.1063/1.3430552. [ Links ]

24. E. K. Lenzi, B. F. de Oliveira, L. R. da Silva, and L. R. Evangelista, Solutions for a Schrödinger equation with a nonlocal term. J. Math. Phys. 49 (2008) 032108, htts://doi.org/10.1063/1.2842069. [ Links ]

25. M. Naber, Time fractional Schrödinger equation. J. Math. Phys. 45 (2004) 3339, htts://doi.org/10.1063/1.1769611. [ Links ]

26. A. Iomin, Fractional-time quantum mechanics. Phys. Rev. E. 80 (2009) 022103, htts://doi.org/10.1103/PhysRevE.80022103. [ Links ]

27. A. Iomin, Fractional-time Schrödinger equation: Fractional dynamics on a comb. Chaos, Solitons & Fractals. 44 (2011) 348, htts://doi.org/10.1016/j.chaos.2011.03.005. [ Links ]

28. J. Dong and M. Xu, Space-time fractional Schrödinger equation with time-independent potentials. J. Math. Phys. 344 (2008) 1005, htts://doi.org/10.1016/j.jmaa.2008.03.061. [ Links ]

29. S. Wang and, M. Xu, Generalized fractional Schrödinger equation with space-time fractional derivatives. J. Math. Phys. 48 (2007) 043502, htts://doi.org/10.1063/1.2716203. [ Links ]

30. A. Patra, Similarity analytical solutions for the Schrödinger equation with the Riesz fractional derivative in quantum mechanics. Mathematical methods in the applied sciences. 43 (2020) 10287, https://doi.org/10.1002/mma.6695. [ Links ]

31. A. Ashyralyev, B. Hicdurmaz, On the numerical solution of fractional Schrödinger differential equations with the Dirichlet condition. Int. J. Comput. Math. 89 (2012) 1927, htts://doi.org/10.1080/00207160.2012.698841. [ Links ]

32. A. Khan, M. I. Syam, A. Zada and, H. Khan, Stability analysis of nonlinear fractional differential equations with Caputo and Riemann-Liouville derivatives. Eur. Phys. J. Plus. 133 (2018) 264, htts://doi.org/10.1140/epjp/i2018-12119-6. [ Links ]

33. A. Khan, J. F. Gómez-Aguilar, T. S. Khan and, H. Khan, Stability analysis and numerical solutions of fractional order HIV/AIDS model. Chaos, Solitons & Fractals. 122 (2019) 119, htts://doi.org/10.1016/j.chaos.2019.03.022. [ Links ]

34. E. C. de Oliveira, F. Silva Costa, J. Vaz, The fractional Schrödinger equation for delta potentials. J. Math. Phys. 51 (2010) 123517, htts://doi.org/10.1063/1.3525976. [ Links ]

35. L. Debnath, Integral Transform and Their Applications, (CRC Press: Boca Raton, 1995). [ Links ]

36. T. R. Prabhakar, A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math. J. 19 (1971) 7. [ Links ]

37. R. Gorenflo, A. A. Kilbas, F. Mainardi, S. V. Rogosin, Mittag-Leffler Functions: Related Topics and Applications, (Springer-Verlag, Berlin-Heidelberg, 2014), https://doi.org/10.1007/978-3-662-43930-2. [ Links ]

38. A. M. Mathai, R. K. Saxena, H. J. Haubold, The H-function: Theory and Applications, (Springer, New York, 2009). [ Links ]

39. A. M. Mathai, R. K. Saxena, The H-function with applications in statistics and other disciplines, (John Wiley and Sons, New York, 1978). htts://doi.org/10.2307/1267776. [ Links ]

40. H. M. Srivastava, K. C. Gupta, S. P. Goyal, The H-function of one and two variables with applications, (South Asian Publishers, New Delhi, 1982). [ Links ]

41. S. Ş. Bayin, Time fractional Schrödinger equation: Fox’s H-functions and the effective potential. J. Math. Phys. 54 (2013) 012103, htts://doi.org/10.1063/1.4773100. [ Links ]

Appendix

A. The Fox’s H-function

According to Refs.38,41, the Fox’s H-function of several variables is defined as

Hp,q,[pr,qr]m,n:[mr,nr]xy(aj;αj1,αjr)1,p:{(cj1,γj1)1,pr}(bj;βj1,βjr)1,q:{(dj1,δj1)1,qr}=1(2πi)rL1Lrϕ(ξ1,,ξr)i=1rθi(ξi)xiξidξi, (A.1)

where

ξ1,,ξr=j=1nΓ1-aj+i=1rαj(i)ξij=1mΓbj-i=1rβj(i)ξij=n+1pΓaj-i=1rαj(i)ξij=m+1qΓ1-bj+i=1rβj(i)ξi, (A.2)

θi(ξi)=j=1niΓ1-cj(i)+γj(i)ξij=1miΓdj(i)-δj(i)ξij=ni+1piΓcj(i)-γj(i)ξij=mi+1qiΓ1-dj(i)+δj(i)ξi. (A.3)

An empty product is always interpreted as unity. Here m, n, p,q, mi, ni, pi, qi, i = 1, …, r are non-negative integers such that 0 ≤ np, 0 ≤ mq, q ≥ 0, 0 mi, ≤ qi and 0 ≤ nipi; the parameters αj(i),βj(i),γj(i),δj(i) are all positive. For ∀i = 1, … r; x1, … xr are complex variables and Lr is a suitable contour of the Mellin-Barnes type running from -i∞ to +i∞, in complex ξi -plane.

For the Fox’s H-function of two variables, by considering r = 2, m = 0 in Eq. A1, we have

Hx,y=Hxy=Hp1,q1:p2,q2;p3,q30,n1:m2,n2;m3,n3xyai;αi,Ai1,p1:ci,γi1,p2;ei,Ei1,p3bj;βj,Bj1,q1:dj,δj1,q2;fj,Fj1,q3=-14π2L1L2ϕ1ξ,ηθ2ξθ3ηxξyηdξdη, (A.4)

where

ϕ1ξ,η=j=1n1Γ1-aj+αjξ+Ajηj=n1+1p1Γaj-αjξ-Ajηj=1q1Γ1-bj+βjξ+Bjη, (A.5)

θ2ξ=j=1n2Γ1-cj+γjξj=1m2Γdj-δjξj=n2+1p2Γcj-γjξj=m2+1q2Γ1-dj+δjξ, (A.6)

θ3η=j=1n3Γ1-ej+Ejηj=1m3Γfj-Fjηj=n3+1p3Γej-Ejηj=m3+1q3Γ1-fj+Fjη (A.7)

Also for the Fox’s H-function of one variable we have

H[x]=Hp,qm,nxaj,αj1,pbj,βj1,q=Hp,qm,nxa1,α1,,ap,αqb1,β1,,bq,βq=12πiLθ(s)xsds=12πiLΠj=1qΓ(bj-βjs)Πj=1qΓ(1-aj+αjs)Πj=m+1qΓ(1-bj+βjs)Πj=m3+1qΓ(aj-αjs). (A.8)

B. Some Useful Relations and an Integral of the Fox’s H-function

Some useful properties of the Fox’s H-function are listed below36.

Property 1: Transformation of Fox’s H-function with argument x to one with argument 1/x.

Hp,qm,nxap,Apbq,Bq=Hq,pn,m1x1-bq,Bq1-ap,Ap. (B.1)

Property 2: If σ ∈ C then we can write

xσHp,qm,nxap,Apbq,Bq=Hp,qm,nxap+σAp,Apbq+σBq,Bq. (B.2)

Property 3: Cosine transform of Fox’s H-function

0xp-1cos(ax)Hp,qm,nbxσ(ap,Ap)(bq,Bq)dx=2p-1πaρHp+2,qm,n+1b(2a)σ2-ρ2,σ2,(ap,Ap),1-ρ2,σ2(bq,Bq), (B.3)

where a, Δ, σ > 0, p, b, ∈ ℂ; |arg b|( < (1/2)πΔ;

R(ρ)+σmin1jmbjBj>0;R(ρ)+σmax1jnaj-1Aj<1, (B.4)

and

Δ=j=1nAj-j=n+1pAj+j=1mBj-j=m+1qBj. (B.5)

Received: August 04, 2021; Accepted: November 06, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License