{"id":"2585adfb-305c-434e-9aae-455c2b187261","arxiv_id":"2501.18207","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A concise, self-contained guide showing how state-based, energy-level, and energy-quantile descriptions of polyatomic internal energy are related and when each should be used.","lead":"This paper lays out three equivalent mathematical descriptions of the internal energy of polyatomic molecules for kinetic theory: internal states, energy levels, and energy quantiles. It shows when each is most useful, with explicit formulas for a rigid-rotor and quantum-oscillator diatomic molecule and for macroscopic quantities.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The worked diatomic quantile formula in Eq. (13) is mathematically wrong on (0,∞), so the advertised DSMC 'cheat sheet' does not yet implement the claimed equivalence.","rationale":"The paper is a pedagogical note whose central value lies in precise, directly usable formulas. The abstract measure/quantile framework is correct, and the equilibrium equivalence between state-, energy-level-, and energy-quantile descriptions is standard once the general inverse-quantile map is accepted. The reader's weakest assumption, that the out-of-equilibrium equivalence requires cross-sections to depend on internal states only through energy, is indeed an assumption and not derived here; for a note whose scope is to relate existing frameworks, this is acceptable rather than a fatal gap. However, the reader's own observation that Eq. (13) is undefined at q = 0 understates the problem: the formula is wrong for every q > 0 in the first interval and is systematically off by roughly one vibrational quantum at each subsequent threshold. Because the paper explicitly advertises itself as a cheat sheet for DSMC, an incorrect energy-quantile map in the main worked diatomic example is a material defect: following Eq. (13) would assign incorrect post-collisional internal energies, thereby changing the sampled physics. The fix is local and algebraic, and the general method of Section 1.2.2 survives, so this should not lead to rejection; it should lead to conditional acceptance with a corrected Eq. (13), verified against the image-measure identity. The Section 3.2 cross-section condition remains a stated assumption, but it does not represent an internal inconsistency in this expository paper.","tokens_in":10619,"tokens_out":14269,"duration_ms":141870,"concrete_test":"Recompute Eq. (13) directly from definition (4). Set J = 2π and Δε = 1, so μ̄ε([0,I)) = ∫₀ᴵ ⌈s⌉ ds. Evaluate at q = 0.5 and q = 1: the correct values are 0.5 and 1, whereas Eq. (13) is undefined at q = 0.5 and gives 2 at q = 1. More generally, invert the piecewise formula to get F←(q) = Δε(q̂/(ℓ+1) + ℓ/2) with ℓ = floor((√(1+8q̂)−1)/2), and verify the image-measure identity ∫ φ(F←(q)) dq = ∫ φ(I) dμ̄ε(I) for test functions on intervals between thresholds.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Eq. (13) in §2.2.2 does not define the total-energy quantile function. From Eqs. (11)–(12), with x = I/Δε and q̂ = J q/(2πΔε), one obtains q̂ = ∫₀ˣ ⌈s⌉ ds. On each interval x ∈ [ℓ, ℓ+1), this gives q̂ = (ℓ+1)x − ℓ(ℓ+1)/2, hence x = q̂/(ℓ+1) + ℓ/2. Eq. (13) instead prints x = q̂/ℓ + (ℓ+1)/2. At q̂ = 0.5, ℓ = 0 and the formula is undefined; at q̂ = 1, it gives x = 2 instead of the correct x = 1; and at every triangular threshold q̂ = ℓ(ℓ+1)/2 it gives x = ℓ+1 instead of x = ℓ. The quantile function is therefore discontinuous and off by about one vibrational quantum. This is load-bearing because the paper's practical claim is precisely that a modeller may adopt the uniform energy-quantile coordinate for DSMC without changing the physics: a code using Eq. (13) would sample the wrong internal energies. The general construction in §1.2.2 is mathematically sound, and the cross-section condition in §3.2 is a clearly labelled assumption rather than a derived molecular fact; the concrete failure is in the worked example meant to make the cheat sheet usable.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript is a pedagogical note on three equivalent ways to encode the internal energy of a polyatomic molecule in kinetic theory: the internal-state description (E, μ, ε), the energy-level description (R+, μ̄ε, ε0), and the energy-quantile description ((0,qmax), Lebesgue, F←). Section 1 develops the measure-theoretic correspondence via the image measure and the quantile transform; Section 2 works out the example of a diatomic molecule with classical rotation and quantum harmonic vibration, giving explicit total and separated energy laws and quantile functions; Section 3 lists macroscopic moment formulas, discusses the mesoscopic equivalence of state-based and energy-based Boltzmann descriptions, and gives equilibrium formulas for the Maxwellian measure, the number of degrees of freedom, and the heat capacity. The central claim is that a modeller may move freely between the descriptions and, in particular, may use the uniform energy-quantile coordinate for DSMC, provided the collision cross-section depends on internal states only through energy (or only through the separate rotational and vibrational energies).","tokens_in":10887,"tokens_out":11445,"duration_ms":108188,"significance":"The note fills a useful niche: it translates the authors' prior framework [5] into a compact set of formulas relating the state, level, and quantile descriptions and explicitly addresses DSMC users. The general construction in Section 1.2 is mathematically clean and parameter-free, and the equilibrium quantities (19)-(23) are standard statistical-mechanics results assembled correctly from the given measures. The pedagogical value is real, and the paper makes no fitted or invented parameters. However, the worked total-energy quantile formula in Eq. (13), which is the advertised 'cheat sheet' item, is incorrect as printed; because the practical DSMC claim rests on this formula, the manuscript needs revision before it can be accepted.","major_comments":[{"comment":"The total-energy quantile function as printed is not the inverse of the cumulative function from Eqs. (11)-(12). With x = I/Δε and q̂ = J q/(2πΔε), Eqs. (11)-(12) give q̂ = ∫₀ˣ ⌈s⌉ ds, and on the interval x ∈ [ℓ, ℓ+1) this yields q̂ = (ℓ+1)x − ℓ(ℓ+1)/2, hence x = q̂/(ℓ+1) + ℓ/2. Eq. (13) instead prints x = q̂/ℓ + (ℓ+1)/2 (with ℓ the floor of the indicated square-root expression). At q̂ = 1 the printed formula gives x = 2 instead of the correct x = 1; at q̂ = 2 it gives x = 3 instead of 1.5; and on q̂ ∈ (0,1) it is undefined because ℓ = 0. This is not a pointwise convention issue: a DSMC code using Eq. (13) would sample internal energies that are systematically too large by about one vibrational quantum. The correct formula is F←(q) = Δε(q̂/(ℓ+1) + ℓ/2) with ℓ = floor((√(1+8q̂)−1)/2), which also matches the breakpoints claimed in §2.4. Please correct Eq. (13), its definition of ℓ, and the corresponding plot in Fig. 2b.","section":"§2.2.2, Eq. (13)"}],"minor_comments":[{"comment":"The sentence 'for any q ≥ 0' conflicts with the domain (0,qmax) of the quantile function defined in Eq. (4); q = 0 should be excluded or treated separately, since F←(0) = 0 is obtained only as an infimum.","section":"§2.2.2, just above Eq. (13)"},{"comment":"With the half-open interval convention in (3)-(4), the quantile function for μ̄εvib = ∑_{n≥0} δ_{nΔε} is Δε(⌈q⌉−1) for q > 0, not Δε⌊q⌋ at integer q; for example F←(1) = 0, not Δε. The discrepancy occurs on a measure-zero set of q for continuous sampling, but since the note is intended as a formula reference, Eq. (16) should be corrected.","section":"§2.3.2, Eq. (16)"},{"comment":"The assertion that the state-based and energy-based Boltzmann models are equivalent if and only if the cross-section depends on internal states only through their energy is stated without proof or a citation in this note; given the pedagogical character, a reference to [5] or a short proof sketch would make the claim verifiable.","section":"§3.2"},{"comment":"The displayed definition of ℓ(q̂) is malformed in the manuscript ('⌊ 4 q̂ √(1+8q̂) + 1 ⌋'); please ensure it is typeset as ℓ(q̂) = ⌊(√(1+8q̂) − 1)/2⌋.","section":"§2.2.2, Eq. (13)"}],"recommendation":"major_revision","confidential_remarks":"This is a short pedagogical note whose value depends entirely on the correctness of its working formulas. The general measure-theoretic construction is sound, and the issue in Eq. (13) is local and fixable; the paper should be accepted only after the quantile formulas are corrected and the figures regenerated. I see no concerns about citation practices or scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a clear, honest cheat sheet, not a new research contribution. The three-way correspondence among state-based, energy-level, and energy-quantile descriptions is already in the authors' earlier CMP paper, and they say so. What is new is the explicit quantile formulation for a diatomic molecule and a very clean statement of when the reduced energy descriptions remain equivalent out of equilibrium.\n\nThe paper does well what it sets out to do. The image-measure and quantile constructions are standard and correctly applied. The diatomic example in Section 2 works out the energy law and the quantile functions; Eqs. (12)-(16) and the equilibrium formulas (22)-(23) are direct computations. Table 4 and the macroscopic formulas are a useful reference. The caveat in Section 3.2, that equivalence holds only if the cross-section depends on internal states via energy, is stated clearly and not hidden.\n\nNow the soft spots. The stress-test claim that Eq. (13) is wrong does not hold up. The derivation from (11)-(12) gives the quantile as Δε( q̂/(ℓ+1) + ℓ/2 ), which is well-behaved at all q>0; the printed expression is easy to misread because the denominator parentheses are missing in the plain text, but the formula itself is correct. A referee might ask the authors to bracket it. The DSMC-relevance argument in Section 1.4 is plausible but not demonstrated: the uniform measure is the right reason, but no actual DSMC code or numerical experiment is shown. That is fine for a note, but it keeps the practical claim at the level of motivation. The equivalence condition in Section 3.2 is an assumption about collision cross-sections; the authors say so, so no marks against them.\n\nThis note is for kinetic theorists who want a compact dictionary before choosing a coordinate for their next model or simulation. It deserves a serious referee, not because it is deep but because it is correct and genuinely useful. I would accept it as a communication after a light revision that brackets Eq. (13) and adds one sentence on the status of the cross-section assumption.","headline":"A sound, self-aware cheat sheet for polyatomic kinetic modelling; the stress-test's Eq. (13) objection is a misreading.","tokens_in":11466,"tokens_out":5041,"would_cite":false,"duration_ms":42363,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B40","82D05","76P05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Three descriptions of a polyatomic molecule's internal structure are equivalent at equilibrium, and out of equilibrium exactly when collisions depend only on internal energy.","keywords":["polyatomic gases","kinetic modelling","internal states","internal energy levels","internal energy quantiles","DSMC","non-polytropic gases","Boltzmann equation"],"falsifier":"Take a polyatomic gas with a collision cross-section that depends on internal states only through total internal energy, run the state-resolved and the total-energy-based Boltzmann equations (or DSMC) from the same non-equilibrium initial condition, and compare the relaxation of the velocity and internal-energy distributions; the paper predicts exact agreement, so any measurable discrepancy would falsify the claimed equivalence. For the converse direction, use a cross-section that distinguishes internal states of equal energy, such as different rotational sub-states, and check whether the reduced energy-only model still matches the full state-based model; a match would refute the 'only if' claim.","tokens_in":10426,"feed_emoji":"⚛️","tokens_out":17798,"duration_ms":153056,"temperature":0.7,"pith_summary":"This paper is a working dictionary for the three standard ways to label a polyatomic molecule's internal degrees of freedom in kinetic theory: the full internal state $\\zeta$, the internal energy level $I$, and the internal energy quantile $q$. It establishes that the two energy-based descriptions are exact changes of variable of the state-based one — the energy law $\\mu_{\\bar\\varepsilon}$ is the image of the state-space measure under the grounded energy map, and the energy quantile function $F^\\leftarrow_{\\mu_{\\bar\\varepsilon}}$ is the generalized inverse of that law — so all single-molecule integrals, macroscopic fields, and equilibrium Maxwellian measures coincide across the three. The paper then states when the reduction is safe at the kinetic level: at equilibrium the descriptions are always equivalent, and away from equilibrium the state-based and energy-based Boltzmann models are equivalent if and only if the cross-section depends on internal states only through their energy; in a separated ro-vibrational model the same criterion applies to the separate rotational and vibrational energies. A complete worked example — a rigid rotor plus quantum harmonic oscillator diatomic molecule — gives the explicit energy density $\\varphi_{\\mathrm{total}}$, the quantile function $F^\\leftarrow$, the temperature-dependent number of internal degrees of freedom, and the heat capacity. The practical payoff for a modeller is that the quantile description turns the internal variables into a one-dimensional space carrying a uniform measure, exactly the coordinate a Direct Simulation Monte-Carlo code wants.","feed_headline":"One energy number can replace full molecular states in kinetic models","feed_subtitle":"Swap detailed molecular states for one uniform energy coordinate in DSMC without changing the physics.","key_machinery":"The central object is the change-of-variable chain connecting a state space $(E,\\mu)$ with internal energy $\\varepsilon$ to the energy line and the quantile line. The links are the push-forward energy law $\\mu_{\\bar\\varepsilon}=\\bar\\varepsilon_\\#\\mu$, the cumulative law $F_{\\mu_{\\bar\\varepsilon}}(I)=\\mu_{\\bar\\varepsilon}([0,I))$, and its left generalized inverse $F^\\leftarrow_{\\mu_{\\bar\\varepsilon}}(q)=\\inf\\{I\\ge 0 : F_{\\mu_{\\bar\\varepsilon}}(I)\\ge q\\}$. The identity $\\int_E \\varphi(\\varepsilon(\\zeta))\\,d\\mu(\\zeta)=\\int_{\\mathbb{R}_+}\\varphi(I+\\varepsilon_0)\\,d\\mu_{\\bar\\varepsilon}(I)=\\int_{(0,q_{\\max})}\\varphi(F^\\leftarrow_{\\mu_{\\bar\\varepsilon}}(q)+\\varepsilon_0)\\,dq$ is what makes the descriptions interchangeable at the level of integrals. At the kinetic level the same machinery carries over as long as the collision cross-section factors through the energy map, which is precisely the condition under which the Boltzmann collision operator agrees in all three representations.","core_discovery":"On the paper's own terms, the discovery is that the state-based description $(E,\\mu,\\varepsilon)$, the energy-level description $(\\mathbb{R}_+,\\mu_{\\bar\\varepsilon})$, and the energy-quantile description $((0,q_{\\max}),\\mathrm{Leb})$ with quantile map $F^\\leftarrow_{\\mu_{\\bar\\varepsilon}}$ are not rival models but exact repackagings of the same single-molecule data. The paper establishes the dictionary through two change-of-variable identities in Section 1, shows that the macroscopic fields and equilibrium Maxwellian measures in Tables 3 and 4 coincide, and then pins down the mesoscopic limitation: at thermodynamic equilibrium the approaches are always equivalent, while out of equilibrium the equivalence holds if and only if the cross-section depends on internal states only through their energy, with the separated ro-vibrational version holding when it depends only on the separate rotational and vibrational energies. The worked diatomic example computes the energy density $\\varphi_{\\mathrm{total}}(I)=\\frac{2\\pi}{J}\\lceil I/\\Delta\\epsilon\\rceil$, the quantile function in (13), and the temperature-dependent degrees of freedom and heat capacity in (22)--(23), all in closed form.","pith_inferences":["The same dictionary suggests a design recipe for coarse-grained Boltzmann kernels: integrate a state-to-state cross-section over all states sharing a given internal energy, and if the result depends only on that energy the reduced kernel is exact; this offers a concrete route from molecular data to energy-only collision models.","The 'if and only if' criterion turns state-resolved scattering measurements into a direct test of whether an energy-only kinetic model is admissible: gases with strong orientation-dependent or mode-specific transition rules should show deviations from the reduced models in non-equilibrium relaxation.","Because the quantile representation is an exact change of variables, the same elementary identity should also carry over to polyatomic BGK and Fokker-Planck models, not only Boltzmann equations, whenever the relaxation operator respects the energy-only condition.","The closed-form diatomic quantile function implies that uniform sampling of $q$ makes the physical internal energy grow like the square root of the quantile at large $q$, so a uniform mesh in $q$ resolves low energies finely and spreads out at high energies; DSMC implementations should sample $q$ uniformly and compute the physical energy from $F^\\leftarrow$."],"forward_implications":["A Direct Simulation Monte-Carlo code may use the internal energy quantile $q$ as the particle's internal coordinate; because the quantile space carries the Lebesgue measure, moving a particle to a new quantile does not change its statistical weight.","Analytic results written in the energy-level language transfer verbatim to the quantile language, since the macroscopic quantities in Table 3 and the Maxwellian measures in Table 4 are identical under the change of variables.","For a diatomic molecule with independent rotation and vibration, the state-based model reduces to the separated two-energy description whenever the cross-section depends separately on rotational and vibrational energy, and to the total-energy description only when it depends on total internal energy.","The temperature-dependent internal degrees of freedom $\\delta(T)$ and heat capacity $c_V(T)$ are expressed through the partition function in a way that is manifestly independent of which of the three descriptions is used, giving a consistency check for any candidate non-polytropic model.","Because the equivalence is conditional on the energy-only cross-section property, the note warns modellers that energy-based or quantile-based reductions hide state-specific physics whenever collisions are sensitive to more than energy."],"supporting_citations":[{"why":"Supplies the general state-based framework for kinetic modelling of polyatomic gases that this note condenses into energy-based descriptions.","marker":"[5]"},{"why":"Its Introduction, Section 1.2 is the source of the state-to-energy correspondence and the equilibrium-equivalence statement used in Section 3.2.","marker":"[4]"},{"why":"Introduces the continuous internal-energy description whose energy law is recovered when the push-forward measure has a density.","marker":"[3]"},{"why":"Shows that the choice of measure on the energy variable controls the temperature dependence of internal energy, which is the non-polytropic setting the note addresses.","marker":"[6]"},{"why":"Introduces the discrete energy-level description for vibrational modes, recovered when the energy law is a discrete measure.","marker":"[12]"},{"why":"Provides the discrete internal-energy kinetic model in the literature that the energy-level description corresponds to.","marker":"[2]"},{"why":"Connects kinetic descriptions of polyatomic gases to temperature-dependent specific heats, motivating the practical interest of the dictionary.","marker":"[11]"}],"fun_headline_variants":["Three polyatomic models become one when scattering is energy-only","Energy-only scattering unifies polyatomic kinetic models","One energy coordinate replaces molecular state lists in DSMC","Kinetic theory: internal states, levels, quantiles are equivalent","Polyatomic gas modeling: three descriptions, one physics core"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole reduction stands on the premise that a molecule's collisions care only about how much internal energy it carries, not about which particular internal state carries it; if real cross-sections depend on orientation, mode, or other quantum labels beyond energy, the energy-level and energy-quantile descriptions lose information away from equilibrium.","fun_headline_variants_meta":{"raw":{"variants":["Three polyatomic models become one when scattering is energy-only","Energy-only scattering unifies polyatomic kinetic models","One energy coordinate replaces molecular state lists in DSMC","Kinetic theory: internal states, levels, quantiles are equivalent","Polyatomic gas modeling: three descriptions, one physics core"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000318,"raw_usage":{"total_tokens":1767,"prompt_tokens":888,"completion_tokens":879,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":504,"completion_tokens_details":{"reasoning_tokens":799}},"tokens_in":504,"tokens_out":879,"duration_ms":9209,"temperature":1.0,"reasoning_tokens":799,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T00:17:32.570733+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a polyatomic gas with a collision cross-section that depends on internal states only through total internal energy, run the state-resolved and the total-energy-based Boltzmann equations (or DSMC) from the same non-equilibrium initial condition, and compare the relaxation of the velocity and internal-energy distributions; the paper predicts exact agreement, so any measurable discrepancy would falsify the claimed equivalence. For the converse direction, use a cross-section that distinguishes internal states of equal energy, such as different rotational sub-states, and check whether the reduced energy-only model still matches the full state-based model; a match would refute the 'only if' claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Its Introduction, Section 1.2 is the source of the state-to-energy correspondence and the equilibrium-equivalence statement used in Section 3.2."},{"cited_title":"Borgnakke and P","cited_arxiv_id":null,"evidence_quote":"Introduces the continuous internal-energy description whose energy law is recovered when the push-forward measure has a density."},{"cited_title":"Desvillettes","cited_arxiv_id":null,"evidence_quote":"Shows that the choice of measure on the energy variable controls the temperature dependence of internal energy, which is the non-polytropic setting the note addresses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the discrete energy-level description for vibrational modes, recovered when the energy law is a discrete measure."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the discrete internal-energy kinetic model in the literature that the energy-level description corresponds to."},{"cited_title":"Pavi ´c- ˇColi´c and S","cited_arxiv_id":null,"evidence_quote":"Connects kinetic descriptions of polyatomic gases to temperature-dependent specific heats, motivating the practical interest of the dictionary."}],"review_version":1}