Out-of-equilibrium dynamical fluctuations in glassy systems

Microscopic dynamics underlying the stress relaxation of arrested soft materials

Arrested soft materials such as gels and glasses exhibit a slow stress relaxation with a broad distribution of relaxation times in response to linear mechanical perturbations. Although this macroscopic stress relaxation is an essential feature in the application of arrested systems as structural materials, consumer products, foods, and biological materials, the microscopic origins of this relaxation remain poorly understood. Here, we elucidate the microscopic dynamics underlying the stress relaxation of such arrested soft materials under both quiescent and mechanically perturbed conditions through X-ray photon correlation spectroscopy. By studying the dynamics of a model associative gel system that undergoes dynamical arrest in the absence of aging effects, we show that the mean stress relaxation time measured from linear rheometry is directly correlated to the quiescent superdiffusive dynamics of the microscopic clusters, which are governed by a buildup of internal stresses during arrest. We also show that perturbing the system via small mechanical deformations can result in large intermittent fluctuations in the form of avalanches, which give rise to a broad non-Gaussian spectrum of relaxation modes at short times that is observed in stress relaxation measurements. These findings suggest that the linear viscoelastic stress relaxation in arrested soft materials may be governed by nonlinear phenomena involving an interplay of internal stress relaxations and perturbation-induced intermittent avalanches.

Analytic form of a two-dimensional critical distribution

This paper explores the possibility of establishing an analytic form of the distribution of the order parameter fluctuations in a two-dimensional critical spin-wave model, or width fluctuations of a two-dimensional Edwards-Wilkinson interface. It is shown that the characteristic function of the distribution can be expressed exactly as a gamma function quotient, while a Charlier series, using the convolution of two Gumbel distributions as the kernel, converges to the exact result over a restricted domain. These results can also be extended to calculate the temperature dependence of the distribution and give an insight into the origin of Gumbel-like distributions in steady-state and equilibrium quantities that are not extreme values.

Universal scaling in the aging of the strong glass former SiO2

We show that the aging dynamics of a strong glass former displays a strikingly simple scaling behavior, connecting the average dynamics with its fluctuations, namely, the dynamical heterogeneities. We perform molecular dynamics simulations of SiO2 with van Beest-Kramer-van Santen interactions, quenching the system from high to low temperature, and study the evolution of the system as a function of the waiting time tw measured from the instant of the quench. We find that both the aging behavior of the dynamic susceptibility χ4 and the aging behavior of the probability distribution P(fs,r) of the local incoherent intermediate scattering function fs,r can be described by simple scaling forms in terms of the global incoherent intermediate scattering function C. The scaling forms are the same that have been found to describe the aging of several fragile glass formers and that, in the case of P(fs,r), have been also predicted theoretically. A thorough study of the length scales involved highlights the importance of intermediate length scales. We also analyze directly the scaling dependence on particle type and on wavevector q and find that both the average and the fluctuations of the slow aging dynamics are controlled by a unique aging clock, which is not only independent of the wavevector q, but is also the same for O and Si atoms.

Fluctuations in the time variable and dynamical heterogeneity in glass-forming systems

We test a hypothesis for the origin of dynamical heterogeneity in slowly relaxing systems, namely that it emerges from soft (Goldstone) modes associated with a broken continuous symmetry under time reparametrizations. We do this by constructing coarse grained observables and decomposing the fluctuations of these observables into transverse components, which are associated with the postulated time-fluctuation soft modes, and a longitudinal component, which represents the rest of the fluctuations. Our test is performed on data obtained in simulations of four models of structural glasses. As the hypothesis predicts, we find that the time reparametrization fluctuations become increasingly dominant as temperature is lowered and timescales are increased. More specifically, the ratio between the strengths of the transverse fluctuations and the longitudinal fluctuations grows as a function of the dynamical susceptibility, χ(4), which represents the strength of the dynamical heterogeneity; and the correlation volumes for the transverse fluctuations are approximately proportional to those for the dynamical heterogeneity, while the correlation volumes for the longitudinal fluctuations remain small and approximately constant.

