An Analytical Study for (2 + 1)-Dimensional Schrödinger Equation

In this paper, the homotopy analysis method has been applied to solve (2 + 1)-dimensional Schrödinger equations. The validity of this method has successfully been accomplished by applying it to find the solution of some of its variety forms. The results obtained by homotopy analysis method have been compared with those of exact solutions. The main objective is to propose alternative methods of finding a solution, which do not require small parameters and avoid linearization and physically unrealistic assumptions. The results show that the solution of homotopy analysis method is in a good agreement with the exact solution.


Introduction
Stationary and time-dependent Schrödinger equation formulated by the Austrian physicist Erwin Schrödinger plays a fundamental role in physics for describing quantum mechanical [1].
In recent years, there exists a considerable number of works dealing with the problem of approximate solutions of the Schrödinger equation by using different methodologies, for example, Adomian decomposition method (ADM) [2], homotopy perturbation method (HPM) [3,4], homotopy analysis method (HAM) [5], the boundary value method [6], and variational iteration method [7].
In this paper, an analytical method called the homotopy analysis method (HAM) will be used to solve (1).
The homotopy analysis method (HAM) was first proposed by Liao [8] in 1992, to get analytic approximations of highly nonlinear equations. Unlike other existing methods, this method is independent of small/large physical parameters, provides us a simple way to ensure the convergence of solution series, and gives us the great freedom to choose proper base functions.
These advantages make the method a powerful and flexible tool in mathematics and engineering, which can be readily distinguished from existing numerical and analytical methods [9,10]. Recently, considerable research have been conducted in applying this method to a class of linear and nonlinear equations [10][11][12][13][14][15].
This paper is arranged in the following manner. In Section 2 the basic idea of standard HAM is illustrated. In Section 3, the implementation of this method on some examples is presented. Finally, conclusions are drawn in Section 4.

Standard Homotopy Analysis Method
Let us consider the differential equation where is a differential operator , and denote independent variables, and ( , , ) is an unknown function. Based on the constructed zero-order deformation equation by Liao [8], we give the following zero-order deformation equation in a similar way: where ℎ ̸ = 0 denotes an auxiliary parameter, ∈ [0, 1] is an embedding parameter, is an auxiliary linear integer-order operator, 0 ( , , ) is an initial guess of unknown function ( , , ), and ( , , ; ) is a kind of mapping, as described later. It is important that one has great freedom to choose auxiliary parameter ℎ in homotopy analysis method. If = 0 and = 1, it holds that Thus as increases from 0 to 1, the solution ( , , ; ) varies from the initial guess 0 ( , , ) to the solution ( , , ). Expanding ( , , ; ) in Taylor series with respect to , one has where If the auxiliary linear operator, the initial guess, and the auxiliary parameter ℎ are so properly chosen, the series (5) converges at = 1, and one has According to the above, the governing equation can be deduced from the zero-order deformation, (3).
Define the vector Differentiation equation (3) times with respect to the embedding parameter and then setting = 0 and finally dividing them by !, we have the so-called th-order deformation equation where The th-order deformation equation (9) is linear and thus can be easily solved, especially by means of symbolic computation software such as Maple.
The th-order component ( , , )( ≥ 1) would be achieved by means of symbolic computation software Maple, Mathematica, and so on.
We still have freedom to choose the auxiliary parameter ℎ. To investigate the influence of ℎ on the solution series, one can consider the convergence of approximation series related to a point in a domain [16]. These curves contain a horizontal line segment. This horizontal line segment denotes the valid region of ℎ which guaranteed the convergence of the related series.

Numerical Examples
We will illustrate the accuracy and efficiency of the homotopy analysis method applied to (1). For this purpose, here we present some numerical examples as considered in [4,6].
The exact solution is given by Starting with 0 ( , , ) = ( , , 0) and by using (15), we now successively obtain by HAM and so forth. The proper value of ℎ = −1 is found from the ℎ-curve shown in Figure 1. Then the series solution expression is obtained by HAM as which clearly converges to the exact solution (18).

Example 2.
Let us have ( , ) = 3 − 2tanh 2 ( ) − 2tanh 2 ( ) in (1) using the following initial condition: The exact solution is In this manner, the rest of components of the standard HAM solution can be found.
Again, the value ℎ = −1 was chosen based on the ℎcurve shown in Figure 2. Then the series solution expression is obtained by HAM as which coincides with the exact solutions (22). Example 3. As another example, let us consider (1) using ( , ) = 1 − 2/ 2 − 2/ 2 , which has the following exact solution: We will solve this example directly by using HAM. We choose the initial approximation Plotting the ℎ-curve similar to what was plotted in Figure 2 suggests that we can take ℎ = −1. Then the series solution expression is obtained by HAM as which converges to the exact solution (25).

Conclusion
In this paper, we have successfully developed homotopy analysis method to obtain the exact solutions of (2 + 1)dimensional Schrödinger equation. It is apparently seen that these method are very powerful and efficient for solving different kinds of problems arising in various fields of science and engineering and present a rapid convergence for the solutions. Mohebbi and Dehghan in [6] reported the computed error for Examples 1-4, and in the present work we have obtained the exact solutions.