Elastically cooperative activated barrier hopping theory of relaxation in viscous fluids. I. General formulation and application to hard sphere fluids

Thickness-Dependent Segmental Dynamics in Supported Thin Films: Insights from a Dynamically Correlated Network Model

A large body of experimental studies shows that the local dynamics in supercooled liquids are significantly altered by spatial nanoconfinement. In a previous study, we proposed a concept of a dynamically correlated network (DCN) model, which assumes that segments in a supercooled liquid undergo cooperative rearrangements within a network-like cluster. We further demonstrated that a model modified for freestanding thin films can predict for the glass transition dynamics in atactic polystyrene (PS) films consistent with experimental results. In this study, we adapted the model to apply it to supported thin films by introducing a layer of virtual vacant segments at the free surface and virtual anchoring segments at the liquid/substrate interface. The latter segments, carrying a finite number of virtual segments, reduce mobility at the interface. We evaluated the cooperative cluster size and distribution with respect to temperature and film thickness, along with the average relaxation time and glass transition temperature Tg for supported thin films of PS. The model predicted that the thickness dependence of Tg for PS becomes stronger with increasing time scale, and this result agreed well with experimental data across different timescales from pseudothermodynamic and dynamic measurements. The results provide insights into the origin of the dynamical decoupling between pseudothermodynamic and dynamic glass transition behaviors.

Microscopic theory of the elastic shear modulus and length-scale-dependent dynamic re-entrancy phenomena in very dense sticky particle fluids

We apply the hybrid projectionless dynamic theory (hybrid PDT) formulation of the elastically collective nonlinear Langevin equation (ECNLE) activated dynamics approach to study dense fluids of sticky spheres interacting with short range attractions. Of special interest is the problem of non-monotonic evolution with short range attraction strength of the elastic modulus ("re-entrancy") at very high packing fractions far beyond the ideal mode coupling theory (MCT) nonergodicity boundary. The dynamic force constraints explicitly treat the bare attractive forces that drive transient physical bond formation, while a projection approximation is employed for the singular hard-sphere potential. The resultant interference between repulsive and attractive forces contribution to the dynamic vertex results in the prediction of localization length and elastic modulus re-entrancy, qualitatively consistent with experiments. The non-monotonic evolution of the structural (alpha) relaxation time predicted by the ECNLE theory with the hybrid PDT approach is explored in depth as a function of packing fraction, attraction strength, and attraction range. Isochronal dynamic arrest boundaries based on activated relaxation display the classic non-monotonic glass melting form. Comparisons of these results with the corresponding predictions of ideal MCT, and also the ECNLE and NLE activated theories based on projection, reveal large qualitative differences. The consequences of stochastic trajectory fluctuations on intra-cage single particle dynamics with variable strength of attractions are also studied. Large dynamical heterogeneity effects in attractive glasses are properly captured. These include a rapidly increasing amplitude of the non-Gaussian parameter with packing fraction and a non-monotonic evolution with attraction strength, in qualitative accord with recent simulations. Extension of the microscopic theoretical approach to treat double yielding in attractive glass nonlinear rheology is possible.

High-density stable glasses formed on soft substrates

Enabled by surface-mediated equilibration, physical vapour deposition can create high-density stable glasses comparable with liquid-quenched glasses aged for millions of years. Deposition is often performed at various rates and temperatures on rigid substrates to control the glass properties. Here we demonstrate that on soft, rubbery substrates, surface-mediated equilibration is enhanced up to 170 nm away from the interface, forming stable glasses with densities up to 2.5% higher than liquid-quenched glasses within 2.5 h of deposition. Gaining similar properties on rigid substrates would require 10 million times slower deposition, taking ~3,000 years. Controlling the modulus of the rubbery substrate provides control over the glass structure and density at constant deposition conditions. These results underscore the significance of substrate elasticity in manipulating the properties of the mobile surface layer and thus the glass structure and properties, allowing access to deeper states of the energy landscape without prohibitively slow deposition rates.

Effect of the nature of the solid substrate on spatially heterogeneous activated dynamics in glass forming supported films

