Relaxation in a glassy binary mixture: mode-coupling-like power laws, dynamic heterogeneity, and a new non-Gaussian parameter

Structural relaxation of polydisperse hard spheres: comparison of the mode-coupling theory to a Langevin dynamics simulation

We analyze the slow glassy structural relaxation as measured through collective and tagged-particle density correlation functions obtained from Brownian dynamics simulations for a polydisperse system of quasi-hard spheres in the framework of the mode-coupling theory (MCT) of the glass transition. Asymptotic analyses show good agreement for the collective dynamics when polydispersity effects are taken into account in a multicomponent calculation, but qualitative disagreement at small q when the system is treated as effectively monodisperse. The origin of the different small-q behavior is attributed to the interplay between interdiffusion processes and structural relaxation. Numerical solutions of the MCT equations are obtained taking properly binned partial static structure factors from the simulations as input. Accounting for a shift in the critical density, the collective density correlation functions are well described by the theory at all densities investigated in the simulations, with quantitative agreement best around the maxima of the static structure factor and worst around its minima. A parameter-free comparison of the tagged-particle dynamics however reveals large quantitative errors for small wave numbers that are connected to the well-known decoupling of self-diffusion from structural relaxation and to dynamical heterogeneities. While deviations from MCT behavior are clearly seen in the tagged-particle quantities for densities close to and on the liquid side of the MCT glass transition, no such deviations are seen in the collective dynamics.

Numerical and theoretical study of a monodisperse hard-sphere glass former

There exists a variety of theories of the glass transition and many more numerical models. But because the models need built-in complexity to prevent crystallization, comparisons with theory can be difficult. We study the dynamics of a deeply supersaturated monodisperse four-dimensional (4D) hard-sphere fluid, which has no such complexity, but whose strong intrinsic geometrical frustration inhibits crystallization, even when deeply supersaturated. As an application, we compare its behavior to the mode-coupling theory (MCT) of glass formation. We find MCT to describe this system better than any other structural glass formers in lower dimensions. The reduction in dynamical heterogeneity in 4D suggested by a milder violation of the Stokes-Einstein relation could explain the agreement. These results are consistent with a mean-field scenario of the glass transition.

MD simulation of concentrated polymer solutions: structural relaxation near the glass transition

We examine by molecular dynamics simulations the relaxation of polymer-solvent mixtures close to the glass transition. The simulations employ a coarse-grained model in which polymers are represented by bead-spring chains and solvent particles by monomers. The interaction parameters between polymer and solvent are adjusted such that mixing is favored. We find that the mixtures have one glass transition temperature T(g) or critical temperature T(c) of mode-coupling theory (MCT). Both T(g) and T(c) (> T(g)) decrease with increasing solvent concentration φ(S). The decrease is linear for the concentrations studied (up to φ(S) = 25%). Above T(c) we explore the structure and relaxation of the polymer-solvent mixtures on cooling. We find that, if the polymer solution is compared to the pure polymer melt at the same T, local spatial correlations on the length scale of the first peak of the static structure factor S(q) are reduced. This difference between melt and solution is largely removed when comparing the S(q) of both systems at similar distance to the respective T(c). Near T(c) we investigate dynamic correlation functions, such as the incoherent intermediate scattering function φ(q)(s)(t), mean-square displacements of the monomers and solvent particles, two non-Gaussian parameters, and the probability distribution P(ln r; t) of the logarithm of single-particle displacements. In accordance with MCT we find, for instance, that φ(q)(s)(t) obeys the time-temperature superposition principle and has α relaxation times τ(q)(s) which are compatible with a power law increase close (but not too close) to T(c). In divergence to MCT, however, the increase of τ(q)(s) depends on the wavelength q, small q values having weaker increase than large ones. This decoupling of local and large-length scale relaxation could be related to the emergence of dynamic heterogeneity at low T. In the time window of the α relaxation an analysis of P(ln r; t) reveals a double-peak structure close to T(c). The first peak corresponds to "slow" particles (monomer or solvent) which have not moved much farther than 10% of their diameter in time t, whereas the second occurs at distances of the order of the particle diameter. These "fast" particles have succeeded in leaving their nearest-neighbor cage in time t. The simulation thus demonstrates that large fluctuations in particle mobility accompany the final structural relaxation of the cold polymer solution in the vicinity of the extrapolated T(c).

