Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

Modified Marrone-Treanor dissociation model: formulation and benchmarking for diatom/atom mixtures

T0 review · 2 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read A two-temperature dissociation model with QCT-derived rates matches direct molecular simulation for nitrogen and oxygen.

desk verdict Solid engineering-model paper; the QCT-derived parameterization is genuinely useful, but the concentration-dependent correction factor is a stated modeling choice that the benchmarks never actually exercise. read the letter →

arxiv 2506.06236 v2 pith:FPEUJ3A2 submitted 2025-06-06 physics.chem-ph

classification physics.chem-ph
keywords two-temperaturemodelMarrone-Treanordissociationkineticsquasiclassicaltrajectoriesnon-Boltzmanndepletiondirectmolecularsimulationhypersonicnonequilibriumairvibrationalrelaxation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that a modest upgrade of the 1963 Marrone-Treanor two-temperature dissociation model, with rate and energy parameters fit to quasiclassical trajectory calculations on ab initio potential energy surfaces, is enough to capture the physics of shock-heated dissociating air. The key claim is that, once the Boltzmann-sampled trajectory data are multiplied by two fixed non-Boltzmann correction factors ($f_k^{\mathrm{NB}} = 0.5$ for the dissociation rate and $f_\epsilon^{\mathrm{NB}} = 0.85$ for vibrational energy removed per dissociation), the model reproduces the main features of direct molecular simulations of N$_2$/N and O$_2$/O mixtures in both isothermal and adiabatic heat baths. The authors also propose a concentration-dependent interpolation of these factors so that the corrections vanish as the mixture approaches chemical equilibrium, which keeps recombination rates consistent with detailed balance. If correct, the work provides a self-consistent, computationally inexpensive two-temperature closure for CFD that ties dissociation rates, vibrational relaxation, and energy removal to the same ab initio surfaces.

What carries the argument

The load-bearing object is the modified Marrone-Treanor two-temperature expression for the dissociation rate coefficient and its companion formula for vibrational energy removal. The rate uses the original Marrone-Treanor nonequilibrium factor $Z(T,T_v)$ built from vibrational partition functions evaluated at a pseudotemperature $T_p$ that depends on $T$, $T_v$, and a characteristic probability temperature $U(T)$; the energy removal uses Knab's simple-harmonic-oscillator formula evaluated at the same pseudotemperature. The paper's main addition is fitting $U$ as a linear function of inverse temperature ($1/U = U_A/T + 1/U^*$, Eq. 15) from QCT data, so the same parameters describe both the rate and the energy removed, and multiplying both by the concentration-dependent non-Boltzmann factors $f_k^{\rm VNB}(\zeta)$ and $f_\epsilon^{\rm VNB}(\zeta)$.

What would settle it

Run DMS isothermal or adiabatic heat-bath calculations for a mixture with $\zeta>1$ (net recombination, e.g. cooling N/N$_2$ or O/O$_2$ starting overpopulated in atoms) and compare the late-time dissociation rate and vibrational energy removal against MMT predictions; if the inferred correction factors depart from the 0.5/0.85 pair by more than the 2-5x depletion band, the universality assumption fails. A simpler check is isothermal DMS for N$_2$+O$_2$ mixtures, whose relaxation and dissociation benchmark is not among the four calibrated reactions.

Watch

Extended reading notes

Core claim

The paper's central discovery is that the classical Marrone-Treanor preferential-dissociation model, equipped with modified parameters extracted from QCT data, forms a complete two-temperature air-dissociation model whose predictions track direct molecular simulation. Concretely, the dissociation rate is written as $k_{\rm diss}^{\rm MMT-NB}(T,T_v) = k_{\rm Arr}(T)\, Z(T,T_v)\, f_k^{\rm NB}$ (Eq. 10), and the average vibrational energy removed per dissociation is $\langle \epsilon_{v,A_2}\rangle_{\rm diss}^{\rm MMT-NB} = \langle \epsilon_{v,A_2}\rangle_{\rm diss}^{\rm Knab}\, f_\epsilon^{\rm NB}$ (Eq. 17), with all Arrhenius parameters $A$, $n$, $\theta_D$, and the two Marrone-Treanor pseudotemperature parameters $U$, $U^*$ fit to QCT data from ab initio PESs. The non-Boltzmann factors $f_k^{\rm NB}=0.5$ and $f_\epsilon^{\rm NB}=0.85$ are calibrated from quasi-steady-state DMS data and are made concentration-dependent through the ratio $\zeta$ of Eq. (30), so that diffusionless heat-bath simulations approach the correct equilibrium dissociation and recombination rates. Benchmarks show close agreement with DMS for N$_2$/N and O$_2$/O in isothermal baths at 5,000-20,000 K and adiabatic baths at two enthalpies each.

Load-bearing premise

The argument assumes that two fixed correction factors, $f_k^{\rm NB}=0.5$ and $f_\epsilon^{\rm NB}=0.85$, calibrated from quasi-steady-state dissociation of four reactions, remain valid for all collision partners, temperatures, and mixture compositions, and that the linear interpolation to unity at chemical equilibrium is a faithful description of how non-Boltzmann effects die out.

Editorial extensions

