Relaxation dynamics and their spatial distribution in a two-dimensional glass-forming mixture

Aging and relaxation in glass-forming systems

We propose that there exists a generic class of glass-forming systems that have competing states (of crystalline order or not) which are locally close in energy to the ground state (which is typically unique). Upon cooling, such systems exhibit patches (or clusters) of these competing states which become locally stable in the sense of having a relatively high local shear modulus. It is in between these clusters where aging, relaxation, and plasticity under strain can take place. We demonstrate explicitly that relaxation events that lead to aging occur where the local shear modulus is low (even negative) and result in an increase in the size of local patches of relative order. We examine the aging events closely from two points of view. On the one hand we show that they are very localized in real space, taking place outside the patches of relative order, and from the other point of view we show that they represent transitions from one local minimum in the potential surface to another. This picture offers a direct relation between structure and dynamics, ascribing the slowing down in glass-forming systems to the reduction in relative volume of the amorphous material which is liquidlike. While we agree with the well-known Adam-Gibbs proposition that the slowing down is due to an entropic squeeze (a dramatic decrease in the number of available configurations), we do not agree with the Adam-Gibbs (or the Volger-Fulcher) formulas that predict an infinite relaxation time at a finite temperature. Rather, we propose that generically there should be no singular crisis at any finite temperature: the relaxation time and the associated correlation length (average cluster size) increase at most superexponentially when the temperature is lowered.

Decoupling of exchange and persistence times in atomistic models of glass formers

With molecular dynamics simulations of a fluid mixture of classical particles interacting with pairwise additive Weeks-Chandler-Andersen potentials, we consider the time series of particle displacements and thereby determine the distributions for local persistence times and local exchange times. These basic characterizations of glassy dynamics are studied over a range of supercooled conditions and were shown to have behaviors, most notably decoupling, similar to those found in kinetically constrained lattice models of structural glasses. Implications are noted.

Mechanical properties of glass forming systems

We address the interesting temperature range of a glass forming system where the mechanical properties are intermediate between those of a liquid and a solid. We employ an efficient Monte Carlo method to calculate the elastic moduli, and show that in this range of temperatures the moduli are finite for short times and vanish for long times, where short and long depend on the temperature. By invoking some exact results from statistical mechanics we offer an alternative method to compute shear moduli using molecular dynamics simulations, and compare those to the Monte Carlo method. The final conclusion is that these systems are not "viscous fluids" in the usual sense, as their actual time-dependence concatenates solid-like materials with varying local shear moduli.

Connection of translational and rotational dynamical heterogeneities with the breakdown of the Stokes-Einstein and Stokes-Einstein-Debye relations in water

We study the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, Dt=kBT/6pietaR and Dr=kBT/8pietaR3, where Dt and Dr are the translational and rotational diffusivity, respectively, T is the temperature, eta the viscosity, kB the Boltzmann constant, and R the "molecular" radius. Our results are based on molecular dynamics simulations of the extended simple point charge model of water. We find that both the SE and SED relations break down at low temperature. To explore the relationship between these breakdowns and dynamical heterogeneities (DHs), we also calculate the SE and SED relations for subsets of the 7% "fastest" and 7% "slowest" molecules. We find that the SE and SED relations break down in both subsets, and that the breakdowns occur on all scales of mobility. Thus these breakdowns appear to be generalized phenomena, in contrast with a view where only the most mobile molecules are the origin of the breakdown of the SE and SED relations, embedded in an inactive background where these relations hold. At low temperature, the SE and SED relations in both subsets of molecules are replaced with "fractional" SE and SED relations, Dt approximately (tau/T)-xit and Dr approximately (tau/T)-xir, where xit approximately 0.84(<1) and xir approximately 0.75(<1). We also find that there is a decoupling between rotational and translational motion, and that this decoupling occurs in both the fastest and slowest subsets of molecules. Further, we find that, the decoupling increases upon cooling, but that the probability of a molecule being classified as both translationally and rotationally fastest also increases. To study the effect of time scale for SE and SED breakdown and decoupling, we introduce a time-dependent version of the SE and SED relations, and a time-dependent function that measures the extent of decoupling. Our results suggest that both the decoupling and SE and SED breakdowns originate at a time scale corresponding to the end of the cage regime, when diffusion starts. This is also the time scale when the DHs are more relevant. Our work also demonstrates that selecting DHs on the basis of translational or rotational motion more strongly biases the calculation of diffusion constants than other dynamical properties such as relaxation times.

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.

