REVIEW 4 major objections 5 minor 31 references
Transport coefficients in hard-sphere fluids: thermodynamic versus kinetic descriptions
T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read For hard-sphere fluids, the ratio of each transport coefficient to its dilute-gas value is fixed by equilibrium thermodynamics alone, yielding parameter-free formulas for thermal conductivity, viscosity, and self-diffusion that match…
desk verdict Self-diffusion formula is Rosenfeld scaling with a prefactor, and the derivation rests on two explicit postulates that the paper itself flags as unverified. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the Onsager-matrix relation $L = (\det H_{\rm id}/\det H)\, L_{\rm id}$, where $H$ is the entropy Hessian for a fluid cell, together with the ansatz that each transported density $x$ (internal energy, transverse momentum, or species-density difference) can be written as $x = c T^a \rho R$ with $R$ depending only on density. The paper assumes the diffusivity of $x$ equals the diffusivity of its ideal part $\tilde{x} = c T^a \rho$, which lets the Onsager ratio become the transport-coefficient ratio. For viscosity and self-diffusion, two additional fluctuation hypotheses fix $a$ and $R$: the transverse-velocity variance is taken proportional to $T Z^2/(m\Gamma)$, and the differential chemical potential of the two species is taken to be set by differential mechanical work during a spontaneous fluctuation.
What would settle it
In a molecular-dynamics simulation of hard spheres, prepare a long-wavelength perturbation in the transverse velocity and separately in the species-density difference, then fit the decay rates of $x$ and of its ideal part $\tilde{x} = cT^a\rho$; if the two fitted diffusivities differ beyond statistical noise, the paper's key assumption is false and the central claim collapses.
Extended reading notes
Core claim
The central result is the formula $\sigma/\sigma_{\rm id} = \bigl((2\Gamma a+3)/(2a+3)\bigr) R$, where $\sigma$ is a transport coefficient, $\sigma_{\rm id}$ its dilute-gas value, $\Gamma = (\rho/T)(\partial\mu/\partial\rho)_T$ the thermodynamic factor, $a$ an exponent identifying the transported density, and $R$ a static density-dependent background factor. For thermal conductivity the paper sets $a=1$, $R=1$, giving $\lambda/\lambda_{\rm id} = 2\Gamma/5 + 3/5$. For viscosity it sets $a=1/2$, $R=Z/\sqrt{\Gamma}$, giving $\eta/\eta_{\rm id} = (\Gamma+3)/4 \cdot Z/\sqrt{\Gamma}$. For self-diffusion it sets $a=5/8$, $R=e^{-f_{\rm ex}/T}$, giving $D/D_{\rm id} = (5\Gamma+12)/17 \cdot e^{-f_{\rm ex}/T}$. Here $Z$ is the compressibility factor and $f_{\rm ex}$ the excess free energy per particle, both obtained from the hard-sphere equation of state. The paper argues that the agreement with simulation data—especially for self-diffusion—demonstrates that transport can be described within a purely thermodynamic framework.
Load-bearing premise
The load-bearing premise is that the density whose transport is computed and its ideal-gas counterpart decay with the same diffusivity, because their ratio is treated as a static background quantity; if that equality fails, the final ratio formula does not follow.
Editorial extensions
If this is right
- Any hard-sphere fluid's three main transport coefficients can be computed directly from its equation of state, with no fitted parameters.
- Self-diffusion, which kinetic theory underestimates by about 20 percent over the fluid range, is captured to about 5 percent by the thermodynamic formula.
- The same Onsager-matrix route should extend to other interaction potentials, because the derivation never uses the detailed collision dynamics of hard spheres beyond the equation of state.
- The theory provides a concrete reason why transport coefficients collapse onto functions of equilibrium thermodynamic state: the ratio $\sigma/\sigma_{\rm id}$ is fixed by $\Gamma$ and $R$.
- Because the formulas are compact and explicit, they can be inserted directly into hydrodynamic or heat-transfer modeling of hard-sphere fluids without solving kinetic equations.
Reading between the lines
- A direct test of the key assumption would be to compare the decay times of a perturbation in $x$ and in its ideal part $\tilde{x}$ in molecular dynamics; the paper does not report such a test, and the formulas stand or fall on it.
- Because $R = e^{-f_{\rm ex}/T}$ enters self-diffusion, the formula hints at a thermodynamic origin for the long-time memory effects that kinetic theory misses; connecting this to excess-entropy scaling would be a natural extension the author leaves implicit.
- The constant-pressure interpretation of the viscosity hypothesis could be tested by computing transverse-velocity fluctuation amplitudes under constant-volume constraints; a discrepancy near freezing would locate the observed high-density underestimate.
- If the framework transfers to soft-sphere or Lennard-Jones fluids, it would provide a parameter-free route to transport coefficients in regimes where kinetic theories are less reliable.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a thermodynamic theory for transport coefficients of hard-sphere fluids, starting from a general relation for the Onsager matrix, Eq. (7), taken from Ref. [18]. The author introduces an ansatz x = c T^a rho R for the energy density, momentum density, and species-density difference, assumes equality of the diffusivities D_x and D_tilde{x} in Eq. (26), and adds two further postulates: a transverse-velocity variance formula, Eq. (31), and a work-difference relation d mu_D = -delta W_D, Eq. (38). The resulting formulas for the thermal conductivity, viscosity, and self-diffusion coefficient, Eqs. (28), (35), and (53), are compared with molecular dynamics simulation data from Refs. [27,28] over the fluid range. The reported mean absolute percentage errors are 3.0%, 10%, and 5.2% for lambda, eta, and D, respectively, with self-diffusion substantially improved over Enskog theory.
Significance. The paper has several strengths: the final formulas are compact, contain no parameters fitted to transport data, and are benchmarked against independent simulations with transparent MAPE reporting. The improved self-diffusion prediction relative to Enskog theory is a genuine quantitative achievement if the underlying assumptions are accepted. However, the central claim that transport coefficients are derived from a purely thermodynamic framework is not fully supported, because the self-diffusion and viscosity results rely on two explicit postulates, Eqs. (31) and (38), that are not derived from equilibrium thermodynamics or from hard-sphere dynamics. The paper is transparent about these postulates and notes in Sec. VII that further verification is required, but the abstract and conclusions state the 'purely thermodynamic' claim without this qualification. The approach is best characterized at present as a thermodynamic-consistent phenomenology with two additional assumptions; validation of those assumptions by direct simulation would make the contribution substantially stronger.
major comments (4)
- [§V, Eq. (38)] The self-diffusion result, Eq. (53), rests entirely on the postulate d mu_D = -delta W_D. This relation is not derived from equilibrium thermodynamics; for a symmetric A/B mixture with equal diameters and interactions, the standard equilibrium fluctuation amplitude is r_RMS ~ 1/sqrt(N), independent of density and temperature, whereas Eq. (51) gives r proportional to T^{5/8} exp(-f_ex/T), which depends exponentially on density. Thus the density dependence of D/D_id is not a consequence of thermodynamics alone but is an external input. Since the manuscript itself acknowledges in Sec. VII that this relation 'requires further verification,' the abstract's claim of a 'purely thermodynamic' framework is stronger than the derivation supports.
- [§III, Eq. (26)] The equality D_x = D_tilde{x} is load-bearing because it converts the Onsager ratio of Eq. (21) into the transport-coefficient ratio of Eq. (27). The argument given, substituting x = R tilde{x} into the diffusion equation, shows only that tilde{x} obeys a diffusion equation with the same coefficient D_x when R is constant; it does not establish that the physical diffusivity of the ideal component tilde{x} equals D_x, since tilde{x} is not an independent conserved field in the real system. This step should be presented as an additional assumption or justified with a more explicit dynamical calculation.
- [§IV, Eq. (31)] The transverse-velocity variance formula, Eq. (31), is an ad hoc postulate. The provided physical rationale (conversion of thermal energy into collective motion at constant pressure) is qualitative and does not uniquely determine the functional form. Because the viscosity prediction, Eq. (35), depends directly on this postulate, and because the MAPE for viscosity (10% over the full fluid range) is notably worse than for thermal conductivity, the paper should include a direct numerical test of Eq. (31) from molecular dynamics trajectories (e.g., measurement of <v_perp^2> as a function of density and temperature) to support the hypothesis.
- [§II, Eq. (7)] The entire theory starts from Eq. (7), taken from Ref. [18], an arXiv preprint by the same research group that is not yet peer-reviewed. Since this relation is the foundation for all subsequent results, the manuscript should either provide a self-contained derivation of Eq. (7) in an appendix or state explicitly the conditions under which it applies to continuous fluids. As written, the reader cannot independently verify the central starting point.
minor comments (5)
- [§IV, Eq. (29)] The integral in Eq. (29) contains a typographical error: it should read v_perp(t) dt, not v_perp(t), dt.
- [Appendix, Eq. (43)] In the sentence preceding Eq. (43), 'It is show' should be 'It is shown'.
- [References] Reference [30] has a typo in the title: 'A Modem Course' should be 'A Modern Course'.
- [§III, Eq. (14)] The function R(rho) is introduced in Eq. (14) without a definition; the text should clarify immediately that R is a density-dependent factor that approaches unity in the ideal-gas limit, as done later.
- [§VI] In the references to simulations, 'Pieprzyket al.' appears without a space; it should be 'Pieprzyk et al.'.
Circularity Check
Load-bearing Onsager formula is imported from a same-group preprint, but the hard-sphere application is benchmarked against independent simulations and not fitted.
-
self citation load bearing
[Sec. II, Eq. (7); conclusions in Sec. VII]
"It was recently demonstrated [18] that the Onsager Matrix is related to its ideal counterpart, L id, through the determinant of the entropy Hessian: L= detH id/detH Lid, (7)."
All three final predictions, Eqs. (28), (35), and (53), are obtained by applying Eq. (7) to the hard-sphere fluid. Eq. (7) is not derived in this paper; it is imported from Ref. [18] (Di Muro and Hoyuelos, arXiv:2605.02579), which overlaps with the present author and is not machine-checked or otherwise independently verified here. The conclusions state: 'The derivation starts from the general expression for the Onsager matrix, Eq. (7), obtained in Ref. [18].' Thus the load-bearing mathematical premise of the derivation reduces to an unverified self-citation; if Eq. (7) were unsupported, the transport-coefficient ratios would not follow. The hard-sphere application and the two additional hypotheses are new, so the circularity is partial rather than total.
full rationale
The paper does not fit any parameter to the transport data; the constants c cancel, and the final formulas are compared with simulation results after derivation. The viscosity factor R=Z/sqrt(Gamma) and the self-diffusion factor R=exp(-f_ex/T) are obtained from explicit hypotheses, Eqs. (31) and (38), which the paper itself labels as requiring further verification (Sec. VII: 'further verification (for example, through direct numerical evaluation of the hypotheses) is still required'). That is a support gap, not a circular reduction by construction. The relation D_x = D_tilde_x in Eq. (26) is a consistency consequence of x = tilde_x R with R a static background quantity, so it is not an independently smuggled assumption. The self-diffusion exponential coincides with known excess-entropy scaling in form, but the prefactor is derived and Rosenfeld's scaling is not used as an input, so this is overlap rather than renaming. The main circularity concern is the load-bearing self-citation of Eq. (7) from the same group's previous work, which justifies a score of 4 rather than 0; the central hard-sphere results still have independent content through the new hypotheses and comparison with independent simulations.
Assumptions & free parameters
assumptions (7)
- ad hoc to paper Onsager matrix relation L = (det H_id / det H) L_id (Eq. 7)
- domain assumption Local thermal equilibrium in cells
- ad hoc to paper Ansatz x = c T^a rho R (Eq. 14) for u, g, rho_D
- ad hoc to paper D_x = D_tilde_x (Eq. 26)
- ad hoc to paper Transverse velocity variance Eq. (31)
- ad hoc to paper Work-difference relation d mu_D = -delta W_D (Eq. 38)
- domain assumption Diffusion-mode dispersion and D_id = lambda_id 8/(25 rho)
Cite this review
Pith. "Pith review of Transport coefficients in hard-sphere fluids: thermodynamic versus kinetic descriptions." pith.science (2026). https://pith.science/paper/4KFYIRNN
@misc{pith2026260806143,
author = {Pith},
title = {Pith review of: Transport coefficients in hard-sphere fluids: thermodynamic versus kinetic descriptions},
year = {2026},
howpublished = {\url{https://pith.science/paper/4KFYIRNN}},
note = {Machine review of arXiv:2608.06143}
}
abstract
A thermodynamic theory for transport coefficients in hard-sphere fluids is developed from a general expression for the Onsager matrix. The theory predicts the ratio $\sigma/\sigma_{\rm id}$, where $\sigma$ is a transport coefficient and $\sigma_{\rm id}$ denotes its dilute-gas value. This ratio depends exclusively on equilibrium thermodynamic properties and can therefore be computed directly from the equation of state. Compact expressions are obtained for the thermal conductivity, viscosity, and self-diffusion coefficient. These expressions quantitatively reproduce simulation data over almost the entire fluid range and, in the case of self-diffusion, significantly improve upon the predictions of Enskog kinetic theory. The results demonstrate that transport coefficients can be accurately described within a purely thermodynamic framework.
Figures
Reference graph
Works this paper leans on
-
[18]
A general thermodynamic approach for diffusion on a lattice
M. A. Di Muro and M. Hoyuelos, “A general thermodynamic approach for diffusion on a lattice,” arXiv:2605.02579 [cond-mat.stat-mech] (2026)
work page Pith review arXiv 2026
-
[1]
S. Chapman and T. G. Cowling,The Mathematical Theory of Non-Uniform Gases, 3rd ed. (Cambridge University Press, 1970)
work page 1970
-
[2]
D. Enskog, “Kinetische theorie der w¨ armeleitung, reibung und selbstdiffusion in gewissen verdichteten gasen und fl¨ ussigkeiten,” Kungl. Svenska Vetensk. Handl.63, 1 (1922)
work page 1922
-
[3]
Modelling of transport properties of hard sphere fluids and related systems, and its applications,
C. M. Silva and H. Liu, “Modelling of transport properties of hard sphere fluids and related systems, and its applications,” inTheory and Simulation of Hard-Sphere Fluids and Related Systems, edited by A. Mulero (Springer, 2008) p. 383
work page 2008
-
[4]
H.J.M. Hanley, R.D. McCarty, and E.G.D. Cohen, “Analysis of the transport coefficients for simple dense fluid: Application of the modified enskog theory,” Physica60, 322–356 (1972)
work page 1972
-
[5]
H. van Beijeren and M. H. Ernst, “The modified enskog equation,” Physica68, 437 (1973)
work page 1973
-
[6]
The role of repulsive forces in determining the equilibrium structure of simple liquids,
J.D. Weeks, D. Chandler, and H.C. Andersen, “The role of repulsive forces in determining the equilibrium structure of simple liquids,” J. Chem. Phys.54, 5237 (1971). 21
work page 1971
-
[7]
Shear viscosity calculated by perturbation theory and molecular dynamics for dense fluids,
F. D. C. Silva, L. A. Coelho, F. W. Tavares, and M. J. E. M. Cardoso, “Shear viscosity calculated by perturbation theory and molecular dynamics for dense fluids,” Inter. J. Quant. Chem.95, 79–87 (2003)
work page 2003
Show all 31 references
-
[8]
Corrected Enskog theory and the transport coefficients of liquids,
J. H. Dymond, “Corrected Enskog theory and the transport coefficients of liquids,” J. Chem. Phys.60, 969 (1974)
1974
-
[9]
Motions of molecules in liquids: Viscosity and diffusivity,
J. H. Hildebrand, “Motions of molecules in liquids: Viscosity and diffusivity,” Science174, 490 (1971)
1971
-
[10]
Investigations of internal friction of fluids,
A. J. Batschinski, “Investigations of internal friction of fluids,” Z. Phys. Chem.84, 643 (1913)
1913
-
[11]
Studies in Newtonian flow. II. The dependence of the viscosity of liquids on free-space,
A. K. Doolittle, “Studies in Newtonian flow. II. The dependence of the viscosity of liquids on free-space,” J. Appl. Phys.22, 1471 (1951)
1951
-
[12]
Molecular transport in liquids and glasses,
M. H. Cohen and D. Turnbull, “Molecular transport in liquids and glasses,” J. Chem. Phys. 31, 1164 (1959)
1959
-
[13]
On the free-volume model of the liquid-glass transition,
D. Turnbull and M. H. Cohen, “On the free-volume model of the liquid-glass transition,” J. Chem. Phys.52, 3038 (1970)
1970
-
[14]
On the relative roles of free volume and activation energy in the viscosity of liquids,
P. B. Macedo and T. A. Litovitz, “On the relative roles of free volume and activation energy in the viscosity of liquids,” J. Chem. Phys.42, 245 (1965)
1965
-
[15]
Relation between the transport coefficients and the internal entropy of simple systems,
Y. Rosenfeld, “Relation between the transport coefficients and the internal entropy of simple systems,” Phys. Rev. A15, 2545 (1977)
1977
-
[16]
A quasi-universal scaling law for atomic transport in simple fluids,
Y. Rosenfeld, “A quasi-universal scaling law for atomic transport in simple fluids,” J. Phys.: Condens. Matter11, 5415 (1999)
1999
-
[17]
Perspective: Excess-entropy scaling,
Jeppe C. Dyre, “Perspective: Excess-entropy scaling,” J. Chem. Phys.149, 210901 (2018)
2018
-
[19]
Dependence on the thermodynamic state of self-diffusion of pseudo-hard-sphere and lennard-jones potentials,
L. Marchioni, M. A. Di Muro, and M. Hoyuelos, “Dependence on the thermodynamic state of self-diffusion of pseudo-hard-sphere and lennard-jones potentials,” Phys. Rev. E107, 014134 (2023)
2023
-
[20]
Viscosity as the product of its ideal low- concentration value and a thermodynamic function,
L. Marchioni, M. A. Di Muro, and M. Hoyuelos, “Viscosity as the product of its ideal low- concentration value and a thermodynamic function,” Phys. Rev. E111, 064128 (2025)
2025
-
[21]
Thermodynamic theory of transport processes in semiconductors,
J. E. Parrott, “Thermodynamic theory of transport processes in semiconductors,” IEEE Trans- actions on electron devices43, 809 (1996)
1996
-
[22]
Linear-nonequilibrium thermodynamics theory for coupled heat 22 and mass transport,
Y. Demirel and S. I. Sandler, “Linear-nonequilibrium thermodynamics theory for coupled heat 22 and mass transport,” International Journal of Heat and Mass Transfer44, 2439 (2001)
2001
-
[23]
Thermodynamic theory of diffusion and thermodiffusion coefficients in mul- ticomponent mixtures,
A. A. Shapiro, “Thermodynamic theory of diffusion and thermodiffusion coefficients in mul- ticomponent mixtures,” J. Non-Equilib. Thermodyn.45, 343 (2020)
2020
-
[24]
Revised enskog theory for mie fluids: Prediction of dif- fusion coefficients, thermal diffusion coefficients, viscosities and thermal conductivities,
V. G. Jervell1 and Ø. Wilhelmsen, “Revised enskog theory for mie fluids: Prediction of dif- fusion coefficients, thermal diffusion coefficients, viscosities and thermal conductivities,” J. Chem. Phys.158, 224101 (2023)
2023
-
[25]
Equation of state for nonattracting rigid spheres,
N. F. Carnahan and K. E. Starling, “Equation of state for nonattracting rigid spheres,” J. Chem. Phys.51, 635 (1969)
1969
-
[26]
Glass transition of dense fluids of hard and compressible spheres,
L. Berthier and T. A. Witten, “Glass transition of dense fluids of hard and compressible spheres,” Phys. Rev. E80, 021502 (2009)
2009
-
[27]
The theoretical predictions correspond to the thermodynamic approach devel- oped here, Eqs
for thermal conductivity and from Pieprzyket al.[28] for viscosity and self-diffusion coefficient. The theoretical predictions correspond to the thermodynamic approach devel- oped here, Eqs. (28), (35), and (53), and are compared with the Enskog kinetic approach, Eqs. (1)–(3)....
-
[28]
A comprehensive study of the thermal conductivity of the hard sphere fluid and solid by molecular dynamics simulation,
S. Pieprzyk, A. C. Bra´ nka, D. M. Heyes, and M. N. Bannerman, “A comprehensive study of the thermal conductivity of the hard sphere fluid and solid by molecular dynamics simulation,” Phys. Chem. Chem. Phys.22, 8834 (2020)
2020
-
[29]
Thermody- namic and dynamical properties of the hard sphere system revisited by molecular dynamics simulation,
S. Pieprzyk, M. N. Bannerman, A. C. Bra´ nka, M. Chudak, and D. M. Heyes, “Thermody- namic and dynamical properties of the hard sphere system revisited by molecular dynamics simulation,” Phys. Chem. Chem. Phys.21, 6886 (2019)
2019
-
[30]
Frenkel and B
D. Frenkel and B. Smit,Understanding Molecular Simulation: From Algorithms to Applica- tions, 2nd ed. (Academic Press, 2002)
2002
-
[31]
L. E. Reichl,A Modem Course in Statistical Physics, 2nd ed. (Wiley, 1998). 23
1998
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.