📄 chemistry

Modified Scarf-type interaction potential for vibrational energy eigenvalues and chemical equilibrium constants of selected diatomic molecules

This paper proposes a modified Scarf-type interaction potential that, when solved via the parametric Nikiforov-Uvarov method, accurately predicts vibrational energy eigenvalues and chemical equilibrium constants for selected diatomic molecules, demonstrating strong agreement with experimental and reference data.

Original authors: E. S. Eyube, M. M. Dahiru, J. E. Ngada, R. D. Musa, C. Joseph, Y. Dahiru, F. M. Sahabo

Published 2026-07-08

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: E. S. Eyube, M. M. Dahiru, J. E. Ngada, R. D. Musa, C. Joseph, Y. Dahiru, F. M. Sahabo

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine two atoms as a pair of dancers holding hands. Sometimes they spin close together, sometimes they stretch out, but they are always connected by an invisible "spring" of force. In the world of physics, scientists need to know exactly how strong that spring is and how the dancers move to predict what will happen in a chemical reaction.

This paper introduces a new, improved "map" (called a Modified Scarf-type interaction potential) to describe that invisible spring for pairs of atoms (diatomic molecules). Here is a breakdown of what the researchers did, using simple analogies:

1. The Problem: Old Maps vs. New Maps

For a long time, scientists have used various mathematical formulas to describe how atoms vibrate and interact. Think of these formulas as different types of maps. Some maps are very detailed but hard to read; others are easy to read but miss important details.

The researchers wanted a map that is both easy to read (mathematically simple) and highly accurate. They created a new map called the Modified Scarf-type (MST) potential.

2. The Solution: Solving the Dance Floor Puzzle

To test their new map, the team had to solve a complex physics equation (the Schrödinger equation) that describes the energy of the dancing atoms.

The Method: They used a specific mathematical tool called the Parametric Nikiforov-Uvarov method. You can think of this as a specialized "calculator" or "key" that unlocks the solution to the equation, turning a messy puzzle into a clean, clear formula.
The Result: They got a precise formula that tells them exactly how much energy the atoms have at different vibration levels (like how fast the dancers are spinning).

3. The Test: Does the Map Match Reality?

A map is only good if it matches the actual terrain. The team tested their new MST map against real-world data for several specific molecules (like Bromine, Phosphorus Oxide, and Sulfur Oxide).

The Comparison: They compared their mathematical predictions against two types of "gold standard" data:
1. RKR Data: This is like a high-resolution satellite photo of the atoms' energy levels.
2. NIST-JANAF Data: This is a massive, trusted library of chemical data used by engineers and scientists worldwide.
The Score: The researchers calculated how far off their map was from the real data. They found the error rate was incredibly low:
- For the energy levels of the atoms, the error was less than 3% (and as low as 0.2% for some molecules).
- For predicting chemical equilibrium (how likely atoms are to stick together or break apart), the error was less than 2% for some molecules and very close for others at higher temperatures.

4. Why This Matters (According to the Paper)

The paper claims that this new MST model is a powerful tool because:

It's Accurate: It predicts how atoms vibrate and interact with almost the same precision as complex, computer-heavy simulations, but without needing a supercomputer.
It's Useful for Chemistry: Because they have a precise formula for how the atoms move, they can easily calculate thermodynamic properties (like heat and energy) and chemical equilibrium constants (which tell us if a reaction will happen).
It Works Well: The model successfully reproduced the behavior of both single-atom pairs (like Nitrogen-Nitrogen) and mixed pairs (like Nitrogen-Oxygen).

Summary

In short, the authors built a new, simpler, and highly accurate mathematical "ruler" for measuring how diatomic molecules vibrate and react. They proved this ruler works by showing it matches real-world measurements almost perfectly. This allows scientists to predict chemical behavior and energy properties more easily and reliably than before.

Technical Summary: Modified Scarf-type Interaction Potential for Vibrational Energy Eigenvalues and Chemical Equilibrium Constants

Problem Statement
Accurate determination of equilibrium thermodynamic properties and chemical equilibrium constants is critical for understanding chemical and physicochemical processes. However, existing computational approaches, such as ab initio, density functional theory, and coupled-cluster methods, often require extensive electronic-structure calculations and highly accurate vibrational frequencies. These methods can be computationally expensive and may suffer from numerical approximations, particularly for weakly interacting systems. While analytically solvable molecular interaction potentials (e.g., Morse, Kratzer, Pöschl-Teller) offer efficient alternatives for deriving vibrational spectra and thermodynamic quantities, systematic investigations utilizing modified Scarf-type frameworks for equilibrium thermodynamic formulations and chemical equilibrium constants remain limited in the literature.

