Abstract

In this work, we apply the parametric Nikiforov-Uvarov method to obtain eigensolutions and total normalized wave function of Schrödinger equation expressed in terms of Jacobi polynomial using Coulomb plus Screened Exponential Hyperbolic Potential (CPSEHP), where we obtained the probability density plots for the proposed potential for various orbital angular quantum number, as well as some special cases (Hellmann and Yukawa potential). The proposed potential is best suitable for smaller values of the screening parameter . The resulting energy eigenvalue is presented in a close form and extended to study thermal properties and superstatistics expressed in terms of partition function and other thermodynamic properties such as vibrational mean energy , vibrational specific heat capacity , vibrational entropy , and vibrational free energy . Using the resulting energy equation and with the help of Matlab software, the numerical bound state solutions were obtained for various values of the screening parameter () as well as different expectation values via Hellmann-Feynman Theorem (HFT). The trend of the partition function and other thermodynamic properties obtained for both thermal properties and superstatistics were in excellent agreement with the existing literatures. Due to the analytical mathematical complexities, the superstatistics and thermal properties were evaluated using Mathematica 10.0 version software. The proposed potential model reduces to Hellmann potential, Yukawa potential, Screened Hyperbolic potential, and Coulomb potential as special cases.

1. Introduction

The approximate analytical solutions of one-dimensional radial Schrödinger equation with a multiple potential function have been studied using a suitable approximation scheme to the centrifugal term within the frame work of the parametric Nikiforov-Uvarov method [1]. The solutions to the wave equations in quantum mechanics and applied physics play a crucial role in understanding the importance of physical systems [2]. The two most important parts in studying Schrödinger equations are the total wave function and energy eigenvalues [3]. The analytic solutions of wave equations for some physical potentials are possible for . For , special approximation schemes like the Greene-Aldrich and Pekeris approximations are employed to deal with the centrifugal barrier in order to obtain approximate bound state solutions [46]. The Greene-Aldrich approximation scheme is mostly applicable for short range potentials [7]. Eigensolutions for both relativistic and nonrelativistic wave equations have been studied with different methods which include the following: Exact quantisation, WKB, Nikiforov-Uvarov method (NU), Laplace transform technique, asymptotic iteration method, proper quantisation, supersymmetric quantum mechanics approach, vibrational approach, formula method, factorisation method, and Shifted 1/N-expansion method [813]. Bound state solutions obtained from the Schrödinger equation has practical applications in investigating tunnelling rate of quantum mechanical systems [14] and mass spectra of quarkonia systems [1519]. Among other goals achieved in this research article is to apply the Hellmann-Feynman Theorem (HFT) to eigenequation of the Schrödinger wave equations to obtain expectation values of , , , and analytically. The Hellmann-Feynman Theorem gives an insight about chemical bonding and other forces existing among atoms of molecules [2025]. To engage HFT in calculating the expectation values, one needs to promote the fixed parameter which appears in the Hamiltonian to be a continuous variable in order to ease the mathematical purpose of taking the derivative [26]. Similarly, the application of Hellmann-Feynman Theorem provides a less mathematical approach of obtaining expectation values of a quantum mechanical systems [27, 28]. Some of the potential models considered within the framework of relativistic and nonrelativistic wave equations are Hulthen-Yukawa Inversely quadratic potential [29], noncentral Inversely quadratic potential [30], Modified Hylleraas potential [31], Yukawa, Hulthen, Eckart, Deng-Fan, Pseudoharmonic, Kratzer, Woods-Saxon, double ring shape, Coulomb, Tietz–Wei, Tietz-Hua, Deng-Fan, Manning-Rosen, trigonometric Rosen-Morse, hyperbolic scalar, and vector potential and exponential type potentials among others [3250]. Coulomb, hyperbolic, and screened exponential type potentials have been of interest to researchers in recent times because of their enormous applications in both chemical and physical sciences. In view of this, Parmar [51] studied ultrageneralized exponential hyperbolic potential where he obtained energy eigenvalues, unnormalized wave function and the partition function. This potential reduces to Yukawa potential, Screened cosine Kratzer potential, Manning-Rosen potential, Hulthen plus Inversely quadratic exponential Mie-type potential, and many others. Diaf et al. [52], in their studies, obtained eigensolutions to the Schrödinger equation with trigonometric Inversely quadratic plus Coulombic hyperbolic potential where they obtained energy eigenvalue and normalized wave function using the Nikiforov-Uvarov method. Onate [53] examined bound state solutions of the Schrödinger equation with second Pöschl-Teller-like potential where he obtained vibrational partition function, mean energy, vibrational specific heat capacity, and mean free energy. In that work, the Pöschl-Teller-like potential was expressed in a hyperbolic form. The practical application of energy eigenvalue of Schrödinger equation in investigating the partition function, thermodynamic properties, and superstatistics arouses the interest of many researchers. Recently, Okon et al. [54] obtained the thermodynamic properties and bound state solutions of the Schrödinger equation using Mobius square plus screened Kratzer potential for two diatomic systems (carbon(II) oxide and scandium fluoride) within the framework of the Nikiforov-Uvarov method. Their results were in agreement to semiclassical WKB among others. They presented energy eigenvalue in a close form in order to obtain partition function and other thermodynamic properties. Omugbe et al. [55] recently studied the unified treatment of the nonrelativistic bound state solutions, thermodynamic properties, and expectation values of exponential-type potentials where they obtained the thermodynamic properties within the framework of semiclassical WKB approach. The authors studied the special cases of the potential as Eckart, Manning-Rosen, and Hulthen potentials. Besides, Oyewumi et al. [56] studied the thermodynamic properties and the approximate solutions of the Schrödinger equation with shifted Deng-Fan potential model within the framework of asymptotic Iteration method where they apply Pekeris-type approximation to centrifugal term to obtain rotational-vibrational energy eigenvalues for selected diatomic systems. A lot of researches have been carried out by Ikot et al. These can be seen in Refs. [5760]. Also, Boumali and Hassanabadi [61] studied thermal properties of a two-dimensional Dirac oscillator under an external magnetic field where they obtained relativistic spin-1\2 fermions subject to Dirac oscillator coupling and a constant magnetic field in both commutative and noncommutative spaces.