Dynamical heterogeneities as fingerprints of a backbone structure in Potts models

We investigate slow nonequilibrium dynamical processes in a two-dimensional q-state Potts model with both ferromagnetic and ±J couplings. Dynamical properties are characterized by means of the mean-flipping time distribution. This quantity is known for clearly unveiling dynamical heterogeneities. Using a two-times protocol we characterize the different time scales observed and relate them to growth processes occurring in the system. In particular we target the possible relation between the different time scales and the spatial heterogeneities originated in the ground-state topology, which are associated to the presence of a backbone structure. We perform numerical simulations using an approach based on graphis processing units (GPUs) which permits us to reach large system sizes. We present evidence supporting both the idea of a growing process in the preasymptotic regime of the glassy phases and the existence of a backbone structure behind this process.

Fluctuation-dissipation relations and field-free algorithms for the computation of response functions

We discuss the relation between the fluctuation-dissipation relation derived by Chatelain and Ricci-Tersenghi [C. Chatelain, J. Phys. A 36, 10739 (2003); F. Ricci-Tersenghi, Phys. Rev. E 68, 065104(R) (2003)] and that by Lippiello-Corberi-Zannetti [E. Lippiello, F. Corberi, and M. Zannetti, Phys. Rev. E 71, 036104 (2005)]. In order to do that, we rederive the fluctuation-dissipation relation for systems of discrete variables evolving in discrete time via a stochastic nonequilibrium Markov process. The calculation is carried out in a general formalism comprising the Chatelain, Ricci-Tersenghi, result and that by Lippiello-Corberi-Zannetti as special cases. The applicability, generality, and experimental feasibility of the two approaches are thoroughly discussed. Extending the analytical calculation to the variance of the response function, we show the advantage of field-free numerical methods with respect to the standard method, where the perturbation is applied. We also show that the signal-to-noise ratio is better (by a factor square root of 2) in the algorithm of Lippiello-Corberi-Zannetti with respect to that of Chatelain-Ricci Tersenghi.

Fluctuation relation and heterogeneous superdiffusion in glassy transport

Current fluctuations and related steady-state fluctuation relation are investigated in simple coarse-grained lattice-gas analogs of a non-Newtonian fluid driven by a constant and uniform force field in two regimes of small entropy production. Non-Gaussian current fluctuations and deviations from fluctuation relation are observed and related to the existence of growing amorphous correlations and heterogeneous anomalous diffusion regimes.

Equilibrium and nonequilibrium fluctuations in a glass-forming liquid

Glass-forming liquids display strong fluctuations-dynamical heterogeneities-near their glass transition. By numerically simulating a binary Weeks-Chandler-Andersen liquid and varying both temperature and time scale, we investigate the probability distributions of two kinds of local fluctuations in the nonequilibrium (aging) regime and in the equilibrium regime, and find them to be very similar in the two regimes and across temperatures. We also observe that, when appropriately rescaled, the integrated dynamic susceptibility is very weakly dependent on temperature and very similar in both regimes.

Spatiotemporal structures in aging and rejuvenating glasses

Complex spatiotemporal structures develop during the process of aging glasses after cooling and of rejuvenating glasses on heating. The key to understanding these structures is the interplay between the activated reconfiguration events that generate mobility and the transport of mobility. These effects are both accounted for by combining the random first-order transition theory of activated events with mode coupling theory in an inhomogeneous setting. The predicted modifications by mobility transport of the time course of the aging regime are modest. In contrast, the rejuvenation process is strongly affected through the propagation of fronts of enhanced mobility originating from the initial reconfiguration events. The structures in a rejuvenating glass resemble flames. An analysis along the lines of combustion theory provides an estimate of the front propagation speed. Heterogeneous rejuvenation naturally should occur for glasses with free surfaces. The analogy with combustion also provides a way of looking at the uptake of diluents by glasses described by case II and super case II diffusion.

