Advances in High Energy Physics
Volume 2015 (2015), Article ID 137038, 8 pages
http://dx.doi.org/10.1155/2015/137038

Research Article

Exact Analytical Solution of the -Dimensional Radial Schrödinger Equation with Pseudoharmonic Potential via Laplace Transform Approach

Tapas Das¹ and Altuğ Arda²

¹Kodalia Prasanna Banga High School (H.S), South 24 Parganas, Sonarpur 700146, India
²Department of Physics Education, Hacettepe University, 06800 Ankara, Turkey

Received 16 September 2015; Revised 1 December 2015; Accepted 7 December 2015

Academic Editor: Ming Liu

Copyright © 2015 Tapas Das and Altuğ Arda. 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. The publication of this article was funded by SCOAP³.

Abstract.

The second-order -dimensional Schrödinger equation with pseudoharmonic potential is reduced to a first-order differential equation by using the Laplace transform approach and exact bound state solutions are obtained using convolution theorem. Some special cases are verified and variations of energy eigenvalues as a function of dimension are furnished. To give an extra depth of this paper, the present approach is also briefly investigated for generalized Morse potential as an example.

The author Tapas Das wishes to dedicate this work to his wife Sonali for her love and care

1. Introduction

Schrödinger equation has long been recognized as an essential tool for the study of atoms, nuclei, and molecules and their spectral behaviors. Much effort has been spent to find the exact bound state solution of this nonrelativistic equation for various potentials describing the nature of bonding or the nature of vibration of quantum systems. A large number of research workers all around the world continue to study the ever fascinating Schrödinger equation, which has wide application over vast areas of theoretical physics. The Schrödinger equation is traditionally solved by operator algebraic method [1], power series method [2, 3], or path integral method [4].

There are various other alternative methods in the literature to solve Schrödinger equation such as Fourier transform method [5–7], Nikiforov-Uvarov method [8], asymptotic iteration method [9], and SUSYQM [10]. The Laplace transformation method is also an alternative method in the list and it has a long history. The LTA was first used by Schrödinger to derive the radial eigenfunctions of the hydrogen atom [11]. Later Englefield used LTA to solve the Coulomb, oscillator, exponential, and Yamaguchi potentials [12]. Using the same methodology, the Schrödinger equation has also been solved for various other potentials, such as pseudoharmonic [13], Dirac delta [14], and Morse-type [15, 16] and harmonic oscillator [17] specially on lower dimensions.

Recently, -dimensional Schrödinger equations have received focal attention in the literature. The hydrogen atom in five dimensions and isotropic oscillator in eight dimensions have been discussed by Davtyan and coworkers [18]. Chatterjee has reviewed several methods commonly adopted for the study of -dimensional Schrödinger equations in the large limit [19], where a relevant expansion can be used. Later Yáñez et al. have investigated the position and momentum information entropies of -dimensional system [20]. The quantization of angular momentum in -dimensions has been described by Al-Jaber [21]. Other recent studies of Schrödinger equation in higher dimension include isotropic harmonic oscillator plus inverse quadratic potential [22], -dimensional radial Schrödinger equation with the Coulomb potential [23]. Some recent works on -dimensional Schrödinger equation can be found in the references list [24–31].

These higher dimension studies facilitate a general treatment of the problem in such a manner that one can obtain the required results in lower dimensions just dialing appropriate . The pseudoharmonic potential is expressed in the form [32]

where is the dissociation energy with the force constant and is the equilibrium constant. The pseudoharmonic potential is generally used to describe the rotovibrational states of diatomic molecules and nuclear rotations and vibrations. Moreover, the pseudoharmonic potential and some kinds of it for -dimensional Schrödinger equation help to test the powerfulness of different analytical methods for solving differential equations. To give an example, the dynamical algebra of the Schrödinger equation in -dimension has been studied by using pseudoharmonic potential [33]. Taşeli and Zafer have tested the accuracy of expanding of the eigenfunction in a Fourier-Bessel basis [34] and a Laguerre basis [35] with the help of different type of polynomial potentials in -dimension.

Motivated by these types of works, in this present paper we discuss the exact solutions of the -dimensional radial Schrödinger equation with pseudoharmonic potential using the Laplace transform approach. To make this paper self-contained we briefly outline Laplace transform method and convolution theorem in the next section. Section 3 is for the bound state spectrum of the potential system. Section 4 is devoted to the results and discussion where we derive some well known results for special cases of the potential. The application of the present method is shown in Section 5 where we briefly show how generalized Morse potential could be solved. Finally the conclusion of the present work is placed in Section 6.

2. Laplace Transform Method and Convolution Theorem