If this is right

  • A single two-temperature model now supplies all three coupled ingredients -- dissociation rate, vibrational energy removed per dissociation, and vibrational relaxation time -- from the same ab initio potential energy surfaces, removing the usual mix-and-match of empirical correlations.
  • In isothermal heat-bath benchmarks, the MMT model reproduces the DMS quasi-steady-state plateau in $T_v$ below the bath temperature, including the counterintuitive reversal at high bath temperatures where rotational relaxation delays vibrational excitation.
  • In adiabatic post-shock-like conditions, the model tracks DMS temperature and composition histories across four decades of time, with the largest early-time discrepancies narrowing to near-exact agreement after roughly 10 microseconds.
  • Because the formulation costs only a few partition-function evaluations per cell per time step, it is practical to deploy in large-scale CFD, unlike state-resolved master equations or DMS.
  • The concentration-dependent factors ensure forward and backward rates obey detailed balance at equilibrium, so the model does not freeze dissociation at a nonphysical depleted state.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial extension: the constant factors $f_k^{\rm NB}=0.5$ and $f_\epsilon^{\rm NB}=0.85$ are calibrated for four reactions in the QSS regime; a natural test is to extract them from DMS for N$_2$+O$_2$ and NO-forming reactions, or at temperatures outside the 4,000-30,000 K fit range, where the paper's own benchmarks show larger deviations.
  • Editorial extension: the capping of the variable factors at unity for recombination-dominated flows ($\zeta>1$) is an admitted modeling choice; if three-body QCT or DMS can quantify preferential recombination into near-threshold levels, the model could be extended to overpopulation factors greater than one, which would matter for nozzle expansions and boundary layers.
  • Editorial extension: the collapse of support-factor and energy-change data when normalized by dissociation energy suggests a possible universal scaling law for diatom dissociation; testing whether the same normalized curves hold for NO or other diatoms would separate species-specific rate parameters from genuinely universal physics.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper presents a modified Marrone-Treanor (MMT) two-temperature dissociation model for N2 and O2 with collision partners N2, O2, N, and O. All rate and vibrational-energy parameters are fitted to quasiclassical trajectory (QCT) data on Minnesota ab initio PESs, and vibrational relaxation times are fitted to DMS data on the same PES family. The model multiplies the conventional MMT rate and Knab vibrational-energy-removal expressions by non-Boltzmann correction factors f_k^NB=0.5 and f_eps^NB=0.85, which are calibrated to QSS DMS data and then made composition-dependent through a piecewise-linear function of the nonequilibrium concentration ratio zeta (Eq. 31). The manuscript benchmarks the CFD implementation against DMS heat-bath calculations for N2/N and O2/O mixtures in isothermal and adiabatic conditions and reports good agreement for dissociation rates, vibrational temperatures, and compositions.

Significance. If the result holds, the paper is significant for hypersonic CFD: it offers a self-consistent two-temperature model in which the dissociation rate, vibrational energy removal, and vibrational relaxation time all trace back to the same ab initio PESs, while remaining cheap enough for large-scale simulations. The QCT fitting is carefully done, with equilibrium Arrhenius fits at R^2>0.9998 and nonequilibrium predictions at 5-22% maximum deviation; the nonequilibrium QCT test set is a genuine predictive check because all parameters are fixed from equilibrium data. The heat-bath benchmarks credibly show that the constant-factor MMT model captures the QSS dissociation regime for N2/N and O2/O. The main weakness is that the novel composition-dependent part of the model, Eq. (31), is never exercised in the benchmarks, and the two cross-reactions N2+O2 and O2+N2 are assigned universal non-Boltzmann factors without QSS DMS validation.

major comments (2)
  1. [Sec. VII and Eq. (31)] The central verification claim for the concentration-dependent non-Boltzmann factors is not supported by the benchmarks. Section VII states explicitly that constant factors f_k^NB=0.5 and f_eps^NB=0.85 are employed throughout all MMT calculations in that section. The isothermal benchmarks therefore test only the QSS dissociation regime (zeta<1), and the adiabatic benchmarks, although they approach equilibrium, cannot validate the interpolation path of Eq. (31) because the same constant factor multiplies both forward and reverse rates in Eq. (28) and cancels at the equilibrium composition. Section VI itself labels Eq. (31) a modeling choice and states that no data exist to constrain the recombining regime. The authors should either rerun the adiabatic cases with the variable factors of Eq. (31) and compare the time-dependent approach to equilibrium, or explicitly restrict the verification claim to the QSS regime and describe Eq. (31) as an unvalidated modeling extension.
  2. [Sec. V, Table 2] The universality of the non-Boltzmann factors is established only for four of the six reactions. The constants f_k^NB=0.5 and f_eps^NB=0.85 are calibrated from DMS QSS data for N2+N2, N2+N, O2+O2, and O2+O (Figs. 11-12), and then applied without further validation to N2+O2 and O2+N2. Since the reaction-specific Arrhenius and U parameters vary markedly with collision partner in Table 2, the assumption that the depletion ratio is partner-independent is not self-evident. The benchmarks in Sec. VII do not exercise N2+O2 or O2+N2, so the six-reaction claim is partly extrapolated. The paper should either add QSS DMS comparisons for the two cross-reactions or explicitly flag those table rows as unverified for the non-Boltzmann factors.
minor comments (4)
  1. [Sec. VI, Fig. 13] The nonlinear alternative from Ref. [66] is shown but not compared quantitatively; a sentence stating whether the linear and nonlinear forms differ significantly over the plotted zeta range would help readers judge the sensitivity to this modeling choice.
  2. [Sec. II.B, Eq. (15)] The manuscript could state the valid temperature range for the linear U fit more prominently, since the fitted parameters in Table 2b imply that 1/U(T) can change sign at temperatures outside the QCT fitting range; this is unlikely to affect the present benchmarks but matters for future five-species-air use.
  3. [Sec. VII] The initial-condition mismatch between CFD and DMS (rotational mode pre-equilibrated in CFD but not in DMS) is discussed but not quantified; adding a short remark on the expected size of the resulting early-time bias in Tv would strengthen the comparison.
  4. [Fig. 6 caption] The caption abbreviates Marrone and Treanor as "M and T (1963)"; using the full author names would be clearer.

Circularity Check

2 steps flagged · score 6.0 of 10