Spontaneous and induced dynamic fluctuations in glass formers. I. General results and dependence on ensemble and dynamics

We study theoretically and numerically a family of multipoint dynamic susceptibilities that quantify the strength and characteristic length scales of dynamic heterogeneities in glass-forming materials. We use general theoretical arguments (fluctuation-dissipation relations and symmetries of relevant dynamical field theories) to relate the sensitivity of averaged two-time correlators to temperature and density to spontaneous fluctuations of the local dynamics. Our theoretical results are then compared to molecular dynamics simulations of the Newtonian, Brownian, and Monte Carlo dynamics of two representative glass-forming liquids, a fragile binary Lennard-Jones mixture, and a model for the strong glass-former silica. We justify in detail the claim made by Berthier et al. [Science 310, 1797 (2005)] that the temperature dependence of correlation functions allows one to extract useful information on dynamic length scales in glassy systems. We also discuss some subtle issues associated with the choice of microscopic dynamics and of statistical ensemble through conserved quantities, which are found to play an important role in determining dynamic correlations.

On the study of collective dynamics in supercooled liquids through the statistics of the isoconfigurational ensemble

The use of the isoconfigurational ensemble to explore structure-dynamic correlations in supercooled liquids is examined. The statistical error of the dynamic propensity and its spatial distribution are determined. The authors present the spatial distribution of the particle non-Gaussian parameter as a measure of the intermittency with which particles exhibit their propensity for motion. The ensemble average of the direction of particle motion is introduced to establish the anisotropy of the dynamic propensity.

Fluorescence lifetime of a single molecule as an observable of meta-basin dynamics in fluids near the glass transition

Using single molecule spectroscopy, we show that the fluorescence lifetime trajectories of single probe molecules embedded in a glass-forming polymer melt exhibit strong fluctuations of a hopping character. Using molecular dynamics simulations targeted to explain these experimental observations, we show that the lifetime fluctuations correlate strongly with the average square displacement function of the matrix particles. The latter observable is a direct probe of the meta-basin transitions in the potential energy landscape of glass-forming liquids. We thus show here that single molecule experiments can provide detailed microscopic information on system properties that hitherto have been accessible via computer simulations only.

Neutron reflectivity measurements of the translational motion of tris(naphthylbenzene) at the glass transition temperature

The translational dynamics of the low molecular weight glass-former tris(naphthylbenzene) have been studied on the length scale of a few nanometers at the glass transition temperature Tg. Neutron reflectivity was used to measure isotopic interdiffusion of multilayer samples created by physical vapor deposition. Deposition with the substrate held at Tg-6 K allows observation of dynamics characterizing the equilibrium supercooled liquid. The diffusion coefficient measured at q = 0.03 A(-1) was determined to be 1x10(-17) cm2/s at 342 K (Tg). The self-part of the intermediate scattering function I(s)(q,t) decays exponentially. Samples deposited well below Tg show a substantial thermal history effect during subsequent translational motion at Tg.

Predicting the long-time dynamic heterogeneity in a supercooled liquid on the basis of short-time heterogeneities

We report that the local Debye-Waller factor in a simulated 2D glass-forming mixture exhibits significant spatial heterogeneities and that these short-time fluctuations provide an excellent predictor of the spatial distribution of the long-time dynamic propensities. In contrast, the potential energy per particle of the inherent structure does not correlate well with the spatially distributed dynamics.

