Advances in Mathematical Physics

Volume 2013, Article ID 239254, 8 pages

http://dx.doi.org/10.1155/2013/239254

## Application of the Asymptotic Taylor Expansion Method to Bistable Potentials

Department of Engineering Physics, Faculty of Engineering, University of Gaziantep, 27310 Gaziantep, Turkey

Received 15 April 2013; Accepted 15 July 2013

Academic Editor: Fedele Lizzi

Copyright © 2013 Okan Ozer et al. 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

A recent method called asymptotic Taylor expansion (ATEM) is applied to determine the analytical expression for eigenfunctions and numerical results for eigenvalues of the Schrödinger equation for the bistable potentials. Optimal truncation of the Taylor series gives a best possible analytical expression for eigenfunctions and numerical results for eigenvalues. It is shown that the results are obtained by a simple algorithm constructed for a computer system using symbolic or numerical calculation. It is observed that ATEM produces excellent results consistent with the existing literature.

#### 1. Introduction

There is no doubt that an interesting problem in fundamental quantum mechanics for lecturers and advanced undergraduate and graduate students in physics and applied mathematics is to obtain the exact solutions of the Schrödinger equation for any type of potential. It is well known that the Schrödinger equation, proposed by Erwin Schrödinger in 1926, is a second-order differential equation that describes how the quantum state of a physical system changes with time [1]. It is as central to quantum mechanics as Newton's laws are to classical mechanics. It is also well known that the quantum mechanics is established on some certain postulates, and in any introductory quantum physics textbook these postulates can be found with the application of the time-independent Schrödinger equation for a single particle in one dimension as a simple and basic example.

Unfortunately, there are not so many potentials that can be solved exactly such as the Coulomb, harmonic oscillator, and Pöschl-Teller potentials. Since one of the source of progress of the science depends on the study of the same problem from different point of view, various methods have been suggested such as numerical calculation [2, 3], the variational [2, 4], the perturbation [5], the WKB [6, 7], the shifted expansion [8, 9], the Nikiforov-Uvarov (NU) [10, 11], the supersymmetry (SUSY) [12, 13], the generalized pseudospectral (GPS) [14], the asymptotic iteration method (AIM) [15], and other methods [16–21] to find the approximate solutions of the potentials that are not exactly solvable.

In this study we will apply a new formalism based on the Taylor series expansion method, namely, asymptotic Taylor expansion method (ATEM) [22], to bistable potentials. These type of potentials have been used in the quantum theory of molecules as a crude model to describe the motion of a particle in the presence of two centers of force [23–29]. It is mentioned in [22] that the taylor series Method [30, 31] is an old one but it has not been fully exploited in the analysis of both physical and mathematical problems in solution. It is also claimed that ATEM can also be easily applied to solve second-order differential equations by introducing a simple Mathematica [32] computer program. Therefore, we focus on the solution of the eigenvalue problems of some type of bistable potentials by using ATEM in this paper.

The organization of the paper is as follows. In Section 2, we present a brief outline of the method ATEM. In Section 3, the eigenvalues of the bistable potentials are determined by using ATEM. Section 4 is devoted to a conclusion.

#### 2. Formalism of Asymptotic Taylor Expansion Method

In this section, we present the solution of the Schrödinger-type equations by modifying Taylor series expansion with the aid of a finite sequence instead of an infinite sequence and its termination possessing the property of quantum mechanical wave function. It is well known in quantum mechanics that the bound-state energy of an atom is quantized and eigenvalues are discrete. Additionally, for each eigenvalues there exist one or more eigenfunctions. If one considers the solution of the Schrödinger equation, discrete eigenvalues of the problem are mainly under investigation. The first main result of this conclusion gives necessary and sufficient conditions for the termination of the Taylor series expansion of the wave function.

