Fluctuating mobility generation and transport in glasses

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.

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.

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.

Dynamical theory of shear bands in structural glasses

The heterogeneous elastoplastic deformation of structural glasses is explored using the framework of the random first-order transition theory of the glass transition along with an extended mode-coupling theory that includes activated events. The theory involves coupling the continuum elastic theory of strain transport with mobility generation and transport as described in the theory of glass aging and rejuvenation. Fluctuations that arise from the generation and transport of mobility, fictive temperature, and stress are treated explicitly. We examine the nonlinear flow of a glass under deformation at finite strain rate. The interplay among the fluctuating fields leads to the spatially heterogeneous dislocation of the particles in the glass, i.e., the appearance of shear bands of the type observed in metallic glasses deforming under mechanical stress.

The role of thermodynamic stability in the characteristics of the devitrification front of vapour-deposited glasses of toluene

Glass stability and molecular shape affect the transformation mechanism of vapour deposited glasses.

The melting of stable glasses is governed by nucleation-and-growth dynamics

We discuss the microscopic mechanisms by which low-temperature amorphous states, such as ultrastable glasses, transform into equilibrium fluids, after a sudden temperature increase. Experiments suggest that this process is similar to the melting of crystals, thus differing from the behaviour found in ordinary glasses. We rationalize these observations using the physical idea that the transformation process takes place close to a "hidden" equilibrium first-order phase transition, which is observed in systems of coupled replicas. We illustrate our views using simulation results for a simple two-dimensional plaquette spin model, which is known to exhibit a range of glassy behaviour. Our results suggest that nucleation-and-growth dynamics, as found near ordinary first-order transitions, is also the correct theoretical framework to analyse the melting of ultrastable glasses. Our approach provides a unified understanding of multiple experimental observations, such as propagating melting fronts, large kinetic stability ratios, and "giant" dynamic length scales. We also provide a comprehensive discussion of available theoretical pictures proposed in the context of ultrastable glass melting.

Vapor-deposited glasses of methyl-m-toluate: How uniform is stable glass transformation?

AC chip nanocalorimetry is used to characterize vapor-deposited glasses of methyl-m-toluate (MMT). Physical vapor deposition can prepare MMT glasses that have lower heat capacity and significantly higher kinetic stability compared to liquid-cooled glasses. When heated, highly stable MMT glasses transform into the supercooled liquid via propagating fronts. We present the first quantitative analysis of the temporal and spatial uniformities of these transformation fronts. The front velocity varies by less than 4% over the duration of the transformation. For films 280 nm thick, the transformation rates at different spatial positions in the film differ by about 25%; this quantity may be related to spatially heterogeneous dynamics in the stable glass. Our characterization of the kinetic stability of MMT stable glasses extends previous dielectric experiments and is in excellent agreement with these results.

Orientational anisotropy in simulated vapor-deposited molecular glasses

Enhanced kinetic stability of vapor-deposited glasses has been established for a variety of glass organic formers. Several recent reports indicate that vapor-deposited glasses can be orientationally anisotropic. In this work, we present results of extensive molecular simulations that mimic a number of features of the experimental vapor deposition process. The simulations are performed on a generic coarse-grained model and an all-atom representation of N,N'-bis(3-methylphenyl)-N,N'-diphenylbenzidine (TPD), a small organic molecule whose vapor-deposited glasses exhibit considerable orientational anisotropy. The coarse-grained model adopted here is found to reproduce several key aspects reported in experiments. In particular, the molecular orientation of vapor-deposited glasses is observed to depend on substrate temperature during deposition. For a fixed deposition rate, the molecular orientation in the glasses changes from isotropic, at the glass transition temperature, Tg, to slightly normal to the substrate at temperatures just below Tg. Well below Tg, molecular orientation becomes predominantly parallel to the substrate. The all-atom model is used to confirm some of the equilibrium structural features of TPD interfaces that arise above the glass transition temperature. We discuss a mechanism based on distinct orientations observed at equilibrium near the surface of the film, which get trapped within the film during the non-equilibrium process of vapor deposition.

Cooling-rate dependence of kinetic and mechanical stabilities of simulated glasses

Glasses created through vapor deposition on a substrate maintained at a proper temperature possess higher kinetic and mechanical stabilities than glasses created by cooling at a constant rate. Molecular dynamics simulations are being increasingly used to understand why vapor deposition improves glasses' stability. There are, however, few detailed molecular dynamics studies of the dependence of the properties of glasses cooled at a constant rate on the rate of cooling. Thus, there is no clear benchmark for comparing ultrastable simulated glasses to simulated glasses prepared through cooling at a constant rate. Here, we examine the dependence of the properties of simulated glasses on the cooling rate used in their preparation. We examine the kinetic stability by measuring the time it takes for a glass to transform back to a liquid upon heating and heterogeneous dynamics during heating. We also examine properties of the energy landscape, and we evaluate mechanical stability by calculating the shear modulus of the glass. The methods outlined here can be used to assess kinetic and mechanical stabilities of simulated glasses generated using specialized algorithms and provide a benchmark for those algorithms.

Thermal stability of vapor-deposited stable glasses of an organic semiconductor

