Abstract

The solutions to the Schrödinger equation with inversely quadratic Yukawa and inversely quadratic Hellmann (IQYIQH) potential for any angular momentum quantum number have been presented using the Nikiforov-Uvarov method. The bound state energy eigenvalues and the corresponding unnormalized eigenfunctions are obtained in terms of the Laguerre polynomials. The NU method is related to the solutions in terms of generalized Jacobi polynomials. In the NU method, the Schrödinger equation is reduced to a generalized equation of hypergeometric type using the coordinate transformation . The equation then yields a form whose polynomial solutions are given by the well-known Rodrigues relation. With the introduction of the IQYIQH potential into the Schrödinger equation, the resultant equation is further transformed in such a way that certain polynomials with four different possible forms are obtained. Out of these forms, only one form is suitable for use in obtaining the energy eigenvalues and the corresponding eigenfunctions of the Schrödinger equation.

1. Introduction

The bound state solutions to the Schrödinger equation (SE) are only possible for some potentials of physical interest [15]. Quite recently, several authors have tried to solve the problem of obtaining exact or approximate solutions to the Schrödinger equation for a number of special potentials [610]. Some of these potentials are known to play very important roles in many fields of physics such as molecular physics, solid state, and chemical physics [8].

The purpose of the present work is to present the solution to the Schrödinger equation with the inversely quadratic Yukawa potential [11] plus inversely quadratic Hellmann potential [12] of the forms The sum of these potentials can be written as where represents the internuclear distance, and are the strengths of the Coulomb and Yukawa potentials, respectively, is the screening parameter, and is the dissociation energy. Equation (2) is then amenable to Nikiforov-Uvarov method. Ita [13] has solved the Schrödinger equation for the Hellman potential and obtained the energy eigenvalues and their corresponding wave functions using expansion method and Nikiforov-Uvarov method. Also, Hamzavi and Rajabi [14] have used the parametric Nikiforov-Uvarov method to obtain tensor coupling and relativistic spin and pseudospin symmetries of the Dirac equation with the Hellmann potential. Kocak et al. [15] solved the Schrödinger equation with the Hellmann potential using asymptotic iteration method and obtained energy eigenvalues and the wave functions. However, not much has been achieved in the area of solving the radial Schrödinger equation for any angular momentum quantum number, , with IQYIQH potential using Nikiforov-Uvarov method in the literature.

2. Overview of the Nikiforov-Uvarov Method

The Nikiforov-Uvarov (NU) method is based on the solutions to a generalized second-order linear differential equation with special orthogonal functions [16]. The Schrödinger equation of the type: could be solved by this method. This can be done by transforming (3) into an equation of hypergeometric type with appropriate coordinate transformation to get To find the exact solution to (4), we write as Substitution of (5) into (4) yields (6) of hypergeometric type In (5), the wave function is defined as the logarithmic derivative [17] with being at most first-order polynomials. Also, the hypergeometric-type functions in (6) for a fixed integer are given by the Rodrigue relation as where is the normalization constant and the weight function must satisfy the condition with In order to accomplish the condition imposed on the weight function it is necessary that the polynomial be equal to zero at some points of an interval and its derivative at this interval at will be negative [17]. That is, The function and the parameter required for the NU method are then defined as [17] The -values in (12) are possible to evaluate if the expression under the square root is square of polynomials. This is possible if and only if its discriminant is zero. Therefore, the new eigenvalue equation becomes [18] A comparison between (13) and (14) yields the energy eigenvalues.

3. The Schrödinger Equation

In spherical coordinate, the Schrödinger equation with the potential is given as [19] Using the common ansatz for the wave function in (9) we get the following set of equations: where and are the separation constants. is the solution to (18) and (19) and their solutions are well known as spherical harmonic functions [19].

4. Solutions to the Radial Equation

Equation (17) is the radial part of the Schrödinger equation which we are interested in solving. Equation (17) together with the potential in (2) and with the transformation yields the following equation: where the radial wave function is and Equation (14) is then compared with (4) and the following expressions are obtained: We then obtain the function by substituting (22) into (12): According to the NU method, the quadratic form under the square-root sign of (23) must be solved by setting the discriminant of this quadratic equation equal to zero; that is, . This discriminant gives a new equation which can be solved for the constant to get the two roots as Thus, we have When the two values of given in (25) are substituted into (23), the four possible forms of are obtained as One of the four values of the polynomial is just proper to obtain the bound state solution since given in (4) must have negative derivative. Therefore, the most suitable expression of is chosen as for . We obtain from (10) and the derivative of this expression would be negative; that is, . From (19) and (20) we obtain When we compare these expressions, , we obtain the energy of the IQYIQH potential as Let us now calculate the radial wave function, . Using and (7) and (9), the following expressions are obtained: Then from (8) one has is a normalization constant. The wave function can be obtained in terms of the generalized Laguerre polynomials as is the normalization constant.

5. Discussion

In summary, we have obtained the energy eigenvalues and the corresponding unnormalized wave function using the NU method for the Schrödinger equation with the inversely quadratic Yukawa plus inversely quadratic Hellmann potential. If we set the parameters and , it is easy to show that (29) reduces to the bound state energy spectrum of a particle in the Coulomb potential; that is, , where is the principal quantum number. Similarly, if we set , (29) results in the bound state energy spectrum of a vibrating-rotating diatomic molecule subject to the inversely quadratic Yukawa potential as follows: This result is similar to that obtained in [11]. Also, if we set , , and , equation reduces to the bound state energy spectrum of a vibrating-rotating diatomic molecule subject to inversely quadratic Hellmann potential as follows: Equation (34) is also similar to the one obtained in [12]. These show the accuracy of our calculations.

6. Conclusion

The bound state solutions to the Schrödinger equation have been obtained for the inversely quadratic Yukawa and inversely quadratic Hellmann potential. Special cases of the potential are also considered and their energy eigenvalues are obtained.

Conflict of Interests

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