Following the notation in [22], one can consider the Taylor series expansion of a function about the point : where is the derivative of the function at . Taylor series specifies the value of a function at one point, , in terms of the value of the function and its derivatives at a reference point . Expansion of the function about the origin () is known as Maclaurin's series and it is given by Here one can develop a method to solve a second-order linear differential equation of the form It is seen that the higher-order derivatives of the can be obtained in terms of the and by differentiating (3). Then, higher-order derivatives of are given by where At this point, one can observe that the eigenvalues and eigenfunctions of the Schrödinger-type equations can efficiently be determined by using ATEM. To this end, the recurrence relations (5) allow one to determine algebraic exact or approximate analytical solution of (3) under some certain conditions. Let us now substitute (5) into (1) to obtain the function that is related to the wave function of the corresponding Hamiltonian: After all, one can now obtain useful formalism of the Taylor expansion method. This form of the Taylor series can also be used to obtain series solution of the second-order differential equations. In the solution of the eigenvalue problems, truncation of the asymptotic expansion to a finite number of terms is useful. If the series optimally truncated at the smallest term then the asymptotic expansion of series is known as superasymptotic [33], and it leads to the determination of eigenvalues with minimum error.

Since the improper sets of boundary conditions may produce nonphysical results, arrangement of the boundary conditions for different problems becomes very important. When only odd or even power of is collected as coefficients of or and vice versa, the series is truncated at ; then an immediate practical consequence of these conditions is obtained for or . In this way, the series truncates at and one of the parameters in the or belongs to the spectrum of the Schrödinger equation. Therefore eigenfunction of the equation becomes a polynomial of degree . Otherwise the spectrum of the system can be obtained as follows. In a quantum mechanical system eigenfunction of the system is discrete. Therefore in order to terminate the eigenfunction we can concisely write that Eliminating and we obtain again one of the parameters in the equation related to the eigenvalues of the problem.

It has been stated that the ATEM reproduces exact solutions to many exactly solvable differential equations and these equations can be related to the Schrödinger equation. It is observed that the process presented in ATEM is iterative and the number of iteration is given by . The method can be applied to the Schrödinger equation with any type of potential as follows. Using a computer program, one first sets up the iteration number , say , to obtain the result. Then, setting , another result is obtained. This procedure is repeated for different values leading to different results. Finally, one can compare the results for each case till desired digits. If the values of the eigenvalue reach their asymptotic values, then one can choose the corresponding value and truncate the iteration for next calculations. For instance, if one can obtain the values of the eigenvalues for , first few of them—first eight eigenvalues, for example—will reach automatically their asymptotic values. The following comment on the function is considerable. For such a solution it may be suitable to take sum of the first eight terms in (6).

#### 3. Applications

We shall illustrate here that (3) with conditions (7) and (8) gives a complete solution for some important Schrödinger-type problems. Through a concrete example we explore the solution of Schrödinger equation (): for the harmonic oscillator potential in one dimension given as In the limit of large , the asymptotic solutions of (10) can be taken as any power of times a decreasing Gaussian function to satisfy the quantum mechanics postulates. With this in mind, one can write the “unnormalized” wave functions as where the functions are to be found by means of the iteration procedure given previously. Substituting (11) into (10), one obtains Comparing (3) and (12) one can deduce that Following the procedure given in (8) yields the exact eigenvalues of the harmonic oscillator potential: and the normalized eigenfunctions of the Schrödinger equation in (10) are shown in Figure 1.

As a second application, we consider the quartic anharmonic oscillator potential that has been a great deal of interest in the analytical and numerical investigations of the one-dimensional anharmonic oscillators because of their importance in molecular vibrations [34] as well as in solid state physics [35, 36] and quantum field theories [37]. Schrödinger equation (9) for the one-dimensional quartic anharmonic potential is written as where . The solution of (15) has been always studied to test accuracy and efficiency of the different methods proposed in [15, 38–41]. Now, we introduce the asymptotic solutions of (15) as and (15) can now be written as Comparing (3) and (17) one can deduce that By the aid of computer program, one can calculate the eigenvalues and the corresponding eigenfunctions for a range of values, changing from to , using number of iterations . The term “asymptotic” means the function approaching to a given value as the iteration number tends to infinity.

We present our results carried out for a range of values in Table 1 with significant digits, and they are compared with those of supersymmetric perturbation approach by [40] and the ones computed numerically by [42]. In our calculations, we set , , and . It is observed that there is remarkable agreement in the whole range of values for all quantum states for different values with results of [42] except .

