Pure vibrational state energies and statistical-mechanical models for the reparameterized scarf oscillator

In this work, the reparameterized Scarf II oscillator was employed to describe the internal vibration of diatomic systems. Analytical equations for bound state pure vibrational energies and canonical partition function were obtained. The equations were used to derive statistical-mechanical models for the prediction of molar entropy, enthalpy, Gibbs free energy and constant pressure (isobaric) heat capacity of gaseous substances. The obtained model equations were used to generate numerical data on bound state energy eigenvalues and, to investigate the thermodynamic properties of the ground states chloroborane (BCl), bromine fluoride (BrF), and bromine chloride (BrCl) molecules. With the aid of the expression for molar entropy of the system, average absolute deviations obtained for the molecules are 0.1878%, 0.1267%, and 0.0586% from experimental data. The isobaric heat capacity model yields average absolute deviation of 2.1608%, 1.8601%, and 1.9805%. The results obtained are in good agreement with available literature data on gaseous molecule. The work could be applicable in the fields of molecular physics, chemical physics, solid-state physics and chemical engineering.


Introduction
Potential energy function (simply known as potential) is a mathematical model used to describe the interaction of a physical system with its environment. One of the problems of this representation is the absence of a universal potential energy function that can model every atomic and molecular interactions. Numerous versions of potential models have been proposed by chemist and physicist to account for observed atomic and molecular phenomena. The list of potential models includes the Morse potential [1], Eckart potential [2], Frost-Musulin potential [3], Rosen-Morse potential [4], Tietz potential [5], Hua potential [6], and Schiöberg potential [7] amongst others.
A potential energy function whose potential parameters are formulated in terms of the spectroscopic constants of a diatomic molecule is referred to as oscillator. The oscillator is a specialized model potential used to describe interactions in diatomic molecules. A diatomic molecule oscillator is required to satisfy the so-called Varshni conditions [8,9]. The Varshni conditions ensure that the potential parameters of an oscillator are expressed in terms of molecular constants such as equilibrium harmonic vibrational frequency (ω e ), rotational-vibrational coupling coefficient (α e ), equilibrium dissociation energy (D e ), anharmonicity constant (ω e x e ), and equilibrium bond length (r e ) [10].
The solution of Schrödinger and other wave equations of quantum mechanics can be obtained analytically or by numerical approach. Quite a number of analytical solution methods have been suggested in the literature, the Nikiforov-Uvarov (NU) method, exact and proper quantization rules, supersymmetric quantum mechanics approach (SUSYQM), asymptotic iteration method (AIM), and the recently introduced Nikiforov-Uvarov functional analysis (NUFA) method are some of the methods. Some illustrative examples where the different solution techniques are used to solve the Schrödinger equation can be found in Ref. [11][12][13][14][15][16] and the references therein. The solution of Schrödinger wave equation has been instrumental in retrieving information regarding the quantum mechanical system of interest. For instance, expectation values, information theoretic, optical, magnetic, electrical, and thermodynamic properties of substances have been investigated through eigen energy levels and eigenfunctions of wave Eqs 17-34.
The need to obtain analytical (or statistical-mechanical) models for the prediction of thermodynamic properties of gaseous substances have recently attracted much attention from the research community. The thermodynamic property of a gaseous molecule can be predicted theoretically with the aid of analytical model equation; it can also be determined by experimental procedures. The statistical-mechanical models are easy to use and are not expensive, only few molecular constants of a diatomic molecule are required to predict the thermal property of the system. In contrast, experimental method is time consuming, usually very expensive and require huge experimental task.
Different oscillator models have been employed to describe the internal vibration of diatomic molecules. Through such representation, canonical partition function is obtained which is then used to derive other useful thermodynamic models for the gaseous substance. Statistical-mechanical models such as Helmholtz free energy (F), mean thermal energy (U), entropy (S), enthalpy (H), Gibbs free energy (G), isobaric heat capacity (C p ), and constant volume heat capacity (C V ) have successfully been used to examine thermal properties of substances [35][36][37][38][39].
This paper is concerned with improved versions of hyperbolical-type oscillators. Available three parameter models are the specialized Pöschl-Teller potential [40], improved generalized Pöschl-Teller potential [41], and improved Scarf II potential energy function (ISPEF) [42]. On the other hand, existing four parameter models are the improved Pöschl-Teller potential [43], improved q-deformed Scarf II oscillator [10], and modified hyperbolical-type potential [44]. Analytical expressions for the prediction of molar entropy, enthalpy, Gibbs free energy, and isobaric heat capacity have been obtained with various formulations of the improved hyperbolical-type oscillators [10,[45][46][47][48]. However, the reparameterized Scarf oscillator has not been considered in the literature. It is against this background that this research is aimed at obtaining the pure vibrational state energies and some statistical-mechanical models for the reparameterized Scarf oscillator (RSO). The paper is organized as follows. In Section 2, potential parameters are developed for the RSO. Equation for pure vibrational state energy is obtained in Section 3. Statistical-mechanical models are derived in Section 4. Results of numerical computations are discussed in Section 5. A brief conclusion of the work is given in Section 6.

Potential parameters of the RSO
In previous study, the Varshni conditions for diatomic molecule potential were used to construct the IqSO. By employing the IqSO, analytical equations for the prediction of molar entropy and Gibbs free energy were derived [10]. However, statistical-mechanical models for the prediction of molar enthalpy and heat capacity have not been reported for the IqSO. In the present work, statistical-mechanical models for the prediction of molar entropy, enthalpy, Gibbs free energy and isobaric heat capacity are obtained using the RSO. The RSO can be obtained by subjecting the ISPEF through the transformation r → r-r 0 , where r 0 is a real constant. Using this transformation on Eq. 1 of Ref.
[42], the RSO is proposed via the following equation where V 1 , V 2 , and α are the potential parameters. The potential Eq. 1 is an oscillator if it satisfies the following constraints (Varshni conditions) where B e = ħ/(4πcμr e 2 ), ħ = h/2π, h being the Planck constant, c represents the speed of light, and f j (j = 1, 2, 3, . . .) is the jth derivative of V (r) at r = r e . Substituting Eq. 1 into Eqs 2-4 yields α π c ω e 2 μ D e 1 2 .
To obtain the parameter r 0 , we first note that f 2 = 2α 2 D e , and f 3 = -6α 3 D e tanhα (r e -r 0 ). Inserting these expressions into Eq. 5 gives 3 Pure vibrational state energy levels for the RSO In this section, the analytical equation for bound state vibrational energy levels is derived using the Nikiforov-Uvarov (NU) solution approach. To ensure continuity of the concept, a brief outline of the NU method is outlined.

An overview of the NU solution method
The NU solution approach is one of the most widely used methods of solving a second order differential equation. With the aid of a convenient transformation, the Schrödinger equation can be converted to a hypergeometric-type differential equation of the form [11].
where σ (s),σ (s) are polynomials at most of second-degree,τ (s) is a first-degree polynomial and R ]J (s) is a function of the hypergeometric-type. Writing R ]J (s) = ψ ]J (s)y ]J (s), by suitably choosing the function ψ ]J (s), expression Eq. 10 assumes the following form where τ s ( ) τ s ( ) + 2 π s ( ), π(s) is a polynomial of degree at most one, it is given by [11].
k is obtained by requiring that the quantity under the square root is a perfect square of a first degree polynomial. λ is deduced via the expression Successively differentiating Eq. 11 ] times (] = 0, 1, 2, . . .) leads to quantum condition given as

