{"id":"5d6e1ebd-226f-4ca3-a0e3-2634bba34c4d","arxiv_id":"2506.06236","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":18,"one_line_summary":"A two-temperature dissociation model, parameterized by ab initio trajectory calculations and non-Boltzmann correction factors, matches direct molecular simulations of N2 and O2 dissociation in heat baths.","lead":"The authors present an updated computer model for how oxygen and nitrogen molecules break apart in extremely hot, fast-moving air. The model's rates come from detailed simulations of individual collisions on accurate energy surfaces, plus two simple corrections for the depletion of highly excited molecules. If correct, it gives hypersonic vehicle designers a cheap way to compute air chemistry that matches much more expensive molecular simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The benchmarks in Sec. VII are run with constant f^NB factors, so the concentration-dependent Eq. (31) advertised in the central claim is never actually tested; the QSS dissociation limit is validated, but the equilibrium/recombination behavior is not.","rationale":"I read the paper as a careful engineering-model contribution. The QCT fitting is thorough: the Arrhenius fits have R^2 > 0.9998, the nonequilibrium factor fits are within 5-22%, and the parameter extraction from equilibrium data only makes the nonequilibrium test set genuinely predictive. The vibrational relaxation fits are openly tied to DMS data, which limits independence but is not a logical flaw. The most load-bearing soft spot is not the universality of the constant factors per se, but that the paper's own benchmarks do not exercise the concentration-dependent Eq. (31) that the abstract advertises. Sec. VII states in plain terms that constant factors are used throughout, so the comparisons validate the QSS dissociation regime only. The reader's weakest_assumption identified the Eq. (31) interpolation as an admitted modeling choice; my read agrees but sharpens the point: the presented DMS benchmarks are structurally incapable of validating that choice. This does not overturn the reader's CONDITIONAL verdict, because the QSS portion of the model is well supported and the paper is candid about the recombination-regime limitation. I would keep the verdict conditional and ask for the Eq. (31) check before the full concentration-dependent claim is accepted.","tokens_in":44847,"tokens_out":6634,"duration_ms":70851,"concrete_test":"Re-run the four adiabatic heat-bath benchmark cases of Sec. VII.B with the zeta-dependent factors of Eq. (31) instead of the constant factors f_k^NB = 0.5 and f_eps^NB = 0.85, keeping all other parameters fixed, and compare the time histories of T, T_v, and species mole fractions against the same DMS reference solutions. If the variable-f results match the constant-f results within the DMS scatter, the benchmark is insensitive to the concentration dependence; if they differ materially, the paper must state which version of the model is actually being benchmarked. A second useful check, if feasible, is to compute f_k^NB and f_eps^NB for zeta > 1 using a state-resolved master-equation treatment with recombination (e.g., following Ref. [65]) to test the validity of capping the factors at unity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract claims that the concentration-dependent functional form of the non-Boltzmann factors ensures depletion effects vanish at chemical equilibrium, and that the resulting model is verified against DMS. However, Sec. VII explicitly states that constant factors f_k^NB = 0.5 and f_eps^NB = 0.85 are employed throughout all MMT calculations presented in that section. This means the benchmarks exercise only the early QSS dissociation regime, where zeta < 1 and the constant factors are appropriate. They never exercise the piecewise-linear interpolation of Eq. (31), nor its transition to f = 1 as zeta approaches unity, nor the cap at unity for zeta > 1. In Sec. VI the authors explicitly call Eq. (31) a modeling choice and state they lack data to constrain the recombining regime. Therefore the benchmark evidence does not support the concentration-dependent part of the central claim. The adiabatic cases do approach equilibrium in principle, but a constant-f run cannot reveal whether the interpolation path or the switch to unity is correct, because the same constant factor multiplies forward and reverse rates and cancels at equilibrium composition, masking errors in the rate of approach. The correct equilibrium composition is guaranteed by detailed balance, but the time-dependent approach to equilibrium is not validated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":45390,"tokens_out":8173,"duration_ms":88215,"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":[{"comment":"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.","section":"Sec. VII and Eq. (31)"},{"comment":"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.","section":"Sec. V, Table 2"}],"minor_comments":[{"comment":"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.","section":"Sec. VI, Fig. 13"},{"comment":"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.","section":"Sec. II.B, Eq. (15)"},{"comment":"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.","section":"Sec. VII"},{"comment":"The caption abbreviates Marrone and Treanor as \"M and T (1963)\"; using the full author names would be clearer.","section":"Fig. 6 caption"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the paper is publishable in principle, but the abstract and conclusions overstate validation of Eq. (31), which is never benchmarked. The recommended revision should either add benchmarks that exercise the concentration-dependent factors or recalibrate the claims. Also note that the DMS benchmark data and the relaxation-time fits come from the same group and PES family; this is not disqualifying, but it means the heat-bath agreement is partly built into the input data, and the genuinely predictive element is the coupling of dissociation with vibrational energy, which is where the paper's real contribution lies."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fairly straightforward take: this is a solid, useful paper for the aerothermodynamics/CFD crowd, and the stress-test note lands. The abstract says the concentration-dependent form of the non-Boltzmann factors is 'verified' against DMS, but Sec. VII explicitly uses constant factors (0.5 and 0.85) throughout. So the interpolation in Eq. (31) – the bridge from QSS to equilibrium – is never tested. The authors themselves call it a modeling choice with no data to constrain the recombining regime. That's an honest admission, but it should be reflected in the abstract.\n\nWhat's genuinely new and good: the consolidated parameterization of the Marrone-Treanor model from the full QCT database, the temperature-dependent U parameter extracted from equilibrium vibrational energy removal, and the demonstration that a two-temperature model with three coupled inputs (rate, energy removal, relaxation time) can track DMS heat-bath results. The equilibrium Arrhenius fits are excellent (R^2 > 0.9998), and the nonequilibrium predictions are within 5-22%. The adiabatic benchmark agreement is reassuring. This is a real step forward compared to older empirical correlations based on shock-tube data.\n\nSoft spots, in order of severity. First, the non-Boltzmann factors are calibrated on four DMS reactions and applied to six. The N2+O2 and O2+N2 pairs get factors by assumption, with no DMS check. Second, as noted, Eq. (31) is unvalidated. Third, the DMS benchmarks come from the same group and the same PES family that supplied the fits; that limits how independent the validation is. It's not disqualifying – the QCT-to-MMT step is a genuine simplification, and agreement is not baked in – but independent DMS or master-equation data would strengthen the case. Fourth, the MMT model assumes rotational equilibrium, which is violated in the high-enthalpy DMS cases; the observed small leads/lags are consistent with that and are handled honestly.\n\nWho gains: people building CFD models for hypersonic entry, especially those who want a self-consistent two-temperature model rather than mismatched empirical pieces. I'd send it to review. The main requested revision would be to recalibrate the abstract's claim, add uncertainty estimates on the trajectory data, and either test or clearly label the recombination-regime interpolation as speculative. The core model is likely to be used regardless.","headline":"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.","tokens_in":45896,"tokens_out":2408,"would_cite":true,"duration_ms":25488,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-temperature dissociation model with QCT-derived rates matches direct molecular simulation for nitrogen and oxygen.","keywords":["two-temperature model","Marrone-Treanor","dissociation kinetics","quasiclassical trajectories","non-Boltzmann depletion","direct molecular simulation","hypersonic nonequilibrium air","vibrational relaxation"],"falsifier":"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.","tokens_in":44647,"feed_emoji":"🔥","tokens_out":5603,"duration_ms":51546,"temperature":0.7,"pith_summary":"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.","feed_headline":"Two-temperature model reproduces heat-bath N2 and O2 dissociation","feed_subtitle":"Rates and energy removal fit to ab initio trajectories, plus two simple correction factors, match direct molecular simulation at CFD cost.","key_machinery":"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)$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the original Marrone-Treanor preferential-dissociation two-temperature rate expression and the characteristic probability temperature concept that the MMT model modifies.","marker":"[5]"},{"why":"Provides the QCT database and the fitted Arrhenius, U, and U* parameters for the six dissociation reactions used in Table 2.","marker":"[25]"},{"why":"Supplies quasi-steady-state DMS dissociation rates and energy removal data used to calibrate the non-Boltzmann factors and to benchmark isothermal heat baths.","marker":"[34]"},{"why":"Provides DMS-derived vibrational relaxation time fits used in the Landau-Teller source term of the two-temperature model.","marker":"[37]"},{"why":"Documents the quasiclassical trajectory method and the improved N4 potential energy surface underlying the nitrogen dissociation database.","marker":"[15]"},{"why":"Demonstrates non-Boltzmann depletion in DMS nitrogen dissociation, motivating the need for correction factors in a two-temperature model.","marker":"[21]"},{"why":"Provides DMS oxygen internal energy relaxation and dissociation data used for the O2/O relaxation times and benchmark comparisons.","marker":"[32]"},{"why":"Introduces the nonequilibrium concentration ratio approach that the paper adapts to make the non-Boltzmann factors vanish at chemical equilibrium.","marker":"[63]"}],"fun_headline_variants":["This modified Marrone-Treanor model matches direct simulations for N2 and O2","QCT-derived parameters make a two-temperature model that reproduces heat-bath data","Non-Boltzmann correction factors let a cheap model mimic DMS for air dissociation","A modified Marrone-Treanor model accurately tracks N2 and O2 dissociation in heat baths","This two-temperature dissociation model is both accurate and CFD-affordable for air"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["This modified Marrone-Treanor model matches direct simulations for N2 and O2","QCT-derived parameters make a two-temperature model that reproduces heat-bath data","Non-Boltzmann correction factors let a cheap model mimic DMS for air dissociation","A modified Marrone-Treanor model accurately tracks N2 and O2 dissociation in heat baths","This two-temperature dissociation model is both accurate and CFD-affordable for air"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000381,"raw_usage":{"total_tokens":2064,"prompt_tokens":1031,"completion_tokens":1033,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":925}},"tokens_in":647,"tokens_out":1033,"duration_ms":9853,"temperature":1.0,"reasoning_tokens":925,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:58:38.796408+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Modeling and Analysis of Chemical Kin etics for Hypersonic Flows in Air,","cited_arxiv_id":null,"evidence_quote":"Provides the QCT database and the fitted Arrhenius, U, and U* parameters for the six dissociation reactions used in Table 2."}],"review_version":1}