In this work, we propose a novel potential called Coulomb plus Screened Exponential Hyperbolic Potential to study bound state solutions, expectation values, superstatistics, and thermal properties within the framework of the parametric Nikiforov-Uvarov method [62]. This article is divided into 9 sections. The introduction is given in Section 1. The parametric Nikiforov-Uvarov method is presented in Section 2. The solutions of the radial Schrödinger equation are presented in Section 3. The application of Hellmann-Feynman Theorem to obtain expectation values is presented in Section 4. The thermodynamic properties and superstatistics formulations are presented in Sections 5 and 6, respectively. Numerical results and discussion are presented Sections 7 and 8, respectively, and the article is concluded in Section 9.

The propose Coulomb plus Screened Hyperbolic Exponential Potential (CPSHEP) is given as where and are the potential depths, is a real constant parameter, and is the adjustable screening parameter. The Pekeris-like approximation to the centrifugal term is given as

The graph of Pekeris approximation to centrifugal term is given in Figure 1.

2. Parametric Nikiforov-Uvarov (NU) Method

The NU method is based on reducing second order linear differential equation to a generalized equation of hypergeometric type and provides exact solutions in terms of special orthogonal functions like Jacobi and Laguerre as well as corresponding energy eigenvalues [6370]. The reference equation for parametric NU method according to Tezcan and Sever [71] is given as

The condition for energy equation is given as [70].

The total wave function is given as

The parametric constants are obtained as follows:

3. The Radial Solution of Schrödinger Wave Equation

The radial Schrödinger wave equation with the centrifugal term is given as

Equation (7) can only be solved analytically to obtain exact solution if the angular orbital quantum number . However, for , equation (7) can only be solve by using the approximations in (2) to the centrifugal term. Substituting equation (1) into (7) gives

By substituting equation (2) into (8) gives the following equation:

By defining and with some simple algebraic simplification, equation (9) can be presented in the form where

Comparing equation (10) to (3), the following polynomials were obtain:

Using equation (6), other parametric constants are obtained as follows:

Using equations (4), (12), and (13) with much algebraic simplification, the energy eigenvalue for the proposed potential is given as

Using equation (5), the total unnormalized wave function is given as

To obtain the normalization constant of equation (15), we employ the normalization condition where

The wave function is assumed to be in bound at and .

Equation (15) reduces to

Let such that the boundary of integration of equation (19) changes from to . Then, equation (19) reduces to

Using the standard integral,