We extend the force-level elastically collective nonlinear Langevin equation theory to treat the spatial gradients of the alpha relaxation time and glass transition temperature, and the corresponding film-averaged quantities, to the geometrically asymmetric case of finite thickness supported films with variable fluid-substrate coupling. The latter typically nonuniversally slows down motion near the solid-liquid interface as modeled via modification of the surface dynamic free energy caging constraints that are spatially transferred into the film and which compete with the accelerated relaxation gradient induced by the vapor interface. Quantitative applications to the foundational hard sphere fluid and a polymer melt are presented. The strength of the effective fluid-substrate coupling has very large consequences for the dynamical gradients and film-averaged quantities in a film thickness and thermodynamic state dependent manner. The interference of the dynamical gradients of opposite nature emanating from the vapor and solid interfaces is determined, including the conditions for the disappearance of a bulk-like region in the film center. The relative importance of surface-induced modification of local caging vs the generic truncation of the long range collective elastic component of the activation barrier is studied. The conditions for the accuracy and failure of a simple superposition approximation for dynamical gradients in thin films are also determined. The emergence of near substrate dead layers, large gradient effects on film-averaged response functions, and a weak non-monotonic evolution of dynamic gradients in thick and cold films are briefly discussed. The connection of our theoretical results to simulations and experiments is briefly discussed, as is the extension to treat more complex glass-forming systems under nanoconfinement.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Interplay between dynamic heterogeneity and interfacial gradients in a model polymer film

Glass-forming liquids exhibit long-lived, spatially correlated dynamical heterogeneity, in which some nm-scale regions in the fluid relax more slowly than others. In the nanoscale vicinity of an interface, glass-formers also exhibit the emergence of massive interfacial gradients in glass transition temperature Tg and relaxation time τ. Both of these forms of heterogeneity have a major impact on material properties. Nevertheless, their interplay has remained poorly understood. Here, we employ molecular dynamics simulations of polymer thin films in the isoconfigurational ensemble in order to probe how bulk dynamic heterogeneity alters and is altered by the large gradient in dynamics at the surface of a glass-forming liquid. Results indicate that the τ spectrum at the surface is broader than in the bulk despite being shifted to shorter times, and yet it is less spatially correlated. This is distinct from the bulk, where the τ distribution becomes broader and more spatially organized as the mean τ increases. We also find that surface gradients in slow dynamics extend further into the film than those in fast dynamics-a result with implications for how distinct properties are perturbed near an interface. None of these features track locally with changes in the heterogeneity of caging scale, emphasizing the local disconnect between these quantities near interfaces. These results are at odds with conceptions of the surface as reflecting simply a higher "rheological temperature" than the bulk, instead pointing to a complex interplay between bulk dynamic heterogeneity and spatially organized dynamical gradients at interfaces in glass-forming liquids.

Penetrant shape effects on activated dynamics and selectivity in polymer melts and networks based on self-consistent cooperative hopping theory

We generalize and apply the microscopic self-consistent cooperative hopping theory for activated penetrant dynamics in polymer melts and crosslinked networks to address the role of highly variable non-spherical molecular shape. The focus is on vastly different shaped penetrants that have identical space filling volume in order to isolate how non-spherical shape explicitly modifies dynamics over a wide range of temperature down to the kinetic glass transition temperature. The theory relates intramolecular and intermolecular structure and kinetic constraints, and reveals how different solvation packing of polymer monomers around variable shaped penetrants impact penetrant hopping. A highly shape-dependent penetrant activated relaxation, including alpha time decoupling and trajectory level cooperativity of the hopping process, is predicted in the deeply supercooled regime for relatively larger penetrants which is sensitive to whether the polymer matrix is a melt or heavily crosslinked network. In contrast, for smaller size penetrants or at high/medium temperatures the effect of isochoric penetrant shape is relatively weak. We propose an aspect ratio variable that organizes how penetrant shape influences the activated relaxation times, leading to a (near) collapse or master curve. The relative absolute values of the penetrant relaxation time (inverse hopping rate) in polymer melts versus in crosslinked networks are found to be opposite when compared at a common absolute temperature versus when they are compared at a fixed value of distance from the glass transition based on the variable Tg/T with Tg the glass transition temperature. Quantitative comparison with recent diffusion experiments on chemically complex molecular penetrants of variable shape but fixed volume in crosslinked networks reveals good agreement, and testable new predictions are made. Extension of the theoretical approach to more complex systems of high experimental interest are discussed, including applications to realize selective transport in membrane separation applications.

Understanding dynamics in coarse-grained models. III. Roles of rotational motion and translation-rotation coupling in coarse-grained dynamics

This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids.

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density-density and hydrostatic pressure-pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. These results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions

Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.

