The Schrödinger wave equation is primarily considered as one of the most commonly used differential equations in non-relativistic quantum mechanics   . 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  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 . Therefore, in most of these potentials, the quasi-exactly solvable potentials are the ones that provide a part of the energy spectrum .
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 , the factorization method  , functional analysis approach (FAA) , supersymmetric (SUSY) quantum mechanics , shape invariance (SI) , the Nikiforov-Uvarov (NU) method , exact quantization rule   and asymptotic iteration method (AIM) .
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  - . 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.
The Nikiforov-Uvarov (NU)  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    , relativistic Dirac, Klien-Gordon equation, etc., can be transformed, using a suitable coordinate transformation s = s(t), into the form:
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  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 
The conditions for the energy eigenvalues and the corresponding eigenstates are, respectively, given as
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 
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
the radial part of Schrödinger equation is given by
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   
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 
Substituting Equation (17) into Equation (16) and using Equation (18) with one obtains Equation (2), where
Comparing Equation (38) with the parameters Equations (5) to (13), one gets
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
3.2. Combination of Manning-Rosen Plus Quadratic Yukawa Potentials
The Manning-Rosen plus quadratic Yukawa potential is given by  
where are potential parameters and is the screening parameter.
Substituting Equation (35) into Equation (16) gives
Substituting Equation (18) in Equation (36) leads
Now, changing of variables using to transform the equation to the form of Equation (2), one obtains
Comparing Equation (38) with the parameters Equations (5) to (13), one gets
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
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
The second case is Maning-Rosen plus quadratic Yukawa potential which can be obtained by setting to zero. The eigenvalues obtained are as follows
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 
The fourth case is Manning-Rosen which can be obtained by substituting and the resulting eigenvalues are given as
which agrees with the results in  and .
The fifth case is Mobius square potential which can be obtained by substituting . The eigenvalues resulting are given by 
The sixth case is inversely quadratic Yukawa potential which can be obtained by setting and the eigenvalues produced are given by 
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 
and the eigenvalues obtained by substituting the parameters in (31) will be
which is the same results as in  and .
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 .
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.
We thank the Editor and the referees for their valuable comments. This generous support is greatly appreciated.
 Hou, C.-F., Zhou, Z.-X. and Li, Y. (1999) Bound States of the Klein-Gordon Equation with Vector and Scalar Wood-Saxon Potentials. Acta Physica Sinica, 8, 561-564.
 Ikhdair, S.M. and Sever, R. (2008) Polynomial Solutions of the Mie-Type Potential in the D-Dimensional Schrödinger Equation. Journal of Molecular Structure: THEOCHEM, 855, 13-17.
 Sever, R., Tezcan, C. and Yesiltas, O. (2008) Exact Solution of Effective Mass Schrödinger Equation for the Hulthen Potential. International Journal of Theoretical Physics, 47, 2243-2248.
 Alhaidari, A.D. (2007) Representation Reduction and Solution Space Contraction in Quasi-Exactly Solvable Systems. Journal of Physics A: Mathematical and Theoretical, 40, 6305-6328.
 Qiang, W.C., Gao, Y. and Zhou, R.S. (2008) Arbitrary/-State Approximate Solutions of the Hulthén Potential through the Exact Quantization Rule. Central European Journal of Physics, 6, 356-362.
 Ikhdair, S.M. and Sever, R. (2009) Exact Quantization Rule to the Kratzer-Type Potentials: An Application to the Diatomic Molecules. Journal of Mathematical Chemistry, 45, Article No. 1137.
 Hamzavi, H., Rajabi, A.A. and Hassanabadi, H. (2010) Exact Pseudospin Symmetry Solution of the Dirac Equation for Spatially-Dependent Mass Coulomb Potential Including a Coulomb-Like Tensor Interaction via Asymptotic Iteration Method. Physics Letters A, 374, 4303-4307.
 Ita, B.I., Louis, H., Orosun, M.M., Nzeata-Ibe, N., Ikeuba, A.I. and Magu, T.O. (2018) Radial Solutions of the Schrodinger’s Equation with Superposition of Yukawa and Manning-Rosen Applying the Nikiforov-Uvarov Method. Communications in Applied Sciences, 6, 1-10.
 Zhang, M.-C. and An, B. (2010) Analytical Solutions of the Manning-Rosen Potential in the Tridiagonal Program. Chinese Physics Letters, 27, Article ID: 110301.
 Okorie, U.S., Ikot, A.N., Onyeaju, M.C. and Chukwuocha, E.O. (2018) Bound State Solutions of Schrödinger Equation with Modified Mobius Square Potential (MMSP) and Its Thermodynamic Properties. Journal of Molecular Modeling, 24, Article No. 289.
 Onate, C.A. and Ojonubah, J.O. (2016) Eigensolutions of the Schrödinger Equation with a Class of Yukawa Potentials via Supersymmetric Approach. Journal of Theoretical and Applied Physics, 10, 21-26.
 Ikot, A.N., Zarrinkamar, S., Ibanga, E.J., Maghsoodi, E. and Hassanabadi, H. (2014) Pseudospin Symmetry of the Dirac Equation for a Möbius Square Plus Mie Type Potential with a Coulomb-Like Tensor Interaction via SUSYQM. Chinese Physics C, 38, Article ID: 013101.
 Maghsoodi, E., Hassanabadi, H., Rahimov, H. and Zarrinkamar, S. (2013) Arbitrary-State Solutions of the Dirac Equation for a Möbius Square Potential Using the Nikiforov-Uvarov Method. Chinese Physics C, 37, Article ID: 043105.
 Yazarloo, B.H., Hassanabadi, H. and Zarrinkamar, S. (2012) Oscillator Strengths Based on the Möbius Square Potential under Schrödinger Equation. The European Physical Journal Plus, 127, Article No. 51.
 Ikot, A.N., Hassanabadi, H., Maghsoodi, E. and Yazarloo, B.H. (2014) Bound and Scattering States of Modified Yukawa Potential under Relativistic Spin and Pseudospin Symmetries with Three Tensor Interactions. The European Physical Journal Plus, 129, Article No. 218.
 Ikot, A.N, Maghsoodi, E., Zarrinkamar, S. and Hassanabadi, H. (2013) Relativistic Spin and Pseudospin Symmetries of Inversely Quadratic Yukawa-Like Plus Mobius Square Potentials Including a Coulomb-Like Tensor Interaction. Few-Body Systems, 54, 2027-2040.
 Ikot, A.N., Obong, H.P. and Hassanabadi, H. (2015) Relativistic Symmetries of (D+1) Dimensional Dirac Equation with Multiparameter Exponential-Type Potentials Using Supersymmetric Quantum Mechanics. Few-Body Systems, 56, 185-196.
 Antia, A.D., Ikot, A.N., Hassanabadi, H. and Maghsoodi, E. (2013) Bound State Solutions of Klein-Gordon Equation with Mobius Square Plus Yukawa Potentials. Indian Journal of Physics, 87, 1133-1139.
 Farout, M., Sever, R. and Ikhdair, S.M. (2020) Approximate Solution to the Time-Dependent Kratzer Plus Screened Coulomb Potential in the Feinberg-Horodecki Equation. Chinese Physics B, 29, Article ID: 060303.
 Farout, M. and Ikhdair, S.M. (2020) Momentum Eigensolutions of Feinberg-Horodecki Equation with Time-Dependent Screened Kratzer-Hellmann Potential. Journal of Applied Mathematics and Physics, 8, 1207-1221.
 Farout, M., Bassalat, A. and Ikhdair, S.M. (2020) Exact Quantized Momentum Eigenvalues and Eigenstates of a General Potential Model. Journal of Applied Mathematics and Physics, 8, 1434-1447.
 Farout, M., Bassalat, A. and Ikhdair, S.M. (2020) Feinberg-Horodecki Exact Momentum States of Improved Deformed Exponential-Type Potential. Journal of Applied Mathematics and Physics, 8, 1496-1506.
 Tezcan, C. and Sever, R. (2009) A General Approach for the Exact Solution of the Schrödinger Equation. International Journal of Theoretical Physics, 48, 337-350.
 Ikot, A.N. (2012) Solutions to the Klein-Gordon Equation with Equal Scalar and Vector Modified Hylleraas Plus Exponential Rosen Morse Potentials. Chinese Physics Letters, 29, Article ID: 060307.
 Ahmaadov, A.I., Demirci, M. and Aslanova, S.M. (2019) Bound State Solutions of the Klein-Fock-Gordon Equation with the Sum of Manning-Rosen Potential and Yukawa Potential within SUSYQM. 26th International Conference on Integrable Systems and Quantum Symmetries, Prague, 8-12 July 2019, Article ID: 012001.
 Edet, C.O., Okoi, P.O. and Chima, S.O. (2020) Analytic Solutions of the Schrödinger Equation with Non-Central Generalized Inverse Quadratic Yukawa Potential. Revista Brasileira de Ensino de Física, 42, Article ID: e20190083.
 Ikot, A.I., Maghsoodi, E., Antia, A.D., Zarrinkamar, S. and Hassanabadi, H. (2013) Approximate κ-State Solutions to the Dirac Mobius Square-Yukawa and Mobius square-Quasi Yukawa Problems under Pseudospin and Spin Symmetry Limits with Coulomb-Like Tensor Interaction. Canadian Journal of Physics, 91, 560-575.
 Antia, A.D., Essien, I.E., Umoren, E.B. and Eze, C.C. (2015) Approximate Solutions of the Non-Relativistic Schrödinger Equation with Inversely Quadratic Yukawa Plus Mobius Square Potential via Parametric Nikiforov-Uvarov Method. Advances in Physics Theories and Applications, 44, 1-13.
 Falaye, B.J., Oyewumi, K.J., Ibrahim, T.T., Punyasena, M.A. and Onate, C.A. (2013) Bound State Solutions of the Manning-Rosen Potential. Canadian Journal of Physics 91, 98-104.
 Ikhdair, S.M. (2011) On the Bound-State Solutions of the Manning-Rosen Potential Including an Improved Approximation to the Orbital Centrifugal Term. Physica Scripta, 83, Article ID: 015010.
 Ita, I. B., Ikeuba, A.I., Louis, H. and Tchoua, P. (2015) Solutions of the Schrödinger Equation with Inversely Quadratic Yukawa plus Attractive Radial Potential Using Nikiforov-Uvarov Method. AIP Conference Proceedings, 1629, 360-362.
 Jia, C.S., Liu, J.Y. and Wang, P.Q. (2008) A New Approximation Scheme for the Centrifugal Term and the Hulthén Potential. Physics Letters A, 372, 4779-4782.
 Ikhdair, S.M. and Abu-Hasna, J. (2011) Quantization Rule Solution to the Hulthén Potential in Arbitrary Dimension with a New Approximate Scheme for the Centrfugal Term. Physica Scripta, 83, Article ID: 025002.