Bound state solutions with a linear combination of Yuakawa plus four-parameter diatomic potentials using path integral approach: Thermodynamic properties

The Gist

Imagine you are trying to understand how two atoms in a molecule hold hands and dance around each other. In the world of quantum physics, this dance is governed by invisible forces and specific rules. This paper is like a detailed map that the authors drew to predict exactly how these atoms move, how much energy they have, and how they behave when the temperature changes.

Here is a simple breakdown of what they did, using everyday analogies:

1. The Problem: A Complicated Dance Floor

In quantum physics, scientists use math equations (like the Schrödinger equation) to describe how particles move. Usually, they look at one specific "force field" or potential at a time. However, real molecules are messy. The force between two atoms isn't just one simple thing; it's a mix of different forces.

The authors decided to study a specific "dance floor" created by mixing two different types of forces:

• The Yukawa Potential: Think of this as a force that gets weaker very quickly as you move away, like a magnet that stops working once you pull it a few inches away.
• The Four-Parameter Potential: This is a more complex force that acts like a custom-made track with specific bumps and dips.

They combined these two into a single, complicated mathematical shape to see how a molecule behaves on this mixed track.

2. The Tool: The "Path Integral" Approach

To solve the math, the authors used a method called the Path Integral approach.

• The Analogy: Imagine you are at a train station and want to get to a destination. A standard map shows you the shortest, straight line. But in the quantum world, a particle doesn't just take one path; it takes every possible path at the same time—some straight, some winding, some looping.
• The authors used this method to add up all these infinite possibilities to find the most likely outcome. It's like calculating the average of every possible route a traveler could take to find the true nature of the journey.

3. The Hurdle: The "Centrifugal" Spin

There was a tricky part of the math called the "centrifugal term."