Importance of Many Particle Correlations to the Collective Debye-Waller Factor in a Single-Particle Activated Dynamic Theory of the Glass Transition

We theoretically study the importance of many body correlations on the collective Debye-Waller (DW) factor in the context of the Nonlinear Langevin Equation (NLE) single-particle activated dynamics theory of glass transition and its extension to include collective elasticity (ECNLE theory). This microscopic force-based approach envisions structural alpha relaxation as a coupled local-nonlocal process involving correlated local cage and longer range collective barriers. The crucial question addressed here is the importance of the deGennes narrowing contribution versus a literal Vineyard approximation for the collective DW factor that enters the construction of the dynamic free energy in NLE theory. While the Vineyard-deGennes approach-based NLE theory and its ECNLE theory extension yields predictions that agree well with experimental and simulation results, use of a literal Vineyard approximation for the collective DW factor massively overpredicts the activated relaxation time. The current study suggests many particle correlations are crucial for a reliable description of activated dynamics theory of model hard sphere fluids.

Thickness Dependence of Segmental Dynamics in Free-Standing Thin Films Predicted by a Dynamically Correlated Network Model

The anomalous dynamics and glass transition behaviors of supercooled liquids under nanoconfinement, such as ultrathin polymer films, have attracted much attention in recent decades. However, a complete elucidation of this mechanism has not yet been achieved. For the dynamics of bulk materials without confinement, we previously proposed a dynamically correlated network (DCN) model, which was found to agree well with the experimental data. The model assumes that segments with thermal fluctuations are dynamically correlated to their neighbors to form string-like clusters, which eventually grow into networks as temperature decreases. In this study, we applied the DCN model to nanoconfined free-standing films by using a simple cubic lattice sandwiched between two free surface layers consisting of virtual "uncorrelated" segments. The average size of DCNs at lower temperatures decreased with decreasing thickness because of confinement. This trend was associated with a decrease in the percolation temperature at which the size of DCN diverges. It was also revealed that the fractal dimension of the generated DCNs exhibits a peak with respect to temperature. The segmental relaxation time for free-standing polystyrene films was evaluated, and the predicted thickness dependence of the glass transition temperature qualitatively agreed with the experimental data. The results suggest that the concept of DCN is compatible with the dynamics of free-standing thin films.

Understanding dynamics in coarse-grained models. II. Coarse-grained diffusion modeled using hard sphere theory

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.

Combined description of pressure-volume-temperature and dielectric relaxation of several polymeric and low-molecular-weight organic glass-formers using SL-TS2 approach

We apply our recently-developed mean-field "SL-TS2" (two-state Sanchez-Lacombe) model to simultaneously describe dielectric α-relaxation time, τ_α, and pressure-volume-temperature (PVT) data in four polymers (polystyrene, poly(methylmethacrylate), poly(vinyl acetate) and poly(cyclohexane methyl acrylate)) and four organic molecular glass formers (ortho-terphenyl, glycerol, PCB-62, and PDE). Previously, it has been shown that for all eight materials, the Casalini-Roland thermodynamical scaling, τ_α = f(Tvγsp) (where T is temperature and v_sp is specific volume) is satisfied (R. Casalini and C. M. Roland, Phys. Rev. E, 2004, 69(6), 62501). It has also been previously shown that the same scaling emerges naturally (for sufficiently low pressures) within the "SL-TS2" framework (V. V. Ginzburg, Soft Matter, 2021, 17, 9094-9106). Here, we fit the ambient pressure curves for the relaxation time and the specific volume as functions of temperature for the eight materials and observe a good agreement between theory and experiment. We then use the Casalini-Roland scaling to convert those results into "master curves", thus enabling predictions of relaxation times and specific volumes at elevated pressures. The proposed approach can be used to describe other glass-forming materials, both low-molecular-weight and polymeric.

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S0(T) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S0(T) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S0(T) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Screening and collective effects in randomly pinned fluids: a new theoretical framework

We propose a theoretical framework for the dynamics of bulk isotropic hard-sphere systems in the presence of randomly pinned particles and apply this theory to supercooled water to validate it. Structural relaxation is mainly governed by local and non-local activated process. As the pinned fraction grows, a local caging constraint becomes stronger and the long range collective aspect of relaxation is screened by immobile obstacles. Different responses of the local and cooperative motions results in subtle predictions for how the alpha relaxation time varies with pinning and density. Our theoretical analysis for the relaxation time of water with pinned molecules quantitatively well describe previous simulations. In addition, the thermal dependence of relaxation for unpinned bulk water is also consistent with prior computational and experimental data.