Low-frequency dynamical heterogeneity in simulated amorphous Ni0.5Zr0.5 below its glass temperature: correlations with cage volume and local order fluctuations

From molecular dynamics simulations results are reported concerning correlations between low-frequency (lf) heterogeneous dynamics in simulated Ni0.5Zr0.5 melts at 700, 760, and 810 K, which means around the Kauzmann temperature of the model, TK approximately 750 K. A method is presented to separate lf dynamics, reflecting the slow relaxation dynamics in the vitrifying melt, and high-frequency (hf) dynamics, characteristic of the thermal fluctuations at the considered temperatures. By means of a suitable quantitative measure of the distribution of heterogeneous lf dynamics in space and time, correlation parameters are evaluated between the spatial distribution of lf dynamics and structural inhomogeneities in the thermodynamically homogeneous melt. Relevant correlations are found between lf dynamics and some involved structure quantities such as the cage volume around Ni atoms, Omega Ni, or the Theta Ni parameter which reflects the geometry of the nearest-neighbor cage around Ni atoms. Further, at 810 K there is a weak correlation between heterogeneous dynamics and fluctuations of the mean potential energy per atom and a comparable weak anticorrelation with the particle density and Ni-atom density inhomogeneities, where these three correlations decrease with decreasing temperature. The present results indicate the existence of long-living regions of enhanced Omega Ni in the structure, which may act as regions of preferential initiation of irreversible lf dynamics and slow relaxation.

Time scale for the onset of Fickian diffusion in supercooled liquids

We propose a quantitative measure of a time scale on which Fickian diffusion sets in for supercooled liquids, and we use Brownian dynamics computer simulations to determine the temperature dependence of this onset time in a Lennard-Jones binary mixture. The time for the onset of Fickian diffusion ranges between 6.5 and 31 times the relaxation time (the relaxation time is the characteristic relaxation time of the incoherent intermediate scattering function). The onset time increases faster with decreasing temperature than the relaxation time. Mean-squared displacement at the onset time increases with decreasing temperature.

Heterogeneity and growing length scales in the dynamics of kinetically constrained lattice gases in two dimensions

We study dynamical heterogeneity and growing dynamical length scales in two kinetically constrained models, namely, the one- and two-vacancy assisted triangular lattice gases. One of the models is a strong glassformer and the other is a fragile glassformer. Both exhibit heterogeneous dynamics with broadly distributed time scales as seen in the distribution of persistence times. We show that the Stokes-Einstein relation is violated, to a greater degree in the fragile glassformer, and show how this violation is related to dynamic heterogeneity. We extract dynamical length scales from structure factors of mobile particles and show, quantitatively, the growth of this length scale as density increases. We comment on how the scaling of lengths and times in these models relates to that in facilitated spin models of glasses.

Heterogeneities in the glassy state

We study heterogeneities in a binary Lennard-Jones system below the glass transition using molecular dynamics simulations. We identify mobile and immobile particles and measure their distribution of vibrational amplitudes. For temperatures near the glass transition the distribution of vibrational amplitudes obeys scaling and compares reasonably well with a mean-field theory for the amorphous solid state. To investigate correlations among the immobile and mobile particles we identify clusters and analyze their size and shape. For a fixed number of immobile particles we observe that the immobile particles cluster more strongly together as the temperature is increased which allows the particles to block each other more effectively and to therefore stay immobile. For the mobile particles, on the other hand, the clustering is most pronounced at low temperatures, indicating that mobility at low temperatures can only be sustained in cooperative motion.

Potential-energy landscape of a supercooled liquid and its resemblance to a collection of traps