Connections of activated hopping processes with the breakdown of the Stokes-Einstein relation and with aspects of dynamical heterogeneities

We develop an extended version of the mode-coupling theory (MCT) for glass transition, which incorporates activated hopping processes via the dynamical theory originally formulated to describe diffusion-jump processes in crystals. The dynamical-theory approach adapted here to glass-forming liquids treats hopping as arising from vibrational fluctuations in the quasiarrested state where particles are trapped inside their cages, and the hopping rate is formulated in terms of the Debye-Waller factors characterizing the structure of the quasiarrested state. The resulting expression for the hopping rate takes an activated form, and the barrier height for the hopping is "self-generated" in the sense that it is present only in those states where the dynamics exhibits a well defined plateau. It is discussed how such a hopping rate can be incorporated into MCT so that the sharp nonergodic transition predicted by the idealized version of the theory is replaced by a rapid but smooth crossover. We then show that the developed theory accounts for the breakdown of the Stokes-Einstein relation observed in a variety of fragile glass formers. It is also demonstrated that characteristic features of dynamical heterogeneities revealed by recent computer simulations are reproduced by the theory. More specifically, a substantial increase of the non-Gaussian parameter, double-peak structure in the probability distribution of particle displacements, and the presence of a growing dynamic length scale are predicted by the extended MCT developed here, which the idealized version of the theory failed to reproduce. These results of the theory are demonstrated for a model of the Lennard-Jones system, and are compared with related computer-simulation results and experimental data.

Exploring the potential energy landscape of glass-forming systems: from inherent structures via metabasins to macroscopic transport

In this review a systematic analysis of the potential energy landscape (PEL) of glass-forming systems is presented. Starting from the thermodynamics, the route towards the dynamics is elucidated. A key step in this endeavor is the concept of metabasins. The relevant energy scales of the PEL can be characterized. Based on the simulation results for some glass-forming systems one can formulate a relevant model system (ideal Gaussian glass-former) which can be treated analytically. The macroscopic transport can be related to the microscopic hopping processes, using either the strong relation between energy (thermodynamics) and waiting times (dynamics) or, alternatively, the concepts of the continuous-time random walk. The relation to the geometric properties of the PEL is stressed. The emergence of length scales within the PEL approach as well as the nature of finite-size effects is discussed. Furthermore, the PEL view is compared to other approaches describing the glass transition.

Glassy dynamics and kinetic vitrification of isotropic suspensions of hard rods

A nonlinear Langevin equation (NLE) theory for the translational center-of-mass dynamics of hard nonspherical objects has been applied to isotropic fluids of rigid rods. The ideal kinetic glass transition volume fraction is predicted to be a monotonically decreasing function beyond an aspect ratio of two. The functional form of the decrease is weaker than the inverse aspect ratio. Vitrification occurs at lower volume fractions for corrugated tangent bead rods compared to their smooth spherocylinder analogs. The ideal glass transition signals a crossover to activated dynamics, which is estimated to be observable before the nematic phase boundary is encountered if the aspect ratio is less than roughly 25. Calculations of the glassy elastic shear modulus and absolute yield stress reveal a roughly exponential growth with volume fraction. The dependence of entropic barriers and mean barrier hopping times on concentration for rods of variable aspect ratios can be collapsed quite well based on a difference volume fraction variable that quantifies the distance from the ideal glass boundary. Full numerical solution of the NLE theory via stochastic trajectory simulation was performed for tangent bead rods, and the results were compared to their hard sphere analogs. With increasing shape anisotropy the characteristic length scales of the nonequilibrium free energy increase and the magnitude of the localization well and entropic barrier curvatures decreases. These changes result in a significant aspect ratio dependence of dynamical properties and time correlation functions including weaker intermediate time subdiffusive transport, stronger two-step decay of the incoherent dynamic structure factor, longer mean alpha relaxation time, and stronger wavevector-dependent decoupling of relaxation times and the self-diffusion constant. The theoretical results are potentially testable via computer simulation, confocal microscopy, and dynamic light scattering.