Solution of schrödinger equation with the RSO by NU method
The radial Schrödinger equation in the presence of any potential field V (r) is given as where ], and J are the vibrational and rotational quantum numbers, respectively. E ]J is the ro-vibrational energy eigenvalue of the quantum state ]J, R ]J is the corresponding radial wave function. where The substitution s = sinh α(r-r 0 ) transforms Eq. 17 to where, for compactness the following notations have been used; Comparing Eqs 10, 18 yieldsτ s, σ = 1 + s 2 ,σ −γ ] s 2 + γ 0 s − γ 1 . Substituting these expressions in Eq. 13 gives Setting the discriminant of the expression under square root to zero gives Eq. 20 leads to Thus, k is a two-valued parameter viz: k = k obtained by choosing the negative square root, and k = k + if the positive square root is chosen. Using Eq. 20 to eliminate γ 0 in Eq. 19, with some algebraic simplifications, one obtains It is evident that for each value of k, π (s) has four possible expressions: π --, π -+ , π +-, and π ++ obtained from all possible combinations of the ± signs in Eq. 22. The next task is to deduce

Molecule
Molecular state Molecular parameter Potential parameter ] with the help of expressions Eqs 12, 23, 24, gives with the help of relationships Eqs 23-25, and recalling thatτ(s) s, σ(s) 1 + s 2 , Eqs 14, 15 gives By equating expressions Eqs 26, 27, and eliminating γ ] , γ 0 , γ 1 in the resulting equation, the vibrational energy levels for the RSO is obtained as The energy spectrum defined by Eq. 28 increases with increasing ], the upper bound vibrational quantum number, ] max at which the eigen energies ceases to increase is given by E ′ ] (] max ) 0, where prime denotes derivative with respect to ν. Inserting Eq. 28 into this equation, we find ]

