# Atomistic mechanism of physical ageing in glassy materials

###### Abstract

Using molecular simulations, we identify microscopic relaxation events of individual particles in ageing structural glasses, and determine the full distribution of relaxation times. We find that the memory of the waiting time elapsed since the quench extends only up to the first relaxation event, while the distribution of all subsequent relaxation times (persistence times) follows a power law completely independent of history. Our results are in remarkable agreement with the well known phenomenological trap model of ageing. A continuous time random walk (CTRW) parametrized with the atomistic distributions captures the entire bulk diffusion behavior and explains the apparent scaling of the relaxation dynamics with during ageing, as well as observed deviations from perfect scaling.

###### pacs:

81.05.Kf###### pacs:

64.70.Q-###### pacs:

05.40.-aGlasses (including metallic glasses) Theory and modeling of the glass transition Fluctuation phenomena, random processes, noise, and Brownian motion

The dynamics of glass forming materials is characterized by slow, spatially heterogeneous and temporally intermittent processes [1, 2, 3, 4, 5, 6]. Long periods of caged particle vibrations are interrupted by rapid, localized structural relaxations. Dynamical heterogeneity near the glass transition leads to stretched exponential relaxation of correlation and response functions, as well as subdiffusive, non-Gaussian transport. In the glassy state, configurational degrees of freedom continue to evolve slowly. This phenomenon is called physical ageing [7] and is a defining and universal feature of nonequilibrium glassy dynamics [8]. Structural relaxations become increasingly sluggish with the waiting time elapsed since vitrification, and almost all material properties change as a result.

Ageing often manifests itself as a simple rescaling of the slow, configurational part of the global correlation functions , where is the ageing exponent. This scaling behaviour has been observed in a wide range of disordered materials [1] and is particularly pronounced in polymer glasses [7], colloidal suspensions [9, 10], and physical gels [11]. Recent simulations of structural glasses show that scaling also applies approximately to the distribution of local correlation functions [12] and the distribution of particle displacements (van Hove function) [13], which superimpose when plotted at times of equal mean. Scaling of local correlations can be proved to hold exactly for certain spin glass models whose dynamical action obeys time reparametrization invariance [14]. Despite the apparently universal nature of this phenomenon [15], the physical mechanism behind it remains elusive.

In this Letter, we directly observe the atomistic dynamics underlying physical ageing on a single particle level through molecular dynamics simulations of two model systems for atomic and polymeric glasses. We consider both the well-known amorphous 80/20 binary Lennard-Jones mixture [16], as well as a bead-spring model for polymer glasses [17, 18], where particles interact via the Lennard-Jones potential and are coupled together by stiff springs to form chains of length 10. Twenty thousand particles each are simulated in a periodic simulation box. For the polymer (80/20) model, the initial ensemble is equilibrated at (4.0), and then quenched rapidly to (0.33) to form a glass. The system is aged at constant pressure for various waiting times and then the dynamics are followed as a function of the measurement time . All results are quoted in Lennard-Jones units.

In both models, the effect of ageing can be immediately observed through the one-dimensional mean squared displacement shown in Fig. 1, which exhibits three regimes: a short time vibrational regime precedes a long plateau characterized by “caged” motion, followed by a cage escape or -relaxation regime. The cage escape time increases with increasing , but -curves for different waiting times can be superimposed in the cage escape regime by rescaling time with . At long times , scaling of the mean squared displacement fails for both models, and curves for the different ages merge. Experimental creep compliance curves show the same behavior, which is a well-known consequence of subageing () [7].

We proceed to analyze the intermittent dynamics in these systems by decomposing single particle trajectories into discrete hopping events. To investigate the hopping dynamics, a subset of five thousand particle trajectories are recorded. A running average of the atom position and its standard deviation is found for each atom over a small time window of 400 (40 individual snapshots of the particle positions) [19]. A hop is identified when the standard deviation is greater than a threshold of [20]. This method differs from the usual technique of identifying hops via a threshold in the hop distance [5, 21, 22], and effectively captures hops at all length scales through their increased activity during the relaxation event. While the correlation factor between hops is nearly zero, approximately of hops are either forward or backward correlated. Removing these hops does not affect the results.

Fig. 2 shows a typical atomic trajectory with the hops identified using this criterion. The initial hop times , the times between all subsequent hops , also called the persistence time, and the displacements were binned to form distribution functions , , and . The initial hop time distribution is shown in Fig. 3(a) for the Lennard-Jones mixture and in (d) for the polymer. For all waiting times, it appears to take the form of two power laws. The exponents are waiting time independent, but the crossover time between the two regimes shifts towards longer times with increasing .