It is analyzed whether the potential energy landscape of a glass-forming system can be effectively mapped on a random model which is described in statistical terms. For this purpose we generalize the simple trap model of Monthus and Bouchaud [J. Phys. A 29, 3847 (1996)] by dividing the total system into M weakly interacting identical subsystems, each being described in terms of a trap model. The distribution of traps in this extended trap model (ETM) is fully determined by the thermodynamics of the glass former. The dynamics is described by two adjustable parameters, one characterizing the common energy level of the barriers, the other the strength of the interaction. The comparison is performed for the standard binary mixture Lennard-Jones system with 65 particles. The metabasins, identified in our previous work, are chosen as traps. Comparing molecular dynamics simulations of the Lennard-Jones system with Monte Carlo calculations of the ETM allows one to determine the adjustable parameters. Analysis of the first moment of the waiting distribution yields an optimum agreement when choosing M approximately 3 subsystems. Comparison with the second moment of the waiting time distribution, reflecting dynamic heterogeneities, indicates that the sizes of the subsystems may fluctuate.

Two-Gaussian excitations model for the glass transition

We develop a modified "two-state" model with Gaussian widths for the site energies of both ground and excited states, consistent with expectations for a disordered system. The thermodynamic properties of the system are analyzed in configuration space and found to bridge the gap between simple two-state models ("logarithmic" model in configuration space) and the random energy model ("Gaussian" model in configuration space). The Kauzmann singularity given by the random energy model remains for very fragile liquids but is suppressed or eliminated for stronger liquids. The sharp form of constant-volume heat capacity found by recent simulations for binary mixed Lennard-Jones and soft-sphere systems is reproduced by the model, as is the excess entropy and heat capacity of a variety of laboratory systems, strong and fragile. The ideal glass in all cases has a narrow Gaussian, almost invariant among molecular and atomic glassformers, while the excited-state Gaussian depends on the system and its width plays a role in the thermodynamic fragility. The model predicts the possibility of first-order phase transitions for fragile liquids. The analysis of laboratory data for toluene and o-terphenyl indicates that fragile liquids resolve the Kauzmann paradox by a first-order transition from supercooled liquid to ideal-glass state at a temperature between T(g) and Kauzmann temperature extrapolated from experimental data. We stress the importance of the temperature dependence of the energy landscape, predicted by the fluctuation-dissipation theorem, in analyzing the liquid thermodynamics.

Dynamics in a two-dimensional core-softened fluid

The dynamical properties of a model one-component core-softened fluid with purely repulsive interactions are found to be very complex. At low temperature the fluid structure exhibits cluster motifs including dimers, stripes, and polygons, depending on density. Single-particle diffusion and the velocity, shear-stress, and wave-vector-dependent current correlation functions have all been calculated using molecular dynamics simulations. The results highlight the presence of well-resolved single-particle and collective motions, which is remarkable for what is essentially a "simple" one-component fluid.

Short-ranged attractions in jammed liquids: how cooling can melt a glass

We demonstrate that an extended picture of kinetic constraints in glass-forming liquids is sufficient to explain dynamic anomalies observed in dense suspensions of strongly attracting colloidal particles. We augment a simple model of heterogeneous relaxation with static attractions between facilitating excitations, in a way that mimics the structural effect of short-ranged interparticle attractions. The resulting spatial correlations among facilitated and unfacilitated regions give rise to relaxation mechanisms that account for nonmonotonic dependence of relaxation times on attraction strength as well as logarithmic decay of density correlations in time. These unusual features are a simple consequence of spatial segregation of kinetic constraints, suggesting an alternative physical perspective on attractive colloids than that suggested by mode-coupling theory. Based on the behavior of our model, we predict a crossover from super-Arrhenius to Arrhenius temperature dependence as attractions become dominant at fixed packing fraction.

Topological lifetimes of polydisperse colloidal hard spheres at a wall

