TY - JOUR

T1 - Integral equation theory of thermodynamics, pair structure, and growing static length scale in metastable hard sphere and Weeks-Chandler-Andersen fluids

AU - Zhou, Yuxing

AU - Mei, Baicheng

AU - Schweizer, Kenneth S.

N1 - Funding Information:
Y.Z. and K.S.S. acknowledge support from the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division. B.M. and K.S.S. acknowledge support from DOE-BES under Grant No. DE-FG02-07ER46471 administered through the Materials Research Laboratory at UIUC.
Publisher Copyright:
© 2020 American Physical Society.

PY - 2020/4

Y1 - 2020/4

N2 - We employ the Ornstein-Zernike integral equation theory with the Percus-Yevick (PY) and modified-Verlet (MV) closures to study the equilibrium structural and thermodynamic properties of metastable monodisperse hard sphere and continuous repulsion Weeks-Chandler-Andersen (WCA) fluids under density and temperature conditions where the system is strongly overcompressed or supercooled, respectively. The theoretical results are compared to crystal-avoiding simulations of these dense monodisperse model one-component fluids. The equation of state (EOS) and dimensionless compressibility are computed using both the virial and compressibility routes. For hard spheres, the MV-based virial route EOS and dimensionless compressibility are in very good agreement with simulation for all packing fractions, much better than the PY analogs. The corresponding MV-based predictions for the static structure factor are also very good. The amplitude of density fluctuations on the local cage scale and in the long wavelength limit, and three technically different measures of the density correlation length, are studied with both closures. All five properties grow in a roughly exponential manner with density in the metastable regime up to packing fractions of 58% with no sign of saturation. The MV-based results are in good agreement with our crystal-avoiding simulations. Interestingly, the density dependences of long and short wavelength quantities are closely related. The MV-based theory is also quite accurate for the thermodynamics and structure of supercooled monodisperse WCA fluids. Overall our findings are also relevant as critical input to microscopic theories that relate the equilibrium pair correlation function or static structure factor to dynamical constraints, barriers, and activated relaxation in glass-forming liquids.

AB - We employ the Ornstein-Zernike integral equation theory with the Percus-Yevick (PY) and modified-Verlet (MV) closures to study the equilibrium structural and thermodynamic properties of metastable monodisperse hard sphere and continuous repulsion Weeks-Chandler-Andersen (WCA) fluids under density and temperature conditions where the system is strongly overcompressed or supercooled, respectively. The theoretical results are compared to crystal-avoiding simulations of these dense monodisperse model one-component fluids. The equation of state (EOS) and dimensionless compressibility are computed using both the virial and compressibility routes. For hard spheres, the MV-based virial route EOS and dimensionless compressibility are in very good agreement with simulation for all packing fractions, much better than the PY analogs. The corresponding MV-based predictions for the static structure factor are also very good. The amplitude of density fluctuations on the local cage scale and in the long wavelength limit, and three technically different measures of the density correlation length, are studied with both closures. All five properties grow in a roughly exponential manner with density in the metastable regime up to packing fractions of 58% with no sign of saturation. The MV-based results are in good agreement with our crystal-avoiding simulations. Interestingly, the density dependences of long and short wavelength quantities are closely related. The MV-based theory is also quite accurate for the thermodynamics and structure of supercooled monodisperse WCA fluids. Overall our findings are also relevant as critical input to microscopic theories that relate the equilibrium pair correlation function or static structure factor to dynamical constraints, barriers, and activated relaxation in glass-forming liquids.

UR - http://www.scopus.com/inward/record.url?scp=85084593502&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85084593502&partnerID=8YFLogxK

U2 - 10.1103/PhysRevE.101.042121

DO - 10.1103/PhysRevE.101.042121

M3 - Article

C2 - 32422713

AN - SCOPUS:85084593502

SN - 2470-0045

VL - 101

JO - Physical Review E

JF - Physical Review E

IS - 4

M1 - 042121

ER -