Heterogeneities in the glassy state

Anisotropic and heterogeneous dynamics in an aging colloidal gel

We investigate the out-of-equilibrium dynamics of a colloidal gel obtained by quenching a suspension of soft polymer-coated gold nanoparticles close to and below its gelation point using X-ray Photon Correlation Spectroscopy (XPCS). A faster relaxation process emergent from the localized motions of the nanoparticles reveals a dynamically-arrested network at the nanoscale as a key signature of the gelation process. We find that the slower network dynamics is hyperdiffusive with a compressed exponential form, consistent with stress-driven relaxation processes. Specifically, we use direction-dependent correlation functions to characterize the anisotropy in dynamics. We show that the anisotropy is greater for the gel close to its gelation point than at lower temperatures, and the anisotropy decreases as the gel ages. We quantify the anisotropic dynamical heterogeneities emergent in such a stress-driven dynamical system using higher order intensity correlations, and demonstrate that the aging phenomenon contributes significantly to the properties evaluated by the fluctuations in the intensity correlations. Our results provide important insights into the structural origin of the emergent anisotropic and cooperative heterogeneous dynamics, and we discuss analogies with previous work on other soft disordered systems.

What does the instantaneous normal mode spectrum tell us about dynamical heterogeneity in glass-forming fluids?

We examine the instantaneous normal mode spectrum of model metallic and polymeric glass-forming liquids. We focus on the localized modes in the unstable part of the spectrum [unstable localized (UL) modes] and find that the particles making the dominant contribution to the participation ratio form clusters that grow upon cooling in a fashion similar to the dynamical heterogeneity in glass-forming fluids, i.e., highly mobile (or immobile) particles form clusters that grow upon cooling; however, a comparison of the UL mode clusters to the mobile and immobile particle clusters indicates that they are distinct entities. We also show that the cluster size provides an alternate method to distinguish localized and delocalized modes, offering a significant practical advantage over the finite-size scaling approach. We examine the trajectories of particles contributing most to the UL modes and find that they have a slightly enhanced mobility compared to the average, and we determine a characteristic time quantifying the persistence time of this excess mobility. This time scale is proportional to the structural relaxation time τ_α of the fluid, consistent with a prediction by Zwanzig [Phys. Rev. 156, 190 (1967)] for the lifetime of collective excitations in cooled liquids. Evidently, these collective excitations serve to facilitate relaxation but do not actually participate in the motion associated with barrier crossing events governing activated transport. They also serve as a possible concrete realization of the "facilitation" clusters postulated in previous modeling of glass-forming liquids.

Effect of Liquid State Organization on Nanostructure and Strength of Model Multicomponent Solids

When a multicomponent liquid composed of particles with random interactions is slowly cooled below the freezing temperature, the fluid reorganizes in order to increase (decrease) the number of strong (weak) attractive interactions and solidifies into a structure composed of domains of strongly and of weakly interacting particles. Using Langevin dynamics simulations of a model system we find that the tensile strength, mode of fracture, and thermal stability of such solids differ from those of one-component solids and that these properties can be controlled by the method of preparation.

Continuous-time random-walk approach to supercooled liquids: Self-part of the van Hove function and related quantities

From equilibrium molecular dynamics (MD) simulations of a bead-spring model for short-chain glass-forming polymer melts we calculate several quantities characterizing the single-monomer dynamics near the (extrapolated) critical temperature [Formula: see text] of mode-coupling theory: the mean-square displacement g₀(t), the non-Gaussian parameter [Formula: see text] and the self-part of the van Hove function [Formula: see text] which measures the distribution of monomer displacements r in time t. We also determine these quantities from a continuous-time random walk (CTRW) approach. The CTRW is defined in terms of various probability distributions which we know from previous analysis. Utilizing these distributions the CTRW can be solved numerically and compared to the MD data with no adjustable parameter. The MD results reveal the heterogeneous and non-Gaussian single-particle dynamics of the supercooled melt near [Formula: see text]. In the time window of the early [Formula: see text] relaxation [Formula: see text] is large and [Formula: see text] is broad, reflecting the coexistence of monomer displacements that are much smaller ("slow particles") and much larger ("fast particles") than the average at time t, i.e. than [Formula: see text]. For large r the tail of [Formula: see text] is compatible with an exponential decay, as found for many glassy systems. The CTRW can reproduce the spatiotemporal dependence of [Formula: see text] at a qualitative to semiquantitative level. However, it is not quantitatively accurate in the studied temperature regime, although the agreement with the MD data improves upon cooling. In the early [Formula: see text] regime we also analyze the MD results for [Formula: see text] via the space-time factorization theorem predicted by ideal mode-coupling theory. While we find the factorization to be well satisfied for small r, both above and below [Formula: see text] , deviations occur for larger r comprising the tail of [Formula: see text]. The CTRW analysis suggests that single-particle "hops" are a contributing factor for these deviations.