Confocal scanning laser microscopy was used to study the behavior of dense suspensions of model colloidal hard spheres at a single wall. Due to the slight polydispersity, our system shows a reentrant melting transition at high densities involving a hexatic structure [Phys. Rev. Lett 92, 195702 (2004)]]. The reentrant melting transition is accompanied by an increase in the mean-squared displacement. The correlation between structure and dynamics was quantitatively analyzed on a single-particle level. In particular, the topological lifetime, being the average time that a particle spends having the same coordination number, is determined for all coordination numbers and as a function of volume fraction. The defective (non-sixfold-coordinated) particles exhibit shorter lifetimes than sixfold-coordinated particles, indicating that the mobility of the system is larger at or close to defective particles. The lifetime itself is a strong function of volume fraction. In particular, the global behavior of the mean-squared displacement is proportional to the hopping frequency (the inverse of the lifetime), showing that particles changing their coordination number contribute most to the local mobility.

Lifetime of dynamic heterogeneities in a binary Lennard-Jones mixture

A four-time correlation function was calculated using a computer simulation of a binary Lennard-Jones mixture. The information content of the four-time correlation function is similar to that of four-time correlation functions measured in NMR experiments. The correlation function selects a subensemble and analyzes its dynamics after some waiting time. The lifetime of the subensemble selected by the four-time correlation function is calculated, and compared to the lifetimes of slow subensembles selected using two different definitions of mobility, and to the alpha relaxation time.

How reproducible are dynamic heterogeneities in a supercooled liquid?

The particle dynamics in a liquid exhibits a transient spatial distribution of dynamic heterogeneities. The relationship between this kinetic structure and the underlying particle configuration remains an outstanding problem. In this Letter, we present a general simulation technique for identifying the features of the dynamic heterogeneity which arise due to a specific configuration, as distinct from the random spatial variation due to the intermittent particle dynamics.

Excitation lines and the breakdown of Stokes-Einstein relations in supercooled liquids

By applying the concept of dynamical facilitation and analyzing the excitation lines that result from this facilitation, we investigate the origin of decoupling of transport coefficients in supercooled liquids. We illustrate our approach with two classes of models. One depicts diffusion in a strong glass former, and the other in a fragile glass former. At low temperatures, both models exhibit violation of the Stokes-Einstein relation, D approximately tau(-1), where D is the self-diffusion constant and tau is the structural relaxation time. In the strong case, the violation is sensitive to dimensionality d, going as D approximately tau(-2/3) for d=1 and as D approximately tau(-0.95) for d=3. In the fragile case, however, we argue that dimensionality dependence is weak, and show that for d=1, D approximately tau(-0.73). This scaling for the fragile case compares favorably with the results of a recent experimental study for a three-dimensional fragile glass former.

Particle dynamics and the development of string-like motion in a simulated monoatomic supercooled liquid

The microscopic details of local particle dynamics is studied in a glass-forming one component supercooled liquid modeled by a Dzugutov potential developed for simple metallic glass formers. Our main goal is to investigate particle motion in the supercooled liquid state, and to ascertain the extent to which this motion is cooperative and occurring in quasi-one-dimesional, string-like paths. To this end we investigate in detail the mechanism by which particles move along these paths. In particular, we show that the degree of coherence--that is, simultaneous motion by consecutive particles along a string--depends on the length of the string. For short strings, the motion is highly coherent. For longer strings, the motion is highly coherent only within shorter segments of the string, which we call "microstrings." Very large strings may contain several microstrings within which particles move simultaneously, but individual microstrings within a given string are temporally uncorrelated with each other. We discuss possible underlying mechanism for this complex dynamical behavior, and examine our results in the context of recent work by Garrahan and Chandler [Phys. Rev. Lett. 89, 035704 (2002)] in which dynamic facilitation plays a central role in the glass transition.

Cooperative dynamics in two dimensions

We report results from molecular dynamics simulations of cooperative motion in a quasi-two-dimensional system of colloid particles. We find that the onset of the deviation of the single-particle displacement distribution from Gaussian form starts in the liquid phase and extends, with increasing magnitude, through the hexatic phase into the crystalline phase. The time for which the deviation is maximum increases exponentially with the density. As the density increases toward the hexatic phase a third dynamical relaxation mode emerges. We argue that the collective motion is generated by superpositions of instantaneous normal mode vibrations, with lifetimes that increase with the density, along paths with strong bond-orientation correlation.