Large-amplitude jumps and non-Gaussian dynamics in highly concentrated hard sphere fluids

Our microscopic stochastic nonlinear Langevin equation theory of activated dynamics has been employed to study the real-space van Hove function of dense hard sphere fluids and suspensions. At very short times, the van Hove function is a narrow Gaussian. At sufficiently high volume fractions, such that the entropic barrier to relaxation is greater than the thermal energy, its functional form evolves with time to include a rapidly decaying component at small displacements and a long-range exponential tail. The "jump" or decay length scale associated with the tail increases with time (or particle root-mean-square displacement) at fixed volume fraction, and with volume fraction at the mean alpha relaxation time. The jump length at the alpha relaxation time is predicted to be proportional to a measure of the decoupling of self-diffusion and structural relaxation. At long times corresponding to mean displacements of order a particle diameter, the volume fraction dependence of the decay length disappears. A good superposition of the exponential tail feature based on the jump length as a scaling variable is predicted at high volume fractions. Overall, the theoretical results are in good accord with recent simulations and experiments. The basic aspects of the theory are also compared with a classic jump model and a dynamically facilitated continuous time random-walk model. Decoupling of the time scales of different parts of the relaxation process predicted by the theory is qualitatively similar to facilitated dynamics models based on the concept of persistence and exchange times if the elementary event is assumed to be associated with transport on a length scale significantly smaller than the particle size.

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.

Fickian crossover and length scales from two point functions in supercooled liquids

Particle motion of a Lennard-Jones supercooled liquid near the glass transition is studied by molecular dynamics simulations. We analyze the wave vector dependence of relaxation times in the incoherent self-scattering function and show that at least three different regimes can be identified and its scaling properties determined. The transition from one regime to another happens at characteristic length scales. The length scale associated with the onset of Fickian diffusion corresponds to the maximum size of heterogeneities in the system, and the characteristic time scale is several times larger than the alpha relaxation time. A second crossover length scale is observed, which corresponds to the typical time and length of heterogeneities, in agreement with results from four point functions. The different regimes can be traced back to the behavior of the van Hove distribution of displacements, which shows a characteristic exponential regime in the heterogeneous region before the crossover to Gaussian diffusion and should be observable in experiments. Our results show that it is possible to obtain characteristic length scales of heterogeneities through the computation of two point functions at different times.

Characterization of the potential energy landscape of an antiplasticized polymer

The nature of the individual transitions on the potential energy landscape (PEL) associated with particle motion are directly examined for model fragile glass-forming polymer melts, and the results are compared to those of an antiplasticized polymer system. In previous work, we established that the addition of antiplasticizer reduces the fragility of glass formation so that the antiplasticized material is a stronger glass former. In the present work, we find that the antiplasticizing molecules reduce the energy barriers for relaxation compared to the pure polymer, implying that the antiplasticized system has smaller barriers to overcome in order to explore its configuration space. We examine the cooperativity of segmental motion in these bulk fluids and find that more extensive stringlike collective motion enables the system to overcome larger potential energy barriers, in qualitative agreement with both the Stillinger-Weber and Adam-Gibbs views of glass formation. Notably, the stringlike collective motion identified by our PEL analysis corresponds to incremental displacements that occur within larger-scale stringlike particle displacement processes associated with PEL metabasin transitions that mediate structural relaxation. These "substrings" nonetheless seem to exhibit changes in relative size with antiplasticization similar to those observed in "superstrings" that arise at elevated temperatures. We also study the effects of confinement on the energy barriers in each system. Film confinement makes the energy barriers substantially smaller in the pure polymer, while it has little effect on the energy barriers in the antiplasticized system. This observation is qualitatively consistent with our previous studies of stringlike motion in these fluids at higher temperatures and with recent experimental measurements by Torkelson and co-workers.