Vapor-deposited organic glasses can show enhanced kinetic stability relative to liquid-cooled glasses. When such stable glasses of model glassformers are annealed above the glass transition temperature Tg, they lose their thermal stability and transform into the supercooled liquid via constant velocity propagating fronts. In this work, we show that vapor-deposited glasses of an organic semiconductor, N,N'-bis(3-methylphenyl)-N,N'-diphenylbenzidine (TPD), also transform via propagating fronts. Using spectroscopic ellipsometry and a new high-throughput annealing protocol, we measure transformation front velocities for TPD glasses prepared with substrate temperatures (TSubstrate) from 0.63 to 0.96 Tg, at many different annealing temperatures. We observe that the front velocity varies by over an order of magnitude with TSubstrate, while the activation energy remains constant. Using dielectric spectroscopy, we measure the structural relaxation time of supercooled TPD. We find that the mobility of the liquid and the structure of the glass are independent factors in controlling the thermal stability of TPD films. In comparison to model glassformers, the transformation fronts of TPD have similar velocities and a similar dependence on TSubstrate, suggesting universal behavior. These results may aid in designing active layers in organic electronic devices with improved thermal stability.

Transformation kinetics of vapor-deposited thin film organic glasses: the role of stability and molecular packing anisotropy

The growth front velocity of indomethacin glasses depends on deposition conditions but is not unambigously determined by its thermodynamic stability when the structure is not completely isotropic.

On the mechanism of activated transport in glassy liquids

We explore several potential issues that have been raised over the years regarding the "entropic droplet" scenario of activated transport in liquids, due to Wolynes and co-workers, with the aim of clarifying the status of various approximations of the random first-order transition theory (RFOT) of the structural glass transition. In doing so, we estimate the mismatch penalty between alternative aperiodic structures, above the glass transition; the penalty is equal to the typical magnitude of free energy fluctuations in the liquid. The resulting expressions for the activation barrier and the cooperativity length contain exclusively bulk, static properties; in their simplest form they contain only the bulk modulus and the configurational entropy per unit volume. The expressions are universal in that they do not depend explicitly on the molecular detail. The predicted values for the barrier and cooperativity length and, in particular, the temperature dependence of the barrier are in satisfactory agreement with observation. We thus confirm that the entropic droplet picture is indeed not only internally consistent but is also fully constructive, consistent with the apparent success of its many quantitative predictions. A simple view of a glassy liquid as a locally metastable, degenerate pattern of frozen-in stress emerges in the present description. Finally, we derive testable relationships between the bulk modulus and several characteristics of glassy liquids and peculiarities in low-temperature glasses.

Role of fragility in the formation of highly stable organic glasses

In situ dielectric spectroscopy has been used to characterize vapor-deposited glasses of methyl-m-toluate (MMT), an organic glass former with low fragility (m = 60). Deposition near 0.84T(g) produces glasses of very high kinetic stability; these materials are comparable in stability to the most stable glasses produced from more fragile glass formers. Highly stable glasses of MMT, when annealed above T(g), transform into the supercooled liquid by a heterogeneous mechanism. A constant velocity propagating front is initiated at the free surface and controls the transformation of thin films. The transition to a bulk-dominated transformation process occurs at 5 μm, the largest length scale reported for any glass. Contrary to recent conclusions, we find that physical vapor deposition can form highly stable organic glasses across the entire range of liquid fragilities.

Dynamical heterogeneity of the glassy state

We compare dynamical heterogeneities in equilibrated supercooled liquids and in the nonequilibrium glassy state within the framework of the random first order transition theory. Fluctuating mobility generation and transport in the glass are treated by numerically solving stochastic continuum equations for mobility and fictive temperature fields that arise from an extended mode coupling theory containing activated events. Fluctuating spatiotemporal structures in aging and rejuvenating glasses lead to dynamical heterogeneity in glasses with characteristics distinct from those found in the equilibrium supercooled liquid. The non-Gaussian distribution of activation free energies, the stretching exponent β, and the growth of characteristic lengths are studied along with the four-point dynamical correlation function. Asymmetric thermodynamic responses upon heating and cooling are predicted to be the result of the heterogeneity and the out-of-equilibrium behavior of glasses below Tg. Our numerical results agree with experimental calorimetry. We numerically confirm the prediction of Lubchenko and Wolynes in the glass that the dynamical heterogeneity can lead to noticeably bimodal distributions of local fictive temperatures during some histories of preparation which explains in a unified way recent experimental observations that have been interpreted as coming from there being two distinct equilibration mechanisms in glasses.

Microscopically based calculations of the free energy barrier and dynamic length scale in supercooled liquids: the comparative role of configurational entropy and elasticity

We compute the temperature-dependent barrier for α-relaxations in several liquids, without adjustable parameters, using experimentally determined elastic, structural, and calorimetric data. We employ the random first order transition (RFOT) theory, in which relaxation occurs via activated reconfigurations between distinct, aperiodic minima of the free energy. Two different approximations for the mismatch penalty between the distinct aperiodic states are compared, one due to Xia and Wolynes (Proc. Natl. Acad. Sci. U. S. A. 2000, 97, 2990), which scales universally with temperature as for hard spheres, and one due to Rabochiy and Lubchenko (J. Chem. Phys. 2013, 138, 12A534), which employs measured elastic and structural data for individual substances. The agreement between the predictions and experiment is satisfactory, given the uncertainty in the measured experimental inputs. The explicitly computed barriers are used to calculate the glass transition temperature for each substance--a kinetic quantity--from the static input data alone. The temperature dependence of both the elastic and structural constants enters the temperature dependence of the barrier over an extended range to a degree that varies from substance to substance. The lowering of the configurational entropy, however, seems to be the dominant contributor to the barrier increase near the laboratory glass transition, consistent with previous experimental tests of the RFOT theory using the XW approximation. In addition, we compute the temperature dependence of the dynamical correlation length, also without using adjustable parameters. These agree well with experimental estimates obtained using the Berthier et al. (Science 2005, 310, 1797) procedure. Finally, we find the temperature dependence of the complexity of a rearranging region is consistent with the picture based on the RFOT theory but is in conflict with the assumptions of the Adam-Gibbs and "shoving" scenarios for the viscous slowing down in supercooled liquids.