A heterogeneous picture of α relaxation for fragile supercooled liquids

Glass spectrum, excess wing phenomenon, and master curves in molecular glass formers: A multi-method approach

The relaxation spectra of glass formers solely displaying an α-peak and excess wing contribution collected by various methods are reanalyzed to pin down their different spectral evolution. We show that master curve construction encompassing both α-peak and emerging excess wing works for depolarized light scattering (DLS) and nuclear magnetic resonance (NMR) relaxometry. It reveals the self-part of the slow dynamics' spectrum. Master curves are to be understood as a result of a more extensive scaling covering all temperatures instead of strict frequency-temperature superposition. DLS and NMR display identical relaxation spectra; yet, comparing different systems, we do not find a generic structural relaxation at variance with recent claims. Dielectric spectroscopy (DS) spectra show particularities, which render master curve construction obsolete. The DS α-peak is enhanced or suppressed with respect to that of DLS or NMR, yet, not correlated to the polarity of the liquid. Attempting to single out the excess wing from the overall spectrum discloses a stronger exponential temperature dependence of its amplitude compared to that below Tg and a link between its exponent and that of the fast dynamics' spectrum. Yet, such a decomposition of α-peak and excess wing appears to be unphysical. Among many different glasses, the amplitude of the excess wing power-law spectrum is found to be identical at Tg, interpreted as a relaxation analog to the Lindemann criterion.

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

A computational approach by an implementation of the principle component analysis (PCA) withK-means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers (WCNs‾) of particles calculated from particles' positions are used as an order parameter to build a low-dimensional representation of feature (structural) space forK-means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states' particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid-liquid phase separation after quenching to a supercooled state.

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

Excess wings and asymmetric relaxation spectra in a facilitated trap model

In a recent computer study, we have shown that the combination of spatially heterogeneous dynamics and kinetic facilitation provides a microscopic explanation for the emergence of excess wings in deeply supercooled liquids. Motivated by these findings, we construct a minimal empirical model to describe this physics and introduce dynamic facilitation in the trap model, which was initially developed to capture the thermally activated dynamics of glassy systems. We fully characterize the relaxation dynamics of this facilitated trap model varying the functional form of energy distributions and the strength of dynamic facilitation, combining numerical results and analytic arguments. Dynamic facilitation generically accelerates the relaxation of the deepest traps, thus making relaxation spectra strongly asymmetric, with an apparent "excess" signal at high frequencies. For well-chosen values of the parameters, the obtained spectra mimic experimental results for organic liquids displaying an excess wing. Overall, our results identify the minimal physical ingredients needed to describe excess processes in the relaxation spectra of supercooled liquids.

How to "measure" a structural relaxation time that is too long to be measured?

It has recently become possible to prepare ultrastable glassy materials characterized by structural relaxation times, which vastly exceed the duration of any feasible experiment. Similarly, new algorithms have led to the production of ultrastable computer glasses. Is it possible to obtain a reliable estimate of a structural relaxation time that is too long to be measured? We review, organize, and critically discuss various methods to estimate very long relaxation times. We also perform computer simulations of three dimensional ultrastable hard spheres glasses to test and quantitatively compare some of these methods for a single model system. The various estimation methods disagree significantly, and non-linear and non-equilibrium methods lead to a strong underestimate of the actual relaxation time. It is not yet clear how to accurately estimate extremely long relaxation times.

Effect of temperature and density fluctuations on the spatially heterogeneous dynamics of glass-forming Van der Waals liquids under high pressure

In this Letter, we show how temperature and density fluctuations affect the spatially heterogeneous dynamics at ambient and elevated pressures. By using high-pressure experimental data for van der Waals liquids, we examine contributions of the temperature and density fluctuations to the dynamics heterogeneity. We show that the dynamic heterogeneity decreases significantly with increasing pressure at a constant structural relaxation time (isochronal condition), while the broadening of the relaxation spectrum remains constant. This observation questions the relationship between spectral broadening and dynamic heterogeneity.

Spatial correlations in the dynamics of glassforming liquids: experimental determination of their temperature dependence

We use recently introduced three-point dynamic susceptibilities to obtain an experimental determination of the temperature evolution of the number of molecules Ncorr that are dynamically correlated during the structural relaxation of supercooled liquids. We first discuss in detail the physical content of three-point functions that relate the sensitivity of the averaged two-time dynamics to external control parameters (such as temperature or density), as well as their connection to the more standard four-point dynamic susceptibility associated with dynamical heterogeneities. We then demonstrate that these functions can be experimentally determined with good precision. We gather available data to obtain the temperature dependence of Ncorr for a large number of supercooled liquids over a wide range of relaxation time scales from the glass transition up to the onset of slow dynamics. We find that Ncorr systematically grows when approaching the glass transition. It does so in a modest manner close to the glass transition, which is consistent with an activation-based picture of the dynamics in glassforming materials. For higher temperatures, there appears to be a regime where Ncorr behaves as a power-law of the relaxation time. Finally, we find that the dynamic response to density, while being smaller than the dynamic response to temperature, behaves similarly, in agreement with theoretical expectations.