QSS heat-bath agreement is largely forced by DMS-calibrated f factors and DMS-fitted relaxation times; independent QCT predictions and adiabatic checks keep circularity partial.

  1. fitted input called prediction [Sec. V (calibration) and Sec. VII (benchmarking)]
    "However, we choose to apply the same factors f_k^NB = 0.5 and f_eps^NB = 0.85 to all four reactions so as to not introduce unnecessary complexity to the MMT model. ... Further note that constant non-Boltzmann factors of f_k^NB = 0.5 and f_eps^NB = 0.85 are employed in Eqs. (10) and (17) respectively throughout all MMT calculations presented in this section."

    The Sec. V factors are calibrated so that MMT's QSS dissociation rate and vibrational-energy removal match DMS data in Figs. 11-12. The Sec. VII isothermal heat baths use those same constants in Eqs. (10) and (17), so during the QSS plateau the MMT mole-fraction slope is k_Arr*Z*0.5 and the energy-removal rate is Knab*0.85, precisely the fitted quantities. The benchmark therefore re-plots the calibration rather than testing a prediction; the Tv plateau is also controlled by the same fitted f_eps and by DMS-fitted relaxation times. The adiabatic cases inherit the same constants in their early QSS phase, so they only partially escape this circularity.

  2. fitted input called prediction [Sec. II.C (relaxation-time fits) and Sec. VII.A (benchmark comparison)]
    "We have derived all these parameters from curve fits to the DMS results of Ref. [37] which relied on the most recent ab initio PESs released by the Minnesota chemists. ... confirming that the values for calculating tau_v^{O2-O2} and tau_v^{O2-O} from Table 3 do an excellent job of reproducing the DMS behavior."

    The pair-wise vibrational relaxation times in Table 3 are obtained by curve-fitting DMS relaxation data from Ref. [37]. The Sec. VII heat-bath benchmarks then use these same tau values inside the Landau-Teller source term, Eq. (6), and compare the resulting Tv evolution against DMS. Since the early Tv rise is largely determined by these fitted relaxation times, the vibrational-excitation agreement in the benchmarks is a reproduction of the fit, not an independent validation. This is secondary to the rate-factor calibration but reinforces that the 'verification' partly reuses its inputs.

full rationale

Two fitted inputs are reused as validation targets. First, the non-Boltzmann factors f_k^NB=0.5 and f_eps^NB=0.85 are calibrated to DMS quasi-steady-state rate coefficients and energy removal (Sec. V, Figs. 11-12), and the same constants are then used in all Sec. VII benchmarks. During the QSS phase of the isothermal heat baths, the dissociation rate and vibrational-energy removal rate are exactly those calibrated quantities, so the close agreement in mole-fraction slope and Tv plateau is a consistency check rather than an independent prediction. Second, the vibrational relaxation times (Table 3) are curve-fits to DMS data from Ref. [37]; the heat-bath Tv evolution is largely governed by these fitted values, so vibrational-excitation agreement is also partly forced. The adiabatic cases provide partial independent evidence because T and Tv evolve together and the late-time approach is not directly calibrated; however, they inherit the same fitted factors in the early QSS phase. The genuinely independent part of the derivation is the QCT-to-MMT construction: Arrhenius parameters and U are fitted to equilibrium QCT data and then successfully reproduce the nonequilibrium QCT test set (Sec. IV.E), including the 22%/0.24 eV error bounds; that step is predictive and not circular. Sec. VI explicitly calls Eq. (31) a modeling choice and says no data constrain the recombination regime, and Sec. VII states that constant factors, not Eq. (31), are used; therefore the advertised concentration-dependent interpolation is never exercised and the 'verification' does not cover it. This is a limitation of the validation, not an additional circular step.

Assumptions & free parameters 18 free parameters · 8 assumptions · 0 invented entities

The model rests on a large set of fitted coefficients: six Arrhenius fits, six MMT U-parameter pairs, two global non-Boltzmann constants, and four relaxation-time fits with four coefficients each, all fit to QCT or DMS data. There are no invented physical entities. The core axioms are standard two-temperature CFD assumptions plus two explicit ad hoc choices: universal non-Boltzmann factors and the Eq. (31) recombination interpolation.

free parameters (18)
  • Arrhenius C,n,T_D for N2+N2 (variable T_D fit, Table 2a) = C=3.5967e18 cm3 mol-1 s-1 K^-n, n=-0.7017, T_D=117529 K
    Fit to QCT equilibrium dissociation rates for N2+N2 via Eq. (11).
  • Arrhenius C,n,T_D for N2+N = C=7.9920e17 cm3 mol-1 s-1 K^-n, n=-0.5625, T_D=113957 K
    Fit to QCT equilibrium dissociation rates for N2+N via Eq. (11).
  • Arrhenius C,n,T_D for N2+O2 = C=5.0420e19 cm3 mol-1 s-1 K^-n, n=-0.9991, T_D=116892 K
    Fit to QCT equilibrium dissociation rates for N2+O2 via Eq. (11).
  • Arrhenius C,n,T_D for O2+O2 = C=3.6932e18 cm3 mol-1 s-1 K^-n, n=-0.7695, T_D=60540 K
    Fit to QCT equilibrium dissociation rates for O2+O2 via Eq. (11).
  • Arrhenius C,n,T_D for O2+O = C=9.2109e17 cm3 mol-1 s-1 K^-n, n=-0.6541, T_D=60552 K
    Fit to QCT equilibrium dissociation rates for O2+O via Eq. (11).
  • Arrhenius C,n,T_D for O2+N2 = C=3.1942e17 cm3 mol-1 s-1 K^-n, n=-0.5430, T_D=62949 K
    Fit to QCT equilibrium dissociation rates for O2+N2 via Eq. (11).
  • MMT U parameters for N2+N2 = a_U=0.3868, U*=254556 K
    Fit to QCT equilibrium vibrational energy removed per dissociation via Eq. (16) with Eq. (15); Table 2a.
  • MMT U parameters for N2+N = a_U=0.3668, U*=478708 K
    Fit to QCT equilibrium vibrational energy removed per dissociation; Table 2a.
  • MMT U parameters for N2+O2 = a_U=0.3001, U*=210253 K
    Fit to QCT equilibrium vibrational energy removed per dissociation; Table 2a.
  • MMT U parameters for O2+O2 = a_U=0.3965, U*=57343 K
    Fit to QCT equilibrium vibrational energy removed per dissociation; Table 2a.
  • MMT U parameters for O2+O = a_U=0.3537, U*=237290 K
    Fit to QCT equilibrium vibrational energy removed per dissociation; Table 2a.
  • MMT U parameters for O2+N2 = a_U=0.3620, U*=385466 K
    Fit to QCT equilibrium vibrational energy removed per dissociation; Table 2a.
  • Non-Boltzmann rate factor f_k^NB = 0.5 (all six reactions)
    Calibrated to DMS QSS dissociation rates in Sec. V; applied uniformly.
  • Non-Boltzmann energy factor f_eps^NB = 0.85 (all six reactions)
    Calibrated to DMS QSS vibrational energy removed per dissociation in Sec. V; applied uniformly.
  • Vibrational relaxation fit parameters for N2-N2 = A_low=221.0 K^(1/3), B_low=-24.83, A_high=33.30, B_high=-17.31
    Fit to DMS relaxation times from Ref. [37]; low-T segment matches Millikan-White; Table 3.
  • Vibrational relaxation fit parameters for N2-N = A_low=239.4 K^(1/3), B_low=-28.52, A_high=32.34, B_high=-18.39
    Fit to DMS relaxation times from Ref. [37]; Table 3.
  • Vibrational relaxation fit parameters for O2-O2 = A_low=129.0 K^(1/3), B_low=-22.29, A_high=-4.384, B_high=-16.70
    Fit to DMS relaxation times from Ref. [37]; Table 3.
  • Vibrational relaxation fit parameters for O2-O = A_low=3.018 K^(1/3), B_low=-17.74, A_high=-121.6, B_high=-15.82
    Fit to DMS relaxation times from Ref. [37]; Table 3.
