Abstract
We construct the approximate solutions of the time-fractional Schrödinger equations, with zero and nonzero trapping potential, by homotopy analysis method (HAM). The fractional derivatives, in the Caputo sense, are used. The method is capable of reducing the size of calculations and handles nonlinear-coupled equations in a direct manner. The results show that HAM is more promising, convenient, efficient and less computational than differential transform method (DTM), and easy to apply in spaces of higher dimensions as well.
1. Introduction
The theory of derivatives of fractional (nonintegers) orders stimulates considerable interest in the areas of mathematics, physics, engineering, and other sciences. Fractional derivatives [1–5] provide an excellent tool for the description of memory and hereditary properties of various material and processes. The beauty of this subject is that a fractional derivative is not a local point property. This considers the history and nonlocal distributed effects. Perhaps, this subject translates the reality of nature better. Application of fractional calculus are found in different areas of sciences such as physics, continuum mechanics, signal processing, electromagnetics, and bioengineering. The electrical properties of nerve cell membranes and the propagation of electrical signals are well characterized by differential equations of fractional order. The fractional differential equations (FDE) [6–11] appear more and more frequently in different research areas and engineering applications. Exact solution of nonlinear partial differential equations has become one of the central themes of perpetual interest in mathematical physics. In order to better understand these phenomena as well as further apply them in the practical life, it is important to seek their exact solutions.
A great deal of effort has also been expanded in attempting to find robust and stable numerical and analytical methods for solving fractional differential equations of physical interest. These methods include Laplace transform method, Fourier transform method, finite difference method [12], fractional linear multistep methods, Adomian decomposition method (ADM) [13], variational iteration method (VIM) [14], homotopy perturbation method (HPM) [15], differential transform method [16], and homotopy analysis method (HAM) [17–20].
The aim of this paper is to solve the Schrödinger equation with fractional order using the homotopy analysis method. By introducing an embedding parameter the nonlinear fractional differential equation is converted to a linear fractional differential equation at . When evolves, the differential equation becomes the original one at . The method has been used in a variety of problems and the details can be found in Liao’s book [17]. The method gives rapidly convergent successive approximations of the exact solution if such solution exists; otherwise, a few approximations can be used for numerical purposes. Wang [12] presented the numerical solution Schrödinger equations by means of finite difference scheme. Khuri [13] applied ADM to obtain the solution of cubic Schrödinger equations. Wazwaz [14] presented the exact solution of the linear and nonlinear one-dimensional Schrödinger equations by VIM. Recently, Ravi Kanth [16] and his colleague presented the exact solution of the linear and nonlinear Schrödinger equations by differential transformation method (DTM).
The aim of this paper is to investigate the approximate solutions of the time-fractional Schrödinger equations, with zero and non-zero trapping potential, by means of HAM. The convergent region is then obtained by looking at the real and imaginary parts of the series in plot.
2. Fractional Schrödinger Equations and Preliminaries
The time-fractional Schrödinger equations (FSE) has the following form: where is the trapping potential and is a real constant. The physical model (2.1) and its generalized forms occur in various areas of physics, including nonlinear optics, plasma physics, superconductivity, and quantum mechanics.
We give some basic definitions, notations, and properties of the fractional calculus theory, which will be used later in this work.
Definition 2.1. The Riemann-Liouville fractional integral operator of order on the usual Lebesgue space is given by This integral operator has the following properties:(i),(ii).
Definition 2.2. The Caputo definition of fractal derivative operator is given by It has the following two basic properties: and The Caputo fractional derivative is considered here because it allows traditional initial and boundary conditions to be included in the formulation of the problem.
3. Basic Idea of HAM
In this paper, we apply the HAM to the linear and nonlinear problems to be discussed. In order to show the basic idea of HAM, consider the following nonlinear fractional differential equation:
where is a nonlinear fractional operator, and denote the independent variables, and is an unknown function. By means of the HAM, we first construct the so-called zeroth-order deformation equation
where is the embedding parameter, is an auxiliary parameter, is an auxiliary linear operator, are unknown functions, and are initial guesses of . It is obvious that for and , (3.3) becomes
respectively. Thus as increases from 0 to 1, the solution varies from the initial guess to the exact solution . Expanding in Taylor series with respect to , one has
where
The convergence of the series (3.4) depends upon the auxiliary parameter . If it is convergent at , one has
which must be one of the solutions of the original nonlinear equations, as proved by Liao [17]. Define the vectors Differentiating the zeroth-order deformation (3.2) -times with respect to and then dividing them by , and finally setting , we get the following th-order deformation equation:
where
It should be emphasized that for is governed by the linear equations (3.8) with boundary conditions that come from the original problem, which can be solved by the symbolic computation software MATHEMATICA. The success of the technique is based on the proper selection of the initial guess.
4. Implementation of the Method
In this section, we introduce the above reliable approach, in a realistic and efficient way, to handle nonlinear Schrödinger equation with time-fractional derivative.
4.1. Nonlinear-Time Fractional Schrödinger Equation (NLFSE)
Setting and in (2.1) leads to the following coupled system of equations
subject to the initial conditions
Equations (4.1) and (4.3) suggest that we define the nonlinear operator as
and the linear operator
where
Obviously, the solution of the th-order deformation (4.8) for m ≥ 1 becomes
Example 4.1. We first consider the one-dimensional NLFSE with zero trapping potential (i.e.,) and
subject to the initial condition .
Solving the above equations
and so on, in this manner the rest of the components can be obtained. Therefore, the approximate solution is
The exact solution of (4.10) for is . When , , the solution obtained by [13–16] is recovered as a special case.
Example 4.2. Consider the one-dimensional NLFSE with trapping potential, that is, and
subject to the initial condition .
Solving above equations we obtain the solution in a series form
The exact solution of (4.13) for is . When , , the solution obtained by [13–16] is recovered as a special case.
Example 4.3. Consider the two dimensional NLFSE with trapping potential Here , , subject to the initial condition . Solving the above equation we obtain the approximate solution in a series form The exact solution of (4.15) for , is .
Example 4.4. Consider the three-dimensional NLFSE with trapping potential where subject to the initial condition Solving the above equations, Finally the approximate solution in the series form is The exact solution of (4.17) for is .
5. Closing Comments
The basic goal of this work has been to extend the works made in the nonlinear physical problem to construct solutions for nonlinear Schrödinger equation with time-fractional derivatives. The goal has been achieved and new solutions of have been derived for nonlinear equations with time fractional derivatives. The proposed approach works successfully in handling nonlinear fractional Schrödinger equations directly with a minimum size of calculations.
Borhanifer and Abazari [11] applied the differential transformation method (DTM) to solving Schrödinger and coupled Schrödinger equations. The major lacks of DTM are that it requires transformation, and the given differential equation and related initial conditions are transformed into a recurrence equation that finally leads to the solution of a system of algebraic equations as coefficients of a power series solution. The main disadvantage of DTM is that it requires transformation, which will be complicated and computational cost will be too much. This emphasizes the fact that the presented approach can be used in a wider class of system of nonlinear fractional differential equations. HAM is a powerful and efficient technique in finding exact and approximate solutions for linear and nonlinear models. HAM provides more realistic solutions that converge very rapidly in real physical problems. The numerical examples show that the solutions are in good agreement with their respective exact solutions for .
In the last, we present the -curves to see the convergent region of as in Figures 1–4. We plot the imaginary and real part of each example in the same figure, and the convergent region is the region of intersection between the convergent regions of the imaginary and real parts. In all figures is in the convergent region.
Acknowledgments
The author N. A. Khan is thankful and grateful to the dean of the Faculty of Sciences, University of Karachi, Pakistan for supporting and facilitating this research work. The author M. Jamil is highly thankful and grateful to the Abdul Salam School of Mathematical Sciences, GC University, Lahore, Pakistan, Department of Mathematics, NED University of Engineering and Technology, Pakistan, and also Higher Education Commission of Pakistan for generous support and facilitating this research work.