Non-Debye response for the structural relaxation in glass-forming liquids: Test of the Avramov model

The experimentally observed characteristic features of the alpha-relaxation process in glass-forming liquids are the non-Arrhenius behavior of the structural relaxation times and the non-Debye character of the macroscopic relaxation function. The Avramov model in which relaxation is considered as an energy activation process of surmounting random barriers in liquid energy landscape was successfully applied to describe the temperature and pressure dependences of the macroscopic relaxation times or viscosity. In this paper, we consider the dielectric spectrum associated with Avramov model. The asymmetrical broadening of the loss spectra was found to be related directly to dispersion of the energy barrier distribution. However, it turns out that temperature dependence of the spectrum broadening as predicted by the Avromov model is at odds to experimental observation in glass-forming liquids.

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.

Spontaneous and induced dynamic correlations in glass formers. II. Model calculations and comparison to numerical simulations

We study in detail the predictions of various theoretical approaches, in particular, mode-coupling theory (MCT) and kinetically constrained models (KCMs), concerning the time, temperature, and wave vector dependence of multipoint correlation functions that quantify the strength of both induced and spontaneous dynamical fluctuations. We also discuss the precise predictions of MCT concerning the statistical ensemble and microscopic dynamics dependence of these multipoint correlation functions. These predictions are compared to simulations of model fragile and strong glass-forming liquids. Overall, MCT fares quite well in the fragile case, in particular, explaining the observed crucial role of the statistical ensemble and microscopic dynamics, while MCT predictions do not seem to hold in the strong case. KCMs provide a simplified framework for understanding how these multipoint correlation functions may encode dynamic correlations in glassy materials. However, our analysis highlights important unresolved questions concerning the application of KCMs to supercooled liquids.

Lamellar order, microphase structures, and glassy phase in a field theoretic model for charged colloids

In this paper we present a detailed analytical study of the phase diagram and of the structural properties of a field theoretic model with a short-range attraction and a competing long-range screened repulsion. We provide a full derivation and expanded discussion and digression on results previously reported briefly in M. Tarzia and A. Coniglio, Phys. Rev. Lett. 96, 075702 (2006). The model contains the essential features of the effective interaction potential among charged colloids in polymeric solutions. We employ the self-consistent Hartree approximation and a replica approach, and we show that varying the parameters of the repulsive potential and the temperature yields a phase coexistence, a lamellar and a glassy phase. Our results suggest that the cluster phase observed in charged colloids might be the signature of an underlying equilibrium lamellar phase, hidden on experimental time scales, and emphasize that the formation of microphase structures may play a prominent role in the process of colloidal gelation.

Breakdown of the Stokes−Einstein Relationship: Role of Interactions in the Size Dependence of Self-Diffusivity

Einstein and others derived the reciprocal dependence of the self-diffusivity D on the solute radius r(u) for large solutes based on kinetic theory. We examine here (a) the range of r(u) over which Stokes-Einstein (SE) dependence is valid and (b) the precise dependence for small solutes outside of the SE regime. We show through molecular dynamics simulations that there are two distinct regimes for smaller solutes: (i) the interaction or Levitation effect (LE) regime for solutes of intermediate size and (ii) the D proportional, variant 1/r(u)(2) for still smaller solutes. We show that as the solute-solvent size ratio decreases, the breakdown in the Stokes-Einstein relationship leading to the LE regime has its origin in dispersion interaction between the solute and the solvent. These results explain reports of enhanced solute diffusion in solvents existing in the literature seen for small solutes for which no explanation exists.

Nature of the breakdown in the Stokes-Einstein relationship in a hard sphere fluid

Molecular dynamics simulations of high density hard sphere fluids clearly show a breakdown of the Stokes-Einstein equation (SE). This result has been conjectured to be due to the presence of mobile particles, i.e., ones which have the propensity to "hop" distances that are integer multiples of the interparticle distance. We conclusively show that the sedentary particles, i.e., ones complementary to the "hoppers," obey the SE relationship to a good approximation, even though the fluid as a whole violates the SE equation at high densities. These results support the notion that the unusual diffusive behavior of supercooled liquids is dominated by the hopping particles.

Pattern formation and glassy phase in the phi4 theory with a screened electrostatic repulsion

We study analytically the structural properties of a system with a short-range attraction and a competing long-range screened repulsion. This model contains the essential features of the effective interaction potential among charged colloids in polymeric solutions and provides novel insights on the equilibrium phase diagram of these systems. Within the self-consistent Hartree approximation and by using a replica approach, we show that varying the parameters of the repulsive potential and the temperature yields a phase coexistence, a lamellar, and a glassy phase. Our results strongly suggest that the cluster phase observed in charged colloids might be the signature of an underlying equilibrium lamellar phase, hidden on experimental time scales.

Direct Experimental Evidence of a Growing Length Scale Accompanying the Glass Transition