assumptions (8)
  • domain assumption A single mixture vibrational temperature describes all vibrational modes, with Landau-Teller relaxation.
    Sec. II.A, Eqs. (5)-(6); this is the standard two-temperature CFD assumption the entire model is built around.
  • domain assumption Quasiclassical trajectories on the University of Minnesota ab initio PESs accurately represent dissociation dynamics.
    All QCT data are from these PESs (Refs. [13]-[20]); the model inherits their accuracy and electronic-ground-state scope.
  • ad hoc to paper Non-Boltzmann correction factors f_k=0.5 and f_eps=0.85 are universal constants across all six reactions and all temperatures.
    Sec. V: the authors note the best factors vary somewhat with reaction and temperature but choose single values to keep the model simple. This is load-bearing for the benchmark agreement.
  • ad hoc to paper The piecewise-linear interpolation of Eq. (31), capped at unity, correctly describes non-Boltzmann factors between fully recombined and equilibrium mixtures.
    Sec. VI: the authors explicitly call Eq. (31) a modeling choice and state they lack first-principles data for recombination-dominated regimes (zeta > 1).
  • domain assumption Detailed balance is used to set recombination rates even in nonequilibrium CFD states.
    Sec. VI, Eqs. (26)-(28); recombination rate is tied to the thermal Arrhenius rate and equilibrium constant, an approximation common to two-temperature CFD models.
  • domain assumption Knab's simple harmonic oscillator approximation accurately represents the vibrational partition function and energy removed per dissociation.
    Sec. IV.C, Eqs. (13) and (16); the authors show the approximation introduces errors up to about 10% of D0, which they accept for CFD efficiency.
  • domain assumption Electronic excited states can be neglected in thermodynamics and kinetics for these benchmark mixtures.
    Sec. VII: thermodynamic properties and DMS reference calculations exclude electronic excited-state contributions.
  • domain assumption Rotational modes are fully equilibrated with translational modes in the CFD model from t=0.
    Sec. VII: CFD-MMT initializes T_rot=T, whereas DMS resolves rotational relaxation, causing the acknowledged early-time lag in vibrational temperature and composition profiles.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Modified Marrone-Treanor dissociation model: formulation and benchmarking for diatom/atom mixtures." pith.science (2026). https://pith.science/paper/FPEUJ3A2

@misc{pith2026250606236,
  author       = {Pith},
  title        = {Pith review of: Modified Marrone-Treanor dissociation model: formulation and benchmarking for diatom/atom mixtures},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FPEUJ3A2}},
  note         = {Machine review of arXiv:2506.06236}
}
read the original abstract

We present a modified Marrone-Treanor model for dissociation with rate parameters derived exclusively from quasiclassical trajectory calculations on ab initio potential energy surfaces. Analysis of the trajectory dataset for reactant O2 and N2 diatoms sampled from Boltzmann internal energy distributions over a wide T,Tv range indicates that a modified version of the classical Marrone-Treanor two-temperature model captures the most relevant physics of shock-heated dissociating diatomic species very well. We find that simple correction factors account for non-Boltzmann depletion effects observed in direct molecular simulations employing the same potentials. The concentration-dependent functional form proposed for these correction factors ensures that depletion effects vanish at chemical equilibrium. Based on comparisons in isothermal and adiabatic heat baths we verify that the resulting two-temperature dissociation model accurately reproduces all major features observed in the direct molecular simulations, while remaining computationally inexpensive enough for large-scale computational fluid dynamics simulations.

Figures

Figures reproduced from arXiv: 2506.06236 by the authors.