The Laplace transform or of a function is defined by [36, 37]

If there is some constant such that for sufficiently large , the integral in (2) will exist for . The Laplace transform may fail to exist because of a sufficiently strong singularity in the function as . In particular

The Laplace transform has the derivative properties

where the superscript denotes the th derivative with respect to for and with respect to for .

The inverse transform is defined as . One of the most important properties of the Laplace transform is that given by the convolution theorem [38]. This theorem is a powerful tool to find the inverse Laplace transform. According to this theorem if we have two transformed functions and , then the product of these two is the Laplace transform of the convolution , where

So the convolution theorem yields

Hence

If we substitute , then we find the important consequence .

3. Bound State Spectrum

The -dimensional time-independent Schrödinger equation for a particle of mass with orbital angular momentum quantum number is given by [39]

where and denote the energy eigenvalues and potential. within the argument of the th state eigenfunctions denotes angular variables . The Laplacian operator in hyperspherical coordinates is written as

where

is known as the hyperangular momentum operator.

We chose the bound state eigenfunctions that are vanishing for and . Applying the separation variable method by means of the solution , (8) provides two separated equations:

where is known as the hyperspherical harmonics, and the hyperradial or in short the “radial” equation

where is the separation constant [40, 41] with .

In spite of taking (1) we take the more general form of pseudoharmonic potential [42]

where , and are three parameters that can take any real value. If we set , , and , (13) converts into the special case which we have given in (1). Taking this into (12) and using the abbreviations

we obtain

We are looking for the bound state solutions for with the following properties: and . Let us assume a solution of type with . Here the term ensures the fact that and is expected to behave like .

Now changing the variable as and taking from (15) we obtain

In order to get an exact solution of the above differential equation we remove the singular term by imposing the condition

So we have

where is taken as the positive solution of , which can be easily found by little algebra from (17). It is worth mentioning here that the condition given by (17) is necessary to get an analytical solution because otherwise only approximate or numerical solution is possible. Introducing the Laplace transform with the boundary condition and using (4), (18) can read

where . This parameter can have integer or noninteger values and there will be integer or noninteger term(s) in energy eigenvalue according to values of which can be seen below.

The solution of the last equation can be written easily as

where is a constant. The term is a positive real number as we restrict ourselves to the choice . Now, since is positive and , then the second factor of (20) could become negative if and thus its power must be a positive integer to get singled valued eigenfunctions. This will also exclude the possibility of getting singularity in the transformation. So we have

Using (21), we have from (20)

where and . In order to find , we find [43]

Therefore using (7) and (23), we have

The integration can be found in [44], which gives

where is the confluent hypergeometric functions. Now using the beta function , can be written:

So we have the radial eigenfunctions

where is the normalization constant. It should be noted that because of the “quantization condition” given by (21), it is possible to write the radial eigenfunctions in a polynomial form of degree , as converges for all finite with and is not a negative integer or zero.

Now the normalization constant can be evaluated from the condition [45]

To evaluate the integration the formula is useful here, where are the Laguerre polynomials. It should be remembered that the formula is applicable only if is a positive integer. Hence identifying ,  , and we have

So using the following formula for the Laguerre polynomials,

we write the normalization constant

Finally, the energy eigenvalues are obtained from (21) along with (17) and (14):

and we write the corresponding normalized eigenfunctions as

or

Finally, the complete orthonormalized eigenfunctions of the -dimensional Schrödinger equation with pseudoharmonic potential can be given by

We give some numerical results about the variation of the energy on the dimensionality obtained from (32) in Figure 1. We summarize the plots for as a function of for a set of physical parameters by taking , especially.

(a) The variation of energy to for

(b) The variation of energy to for

(c) The variation of energy to for

(d) The variation of energy to for

Figure 1: The dependence of on . Parameters: ,  , and  .

4. Results and Discussion

In this section we have shown that the results obtained in Section 3 are very useful in deriving the special cases of several potentials for lower as well as for higher dimensional wave equation.

4.1. Isotropic Harmonic Oscillator

(1) Three Dimensions ( ). For this case and which gives from (32)

where is the circular frequency of the particle. From (14) and (17) we obtain . This makes and we get radial eigenfunctions from (34). The result agrees with those obtained in [22].

(2) Arbitrary -Dimensions. Here is an arbitrary constant and as before ,  . We have the energy eigenvalues from (32):

Solving (17) with the help of (14) we have and one can easily obtain the normalization constant from (31):

where . The radial eigenfunctions are given in (34) with the above normalization constant. The results obtained here agree with those found in some earlier works [22, 46, 47].

4.2. Isotropic Harmonic Oscillator Plus Inverse Quadratic Potential