Freezing transition and correlated motion in a quasi-two-dimensional colloid suspension

Recent experiments have demonstrated that the deviation of the single-particle displacement distribution from Gaussian form in a dense quasi-two-dimensional colloid suspension is a result of heterogenous dynamics that involves cooperative motions of neighboring colloid particles [J. Chem. Phys. 47, 9142 (2001)]. In this paper, we report the results of molecular dynamics (MD) simulations of a quasi-two-dimensional assembly of nearly hard-sphere colloid particles. The colloid-colloid interaction we use is short ranged and everywhere repulsive; it is related to the Marcus-Rice (MR) and modified MR interactions used in a previous study [Phys. Rev. E 58, 7529 (1998)]. As is the case for those systems, the one we study supports liquid, hexatic, and solid phases. Our calculations show that the deviation of the single-particle displacement distribution from Gaussian form is present in the liquid phase, and that a sharp increase in its magnitude occurs at the liquidus density and extends into the crystalline phase. For densities greater than the liquidus density we find three dynamical relaxation processes that include, at intermediate times, a slowing down in the rate of growth of the diffusive displacement of a particle due to the cage effect. As the density increases toward the solidus density, the dependence of the mean squared displacement on time, at intermediate times, changes from sublinear to zero. The onset of the long-time relaxation mode corresponds to the time at which the deviation of the particle displacement distribution from Gaussian form is a maximum. At this time, which increases exponentially with the density, the self-part of the van Hove function exhibits multiple maxima with respect to r while the distinct part of the van Hove function is a maximum at the origin, thereby signaling jump dynamics. At long times the particle mean square displacement has diffusive character at all densities including solid phase densities. A remarkable feature of our findings is the continuity of character of the particle displacement from the liquid phase through the hexatic phase and into the solid phase. Cooperative jumps that lead to diffusive process in crystals can be explained by a mechanism that involves many such correlated hops in random locations and random directions (but along the crystallographic axes) thereby generating effective random walk behavior. We argue that the collective motion we have found is generated by superpositions of instantaneous normal mode vibrations along diffusive paths. The diffusive paths are along the directions with strong bond orientation correlation, and start to grow in amplitude rapidly on entry into the hexatic phase.

Real space origin of temperature crossovers in supercooled liquids

We show that the various crossovers between dynamical regimes observed in experiments and simulations of supercooled liquids can be explained in simple terms from the existence and statistical properties of dynamical heterogeneities. We confirm that dynamic heterogeneity is responsible for the slowing down of glass formers at temperatures well above the dynamic singularity Tc predicted by mode-coupling theory. Our results imply that activated processes govern the long-time dynamics even in the temperature regime where they are neglected by mode-coupling theory. We show that alternative interpretations based on topographic properties of the potential energy landscape are inefficient ways of describing simple physical features which are naturally accounted for within our approach. We show in particular that the reported links between mode coupling and landscape singularities do not exist.

Diffusion and viscosity in a supercooled polydisperse system

We have carried out extensive molecular dynamics simulations of a supercooled polydisperse Lennard-Jones liquid with large variations in temperature at a fixed pressure. The particles in the system are considered to be polydisperse in both size and mass. The temperature dependence of dynamical properties such as the viscosity (eta) and the self-diffusion coefficients (D(i)) of different size particles is studied. Both viscosity and diffusion coefficients show super-Arrhenius temperature dependence and fit well to the well-known Vogel-Fulcher-Tammann equation. Within the temperature range investigated, the value of Angell's fragility parameter (D approximately 1.4) classifies the present system as a very fragile liquid. The critical temperature for diffusion (T(D(i))(o)) increases with the size of the particles. The critical temperature for viscosity (T(eta)(o)) is larger than that for diffusion, and sizable deviations appear for the smaller size particles, implying a decoupling of translational diffusion from viscosity in deeply supercooled liquids. Indeed, the diffusion shows markedly non-Stokesian behavior at low temperatures where a highly nonlinear dependence on size is observed. An inspection of the trajectories of the particles shows that at low temperatures the motions of both the smallest and largest size particles are discontinuous (jump type). However, the crossover from continuous Brownian to large length hopping motion takes place at shorter time scales for the smaller size particles.