We also present and compare our results for with those of Bacus et al. [39] in Table 2 with 20 significant digits. The function , for , state for different values of is found as follow:
For the first six states, the plot of the normalized wave functions for are given in Figure 2. As a last example, we consider the symmetric bistable potential defined in [43] as
Inserting the potential function into Schrödinger equation in (9) and using the *ansatz* wave function of the form defined in (16), one obtains
Following the procedure given in (8) for , , and yields the eigenvalues of the symmetric bistable potential in (20). We note here the selection of adjustable parameters and . Since we choose the *ansatz* wave function given in (16), the form of the superpotential for the partner potentials by supersymmetric quantum mechanics (SUSYQM) [44] is found to be
and the partner potentials are defined as
Since the superpotential obtained from (16) is
then the partner potential is obtained as
If , the potential (25) has one minimum. Investigating this one-minimum case one can select the adjustable parameters and that lead to an increase in the iteration number . Therefore, we select and set the parameters as and in our calculations.

Our results are presented and compared with the values found by the variational supersymmetric method () [43], state-dependent diagonalization method () [45], and by direct numerical integration () [43], in Table 3. The ATEM results are in a very good agreement, by low percent errors, for all values of energies with the ones obtained by numerical calculation. For lower percent errors, the iteration number must be increased.

#### 4. Conclusion

An approximate method based on the asymptotic Taylor series expansion of a function and its fundamental features are presented. It is observed that the method is applicable for obtaining both eigenvalues and eigenfunctions of the Schrödinger-type equations. After applying the method to the one-dimensional harmonic oscillator potential, it is shown that the approach gives accurate results for eigenvalue problems of some certain type of bistable potentials. It is thought that the approach opens the way to the treatment of the Schrödinger equation including large class of potentials of practical interest. As a future study, the method can be developed and applied to the non-Hermitian systems and QES potentials. The direct application of the method to the Fokker-Planck equation for quasi exactly solvable bistable potentials when the drift coefficient has a nonpolynomial nature would also be interesting. Studies along this line are in progress.

#### Acknowledgments

This work is supported by the Research Fund Unit of Gaziantep University. Okan Ozer would like to thank Dr. G. Lévai (ATOMKI, Debrecen-Hungary) for supporting some papers used in this study, and it is a pleasure to thank Dr. R. Koc for valuable conversations and comments. The authors are also very much indebted to the referee(s) for useful comments and suggestions.

#### References