Understanding glass formation is a challenge, because the existence of a true glass state, distinct from liquid and solid, remains elusive: Glasses are liquids that have become too viscous to flow. An old idea, as yet unproven experimentally, is that the dynamics becomes sluggish as the glass transition approaches, because increasingly larger regions of the material have to move simultaneously to allow flow. We introduce new multipoint dynamical susceptibilities to estimate quantitatively the size of these regions and provide direct experimental evidence that the glass formation of molecular liquids and colloidal suspensions is accompanied by growing dynamic correlation length scales.

Numerical Study of a Fragile Three-Dimensional Kinetically Constrained Model

We numerically study the three-dimensional generalization of the kinetically constrained east model, the north-or-east-or-front (NEF) model. We characterize the equilibrium behavior of the NEF model in detail, measuring the temperature dependence of several quantities: alpha-relaxation time, distributions of relaxation times, dynamic susceptibility, dynamic correlation length, and four-point susceptibility. We show that the NEF model describes quantitatively experimental observations over an exceptionally wide range of time scales. We illustrate this by fitting experimental data obtained both in the mildly supercooled regime by optical Kerr effect and close to the glass transition by dielectric spectroscopy.

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.

Coarse-grained microscopic model of glass formers

We introduce a coarse-grained model for atomic glass formers. Its elements are physically motivated local microscopic dynamical rules parameterized by observables. Results of the model are established and used to interpret the measured behaviors of supercooled fluids approaching glass transitions. The model predicts the presence of a crossover from hierarchical super-Arrhenius dynamics at short length scales to diffusive Arrhenius dynamics at large length scales. This prediction distinguishes our model from other theories of glass formers and can be tested by experiment.

Possible resolution of the Kauzmann paradox in supercooled liquids

Generally, the entropy of the supercooled liquid decreases more rapidly than that of the crystal. Thus, the former, if we extrapolate it smoothly below the glass-transition temperature T(g), becomes equal to the latter at the so-called Kauzmann temperature T(K). Further extrapolation below T(K) leads to the unphysical situation that the entropy of disordered liquid is lower than the ordered crystal, which results in the violation of the third law of thermodynamics. This is known as the "Kauzmann paradox" which has been the key problem of liquid-glass transition for a long time. Here we propose a simple resolution of the Kauzmann paradox by answering a fundamental question of how deeply we can supercool a liquid. We argue that the lower metastable limit T(LML), below which a liquid should crystallize before its structural relaxation, is located above the Kauzmann temperature T(K). Thus, the entropy crisis at T(K) is naturally avoided by crystallization. We suggest that it is dynamic heterogeneity that destabilizes a deeply supercooled "equilibrium" liquid state as well as a glassy state against crystallization. This may have a significant implication on the stability of a glassy state, which is of industrial importance in relation to the storage of glassy material.

Langevin dynamics of the Coulomb frustrated ferromagnet: A mode-coupling analysis

We study the Langevin dynamics of the soft-spin, continuum version of the Coulomb-frustrated Ising ferromagnet. By using the dynamical mode-coupling approximation, supplemented by reasonable approximations for describing the equilibrium static correlation function, and the somewhat improved dynamical self-consistent screening approximation, we find that the system displays a transition from an ergodic to a nonergodic behavior. This transition is similar to that obtained in the idealized mode-coupling theory of glass-forming liquids and in the mean-field generalized spin glasses with one-step replica symmetry breaking. The significance of this result and the relation to the appearance of a complex free-energy landscape are also discussed.

Potential energy landscape and mechanisms of diffusion in liquids

The mechanism of diffusion in supercooled liquids is investigated from the potential energy landscape point of view, with emphasis on the crossover from high- to low-T dynamics over the range T(A) > or =T > or =T(c). Molecular dynamics simulations with a time dependent mapping to the associated local minimum or inherent structure (IS) are performed on unit-density Lennard-Jones. Dynamical quantities introduced include r2(is)(t), the mean-square displacement (MSD) within a basin of attraction of an IS, R2(t), the MSD of the IS itself, and g(t), the distribution of IS waiting times. The configuration space is treated as a composite of the contributions of cooperative local regions, and a method is given to obtain the physically meaningful g(loc)(t) and mean waiting time tau(loc) from g(t). An understanding of the crossover is obtained in terms of r2(is)(t) and tau(loc). At intermediate T, r2(is)(t) possesses an interval of linear t dependence allowing calculation of an intrabasin diffusion constant D(is). Near T(c), where intrabasin diffusion is well established for t<tau(loc), diffusion is intrabasin dominated with D=D(is); D may be calculated within a basin. Below T(c), tau(loc) exceeds the time tau(pl) needed for the system to explore the basin, indicating the action of barriers at the border; tau(loc)=tau(pl) is a criterion for transition to activated hopping. Intrabasin diffusion provides a means of confinement not involving barriers and plays a key role in the dynamics above T(c). The distinction is discussed between motion among the IS (IS dynamics) below T(c) and saddle or border dynamics above T(c), where the system is always close to one of the saddle barriers connecting the basins and IS boundaries are closely spaced and easily crossed. A border index is introduced based upon the relation of R2(t) to the conventional MSD, and shown to vanish at T approximately T(c). It is proposed that intrabasin diffusion is a manifestation of saddle dynamics.