Self-diffusion of tris-naphthylbenzene near the glass transition temperature

We present a direct measurement of self-diffusion of a single-component glass-forming liquid at the glass transition temperature. Forward recoil spectrometry is used to measure the concentration profiles of deuterio and protio 1,3-bis-(1-naphthyl)-5-(2-naphthyl)benzene (TNB) following annealing-induced diffusion in a vapor-deposited bilayer. These experiments extend the range of measured diffusion coefficients in TNB by 6 orders of magnitude. The results indicate a decoupling of translational diffusion coefficients from viscosity or rotation. At T(g), D(T) is 400 times larger than expected from the Stokes-Einstein equation.

Properties of cage rearrangements observed near the colloidal glass transition

We use confocal microscopy to study particle motion in colloidal systems. Near the glass transition, motion is inhibited, as particles spend time trapped in transient "cages" formed by neighboring particles. We measure the cage sizes and lifetimes, which, respectively, shrink and grow as the glass transition approaches. Cage rearrangements are more prevalent in regions with lower concentrations and higher disorder. Neighboring rearranging particles typically move in parallel directions, although a nontrivial fraction moves in antiparallel directions, usually from particle pairs with initial separations corresponding to local maxima and minima of the pair correlation function g(r), respectively.

Spatially correlated dynamics in a simulated glass-forming polymer melt: analysis of clustering phenomena

In recent years, experimental and computational studies have demonstrated that the dynamics of glass-forming liquids are spatially heterogeneous, exhibiting regions of temporarily enhanced or diminished mobility. Here we present a detailed analysis of dynamical heterogeneity in a simulated "bead-spring" model of a low-molecular-weight polymer melt. We investigate the transient nature and size distribution of clusters of "mobile" chain segments (monomers) as the polymer melt is cooled toward its glass transition. We also explore the dependence of this clustering on the way in which the mobile subset is defined. We show that the mean cluster size is time dependent with a peak at intermediate time, and that the mean cluster size at the peak time grows with decreasing temperature T. We show that for each T a particular fraction of particles maximizes the mean cluster size at some characteristic time, and this fraction depends on T. The growing size of the clusters demonstrates the growing range of correlated motion, previously reported for this same system [C. Beneman et al. Nature (London) 399, 246 (1999)]. The distribution of cluster sizes approaches a power law near the mode-coupling temperature, similar to behavior reported for a simulated binary mixture and a dense colloidal suspension, but with a different exponent. We calculate the correlation length of the clusters, and show that it exhibits similar temperature- and time-dependent behavior as the mean cluster size, with a maximum at intermediate time. We show that the characteristic time of the maximum cluster size follows the scaling predicted by mode-coupling theory (MCT) for the beta time scale, revealing a possible connection between spatially heterogeneous dynamics and MCT.

Spatially heterogeneous dynamics in supercooled liquids

Although it has long been recognized that dynamics in supercooled liquids might be spatially heterogeneous, only in the past few years has clear evidence emerged to support this view. As a liquid is cooled far below its melting point, dynamics in some regions of the sample can be orders of magnitude faster than dynamics in other regions only a few nanometers away. In this review, the experimental work that characterizes this heterogeneity is described. In particular, the following questions are addressed: How large are the heterogeneities? How long do they last? How much do dynamics vary between the fastest and slowest regions? Why do these heterogeneities arise? The answers to these questions influence practical applications of glass-forming materials, including polymers, metallic glasses, and pharmaceuticals.