Growth of spatial correlations in the aging of a simple structural glass

We present a detailed numerical study of dynamical heterogeneities in the aging regime of a simple binary Lennard-Jones glass former. For most waiting times t_{w} and final times t , both the dynamical susceptibility chi_{4}(t,t_{w}) and the dynamical correlation length xi_{4}(t,t_{w}) can be approximated as products of two factors: (i) a waiting-time-dependent scale that grows as a power of t_{w} , and (ii) a scaling function dependent on t,t_{w} only through the value of the intermediate scattering function C(t,t_{w}) . We find that chi_{4}(t,t_{w}) is determined only in part by the correlation volume.

Connecting microscopic simulations with kinetically constrained models of glasses

Kinetically constrained spin models are known to exhibit dynamical behavior mimicking that of glass forming systems. They are often understood as coarse-grained models of glass formers, in terms of some "mobility" field. The identity of this "mobility" field has remained elusive due to the lack of coarse-graining procedures to obtain these models from a more microscopic point of view. Here we exhibit a scheme to map the dynamics of a two-dimensional soft disk glass former onto a kinetically constrained spin model, providing an attempt at bridging these two approaches.

Strong dynamical heterogeneities in the violation of the fluctuation-dissipation theorem in spin glasses

We analyze numerically the violation of the fluctuation-dissipation theorem (FDT) in the +/-J Edwards-Anderson (EA) spin-glass model. Using single spin probability densities we reveal the presence of strong dynamical heterogeneities, which correlate with ground-state information. The physical interpretation of the results shows that the spins can be divided into two sets. In 3D, one set forms a compact structure which presents a coarseninglike behavior with its characteristic violation of the FDT, while the other asymptotically follows the FDT. Finally, we compare the dynamical behavior observed in 3D with 2D.

Space-time thermodynamics and subsystem observables in a kinetically constrained model of glassy materials

In a recent article [M. Merolle et al., Proc. Natl. Acad. Sci. U.S.A. 102, 10837 (2005)], it was argued that dynamic heterogeneity in d-dimensional glass formers is a manifestation of an order-disorder phenomenon in the d+1 dimensions of space time. By considering a dynamical analog of the free energy, evidence was found for phase coexistence between active and inactive regions of space time, and it was suggested that this phenomenon underlies the glass transition. Here we develop these ideas further by investigating in detail the one-dimensional Fredrickson-Andersen (FA) model, in which the active and inactive phases originate in the reducibility of the dynamics. We illustrate the phase coexistence by considering the distributions of mesoscopic space-time observables. We show how the analogy with phase coexistence can be strengthened by breaking microscopic reversibility in the FA model, leading to a nonequilibrium theory in the directed percolation universality class.

Aging and intermittency in a p-spin model of a glass

We numerically analyze the statistics of the heat flow between an aging system and its thermal bath, following a method proposed and tested for a spin-glass model in a recent paper [P. Sibani and H. J. Jensen, Europhys. Lett. 69, 563 (2005)]. The present system, which lacks quenched randomness, consists of Ising spins located on a cubic lattice, with each plaquette contributing to the total energy the product of the four spins located at its corners. Similarly to our previous findings, energy leaves the system in rare but large, so-called intermittent, bursts which are embedded in reversible and equilibriumlike fluctuations of zero average. The intermittent bursts, or quakes, dissipate the excess energy trapped in the initial state at a rate which falls off with the inverse of the age. This strongly heterogeneous dynamical picture is explained using the idea that quakes are triggered by energy fluctuations of record size, which occur independently within a number of thermalized domains. From the temperature dependence of the width of the reversible heat fluctuations we surmise that these domains have an exponential density of states. Finally, we show that the heat flow consists of a temperature independent term and a term with an Arrhenius temperature dependence. Microscopic dynamical and structural information can thus be extracted from numerical intermittency data. This type of analysis seems now within the reach of time resolved microcalorimetry techniques.

Signature of the ground-state topology in the low-temperature dynamics of spin glasses