(1) Two Dimensions ( ). Here we have ,  , and , where is the circular frequency of the particle. So from (31) we obtain

where . can be obtained from (14) and (17) as . Hence (32) gives the energy eigenvalues of the system

This result has already been obtained in [22, 48].

(2) Three Dimensions ( ). Here , , and which gives the energy eigenvalues as

Solving (17) we get and hence . The radial eigenfunction can hence be obtained from (34) which also corresponds to the result obtained in [22].

4.3. 3-Dimensional Schrödinger Equation with Pseudoharmonic Potential

Here (17) gives and this makes . So (31) provides

and hence the normalized eigenfunctions become

with the energy eigenvalues

This result corresponds exactly to the ones given in [13].

5. Short Review of Generalized Morse Potential: An Example

The generalized Morse potential [49] in terms of four parameters , and is

where describes the characteristic range of the potential, is the equilibrium molecular separation, and , and are related to the potential depth. In this short section we will only show Laplace transformable differential equation like (18) can also be achieved if the exponential and singular terms and are properly handled. Looking back to the solution given by (34) we predetermine the solution for (8) for the potential given by (45) as

This substitution facilitates following easier differential equation (without the term) similar to (12) of Section 3:

where .

Let us introduce a new variable . Hence inserting (45) into (47) we have

where .

It is possible to expand the term in power series as within the molecular point of view . So

This expansion can also be rewritten as

By comparing these two it is not hard to justify that ;  ;  .

Inserting (50) into (48) and changing the variable as one can easily get

Now further assuming the solution of the above differential equation , as we did in the previous section with proper requirement for bound state scenario, finally we can construct the following differential equation just like (16) for the perfect platform for Laplace transformation:

where

Imposing the condition as previously and approaching the same way as we did in Section 3 one can obtain the energy eigenvalues and bound state wave functions in terms of confluent hypergeometric function. We have investigated that the results are exact match with [15] for , , and , where describes the depth of the potential.

6. Conclusions

We have investigated some aspects of -dimensional hyperradial Schrödinger equation for pseudoharmonic potential by Laplace transformation approach. It is found that the energy eigenfunctions and the energy eigenvalues depend on the dimensionality of the problem. In this connection we have furnished few plots of the energy spectrum as a function of for a given set of physical parameters , and keeping . The general results obtained in this paper have been verified with earlier reported results, which were obtained for certain special values of potential parameters and dimensionality.

The Laplace transform is a powerful, efficient, and accurate alternative method of deriving energy eigenvalues and eigenfunctions of some spherically symmetric potentials that are analytically solvable. It may be hard to predict which kind of potentials is solvable analytically by Laplace transform, but in general prediction of the eigenfunctions with the form like always opens the gate of the possibility of closed form solutions for a particular potential model. This kind of substitution is called Universal Laplace transformation scheme. In this connection one might go through [50] to check out how Laplace transformation technique behaves over different potentials, specially for Schrödinger equation in lower dimensional domain. It is also true that the technique of Laplace transformation is useful if, inserting the potential into the Schrödinger equation and using Universal Laplace transformation scheme via some suitable parametric restrictions, one is able to get a differential equation with variable coefficient   . This is not easy to achieve every time. However, for a given potential if there is no such achievement, iterative approach facilitates a better way to overcome the situation [51]. The results are sufficiently accurate for such special potentials at least for practical purpose.

Before concluding we want to mention here that we have not succeeded in developing the scattering state solution for the pseudoharmonic potential. If we could develop those solutions using the LTA that would have been a remarkable achievement. The main barrier of this success lies on the realization of complex index in (20), because at scattering state situation, that is, , should be replaced with , where . Maybe this cumbersome situation would attract the researchers to study the subject further and we are looking forward to it.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The author Altuğ Arda thanks Dr. Andreas Fring from City University London and the Department of Mathematics for hospitality where the final version of this work has been completed.

References

J. J. Sakurai, Modern Quantum Mechanics, Addison-Wesley, New York, NY, USA, 1967.

A. Sommerfeld, Wave-Mechanics, Methuen, London, UK, 1930.

D. J. Griffiths, Introduction to Quantum Mechanics, Pearson Prentice Hall, Upper Saddle River, NJ, USA, 2nd edition, 2005.

R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, New York, NY, USA, 1965.

S. A. Ponomarenko, “Quantum harmonic oscillator revisited: a fourier transform approach,” American Journal of Physics, vol. 72, article 1259, 2004.

A. Engel, “Comment on quantum harmonic oscillator revisited: a fourier transform approach,” American Journal of Physics, vol. 74, no. 9, article 837, 2006.