Frequency-dependent Stokes-Einstein relation in supercooled liquids

We investigate by molecular dynamics simulations the validity of the frequency-dependent Stokes-Einstein (SE) relation in supercooled liquids at different temperatures. The results indicate that the SE relation holds at intermediate frequencies that correspond to the beta -relaxation and the onset of the alpha -relaxation regimes. Large deviations, which increase as the temperature decreases, are observed at frequencies well below the frequency at which the non-Gaussian parameter alpha2 is maximum. We argue that the breakdown of the SE relation in supercooled liquids arises from underestimation of the diffusion coefficient due to neglect of correlated motions.

Probe particles alter dynamic heterogeneities in simple supercooled systems

The authors present results from molecular dynamics simulations on the effect of smooth and rough probes on the dynamics of a supercooled Lennard-Jones (LJ) mixture. The probe diameter was systematically varied from one to seven times the diameter of the large particles of the LJ mixture. Mean square displacements show that in the presence of a large smooth probe the supercooled liquid speeds up, while in the presence of a large rough probe, the supercooled liquid slows down. Non-Gaussian parameters indicate that with both smooth and rough probes, the heterogeneity of the supercooled system increases. From the analysis of local Debye-Waller factors, it is evident that the change in the dynamics of the LJ system is heterogeneous, with the largest perturbations close to the probes. Large smooth and rough probes appear to set up heterogeneities in these supercooled systems that would otherwise not occur, and these heterogeneities persist for long times.

Ideal glass transitions, shear modulus, activated dynamics, and yielding in fluids of nonspherical objects

An extension of naive ideal mode coupling theory (MCT) and its generalization to treat activated barrier hopping and glassy dynamics in fluids and suspensions composed of nonspherical hard core objects is proposed. An effective center-of-mass description is adopted. It corresponds to a specific type of pre-averaging of the dynamical consequences of orientational degrees of freedom. The simplest case of particles composed of symmetry-equivalent interaction sites is considered. The theory is implemented for a homonuclear diatomic shape of variable bond length. The naive MCT glass transition boundary is predicted to be a nonmonotonic function of the length-to-width or aspect ratio and occurs at a nearly unique value of the dimensionless compressibility. The latter quantifies the amplitude of long wavelength thermal density fluctuations, thereby (empirically) suggesting a tight connection between the onset of localization and thermodynamics. Localization lengths and elastic shear moduli for different aspect ratio and volume fraction systems approximately collapse onto master curves based on a reduced volume fraction variable that quantifies the distance from the ideal glass transition. Calculations of the entropic barrier height and hopping time, maximum restoring force, and absolute yield stress and strain as a function of diatomic aspect ratio and volume fraction have been performed. Strong correlations of these properties with the dimensionless compressibility are also found, and nearly universal dependences have been numerically identified based on property-specific nondimensionalizations. Generalization of the approach to rigid rods, disks, and variable shaped molecules is possible, including oriented liquid crystalline phases.

Non-Gaussian effects, space-time decoupling, and mobility bifurcation in glassy hard-sphere fluids and suspensions