We numerically address the issue of how the ground-state topology is reflected in the finite temperature dynamics of the +/-J Edwards-Anderson spin glass model. In this system a careful study of the ground-state configurations allows us to classify spins into two sets: solidary and nonsolidary spins. We show that these sets quantitatively account for the dynamical heterogeneities found in the mean flipping time distribution at finite low temperatures. The results highlight the relevance of taking into account the ground-state topology in the analysis of the finite temperature dynamics of spin glasses.

Time-resolved-correlation measurements of temporally heterogeneous dynamics

Time resolved correlation (TRC) is a recently introduced light scattering technique that allows one to detect and quantify dynamic heterogeneities. The technique is based on the analysis of the temporal evolution of the speckle pattern generated by the light scattered by a sample, which is quantified by cI(t, tau), the degree of correlation between speckle images recorded at time t and t + tau. Heterogeneous dynamics results in significant fluctuations of cI(t,tau) with time t. We describe how to optimize TRC measurements and how to detect and avoid possible artifacts. The statistical properties of the fluctuations of cI are analyzed by studying their variance, probability distribution function, and time autocorrelation function. We show that these quantities are affected by a noise contribution due to the finite number N of detected speckles. We propose and demonstrate a method to correct for the noise contribution, based on a N--> infinity extrapolation scheme. Examples from both homogeneous and heterogeneous dynamics are provided. Connections with recent numerical and analytical works on heterogeneous glassy dynamics are briefly discussed.

Space-time thermodynamics of the glass transition

We consider the probability distribution for fluctuations in dynamical action and similar quantities related to dynamic heterogeneity. We argue that the so-called "glass transition" is a manifestation of low action tails in these distributions where the entropy of trajectory space is subextensive in time. These low action tails are a consequence of dynamic heterogeneity and an indication of phase coexistence in trajectory space. The glass transition, where the system falls out of equilibrium, is then an order-disorder phenomenon in space-time occurring at a temperature T(g), which is a weak function of measurement time. We illustrate our perspective ideas with facilitated lattice models and note how these ideas apply more generally.

Dynamical susceptibility of glass formers: contrasting the predictions of theoretical scenarios

We compute analytically and numerically the four-point correlation function that characterizes nontrivial cooperative dynamics in glassy systems within several models of glasses: elastoplastic deformations, mode-coupling theory (MCT), collectively rearranging regions (CRR's), diffusing defects, and kinetically constrained models (KCM's). Some features of the four-point susceptibility chi(4) (t) are expected to be universal: at short times we expect a power-law increase in time as t(4) due to ballistic motion (t(2) if the dynamics is Brownian) followed by an elastic regime (most relevant deep in the glass phase) characterized by a t or sqrt[t] growth, depending on whether phonons are propagative or diffusive. We find in both the beta and early alpha regime that chi(4) approximately t(mu), where mu is directly related to the mechanism responsible for relaxation. This regime ends when a maximum of chi(4) is reached at a time t= t(*) of the order of the relaxation time of the system. This maximum is followed by a fast decay to zero at large times. The height of the maximum also follows a power law chi(4) (t(*)) approximately t(*lambda). The value of the exponents mu and lambda allows one to distinguish between different mechanisms. For example, freely diffusing defects in d=3 lead to mu=2 and lambda=1 , whereas the CRR scenario rather predicts either mu=1 or a logarithmic behavior depending on the nature of the nucleation events and a logarithmic behavior of chi(4) (t(*)) . MCT leads to mu=b and lambda=1/gamma , where b and gamma are the standard MCT exponents. We compare our theoretical results with numerical simulations on a Lennard-Jones and a soft-sphere system. Within the limited time scales accessible to numerical simulations, we find that the exponent mu is rather small, mu<1 , with a value in reasonable agreement with the MCT predictions, but not with the prediction of simple diffusive defect models, KCM's with noncooperative defects, and CRR's. Experimental and numerical determination of chi(4) (t) for longer time scales and lower temperatures would yield highly valuable information on the glass formation mechanism.