Methodology
The study proposes a Modified Scarf-type (MST) interaction potential to address these gaps. The methodology proceeds through the following analytical and computational steps:

Potential Definition: A modified Scarf-type potential, $U(\rho)$ , is defined using hyperbolic functions involving strength parameters ( $U_0, U_1, U_2$ ) and a screening parameter ( $\alpha$ ). The potential parameters are constrained by molecular boundary conditions, ensuring the model satisfies physical requirements for dissociation energy ( $D_e$ ), equilibrium bond length ( $\rho_e$ ), stiffness constant ( $k$ ), and rotation-vibration coupling ( $\alpha_e$ ).
Analytical Solution: The radial Schrödinger equation for the MST potential is solved analytically using the Parametric Nikiforov-Uvarov (PNU) method. This approach transforms the differential equation into a hypergeometric-type form to derive an explicit expression for the vibrational energy eigenvalues ( $E_n$ ).
Thermodynamic Formulation: The derived energy eigenvalues are used to construct the vibrational partition function ( $Q_A$ ). To improve agreement with reference data, a fitting parameter ( $\xi$ ) is introduced at the partition function level. The total molecular partition function is obtained by combining vibrational, rotational, and translational components.
Property Derivation: Using the total partition function, the study derives expressions for molar enthalpy ( $H$ ), Gibbs free energy ( $G$ ), and the chemical equilibrium constant ( $K_p$ ) for reversible reactions forming heteronuclear diatomic molecules.
Validation: The model is validated against experimental and reference data for ten diatomic molecules (Br $_2$ $_{2}$ , F $_2$ $_{2}$ , N $_2$ $_{2}$ , O $_2$ $_{2}$ , P $_2$ $_{2}$ , S $_2$ $_{2}$ , BrF, NO, PO, SO). Validation metrics include:
- Comparison of potential energy curves and vibrational energy eigenvalues against Rydberg-Klein-Rees (RKR) data.
- Comparison of chemical equilibrium constants against NIST-JANAF thermochemical data.
- Calculation of Mean Percentage Absolute Deviation (MPAD) to quantify accuracy.

Key Contributions and Results
The paper presents the following specific findings:

Energy Eigenvalues: The analytical expression for vibrational energy eigenvalues successfully reproduces experimental RKR data. For Br $_2$ , PO, and SO, the MPAD values for vibrational energies are 2.669%, 0.547%, and 1.86%, respectively.
Potential Energy Surfaces: The MST potential accurately models the potential energy profiles. MPAD values for the potential energy curves are 2.009% (Br $_2$ ), 0.245% (PO), and 2.972% (SO), with PO showing the highest agreement.
Chemical Equilibrium Constants: The model demonstrates strong capability in predicting chemical equilibrium constants ( $K_p$ $K_{p}$ ) across the temperature range of 300–6000 K.
- For BrF and NO, the theoretical results closely match NIST-JANAF data with MPAD values of 1.33% and 1.93%, respectively.
- For PO and SO, deviations are observed at low temperatures, but the model converges with reference data at moderate and high temperatures, maintaining favorable agreement.
Thermodynamic Consistency: The study establishes a consistent framework linking the MST interaction potential to thermodynamic quantities (enthalpy, Gibbs free energy) and equilibrium constants without relying on high-cost electronic-structure calculations.

Significance and Claims
The authors claim that the Modified Scarf-type interaction potential serves as a useful and reliable framework for molecular spectroscopy, thermodynamic modeling, and chemical equilibrium studies of diatomic molecules. The work demonstrates that analytically tractable models can provide accurate descriptions of vibrational spectra and equilibrium properties while avoiding the computational costs associated with large-scale electronic-structure calculations. The study concludes that the MST potential is particularly effective for systems where moderate to high-temperature equilibrium behavior is of interest, offering a transparent theoretical approach for investigating molecular interactions and thermodynamic properties.

1. The Problem: Old Maps vs. New Maps

2. The Solution: Solving the Dance Floor Puzzle

3. The Test: Does the Map Match Reality?

4. Why This Matters (According to the Paper)

Summary

More like this