• The Analogy: Imagine a child spinning on a merry-go-round. If they spin too fast, they want to fly off. In atoms, if the electron or nucleus has "angular momentum" (it's spinning or orbiting), it creates a force trying to push it away from the center.
• This force made the math impossible to solve exactly. So, the authors used a clever approximation (a smart guess) to simplify this spinning force so it looked like the rest of the track. This allowed them to solve the puzzle.

4. The Results: The Energy Map and the Wave

Once they solved the math, they found two main things:

• The Energy Spectrum: This is like a ladder. The atoms can only stand on specific rungs of the ladder, not in between. The authors calculated exactly how high each rung is. They found that the height of these rungs changes depending on how "stretched" or "squished" the molecule is (controlled by parameters like the screening parameter $\alpha$ and deformation parameter $q$).
• The Wave Functions: These describe the "shape" of the atom's dance. The authors figured out the exact shape of the dance for every rung on the ladder.

5. The Heat: Thermodynamics

After mapping the energy levels, they asked: "What happens when we heat this molecule up?"

• They calculated the Partition Function, which is essentially a scorecard that tells you how many different ways the molecule can vibrate at a certain temperature.
• From this scorecard, they derived other properties:
  • Free Energy: How much "work" the molecule can do.
  • Heat Capacity: How much heat the molecule can soak up before getting hotter.
  • Entropy: A measure of disorder or chaos. As the molecule gets hotter, it vibrates more wildly, increasing its chaos.

6. Testing the Theory: Real Molecules

To make sure their math wasn't just theory, they plugged in real numbers for actual molecules like Hydrogen ($H_2$), Carbon Monoxide ($CO$), and Iodine ($I_2$).

• They found that for heavy molecules (like Iodine), the energy levels are very close together, like steps on a staircase that are barely visible.
• For lighter molecules (like Hydrogen), the steps are wider apart.
• They also discovered that changing the "shape" of the force (the deformation parameter) changes the energy levels, but the effect is different for different molecules. For example, the force affects Hydrogen and Iodine very differently.

Summary

In short, this paper is a mathematical recipe. The authors mixed two different force models, used a complex "sum-of-all-paths" technique to solve the resulting equation, and created a new map of energy levels and heat behaviors for diatomic molecules. They then tested this map against real-world molecules to show that their recipe works and gives consistent results.

Technical Summary

Technical Summary: Bound State Solutions with a Linear Combination of Yukawa plus Four-Parameter Diatomic Potentials using Path Integral Approach

Problem Statement
The paper addresses the challenge of obtaining approximate analytical bound state solutions for a non-relativistic quantum system governed by a linear combination of the Yukawa potential and a generalized four-parameter diatomic potential. The specific potential under investigation is given by:
$$V(r) = \frac{a}{(e^{2\alpha r} - q)^2} - \frac{b}{e^{2\alpha r} - q} - \frac{ce^{-\alpha r}}{r}$$
where $a, b, c$ are constants related to the potential depth ($D_e$) and equilibrium distance ($r_e$), while $q$ and $\alpha$ represent deformation and screening parameters, respectively. A primary difficulty in solving the radial Schrödinger equation for such potentials arises from the centrifugal term ($l(l+1)/r^2$) for non-zero angular momentum ($l \neq 0$), which prevents exact analytical solutions due to the mismatch between the functional forms of the potential and the centrifugal barrier.

Methodology
The authors employ the Feynman path integral formalism to derive the energy spectrum and wave functions. The methodology proceeds through the following steps:

1. Approximation of the Centrifugal Term: To overcome the analytical intractability of the centrifugal term, the authors introduce a specific approximation scheme valid for $q \geq 1$:
   $$\frac{1}{r} \approx \frac{2\alpha e^{-\alpha r}}{1 - qe^{-2\alpha r}}, \quad \frac{1}{r^2} \approx \frac{4\alpha^2 e^{-2\alpha r}}{(1 - qe^{-2\alpha r})^2}$$
   This transforms the Hamiltonian into a form compatible with the potential's structure.

2. Path Integral Formulation: The Green's function is expressed as a path integral. Due to singularities in the transformed potential at $r_0 = \frac{1}{2\alpha} \ln q$, the authors restrict the analysis to the region $r > r_0$. A regulating function $f(r)$ is introduced to handle the singularity and define a discrete path integral action.

3. Coordinate Transformation: A change of variables is performed to map the problem onto a known solvable model. Specifically, the radial coordinate is transformed using a relation involving deformed hyperbolic functions (introduced by Arai), converting the effective potential into a Rosen-Morse potential (a generalization of the modified Pöschl-Teller potential).

4. Derivation of Spectral Properties: The radial Green's function for the resulting Rosen-Morse potential is identified. The energy spectrum is extracted from the poles of the Gamma function appearing in the Green's function's explicit expression. The normalized wave functions are derived from the residues at these poles, expressed in terms of Wigner functions and subsequently converted to hypergeometric functions.

5. Thermodynamic Analysis: Using the derived energy spectrum, the vibrational partition function is calculated via the Poisson summation formula. This allows for the analytical derivation of thermodynamic properties, including vibrational free energy, mean energy, heat capacity, and entropy.

Key Contributions and Results

• Analytical Solutions: The paper provides explicit analytical expressions for the energy eigenvalues ($E_{n,l}$) and the normalized radial wave functions ($\chi_{n,l}(r)$) for the combined Yukawa and four-parameter potential. The energy spectrum is shown to depend on the principal quantum number $n$, angular momentum $l$, and the potential parameters ($\alpha, q, a, b, c$).
• Thermodynamic Formulas: The authors derive closed-form expressions for the vibrational partition function and subsequent thermodynamic quantities. These formulas involve complex error functions ($\text{erfi}$) and depend on the maximum vibrational quantum number $\lambda_{\text{max}}$.
• Numerical Validation: The theoretical framework is applied to several diatomic molecules ($H_2, I_2, LiH, CO, HCl, NO$). Numerical calculations demonstrate the behavior of energy levels as functions of the screening parameter $\alpha$, deformation parameter $q$, and quantum numbers.
  • Results indicate that energy levels decrease as the quantum number $n$ or the deformation parameter $q$ increases.
  • The influence of $q$ varies by molecule; for massive molecules like $I_2$, the influence is negligible, whereas for lighter molecules like $H_2$ and $HCl$, the deformation parameter significantly affects the energy variation.
• Thermodynamic Behavior: The study illustrates the behavior of thermodynamic functions against the inverse temperature $\beta$.
  • The partition function increases with $\beta$ and $\lambda_{\text{max}}$ but decreases with increasing $q$.
  • The vibrational heat capacity ($C_{\text{vib}}$) exhibits a maximum value, separating low-temperature (increasing) and high-temperature (decreasing) regimes. The position of this maximum varies depending on the specific molecular parameters.

Significance and Claims
The paper claims to successfully extend the path integral approach to a complex linear combination of diatomic potentials that had not been treated in this specific manner. By utilizing the approximation for the centrifugal term and the path integral formalism, the authors provide a unified method to obtain bound state solutions and their associated thermodynamic properties.

The authors emphasize that their numerical results for various diatomic molecules demonstrate the coherence of the model. They highlight that the model captures the sensitivity of energy levels and thermodynamic properties to the deformation parameter $q$ and the screening parameter $\alpha$. The work serves as a theoretical tool for understanding the vibrational characteristics of diatomic molecules under these specific potential interactions, offering explicit formulas that can be used to assess thermodynamic properties without resorting to purely numerical diagonalization of the Hamiltonian. The paper concludes by noting that the Green's function for this potential cannot be evaluated in a unified manner for any deformation parameter without the specific approximations employed, validating the necessity of their approach.