{"id":"498e9c67-a4d3-474a-9e65-f0783f20ec87","arxiv_id":"2501.01626","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A compact H2 dissociation source-term expression based on master equation QSS theory is derived, calibrated with a pre-QSS correction factor, and paired with new rate constant fits valid from 200 to 20,000 K.","lead":"This paper derives a simplified one-temperature model for hydrogen dissociation that captures quasi-steady-state and early-time non-equilibrium effects using only the translational temperature and the degree of dissociation. The authors also compile and refit hydrogen dissociation rate constants from shock tubes, flow tubes, and master equation simulations for use in planetary entry modeling.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The pre-QSS correction η(Tt) is fitted and validated on the same isothermal/isochoric master-equation cases, all with identical initial conditions; its transferability to other densities, initial internal temperatures, or mixtures is the load-bearing unproven step.","rationale":"The QSS source-term equivalence in Eq. (41) is a genuine analytic result: substituting φA^2/χ and kr = kd,th/Keq cancels the kd,th terms exactly, and the low-temperature number-density profiles in Figure 6 verify this without invoking the pre-QSS correction. The concern is therefore not with the QSS core but with the pre-QSS extension, which is the part needed for the paper's stronger claim of a one-temperature model capturing both regimes. The paper fits η on the same master-equation data it uses for validation, and all cases share the same initial internal temperature and similar number densities, so no out-of-sample evidence exists for density, initial-condition, or mixture transferability. This is precisely the reader's weakest-assumption concern, and it justifies the CONDITIONAL verdict. The proposed out-of-sample master-equation test would settle whether η(Tt) is genuinely a function of Tt alone or whether it must depend on additional state variables such as initial internal temperature or number density.","tokens_in":37434,"tokens_out":7648,"duration_ms":78256,"concrete_test":"Run a new master-equation simulation for M = H2 at Tt = 10,000 K with initial Tr,0 = Tv,0 = 300 K and nH2,0 lowered to 1e16 cm^-3, holding Tt and the total number density fixed, using the same state-specific rates as Kim and Boyd. Then integrate Eq. (78) with the published η(Tt) Arrhenius fit, ε = 10^-3, and the published kd,nr(Tt) fit, without refitting η. Compare the predicted nH2(t) and kd(φH) against the master-equation result. If the deviation exceeds a factor of two, η is not a transferable function of Tt alone; if it agrees, the transferability concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central QSS result, Eq. (41), is an exact algebraic consequence of the QSS distribution and is well supported by the low-temperature master-equation comparisons. The load-bearing weakness is the pre-QSS correction. In Section III A, η(Tt) is obtained by least-squares fitting Eq. (75) to the master-equation kd(φH) curves that are then used as the validation data in Figures 3, 6, and 7, so the agreement is partly circular. All fitted/validated cases share a single initial rovibrational condition (Tr,0 = Tv,0 = 1000 K) and nearly fixed number densities (nH2,0 = 1e18 cm^-3; nH,0 = nHe,0 = 5e17 cm^-3), with only Tt varied. The derivation also assumes an isothermal, isochoric, inert bath, a Boltzmann pre-shock distribution, decay dominated by a single eigenmode, and α ≈ φA. Consequently, the claim that kd,pre-QSS depends only on Tt and φA, and that a one-temperature bulk-species model captures both QSS and pre-QSS dissociation, is not independently established beyond the fitted test set. The reader's weakest-assumption identification is correct: the generality of η(Tt) is the key unverified premise.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives, from the rovibrational master equations, rate expressions for diatomic dissociation in thermal non-equilibrium, and applies them to H2. The central QSS result is Eq. (41), which states that when the QSS assumption holds the chemical source term can be written as dnA/dt = 2 nA2 nM kd,nr - 2 nA^2 nM kd,nr/Keq, with kd,nr a function of Tt alone. The paper then adds a pre-QSS correction leading to Eq. (75), making kd a function of Tt and the dissociation fraction phi_A only. The model is compared with master-equation simulations of 0-D isothermal/isochoric H2 dissociation with H2, H, and He as third bodies, and an extensive literature review yields Arrhenius fits for kd,nr/Keq from 200 to 20,000 K with claimed uncertainty better than a factor of two.","tokens_in":37819,"tokens_out":2791,"duration_ms":29211,"significance":"If the QSS source-term equivalence of Eq. (41) is accepted, it is a significant and elegant result: a one-temperature bulk-species model captures the full QSS dissociation/recombination source term without fitting a recombination rate separately. The paper's algebraic derivation of Eq. (41) from the QSS distribution is internally consistent and is verified against master-equation number densities at 4,000 and 6,000 K, where pre-QSS effects are negligible. The literature review of H2 dissociation rate constants, with new fits and a factor-of-two uncertainty estimate, is also a useful contribution for ice-giant entry modeling. The pre-QSS correction, however, is the load-bearing novelty for high-temperature flows, and its validation is currently tied to the fitting procedure and to a narrow set of test conditions, so the broader claim that kd,pre-QSS depends only on Tt and phi_A is not yet independently established.","major_comments":[{"comment":"The pre-QSS correction is fitted and validated on the same data. Section III A states that the eta values are computed by a least-squares fit to the master equation results, and the resulting corrected kd curves are then shown to reproduce those same master equation curves in Fig. 3 and, after integration, in Figs. 6 and 7. The agreement in the pre-QSS region is therefore partly circular. The only genuinely independent checks are the 4,000 and 6,000 K number-density and energy profiles, where the pre-QSS contribution is small. To support the central claim that Eq. (75) is a predictive one-temperature source term, the authors should validate eta(Tt) on master-equation cases not used in the fit, e.g., different initial Tr/Tv values, different number densities, or post-shock (non-isothermal) conditions.","section":"III A, Eq. (75), Figs. 3, 6, 7"},{"comment":"Several uncontrolled approximations enter the pre-QSS derivation: the slow-mode dominance of Eq. (66), the identification of psi_infinity with psi_nr, the first-order Taylor expansion used to derive Eq. (73) while the exponential form is retained in Eq. (74), and the alpha approx phi_A substitution. These approximations are stated but not quantified, and no sensitivity analysis is provided. Since Eq. (75) is the basis for the claimed extension beyond QSS, the manuscript should either bound the error of these approximations against full master-equation solutions or clearly limit the claimed validity range to the conditions that are actually tested.","section":"II D 2, Eqs. (66), (73), (74), (75)"},{"comment":"All fitted and validated master-equation cases share nearly identical initial conditions: Tr,0 = Tv,0 = 1000 K, nH2,0 = 1e18 cm^-3, and nH,0 = nHe,0 = 5e17 cm^-3, with only Tt varied. The derivation of Eq. (75) removes the nM dependence present in Eq. (70), but this removal is not empirically verified over a range of densities or mixture compositions. Consequently, the transferability of eta(Tt) to other number densities, initial internal temperatures, and H2/H/He mixtures is the key unproven premise for CFD application, and the linear-mixture rule of Appendix A does not cure this because no mixed-bath master-equation data are presented.","section":"III A / III C, Figs. 1-7"},{"comment":"The claim that a one-temperature bulk-species model captures both QSS and pre-QSS dissociation in practical flows also relies on the source-term form of Eq. (78), in which the recombination term is written as kd,pre-QSS/Keq. The paper correctly notes in Section II C 3 that this term is not the physical kr, but the pre-QSS version of this equivalence has not been independently verified: the integrated number-density profiles in Fig. 6 at high Tt use the same fitted eta and the same isothermal, isochoric, inert-bath assumptions as the fit. A test where recombination dominates, or where the bath is not inert, would be needed before the source term is asserted to be general for entry-flow conditions.","section":"Section III C / Conclusion (Eq. (78))"}],"minor_comments":[{"comment":"The table heading reads 'Experimental Sources for High Temperature Rate Constants' but the table lists discharge-flow tube studies at Tt <= 350 K; the heading should be 'Low Temperature'.","section":"Table II"},{"comment":"The horizontal axis is labeled only with numeric tick values (0.05, 0.1, 0.15, 0.2); the axis label, presumably 10^4/Tt or similar, is missing.","section":"Fig. 4"},{"comment":"The caption says 'dashed blue and dash-dotted red lines' but the text refers to 'dash-dotted red and dashed blue lines'; the order should be made consistent so the reader knows which line is QSS and which is pre-QSS.","section":"III A, Fig. 3 caption"},{"comment":"In the paragraph describing the three Leibowitz rate sets, 'Leibowitz and Kuo128' should presumably be 'Leibowitz and Kuo129'; the current reference numbering is inconsistent with the bibliography.","section":"Appendix C"},{"comment":"The symbol QA2 is used both for the total internal partition function and, in context, for the bulk quantity; near Eq. (10) the distinction between QA2 and QA2(J,nu) should be stated explicitly to avoid confusion.","section":"Eq. (10) and Eq. (15)"}],"recommendation":"major_revision","confidential_remarks":"The QSS part of the paper is solid and the literature review is useful. My main concern is that the pre-QSS claim, which is the distinctive high-temperature contribution, is validated only on the fitting set with a single initial internal temperature and nearly fixed number densities. This is fixable with additional master-equation comparisons or by restricting the claims; I therefore do not recommend rejection, but the revision needs to address the transferability issue head-on."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a careful read. The paper does two genuinely useful things. First, it makes explicit that under the QSS assumption the dissociation/recombination source term collapses to Eq. (41), a one-temperature expression in kd,nr and Keq alone. That is a clean analytic result, and the master-equation comparisons at 4,000 and 6,000 K back it up. Second, the H2 dissociation rate compilation and the new kd,nr/Keq fits are a practical contribution; the paper makes a decent case that the old Leibowitz rates are badly off at low temperature and that a factor-of-two uncertainty bound is more honest than the order-of-magnitude or 0.25x-3x bounds used in prior entry studies. The discussion of the three inconsistent Leibowitz model variants is a nice detail.\n\nThe soft spot is the pre-QSS correction, just as the stress-test says. In Section III A, eta(Tt) is obtained by least-squares fitting Eq. (75) to the very same master-equation kd(phi_H) cases that are then shown as the validating comparisons in Figures 3, 6, and 7. All of those cases share one initial rovibrational condition (Tr,0 = Tv,0 = 1000 K) and nearly fixed number densities; only Tt varies. So the claim that kd,pre-QSS depends only on Tt and phi_A, and that a one-temperature bulk model captures both QSS and pre-QSS behavior, is not independently established beyond that fitted test set. The derivation itself also leans on several uncontrolled approximations: single-slowest-mode dominance, psi_inf = psi_nr, and alpha = phi_A. None of this makes the QSS result wrong, but it means the pre-QSS expression is presently a calibrated correlation, not a validated predictive model.\n\nTwo smaller gripes: no code or data are shipped, and the factor-of-two uncertainty on the rate fits is based on visual scatter rather than formal UQ. Both are fixable.\n\nWho is this for? Anyone doing CFD for hydrogen-helium planetary entry, and anyone who wants a compact source term for QSS dissociation. The paper deserves a serious referee: the QSS equivalence and the rate review are solid enough that the manuscript should not be desk-rejected, but the referee should push for an out-of-sample test of eta (different densities, initial temperatures, or mixtures) and for release of fitting scripts and extracted data.","headline":"The QSS source-term equivalence is a real and useful result; the pre-QSS correction is a fitted correlation that needs independent validation before it carries the load.","tokens_in":38312,"tokens_out":1539,"would_cite":true,"duration_ms":17206,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that non-equilibrium dissociation of a diatomic gas can be reduced to a one-temperature source term with two fitted rates.","keywords":["non-equilibrium dissociation","quasi-steady-state","master equations","pre-QSS correction","hydrogen dissociation","one-temperature model","rate constant fits","thermochemical nonequilibrium"],"falsifier":"Set up a 0-D isothermal, isochoric master-equation simulation of $\\rm H_2$ dissociation in a $\\rm H_2$/\\rm H$/\\rm He$ mixture at a number density outside the fitted range, say $n_{\\rm M}=10^{17}\\ \\mathrm{cm^{-3}}$ at $T_{\\rm t}=12{,}000\\ \\mathrm{K}$, and compare the predicted $n_{\\rm H}(t)$ and rovibrational energy history against Eq. (78) using the paper's $\\eta$ and $k_{\\rm d,nr}$ fits; systematic disagreement beyond the factor-of-two uncertainty would falsify the transferability of $\\eta$.","tokens_in":37228,"feed_emoji":"⚗️","tokens_out":9222,"duration_ms":78815,"temperature":0.7,"pith_summary":"This paper tries to establish that the full state-resolved chemistry of a dissociating diatomic can be collapsed, under the quasi-steady-state assumption, into a chemical source term that depends only on the translational temperature $T_{\\rm t}$ and the fraction of dissociation $\\phi_{\\rm A}$. It derives the exact QSS source term from the master equations, then adds a simple pre-QSS correction factor $\\eta(T_{\\rm t})$. For $\\rm H_2$ dissociation with the third bodies $\\rm H_2$, $\\rm H$, and $\\rm He$, the resulting two-parameter expression reproduces the number-density and rovibrational-energy profiles of detailed master-equation simulations. If correct, this gives computational fluid dynamics a cheap way to include non-equilibrium dissociation without solving additional internal-energy equations.","feed_headline":"One-temperature chemistry captures non-equilibrium H2 dissociation","feed_subtitle":"Source term depends only on gas temperature and dissociation progress, matching full master-equation simulations of H2.","key_machinery":"The central object is the quasi-steady-state decomposition of the rovibrational distribution, $$\\vec{\\psi}_{\\rm A_2} = \\vec{\\psi}_{\\rm A_2,nr}\\left(1 - \\frac{\\phi_{\\rm A}^2}{\\chi}\\right) + \\frac{\\phi_{\\rm A}^2}{\\chi}\\,\\vec{1},$$ which splits the QSS solution into a non-recombining part and a recombining part, each a function of $T_{\\rm t}$ alone. Feeding this into the aggregate rate expression makes $k_{\\rm d}$ a convex combination of $k_{\\rm d,nr}$ and $k_{\\rm d,th}$, and substituting into the source term cancels $k_{\\rm d,th}$ exactly, producing Eq. (41). For the pre-QSS region the machinery is the slowest eigenmode of the relaxation matrix: $\\eta(T_{\\rm t}) \\equiv k_{\\rm d,nr}/(2(\\lambda_1-\\lambda_0))$ converts the time-dependent relaxation into a function of $\\phi_{\\rm A}$, removing the dependence on third-body number density.","core_discovery":"On the paper's own terms, the central discovery is that the QSS chemical source term for dissociation, $$\\frac{dn_{\\rm A}}{dt} = 2n_{\\rm A_2} n_{\\rm M} k_{\\rm d,nr} - 2n_{\\rm A}^2 n_{\\rm M} \\frac{k_{\\rm d,nr}}{K_{\\rm eq}},$$ depends only on the non-recombining dissociation rate constant $k_{\\rm d,nr}(T_{\\rm t})$; the thermal-limit rate constant drops out exactly through the progress variable $\\phi_{\\rm A}^2/\\chi$. The pre-QSS extension replaces $k_{\\rm d,nr}$ by $$k_{\\rm d,pre\\text{-}QSS} = k_{\\rm d,nr}\\left[1 - (1-\\varepsilon)\\exp\\left(-\\sqrt{-\\ln(1-\\phi_{\\rm A})/\\eta}\\right)\\right],$$ so the complete source term is a function of $T_{\\rm t}$ and $\\phi_{\\rm A}$ alone, with inputs $k_{\\rm d,nr}(T_{\\rm t})$ and $\\eta(T_{\\rm t})$. The authors show that this reproduction holds for the majority of the tested master-equation cases for $\\rm H_2$ with the third bodies $\\rm H_2$, $\\rm H$, and $\\rm He$, and they use the same framework to reinterpret reported shock-tube rates as $k_{\\rm d,nr}$ data.","pith_inferences":["A genuinely independent test would run master-equation simulations at number densities and in $\\rm H_2$/\\rm H$/\\rm He$ mixtures outside the cases used to fit $\\eta$, since the paper's $\\eta$ is fitted to the same 0-D isothermal, isochoric reactor data used for validation.","The exactness of Eq. (41) is structural, not $\\rm H_2$-specific, so the same QSS source-term identity should transfer to other diatomics such as $\\rm N_2$ or $\\rm O_2$; only the fits of $k_{\\rm d,nr}$ and $\\eta$ would need to be re-established.","If the linear mixture rule is applied, the per-third-body $\\eta(T_{\\rm t})$ fits could be summed over species, but the paper notes the mixture rovibrational distribution need not match any single-bath distribution, so this is an extrapolation rather than a proven result.","For recombination-dominated flows the model may lose accuracy, because those flows overpopulate excited states while the QSS assumption here is tied to the underpopulated, dissociation-dominated distribution; the authors flag this as future work."],"forward_implications":["A one-temperature bulk-species model can capture both QSS and pre-QSS non-equilibrium dissociation for $\\rm H_2$ in the tested regimes, without separate vibrational or rotational temperatures.","The same $k_{\\rm d,nr}$ fit must be used for both the production and consumption terms in the QSS source term; the apparent recombination constant is $k_{\\rm d,nr}/K_{\\rm eq}$, not the physical $k_{\\rm r}$.","The transition from the non-recombining QSS limit to thermal equilibrium is captured implicitly by Eq. (41), so no separate fit of $k_{\\rm d,th}$ is needed in the source term.","The new fits of $k_{\\rm d,nr}/K_{\\rm eq}$ for $\\rm H_2$, $\\rm H$, and noble-gas third bodies are claimed valid from 200 to 20,000 K with uncertainty under a factor of two.","The fraction of dissociation that happens in the pre-QSS region is a function of temperature alone once $\\eta(T_{\\rm t})$ is known, independent of number density."],"supporting_citations":[{"why":"Supplies the master-equation number-density, temperature, and rate-constant data for $\\rm H_2$+$\\rm H_2$ that the QSS and pre-QSS predictions are compared against.","marker":"Kim and Boyd30"},{"why":"Supplies the master-equation simulation data for the third bodies $\\rm H$ and $\\rm He$ used to extract $k_{\\rm d}$ and to fit $\\eta(T_{\\rm t})$.","marker":"Kim68"},{"why":"Supplies master-equation data and $k_{\\rm d,nr}$ values for $\\rm H$ and $\\rm He$ at 10,000 and 16,000 K, extending the validation temperature range.","marker":"Kim et al.29"},{"why":"Provides the thermal dissociation rate constants $k_{\\rm d,th}$ for $\\rm H$ and $\\rm He$ third bodies needed to extract $k_{\\rm d}$ from the master-equation number-density profiles.","marker":"Vargas et al.71"},{"why":"Establishes the QSS formulation and the non-recombining/recombining decomposition that this paper re-derives and extends to the exact source-term identity.","marker":"Park38"},{"why":"Contributes master-equation-derived $k_{\\rm d,nr}$ data for $\\rm H_2$ and $\\rm H$ that anchor the proposed rate-constant fits at high temperature.","marker":"Schwenke89"},{"why":"Supplies master-equation $k_{\\rm d,nr}$ values for $\\rm H$+$\\rm H_2$ over a wide temperature range, constraining the uncertain $\\rm M=\\rm H$ fit.","marker":"Martin et al.91"},{"why":"Defines the baseline $\\rm H_2$ dissociation rates used in entry-flow studies that the new fits are shown to improve upon.","marker":"Leibowitz58"}],"fun_headline_variants":["H2 dissociation rates reduced to one temperature","Non-equilibrium H2 chemistry: one simple formula","H2 dissociation: master equations captured by two inputs","Single expression for H2 dissociation across regimes","One-temperature source term for H2 dissociation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the pre-QSS correction factor $\\eta(T_{\\rm t})$, fitted to a small set of 0-D isothermal and isochoric master-equation cases for three third bodies, is a function of translational temperature alone and transfers to other number densities, mixtures, and flow conditions.","fun_headline_variants_meta":{"raw":{"variants":["H2 dissociation rates reduced to one temperature","Non-equilibrium H2 chemistry: one simple formula","H2 dissociation: master equations captured by two inputs","Single expression for H2 dissociation across regimes","One-temperature source term for H2 dissociation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000292,"raw_usage":{"total_tokens":1836,"prompt_tokens":1207,"completion_tokens":629,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":823,"completion_tokens_details":{"reasoning_tokens":559}},"tokens_in":823,"tokens_out":629,"duration_ms":6621,"temperature":1.0,"reasoning_tokens":559,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:23:45.970005+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Set up a 0-D isothermal, isochoric master-equation simulation of $\\rm H_2$ dissociation in a $\\rm H_2$/\\rm H$/\\rm He$ mixture at a number density outside the fitted range, say $n_{\\rm M}=10^{17}\\ \\mathrm{cm^{-3}}$ at $T_{\\rm t}=12{,}000\\ \\mathrm{K}$, and compare the predicted $n_{\\rm H}(t)$ and rovibrational energy history against Eq. (78) using the paper's $\\eta$ and $k_{\\rm d,nr}$ fits; systematic disagreement beyond the factor-of-two uncertainty would falsify the transferability of $\\eta$.","supporting_citations":[],"review_version":1}