Temperature Dependence of Structural Relaxation in Glass-Forming Liquids and Polymers

Understanding the microscopic mechanism of the transition of glass remains one of the most challenging topics in Condensed Matter Physics. What controls the sharp slowing down of molecular motion upon approaching the glass transition temperature Tg, whether there is an underlying thermodynamic transition at some finite temperature below Tg, what the role of cooperativity and heterogeneity are, and many other questions continue to be topics of active discussions. This review focuses on the mechanisms that control the steepness of the temperature dependence of structural relaxation (fragility) in glass-forming liquids. We present a brief overview of the basic theoretical models and their experimental tests, analyzing their predictions for fragility and emphasizing the successes and failures of the models. Special attention is focused on the connection of fast dynamics on picosecond time scales to the behavior of structural relaxation on much longer time scales. A separate section discusses the specific case of polymeric glass-forming liquids, which usually have extremely high fragility. We emphasize the apparent difference between the glass transitions in polymers and small molecules. We also discuss the possible role of quantum effects in the glass transition of light molecules and highlight the recent discovery of the unusually low fragility of water. At the end, we formulate the major challenges and questions remaining in this field.

An alternative, dynamic density functional-like theory for time-dependent density fluctuations in glass-forming fluids

We propose an alternative theory for the relaxation of density fluctuations in glass-forming fluids. We derive an equation of motion for the density correlation function that is local in time and is similar in spirit to the equation of motion for the average non-uniform density profile derived within the dynamic density functional theory. We identify the Franz-Parisi free energy functional as the non-equilibrium free energy for the evolution of the density correlation function. An appearance of a local minimum of this functional leads to a dynamic arrest. Thus, the ergodicity breaking transition predicted by our theory coincides with the dynamic transition of the static approach based on the same non-equilibrium free energy functional.

Local-Average Free Volume Correlates with Dynamics in Glass Formers

Glass formers exhibit a pronounced slowdown in dynamics, accompanied by progressive heterogeneity as they approach the glass transition. There is intense debate over whether the dramatic slowdown is caused by dynamical heterogeneity and whether the enhanced dynamical heterogeneity originates from structural causes. However, the connection between dynamical heterogeneity and the spatial distribution of the single-particle free volume (a purely static structural quantity) was found to be rather weak, which raises the question of whether dynamic heterogeneity has a purely structural origin. Here, by introducing the concept of local-average free volume, we present numerical evidence that long-time dynamic heterogeneity shows significantly enhanced correlation with the average local free volume over a length scale of a few neighboring shells. Our results resolve the long-standing controversy about whether free volume plays an important role in particle rearrangements associated with the activated hopping relaxation. The concept of "local average" can be applied to other local structural descriptors to better correlate with dynamic heterogeneity in glass-forming liquids.

The dimensional evolution of structure and dynamics in hard sphere liquids

The formulation of the mean-field infinite-dimensional solution of hard sphere glasses is a significant milestone for theoretical physics. How relevant this description might be for understanding low-dimensional glass-forming liquids, however, remains unclear. These liquids indeed exhibit a complex interplay between structure and dynamics, and the importance of this interplay might only slowly diminish as dimension d increases. A careful numerical assessment of the matter has long been hindered by the exponential increase in computational costs with d. By revisiting a once common simulation technique involving the use of periodic boundary conditions modeled on D_d lattices, we here partly sidestep this difficulty, thus allowing the study of hard sphere liquids up to d = 13. Parallel efforts by Mangeat and Zamponi [Phys. Rev. E 93, 012609 (2016)] have expanded the mean-field description of glasses to finite d by leveraging the standard liquid-state theory and, thus, help bridge the gap from the other direction. The relatively smooth evolution of both the structure and dynamics across the d gap allows us to relate the two approaches and to identify some of the missing features that a finite-d theory of glasses might hope to include to achieve near quantitative agreement.

Temperature dependence of aging dynamics in highly non-equilibrium model polymer glasses

