Surface-Bulk Interplay in Vapor-Deposited Glasses: Crossover Length and the Origin of Front Transformation

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.

Two-step devitrification of ultrastable glasses

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.

Surface premelting and melting of colloidal glasses

The nature of liquid-to-glass transition is a major puzzle in science. A similar challenge exists in glass-to-liquid transition, i.e., glass melting, especially for the poorly investigated surface effects. Here, we assemble colloidal glasses by vapor deposition and melt them by tuning particle attractions. The structural and dynamic parameters saturate at different depths, which define a surface liquid layer and an intermediate glassy layer. The power-law growth of both layers and melting front behaviors at different heating rates are similar to crystal premelting and melting, suggesting that premelting and melting can be generalized to amorphous solids. The measured single-particle kinetics reveal various features and confirm theoretical predictions for glass surface layer.

The role of intramolecular relaxations on the structure and stability of vapor-deposited glasses

Stable glasses (SGs) are formed through surface-mediated equilibration (SME) during physical vapor deposition (PVD). Unlike intermolecular interactions, the role of intramolecular degrees of freedom in this process remains unexplored. Here, using experiments and coarse-grained molecular dynamics simulations, we demonstrate that varying dihedral rotation barriers of even a single bond, in otherwise isomeric molecules, can strongly influence the structure and stability of PVD glasses. These effects arise from variations in the degree of surface mobility, mobility gradients, and mobility anisotropy, at a given deposition temperature ( T dep). At high T dep, flexible molecules have access to more configurations, which enhances the rate of SME, forming isotropic SGs. At low T dep, stability is achieved by out of equilibrium aging of the surface layer. Here, the poor packing of rigid molecules enhances the rate of surface-mediated aging, producing stable glasses with layered structures in a broad range of T dep. In contrast, the dynamics of flexible molecules couple more efficiently to the glass layers underneath, resulting in reduced mobility and weaker mobility gradients, producing unstable glasses. Independent of stability, the flattened shape of flexible molecules can also promote in-plane orientational order at low T dep. These results indicate that small changes in intramolecular relaxation barriers can be used as an approach to independently tune the structure and mobility profiles of the surface layer and, thus, the stability and structure of PVD glasses.

Vapor deposition rate modifies anisotropic glassy structure of an anthracene-based organic semiconductor

We control the anisotropic molecular packing of vapor-deposited glasses of ABH113, a deuterated anthracene derivative with promise for future organic light emitting diode materials, by changing the deposition rate and substrate temperature at which they are prepared. We find that at substrate temperatures from 0.65 T_g to 0.92 T_g, the deposition rate significantly modifies the orientational order in the vapor-deposited glasses as characterized by x-ray scattering and birefringence. Both measures of anisotropic order can be described by a single deposition rate-substrate temperature superposition (RTS). This supports the applicability of the surface equilibration mechanism and generalizes the RTS principle from previous model systems with liquid crystalline order to non-mesogenic organic semiconductors. We find that vapor-deposited glasses of ABH113 have significantly enhanced density and thermal stability compared to their counterparts prepared by liquid-cooling. For organic semiconductors, the results of this study provide an efficient guide for using the deposition rate to prepare stable glasses with controlled molecular packing.

Microscopic Structural Evolution during Ultrastable Metallic Glass Formation

By decreasing the rate of physical vapor deposition, ZrCuAl metallic glasses with improved stability and mechanical performances can be formed, while the microscopic structural mechanisms remain unclear. Here, with scanning transmission electron microscopy and high-energy synchrotron X-ray diffraction, we found that the metallic glass deposited at a higher rate exhibits a heterogeneous structure with compositional fluctuations at a distance of a few nanometers, which gradually disappear on decreasing the deposition rate; eventually, a homogeneous structure is developed approaching ultrastability. This microscopic structural evolution suggests the existence of the following two dynamical processes during ultrastable metallic glass formation: a faster diffusion process driven by the kinetic energy of the depositing atoms, which results in nanoscale compositional fluctuations, and a slower collective relaxation process that eliminates the compositional and structural heterogeneity, equilibrates the deposited atoms, and strengthens the local atomic connectivity.

Glass Dynamics Deep in the Energy Landscape

When a liquid is cooled, progress down the energy landscape is arrested near the glass transition temperature T_g. In principle, lower energy states can be accessed by waiting for further equilibration, but the rough energy landscape of glasses quickly leads to kinetics on geologically slow time scales below T_g. Over the past decade, progress has been made probing deeper into the energy landscape via several techniques. By looking at bulk and surface diffusion, using layered deposition that promotes equilibration, imaging glass surfaces with faster dynamics below T_g, and optically exciting glasses, experiments have moved into a regime of ultrastable, low energy glasses that was difficult to access in the past. At the same time, both simulations and energy landscape theory based on a random first order transition (RFOT) have tackled systems that include surfaces, optical excitation, and interfacial dynamics. Here we review some of the recent experimental work, and how energy landscape theory illuminates glassy dynamics well below the glass transition temperature by making direct connections between configurational entropy, energy landscape barriers, and the resulting dynamics.

Spatial Heterogeneities in Structural Temperature Cause Kovacs' Expansion Gap Paradox in Aging of Glasses

Volume and enthalpy relaxation of glasses after a sudden temperature change has been extensively studied since Kovacs' seminal work. One observes an asymmetric approach to equilibrium upon cooling versus heating and, more counterintuitively, the expansion gap paradox, i.e., a dependence on the initial temperature of the effective relaxation time even close to equilibrium when heating. Here, we show that a distinguishable-particle lattice model can capture both the asymmetry and the paradox. We quantitatively characterize the energetic states of the particle configurations using a physical realization of the fictive temperature called the structural temperature, which, in the heating case, displays a strong spatial heterogeneity. The system relaxes by nucleation and expansion of warmer mobile domains having attained the final temperature, against cooler immobile domains maintained at the initial temperature. A small population of these cooler regions persists close to equilibrium, thus explaining the paradox.

Nucleation and Growth of the Supercooled Liquid Phase Control Glass Transition in Bulk Ultrastable Glasses

We report the anomalous bulk transformation of vapor deposited stable glasses into the liquid state. The transformation proceeds through two competing parallel processes: partial rejuvenation of the stable glass and nucleation and growth of liquid patches within the glass. The kinetics of the transformation extracted from heat capacity curves after isothermal runs is dominated by the heterogeneous nucleation and growth process that initiates at preexisting seeds and propagates radially at a velocity proportional to the alpha relaxation time. Remarkably, the distance between the activation seeds is independent of temperature within experimental uncertainty and amounts to several micrometers, a value in close agreement with the crossover length for TPD glasses. We speculate the initiation sites for the transformation of the glass into the supercooled liquid are localized regions of lower stability (or density).

Front-Mediated Melting of Isotropic Ultrastable Glasses

Ultrastable vapor-deposited glasses display uncommon material properties. Most remarkably, upon heating they are believed to melt via a liquid front that originates at the free surface and propagates over a mesoscopic crossover length, before crossing over to bulk melting. We combine swap Monte Carlo with molecular dynamics simulations to prepare and melt isotropic amorphous films of unprecedendtly high kinetic stability. We are able to directly observe both bulk and front melting, and the crossover between them. We measure the front velocity over a broad range of conditions, and a crossover length scale that grows to nearly 400 particle diameters in the regime accessible to simulations. Our results disentangle the relative roles of kinetic stability and vapor deposition in the physical properties of stable glasses.