On the modelling of polyatomic molecules in kinetic theory

Cermics); LJLL (UMR\_7598)); Maria Groppi (UNIPR; Marzia Bisi (UNIPR); Thomas Borsoni (LJLL (UMR\_7598)

REVIEW 1 major objections 4 minor 1 cited by

On the modelling of polyatomic molecules in kinetic theory

T0 review · 1 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict A sound, self-aware cheat sheet for polyatomic kinetic modelling; the stress-test's Eq. (13) objection is a misreading. read the letter →

arxiv 2501.18207 v1 pith:6GKBE7UC submitted 2025-01-30 math.AP

Marzia Bisi (UNIPR) , Thomas Borsoni (LJLL (UMR_7598) , CERMICS) , Maria Groppi (UNIPR , LJLL (UMR_7598)) This is my paper · ORCID

classification math.AP MSC 82B4082D0576P05

keywords polyatomicgaseskineticmodellinginternalstatesenergylevelsquantilesDSMCnon-polytropicBoltzmannequation

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 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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

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

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$.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the 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.

The 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.

Now 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.

This 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.

Referee Report

1 major / 4 minor

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).

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 (1)

[§2.2.2, Eq. (13)] 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.

minor comments (4)

[§2.2.2, just above Eq. (13)] 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.
[§2.3.2, Eq. (16)] 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.
[§3.2] 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.
[§2.2.2, Eq. (13)] 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⌋.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the level and quantile descriptions are linked by explicit measure-theoretic change-of-variable identities, and no fitted quantity is relabelled as a prediction.

full rationale

The paper is a dictionary between three equivalent mathematical encodings of the same molecular internal-state information. Equation (2) defines the energy-level measure as the image of the state measure under the grounded energy function, and Eq. (4) defines the quantile function as the generalized inverse of the cumulative energy law; the integral equalities in Sections 1.2.1 and 1.2.2 are change-of-variable identities and are explicitly presented as correspondences, not as empirical predictions. The equilibrium Maxwellian forms in Table 4 and the thermodynamic quantities in Eqs. (18)-(23) are standard statistical-mechanics computations from the stated partition function and the definition of the Maxwellian measure, not outputs fitted to data. The mesoscopic equivalence criterion in Section 3.2, namely that energy-based reduction is lossless when cross-sections depend on internal states only through the relevant energy variables, is the only load-bearing external ingredient; it belongs to the authors' prior framework [5], which is an independent published mathematical work with its own derivations and stated assumptions. Under the reviewing rules, such a citation is real evidence rather than circularity, because the theorem is not parameterized by the present paper's outputs and is not being used to disguise a fitted quantity as a prediction. No parameter is fitted, no uniqueness claim is imported to forbid alternatives by self-citation alone, and no empirical result is renamed as a new derivation. The algebraic error in the worked quantile formula (13) noted by the skeptical reading is a correctness defect in an illustrative formula, not a circularity, and it does not change the verdict.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No parameters are fitted to data; J, Δε, and the molecular mass are physical inputs of the molecular model. The ledger contains four modelling assumptions from Sections 1.1, 1.5.1, 2.1, and 3.2. No new physical entities are introduced; the quantile variable q is a mathematical reparametrization of the same state space.

assumptions (4)

domain assumption The internal state of a molecule is represented by a measure space (E, μ) equipped with a measurable internal energy function ε with finite minimum ε0.
Section 1.1 opens the modelling framework; every derived description, including the energy law and quantile function, is built on this representation.
standard math The measures involved (μ, μ̄ε, Lebesgue on quantiles) are σ-finite so that image measures, product measures, and convolutions are well defined.
Needed for equations (2), (10), and the separated descriptions; the paper notes the σ-finite requirement in Section 1.5.1.
domain assumption In the diatomic example, rotation is a classical rigid rotor and vibration is an independent quantum harmonic oscillator.
Section 2.1; this independence and the choice of energy functions produce the explicit formulas (12)-(16) and (22)-(23).
domain assumption State-based and energy-based mesoscopic models are equivalent only when collision cross-sections depend on internal states only through energy.
Section 3.2 states this as the condition under which energy-based descriptions lose no information at kinetic level; it is the load-bearing assumption for the equivalence claims.

how reviews work

0 comments

Cite this review

Pith. "Pith review of On the modelling of polyatomic molecules in kinetic theory." pith.science (2026). https://pith.science/paper/6GKBE7UC

@misc{pith2026250118207,
  author       = {Pith},
  title        = {Pith review of: On the modelling of polyatomic molecules in kinetic theory},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6GKBE7UC}},
  note         = {Machine review of arXiv:2501.18207}
}

read the original abstract

This communication is both a pedagogical note for understanding polyatomic modelling in kinetic theory and a ''cheat sheet'' for a series of corresponding concepts and formulas. We explain, detail and relate three possible approaches for modelling the polyatomic internal structure, that are: the internal states approach, well suited for physical modelling and general proofs, the internal energy levels approach, useful for analytic studies and corresponding to the common models of the literature, and the internal energy quantiles approach, less known while being a powerful tool for particle-based numerical simulations such as Direct Simulation Monte-Carlo (DSMC). This note may in particular be useful in the study of non-polytropic gases.

