Approximate Bound State Solutions for Certain Molecular Potentials ()
1. Introduction
The Schrödinger wave equation is primarily considered as one of the most commonly used differential equations in non-relativistic quantum mechanics [1] [2] [3]. However, since the early times of quantum mechanics, the exact solutions of the Schrödinger equation with some particular physical potentials are of much interest. Such solutions provide profound conceptual understanding to physical models and certainly lead to a strong judgment supporting the correctness of quantum theory. The exact solutions of central and non-central potentials find their applications in various branches of physics such as molecular, solid-state and chemical physics [4] and so forth. Our choice for the real potential gives the bound state energy eigenvalues and wave functions of the Schrödinger wave equation that might describe essentially the particle dynamics in non-relativistic quantum mechanics. Moreover, these solutions are used in checking and try to improve models under study and then also finding methods in solving complicated physical models.
Since the early times of quantum mechanics the number of exactly solvable physical problems is very limited. Several authors have paid many efforts toward studying the exactly solvable physical problems by which one can determine the whole energy spectrum analytically for wide range values of potential parameters [5]. Therefore, in most of these potentials, the quasi-exactly solvable potentials are the ones that provide a part of the energy spectrum [6].
Recently, various methods are introduced and employed in quantum mechanics in solving the wave equations with a particular given solvable potential. We mention few among the many methods: the group theoretical technique [7], the factorization method [8] [9], functional analysis approach (FAA) [10], supersymmetric (SUSY) quantum mechanics [11], shape invariance (SI) [12], the Nikiforov-Uvarov (NU) method [13], exact quantization rule [14] [15] and asymptotic iteration method (AIM) [16].
The Manning-Rosen, the quadratic Yukawa and the Mobius square potentials have been intensively considered and studied in non-relativistic and relativistic wave equations in recent years [17] - [27]. Therefore, the main motivation of the present work is to give approximate solution to the non-relativistic Schrödinger equation with the superposition of Manning-Rosen plus inversely Mobius square plus Yukawa potential models. Hence we need to treat the centrifugal term with Greene-Aldrich approximation to enable for analytical solution of the Schrödinger equation for any angular momentum quantum number l. This would provide us the bound state energy spectrum for any angular momentum quantum number l and the corresponding wave functions by simply applying the parametric Nikiforov-Uvarov (pNU) method.
The structure of the present work is as follows. In Section 2, we present the brief methodology. In Section 3, we apply this method to derive the bound state energy and wave functions for the Schrödinger equation with the present potential model. Section 4 presents our results and discussion. Finally, in Section 5, we give our conclusion.
2. Methodology
The Nikiforov-Uvarov (NU) [13] method is an efficient tool which is usually used to reduce the second-order differential equation into a general form of a hypergeometric type. In that sense, any second order differential equation, i.e. Schrödinger, Fienberg-Horodecki [28] [29] [30] [31], relativistic Dirac, Klien-Gordon equation, etc., can be transformed, using a suitable coordinate transformation s = s(t), into the form:
(1)
where
and
are polynomials, at most second-degree, and
is a first-degree polynomial. The method is noted to be tiresome and time-consuming. Therefore, Tezcan and Sever [32] derived a parametric form of the NU method to popularize the method more. The parametric NU method is straightforward, simpler and more accurate for the determination of the energy eigenvalues and the corresponding eigenstates. To apply the parametric NU method, the differential equation must be set into the general form given by [33]
(2)
The conditions for the energy eigenvalues and the corresponding eigenstates are, respectively, given as
(3)
(4)
where
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
where
is the normalisation constant and
is the orthogonal Jacobi polynomial.
3. Solution of the Schrödinger Equation with Two Molecular Potential Models
The Schrödinger equation in spherical coordinates is given as [33]
(14)
where
is the reduced Plank constant,
is the reduced mass, E is the energy eigenvalues, and
is the wave function of the particle. If we define the wave function as
(15)
the radial part of Schrödinger equation is given by
(16)
where n and l are the radial and the angular momentum quantum numbers, respectively.
We shall solve the Schrodinger equation for the following two molecular potential models:
3.1. Combination of Manning-Rosen Plus Mobius Square Plus Quadratic Yukawa Potentials
The general potential is given as [34] [35] [36]
(17)
where
are potential parameters and
is the screening parameter.
It is obvious that Equation (16) cannot be solved analytically due to the quadratic Yukawa and the centrifugal terms. However, this can be addressed using the Green-Aldrich approximation [37]
(18)
Substituting Equation (17) into Equation (16) and using Equation (18) with
one obtains Equation (2), where
(19)
(20)
(21)
Comparing Equation (38) with the parameters Equations (5) to (13), one gets
(22)
(23)
(24)
(25)
(26)
(27)
(28)
(29)
(30)
Substituting the values of the parametric constants Equations (22) to (30) into Equations (3) and (4), respectively, one gets the energy eigenvalues and the corresponding unnormalized radial eigenstates as
(31)
and
(32)
where
(33)
(34)
and
.
3.2. Combination of Manning-Rosen Plus Quadratic Yukawa Potentials
The Manning-Rosen plus quadratic Yukawa potential is given by [34] [35]
(35)
where
are potential parameters and
is the screening parameter.
Substituting Equation (35) into Equation (16) gives
(36)
Substituting Equation (18) in Equation (36) leads
(37)
Now, changing of variables using
to transform the equation to the form of Equation (2), one obtains
(38)
where
(39)
(40)
(41)
Comparing Equation (38) with the parameters Equations (5) to (13), one gets
(42)
(43)
(44)
(45)
(46)
(47)
(48)
(49)
(50)
Substituting the values of the parametric constants Equations (42) to (50) into Equations (3) and (4), respectively, one gets the energy eigenvalues and the corresponding unnormalized radial eigenstates as
(51)
and
(52)
where
(53)
(54)
and
.
3.3. Special Cases
To get special cases some parameters should be set to zero. The first case is Manning-Rosen plus Mobius square which can be obtained by setting
to zero and the energy eigenvalues will be
(55)
where
(56)
The second case is Maning-Rosen plus quadratic Yukawa potential which can be obtained by setting
to zero. The eigenvalues obtained are as follows
(57)
where
(58)
which is the same results as in (51) and (53).
The third case is Mobius square plus inversely quadratic Yukawa potential, which results from substituting
in (17). The eigenvalues resulting from substituting these parameters in (55) are given by [38]
(59)
where
(60)
The fourth case is Manning-Rosen which can be obtained by substituting
and the resulting eigenvalues are given as
(61)
where
(62)
which agrees with the results in [39] and [40].
The fifth case is Mobius square potential which can be obtained by substituting
. The eigenvalues resulting are given by [19]
(63)
where
(64)
The sixth case is inversely quadratic Yukawa potential which can be obtained by setting
and the eigenvalues produced are given by [41]
(65)
where
(66)
The seventh case is the Hulthèn potential which can be obtained by substituting
in the Manning-Rosen Potential and the resulting potential will be [42]
(67)
and the eigenvalues obtained by substituting the parameters in (31) will be
(68)
which is the same results as in [43] and [44].
4. Numerical Results and Discussion
In this work, we have studied the solution of the Schrödinger wave equation with two sets of potentials. Here we tend to explain our results by commenting on the plotted Figures. In Figure 1, we showed the variation in the vibrational energy levels against the screening parameter
. It is noted that as
increases, the energy levels of the system decrease monotonically from zero. It is equally seen that rotational energy levels of the system decrease as the screening parameter increases as shown in Figure 2. Figure 3 shows the energy states decrease as the principal quantum number increases for various values of the screening parameter.
Figure 4 indicates the variation of the vibrational energy levels against the strength parameter D. A decrease in the strength D results in an increase in the energy. Moreover, as the potential strength D decreases beyond some value it results in a sharp decrease in the energy.
In Figure 5, we examined the variation in the energy against the potential strength V1. It is seen that the energy of the system decreases monotonically from zero as the potential strength increases for various values of n. A reverse case is shown in Figure 6 when energy is plotted against the potential strength C. It is obvious that when a particle is subjected to this system, the particle exhibits different features of V1 and C for various values of screening parameters; namely,
,
. and
. However, when the strength parameter C gets a large value, the energy drops sharply for
.
Figure 1. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs
for different values of n.
,
,
and
.
Figure 2. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs
.
,
, and
for
(
),
(
,
), and
(
,
,
).
Figure 3. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs n.
,
,
and
for different values of
.
Figure 4. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs D.
,
,
and
for different values of n.
Figure 5. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs
.
,
,
and
for different values of n.
Figure 6. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs C.
,
,
and
for different values of
.
Figure 7 shows a plot of the variation in the energy against the strength parameter C for various values of n. It is seen that the energy decreases as the strength parameter C increases. It is equally seen that the vibrational energy of the system decreases as the screening parameter increases for various values of n as shown in Figure 8. A similar behavior to Figure 8 appears for the rotational energy levels when plotted against the screening parameters are shown in Figure 9.
Finally, Figure 10 shows the decrease of vibrational energy levels as the quantum number n increases for various values of screening parameter
.
Figure 7. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen plus quadratic Yukawa potential vs C.
,
and
for different values of n.
Figure 8. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen Mobius square plus quadratic Yukawa potential vs
.
,
,
and
for different values of n.
Figure 9. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen Mobius square plus quadratic Yukawa potential vs
.
,
and
for different values of n.
Figure 10. The quantized energy eigenvalues of Schrödinger equation with Manning-Rosen Mobius square plus quadratic Yukawa potential vs n.
,
,
and
for different values of
.
5. Conclusion
In this work, we have analytically solved the non-relativistic wave equation with the combination of three important potentials via the parametric Nikiforov-Uvarov method. We have obtained the energy equation and the corresponding non-normalized wave functions of the combination set of Manning-Rosen plus Mobius square plus quadratic Yukawa potential and their subset of potentials. We have obtained in detail the energy eigenvalues and the corresponding wave function for subset of potentials. These results could find their applications in atomic as well as molecular physics. The effects of the strength parameters as well as screening parameter on the vibrational and rotational energy levels were also studied.
Acknowledgements
We thank the Editor and the referees for their valuable comments. This generous support is greatly appreciated.