Fractional Schrödinger Equation in the Presence of the Linear Potential

with the Riesz space-fractional derivative of order 0 < α ≤ 2 in the presence of the linear potential V(x) = βx. The wave function to the one-dimensional Schrödinger equation in momentum space is given in closed form allowing the determination of other measurable quantities such as the mean square displacement. Analytical solutions are derived for the relevant case of α = 1, which are useable for studying the propagation of wave packets that undergo spreading and splitting. We furthermore address the two-dimensional space-fractional Schrödinger equation under consideration of the potential V(ρ) = F · ρ including the free particle case. The derived equations are illustrated in different ways and verified by comparisons with a recently proposed numerical approach.


Introduction
The time-dependent Schrödinger equation is a fundamental equation in quantum mechanics for studying the dynamics and evolution of wave packets over time.Several years ago, the classical Schrödinger equation has been generalized to a fractional partial differential equation that takes into account the Riesz space-fractional derivative instead of the conventional Laplacian [1,2].Apart from quantum mechanics, there are many other equations occurring in science that have been reconsidered in terms of fractional derivatives such as the diffusion-wave equation [3][4][5][6], the Langevin equation [7] or the radiative transport equation [8].Analytical solutions to fractional differential equations are, in general, not available in terms of elementary functions.As an example, the fundamental solution to the Schrödinger equation for a free particle can only be written in closed form under consideration of the Fox H-function [9].Very recently, the space-fractional Schrödinger equation has been employed for studying the propagation dynamics of wave packets in the presence of the harmonic potential [10] as well as of the free particle [11].In addition, the fractional Schrödinger equation subject to a periodic P T -symmetric potential has been used to investigate the conical diffraction of a light beam [12].An optical realization of the space-fractional Schrödinger equation, based on transverse light dynamics in aspherical optical cavities, has been recently proposed in the study [13].Besides the free particle and the harmonic potential, the case of a linear potential is also a fundamental problem in quantum mechanics that can be treated and solved analytically [14,15].In this context, we refer to the paper [16] dealing with the nonlinear Schrödinger equation in the presence of uniform acceleration.We additionally note on the time-dependent Schrödinger equation with a nonlocal term that has been analyzed in the publication [17].
In this article, we consider the time-dependent Schrödinger equation with the Riesz space-fractional derivative of order 0 < α ≤ 2 in the presence of the linear potential.In the one-dimensional (1D) case, the wave function in momentum space is derived for an arbitrary fractional order α in closed form, which allows the determination of other important quantities such as the mean square displacement (MSD).For the case of α = 1, we provide some explicit analytical solutions useable for studying the time evolution of wave packets that undergo both spreading and splitting.We additionally address the space-fractional Schrödinger equation in two spatial dimensions.The obtained equations are illustrated in different ways and partly verified by comparisons with a recently proposed numerical approach.

