Bound state solutions and thermodynamic properties of modified exponential screened plus Yukawa potential

Introduction In quantum mechanics, the study of relativistic and nonrelativistic wave equations arouse the interest of different researchers [1–5]. Schrӧdinger equation is the nonrelativistic wave equation, while Dirac and Klein–Gordon equations are relativistic wave equation describing spin-half and spinless particles, respectively [6–10]. The total wave function provides implicitly the information about the quantum mechanical system. [11]. Researchers have adopted many methods in providing solutions to both relativistic and nonrelativistic wave equations. Among them are Nikiforov–Uvarov method [12–16], supersymmetric quantum mechanics Abstract

approach [17], exact proper quantization [18] asymptotic iteration method [19], Wentzel-Kramers-Brillouin (WKB) approach [20] and many others [21][22][23][24][25][26][27][28][29]. These techniques have been used to solve some quantum mechanical potentials, like Hulthen, Yukawa, Poschl-Teller, Tietz-Wei, Tietz-Hua, exponential-type potentials, hyperbolic potentials, Kratzer potentials, screened-Kratzer potential, Mobius square, Hellmann, coulomb, Cornel, Killingbeck, Woods-Saxon, Deng-Fan, Hylleraas, Eckart, pseudoharmonic, Poschl-Teller, modified Yukawa potential and many others [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48]. Yukawa potential which is otherwise known as screened Coulomb potential is a short-range potential with application in particle, high energy and molecular physics which is basically used for the description of interaction existing between atoms of diatomic molecules [49]. A lot of research work has been carried out on thermodynamic properties of some considerable potentials. Okon et.al [50] studied spin and pseudospin solutions of Dirac equation and its thermodynamic properties using hyperbolic Hulthen plus hyperbolic exponential inversely quadratic potential where they obtained numerical bound state solutions for both spin and pseudospin symmetries. Using nonrelativistic limit, they also obtain nonrelativistic energy eigenequation presented in a close form to study partition function and other thermodynamic properties. Okon et al. [51] studied thermodynamic properties and bound state solutions of Schrӧdinger equation using Mobius square plus screened-Kratzer potential within the framework of parametric Nikiforov-Uvarov (NU) method. In their study, they applied NU and semi-classical WKB to obtain bound state solutions and thermodynamic properties for two diatomic molecules (Carbon (II) oxide and Scandium Fluoride). Okorie et al. [52] studied energy spectra and thermodynamic properties of hyperbolic Poschl-Teller potential model. In that work, they solved Dirac equation using modified factorization method to obtain both relativistic and nonrelativistic ro-vibrational energy spectra and thermodynamic properties as applied to some diatomic molecules: Hydrogen chloride (HCl), chlorine (Cl 2 ), carbon (II) oxide and lithium hydride (LiH). The potential they used is reduced to screened-Kratzer and Kratzer potential as special cases. Onate et al. [53] studied bound state solutions and thermal properties of modified Tietz-Hua potential using supersymmetric quantum mechanics approach where they study ro-vibrational energy spectra and thermodynamic properties as applicable to some diatomic properties. Purohit et al. [54] studied eigensolutions and various properties of the screened cosine Kratzer potential in D-dimension via relativistic and nonrelativistic treatment. Here, the screened-Kratzer potential model was extended to study rotational and vibrational energies for few heterogeneous diatomic molecules. They authors also extend the potential model to study the partition function as well as information theoretic measures like Tsallis, Renyi, Shannon and Fisher information entropies. Recently, researchers have studied potential model under the influence of Aharonov-Bohm flux and external magnetic field. Purohit et al. [55] studied the thermomagnetic properties of the screened-Kratzer potential under the influence of Aharonov-Bohm flux and external magnetic field where the obtained thermomagnetic properties as well as persistent current, magnetization and magnetic susceptibility. Okon et al. [56] investigated the effect of Aharonov-Bohm and external magnetic field on Hellmann plus screened-Kratzer potential where the energy equation was presented in a closed form and applied to study thermomagnetic properties as applicable to some diatomic molecules. Here, the authors obtained a normalized wave function expressed in terms of Jacobi polynomials as well as the wave function and probability density plots for some selected diatomic molecules. They also obtained other thermomagnetic properties including partition function, vibrational mean energy, vibrational heat capacity, magnetization, persistent current and magnetic susceptibility. Based on this motivation, in this paper, we study an approximate bound state solutions of Schrӧdinger equation and thermodynamic properties of newly proposed potential called modified exponential screened plus Yukawa potential (MESPYP) within the framework of parametric Nikiforov-Uvarov method. This is a new potential model that has not yet been studied to the best of our knowledge. This article is divided into seven sections. Section one gives the brief introduction of the article. Parametric Nikiforov-Uvarov method is presented in section two. The radial solution of the proposed potential is presented in section three. Numerical computation of energy eigenvalues is carried out in section four. Thermodynamic properties are presented in section five. Results and discussion are presented in section 6, while the article is concluded in section seven.