G. Palma and U. Raff, “A novel application of a Fourier integral representation of bound states in quantum mechanics,” American Journal of Physics, vol. 79, no. 2, pp. 201–205, 2011.

A. F. Nikiforov and V. B. Uvarov, Special Functions of Mathematical Physics, Birkhaüser, Bazel, Switzerland, 1988.

O. Bayrak, I. Boztosun, and H. Ciftci, “Exact analytical solutions to the kratzer potential by the asymptotic iteration method,” International Journal of Quantum Chemistry, vol. 107, no. 3, pp. 540–544, 2007.

H. Hassanabadi, L. L. Lu, S. Zarrinkamar, G. H. Liu, and H. Rahimov, “Approximate solutions of schrödinger equation under manning-rosen potential in arbitrary dimension via SUSYQM,” Acta Physica Polonica A, vol. 122, no. 4, pp. 650–654, 2012.

E. Schrödinger, “Quantisierung als eigenwertproblem,” Annalen der Physik, vol. 384, no. 4, pp. 361–376, 1926.

M. J. Englefield, “Solution of Coulomb problem by Laplace transform,” Journal of Mathematical Analysis and Applications, vol. 48, pp. 270–275, 1974.

A. Arda and R. Sever, “Exact solutions of the Schrödinger equation via Laplace transform approach: pseudoharmonic potential and Mie-type potentials,” Journal of Mathematical Chemistry, vol. 50, no. 4, pp. 971–980, 2012.

A. S. de Castro, “Bound states in a double delta potential via laplace transform,” Revista Brasileira de Ensino de Fisica, vol. 34, article 4301, 2012.

S. Miraboutalebi and L. Rajaei, “Solutions of N-dimensional Schrödinger equation with Morse potential via Laplace transforms,” Journal of Mathematical Chemistry, vol. 52, no. 4, pp. 1119–1128, 2014.

A. Arda and R. Sever, “Exact solutions of the morse-like potential, step-up and step-down operators via laplace transform approach,” Communications in Theoretical Physics, vol. 58, no. 1, pp. 27–30, 2012.

D. R. M. Pimentel and A. S. de Castro, “A Laplace transform approach to the quantum harmonic oscillator,” European Journal of Physics, vol. 34, no. 1, pp. 199–204, 2013.

L. S. Davtyan, L. G. Mardoyan, G. S. Pogosyan, A. N. Sissakian, and V. M. Ter-Antonyan, “Generalised KS transformation: from five-dimensional hydrogen atom to eight-dimensional isotrope oscillator,” Journal of Physics A: Mathematical and General, vol. 20, no. 17, pp. 6121–6125, 1987.

A. Chatterjee, “Large-N expansions in quantum mechanics, atomic physics and some O(N) invariant systems,” Physics Report, vol. 186, no. 6, pp. 249–370, 1990.

R. J. Yáñez, W. Van Assche, and J. S. Dehesa, “Position and momentum information entropies of the D-dimensional harmonic oscillator and hydrogen atom,” Physical Review A, vol. 50, no. 4, article 3065, 1994.

S. M. Al-Jaber, “Quantization of angular momentum in the N-dimensional space,” Il Nuovo Cimento B, vol. 110, no. 8, pp. 993–995, 1995.

K. J. Oyewumi and E. A. Bangudu, “Isotropic harmonic oscillator plus inverse quadratic potential in N-dimensional spaces,” The Arabian Journal for Science and Engineering, vol. 28, no. 2, pp. 173–182, 2003.

C. Gang, “Exact solutions of N-dimensional harmonic oscillator via Laplace transformation,” Chinese Physics, vol. 14, no. 6, article 1075, 2005.

H. Hassanabadi, S. Zarrinkamar, and A. A. Rajabi, “Exact solutions of D-dimensional Schrödinger equation for an energy-dependent potential by NU method,” Communications in Theoretical Physics, vol. 55, no. 4, pp. 541–544, 2011.

D. Agboola, “Solutions to the modified pöschl–teller potential in D-dimensions,” Chinese Physics Letters, vol. 27, no. 4, Article ID 040301, 2010.

G. Chen, “Exact solutions of the N-dimensional radial schrödinger equation with the coulomb potential via the laplace tranform,” Zeitschrift für Naturforschung A, vol. 59, no. 11, pp. 875–876, 2004.

S.-H. Dong and G.-H. Sun, “The series solutions of the non-relativistic equation with the Morse potential,” Physics Letters. A, vol. 314, no. 4, pp. 261–266, 2003.

S.-H. Dong, “On the solutions of the Schrödinger equation with some anharmonic potentials: wave function ansatz,” Physica Scripta, vol. 65, no. 4, pp. 289–295, 2002.

