Superelastic Heating in Treanor-Gordiets Plasmas: A Unified Analytic Closure

This paper presents a unified, thermodynamically consistent analytic closure that corrects conventional harmonic models' mispredictions of superelastic electron heating in non-equilibrium Treanor-Gordiets plasmas by introducing an anharmonic correction factor to accurately capture the kinetic competition between vibrational up-pumping and relaxation.

The Gist

Imagine a bustling dance floor inside a plasma (a super-hot, electrically charged gas). On this floor, you have two main groups of dancers: electrons (tiny, fast, energetic particles) and molecules (larger, slower particles that can vibrate like springs).

Usually, these two groups dance separately. But sometimes, they bump into each other. When a fast electron hits a vibrating molecule, two things can happen:

1. Cooling: The electron gives energy to the molecule, making the molecule vibrate faster and the electron slow down.
2. Superelastic Heating: The electron hits a molecule that is already vibrating wildly. The molecule gives its energy back to the electron, making the electron zoom even faster. This is the "heating" the paper focuses on.

The Problem: The "Perfect Spring" Mistake

For a long time, scientists modeled these molecules as perfect springs. They assumed that no matter how much a molecule vibrates, the energy steps between "low vibration" and "high vibration" are always exactly the same size.

The author of this paper says: "That's wrong."

Real molecules are more like rubber bands. As you stretch a rubber band further, it gets harder to stretch, and the energy steps between vibrations change.

• When the gas is cold but the molecules are vibrating wildly (a common situation in plasma engines or combustion), the "rubber band" effect causes a massive pile-up of molecules in the high-energy states.
• The old "perfect spring" models missed this pile-up. They thought there were fewer high-energy molecules than there actually were.
• The Result: The old models predicted that the electrons would get heated up much less than they actually do. In some cases, they were off by a factor of five. It's like trying to predict how much a car will speed up by ignoring that the driver is actually pressing the gas pedal harder than you thought.

The Solution: A New "Unified" Rulebook

The author created a new mathematical formula (a "closure") that fixes this. Think of it as upgrading the dance floor's rulebook to account for the "rubber band" nature of molecules.

This new rulebook does three clever things:

1. It tracks the "Rubber Band" distortion: It calculates exactly how the energy steps change as molecules vibrate harder.
2. It finds the "Traffic Jam" (The Treanor Minimum): In these plasmas, molecules pile up at a certain high-energy level before they start falling back down. The new math finds exactly where this traffic jam happens.
3. It balances the books: It ensures that if the system reaches a perfect, calm equilibrium (where everything is the same temperature), the heating and cooling cancel out perfectly, obeying the laws of thermodynamics.

The "Magic" Shortcut

Calculating every single collision between every single electron and every single molecule is like trying to count every grain of sand on a beach. It's too slow for computer simulations of real-world engines or spacecraft.

The author didn't just fix the math; they found a shortcut.

• Instead of tracking every single grain of sand, they created a "representative average" grain.
• By using this average, they reduced the amount of computer work needed by 40 to 70 times.
• This means scientists can now run fast, accurate simulations of complex systems (like plasma-assisted combustion or hypersonic flight) without needing supercomputers to do the heavy lifting.

Why This Matters (According to the Paper)

The paper specifically mentions that this new model helps us understand and predict what happens in:

• Hypersonic flight: When spacecraft re-enter the atmosphere and create shockwaves.
• Plasma-assisted combustion: Using plasma to help engines burn fuel more efficiently.
• Laser-induced plasmas: Creating plasma with lasers for various industrial or scientific uses.

In short, the paper says: "We found a way to stop underestimating how much energy electrons gain from vibrating molecules. We fixed the math to account for real-world 'rubber band' physics, and we made it fast enough to use in real engineering simulations."

Technical Summary

Technical Summary: Superelastic Heating in Treanor-Gordiets Plasmas

Problem Statement
In thermally non-equilibrium plasmas, conventional harmonic models for electron-vibrational (e-V) energy coupling significantly mispredict superelastic heating rates ($Q_{v-e}$). These models fail when the vibrational temperature ($T_v$) exceeds the gas temperature ($T_g$), a condition common in the inter-pulse stages of plasma-assisted combustion (PAC), hypersonic re-entry wakes, and non-equilibrium chemical synthesis discharges. Under these conditions, the vibrational manifold follows a Treanor distribution characterized by the overpopulation of high-lying states. Conventional harmonic models, which assume equal energy spacing between levels, neglect the exponential heating contribution from these overpopulated anharmonic states, leading to underestimations of heating rates by factors of approximately five. Conversely, when $T_g > T_v$, harmonic models overestimate heating by roughly 15%. Furthermore, existing thermodynamically consistent models have been limited to low-temperature regimes or mono-quantum transitions, failing to capture the multi-quantum overtone transitions and anharmonic distortions present in high-energy discharges.

