Non-Relativistic Bound State Solutions of Modified Quadratic Yukawa plus q-Deformed Eckart Potential ()
1. Introduction
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 [1], while others can only be solved approximately [2] [3], with the use of different approximation schemes [4] [5]. 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 [6], functional analysis approach [7] [8] [9], supersymmetry quantum mechanics (SUSYQM) [10], Nikiforov-Uvarov method (NU) [11] [12] 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 [13] [14] [15]. Recently, [16] 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 [17] [18]. The screened Coulomb potential is used mostly in short-range interactions [19] [20]. 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 [21]. 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 [22] [23]. 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:
, (1)
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 [2]. 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
(2)
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
(3)
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
(4)
The energy equation is given as
(5)
The total wave function is given by
(6)
3. The Radial Solutions of Schrödinger Equation of the Proposed Potential
The Schrödinger equation is given by
(7)
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 [7] [14].
Substituting Equation (1) into (7) gives
(8)
Let’s define suitable approximation to the centrifugal term as
(9)
Substituting Equation (9) into (8) gives
(10)
In order to present Equation (10) in a form that is conformable to Equation (3) then, there is a need for the transformation
and
(11)
Substituting (11) into (10) and simplifying gives
(12)
where
and
(13)
Comparing with parametric NU
(14)
(15)
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
(16)
Factorizing and rearranging Equation (16)
(17)
(18)
3.2. Deductions from the Energy Equation
Case 1:
If we replace parameter
and
, then the potential function reduces to Poschl-Teller potential which is written as:
(19)
The energy equation for Poschl-Teller potential becomes:
(20)
Case 2:
If we also replace parameter
we obtain exponential inversely quadratic potential which is expressed as
(21)
Also, the energy equation becomes:
(22)
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,
(21)
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
.
5. Discussion
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.
6. Conclusion
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.
Acknowledgements
The authors are very grateful to the reviewer for the significant suggestions which have led to the improvement of this article.