A central feature of the non-equilibrium glassy “state” is its tendency to age toward equilibrium, obeying signatures identified by Kovacs over 50 years ago. The origin of these signatures, their fate far from equilibrium and at high temperatures, and the underlying nature of the glassy “state” far from equilibrium remain unsettled. Here, we simulate physical aging of polymeric glasses, driven much farther from equilibrium and at much higher temperatures than possible in experimental melt-quenched glasses. While these glasses exhibit Kovacs’ signatures of glassy aging at sufficiently low temperatures, these signatures disappear above the onset TA of non-Arrhenius equilibrium dynamics, suggesting that TA demarcates an upper bound to genuinely glassy states. Aging times in glasses after temperature up-jumps are found to obey an Arrhenius law interpolating between equilibrium dynamics at TA and at the start of the temperature up-jump, providing a zero-parameter rule predicting their aging behavior and identifying another unrecognized centrality of TA to aging behavior. This differs qualitatively from behavior of our glasses produced by temperature down-jumps, which exhibit a fractional power law decoupling relation with equilibrium dynamics. While the Tool–Narayanaswamy–Moynihan model can predict the qualitative single-temperature behavior of these systems, we find that it fails to predict the disappearance of Kovacs signatures above TA and the temperature dependence of aging after large temperature up-jumps. These findings highlight a need for new theoretical insights into the aging behavior of glasses at ultra-high fictive temperatures and far from equilibrium.

Confinement Effects on the Spatially Inhomogeneous Dynamics in Metallic Glass Films

This work develops the elastically collective nonlinear Langevin equation theory to investigate, for the first time, the glassy dynamics in capped metallic glass thin films. Finite-size effects on the spatial gradient of structural relaxation time and glass transition temperature (T_g) are calculated at different temperatures and vitrification criteria. Molecular dynamics is significantly slowed down near rough solid surfaces, and the dynamics at location far from the interfaces is sped up. In thick films, the mobility gradient normalized by the bulk value obeys the double-exponential form since interference effects between two surfaces are weak. Reducing the film thickness induces a strong dynamic coupling between two surfaces and flattens the relaxation gradient. The normalized gradient of the glass transition temperature is independent of vitrification time scale criterion and can be fitted by a superposition function as the films are not ultrathin. The local fragility is found to remain unchanged with location. This finding suggests that one can use Angell plots of bulk relaxation time and the T_g spatial gradient to characterize glassy dynamics in metallic glass films. Our computational results agree well with experimental data and simulation.

Design of a homologous series of molecular glassformers

We design and synthesize a set of homologous organic molecules by taking advantage of facile and tailorable Suzuki cross coupling reactions to produce triarylbenzene derivatives. By adjusting the number and the arrangement of conjugated rings, the identity of heteroatoms, lengths of fluorinated alkyl chains, and other interaction parameters, we create a library of glassformers with a wide range of properties. Measurements of the glass transition temperature (T_g) show a power-law relationship between T_g and molecular weight (M_W), with of the molecules, with an exponent of 0.3 ± 0.1, for T_g values spanning a range of 300-450 K. The trends in indices of refraction and expansion coefficients indicate a general increase in the glass density with M_W, consistent with the trends observed in T_g variations. A notable exception to these trends was observed with the addition of alkyl and fluorinated alkyl groups, which significantly reduced T_g and increased the dynamical fragility (which is otherwise insensitive to M_W). This is an indication of reduced density and increased packing frustrations in these systems, which is also corroborated by the observations of the decreasing index of refraction with increasing length of these groups. These data were used to launch a new database for glassforming materials, glass.apps.sas.upenn.edu.

Construction of a quantitative relation between structural relaxation and dynamic heterogeneity by vibrational dynamics in glass-forming liquids and polymers

The structural relaxation slows down drastically upon approaching the glass transition, accompanied by the significant growth of dynamic heterogeneity. The fundamental question of elusiveness and interest is whether there exists an underlying quantitative relationship between structural relaxation and dynamic heterogeneity. Here, we reveal that b̃ which is related to the reduced mean square displacements to overcome the energy barriers of activated jumps, instead of the kinetic fragility m, is the genuine key parameter connecting dynamic heterogeneity with structural relaxation for varying types of glass formers. Furthermore, based on the dependence of dynamic heterogeneity on the Debye-Waller factor we obtained a direct quantitative relation between dynamic heterogeneity and structural relaxation is built for different glass-forming liquids. More importantly, a scaling collapse of structural relaxation and dynamic heterogeneity is achieved by the important parameter b̃. These results are of fundamental and critical importance for developing a unified theory of glassy dynamics.

Long Wavelength Thermal Density Fluctuations in Molecular and Polymer Glass-Forming Liquids: Experimental and Theoretical Analysis under Isobaric Conditions