S. M. Ikhdair and R. Sever, “An alternative simple solution of the sextic anharmonic oscillator and perturbed Coulomb problems,” International Journal of Modern Physics C, vol. 18, no. 10, pp. 1571–1581, 2007.

S. M. Ikhdair and R. Server, “Exact solutions of the modified kratzer potential plus ring-shaped potential in the D-dimensional schrödinger equation by the nikiforov–uvarov method,” International Journal of Modern Physics C, vol. 19, no. 2, pp. 221–235, 2008.

L.-Y. Wang, X.-Y. Gu, Z.-Q. Ma, and S.-H. Dong, “Exact solutions to D-dimensional Schrödinger equation with a pseudoharmonic oscillator Found,” Foundations of Physics Letters, vol. 15, no. 6, pp. 569–576, 2002.

I. Goldman and V. Krivchenkov, Problems in Quantum Mechanics, Pergamon Press, London, UK, 1960.

D. Martínez, J. C. Flores-Urbina, R. D. Mota, and V. D. Granados, “The su(1, 1) dynamical algebra from the Schrödinger ladder operators for N-dimensional systems: hydrogen atom, Mie-type potential, harmonic oscillator and pseudo-harmonic oscillator,” Journal of Physics A: Mathematical and Theoretical, vol. 43, no. 13, Article ID 135201, 2010.

H. Taşeli and A. Zafer, “A Fourier-Bessel expansion for solving radial Schrödinger equation in two dimensions,” International Journal of Quantum Chemistry, vol. 61, no. 5, pp. 759–768, 1997.

H. Taşeli, “An alternative series solution to the isotropic quartic oscillator in N dimensions,” Journal of Mathematical Chemistry, vol. 20, no. 3-4, pp. 235–245, 1996.

G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists, Academic Press, San Diego, Calif, USA, 5th edition, 2001.

P. K. Chattopadhyay, Mathematical Physics, New Age International Publishers, New Delhi, India, 1st edition, 1990.

M. R. Spiegel, Schaum's Outline of Theory and Problems of Laplace Transforms, McGraw-Hill, New York, NY, USA, 1965.

U. A. Deta, “Approximate solution of Schrödinger equation in dimensions for scarf hyperbolic potential using Nikiforov-Uvarov method,” Advanced Studies in Theoretical Physics, vol. 7, p. 647, 2013.

J. Avery, Hyperspherical Harmonics: Applications in Quantum Theory, Kluwer Academic, Dordrecht, The Netherlands, 1989.

N. Shimakura, Partial Differential Operator of Elliptic Type, American Mathematical Society, Providence, RI, USA, 1992.

S. H. Patil and K. D. Sen, “Uncertainty relations for modified isotropic harmonic oscillator and Coulomb potentials,” Physics Letters A, vol. 362, no. 2-3, pp. 109–114, 2007.

J. L. Schiff, The Laplace Transform: Theory and Applications, Undergraduate Texts in Mathematics, Springer, New York, NY, USA, 1999.

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Academic Press, New York, NY, USA, 7th edition, 2007.

T. Das, “A Laplace transform approach to find the exact solution of the N-dimensional Schrödinger equation with Mie-type potentials and construction of Ladder operators,” Journal of Mathematical Chemistry, vol. 53, no. 2, pp. 618–628, 2015.

V. L. Bakhrakh, S. I. Vetchinkin, and S. V. Khristenko, “Green's function of a multidimensional isotropic harmonic oscillator,” Theoretical and Mathematical Physics, vol. 12, no. 2, pp. 776–778, 1972.

C. Grosche and F. Steiner, Handbook of Feynman Path Integrals, Springer, New York, NY, USA, 1998.

B. K. Cheng and M. G. da Luz, “Exact propagator for a two-dimensional inverse quadratic oscillator interacting with a wedge,” Journal of Physics A: Mathematical and General, vol. 25, no. 7, pp. 2033–2042, 1992.

M. Chabab, R. Jourdani, and M. Oulne, “Exact solutions for vibrational states with generalized q-deformed Morse potential within the asymptotic iteration method,” International Journal of Physical Sciences, vol. 7, no. 8, pp. 1150–1154, 2012.

G.-Y. Tsaur and J. Wang, “A universal Laplace-transform approach to solving Schrödinger equations for all known solvable models,” European Journal of Physics, vol. 35, no. 1, Article ID 015006, 2014.

T. Barakat, “The asymptotic iteration method for the eigenenergies of the anharmonic oscillator potential ,” Physics Letters A, vol. 344, no. 6, pp. 411–417, 2005.