Researchers have put on their interest over the years with the aim of investigating the bound state solutions of relativistic and nonrelativistic wave equations for different potentials. A few of these potentials have been solved exactly , while others can only be solved approximately  , with the use of different approximation schemes  . Subsequently, various methods have been applied to obtain the solutions of the nonrelativistic wave equations with a chosen potential model. Methods include the factorization method , functional analysis approach   , supersymmetry quantum mechanics (SUSYQM) , Nikiforov-Uvarov method (NU)   and others.
The Eckart potential is one of the most important exponential-type potentials in the literature and it has been widely used in physics and chemical physics   . Recently,  has employed the asymptotic iteration method to calculate any ℓ-state solutions of the Schrödinger equation with the Eckart potential by proper approximation of the centrifugal term. Energy eigenvalues and corresponding eigen functions are obtained explicitly. The energy eigenvalues are calculated numerically for some values of ℓ and n.
The Yukawa potential, which is also known as the screened Coulomb potential, is greatly important with applications cutting across nuclear physics and condensed matter physics  . The screened Coulomb potential is used mostly in short-range interactions  . The Yukawa potential is known to be the potential of a charged particle in a weakly non-ideal plasma. It also describes the charged particle effects in a sea of conduction electrons in solid-state physics.
With approximate analytical solution of the Yukawa potential with arbitrary momenta using the Nikiforov-Uvarov method, they obtained approximate analytical solutions of the radial Schrödinger equation for the Yukawa potential, and the energy eigenvalues and the corresponding eigen functions are calculated in closed forms . Some numerical results are presented and show that these results are in good agreement with those obtained previously by other methods. As the screening parameter tends to zero, the energy levels of Yukawa potential are equivalent to that of the familiar pure Coulomb potential. Eckart potential is one of the most important exponential type potentials for the study of nuclear forces  . Basically, many properties of nuclear forces are explained using Yukawa potential and to predict the existence of pi-meson. This is the main reason why the study of the combination of the two potential models is very significant and highly essential. However, after the above mentioned studies on these different potentials and their lofty importance, we seek to investigate the bound state solutions of the Schrodinger equation with the modified quadratic Yukawa plus q-deformed Eckart potential of the form:
where are the potential depth, q is the deformation parameter, r is the inter-nuclear distance and is the screening parameter. It can be deduced that when and , the above combined potential reduces to the Poschl-Teller potential. When , the potential reduces to Exponential inversely quadratic potential. Using the parametric NU method, we derive the energy bound state solutions and their wave functions of the Schrodinger equation for the modified quadratic Yukawa plus q-deformed Eckart potential, analytically and numerically. The rest of the manuscript is organized as follows: In Section 2, a review of Nikiforov-Uvarov method is presented. The radial solutions of Schrodinger equation are presented in Section 3. The numerical solutions are presented in Section 4 while discussion and conclusion are presented in Sections 5 and 6 respectively.
2. Nikiforov-Uvarov Method: Parametric Method
The parametric form is simply using parameters to obtain explicitly energy eigenvalues and it is based on the solutions of a generalized second order linear differential equation with special orthogonal functions . The hypergeometric NU method has shown high utility in calculating the exact energy levels of all bound states for some solvable quantum systems.
Given a second order differential equation of the form
where and are polynomials at most degree and is first degree polynomials. The parametric generalization of the N-U method is given by the generalized hypergeometric-type equation
Thus Equation (2) can be solved by comparing it with Equation (3) and the following polynomials are obtained
The parameters obtainable from Equation (3) serve as an important tool for finding the energy eigenvalue and eigenfunctions.
Where and are parametric constants, other parametric constants are given below
The energy equation is given as
The total wave function is given by
3. The Radial Solutions of Schrödinger Equation of the Proposed Potential
The Schrödinger equation is given by
where is the reduced mass, is the rotational vibrational energy spectra of the diatomic molecules, is the reduced Planck’s constant, n and l are the radial and orbital angular momentum quantum numbers, respectively  .
Substituting Equation (1) into (7) gives
Let’s define suitable approximation to the centrifugal term as
Substituting Equation (9) into (8) gives
In order to present Equation (10) in a form that is conformable to Equation (3) then, there is a need for the transformation
Substituting (11) into (10) and simplifying gives
Comparing with parametric NU
3.1. Energy of the Proposed Potential
The energy can be calculated by substituting Equation (14) and Equation (15) into (5) with simple mathematical algebraic simplification. Hence, the energy eigen equation become
Factorizing and rearranging Equation (16)
3.2. Deductions from the Energy Equation
If we replace parameter and , then the potential function reduces to Poschl-Teller potential which is written as:
The energy equation for Poschl-Teller potential becomes:
If we also replace parameter we obtain exponential inversely quadratic potential which is expressed as
Also, the energy equation becomes:
3.3. The Wave Function of the Proposed Potential
The wave function for the proposed potential can be calculated using equation (6), expressed in terms of Jacobi polynomial and substituting the necessary parameters. Thus,
where is the normalization constant.
4. The Numerical Computation for Energy Equation of the Proposed Potential Using Schrödinger Equation
The numerical computation was carried out using MATLAB 8.0 software while graphical plot was carried out using origin software. The computation parameters have been incorporated in the numerical table. The energy equation for the proposed potential is given in Equation (18). The computation was carried for various range of values of the screening i.e. for large and small values of .
The results of the numerical computation that were carried out for the screening parameter α = (0.1 to 0.5) with the potential depths V0 = 0.3, V1 = 0.1, V2 = 0.2. Tables 1-5 show that for α = (0.1 to 0.5), the numerical bound state energy decreases with an increase in quantum state. It is also observed that the energy increases with a decrease in screening parameter α. However, the energy spectral diagram as shown in Figures 1-5 shows unique quantization of different energy levels with respect to the quantum state. The numerical bound state energies obtained for the various values of α are predominantly negative which the sufficient and necessary condition for bound state solutions. Also the negative bound state eigenvalues obtained show that the potential is suitable for describing particle state and their bound state properties.
We used the Nikiforov-Uvarov method to compute bound state solution of non-relativistic wave equation. We developed a new potential model which is the
Table 1. Numerical computation for .
Table 2. Numerical computation for .
Table 3. Numerical computation for .
Table 4. Numerical computation for .
Figure 1. Energy spectral diagram with Schrödinger equation for .
Table 5. Numerical computation for .
Figure 2. Energy spectral diagram with Schrödinger equation for .
Figure 3. Energy spectral diagram with Schrödinger equation for .
Figure 4. Energy spectral diagram with Schrödinger equation for .
Figure 5. Energy spectral diagram with Schrödinger equation for .
modified quadratic Yukawa plus q-deformed Eckart potential. This potential reduces to a well-known potential called the Poschl-Teller potential and the exponential inversely quadratic potential. We also implement some algorithms to carry out numerical computation for the resulting energy equations for various values of screening parameters. Our analytic and computational results for the proposed potential are in agreement with the existing literature.
The authors are very grateful to the reviewer for the significant suggestions which have led to the improvement of this article.
 Ikot, A., Akpabio, E. and Umoren, B. (2011) Exact Solutions of Schrodinger Equation with Inverted Woods-Saxon and Manning-Rosen Potential. Journal of Scientific Research, 3, 25-33.
 Okon, I., Popoola, O. and Ituen, E. (2016) Bound State Solution to Schrödinger Equation with Hulthen plus Exponential Coulombic Potential with Centrifugal Potential Barrier Using Parametric Nikiforovuvarov Method. International Journal of Recent advances in Physics, 2, 5-7.
 Louis, H., Ita, B., Amos, P. and Phillip, P. (2018) Bound State Solutions of Klein-Gordon Equation with Manning-Rosen Plus a Class of Yukawa Potential Using Pekeris-Like Approximation of the Coulomb Term and Parametric Nikiforov-Uvarov. International Journal of Chemical Physical Science, 1, 7.
 Berkdemir, C. and Han, J. (2005) Any l-State Solutions of Morse Potential through the Pekeris Approximation and Nikiforov-Uvarov Method. Chemical Physics Letters, 409, 203-207.
 Qiang, W. and Dong, S. (2007) Arbitrary l-State Solutions of the Rotating Morse Potential through the Exact Quantization Rule Method. Physics Letters A, 363, 169.
 Tang, H., Liang, C., Zhang, H., Zhao, F. and Jia, C. (2014) Diatomic Molecule Energies of the Modified Rosen-Morse Potential Energy Model. Canadian Journal of Chemistry, 92, 341.
 Miftahul, M., Naufal, I., Suparmi, A. and Cari, C. (2019) Exact Solution of Schrödinger Equation in Cylindrical Coordinates for Double Ring-Shaped Coulomb Oscillator Potential Using SUSYQM Method. Journal of Physics: Conference Series, 1, 1127.
 Antia, A., Christian, C. and Akpabio, L. (2018) Solutions of Schrödinger Equation with the Harmonic Oscillator Potential (HOP) in Cylindrical Basis. Physics and Astronomy International Journal, 2, 3.
 Parmar, H. (2019) Generalized Improved Non-Central Potential and Solution of Schrödinger Equation with Extended Ring-Shaped Potential through Nikiforov-Uvarov Method. The European Physical Journal Plus, 3, 134.
 Edet, O., Okorie, S., Ngiangia, T. and Ikot, A. (2019) Bound State Solutions of the Schrodinger Equation for the Modified Kratzer Potential plus Screened Coulomb Potential. Indian Journal of Physics, 93, 1171-1179.
 Dong, S., Sun, H. and Dong, H. (2013) Arbitrary l-Wave Solution of the Schrodinger Equation for the Screen Coulomb Potential. International Journal of Modern Physics E, 22, Article ID: 1350036.
 Ikhdair, S. and Sever, R. (2007) Polynomial Solutions of the Mie-Type Potential in the Dimensional Schrodinger Equation. Journal of Molecular Structure, 855, 13-17.
 Liverts, Z., Drukarev, G., Krivec, R. and Mandelzweig, B. (2008) Analytic Presentation of a Solution of the Schrodinger Equation. Few-Body Systems, 44, 367.
 Maghsoodi, E., Hassanabadi, H. and Aydodu, O. (2012) Dirac Particles in the Presence of the Yukawa Potential plus a Tensor Interaction in SUSYQM Frame Work. Physica Scripta, 86, Article ID: 015005.
 Ikhdair, M. and Sever, R. (2009) Exact Quantization Rule to the Kratzer-Type Potentials: An Application to the Diatomic Molecules. Journal of Materials Chemistry, 45, 1137.
 Taskin, F. and Kocak, G. (2010) Approximate Solutions of Schrodinger Equation for Eckart Potential with Centrifugal Term. Chinese Physics B, 9, Article ID: 090314.
 Hamzavi, M.,Thylwe, K. and Rajabi, A. (2012) Approximate Analytical Solution of the Yukawa Potential with Arbitrary Angular Momenta. Chinese Physics Letters, 29, Article ID: 080302.