Brownian trajectory simulation methods are employed to fully establish the non-Gaussian fluctuation effects predicted by our nonlinear Langevin equation theory of single particle activated dynamics in glassy hard-sphere fluids. The consequences of stochastic mobility fluctuations associated with the space-time complexities of the transient localization and barrier hopping processes have been determined. The incoherent dynamic structure factor was computed for a range of wave vectors and becomes of an increasingly non-Gaussian form for volume fractions beyond the (naive) ideal mode coupling theory (MCT) transition. The non-Gaussian parameter (NGP) amplitude increases markedly with volume fraction and is well described by a power law in the maximum restoring force of the nonequilibrium free energy profile. The time scale associated with the NGP peak becomes much smaller than the alpha relaxation time for systems characterized by significant entropic barriers. An alternate non-Gaussian parameter that probes the long time alpha relaxation process displays a different shape, peak intensity, and time scale of its maximum. However, a strong correspondence between the classic and alternate NGP amplitudes is predicted which suggests a deep connection between the early and final stages of cage escape. Strong space-time decoupling emerges at high volume fractions as indicated by a nondiffusive wave vector dependence of the relaxation time and growth of the translation-relaxation decoupling parameter. Displacement distributions exhibit non-Gaussian behavior at intermediate times, evolving into a strongly bimodal form with slow and fast subpopulations at high volume fractions. Qualitative and semiquantitative comparisons of the theoretical results with colloid experiments, ideal MCT, and multiple simulation studies are presented.

Relaxation scenarios in a mixture of large and small spheres: Dependence on the size disparity

We present a computational investigation on the slow dynamics of a mixture of large and small soft spheres. By varying the size disparity at a moderate fixed composition different relaxation scenarios are observed for the small particles. For small disparity density-density correlators exhibit moderate stretching. Only small quantitative differences are observed between dynamic features for large and small particles. On the contrary, large disparity induces a clear time scale separation between the large and small particles. Density-density correlators for the small particles become extremely stretched and display logarithmic relaxation by properly tuning the temperature or the wave vector. Self-correlators decay much faster than density-density correlators. For very large size disparity, a complete separation between self- and collective dynamics is observed for the small particles. Self-correlators decay to zero at temperatures where density-density correlations are frozen. The dynamic picture obtained by varying the size disparity resembles features associated with mode coupling transition lines of the types B and A at, respectively, small and very large size disparities. Both lines might merge, at some intermediate disparity, at a higher-order point, to which logarithmic relaxation would be associated. This picture resembles predictions of a recent mode coupling theory for fluids confined in matrices with interconnected voids [V. Krakoviack, Phys. Rev. Lett. 94, 065703 (2005)].

Cooperativity beyond caging: Generalized mode-coupling theory

The validity of mode-coupling theory (MCT) is restricted by an uncontrolled factorization approximation of density correlations. The factorization can be delayed and ultimately avoided, however, by explicitly including higher order correlations. We explore this approach within a microscopically motivated schematic model. Analytic tractability allows us to discuss in great detail the impact of factorization at arbitrary order, including the limit of avoided factorization. Our results indicate a coherent picture for the capabilities as well as limitations of MCT. Moreover, including higher order correlations systematically defers the transition and ultimately restores ergodicity. Power-law divergence of the relaxation time is then replaced by continuous but exponential growth.

Activated hopping and dynamical fluctuation effects in hard sphere suspensions and fluids

Single particle Brownian dynamics simulation methods are employed to establish the full trajectory level predictions of our nonlinear stochastic Langevin equation theory of activated hopping dynamics in glassy hard sphere suspensions and fluids. The consequences of thermal noise driven mobility fluctuations associated with the barrier hopping process are determined for various ensemble-averaged properties and their distributions. The predicted mean square displacements show classic signatures of transient trapping and anomalous diffusion on intermediate time and length scales. A crossover to a stronger volume fraction dependence of the apparent nondiffusive exponent occurs when the entropic barrier is of order the thermal energy. The volume fraction dependences of various mean relaxation times and rates can be fitted by empirical critical power laws with parameters consistent with ideal mode-coupling theory. However, the results of our divergence-free theory are largely a consequence of activated dynamics. The experimentally measurable alpha relaxation time is found to be very similar to the theoretically defined mean reaction time for escape from the barrier-dominated regime. Various measures of decoupling have been studied. For fluid states with small or nonexistent barriers, relaxation times obey a simple log-normal distribution, while for high volume fractions the relaxation time distributions become Poissonian. The product of the self-diffusion constant and mean alpha relaxation time increases roughly as a logarithmic function of the alpha relaxation time. The cage scale incoherent dynamic structure factor exhibits nonexponential decay with a modest degree of stretching. A nearly universal collapse of the different volume fraction results occurs if time is scaled by the mean alpha relaxation time. Hence, time-volume fraction superposition holds quite well, despite the presence of stretching and volume fraction dependent decoupling associated with the stochastic barrier hopping process. The relevance of other origins of dynamic heterogeneity (e.g., mesoscopic domains), and comparison of our results with experiments, simulations, and alternative theories, is discussed.