Spatial correlations of elementary relaxation events in glass-forming liquids

The dynamical facilitation scenario, by which localized relaxation events promote nearby relaxation events in an avalanche process, has been suggested as the key mechanism connecting the microscopic and the macroscopic dynamics of structural glasses. Here we investigate the statistical features of this process via numerical simulations of a model structural glass. First we show that the relaxation dynamics of the system occurs through particle jumps that are irreversible, and that cannot be decomposed in smaller irreversible events. Then we show that each jump does actually trigger an avalanche. The characteristics of this avalanche change upon cooling, suggesting that the relaxation dynamics crossovers from a noise dominated regime, where jumps do not trigger other relaxation events, to a regime dominated by the facilitation process, where a jump triggers more relaxation events.

Continuous-time random-walk approach to supercooled liquids. II. Mean-square displacements in polymer melts

The continuous-time random walk (CTRW) describes the single-particle dynamics as a series of jumps separated by random waiting times. This description is applied to analyze trajectories from molecular dynamics (MD) simulations of a supercooled polymer melt. Based on the algorithm presented by Helfferich et al. [Phys. Rev. E 89, 042603 (2014)], we detect jump events of the monomers. As a function of temperature and chain length, we examine key distributions of the CTRW: the jump-length distribution (JLD), the waiting-time distribution (WTD), and the persistence-time distribution (PTD), i.e., the distribution of waiting times for the first jump. For the equilibrium (polymer) liquid under consideration, we verify that the PTD is determined by the WTD. For the mean-square displacement (MSD) of a monomer, the results for the CTRW model are compared with the underlying MD data. The MD data exhibit two regimes of subdiffusive behavior, one for the early α process and another at later times due to chain connectivity. By contrast, the analytical solution of the CTRW yields diffusive behavior for the MSD at all times. Empirically, we can account for the effect of chain connectivity in Monte Carlo simulations of the CTRW. The results of these simulations are then in good agreement with the MD data in the connectivity-dominated regime, but not in the early α regime where they systematically underestimate the MSD from the MD.

Continuous-time random-walk approach to supercooled liquids. I. Different definitions of particle jumps and their consequences

Single-particle trajectories in supercooled liquids display long periods of localization interrupted by "fast moves." This observation suggests a modeling by a continuous-time random walk (CTRW). We perform molecular dynamics simulations of equilibrated short-chain polymer melts near the critical temperature of mode-coupling theory Tc and extract "moves" from the monomer trajectories. We show that not all moves comply with the conditions of a CTRW. Strong forward-backward correlations are found in the supercooled state. A refinement procedure is suggested to exclude these moves from the analysis. We discuss the repercussions of the refinement on the jump-length and waiting-time distributions as well as on characteristic time scales, such as the average waiting time ("exchange time") and the average time for the first move ("persistence time"). The refinement modifies the temperature (T) dependence of these time scales. For instance, the average waiting time changes from an Arrhenius-type to a Vogel-Fulcher-type T dependence. We discuss this observation in the context of the bifurcation of the α process and (Johari) β process found in many glass-forming materials to occur near Tc. Our analysis lays the foundation for a study of the jump-length and waiting-time distributions, their temperature and chain-length dependencies, and the modeling of the monomer dynamics by a CTRW approach in the companion paper [J. Helfferich et al., Phys. Rev. E 89, 042604 (2014)].

The relationship of dynamical heterogeneity to the Adam-Gibbs and random first-order transition theories of glass formation