We establish via an in-depth analysis of experimental data that the dimensionless compressibility (proportional to the dimensionless amplitude of long wavelength thermal density fluctuations) of one-component normal and supercooled liquids of chemically complex nonpolar and weakly polar molecules and polymers follows extremely well a surprisingly simple and general temperature dependence over an exceptionally wide range of pressures and temperatures. A theoretical basis for this behavior is shown to exist in the venerable van der Waals model and its more modern interpretations. Although associated hydrogen-bonding (and to a lesser degree strongly polar) liquids display modestly more complex behavior, rather simple temperature and pressure dependences are also discovered. A new approach to collapse the temperature- and pressure-dependent dimensionless compressibility data onto a master curve is formulated that differs from the empirical thermodynamic scaling approach. As a practical matter, we also find that the dimensionless compressibility scales well as an inverse power law with temperature with an exponent that is system dependent and decreases with pressure. At very high pressures and low temperatures, the thermal liquid behavior appears to approach (but not reach) a repulsion-dominated random close packing limit. All these findings are relevant to our recent theoretical work on the problem of activated relaxation and vitrification of supercooled molecular and polymeric liquids.

Activation free energy gradient controls interfacial mobility gradient in thin polymer films

We examine the mobility gradient in the interfacial region of substrate-supported polymer films using molecular dynamics simulations and interpret these gradients within the string model of glass-formation. No large gradients in the extent of collective motion exist in these simulated films, and an analysis of the mobility gradient on a layer-by-layer basis indicates that the string model provides a quantitative description of the relaxation time gradient. Consequently, the string model indicates that the interfacial mobility gradient derives mainly from a gradient in the high-temperature activation enthalpy ΔH₀ and entropy ΔS₀ as a function of depth z, an effect that exists even in the high-temperature Arrhenius relaxation regime far above the glass transition temperature. To gain insight into the interfacial mobility gradient, we examined various material properties suggested previously to influence ΔH₀ in condensed materials, including density, potential and cohesive energy density, and a local measure of stiffness or u²(z)^-3/2, where u²(z) is the average mean squared particle displacement at a caging time (on the order of a ps). We find that changes in local stiffness best correlate with changes in ΔH₀(z) and that ΔS₀(z) also contributes significantly to the interfacial mobility gradient, so it must not be neglected.

Combined description of polymer PVT and relaxation data using a dynamic "SL-TS2" mean-field lattice model

We develop a combined model to describe the pressure-volume-temperature (PVT) thermodynamics and the α- and β-relaxation time dynamics in glass-forming amorphous materials. The PVT results are described using a two-state modification of the Sanchez-Lacombe equation of state (SL-EoS). The minimization of the Sanchez-Lacombe free energy expression allows one to calculate the density (or specific volume) and the "solid fraction" as a function of temperature and pressure. The solid fraction, ψ(T,P), is then substituted into the TS2 mean-field model (V. V. Ginzburg, Soft Matter, 2020, 16, 810), and the equilibrium relaxation times and the glass transition temperature are calculated. Finally, we use a dynamic model of a Tool-Narayanaswami-Moynihan (TNM)-type to describe the density relaxation as a function of the cooling rate. We applied the new framework to describe experimental data for polystyrene and poly(methylmethacrylate) and found a good qualitative and quantitative agreement.

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical "slaving onset" stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

Nature of dynamic gradients, glass formation, and collective effects in ultrathin freestanding films

Molecular, polymeric, colloidal, and other classes of liquids can exhibit very large, spatially heterogeneous alterations of their dynamics and glass transition temperature when confined to nanoscale domains. Considerable progress has been made in understanding the related problem of near-interface relaxation and diffusion in thick films. However, the origin of "nanoconfinement effects" on the glassy dynamics of thin films, where gradients from different interfaces interact and genuine collective finite size effects may emerge, remains a longstanding open question. Here, we combine molecular dynamics simulations, probing 5 decades of relaxation, and the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory, addressing 14 decades in timescale, to establish a microscopic and mechanistic understanding of the key features of altered dynamics in freestanding films spanning the full range from ultrathin to thick films. Simulations and theory are in qualitative and near-quantitative agreement without use of any adjustable parameters. For films of intermediate thickness, the dynamical behavior is well predicted to leading order using a simple linear superposition of thick-film exponential barrier gradients, including a remarkable suppression and flattening of various dynamical gradients in thin films. However, in sufficiently thin films the superposition approximation breaks down due to the emergence of genuine finite size confinement effects. ECNLE theory extended to treat thin films captures the phenomenology found in simulation, without invocation of any critical-like phenomena, on the basis of interface-nucleated gradients of local caging constraints, combined with interfacial and finite size-induced alterations of the collective elastic component of the structural relaxation process.