FIGURE 2
Variation of bound state energy eigenvalues versus vibrational quantum number for the ground state BCl, BrF, and BrCl molecules.

Some statistical-mechanical models for the RSO
The partition function is the master key connecting the microscopic model of a system with its macroscopic property. Other thermodynamic quantities can easily be expressed in terms of the partition function. Therefore, explicit expression of the partition function of a gas system is vital in this study. Statisticalmechanical models considered in this paper include molar entropy, enthalpy, Gibbs free energy, and isobaric specific heat capacity.

The canonical partition function
The canonical partition function of a gaseous molecule is composed of the translation (Q tra ), vibrational (Q vib ), and rotational (Q rot ) [49] components. The canonical partition function is given by the product Q (T) = Q tra Q vib Q rot , T being the temperature of the enclosed gas. The vibrational partition function depends on the vibrational energy levels of the diatomic system which in turn depends on the diatomic oscillator used to describe the internal vibration of the gas molecules. The vibrational partition function is given by [50].
where g ] is the factor of the degeneration of the spectrum, β = (k B T) −1 , k B is the Boltzmann constant. Given a non-degenerate system of gas molecules, g ] = 1. Eqs 28-30 leads to where for a finite series with an upper bound ] max , the Poisson series formula can be written [51].
by evaluating the definite integral, Eq. 34 yields x 0 exp(t 2 )d t is the imaginary error function of parameter x. Eq. 35 contains quantum correction terms, these are terms with j ≠ 0, contributions from these terms are significant if T is small. However, if T is large, the contributions from the quantum correction terms are relatively small. The lowest order approximation has previously been used to obtain approximate expression for the high temperature vibrational partition function of diatomic molecule oscillators [52,53]. In the lowest order approximation, only contribution from the term with j = 0 is considered, the contributions from quantum correction terms are ignored. Due to the temperature range of the diatomic systems considered in this work, with the help of the lowest order approximation model, expression Eq. 35 is reduced to If the molecules of a diatomic gas are visualized as rigid rotors, neglecting molecular interactions, the translational and rotational partition functions are given as [33,47].
where σ is assigned the value 1, 2 for heteronuclear and homonuclear diatomic molecules, respectively. Θ rot = ħ 2 /(2μr e 2 k B ) is the characteristic temperature, m is the mass of the molecules which make up gas, p is the gas pressure.

Molar entropy model for the RSO
The molar entropy for the system is evaluated from the expression [48].
where R = N A k B is the molar gas constant, N A is the Avogadro number, V is the volume of gas enclosed. Substituting Q = Q tra Q vib Q rot into the second term in Eq. 39 gives where, for brevity we have employed the representation