We carefully examine common measures of dynamical heterogeneity for a model polymer melt and test how these scales compare with those hypothesized by the Adam and Gibbs (AG) and random first-order transition (RFOT) theories of relaxation in glass-forming liquids. To this end, we first analyze clusters of highly mobile particles, the string-like collective motion of these mobile particles, and clusters of relative low mobility. We show that the time scale of the high-mobility clusters and strings is associated with a diffusive time scale, while the low-mobility particles' time scale relates to a structural relaxation time. The difference of the characteristic times for the high- and low-mobility particles naturally explains the well-known decoupling of diffusion and structural relaxation time scales. Despite the inherent difference of dynamics between high- and low-mobility particles, we find a high degree of similarity in the geometrical structure of these particle clusters. In particular, we show that the fractal dimensions of these clusters are consistent with those of swollen branched polymers or branched polymers with screened excluded-volume interactions, corresponding to lattice animals and percolation clusters, respectively. In contrast, the fractal dimension of the strings crosses over from that of self-avoiding walks for small strings, to simple random walks for longer, more strongly interacting, strings, corresponding to flexible polymers with screened excluded-volume interactions. We examine the appropriateness of identifying the size scales of either mobile particle clusters or strings with the size of cooperatively rearranging regions (CRR) in the AG and RFOT theories. We find that the string size appears to be the most consistent measure of CRR for both the AG and RFOT models. Identifying strings or clusters with the "mosaic" length of the RFOT model relaxes the conventional assumption that the "entropic droplets" are compact. We also confirm the validity of the entropy formulation of the AG theory, constraining the exponent values of the RFOT theory. This constraint, together with the analysis of size scales, enables us to estimate the characteristic exponents of RFOT.

Distribution of local relaxation events in an aging three-dimensional glass: spatiotemporal correlation and dynamical heterogeneity

We investigate the spatiotemporal distribution of microscopic relaxation events, defined through particle hops, in a model polymer glass using molecular dynamics simulations. We introduce an efficient algorithm to directly identify hops during the simulation, which allows the creation of a map of relaxation events for the whole system. Based on this map, we present density-density correlations between hops and directly extract correlation scales. These scales define collaboratively rearranging groups of particles and their size distributions are presented as a function of temperature and age. Dynamical heterogeneity is spatially resolved as the aggregation of hops into clusters, and we analyze their volume distribution and growth during aging. A direct comparison with the four-point dynamical susceptibility χ(4) reveals the formation of a single dominating cluster prior to the χ(4) peak, which indicates maximally correlated dynamics. An analysis of the fractal dimension of the hop clusters finds slightly noncompact shapes in excellent agreement with independent estimates from four-point correlations.

The physics of the colloidal glass transition

As one increases the concentration of a colloidal suspension, the system exhibits a dramatic increase in viscosity. Beyond a certain concentration, the system is said to be a colloidal glass; structurally, the system resembles a liquid, yet motions within the suspension are slow enough that it can be considered essentially frozen. For several decades, colloids have served as a valuable model system for understanding the glass transition in molecular systems. The spatial and temporal scales involved allow these systems to be studied by a wide variety of experimental techniques. The focus of this review is the current state of understanding of the colloidal glass transition, with an emphasis on experimental observations. A brief introduction is given to important experimental techniques used to study the glass transition in colloids. We describe features of colloidal systems near and in glassy states, including increases in viscosity and relaxation times, dynamical heterogeneity and ageing, among others. We also compare and contrast the glass transition in colloids to that in molecular liquids. Other glassy systems are briefly discussed, as well as recently developed synthesis techniques that will keep these systems rich with interesting physics for years to come.

Evidence for compact cooperatively rearranging regions in a supercooled liquid

We examine structural relaxation in a supercooled glass-forming liquid simulated by constant-energy constant-volume (NVE) molecular dynamics. Time correlations of the total kinetic energy fluctuations are used as a comprehensive measure of the system's approach to the ergodic equilibrium. We find that, under cooling, the total structural relaxation becomes delayed as compared with the decay of the component of the intermediate scattering function corresponding to the main peak of the structure factor. This observation can be explained by collective movements of particles preserving many-body structural correlations within compact three-dimensional (3D) cooperatively rearranging regions.

Low-frequency transient dynamic clusters in simulated amorphous Ni0.5Zr0.5 around the glass temperature

Molecular dynamics results are reported concerning cooperatively rearranging regions in simulated Ni0.5Zr0.5 melts down to 700 K . Emphasis is laid on discriminating between clusters of mobile atoms (CMA) from low-frequency dynamics and the all-frequency case, where the former characterize fluctuations and relaxations on the scale of the late beta regime and alpha decay, while the latter include, in addition, reversible high-frequency vibrations. Separation of the low-frequency part of the dynamics is carried out by low-pass filtering, exploiting the separation of time scales below the critical temperature T{c} of the mode-coupling theory. With increasing temperature, the low-frequency and all-frequency dynamics merge in the range of T{c} when the separation of time scales disappears. In the low-frequency CMA, the average size of correlated clusters of connected atoms turns out to be nearly one order of magnitude larger than in the all-frequency CMA. The low-frequency CMA appear as local clusters propagating extremely slowly in space with characteristic time scale of mus at 700 K , the scale of the onset of alpha decay.

