Theoretical analysis of the structure, thermodynamics, and shear elasticity of deeply metastable hard sphere fluids

The structure, thermodynamics, and slow activated dynamics of the equilibrated metastable regime of glass-forming fluids remain a poorly understood problem of high theoretical and experimental interest. We apply a highly accurate microscopic equilibrium liquid state integral equation theory, in conjunction with naïve mode coupling theory of particle localization, to study in a unified manner the structural correlations, thermodynamic properties, and dynamic elastic shear modulus in deeply metastable hard sphere fluids. Distinctive behaviors are predicted including divergent inverse critical power laws for the contact value of the pair correlation function, pressure, and inverse dimensionless compressibility, and a splitting of the second peak and large suppression of interstitial configurations of the pair correlation function. The dynamic elastic modulus is predicted to exhibit two distinct exponential growth regimes with packing fraction that have strongly different slopes. These thermodynamic, structural, and elastic modulus results are consistent with simulations and experiments. Perhaps most unexpectedly, connections between the amplitude of long wavelength density fluctuations, dimensionless compressibility, local structure, and the dynamic elastic shear modulus have been theoretically elucidated. These connections are more broadly relevant to understanding the slow activated relaxation and mechanical response of colloidal suspensions in the ultradense metastable region and deeply supercooled thermal liquids in equilibrium.