Molar enthalpy model for the RSO
The molar enthalpy of the system is defined by the equation [48].
with the aid of the substitution Q = Q tra Q vib Q rot , expression Eq. 42 and the relations Eqs 36-38 gives For the purpose of comparing theoretical results to experimental data, the reduced molar enthalpy is considered, it is given as [46][47][48].
where H 298.15 is the molar enthalpy computed at p = 1 bar, and T = 298.15 K. Therefore, putting Eq. 43 in Eq. 44, we have

Molar Gibbs free energy model for the RSO
The molar Gibbs free energy of the system can be deduced from the expression [48]. TABLE 6 Computed data on molar entropy (J mol −1 K −1 ), reduced enthalpy (kJ mol −1 ), reduced Gibbs free energy, and isobaric specific heat capacity for the BrCl (X 1 Σ + ) molecule. N p = 60. Frontiers in Physics frontiersin.org 08

Frontiers in Physics
frontiersin.org 09 Inserting the equation of canonical partition function in to the first term of Eq. 46, we obtain the molar Gibbs free energy in compact form as The reduced or scaled Gibbs free energy is often required for comparison with experimental data, it is given as G red = (G-H′ 298.15 )/T. Inserting Eq. 47 into this equation yields where the parameter Γ vib is given by

Results and discussion
In this section, the analytical equations developed for the reparameterized Scarf II oscillator are analyzed on some selected diatomic molecules. The experimental values of the relevant molecular constants D e , r e , ω e , α e of the ground state BCl, BrF, and BrCl molecules are taken from Refs. [54][55][56][57][58][59]. The molecular constants and computed potential parameters of the molecules are shown in Table 1. With the aid of Eq. 1, potential energy data, V (r min ) = V min , and V (r max ) = V max are obtained. Due to unavailability of literature data on BrCl molecule, only results for BCl and BrF molecules are obtained. Tables 2, 3 summarizes results of numerical computations, and available literature data on vibrational energies of the molecules. The literature data were those obtained by multireference configuration interaction (MRCI), and Rydberg-Klein-Rees (RKR) method [54,59]. Figure 1 shows graphical fitting of the RSO to (A) MRCI data points of BCl molecule, and (B) RKR data points of BrF molecule. The plots reveal that the RSO could model the internuclear potential energy curves of the BCl and BrF molecules. However, graphical plots only give an idea of the agreement between predicted data and experimental results. The average absolute deviation from experimental data (σ ave ) is one of the most widely used goodness-of-fit indicators to evaluate the accuracy of an empirical model. Previously, the average absolute deviation has been used to substantiate the accuracy of proposed model Eqs 38,41,43,48. The average absolute deviation can be written as where N p is the number of experimental data points, (X, Y, Z) ≡ (V, E, S, H, G, C p ). Average absolute deviations less than 1% shows good agreement between predicted and observed data, the smaller the σ ave , better the model. Using Eq. 52 and the data in Tables 2, 3, average absolute deviations of 5.3976% and 1.6790% are obtained from the experimental data of the BCl and BrF molecules, respectively. The results show that data predicted by the RSO are relatively high for the BCl molecule. To within 2% error limit of the RKR data, the RSO can approximately reproduced the internuclear potential energy curve of the ground state BrF molecule. Figure 2 shows the variation of energy levels Eq. 28 of the molecules against vibrational quantum number, ]. The plot reveal that as ] is gradually increased from zero, the energy of the molecules also increases. The upper bound vibrational quantum numbers obtained for the BCl, BrF and BrCl molecules are 102, 60 and 82, respectively. As ] is further increased beyond ] max , the bound state energy of the molecules begins to decrease, leading to degenerate energy levels of the molecules. Eq. 28 is also used to generate numerical data on bound state pure vibrational energies of the molecules. As there are no available literature data on energies of BrCl molecule to allow comparison, only the results for BCl and BrF molecules are displayed in Tables 2, 3. Using Eq. 52 and the energy data in Tables 2, 3, average absolute deviation obtained are 1.1949% for BCl molecule, and 1.8353% for BrF molecule. Thus, it can be seen that the energy values predicted by expression Eq. 28 are in good agreement with existing experimental data on the diatomic molecules.
To authenticate the applicability of the statistical-mechanical models proposed in this study, Eq. 43 is used to compute the value of H′ 298.15 at T = 298.15 K, p = 0.1 MPa to yield; 14.502 kJmol −1 for BCl, 13.626 kJmol −1 for BrF and 12.533 kJmol −1 for BrCl molecules. Using these results and the data in Table 1, the expressions for molar entropy Eq. 40, reduced enthalpy Eq. 45, reduced Gibbs free energy Eq. 48, and constant pressure specific heat capacity Eq. 51 are used to obtain numerical data for the diatomic molecules. The computations are carried out at p = 0.1 MPa and temperature in the range 300 K-6,000 K. The results of computations are shown in Tables 4-6. Also included in the tables are literature data on molar entropy (S NIST ), reduced molar enthalpy (H NIST ), reduced molar Gibbs free energy (G NIST ) and constant pressure specific heat capacity (C pNIST ). The data are those reported in the National Institute of Standards and Technology (NIST) database [60]. Choosing X ≡ Z ≡ S NIST to represent experimental data on molar entropy and setting Y ≡ S, average absolute deviations obtained for the BCl, BrF and BrCl molecules are 0.2011%, 0.1224% and 0.5323%, respectively. Thus, it is obvious that within error limit of 1% of the NIST data, entropy equation proposed in this work can accurately predict the experimental data on molar entropy of the diatomic molecules.
To confirm the suitability of expression Eq. 45 to predict molar entropy of the diatomic molecules, we let X ≡ Y ≡ H NIST to denote the experimental data on reduced molar entropy. Choosing Y ≡ H red . The average absolute deviation deduced for the diatomic molecules are 1.5346% for BCl, 0.688% for BrF and 0.5323% for BrCl molecules. The results reveal that statistical-mechanical model proposed for the RSO can accurately predict the NIST data on molar enthalpy of the examined molecules.
To show the relevance of expression Eq. 51 to model experimental data on molar Gibbs free energy, average absolute deviation is computed for each of the molecules. By appropriately choosing the parameters in Eq. 52 with respect to experimental data, such that X ≡ Y ≡ G NIST , and Z ≡ G red , average absolute deviations obtained for the BCl, BrF and BrCl diatomic molecules are 0.1033%, 0.0903%, and 0.0367%, respectively. The results show that on a scale of 1% of the NIST data, Gibbs free energy expression proposed by the RSO is a near perfect model to accurately predict experimental results of the examined diatomic molecules.
The analytical expression of constant pressure specific heat capacity is also analyzed for the diatomic molecules. Average absolute deviation deduced for the BCl, BrF, and BrCl molecules are 2.1565%, 1.9731%, and 1.9805%, respectively. The results are in good agreement with NIST data particularly in the low to moderate temperature regions. However, in the high temperature regime, the predicted results are relatively higher than observed data. The discrepancy arises as a result of rigid-rotor approximation of the diatomic molecules and also, the quantum correction terms excluded in the expression of vibrational partition function, which by extension are also excluded in all the statisticalmechanical models of the system.

Conclusion
In this paper, conditions to be satisfied by a diatomic molecule potential are used to construct the reparameterized Scarf II oscillator (RSO), suitable for application to diatomic molecules. Using the Nikiforov-Uvarov method to solve the radial Schrödinger equation for the RSO, analytical expression of bound state pure vibrational energy is derived for the system. With the aid of the formula for energy eigenvalues, analytical equations representing canonical partition function and other relevant statistical-mechanical models are obtained, including molar entropy, enthalpy, Gibbs Frontiers in Physics frontiersin.org free energy and isobaric specific heat capacity. The obtained equations were used to study thermodynamic properties of three diatomic molecules viz BCl, BrF, and BrCl. Average absolute deviations of 1.5364%, 0.688%, and 0.5323%, respectively are obtained using the expression of reduced molar enthalpy. The equation of reduced Gibbs free energy yields average absolute deviations of 0.1033%, 0.0903%, and 0.0367% for the diatomic molecules. The results are in excellent agreement with existing literature data on the diatomic molecules. The developed statistical-mechanical models could be useful in scientific and engineering researches involving thermochemical processes.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.