Review of parametric Nikiforov-Uvarov method
The NU method was proposed to solve the second-order linear differential equation by reducing it to a generalized equation of hypergeometric type of the form where σ (s) and σ 2 (s) are polynomials at most second degree, and τ (s) is a first degree polynomial.
According to Tezcan and Sever [57], the parametric Nikiforov-Uvarov method is given as The total wave function is given as , while the total energy eigenequation is given as where the parametric constants can be obtained as follows: The Jacobi polynomial can be expressed in terms of Laguerre polynomial for a special case where C 3 = 0 [32] and Hence, the solution given by Eq. (3) becomes

Radial solution of Schrӧdinger equation with modified exponential screened plus Yukawa potential
The modified exponential screened plus Yukawa potential is given as , D e is the dissociation energy that describes the depth of the potential well, α is the screening parameter which characterizes the strength of the potential. D 1 is a real constant which also serve as a control parameter for the potential model, while r is the internuclear distance between the atoms of diatomic molecules.

Numerical computation of energy eigenvalues
Using Eq. (17), the numerical bound state solutions were carried out for fixed principal quantum number (n) with varying orbital angular quantum number l = 0, 1, 2 and 3.

Thermodynamic properties
In this section, we present the thermodynamic properties for the potential model. The thermodynamic properties of quantum systems can be obtained from the exact partition function given by [53,54] where is an upper bound of the vibrational quantum number obtained from the numerical solution of dE n dn = 0 , β = 1 kT where k and T are Boltzmann constant and absolute temperature, respectively. In the classical limit, the summation in Eq. (19) can be replaced with an integral: In order to obtain the partition function, energy equation (17) can be presented in a close and compact form as The energy equation of (21) can be presented in a compact form as Using Eq. (20), the partition function can be expressed as Using Maple 10.0 version, the partition function of equation (25) can be evaluated as Other thermodynamic properties can be obtained using the partition function.   (c) Vibrational Entropy: The vibrational entropy is given as where (d) Vibrational free energy: The vibrational free energy F (β) is given as       0.226 and 0.77 are spectroscopic constants for ZnH, CdH, HBr, and HF molecules, respectively, taken from Refs. [38,48] for different values of the reduced mass ( µ ) and dissociation energy ( D e ). Throughout our computation, D 1 is a control and optimizing parameter which is fixed to a value of 0.5, while takes a unitary value of 1. Figure 1a, b shows the variation of partition function ( Z( , β) ) which increases exponentially with respect to β and decreases to -0.5 with respect to . Figure 2a, b shows variation of vibrational mean energy ( U ( , β) ) with respect to β in Fig. 2a and variation with respect to as shown in Fig. 2b. In Fig. 2a, the vibrational mean energy increases with an increase in energy-dependent parameter but monotonically decreases with respect to . The variation of vibrational specific heat capacity C( , β) with respect to β is shown in Fig. 3a where there are several minimum turning points between β = 0.1 and 0.2. In Fig. 3b, the vibrational specific heat capacity with respect to reproduces the same characteristics of vibrational mean energy with respect to . In Fig. 4a, the vibrational entropy increases exponentially with respect to β but display an hyperbolic nature with respect as shown in Fig. 4b. The vibrational free energy decreases monotonically to a constant value of 0.08 with respect to β as shown in Fig. 5a but exhibits an hyperbolic nature with multiple maximum points at = −1 with respect to as shown in Fig. 5b. The trend of partition function and other thermodynamic properties obtained in this work is in excellent agreement with researches of refs [50][51][52][53][54] which affirm the accuracy of the work.

Conclusion
In this work, we study an approximate bound state solutions of Schrӧdinger wave equation with MESPYP potential model using parametric Nikiforov-Uvarov method where we obtain energy eigenequation and total unnormalized wave function expressed in terms of associated Jacobi polynomial. We obtained numerical bound state energies for five selected diatomic molecules. The resulting energy eigenequation was presented in a compact form and employed to study partition function and other thermodynamic properties. The numerical bound state solutions were carried out with a fixed control and optimizing parameter. The thermodynamic curves obtained are in excellent agreement to work of an existing literature.