Understanding in a unified manner the generic and chemically specific aspects of activated dynamics in diverse glass-forming liquids over 14 or more decades in time is a grand challenge in condensed matter physics, physical chemistry, and materials science and engineering. Large families of conceptually distinct models have postulated a causal connection with qualitatively different “order parameters” including various measures of structure, free volume, thermodynamic properties, short or intermediate time dynamics, and mechanical properties. Construction of a predictive theory that covers both the noncooperative and cooperative activated relaxation regimes remains elusive. Here, we test using solely experimental data a recent microscopic dynamical theory prediction that although activated relaxation is a spatially coupled local–nonlocal event with barriers quantified by local pair structure, it can also be understood based on the dimensionless compressibility via an equilibrium statistical mechanics connection between thermodynamics and structure. This prediction is found to be consistent with observations on diverse fragile molecular liquids under isobaric and isochoric conditions and provides a different conceptual view of the global relaxation map. As a corollary, a theoretical basis is established for the structural relaxation time scale growing exponentially with inverse temperature to a high power, consistent with experiments in the deeply supercooled regime. A criterion for the irrelevance of collective elasticity effects is deduced and shown to be consistent with viscous flow in low-fragility inorganic network-forming melts. Finally, implications for relaxation in the equilibrated deep glass state are briefly considered.

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four interrelated questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus−Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics−structure−dynamics "correlations" based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. The predicted thermodynamics−structure−dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

Suppression of the glass transition temperature, T g, in polymer thin films is of great practical importance and theoretical significance. It is widely believed that such T g suppression results from enhanced segmental mobility at the free surface. To investigate these effects, we carry out united-atom molecular dynamics simulations on free-standing polystyrene thin films. Care has been taken to ensure consistent behavior between thin films and the bulk. The dilatometric T g inferred from the density versus temperature shows substantial reduction in thin films compared to the bulk even at high cooling rates. Furthermore, we find that dynamical T g shifts, obtained by collapsing temperature-dependent short-time dynamical properties onto a master curve, vary with film thickness just like the dilatometric T g. We apply the same data collapse procedure to dynamics of segments within a given distance from the free surface to obtain the local T g(z), which reveals a mobile surface layer of about 4 nm.

In the IR spectra of the silica and silicate glasses, the shifts of the maximum intensity position of the m Si-O-Si,as band upon heating or applying mechanical stress could be attributed to changes in the distribution of bond parameters such as bond length and bond angle. Upon heating, isotropic expansion occurs and the density changes; upon applying mechanical stress, anisotropic strain is induced and a significant change in the Si-O-Si bond angle is observed. From molecular dynamics simulations of a silica glass, we show that the peak position shift correlates better with the asymmetric change in the Si-O bond length distribution, rather than the Si-O-Si bridge angle, the O-Si-O tetrahedral angle, or the density change. This new finding provides an insight into how and why the m Si-O-Si,as IR peak of soda lime silica (SLS) glass shifts upon chemical strengthening via ion exchange and thermal tempering.

We employ 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 that the system is strongly over-compressed or supercooled, respectively. The theoretical results are compared to new 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 MVbased 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 new 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.