1D Space-Fractional Schrödinger Equation
The 1D space-fractional Schrödinger equation in the presence of the linear potential V(x) = βx is given by where ψ(x, t) is the wave function, 0 < α ≤ 2 is the order of the Riesz space-fractional derivative and β ∈ R. The nonlocal operator occurring in Equation ( 1) is defined in the sense of [18] F {(−∆) with Alternatively, it can also be defined directly in real space in terms of the left-and right-sided Riemann-Liouville derivatives [19].In this section, Equation ( 1) is considered under the initial and boundary conditions where ψ 0 (x) is the initial wave packet that is normalized according to The corresponding Schrödinger equation in momentum space is given by where ψ(ω, 0 + ) = ψ 0 (ω).In order to solve this equation, we introduce a wave function of the form ψ(ω, t) = e iϕ(ω) u(ω + βt), where ϕ and u are unknown.Inserting this function into Equation ( 4) results in where u denotes the first derivative of the unknown function u.Comparing both sides of Equation ( 5) leads to the equation which is solved by the function with C being an arbitrary constant.The unknown quantity u can be found by employing the initial condition according to Therefore, the general solution of the Schrödinger Equation ( 4) in momentum space is found to be , which is a measure for the probability density in momentum space, is independent of the fractional order.The wave function in real space is formally given by In general, this integral must be carried out numerically using one of the established quadrature rules.For the classical case of α = 2, the wave function can also be given in terms of the convolution The probability density function (pdf) associated with the position of the particle is defined as ρ(x, t) = |ψ(x, t)| 2 .Furthermore, in view of experimental activities, the MSD is known to be an important quantity that can be computed by means of the pdf as x 2 (t) = ∞ −∞ x 2 ρ(x, t) dx.In our case, when ρ(x, t) is not given, the MSD can be determined from the wave function in momentum space according to Using Equation (9), one obtains the MSD for an arbitrary initial stage as In the following, we also want to address the case of the free particle which is obtained for β = 0.The corresponding solution in momentum space Equation (9) reduces to However, also in this case, the wave function in real space is not available in elementary form.Instead, we have with G(x, t) being Green's function, which is given for 1 < α ≤ 2 by The Fox H-function occurring above is defined in terms of the Mellin-Barnes integral [9] H m,n p,q z where Γ(s) is the Gamma function and L denotes an appropriate integration path in the complex plane.The fundamental solution Equation ( 16) can be evaluated by means of the series where G(0, t) = (−2i/t) 1/α Γ(1/α)/(απ).The corresponding MSD for β = 0 becomes where sign(x) is the sign function.It can be verified that the MSD Equation ( 19) as function of time describes for all fractional orders a parabola.As an example, for the commonly applied Gaussian wave packet one obtains in the free particle case the MSD which exists for all 1/2 < α ≤ 2. Figure 1 displays the MSD Equation ( 21) for different values of the fractional order.

Explicit Solutions for the Square Root of the Laplacian (−∆)
In this subsection, we provide some explicit analytical solutions of Equation ( 1) for the order α = 1.In this case, the fractional Laplacian can be represented as [20] (−∆) where H is the Hilbert transform that is defined by and P denotes the Cauchy principal value.It should be noted that, in recent time, the case of α = 1 has attracted attention in the context of the propagation dynamics of wave packets under the use of the space-fractional Schrödinger equation [10,11,13].In this study, we investigate the time evolution of the Cauchy initial wave packet where ε > 0 is a measure for the beam width.We afore have to note that the time evolution of the Gaussian wave packet and the Cauchy wave packet shows a similar behavior.However, in the latter case, we are able to derive simple analytical solutions that are easy to implement and reproduce.Our starting point is the evaluation of the Fourier integral Equation ( 10) under consideration of the wave function Equation (9).For this, we also need the Cauchy wave packet Equation ( 24) in momentum space, that is ψ 0 (ω) = √ 2πε exp(−ε|ω|).Splitting the integrand into three sections and assuming that β ≥ 0 enables the derivation of the following closed-form expression where erf(z) is the error function.Note that the case β < 0 can be treated in the same manner.The associated MSD belonging to this process can be obtained from Equation ( 13) and summarized as In this context, we note on the limits In the case of the free particle, there is the possibility to provide the general solution for an arbitrary initial stage.In doing so, we need the inverse Fourier transform of the impulse response G(ω, t) = exp(−it|ω|/2), which can be written as Inserting this finding into Equation ( 15) leads to the general solution of the Schrödinger equation In the actual case of the Cauchy initial wave packet Equation (24), we find the time evolution which can be confirmed via Equation ( 25) by setting β = 0.The resulting pdf belonging to this process is In addition, the cumulative distribution function F(x, t) = x −∞ ρ(y, t) dy is found to be Within the classical quantum mechanics, it is known that a pdf belonging to the free particle solution typically spreads out but it does not split.However, as has been very recently shown, the pdf belonging to the space-fractional Schrödinger equation additionally undergoes splitting [11].For the Cauchy wave packet, there is the possibility to provide some more details concerning this feature.Taking the first derivative of the pdf Equation (31) with respect to the space and setting ρ x (x, t) = 0, yields three critical positions, namely x = 0 or x = ± √ t 2 /4 − ε 2 .From the behavior of ρ(x, t), it follows that x = 0 is an absolute maximum as long as 0 ≤ t ≤ 2ε.On the other side, for time values t > 2ε, the position x = 0 becomes a relative minimum.At the same time, we obtain two absolute maxima located at x = ± √ t 2 /4 − ε 2 .Consequently, the probability density spreads out for t > 0, and additionally, starts to split at t = 2ε.For illustration purposes, Figure 2 shows the time evolution of a Cauchy initial wave packet according to Equation (31).A similar behavior can be observed in the case of the Gaussian initial wave packet Equation (20), where the critical time regarding the splitting process must be determined numerically from the wave function where erfi(x) is the imaginary error function and ψ 0 (x) is given by Equation (20). Figure 3 displays the same situation as shown in Figure 2 for the the Gaussian initial wave packet Equation (20).Concerning an experimental verification/realization of the derived theoretical results we refer to the recent publication [11], which proposes an optical system composed of two convex lenses and a phase mask.We furthermore note on the free particle solution to the classical Schrödinger equation.In the case of the Cauchy wave packet Equation ( 24), we have where w(z) = exp(−z 2 ) erfc(−iz) is the Faddeeva function.The corresponding MSD belonging to this process is x 2 (t) = t 2 /(2ε 2 ) + ε 2 .Figure 4 displays the time evolution of the wave function Equation (34) for the same time values considered in Figure 2. It can can be seen that the probability density indeed spreads out but it does not split.

2D Space-Fractional Schrödinger Equation
The 2D space-fractional Schrödinger equation in the presence of the potential where F = (F 1 , F 2 ) and F = |F|.Again, the 2D fractional Laplacian is most convenient defined by where |q| = q denotes the length of the wave vector q = (q 1 , q 2 ) that belongs to the Fourier transform In addition, for α = 1, the 2D fractional Laplacian can be represented as where R denotes the Riesz transform [21], that maps a scalar function f (ρ) into a vector field according to We note that relation Equation ( 37) is a generalization of the 1D differential operator Equation (22).Concerning the initial and boundary conditions in 2D, we have where ρ = |ρ|.The Schrödinger Equation (35) in momentum space becomes ∂ψ(q, t) ∂t where the Del operator has to be applied with regard to the Fourier variables.Similar as in 1D, we seek a wave function of the form ψ(q, t) = exp[iϕ(q)]u(q + Ft).Inserting this ansatz into Equation (40), leads to the first order partial differential equation In contrast to the 1D case, we are not able to solve this equation for all fractional orders in elementary form.However, for the relevant case of α = 1, an exact solution can be derived and summarized as The solution belonging to the classical case with α = 2 is less complicated and given by ϕ(q) = q 3  1 /(3F 1 ) + q 3 2 /(3F 2 ).Again, the unknown function u can be recovered under consideration of the initial condition such as Therefore, we obtain the wave function in momentum space according to ψ(q, t) = ψ 0 (q + Ft) exp[i(ϕ(q) − ϕ(q + Ft))] (44) In the case of the free particle with F = 0, the solution in momentum space can be given for all orders as ψ(q, t) = ψ 0 (q) exp(−itq α /2), 0 < α ≤ 2 (45) Moreover, if the initial wave packet exhibits rotational symmetry which implies that ψ 0 (q) = ψ 0 (q), the solution in real space can be evaluated by means of the inverse Hankel transform with J 0 (x) being the zero order Bessel function.For ψ 0 (q) = 1 and 1 < α ≤ 2, Equation (46) yields the 2D Green's function which can be evaluated by means of the series where G(0, t) = (−2i/t) 2/α Γ(2/α)/(2απ).In the case of α = 2, this relation coincides with the known result G(ρ, t) = exp(iρ 2 /(2t))/(2πit).On the other side, Green's function for α = 1 and ρ = t/2 is found to be Based on this relation, we obtain the wave function which is the general solution of Equation ( 35) for α = 1 and F = 0.In the following, we investigate the time evolution of the rotationally symmetric Cauchy wave packet which is normalized according to |ψ 0 (ρ)| 2 dρ = 1.Similar as in 1D, we consider the case of α = 1 in detail.The desired wave function can be obtained either via the convolution Equation ( 50) or by means of the inverse Hankel transform Equation ( 46).In the latter case, we need the initial wave packet Equation (51) in momentum space, that is ψ 0 (q) = 2ε √ 2π exp(−εq).The required inverse Hankel transform can be taken exactly and summarized as The associated probability density to this wave function is given by In accordance with the 1D case, the probability density in 2D also undergoes spreading and splitting.Remarkable, the probability density Equation (53) starts to split at t = 2ε, which is exactly the same time value as for the 1D pendant Equation (31). Figure 5 displays the time evolution of a rotationally symmetric Cauchy initial wave packet for different time values.We note that the solution for the rotationally symmetric Gaussian wave packet ψ 0 (ρ) = 2 exp(−2ρ 2 /σ 2 )/( √ πσ) can be derived in a similar manner.The required Hankel transform for evaluation of the wave function Equation ( 46) is given by ψ 0 (q) = √ πσ exp(−σ 2 q 2 /8).

Numerical Solution of the Fractional Schrödinger Equation
In this section, we describe a recently proposed matrix approach [22] that has been used for solving numerically the 1D space-fractional diffusion equation.In this context, we additionally refer to the recent publication [23], which takes into account the so-called Adomian decomposition method.
We adopt the matrix approach in order to solve the space-fractional Schrödinger Equation ( 1) subject to an arbitrary potential function V(x).Within this framework, we consider the computational domain |x| ≤ L/2 together with the uniform grid x j = −L/2 + jh for j = 0, . . ., N as well as the node spacing h = L/N.Furthermore, the fractional Laplacian has to be restricted to a bounded domain.If 0 < α < 1, we have the representation whereas for 1 < α < 2 it is given by In particular, the square root of the negative Laplacian on a bounded domain becomes For our applications, the quantity L has to be chosen relatively large in order to obtain a reasonable model for the infinite domain.On the other side, the matrix approach is designed for solving the fractional Schrödinger on a bounded domain, where analytical solutions are typically not available.
The key point here is the approximation of the fractional Laplacian by means of the fractional centered difference [22] (−∆) which can be applied for all 0 < α ≤ 2. The corresponding weights are given by [22] We note that within the definition Equation (57), we have taken into account the boundary data f (x 0 ) = f (x N ) = 0.The weights given in Equation (58) can be computed efficiently by means of the recursion [22] Applying the fractional centered difference scheme Equation (57) to the Schrödinger Equation (1) leads to the following system of ordinary differential equations where j = 1, . . ., N − 1, ψ j (t) = ψ(x j , t) and V j = V(x j ).The system Equation (60) in matrix notation becomes where Φ(t) = (ψ 1 (t), ψ 2 (t), . . ., ψ N−1 (t)).The associated solution can be directly given in terms of the matrix exponential according to where Φ 0 = (ψ 0 (x 1 ), ψ 0 (x 2 ), . . ., ψ 0 (x N−1 )) is a column vector that contains the sampled initial wave packet ψ 0 (x).If MATLAB is used, the matrix exponential can be evaluated via the function expm.For illustration purposes, the corresponding matrix M in the case of N = 5 looks like Now, the structure of the matrix M can be readily continued for an arbitrary N. In the following, we use the above described numerical scheme in order to solve the space-fractional Schrödinger Equation (1).We take into account the Cauchy initial wave packet Equation ( 24) with ε = 1 and evaluate the resulting probability density for α = 1 and t = 5.The size of the computational domain for approximation of the infinite domain is set to L = 30.Figure 6 displays the comparison between the analytical solution Equation (31) and the matrix approach Equation (62) for the case of a free particle.Next, we consider the case of the linear potential with β = 0.5.The resulting probability density belonging to this process is shown in Figure 7.In contrast to the free particle, the probability density is no longer symmetric and additionally exhibits some oscillations.(25) and the matrix approach Equation (62) for the case of a liner potential with β = 0.5.
We repeat the last numerical experiment for the case of a stronger accelerated particle with β = 4.The resulting probability density, which is depicted in Figure 8, is given by a smooth curve without oscillations.Furthermore, in contrast to the free particle case, the probability density has not been splitted.The numerical experiments outlined above confirm both the correctness of the analytical solutions as well as the good quality of the matrix approach.Besides these comparisons, we additionally have performed numerical experiments for other fractional orders α.As a result, we observed in all cases the same good agreement between the analytical solution Equation ( 10) and the matrix approach Equation (62).It should also be noted that the above matrix approach can be readily extended for solving the space-time fractional Schrödinger equation.More precisely, if we replace the classical time derivative occurring in Equation ( 1) by the Caputo derivative [24] of order 0 < ν ≤ 1 one obtains the following system of fractional differential equations The solution to this system is formally given by Φ is the Mittag-Leffler matrix function.It can be computed under consideration of the eigenvalue decomposition A = UΛU −1 , where U is a matrix containing the eigenvectors of A. For A ∈ C n×n with eigenvalues Λ = diag(λ 1 , λ 2 , . . ., λ n ), we have E ν (A) = UE ν (Λ)U −1 together with the diagonal matrix As a final numerical experiment, we consider the space-time fractional Schrödinger equation for the joint fractional order α ∈ (0, 1).In this case, the wave function in momentum space can be written as ψ(ω, t) = ψ 0 (ω)G(ω, t), where G(ω, t) = E α (−i|ω| α t α /2).The corresponding Green's function in space domain is available in elementary form as We note that this result can be derived in a similar manner as the fundamental solution to the neutral-fractional wave equation [25].The response to an arbitrary initial stage is given by the convolution ψ(x, t) = it α sin(απ/2) 2απ ∞ 0 ψ 0 (x − u 1/α ) + ψ 0 (x + u 1/α ) u 2 + iut α cos(απ/2) − t 2α /4 du (69) In most cases, this integral must be evaluated numerically.However, for the Cauchy initial wave packet Equation ( 24) the integration can be carried out analytically and summarized as ψ(x, t) = ε 2π In Figure 9 we have evaluated the solution of the space-time fractional Schrödinger equation for two different fractional orders subject to the initial stage Equation (20).The analytical solution Equation ( 69) is denotes by filled dots, whereas the solid line is the numerical solution obtained by means of the Mittag-Leffler matrix function.In both cases the wave function is evaluated for the time value t = 5.

Conclusions and Discussion
In this article, we considered the space-fractional Schrödinger equation with the Riesz space-fractional derivative in the presence of the linear potential.The general solution to the 1D Schrödinger equation in momentum space has been obtained in closed form, which allows the determination of other measurable quantities such as the spatial moments.Furthermore, we investigated the time evolution of a Cauchy wave packet in detail.Besides the 1D Schrödinger equation we also addressed the case of two spatial dimensions.The obtained equations have been illustrated and verified by comparisons with an independent recently proposed numerical scheme.Consequently, the derived equations are therefore also useful for verification of other numerical approaches which are designed for solving the fractional Schrödinger equation.Besides the space-fractional time-dependent Schrödinger equation solutions to the associated time-independent counterpart are also of importance.Concerning the relationship between the space-fractional time-dependent and time-independent Schrödinger equation we refer to the article [26].

Figure 1 .
Figure 1.Illustration of the mean square displacement (MSD) Equation (21) for different values of the fractional order α under consideration of the Gaussian initial wave packet Equation (20) with σ = 1.

Figure 2 .
Figure 2. Time evolution of the Cauchy initial wave packet Equation (24) for the parameter ε = 1.

Figure 3 .
Figure 3.Time evolution of the Gaussian initial wave packet Equation (20) for the parameter σ = 1.

Figure 4 .
Figure 4. Time evolution of the Cauchy initial wave packet Equation (24) for ε = 1 within the classical quantum mechanics.

Figure 5 .
Figure 5.Time evolution of the rotationally symmetric Cauchy initial wave packet Equation (51) for the parameter ε = 1.

Figure 6 .
Figure 6.Comparison between the analytical solution Equation (31) and the matrix approach Equation (62) for the case of a free particle.

Figure 7 .
Figure 7.Comparison between the analytical solution Equation(25) and the matrix approach Equation (62) for the case of a liner potential with β = 0.5.

Figure 8 .
Figure 8.Comparison between the analytical solution Equation(25) and the matrix approach Equation (62) for the case of a liner potential with β = 4.

Figure 9 .
Figure 9.Comparison between the analytical solution Equation (69) and the matrix approach for the case of the Gaussian initial wave packet Equation (20) with σ = 0.5.