Figure 3
Figure 3. of Ref. [25] shows a contour plot of dissociation [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figure 4
Figure 4. Although neither of these molecules is quasibound, t [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. We find that the remainder energy is a relatively good p [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Modified Marrone-Treanor model: parameterization and benchmarking for five-species air

    physics.chem-ph 2025-06 conditional novelty 6.0 of 10

    A Modified Marrone-Treanor model for five-species air is parameterized from ab initio trajectories and reproduces direct molecular simulation heat-bath benchmarks while predicting more NO and slower N2 dissociation th...

Reference graph

Works this paper leans on

66 extracted references · 26 canonical work pages · cited by 1 Pith paper

  1. [1]

    Planetary-Entry Gas Dynamics,

    Gnoffo, P. A., “Planetary-Entry Gas Dynamics,” Annual Review of Fluid Mechanics , Vol. 31, No. 1, 1999, pp. 459–494. doi:10.1146/annurev.fluid.31.1.459

  2. [2]

    Rate Effects in Hypersonic Flows,

    Candler, G. V ., “Rate Effects in Hypersonic Flows,” Annual Review of Fluid Mechanics , Vol. 51, No. 1, 2019, pp. 379 – 402. doi:10.1146/annurev-fluid-010518-040258

  3. [3]

    Theory of Ra diation from Luminous Shock Waves in Nitrogen,

    Hammerling, P., Teare, J. D., and Kivel, B., “Theory of Ra diation from Luminous Shock Waves in Nitrogen,” The Physics of Fluids, Vol. 2, No. 4, 1959, pp. 422–426. doi:10.1063/1.1724413

  4. [4]

    Effect of Dissociatio n on the Rate of Vibrational Relaxation,

    Treanor, C. E., and Marrone, P. V ., “Effect of Dissociatio n on the Rate of Vibrational Relaxation,” The Physics of Fluids , Vol. 5, No. 9, 1962, pp. 1022–1026. doi:10.1063/1.1724467

  5. [5]

    Chemical Relaxation with Preferential Dissociation from Excited Vibrational Levels,

    Marrone, P. V ., and Treanor, C. E., “Chemical Relaxation with Preferential Dissociation from Excited Vibrational Levels,” The Physics of Fluids, Vol. 6, No. 9, 1963, pp. 1215–1221. doi:10.1063/1.1706888

  6. [6]

    Conservatio n Equations and Physical Models for Hypersonic Air Flows in Thermal and Chemical Nonequilibrium,

    Gnoffo, P. A., Gupta, R. N., and Shinn, J. L., “Conservatio n Equations and Physical Models for Hypersonic Air Flows in Thermal and Chemical Nonequilibrium,” Tech. Rep. TP-2867, NASA, February 1989

  7. [7]

    A Study of the Excitat ion of Vibrations and Dissociation of Oxygen Molecules at Hi gh Temperatures,

    Losev, S. A., and Generalov, N. A., “A Study of the Excitat ion of Vibrations and Dissociation of Oxygen Molecules at Hi gh Temperatures,” Soviet Physics Doklady, Vol. 6, 1962, p. 1081

  8. [8]

    Assessment of a T wo-Temperature Kinetic Mode l for Dissociating and Weakly Ionizing Nitrogen,

    Park, C., “Assessment of a T wo-Temperature Kinetic Mode l for Dissociating and Weakly Ionizing Nitrogen,” Journal of Thermophysics and Heat Transfer, Vol. 2, No. 1, 1988, pp. 8–16. doi:10.2514/3.55

Show all 66 references
  1. [9]

    Theo ry and Validation of the Physically Consistent Coupled Vibration- Chemistry-Vibration Model,

    Knab, O., Fruehauf, H. H., and Messerschmid, E. W., “Theo ry and Validation of the Physically Consistent Coupled Vibration- Chemistry-Vibration Model,” Journal of Thermophysics and Heat Transfer , Vol. 9, No. 2, 1995, pp. 219–226. doi:10.2514/3. 649

  2. [10]

    Assessment of Classical Impulsive Models of Dissociation in Thermochemical Nonequilibrium,

    Luo, H., Alexeenko, A. A., and Macheret, S. O., “Assessment of Classical Impulsive Models of Dissociation in Thermochemical Nonequilibrium,” Journal of Thermophysics and Heat Transfer , Vol. 32, No. 4, 2018, pp. 861–868. doi:10.2514/1.T5375

  3. [11]

    Systematics of Vibra tional Relaxation,

    Millikan, R. C., and White, D. R., “Systematics of Vibra tional Relaxation,” The Journal of Chemical Physics , Vol. 39, No. 12, 1963, pp. 3209–3213. doi:10.1063/1.1734182

  4. [12]

    Reactive scattering c ross sections III: Quasiclassical and semiclassical metho ds,

    Truhlar, D., and Muckerman, J., “Reactive scattering c ross sections III: Quasiclassical and semiclassical metho ds,” Atom- Molecule Collision Theory , Springer, 1979, pp. 505–566

  5. [13]

    G lobal ab initio Ground-State Potential Energy Surface of N 4,

    Paukku, Y ., Yang, K. R., Varga, Z., and Truhlar, D. G., “G lobal ab initio Ground-State Potential Energy Surface of N 4,” The Journal of Chemical Physics , Vol. 139, No. 4, 2013, p. 044309. doi:10.1063/1.4811653

  6. [14]

    E rratum: “Global ab initio Ground-State Potential Energy Surface of N 4

    Paukku, Y ., Yang, K. R., Varga, Z., and Truhlar, D. G., “E rratum: “Global ab initio Ground-State Potential Energy Surface of N 4” [J. Chem. Phys. 139, 044309 (2013)],” The Journal of Chemical Physics , Vol. 140, No. 1, 2014, p. 019903. doi: 10.1063/1.4861562. 46

  7. [15]

    An Improved Potential Energy Surface and Multi-Temperature Q uasiclassical Trajectory Calculations of N 2 + N 2 Dissociation Reactions,

    Bender, J. D., Valentini, P., Nompelis, I., Paukku, Y ., Varga, Z., Truhlar, D. G., Schwartzentruber, T., and Candler, G. V ., “An Improved Potential Energy Surface and Multi-Temperature Q uasiclassical Trajectory Calculations of N 2 + N 2 Dissociation Reactions,” The Journal ...

  8. [16]

    Potential energy surface f or high-energy N + N2 collisions,

    Varga, Z., and Truhlar, D., “Potential energy surface f or high-energy N + N2 collisions,” Phys. Chem. Chem. Phys. , Vol. 23, 2021, pp. 26273–26284. doi:10.1039/D1CP04373K

  9. [17]

    Potential Energy Surfaces of Quintet and Singlet O4,

    Paukku, Y ., Yang, K. R., Varga, Z., Song, G., Bender, J. D ., and Truhlar, D. G., “Potential Energy Surfaces of Quintet and Singlet O4,” The Journal of Chemical Physics , Vol. 147, No. 3, 2017, p. 034301. doi:10.1063/1.4993624

  10. [18]

    Potential En ergy Surface of Triplet O 4,

    Paukku, Y ., Varga, Z., and Truhlar, D. G., “Potential En ergy Surface of Triplet O 4,” The Journal of Chemical Physics , Vol. 148, No. 12, 2018, p. 124314. doi:10.1063/1.5017489

  11. [19]

    Potential Energy Surfaces for O + O2 Collisions,

    Varga, Z., Paukku, Y ., and Truhlar, D. G., “Potential Energy Surfaces for O + O2 Collisions,” The Journal of Chemical Physics, Vol. 147, No. 15, 2017, p. 154312. doi:10.1063/1.4997169

  12. [20]

    Potential Energy Surface of Triplet N 2O2,

    Varga, Z., Meana-Pañeda, R., Song, G., Paukku, Y ., and T ruhlar, D. G., “Potential Energy Surface of Triplet N 2O2,” The Journal of Chemical Physics , Vol. 144, No. 2, 2016, p. 024310. doi:10.1063/1.4939008

  13. [21]

    Dynamics of Nitrogen Dissociation from Direct Molecular Simulation,

    Valentini, P., Schwartzentruber, T. E., Bender, J. D., and Candler, G. V ., “Dynamics of Nitrogen Dissociation from Direct Molecular Simulation,” Physical Review Fluids, Vol. 1, No. 4, 2016, p. 043402. doi:10.1103/PhysRevFluids .1.043402

  14. [22]

    Quasiclassical Trajectory Analysis o f Nitrogen for High-Temperature Chemical Kinetics,

    Chaudhry, R. S., Bender, J. D., Schwartzentruber, T. E. , and Candler, G. V ., “Quasiclassical Trajectory Analysis o f Nitrogen for High-Temperature Chemical Kinetics,” Journal of Thermophysics and Heat Transfer , Vol. 32, No. 4, 2018, pp. 833–845. doi:10.2514/1.T5484

  15. [23]

    Analysis of Dissociation and Internal Energy Transfer in High-Energy N 2+O2 Collisions using the Quasiclassical Trajectory Method,

    Chaudhry, R. S., Bender, J. D., Valentini, P., Schwartzentruber, T. E., and Candler, G. V ., “Analysis of Dissociation and Internal Energy Transfer in High-Energy N 2+O2 Collisions using the Quasiclassical Trajectory Method,” 46th AIAA Thermophysics Conference, AIAA Paper 2016...

  16. [24]

    Quasiclassical Traject ory Analysis of Oxygen Dissociation via O 2, O, and N 2,

    Chaudhry, R. S., Grover, M. S., Bender, J. D., Schwartze ntruber, T. E., and Candler, G. V ., “Quasiclassical Traject ory Analysis of Oxygen Dissociation via O 2, O, and N 2,” 2018 AIAA Aerospace Sciences Meeting, AIAA Paper 2018-0237, 2018. doi:10.2514/6.2018-0237

  17. [25]

    Modeling and Analysis of Chemical Kin etics for Hypersonic Flows in Air,

    Chaudhry, R. S., “Modeling and Analysis of Chemical Kin etics for Hypersonic Flows in Air,” Ph.D. thesis, Universit y of Minnesota, Minneapolis, MN, November 2018. URL https://hdl.handle.net/11299/201709

  18. [26]

    Rovibrational Internal Energy Transfer and Dissociatio n of /u1D4412 (1Σ+ /u1D454) – /u1D441(4/u1D446/u1D462) System in Hypersonic Flows,

    Panesi, M., Jaffe, R. L., Schwenke, D. W., and Magin, T. E. , “Rovibrational Internal Energy Transfer and Dissociatio n of /u1D4412 (1Σ+ /u1D454) – /u1D441(4/u1D446/u1D462) System in Hypersonic Flows,” The Journal of Chemical Physics , Vol. 138, No. 4, 2013, p. 044312. doi:10.1...

  19. [27]

    State-Resolved Master Equa tion Analysis of Thermochemical Nonequilibrium of Nitroge n,

    Kim, J. G., and Boyd, I. D., “State-Resolved Master Equa tion Analysis of Thermochemical Nonequilibrium of Nitroge n,” Chemical Physics, Vol. 415, 2013, pp. 237–246. doi:10.1016/j.chemphys.201 3.01.027. 47

  20. [28]

    First Princi ples Calculation of Heavy Particle Rate Coefficients,

    Richard L. Jaffe, M. P., David W. Schwenke, “First Princi ples Calculation of Heavy Particle Rate Coefficients,” Hypersonic Nonequilibrium Flows: Fundamentals and Recent Advances , Vol. 247, edited by E. Josyula, Progress in Aerospace Scien ces, AIAA, New Y ork, 2015, pp. 103–158...

  21. [29]

    Direct Molecular Simulation of Nonequilibrium Dilute Gases,

    Schwartzentruber, T. E., Grover, M. S., and Valentini, P., “Direct Molecular Simulation of Nonequilibrium Dilute Gases,” Journal of Thermophysics and Heat Transfer , Vol. 32, No. 4, 2017, pp. 892–903. doi:10.2514/1.T5188

  22. [30]

    Direct Molecular Simulat ion of Nitrogen Dissociation Based on an ab initio Potential Energy Surface,

    Valentini, P., Schwartzentruber, T. E., Bender, J. D., Nompelis, I., and Candler, G. V ., “Direct Molecular Simulat ion of Nitrogen Dissociation Based on an ab initio Potential Energy Surface,” Physics of Fluids , Vol. 27, No. 8, 2015, p. 086102. doi:10.1063/1.4929394

  23. [31]

    Vibrational Energy Transfer and Collision-Induced Dissociation in O + O2 Collisions,

    Grover, M., Schwartzentruber, T., Varga, Z., and Truhlar, D., “Vibrational Energy Transfer and Collision-Induced Dissociation in O + O2 Collisions,” Journal of Thermophysics and Heat Transfer, Vol. 33, No. 3, 2019, pp. 797–807. doi:10.2514/1.T5551

  24. [32]

    Dir ect molecular simulation of internal energy relaxation and dissociation in oxygen,

    Grover, M., Torres, E., and Schwartzentruber, T., “Dir ect molecular simulation of internal energy relaxation and dissociation in oxygen,” Physics of Fluids, Vol. 31, 2019, p. 076107. doi:10.1063/1.5108666

  25. [33]

    Direct molecul ar simulation of nitrogen dissociation under adiabatic pos t-shock condi- tions,

    Torres, E., and Schwartzentruber, T., “Direct molecul ar simulation of nitrogen dissociation under adiabatic pos t-shock condi- tions,” Journal of Thermophysics and Heat transfer , Vol. 34, No. 4, 2020, pp. 801–815. doi:10.2514/1.T5970

  26. [34]

    High-Temperature Nonequilibrium Air Chemistry from F irst Principles,

    Torres, E., Geistfeld, E. C., and Schwartzentruber, T. E., “High-Temperature Nonequilibrium Air Chemistry from F irst Principles,” Journal of Thermophysics and Heat Transfer , Vol. 38, No. 2, 2024, pp. 260–291. doi:10.2514/1.T6863

  27. [35]

    Consistent kine tic–continuum dissociation model I. Kinetic formulation,

    Singh, N., and Schwartzentruber, T., “Consistent kine tic–continuum dissociation model I. Kinetic formulation, ” The Journal of Chemical Physics, Vol. 152, No. 22, 2020, p. 224302. doi:10.1063/1.5142752

  28. [36]

    Consistent kine tic-continuum dissociation model. II. Continuum formulat ion and verification,

    Singh, N., and Schwartzentruber, T., “Consistent kine tic-continuum dissociation model. II. Continuum formulat ion and verification,” The Journal of Chemical Physics , Vol. 152, No. 22, 2020, p. 224303. doi:10.1063/1.5142754

  29. [37]

    Characteris tic vibrational and rotational relaxation times for air spe cies from first- principles calculations,

    Torres, E., and Schwartzentruber, T. E., “Characteris tic vibrational and rotational relaxation times for air spe cies from first- principles calculations,” Journal of Thermophysics and Heat Transfer, Vol. 39, No. 2, 2025, pp. 223–249. doi:10.2514/1.T7042

  30. [38]

    POTLIB: An Online Library of Potential Energ y Surfaces,

    Varga, Z., Meana-Pañeda, R., Song, G., Paukku, Y ., and T ruhlar, D. G., “POTLIB: An Online Library of Potential Energ y Surfaces,” , 2021. URL comp.chem.umn.edu/potlib/

  31. [39]

    Multiscale Computational Analysis of N itrogen and Oxygen Gas-Phase Thermochemistry in Hypersonic Flows,

    Bender, J. D., “Multiscale Computational Analysis of N itrogen and Oxygen Gas-Phase Thermochemistry in Hypersonic Flows,” Ph.D. thesis, University of Minnesota, Minneapolis, MN, February 2016. URL https://hdl.handle.net/11299/178972

  32. [40]

    Investigation of Oxygen Dissociation and Vibrational Relaxation at Tempe ratures 4000–10800 K,

    Ibraguimova, L. B., Sergievskaya, A. L., Levashov, V . Y ., Shatalov, O. P., Tunik, Y . V ., and Zabelinskii, I. E., “Investigation of Oxygen Dissociation and Vibrational Relaxation at Tempe ratures 4000–10800 K,” The Journal of Chemical Physics , Vol. 139, No. 3, 2013, p. 034...

  33. [41]

    Vibrati onal relaxation time measurements in shock-heated oxygen and air from 2000 K to 9000 K using ultraviolet laser absorption,

    Streicher, J. W., Krish, A., and Hanson, R. K., “Vibrati onal relaxation time measurements in shock-heated oxygen and air from 2000 K to 9000 K using ultraviolet laser absorption,”Physics of Fluids, Vol. 32, No. 8, 2020, p. 086101. doi:10.1063/5.0015890. 48

  34. [42]

    Coupled vibr ation-dissociation time-histories and rate measurements in shock- heated, nondilute O 2 and O 2 − Ar mixtures from 6000 to 14000 K,

    Streicher, J., Krish, A., and Hanson, R., “Coupled vibr ation-dissociation time-histories and rate measurements in shock- heated, nondilute O 2 and O 2 − Ar mixtures from 6000 to 14000 K,” Physics of Fluids , Vol. 33, No. 5, 2021, p. 056107. doi:10.1063/5.0048059

  35. [44]

    High-te mperature vibrational relaxation and decomposition of sho ck-heated nitric oxide: II. Nitrogen dilution from 1900 to 8200 K,

    Streicher, J. W., Krish, A., and Hanson, R. K., “High-te mperature vibrational relaxation and decomposition of sho ck-heated nitric oxide: II. Nitrogen dilution from 1900 to 8200 K,” Physics of Fluids , Vol. 34, No. 11, 2022, p. 116123. doi: 10.1063/5.0122787

  36. [45]

    Review of Chemical-Kinetic Problems of Futu re NASA Missions, I: Earth Entries,

    Park, C., “Review of Chemical-Kinetic Problems of Futu re NASA Missions, I: Earth Entries,” Journal of Thermophysics and Heat Transfer, Vol. 7, No. 3, 1993, pp. 385–398. doi:10.2514/3.431

  37. [46]

    Nikitin, E., Theory of Elementary Atomic and Molecular Processes in Gase s, Clarendon Press, Oxford, 1974

  38. [47]

    Rate Constants for Chemical Reactions in H igh-Temperature Nonequilibrium Air,

    Jaffe, R. L., “Rate Constants for Chemical Reactions in H igh-Temperature Nonequilibrium Air,” Thermophysical Aspects of Re-Entry Flows, Vol. 103, edited by J. N. Moss and C. D. Scott, Progress in Aer ospace Sciences, AIAA, New Y ork, 1986, pp. 123–151. doi:10.2514/5.9781600865...

  39. [48]

    F., Hobson, M

    Riley, K. F., Hobson, M. P., and Bence, S. J., Mathematical Methods for Physics and Engineering , 3 rd ed., Cambridge University Press, Cambridge, UK, 2006

  40. [49]

    Measurement of the Rate of Dissociation o f Oxygen,

    Byron, S. R., “Measurement of the Rate of Dissociation o f Oxygen,” The Journal of Chemical Physics , Vol. 30, No. 6, 1959, pp. 1380–1392. doi:10.1063/1.1730209

  41. [50]

    Review of Finite-Rate Chemistry Models for A ir Dissociation and Ionization,

    Park, C., “Review of Finite-Rate Chemistry Models for A ir Dissociation and Ionization,” Molecular Physics and Hypersonic Flows, NATO ASI Series (Series C: Mathematical and Physical Scien ces), Vol. 482, edited by M. Capitelli, Kluwer Academic Publishers, Dordrecht, 1996, pp. ...

  42. [51]

    Chemistry Models for Air Dissociation ,

    Losev, S., Kovach, E., Makarov, V ., Pogosbekjan, M., an d Sergievskaya, A., “Chemistry Models for Air Dissociation ,” Molecular Physics and Hypersonic Flows, NATO ASI Series (Series C: Mathematical and Physical Sciences), Vol. 482, edited by M. Capitelli, Kluwer Academic Publ...

  43. [52]

    T wo-Temperature Interpretation of Dissociation Rate Data for N2 and O2,

    Park, C., “T wo-Temperature Interpretation of Dissociation Rate Data for N2 and O2,” 26th Aerospace Sciences Meeting, AIAA Paper 1988-458, 1988. doi:10.2514/6.1988-458

  44. [53]

    Park, C., Nonequilibrium Hypersonic Aerothermodynamics, John Wiley & Sons, Inc., New Y ork, 1990

  45. [54]

    State-to-State Master Equation and Direct Mole cular Simulation Study of Energy Transfer and Dissociation for th e N 2 − N System,

    Macdonald, R., Torres, E., Schwartzentruber, T., and P anesi, M., “State-to-State Master Equation and Direct Mole cular Simulation Study of Energy Transfer and Dissociation for th e N 2 − N System,” The Journal of Physical Chemistry A , 2020. doi:10.1021/acs.jpca.0c04029. 49

  46. [55]

    Non-Boltzmann v ibrational energy distributions and coupling to dissociat ion rate,

    Singh, N., and Schwartzentruber, T., “Non-Boltzmann v ibrational energy distributions and coupling to dissociat ion rate,” The Journal of Chemical Physics , Vol. 152, No. 22, 2020, p. 224301. doi:10.1063/1.5142732

  47. [56]

    Nonequilibri um Internal Energy Distributions during Dissociation,

    Singh, N., and Schwartzentruber, T. E., “Nonequilibri um Internal Energy Distributions during Dissociation,” Proceedings of the National Academy of Sciences , Vol. 115, No. 1, 2018, pp. 47–52. doi:10.1073/pnas.171384 0115

  48. [57]

    A simple model for state-specific diatomic dissociation,

    Gonzales, D. A., and Varghese, P. L., “A simple model for state-specific diatomic dissociation,” The Journal of Physical Chemistry, Vol. 97, No. 29, 1993, pp. 7612–7622. doi:10.1021/j100131 a034

  49. [58]

    Evaluation of sim ple rate expressions for vibration-dissociation coupling ,

    Gonzales, D. A., and Varghese, P. L., “Evaluation of sim ple rate expressions for vibration-dissociation coupling ,” Journal of Thermophysics and Heat Transfer, Vol. 8, No. 2, 1994, pp. 236–243. doi:10.2514/3.529

  50. [59]

    Vibrational rela xation models for dilute shock heated gases,

    Gonzales, D. A., and Varghese, P. L., “Vibrational rela xation models for dilute shock heated gases,” Chemical Physics, Vol. 195, No. 1, 1995, pp. 83–91. doi:10.1016/0301-0104(95)000 78-3

  51. [60]

    Direct molecul ar simulation of oxygen dissociation across normal shocks,

    Torres, E., and Schwartzentruber, T., “Direct molecul ar simulation of oxygen dissociation across normal shocks, ” Theoretical and Computational Fluid Dynamics , Vol. 36, No. 1, 2022, pp. 41–80. doi:10.1007/s00162-021-0 0596-6

  52. [61]

    Quasi-classical trajectory analysis of three-body coll ision induced recombination in neutral nitrogen and oxygen,

    Geistfeld, E. C., Torres, E., and Schwartzentruber, T. , “Quasi-classical trajectory analysis of three-body coll ision induced recombination in neutral nitrogen and oxygen,” The Journal of Chemical Physics , Vol. 159, No. 15, 2023, p. 154111. doi: 10.1063/5.0163942

  53. [62]

    Coefficients for C alculating Thermodynamic and Transport Properties of Indi vidual Species,

    McBride, B., Gordon, S., and Reno, M., “Coefficients for C alculating Thermodynamic and Transport Properties of Indi vidual Species,” Tech. rep., NASA Technical Memorandum 4513, 1993

  54. [63]

    Nonequilibri um Dissociation and Recombination Models for Hypersonic Fl ows,

    Singh, N., and Schwartzentruber, T. E., “Nonequilibri um Dissociation and Recombination Models for Hypersonic Fl ows,” AIAA Journal, Vol. 60, No. 5, 2022, pp. 2810–2825. doi:10.2514/1.J06115 4

  55. [64]

    Objective molecular dynamics investigation of d issociation and recombination kinetics in high-temperature nitrogen,

    Pahlani, G., Torres, E., Schwartzentruber, T., and Jam es, R. D., “Objective molecular dynamics investigation of d issociation and recombination kinetics in high-temperature nitrogen, ” Physics of Fluids , Vol. 35, No. 6, 2023, p. 067111. doi:10.1063/5. 0150492

  56. [65]

    State-to-state study of non-equili brium recombination of oxygen and nitrogen molecules,

    Macdonald, R. L., “State-to-state study of non-equili brium recombination of oxygen and nitrogen molecules,” The Journal of Chemical Physics, Vol. 160, No. 13, 2024, p. 134307. doi:10.1063/5.0195238

  57. [66]

    Vehicle-Scal e Simulations of Hypersonic Flows using the MMT Chemical Kin etics Model,

    Chaudhry, R., Boyd, I., and Candler, G., “Vehicle-Scal e Simulations of Hypersonic Flows using the MMT Chemical Kin etics Model,” AIAA Aviation 2020 Forum, 2020. doi:10.2514/6.2020-3272, AIAA 2020-3272

  58. [67]

    Direct molecul ar simulation of dissociating oxygen under adiabatic and no rmal shock wave conditions,

    Torres, E., and Schwartzentruber, T., “Direct molecul ar simulation of dissociating oxygen under adiabatic and no rmal shock wave conditions,” AIAA Scitech 2021 Forum, 2021. doi:10.2514/6.2021-0318, AIAA 2021-0318. 50

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.