Let . Then, using equation (20), the normalization constant can be obtained as

Hence, the total normalized wave function is given as

4. Expectation Values Using Hellmann-Feynman Theorem

In this section, some expectation values are obtain using Hellmann-Feynman Theorem (HFT). According to the Hellmann-Feynman Theorem, the Hamiltonian for a particular quantum mechanical system is expressed as a function of some parameters . Let and be the eigenvalues and eigenfunction of the Hamiltonian. Then,

For the purpose of clarity, the Hamiltonian for the propose potential using HFT is

4.1. Expectation Value for

To obtain the expectation value for , we set to have where

4.2. Expectation Value for

Setting , then, the expectation value for becomes

4.3. Expectation Value for

Setting , then,

4.4. Expectation Value for

The relationship between and is given as , therefore

5. 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 where is an upper bound of the vibrational quantum number obtained from the numerical solution of , given as , , where and are Boltzmann constant and absolute temperature, respectively. In the classical limit, the summation in equation (31) can be replaced with an integral:

In order to obtain the partition function, the energy equation (14) can be presented in a close and compact form as

Equation (33) can further be simplified to where and equation (34) can be represented in the form where (i)Partition function is obtained as follows:

Substituting equation (36) into equation (32) taking note of changes in the integration boundaries using equation (37) gives the partition function as

Using Mathematica 10.0 version, the partition function of equation (38) is given as (ii)Vibrational mean energy is given as(iii)Vibrational entropy is given as(iv)Vibrational free energy is given as(v)Vibrational specific heat capacity is given aswhere

6. Superstatistics Formulation

Superstatistics is the superposition of two different statistics which is applicable to driven nonequilibrium systems to statistical intensive parameter () fluctuation [72]. This intensive parameter which undergoes spatiotemporal fluctuations includes chemical potential and energy fluctuation which is basically describe in terms of effective Boltzmann factor [73]. According to Edet et.al [74], the effective Boltzmann factor is given as where is the Dirac delta function. However, the generalized Boltzmann factor expressed in terms of deformation parameter is given as

The partition function for superstatistics formalism is then given as

Substituting equation (34) into equation (46) gives the generalized Boltzmann factor equation as

Using equation (47), the superstatistics partition function equation is given as

Using Mathematica 10.0 version, the partition obtain from equation (47) is

We use the same procedure of thermodynamics (Section 5) to obtain superstatistics vibrational mean energy (), vibrational specific heat capacity (), vibrational entropy (), and vibrational free energy () from the partition equation (41). However, this solution is not included in the article because of the lengthy and bulky analytical equations.

7. Numerical Results

Using Matlab 10.0 version, the numerical bound state solutions for the proposed potential were calculated using (14) for different quantum state. Also, using equations (18), (20), (21), and (22), the expectation values for , , , and , respectively, were calculated as shown in Tables 15.

7.1. Special Cases

(a)Hellmann potential: substituting into equation (1), then, the potential reduces to Hellmann potential

The required energy equation is

The result of equation (52) agrees excellently with Ref. [75]. (b)Yukawa potential: if , then, equation (1) reduces to Yukawa potential

The corresponding energy eigenequation is

Equation (54) is consistent with result obtain in Ref. [76] in order to prove the mathematical accuracy of our analytical calculation. (c)Screened-hyperbolic inversely quadratic potential: substituting into equation (1), the potential reduces to screened-hyperbolic inversely quadratic potential

The resulting energy eigenvalue is (d)Coulomb potential: substituting into equation (53), then, the potential reduces to Coulomb potential

By substituting into equation (54), it gives the corresponding energy eigenvalue for Coulomb’s potential as

Equation (58) agrees we with result obtain in Ref [76].

8. Discussion