Signatures of dynamical heterogeneity in the structure of glassy free-energy minima

From numerical minimization of a model free-energy functional for a system of hard spheres, we show that the width of the local peaks of the time-averaged density field at a glassy free-energy minimum exhibits large spatial variation, similar to that of the "local Debye-Waller factor" in simulations of dynamical heterogeneity. Molecular dynamics simulations starting from a particle configuration generated from the density distribution at a glassy free-energy minimum show similar spatial heterogeneity in the degree of localization, implying a direct connection between dynamical heterogeneity and the structure of glassy free-energy minima.

Simulations of aging and plastic deformation in polymer glasses

We study the effect of physical aging on the mechanical properties of a model polymer glass using molecular dynamics simulations. The creep compliance is determined simultaneously with the structural relaxation under a constant uniaxial load below yield at constant temperature. The model successfully captures universal features found experimentally in polymer glasses, including signatures of mechanical rejuvenation. We analyze microscopic relaxation time scales and show that they exhibit the same aging characteristics as the macroscopic creep compliance. In addition, our model indicates that the entire distribution of relaxation times scales identically with age. Despite large changes in mobility, we observe comparatively little structural change except for a weak logarithmic increase in the degree of short-range order that may be correlated with an observed decrease in aging with increasing load.

Avalanche size scaling in sheared three-dimensional amorphous solid

We study the statistics of plastic rearrangement events in a simulated amorphous solid at T=0. Events are characterized by the energy release and the "slip volume", the product of plastic strain and system volume. Their distributions for a given system size L appear to be exponential, but a characteristic event size cannot be inferred, because the mean values of these quantities increase as Lalpha with alpha approximately 3/2. In contrast with results obtained in 2D models, we do not see simply connected avalanches. The exponent suggests a fractal shape of the avalanches, which is also evidenced by the mean fractal dimension and participation ratio.

Aging in attraction-driven colloidal glasses

Aging in an attraction-driven colloidal glass is studied by computer simulations. The system is equilibrated without attraction and instantaneously "quenched," at constant colloid volume fraction, to one of two states beyond the glass transition; one is close to the transition, and the other one deep in the glass. The evolution of structural properties shows that bonds form in the system, increasing the local density, creating density deficits (holes) elsewhere. This process slows down with the time elapsed since the quench. As a consequence of bond formation, there is a slowing down of the dynamics, as measured by the mean-squared displacement and the density, bond, and environment correlation functions. The density correlations can be time rescaled to collapse their long time (structural) decay. The time scale for structural relaxation shows for both quenches a superlinear dependence on waiting time; it grows faster than the bond lifetime, showing the collective origin of the transition. At long waiting times and high attraction strength, we observe arrested dynamics for more than three decades in time, although individual bonds are not permanent on this time scale. The localization length decreases as the state moves deeper in the glass; the nonergodicity parameter oscillates in phase with the structure factor. Our main results are obtained for systems with a barrier in the pair potential that inhibits phase separation. However, when this barrier is removed for the case of a deep quench, we find changes in the static structure but almost none in the dynamics. Hence our results for the aging behavior remain relevant to experiments in which the glass transition competes with phase separation.

Influence of fictive temperature and composition of silica glass on anomalous elastic behaviour

In order to determine the influence of the thermal history (fictive temperature) and OH content on the elastic properties of silica glass, we have investigated high resolution in situ Brillouin experiments on SiO(2) glass from room temperature to the supercooled liquid at 1773 K across the glass transition. The well known anomalous increase of elastic modulus in the glassy state and in the supercooled liquid regime is observed. No change in the slope of the elastic moduli of silica appears as a characteristic of the glass transition, in contrast to what happens in various other glasses. We show that thermal history has a weak effect on elastic moduli in the glass transition regime for silica glass. The effect of the water content in silica glass is greater than the fictive temperature effect and gives larger changes in the amplitude of the elastic modulus for the same thermal dependence. A singular decrease above 1223 K is also observed in the shear moduli for hydrated samples. Different models explaining the temperature dependence of the elastic properties in relationship with frozen-in density fluctuations or with the structure are discussed.