Figures

Figures reproduced from arXiv: 2501.18207 by the authors.

**Figure 2.** Plots of energy laws and quantile functions associated with the discussed [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗

**Figure 3.** Plot of the total number of degrees of freedom and heat capacity at constant [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗

Discussion (0). Continue with ORCID 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

A kinetic model for polyatomic gas with quasi-resonant collisions leading to bi-temperature relaxation processes
math.AP 2025-06 conditional novelty 6.0 of 10

A new Boltzmann-type model for polyatomic gases with quasi-resonant collisions yields a Landau-Teller relaxation between kinetic and internal temperatures.

Reference graph

Works this paper leans on

12 extracted references · 11 canonical work pages · cited by 1 Pith paper

[5]

Borsoni, M

T. Borsoni, M. Bisi, and M. Groppi. A general framework for the kinetic modelling of polyatomic gases. Commun. Math. Phys., 393(1):215–266, 2022

2022
[1]

Arima, T

T. Arima, T. Ruggeri, and M. Sugiyama. Extended thermodynamics of rarefied polyatomic gases: 15-field theory incorporating relaxation processes of molecular rotation and vibration. Entropy, 20(4):301, 2018

work page 2018
[2]

M. Bisi, M. Groppi, and G. Spiga. A kinetic model for bimolecular chemical reactions. In Kinetic Methods for Nonconservative and Reacting Systems (G. Toscani, ed.), Quaderni di Matematica, vol. 16, pp. 1-143, Aracne Editrice, 2005

work page 2005
[3]

Borgnakke and P

C. Borgnakke and P. S. Larsen. Statistical collision model for Monte Carlo simulation of polyatomic gas mixture. J. Comput. Phys., 18:405–420, 1975

work page 1975
[4]

T. Borsoni. Contributions autour de l’équation de Boltzmann et certaines de ses variantes . PhD thesis, Sorbonne Université, 2024

work page 2024
[6]

Desvillettes

L. Desvillettes. Sur un modèle de type Borgnakke-Larsen conduisant à des lois d’énergie non linéaires en température pour les gaz parfaits polyatomiques. Ann. Fac. Sci. Toulouse Math., 6:257–262, 1997

work page 1997
[7]

I. M. Gamba and M. Pavi´c- ˇColi´c. On the Cauchy problem for Boltzmann equation modelling a polyatomic gas. J. Math. Phys., 64(1):013303, 2023

work page 2023
[8]

Giovangigli

V . Giovangigli. Multicomponent Flow Modeling. Birkhäuser, Boston, 1999

work page 1999

Show all 12 references

[9]

Groppi and G

M. Groppi and G. Spiga. Kinetic approach to chemical reactions and inelastic transitions in a rarefied gas. J. Math. Chem., 26:197–219, 1999

1999
[10]

A.R. Leach. Molecular modelling: principles and applications. Pearson education, 2001

2001
[11]

Pavi ´c- ˇColi´c and S

M. Pavi ´c- ˇColi´c and S. Simi ´c. Kinetic description of polyatomic gases with temperature- dependent specific heats. Phys. Rev. Fluids, 7(8):083401, 2022

2022
[12]

C. S. Wang Chang and G. E. Uhlenbeck. Transport phenomena in polyatomic gases. Research Rep. CM-681. University of Michigan Engineering, 1951

1951

Pith tools

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

[1] [5]

Borsoni, M

T. Borsoni, M. Bisi, and M. Groppi. A general framework for the kinetic modelling of polyatomic gases. Commun. Math. Phys., 393(1):215–266, 2022

2022

[2] [1]

Arima, T

T. Arima, T. Ruggeri, and M. Sugiyama. Extended thermodynamics of rarefied polyatomic gases: 15-field theory incorporating relaxation processes of molecular rotation and vibration. Entropy, 20(4):301, 2018

work page 2018

[3] [2]

M. Bisi, M. Groppi, and G. Spiga. A kinetic model for bimolecular chemical reactions. In Kinetic Methods for Nonconservative and Reacting Systems (G. Toscani, ed.), Quaderni di Matematica, vol. 16, pp. 1-143, Aracne Editrice, 2005

work page 2005

[4] [3]

Borgnakke and P

C. Borgnakke and P. S. Larsen. Statistical collision model for Monte Carlo simulation of polyatomic gas mixture. J. Comput. Phys., 18:405–420, 1975

work page 1975

[5] [4]

T. Borsoni. Contributions autour de l’équation de Boltzmann et certaines de ses variantes . PhD thesis, Sorbonne Université, 2024

work page 2024

[6] [6]

Desvillettes

L. Desvillettes. Sur un modèle de type Borgnakke-Larsen conduisant à des lois d’énergie non linéaires en température pour les gaz parfaits polyatomiques. Ann. Fac. Sci. Toulouse Math., 6:257–262, 1997

work page 1997

[7] [7]

I. M. Gamba and M. Pavi´c- ˇColi´c. On the Cauchy problem for Boltzmann equation modelling a polyatomic gas. J. Math. Phys., 64(1):013303, 2023

work page 2023

[8] [8]

Giovangigli

V . Giovangigli. Multicomponent Flow Modeling. Birkhäuser, Boston, 1999

work page 1999