A comprehensive study of the thermal conductivity of the hard sphere fluid and solid by Molecular Dynamics simulation

Slawomir Pieprzyk* (Corresponding Author), Arkadiusz C. Branka, David M. Heyes, Marcus Campbell Bannerman

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

14 Citations (Scopus)
6 Downloads (Pure)


This work reports a new set of hard sphere (HS) thermal conductivity coefficient, λ, data obtained by Molecular Dynamics (MD) computer simulation, over a density range covering the dilute fluid to near the close-packed solid, and for a large number of particles (up to N = 13 1072) and long simulation times. The N-dependence of the thermal conductivity is shown to be proportional to N−2/3 to a good approximation over a wide range of system sizes, which enabled λ values in the thermodynamic limit to be predicted accurately. The fluid and solid λ can be represented well by the Enskog theory (ET) formula, λE, times a density-dependent correction term, which is close to unity for the fluid and practically constant for the solid. The convergence of the MD λ data back towards ET in the metastable fluid starts just above the freezing density. For the HS solid and dense fluid it was found that the thermal conductivity is nearly linear in pressure, as has been observed experimentally for a number of solids. Simple excess entropy scaling over the higher density fluid phase region was found, and Rosenfeld's exponential relationship can be fitted to the simulation data for the solid to a high degree of accuracy. The simulation analysis has revealed a number of new trends in the behaviour of the HS thermal conductivity which could be useful in building more accurate models for heat conduction in experimental systems.

Original languageEnglish
Pages (from-to)8834-8845
Number of pages12
JournalPhysical Chemistry Chemical Physics
Issue number16
Early online date14 Apr 2020
Publication statusPublished - 28 Apr 2020




Dive into the research topics of 'A comprehensive study of the thermal conductivity of the hard sphere fluid and solid by Molecular Dynamics simulation'. Together they form a unique fingerprint.

Cite this