Methodology
To resolve these discrepancies, the author derives a unified, thermodynamically consistent macroscopic closure based on the principle of detailed balance and a second-order Dunham expansion. The methodology proceeds through the following steps:

1. Anharmonic Energy Formulation: The vibrational energy levels are described using a truncated second-order Dunham expansion. This yields an analytical expression for the energy gap ($\Delta E_{n,m}$) between vibrational levels $n$ and $n+m$, incorporating an anharmonic defect function $\delta(n,m)$ that captures deviations from harmonicity up to the second order.
2. Kinetic Modeling of Populations: The model assumes a Treanor distribution for the vibrational states to account for overpopulation when $T_v > T_g$. To address the physical divergence of the pure Treanor distribution at high levels, the formulation integrates a flux-corrected Gordiets ratio ($G(n,m)$). This ratio models the plateau region beyond the Treanor minimum ($n^*$), where Vibrational-Translational (V-T) relaxation overtakes Vibrational-Vibrational (V-V) up-pumping.
3. Unified Population Ratio: A continuous, piecewise-free population ratio is constructed by dynamically truncating the Treanor regime at the critical crossover point $n^*$. This involves an effective quantum jump ($m^*$) that transitions the system from the Treanor curve to the Gordiets plateau.
4. Derivation of the Correction Factor: By applying detailed balance to the anharmonic energy gaps and the unified population ratios, the author derives a closed-form, generalized anharmonic correction factor ($\Phi_{anh}^{(n,m)}$). This factor modifies the standard harmonic heating term to account for the kinetic competition between V-V up-pumping and V-T relaxation.
5. Macroscopic Decoupling: To enable efficient implementation in fluid solvers, a decoupled formulation is derived. This assumes the forward electron-impact rate coefficient is independent of the initial vibrational level, allowing the separation of collision cross-sections from the thermodynamic state of the vibrational manifold. This reduces the computational burden from a full state-to-state (StS) matrix to a set of channel-specific rates.

Key Contributions

• Unified Analytic Closure: The paper presents a single governing equation for electron-vibrational heat exchange that remains valid across both low-temperature (Boltzmann-like) and high-temperature (Treanor-Gordiets) regimes.
• Thermodynamic Consistency: The formulation strictly enforces the second law of thermodynamics. It is proven that when $T_e = T_v = T_g$, the net energy transfer vanishes, and the heating rate exactly equals the cooling rate ($Q_{v-e} = Q_{e-v}$), a condition often violated by previous models.
• Anharmonic Correction Factor: The introduction of $\Phi_{anh}$ explicitly captures the physics of anharmonic distortion, correcting the underestimation of heating in $T_v > T_g$ regimes and the overestimation in $T_g > T_v$ regimes.
• Computational Efficiency: The decoupled macroscopic model achieves a 40- to 70-fold reduction in rate evaluations compared to full StS benchmarks while maintaining high fidelity (deviating by at most 15% from exact StS solutions).

Results
The model was validated against direct state-to-state (StS) summation benchmarks for nitrogen ($N_2$) and carbon monoxide ($CO$) plasmas across two distinct regimes:

1. Cold-Gas Regime ($T_g = 300$ K, $T_e = 2.0$ eV): In this regime, characteristic of PAC inter-pulse stages, $T_v$ can significantly exceed $T_g$. The harmonic model underestimated the heating rate by a factor of five. The proposed anharmonic model, with $\Phi_{anh} > 1$, accurately recovered the StS benchmark by accounting for the superelastic surplus from overpopulated high-lying states.
2. Hot-Gas Regime ($T_g = 20,000$ K, $T_e = 0.4$ eV): In this regime, typical of hypersonic shock heating, $T_g > T_v$. The harmonic model overpredicted heating by approximately 15%. The generalized model correctly scaled down the heating flux ($\Phi_{anh} < 1$) to match the StS benchmark.
3. Decoupled Performance: The decoupled macroscopic model, utilizing dynamically weighted representative levels, closely matched the exact summation results, demonstrating that the anharmonic physics can be captured without resolving every individual state-to-state transition.

Significance
The paper claims that this unified framework provides a governing equation for heat exchange between electrons and excited states in non-equilibrium environments, including plasma-assisted combustion and hypersonic flows. By resolving the limitations of conventional harmonic models and ensuring thermodynamic consistency, the work enables the development of accurate, rate-limited reduced-order models for macroscopic fluid solvers. This is particularly significant for applications where full state-to-state kinetic data is incomplete or computationally prohibitive, offering a robust path toward simulating complex plasma-gas interactions with high fidelity and reduced computational cost.