Experimental test of a predicted dynamics-structure-thermodynamics connection in molecularly complex glass-forming liquids

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.

Vitrification is a spontaneous non-equilibrium transition driven by osmotic pressure

Persistent dynamics in colloidal glasses suggest the existence of a non-equilibrium driving force for structural relaxation during glassy aging. But the implicit assumption in the literature that colloidal glasses form within the metastable state bypasses the search for a driving force for vitrification and glassy aging and its connection with a metastable state. The natural relation of osmotic pressure to number-density gradients motivates us to investigate the osmotic pressure as this driving force. We use dynamic simulation to quench a polydisperse hard-sphere colloidal liquid into the putative glass region while monitoring structural relaxation and osmotic pressure. Following quenches to various depths in volume fractionϕ(whereϕ_RCP≈ 0.678 for 7% polydispersity), the osmotic pressure overshoots its metastable value, then decreases with age toward the metastable pressure, driving redistribution of coordination number and interparticle voids that smooths structural heterogeneity with age. For quenches to 0.56 ⩽ϕ⩽ 0.58, accessible post-quench volume redistributes with age, allowing the glass to relax into a strong supercooled liquid and easily reach a metastable state. At higher volume fractions, 0.59 ⩽ϕ< 0.64, this redistribution encounters a barrier that is subsequently overcome by osmotic pressure, allowing the system to relax toward the metastable state. But forϕ⩾ 0.64, the overshoot is small compared to the high metastable pressure; redistribution of volume stops as particles acquire contacts and get stuck, freezing the system far from the metastable state. Overall, the osmotic pressure drives structural rearrangements responsible for both vitrification and glassy age-relaxation. The connection of energy, pressure, and structure identifies the glass transition, 0.63 <ϕ_g⩽ 0.64. We leverage the connection of osmotic pressure to energy density to put forth the mechanistic view that relaxation of structural heterogeneity in colloidal glasses occurs via individual particle motion driven by osmotic pressure, and is a spontaneous energy minimization process that drives the glass off and back to the metastable state.

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

We employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant-matrix pair correlation function and the penetrant-matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Theory of Pressure-Induced Rejuvenation and Strain Hardening in Metallic Glasses

We theoretically investigate high-pressure effects on the atomic dynamics of metallic glasses. The theory predicts compression-induced rejuvenation and the resulting strain hardening that have been recently observed in metallic glasses. Structural relaxation under pressure is mainly governed by local cage dynamics. The external pressure restricts the dynamical constraints and slows down the atomic mobility. In addition, the compression induces a rejuvenated metastable state (local minimum) at a higher energy in the free-energy landscape. Thus, compressed metallic glasses can rejuvenate and the corresponding relaxation is reversible. This behavior leads to strain hardening in mechanical deformation experiments. Theoretical predictions agree well with experiments.

MD simulations of charged binary mixtures reveal a generic relation between high- and low-temperature behavior

Experimental studies of the glassy slowdown in molecular liquids indicate that the high-temperature activation energy E_∞ of glass-forming liquids is directly related to their glass transition temperature T_g. To further investigate such a possible relation between high- and low-temperature dynamics in glass-forming liquids, we analyze the glassy dynamics of binary mixtures using molecular dynamics simulations. We consider a binary mixture of charged Lennard-Jones particles and vary the partial charges of the particles and, thus, the high-temperature activation energy and the glass transition temperature of the system. Based on previous results, we introduce a phenomenological model describing relaxation times over the whole temperature regime from high temperatures to temperatures well inside the supercooled regime. By investigating the dynamics of both particle species on molecular and diffusive length scales along isochoric and isobaric pathways, we find a quadratic charge dependence of both E_∞ and T_g, resulting in an approximately constant ratio of both quantities independent of the underlying observable, the thermodynamic ensemble, and the particle species, and this result is robust against the actual definition of T_g. This generic relation between the activation energy and the glass transition temperature indicates that high-temperature dynamics and the glassy slowdown are related phenomena, and the knowledge of E_∞ may allow us to approximately predict T_g.