Figure 1 is the graph of Pekeris approximation against the screening parameter . The nature of the graph shows that the approximation is suitable for the proposed potential. Variation of the probability density against the internuclear separation at various quantum state for and , respectively, is shown in Figures 2(a) and 2(b), respectively. In Figure 2(a), the probability density curves produce thermal curves with regular peaks compacted close to the origin with uneven peaks at various internuclear distance. This curve shows that for orbital angular quantum number there is more concentration of the electron density at the origin for all the quantum state studied. The same situation is also observed for as shown in Figure 2(b). It can also be observed that at every value of the internuclear distance, the probability density for is higher than the probability density for . Figure 3 shows variation of the probability density against the internuclear separation at various quantum state for and for Hellmann potential presented in Figures 3(a) and 3(b), respectively. More concentration of the electron density is observed at the origin in both cases. It is also seen that the probability density obtained for is lower than the probability density obtained for In Figure 4, we presented the variation of the probability density against the internuclear separation at various quantum state for and for Yukawa potential as shown in Figures 4(a) and 4(b), respectively. Here, there are more concentration and localization of electron density at the origin for with uneven distribution of peak curves as shown in Figure 4(a), but this situation is not the same for as presented in Figure 4(b). The probability density at the second excited state for remains constant with uniform distribution of the density curves for all values of the internuclear separation.

In Figure 5, the variation of partition function for thermodynamic properties and superstatistics with the temperature parameter were observed. Here, the partition function increases nonlinearly with . However, in Figure 5(a), the partition function diverged as increases, but later converges for the superstatistics as shown in Figure 5(b). It can also be observed that the partition function converges as becomes positive in the superstatistics in Figure 5(b). In Figure 6, we presented the variation of vibrational mean energy against for thermodynamic properties and superstatistics, respectively, as shown in Figures 6(a) and 6(b). For the thermal property, the mean energy increases as goes up for all values of as shown in Figure 6(a), and at higher values of the mean energy for various tends to converge. The superstatistics mean energy rises as the temperature of the system decreases as shown in Figure 6(b). However, at a certain absolute temperature, the superstatistics mean energy increases vertically for various values of the deformed parameter (). The variation is opposite in Figure 7 which is the variation of the heat capacity against for different values of and deformed parameter is shown in Figures 7(a) and 7(b), respectively. In each case, the heat capacity decreases monotonically with an increasing for the nonsuperstatistics (thermodynamic property). At zero value of , the heat capacity for various converged and diverges as increases gradually. For the superstatistics, the heat capacity for various deformed parameter rises while the temperature cools down. The specific heat capacity has a turning point when equals -150 as shown in Figure 7(b) while the specific heat capacity for the nonsuperstatistics has a maximum turning point at about 0.15 as shown in Figure 7(a). In Figure 8, the vibrational entropy decreases and diverged while the temperature of the system decreases ( increases) in the nonsuperstatistics as shown in Figure 8(a). This decrease is sharper for negative values of the entropy when is almost constant as shown in Figure 8(b). The superstatistics entropy varies inversely with (directly with temperature). This means that when the temperature of the system is raised, the disorderliness of the system also increases for every value of the deformed parameter. The entropy for the superstatistics converged as the temperature parameter tends to zero. In Figure 9, the variation of free energy against the is seen to be two different steps in the case of a nonsuperstatistics as shown in Figure 9(a). Here, between 0 and 60 values of β, the free energy increases steadily for the various values of but beyond this range; the free energy increases sharply at constant . For the superstatistics, the free energy increases monotonically as the temperature of the system reduces gradually. The free energy is always higher when the deformed parameter is increased as shown in Figure 9(b).

Table 1 is the numerical bound state solutions for CPSEHP. The numerical bound state solutions increases with an increase in quantum state, but decreases with an increase in the screening parameter. Tables 25 are expectation values for , , , and , respectively. Here, the numerical values in all cases decrease with an increase in the screening parameter.

9. Conclusion

In this work, we apply the parametric Nikiforov-Uvarov method to obtain the bound state solutions of Coulomb plus Screened Exponential Hyperbolic Potential. The resulting energy eigenvalues were presented in a close and compact form. The research work was extended to study thermal properties, superstatistics, and various expectation values. The proposed potential also reduced to Hellmann potential, Yukawa potential, Screened Hyperbolic potential, and Coulomb potential as special cases. The normalized wave function for the mother potential and that of the Hellmann potential are similar, but the normalized wave function of the Yukawa potential seams different. The trend of the thermodynamic and superstatistics curves is in agreement to the results of an existing literature. The results of the thermodynamic properties and superstatistics revealed that the effect of the temperature on the thermodynamic properties and the superstatistics are similar. Finally, this research work has practical applications in physical and chemical sciences.

Data Availability

The numerical data used for this work are generated using Matlab programme.

Disclosure

An earlier version of this article has been deposited in arxiv in Cornell University which is accessible at https://arxiv.org/abs/2110.09896.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The article processing charge will be funded by SCOAP3.