The Scientific World Journal

Volume 2014 (2014), Article ID 438345, 5 pages

http://dx.doi.org/10.1155/2014/438345

Research Article

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

Behzad Ghanbari

Department of Basic Sciences, Kermanshah University of Technology, Kermanshah, Iran

Received 30 August 2013; Accepted 30 October 2013; Published 27 January 2014

Academic Editors: N. Damil and X.-w. Gao

Copyright © 2014 Behzad Ghanbari. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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.

1. 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].

We consider the linear ()-dimensional Schrödinger equation with variable coefficients of the form [6]

where is an arbitrary potential function and .

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–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.

2. 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 denotes an auxiliary parameter, is an embedding parameter, is an auxiliary linear integer-order operator, 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 and , it holds that

Thus as increases from 0 to 1, the solution varies from the initial guess 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 , 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 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.

3. Application

To solve (1) by means of the standard HAM, we choose

According to (10)-(11) and (12), one has

where

Substituting (11) and (14) into (13), we have

where the constant of integration is determined by the initial conditions

The th-order component 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 . These curves contain a horizontal line segment. This horizontal line segment denotes the valid region of which guaranteed the convergence of the related series.

4. 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].

Example 1. Consider (1) with and the following initial condition:

The exact solution is given by

Starting with and by using (15), we now successively obtain by HAM

and so forth.

The proper value of 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).

Figure 1: The -curves for the 7th-order of HAM approximation of for Example 1; dotted line: real part of approximation; solid line: imaginary part of approximation.

Example 2. Let us have in (1) using the following initial condition:

The exact solution is

Starting with in recursive scheme (15), the following components are obtained:

In this manner, the rest of components of the standard HAM solution can be found.

Again, the value 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).

Figure 2: The -curves for the 7th-order of HAM approximation of for Example 2; dotted line: real part of approximation; solid line: imaginary part of approximation.

Example 3. As another example, let us consider (1) using , which has the following exact solution:

We will solve this example directly by using HAM. We choose the initial approximation

Using such starting with , the following solutions are obtained:

Plotting the -curve similar to what was plotted in Figure 2 suggests that we can take . Then the series solution expression is obtained by HAM as

which converges to the exact solution (25).

Example 4. As the last example, let us try to solve (1) with and the following initial condition:

Its exact solution reads

Starting with in HAM procedure, we successively obtain

and so forth. The value (which can be obtained by plotting the same -curve as was plotted in Figure 2) yields

which converge to the exact solutions (30).

5. Conclusion

In this paper, we have successfully developed homotopy analysis method to obtain the exact solutions of ()-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.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

A. Arnold, “Numerically absorbing boundary conditions for quantum evolution equations,” VLSI Design, vol. 6, no. 1–4, pp. 313–319, 1998.

S. A. Khuri, “A new approach to the cubic Schrodinger equation: an application of the decomposition technique,” Applied Mathematics and Computation, vol. 97, pp. 251–254, 1998.

J. Biazar and H. Ghazvini, “Exact solutions for non-linear Schrödinger equations by He's homotopy perturbation method,” Physics Letters A, vol. 366, no. 1-2, pp. 79–84, 2007.

J. Biazar and B. Ghanbari, “The homotopy perturbation method for solving neutral functional-differential equations with proportional delays,” Journal of King Saud University-Science, vol. 24, no. 1, pp. 33–37, 2012.

A. K. Alomari, M. S. M. Noorani, and R. Nazar, “Explicit series solutions of some linear and nonlinear Schrodinger equations via the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1196–1207, 2009.

A. Mohebbi and M. Dehghan, “The use of compact boundary value method for the solution of two-dimensional Schrödinger equation,” Journal of Computational and Applied Mathematics, vol. 225, no. 1, pp. 124–134, 2009.

N. H. Sweilam and R. F. Al-Bar, “Variational iteration method for coupled nonlinear Schrödinger equations,” Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 993–999, 2007.

S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.

S. J. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.

S. J. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman & Hall, CRC Press, Boca Raton, Fla, USA, 2003.

S. J. Liao, “Notes on the homotopy analysis method: some definitions and theorems,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 983–997, 2009.

S. J. Liao, “A new branch of solutions of boundary-layer flows over an impermeable stretched plate,” International Journal of Heat and Mass Transfer, vol. 48, no. 12, pp. 2529–2539, 2005.

B. Ghanbari, L. Rada, and K. Chen, “A restarted iterative homotopy analysis method for twononlinear models from image processing,” International Journal of Computer Mathematics, 2013.

J. Biazar and B. Ghanbari, “HAM solution of some initial value problems arising in heat radiation equations,” Journal of King Saud University-Science, vol. 24, no. 2, pp. 161–165, 2012.

S. Liang and D. J. Jeffrey, “Comparison of homotopy analysis method and homotopy perturbation method through an evolution equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 12, pp. 4057–4064, 2009.