Figures 3(b) and (e) show that the distributions also decay as a power law: (Lennard-Jones mixture) and (polymer). By contrast, is independent of and shows no ageing effects. The one-dimensional hop displacement distribution (Fig. 3(c) and (f)) is sharply peaked at zero displacement, and is also independent of waiting time. This is perhaps not surprising, as it has been shown that the density and the dynamical correlation lengths vary little during ageing [6]. has a purely exponential form for the mixture, whereas in the polymer model it decays more rapidly in the tails due to chain connectivity effects. Both the first hop time distribution, and the displacement distribution are in qualitative agreement with results obtained by analyzing hops between adjacent inherent structures [23, 24, 25].

The measured distributions have much in common with the well-known phenomenological trap model of ageing [26]. In the trap model, as in our results, is constant and continuously evolves towards longer times in the ageing regime. In fact, even the form of our measured distributions is similar. The distributions and decay as weak power laws at long times with an exponent , where is the noise temperature which governs activated hopping between states [27]. In this formulation, is either a constant (exponential density of states), or more realistically [23, 24, 25], a slowly varying function of time (Gaussian density of states) [26]. We do not observe any curvature in , but cannot rule out its existence at longer timescales. While our results clearly show subageing behavior, the ageing exponent in the trap model is strictly unity [26]. However, this holds true only in the limit . We have observed that depending on the quench conditions, there can be a long lived transient regime where correlation functions approximately scale with an ageing exponent .

To confirm that our analysis captures all important physics, we recompute the ageing dynamics through stochastic realizations of a continuous time random walk (CTRW) using the measured hop distributions , and as input. For the polymer model, a pure random walk is insufficient as the chain connectivity induces subdiffusive Rouse dynamics. The mean squared displacement of the atoms as a function of the number of hops is found to be ( for the binary mixture). In the CTRW, the Rouse behaviour is accounted for by choosing a new particle position from a Gaussian distribution widening as at every step. This procedure creates the proper displacement statistics. In both systems, the initial displacement is chosen such that the mean squared displacement coincides with the MD data at the shortest measured relaxation time. This is slightly higher than the actual vibrational amplitude, but accounts for relaxations which occur at times down to the vibrational timescale, but could not be captured otherwise.

The full distribution of displacements (van Hove function)

(1) |

as well the as the mean squared displacements are compared with MD data in Fig. 4 for the polymer model. As time progresses, a Gaussian caged peak at small displacements loses particles to widening exponential tails. The origin of the exponential tails has recently generated intense interest, as they appear to be a universal feature of structural glasses [28, 22, 29, 30, 31]. In our calculation, they arise naturally from the wide distribution of persistence times. The CTRW reproduces the true ageing dynamics exceedingly well, with no fit parameters. The same conclusion holds for the binary mixture. As in the trap model, the power law form of is sufficient to generate ageing. Distributions evolved via the CTRW from the shortest waiting time alone closely track measured distributions. Ageing emerges here as a self-generating, dynamical phenomenon.

Our results clearly show that there are two relevant timescales in the relaxation dynamics of the glass: the first hop time which experiences ageing, and the subsequent hop times which are independent. This would seem to refute the hypothesis that the dynamics are universally reparametrized with waiting time as has been suggested for spin glasses [12], and may serve to explain some deviations from local dynamical scaling that have been observed. In ref. [12], the authors note that the distribution of local correlation functions for a binary Lennard-Jones glass superimposes to first order at times where the global correlation function is equal. However, good collapse was found only for large and small global correlation with significant deviations for and .

We further investigate the local dynamical scaling by using the CTRW model to calculate the van Hove distribution for three waiting times and three different values of the mean in Fig. 5. Note that here is generated by evolving the system directly in the ageing CTRW. The superposition of the van Hove functions displays the same behavior reported in ref. [12]. At short times (large ), the van Hove functions closely superimpose at times of equal mean. In this regime, few atoms have hopped and their mean squared displacement and van Hove function are dominated by the distribution , which scales with . For the mid-range of the correlation functions, the overlap of the van Hove functions becomes weaker. The long waiting time distributions have a sharper trapped peak and wider tails than the short function. This is due to the fact that the width of the tails is controlled by , whereas the height of the caged peak is controlled by . At very long times (small ), the mean number of hops per atom is much greater than one, and the waiting time independent behavior begins to dominate the dynamics. In this case, the mean squared displacement curves merge and there is trivial superposition of the van Hove functions.

