The discovery of ultrastable glasses raises novel challenges about glassy systems. Recent experiments studied the macroscopic devitrification of ultrastable glasses into liquids upon heating but lacked microscopic resolution. We use molecular dynamics simulations to analyze the kinetics of this transformation. In the most stable systems, devitrification occurs after a very large time, but the liquid emerges in two steps. At short times, we observe the rare nucleation and slow growth of isolated droplets containing a liquid maintained under pressure by the rigidity of the surrounding glass. At large times, pressure is released after the droplets coalesce into large domains, which accelerates devitrification. This two-step process produces pronounced deviations from the classical Avrami kinetics and explains the emergence of a giant lengthscale characterizing the devitrification of bulk ultrastable glasses. Our study elucidates the nonequilibrium kinetics of glasses following a large temperature jump, which differs from both equilibrium relaxation and aging dynamics, and will guide future experimental studies.
more »
« less
Front propagation in ultrastable glasses is dynamically heterogeneous
Upon heating, ultrastable glassy films transform into liquids via a propagating equilibration front, resembling the heterogeneous melting of crystals. A microscopic understanding of this robust phenomenology is, however, lacking because experimental resolution is limited. We simulate the heterogeneous transformation kinetics of ultrastable configurations prepared using the swap Monte Carlo algorithm, thus allowing a direct comparison with experiments. We resolve the liquid–glass interface both in space and in time as well as the underlying particle motion responsible for its propagation. We perform a detailed statistical analysis of the interface geometry and kinetics over a broad range of temperatures. We show that the dynamic heterogeneity of the bulk liquid is passed on to the front that propagates heterogeneously in space and intermittently in time. This observation allows us to relate the averaged front velocity to the equilibrium diffusion coefficient of the liquid. We suggest that an experimental characterization of the interface geometry during the heterogeneous devitrification of ultrastable glassy films could provide direct experimental access to the long-sought characteristic length scale of dynamic heterogeneity in bulk supercooled liquids.
more »
« less
- Award ID(s):
- 2153944
- PAR ID:
- 10511014
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 159
- Issue:
- 11
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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.more » « less
-
Understanding the underlying nature of dynamical correlations believed to drive the bulk glass transition is a long-standing problem. Here we show that the form of spatial gradients of the glass transition temperature and structural relaxation time near an interface indeed provide signatures of the nature of relaxation in bulk glass forming liquids. We report results of long-time, large-system molecular dynamics simulations of thick glass-forming polymer films with one vapor interface, supported on a dynamically neutral substrate. We find that gradients in the glass transition temperature and logarithm of the structural relaxation time nucleated at a vapor interface exhibit two distinct regimes: a medium-ranged, large amplitude exponential gradient, followed by a long-range slowly decaying tail that can be described by an inverse power law. This behavior disagrees with multiple proposed theories of glassy dynamics but is predicted by the Elastically Collective Nonlinear Langevin Equation theory as a consequence of two coupled mechanisms: a medium-ranged interface-nucleated gradient of surface modified local caging constraints, and an interfacial truncation of a long-ranged collective elastic field. These findings support a coupled spatially local-nonlocal mechanism of activated glassy relaxation and kinetic vitrification.in both the isotropic bulk and in broken symmetry films.more » « less
-
The authors recently reported that undercooled liquid Ag and Ag–Cu alloys both exhibit a first order phase transition from the homogeneous liquid (L-phase) to a heterogeneous solid-like G-phase under isothermal evolution. Here, we report a similar L–G transition and heterogenous G-phase in simulations of liquid Cu–Zr bulk glass. The thermodynamic description and kinetic features (viscosity) of the L-G-phase transition in Cu–Zr simulations suggest it corresponds to experimentally reported liquid–liquid phase transitions in Vitreloy 1 (Vit1) and other Cu–Zr-bearing bulk glass forming alloys. The Cu–Zr G-phase has icosahedrally ordered cores versus fcc/hcp core structures in Ag and Ag–Cu with a notably smaller heterogeneity length scale Λ . We propose the L–G transition is a phenomenon in metallic liquids associated with the emergence of elastic rigidity. The heterogeneous core–shell nano-composite structure likely results from accommodating strain mismatch of stiff core regions by more compliant intervening liquid-like medium.more » « less
-
Rotational–translational decoupling in systems near T g , in which translational diffusion is apparently enhanced relative to rotation, has been observed in ensemble and single molecule experiments and has been linked to dynamic heterogeneity. Here, simulations of single molecules experiencing homogeneous diffusion and static and dynamic heterogeneous diffusion are performed to clarify the contributions of heterogeneity to such enhanced translational diffusion. Results show that time-limited trajectories broaden the distribution of diffusion coefficients in the presence of homogeneous diffusion but not when physically reasonable degrees of static heterogeneity are present. When dynamic heterogeneity is introduced, measured diffusion coefficients uniformly increase relative to input diffusion coefficients, and the widths of output distributions decrease, providing support for the idea that dynamic heterogeneity can drive apparent translational enhancement. Among simulations with dynamic heterogeneity, when the frequency of dynamic exchange is correlated with the initial diffusion coefficient, the measured diffusion coefficient behavior as a function of observation time matches that seen experimentally, the only set of simulations explored in which this occurs. Taken together with experimental results, this suggests that enhanced translational diffusion in glassy systems occurs through dynamic exchange consistent with wide underlying distributions of diffusion coefficients and exchange coupled to local spatiotemporal dynamics.more » « less
An official website of the United States government