Hybrid Monte Carlo simulation of a glass-forming binary mixture

We propose a different, hybrid Monte Carlo algorithm that combines configurational bias particle swaps with parallel tempering. We use this algorithm to simulate a standard model of a glass-forming binary mixture above and below the so-called mode-coupling temperature T(MCT). We find that an ansatz that was used previously to extrapolate thermodynamic quantities to temperatures below T(MCT) breaks down in the vicinity of the mode-coupling temperature. We also find that the Adam-Gibbs relations D(alpha)exp(-a/TS(c)) and tau(alpha)exp(b/TS(c)), which connect the diffusion coefficient D and the relaxation time tau with the configurational entropy S(c), are valid for all temperatures for which the configurational and vibrational contributions to the free energy decouple.

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.

Derivation of a microscopic theory of barriers and activated hopping transport in glassy liquids and suspensions

A recently proposed microscopic activated barrier hopping theory [K. S. Schweizer and E. J. Saltzman, J. Chem. Phys. 119, 1181 (2003)] of slow single-particle dynamics in glassy liquids, suspensions, and gels is derived using nonequilibrium statistical mechanics. Fundamental elements underlying the stochastic nonlinear Langevin equation description include an inhomogeneous liquid or locally solid-state perspective, dynamic density-functional theory (DDFT), a local equilibrium closure, and a coarse-grained free-energy functional. A dynamic Gaussian approximation is not adopted which is the key for avoiding a kinetic ideal glass transition. The relevant excess free energy is of a nonequilibrium origin and is related to dynamic force correlations in the fluid. The simplicity of the approach allows external perturbations to be rather easily incorporated. Dynamic heterogeneity enters naturally via mobility fluctuations associated with the stochastic barrier-hopping process. The derivation both identifies the limitations of the theory and suggests new avenues for its systematic improvement. Comparisons with ideal mode-coupling theory, alternative DDFT approaches and a field theoretic path-integral formulation are presented.

Relaxation in a glassy binary mixture: comparison of the mode-coupling theory to a Brownian dynamics simulation

We solved the mode-coupling equations for the Kob-Andersen binary mixture using structure factors calculated from Brownian dynamics simulations of the same system. We found, as was previously observed, that the mode-coupling temperature T(c) inferred from simulations is about two times greater than that predicted by the theory. However, we find that many time-dependent quantities agree reasonably well with the predictions of the mode-coupling theory if they are compared at the same reduced temperature epsilon = (T - T(c))/T(c), and if epsilon is not too small. Specifically, the simulation results for the incoherent intermediate scattering function, the mean square displacement, the relaxation time, and the self-diffusion coefficient agree reasonably well with the predictions of the mode-coupling theory. We find that there are substantial differences for the non-Gaussian parameter. At small reduced temperatures the probabilities of the logarithm of single particle displacements demonstrate that there is hopping-like motion present in the simulations, and this motion is not predicted by the mode-coupling theory. The wave-vector-dependent relaxation time is shown to be qualitatively different from the predictions of the mode-coupling theory for temperatures where hopping-like motion is present.