- E. Schrödinger, “Quantisierung als Eigenwertproblem,”
*Annalen der Physik*, vol. 79, no. 4, pp. 361–376, 1926. View at Publisher · View at Google Scholar - Y. P. Varshni, “Eigenenergies and oscillator strengths for the Hulthén potential,”
*Physical Review A*, vol. 41, no. 9, pp. 4682–4689, 1990. View at Publisher · View at Google Scholar · View at Scopus - M. A. Nunez, “Accurate computation of eigenfunctions for Schrödinger operators associated with Coulomb-type potentials,”
*Physical Review A*, vol. 47, no. 5, pp. 3620–3631, 1993. View at Publisher · View at Google Scholar - C. Stubbins, “Bound states of the Hulthén and Yukawa potentials,”
*Physical Review A*, vol. 48, no. 1, pp. 220–227, 1993. View at Publisher · View at Google Scholar - P. Matthys and H. de Meyer, “Dynamical-group approach to the Hulthén potential,”
*Physical Review A*, vol. 38, no. 3, pp. 1168–1171, 1988. View at Publisher · View at Google Scholar · View at Scopus - Y. P. Varshni, “Relative convergences of the WKB and SWKB approximations,”
*Journal of Physics A*, vol. 25, no. 21, pp. 5761–5777, 1992. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. A. Dobrovolsky and R. S. Tutik, “Regularization of the WKB integrals,”
*Journal of Physics A*, vol. 33, no. 37, pp. 6593–6599, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Z. Tang and F. T. Chan, “Shifted
*1/N*expansion for the Hulthén potential,”*Physical Review A*, vol. 35, no. 2, pp. 911–814, 1987. View at Publisher · View at Google Scholar - R. K. Roychoudhury and Y. P. Varshni, “Shifted
*1/N*expansion and exact solutions for the potential $V(r)=-Z/r+gr+{\lambda r}^{2}$,”*Journal of Physics A*, vol. 21, no. 13, pp. 3025–3034, 1988. View at Publisher · View at Google Scholar · View at MathSciNet - A. F. Nikiforov and V. B. Uvarov,
*Special Functions of Mathematical Physics*, Birkhäuser, Basel, Switzerland, 1988. View at MathSciNet - B. Gönül and K. Köksal, “A search on the Nikiforov-Uvarov formalism,”
*Physica Scripta*, vol. 75, no. 5, pp. 686–690, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Gönül, O. Özer, Y. Cançelik, and M. Koçak, “Hamiltonian hierarchy and the Hulthén potential,”
*Physics Letters A*, vol. 275, no. 4, pp. 238–243, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. W. Qian, B. W. Huang, and Z. Y. Gu, “Supersymmetry and shape invariance of the effective screened potential,”
*New Journal of Physics*, vol. 4, article 13, 2002. View at Publisher · View at Google Scholar - A. K. Roy, “The generalized pseudospectral approach to the bound states of the Hulthén and the Yukawa potentials,”
*Pramana Journal of Physics*, vol. 65, no. 1, pp. 1–15, 2005. View at Publisher · View at Google Scholar - H. Ciftci, R. L. Hall, and N. Saad, “Asymptotic iteration method for eigenvalue problems,”
*Journal of Physics A*, vol. 36, no. 47, pp. 11807–11816, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. S. Lai and W. C. Lin, “Energies of the Hulthén potential for
*l**≠*0,”*Physics Letters A*, vol. 78, no. 4, pp. 335–337, 1980. View at Publisher · View at Google Scholar - S. H. Patil, “Energy levels of screened Coulomb and Hulthen potentials,”
*Journal of Physics A*, vol. 17, no. 3, article 575, 1984. View at Publisher · View at Google Scholar - B. Roy and R. Roychoudhury, “The shifted
*1/N*expansion and the energy eigenvalues of the Hulthen potential for l not = 0,”*Journal of Physics A*, vol. 20, no. 10, p. 3051, 1987. View at Publisher · View at Google Scholar - R. L. Greene and C. Aldrich, “Variational wave functions for a screened Coulomb potential,”
*Physical Review A*, vol. 14, no. 6, pp. 2363–2366, 1976. View at Publisher · View at Google Scholar · View at Scopus - U. Myhrman, “A recurrence formula for obtaining certain matrix elements in the base of eigenfunctions of the Hamiltonian for a particular screened potential,”
*Journal of Physics A*, vol. 16, no. 2, pp. 263–270, 1983. View at Publisher · View at Google Scholar · View at MathSciNet - A. Bechlert and W. Bhring, “Analytic model of atomic screening based on a generalised Hulthen potential,”
*Journal of Physics B*, vol. 21, no. 5, p. 817, 1988. View at Publisher · View at Google Scholar - R. Koç and S. Sayın, “Remarks on the solution of the position-dependent mass Schrödinger equation,”
*Journal of Physics A*, vol. 43, no. 45, Article ID 455203, 8 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Razavy, “An exactly soluble Schrödinger equation with a bistable potential,”
*American Journal of Physics*, vol. 48, no. 4, p. 285, 1980. View at Publisher · View at Google Scholar - Q.-T. Xie, “New quasi-exactly solvable double-well potentials,”
*Journal of Physics A*, vol. 45, no. 17, Article ID 175302, 7 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. U. Felderhof, “Diffusion in a bistable potential,”
*Physica A*, vol. 387, no. 21, pp. 5017–5023, 2008. View at Publisher · View at Google Scholar · View at Scopus - H. Ciftci, O. Ozer, and P. Roy, “Asymptotic iteration approach to supersymmetric bistable potentials,”
*Chinese Physics B*, vol. 21, no. 1, Article ID 010303, 2012. View at Publisher · View at Google Scholar - T. Barakat, “The asymptotic iteration method for the eigenenergies of the anharmonic oscillator potential $V(x)=Ax2\alpha +{Bx}^{2}$,”
*Physics Letters A*, vol. 344, no. 6, pp. 411–417, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. J. Sous, “Solution for the eigenenergies of sextic anharmonic oscillator potential $V(x)={A}_{6}{x}^{6}+{A}_{4}{x}^{4}+{A}_{2}{x}^{2}$,”
*Modern Physics Letters A*, vol. 21, no. 21, pp. 1675–1682, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Barakat and O. M. Al-Dossary, “The asymptotic iteration method for the eigenenergies of the asymmetrical quantum anharmonic oscillator potentials $V\left(x\right)={\sum}_{j=2}^{2\alpha}{A}_{j}{x}^{j}$,”
*International Journal of Modern Physics A*, vol. 22, no. 1, pp. 203–212, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Brook,
*Methodus Incrementorum Directa et Inversa (Direct and Reverse Methods of Incrementation)*, impensis Gulielmi Innys, London, UK, 1715. - D. J. Struik,
*A Source Book in Mathematics*, Translated into English, University Press, Cambridge, Mass, USA, 1969. - Wolfram Research,
*Mathematica, Version 8. 0*, Wolfram Research, Champaign, Ill, USA, 2010. - J. P. Boyd, “The devil's invention: asymptotic, superasymptotic and hyperasymptotic series,”
*Acta Applicandae Mathematicae*, vol. 56, no. 1, pp. 1–98, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. S. Hsue and J. L. Chern, “Two-step approach to one-dimensional anharmonic oscillators,”
*Physical Review D*, vol. 29, no. 4, pp. 643–647, 1984. View at Publisher · View at Google Scholar · View at MathSciNet - G. P. Flessas, R. R. Whitehead, and A. Rigas, “On the ${\alpha x}^{2}+{\beta x}^{4}$ interaction,”
*Journal of Physics A*, vol. 16, no. 1, pp. 85–97, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. A. Bonham and L. S. Su, “Use of Hellmann-Feynman and hypervirial theorems to obtain anharmonic vibration-rotation expectation values and their application to gas diffraction,”
*Journal of Chemical Physics*, vol. 45, no. 8, p. 2827, 1996. View at Publisher · View at Google Scholar - M. Reed and B. Simon,
*Methods of Modern Mathematical Physics, IV: Analysis of Operators*, Academic Press, New York, NY, USA, 1978. View at MathSciNet - N. Aquino, “Asymptotic series for the quantum quartic anharmonic oscillator,”
*Journal of Mathematical Chemistry*, vol. 18, no. 2, pp. 349–357, 1995. View at Publisher · View at Google Scholar - B. Bacus, Y. Meurice, and A. Soemadi, “Precise determination of the energy levels of the anharmonic oscillator from the quantization of the angle variable,”
*Journal of Physics A*, vol. 28, no. 14, pp. L381–L385, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Gönül, N. Çelik, and E. Olğar, “A new algebraic approach to perturbation theory,”
*Modern Physics Letters A*, vol. 20, no. 22, pp. 1683–1694, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Koc and E. Olgar, 2010, Preprint math-ph/1008.0697.
- K. Banerjee, “General anharmonic oscillators,”
*Proceedings of the Royal Society of London A*, vol. 364, no. 1717, pp. 265–275, 1978. View at Publisher · View at Google Scholar · View at MathSciNet - G. R. P. Borges, E. D. Filho, and R. M. Ricotta, “Variational supersymmetric approach to evaluate Fokker-Planck probability,”
*Physica A*, vol. 389, no. 18, pp. 3892–3899, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Khare and Y. P. Varshni, “Is shape invariance also necessary for lowest order supersymmetric WKB to be exact?”
*Physics Letters A*, vol. 142, no. 1, pp. 1–4, 1989. View at Google Scholar · View at Scopus - F. So and K. L. Liu, “Study of the Fokker-Planck equation of bistable systems by the method of state-dependent diagonalization,”
*Physica A*, vol. 277, no. 3, pp. 335–348, 2000. View at Publisher · View at Google Scholar · View at Scopus