To summarize, we have presented a complete characterization of particle diffusion in ageing structural glasses. We find that on a microscopic level, time is not simply reparameterized during ageing. Only the first hop time depends on the waiting time, and the particle subsequently “forgets” its age. The evolution of the first hop time distributions is remarkably close to that assumed in the trap model with annealed disorder. Although this model is a rather abstract picture of glassy dynamics, our results indicate that it may be closer to physical reality on a single particle level than previously thought. An explicit mapping of the MD data onto a CTRW with measured distributions of hop times and displacements is able to completely reproduce all dynamical features of the relaxing glass, and provides an explanation for the scaling with and small deviations thereof.

We thank the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canadian Foundation for Innovation (CFI) for financial support.

## References

- [1] \NameBarrat J.-L., Feigelman M. Kurchan J. (Editors) \BookSlow relaxations and nonequilibrium dynamics in condensed matter (Springer, Berlin) 2003.
- [2] \NameEdiger M. D. \REVIEWAnnu. Rev. Phys. Chem. 51200099.
- [3] \NameKob W., Donati C., Plimpton S. J., Poole P. H. Glotzer S. C. \REVIEWPhys. Rev. Lett. 7919972827.
- [4] \NameShintani H. Tanaka H. \REVIEWNature Physics 22006200.
- [5] \NameVollmayr-Lee K. Zippelius A. \REVIEWPhys. Rev. E 722005041507.
- [6] \NameCourtland R. E. Weeks E. R. \REVIEWJ. Phys.:Condens. Matter 152003S359.
- [7] \NameStruik L. C. E. \BookPhysical Aging in Amorphous Polymers and Other Materials (Elsevier, Amsterdam) 1978.
- [8] \NameCipelletti L. Ramos L. \REVIEWJ. Phys.: Condens. Matter 172005R253.
- [9] \NameWang P., Song C. Makse H. \REVIEWNature Physics 22006526.
- [10] \NameMartinez V. A., Bryant G. van Megen W. \REVIEWPhys. Rev. Lett. 1012008135702.
- [11] \NameCloitre M., Borrega R. Leibler L. \REVIEWPhys. Rev. Lett. 8520004819.
- [12] \NameCastillo H. E. Parsaeian A. \REVIEWNature Physics 3200726.
- [13] \NameWarren M. Rottler J. \REVIEWPhys. Rev. E 762007031802.
- [14] \NameCastillo H., Chamon C., Cugliandolo L., Iguain J. Kennett M. \REVIEWPhys. Rev. B 682003134442.
- [15] \NameParsaeian A. Castillo H. \REVIEWPhys. Rev. Lett. 1022009055704.
- [16] \NameKob W. Andersen H. C. \REVIEWPhys. Rev. E 5119954626.
- [17] \NameBennemann C., Paul W., Binder K. Dünweg B. \REVIEWPhys. Rev. E 571998843.
- [18] \NameBaschnagel J., Bennemann C., Paul W. Binder K. \REVIEWJ. Phys.: Condens. Matter 1220006365.
- [19] The choice of time window affects the width of the resulting displacement distribution and the normalization constant of the persistence time distribution, but the slope of the persistence time and the distribution of first hop times are unaffected.
- [20] Increasing the threshold broadens both hop time and displacement distributions; however, thresholds in a range of 1.5 to 2.5 give self-consistent results for the mean-squared displacement in the continuous time random walk. A threshold of two was chosen as a compromise between resolving ageing and reducing noise.
- [21] \NameHedges L., Maibaum L., Chandler D. Garrahan J. \REVIEWJ. Chem. Phys 1272007211101.
- [22] \NameChaudhuri P., Gao Y., Berthier L., Kilfoil M. Kob W. \REVIEWJ. Phys.: Condens. Matt. 202008244126.
- [23] \NameDoliwa B. Heuer A. \REVIEWPhys. Rev. Lett. 912003235501.
- [24] \NameRubner O. Heuer A. \REVIEWPhys. Rev. E 782008011504.
- [25] \NameHeuer A. \REVIEWJ. Phys.: Condens. Matt. 202008373101.
- [26] \NameMonthus C. Bouchaud J.-P. \REVIEWJ. Phys. A: Math. Gen. 2919963847.
- [27] \NameFielding S. M., Sollich P. Cates M. E. \REVIEWJ. Rheology 442000329.
- [28] \NameChaudhuri P., Berthier L. Kob W. \REVIEWPhys. Rev. Lett. 992007060604.
- [29] \NameLanger J. S. Swagatam M. \REVIEWPhys. Rev. E 772008061505.
- [30] \NameLanger J. S. \REVIEWPhys. Rev. E 782008051115.
- [31] \NameSaltzman E. Schweizer K. S. \REVIEWPhys. Rev. E 772008051504.