Pith. sign in

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 →

arxiv 2608.06143 v1 pith:4KFYIRNN submitted 2026-08-06 cond-mat.stat-mech

classification cond-mat.stat-mech MSC 82C4082C70
keywords hard-spherefluidstransportcoefficientsOnsagermatrixequationofstateself-diffusionviscositythermalconductivitythermodynamictheory
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

The paper claims that transport coefficients of a hard-sphere fluid can be computed from equilibrium thermodynamics alone, without solving kinetic equations. Starting from a general relation for the Onsager matrix, it derives a single formula for the ratio of each transport coefficient to its dilute-gas value, with all inputs coming from the equation of state. The resulting closed expressions reproduce molecular-dynamics simulation data over nearly the entire fluid density range. For self-diffusion, the thermodynamic formula reduces the mean error from about 20 percent to about 5 percent, a clear improvement over standard kinetic theory. If correct, this shows that the density dependence of transport is already encoded in equilibrium thermodynamics.

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.

Watch

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

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

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

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)
  1. [§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.
  2. [§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.
  3. [§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.
  4. [§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)
  1. [§IV, Eq. (29)] The integral in Eq. (29) contains a typographical error: it should read v_perp(t) dt, not v_perp(t), dt.
  2. [Appendix, Eq. (43)] In the sentence preceding Eq. (43), 'It is show' should be 'It is shown'.
  3. [References] Reference [30] has a typo in the title: 'A Modem Course' should be 'A Modern Course'.
  4. [§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.
  5. [§VI] In the references to simulations, 'Pieprzyket al.' appears without a space; it should be 'Pieprzyk et al.'.

Circularity Check

1 steps flagged · score 4.0 of 10

Load-bearing Onsager formula is imported from a same-group preprint, but the hard-sphere application is benchmarked against independent simulations and not fitted.

  1. 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 0 free parameters · 7 assumptions · 0 invented entities

The central formulas contain no fitted transport parameters, but they rest on the Onsager relation of Ref. [18], three structural assumptions about fluctuation amplitudes, and a kinetic-theory input to fix the self-diffusion exponent.

assumptions (7)
  • ad hoc to paper Onsager matrix relation L = (det H_id / det H) L_id (Eq. 7)
    Taken from Ref. [18], a preprint by the same group; the entire derivation depends on it and it is not rederived here.
  • domain assumption Local thermal equilibrium in cells
    Sec. III assumes cells large enough for many particles and small against perturbation wavelength.
  • ad hoc to paper Ansatz x = c T^a rho R (Eq. 14) for u, g, rho_D
    The fluctuation amplitudes of momentum and composition are assumed to have this scaling; c cancels.
  • ad hoc to paper D_x = D_tilde_x (Eq. 26)
    Asserts the transported density and its ideal part decay with the same diffusivity; not derived from dynamics.
  • ad hoc to paper Transverse velocity variance Eq. (31)
    Postulates <v_perp^2> proportional to T Z^2/(m Gamma), which fixes viscosity exponents.
  • ad hoc to paper Work-difference relation d mu_D = -delta W_D (Eq. 38)
    Phenomenological postulate for self-diffusion fluctuations; paper says it is not a general identity.
  • domain assumption Diffusion-mode dispersion and D_id = lambda_id 8/(25 rho)
    Appendix projects hydrodynamic equations onto the diffusion mode using kinetic theory; fixes exponent 5/8.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2608.06143 by the authors.

Figure 1
Figure 1. FIG. 1. Transport coefficients over their ideal value against particle density, [PITH_FULL_IMAGE:figures/full_fig_p015_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Relative error of theoretical transport coefficients, relative to numerical values, against [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

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

  2. [1]

    Chapman and T

    S. Chapman and T. G. Cowling,The Mathematical Theory of Non-Uniform Gases, 3rd ed. (Cambridge University Press, 1970)

  3. [2]

    Kinetische theorie der w¨ armeleitung, reibung und selbstdiffusion in gewissen verdichteten gasen und fl¨ ussigkeiten,

    D. Enskog, “Kinetische theorie der w¨ armeleitung, reibung und selbstdiffusion in gewissen verdichteten gasen und fl¨ ussigkeiten,” Kungl. Svenska Vetensk. Handl.63, 1 (1922)

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

  5. [4]

    Analysis of the transport coefficients for simple dense fluid: Application of the modified enskog theory,

    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)

  6. [5]

    The modified enskog equation,

    H. van Beijeren and M. H. Ernst, “The modified enskog equation,” Physica68, 437 (1973)

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

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

Show all 31 references
  1. [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)

  2. [9]

    Motions of molecules in liquids: Viscosity and diffusivity,

    J. H. Hildebrand, “Motions of molecules in liquids: Viscosity and diffusivity,” Science174, 490 (1971)

  3. [10]

    Investigations of internal friction of fluids,

    A. J. Batschinski, “Investigations of internal friction of fluids,” Z. Phys. Chem.84, 643 (1913)

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

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

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

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

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

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

  10. [17]

    Perspective: Excess-entropy scaling,

    Jeppe C. Dyre, “Perspective: Excess-entropy scaling,” J. Chem. Phys.149, 210901 (2018)

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

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

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

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

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

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

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

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

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

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

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

  22. [30]

    Frenkel and B

    D. Frenkel and B. Smit,Understanding Molecular Simulation: From Algorithms to Applica- tions, 2nd ed. (Academic Press, 2002)

  23. [31]

    L. E. Reichl,A Modem Course in Statistical Physics, 2nd ed. (Wiley, 1998). 23

Pith tools

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