Phonon gas model for thermal conductivity of dense, strongly interacting liquids
Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.