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

Evidence of ferroelectric features in low-density supercooled water from ab initio deep neural-network simulations

Over the last decade, an increasing body of evidence has emerged, supporting the existence of a metastable liquid-liquid critical point in supercooled water whereby two distinct liquid phases of different densities coexist. Analyzing long molecular dynamics simulations performed using deep neural-network force fields trained to accurate quantum mechanical data, we demonstrate that the low-density liquid phase displays a strong propensity toward spontaneous polarization, as witnessed by large and long-lived collective dipole fluctuations. Our findings suggest that the dynamical stability of the low-density phase, and hence the transition from high-density to low-density liquid, is triggered by a collective process involving an accumulation of rotational angular jumps, which could ignite large dipole fluctuations. This dynamical transition involves subtle changes in the electronic polarizability of water molecules which affects their rotational mobility within the two phases. These findings hold the potential for catalyzing activity in the search for dielectric-based probes of the putative second critical point.

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

A method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob-Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density-density and hydrostatic pressure-pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. These results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

Predicting Dynamic Heterogeneity in Glass-Forming Liquids by Physics-Inspired Machine Learning

We introduce GlassMLP, a machine learning framework using physics-inspired structural input to predict the long-time dynamics in deeply supercooled liquids. We apply this deep neural network to atomistic models in 2D and 3D. Its performance is better than the state of the art while being more parsimonious in terms of training data and fitting parameters. GlassMLP quantitatively predicts four-point dynamic correlations and the geometry of dynamic heterogeneity. Transferability across system sizes allows us to efficiently probe the temperature evolution of spatial dynamic correlations, revealing a profound change with temperature in the geometry of rearranging regions.

Extending the analogy between intracellular motion in mammalian cells and glassy dynamics

How molecules, organelles, and foreign objects move within living cells has been studied in organisms ranging from bacteria to human cells. In mammalian cells, in particular, cellular vesicles move across the cell using motor proteins that carry the vesicle down the cytoskeleton to their destination. We have recently noted several similarities between the motion of such vesicles and that in disordered, "glassy", systems, but the generality of this observation remains unclear. Here we follow the motion of mitochondria, the organelles responsible for cell energy production, in mammalian cells over timescales from 50 ms to 70 s. Qualitative observations show that single mitochondria remain within a spatially limited region for extended periods of time, before moving longer distances relatively quickly. The displacement distribution is roughly Gaussian for shorter distances (≲0.05 μm) but exhibits exponentially decaying tails at longer distances (up to 0.40 μm). This behaviour is well-described by a model developed to describe the motion in glassy systems. These observations are extended to in total 3 different objects (mitochondria, lysosomes and nano-sized beads enclosed in vesicles), 3 different mammalian cell types (HEK 293, HeLa, and HT22), from 2 different organisms (human and mouse). Further evidence that supports glass-like characteristics of the motion is a difference between the time it takes to move a longer distance for the first time and subsequent times, as well as a weak ergodicity breaking of the motion. Overall, we demonstrate the ubiquity of glass-like motion in mammalian cells, providing a different perspective on intracellular motion.

Choice of diameters in a polydisperse model glassformer: Deterministic or stochastic?

In particle-based computer simulations of polydisperse glassforming systems, the particle diameters σ=σ_{1},⋯,σ_{N} of a system with N particles are chosen with the intention to approximate a desired distribution density f with the corresponding histogram. One method to accomplish this is to draw each diameter randomly and independently from the density f. We refer to this stochastic scheme as model S. Alternatively, one can apply a deterministic method, assigning an appropriate set of N values to the diameters. We refer to this method as model D. We show that, for sample-to-sample fluctuations, especially for the glassy dynamics at low temperatures, it matters whether one chooses model S or model D. Using molecular dynamics computer simulations, we investigate a three-dimensional polydisperse nonadditive soft-sphere system with f(s)∼s^{-3}. The swap Monte Carlo method is employed to obtain equilibrated samples at very low temperatures. We show that for model S the sample-to-sample fluctuations due to the quenched disorder imposed by the diameters σ can be explained by an effective packing fraction. Dynamic susceptibilities in model S can be split into two terms: one that is of thermal nature and can be identified with the susceptibility of model D, and another one originating from the disorder in σ. At low temperatures the latter contribution is the dominating term in the dynamic susceptibility. Our study clarifies the pros and cons of the use of models S and D in practice.

Dynamical susceptibilities near ideal glass transitions

Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalized mode-coupling theory of supercooled liquids to inhomogeneous environments. This provides a first-principles-based, systematic, and rigorous way of deriving high-point dynamical susceptibilities from variations of the many-body dynamic structure factors with respect to their conjugate field. This framework allows for a fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition. The behavior of these dynamical susceptibilities is then studied in the context of simplified self-consistent relaxation models.

Beyond the Average: Spatial and Temporal Fluctuations in Oxide Glass-Forming Systems

Atomic structure dictates the performance of all materials systems; the characteristic of disordered materials is the significance of spatial and temporal fluctuations on composition-structure-property-performance relationships. Glass has a disordered atomic arrangement, which induces localized distributions in physical properties that are conventionally defined by average values. Quantifying these statistical distributions (including variances, fluctuations, and heterogeneities) is necessary to describe the complexity of glass-forming systems. Only recently have rigorous theories been developed to predict heterogeneities to manipulate and optimize glass properties. This article provides a comprehensive review of experimental, computational, and theoretical approaches to characterize and demonstrate the effects of short-, medium-, and long-range statistical fluctuations on physical properties (e.g., thermodynamic, kinetic, mechanical, and optical) and processes (e.g., relaxation, crystallization, and phase separation), focusing primarily on commercially relevant oxide glasses. Rigorous investigations of fluctuations enable researchers to improve the fundamental understanding of the chemistry and physics governing glass-forming systems and optimize structure-property-performance relationships for next-generation technological applications of glass, including damage-resistant electronic displays, safer pharmaceutical vials to store and transport vaccines, and lower-attenuation fiber optics. We invite the reader to join us in exploring what can be discovered by going beyond the average.

Revisiting the single-saddle model for the β-relaxation of supercooled liquids

The dynamics of glass-forming liquids display several outstanding features, such as two-step relaxation and dynamic heterogeneities, which are difficult to predict quantitatively from first principles. In this work, we revisit a simple theoretical model of the β-relaxation, i.e., the first step of the relaxation dynamics. The model, first introduced by Cavagna et al. [J. Phys. A: Math. Gen. 36, 10721 (2003)], describes the dynamics of the system in the neighborhood of a saddle point of the potential energy surface. We extend the model to account for density–density correlation functions and for the four-point dynamic susceptibility. We obtain analytical results for a simple schematic model, making contact with related results for p-spin models and with the predictions of inhomogeneous mode-coupling theory. Building on recent computational advances, we also explicitly compare the model predictions against overdamped Langevin dynamics simulations of a glass-forming liquid close to the mode-coupling crossover. The agreement is quantitative at the level of single-particle dynamic properties only up to the early β-regime. Due to its inherent harmonic approximation, however, the model is unable to predict the dynamics on the time scale relevant for structural relaxation. Nonetheless, our analysis suggests that the agreement with the simulations may be largely improved if the modes’ spatial localization is properly taken into account.

Glass-like characteristics of intracellular motion in human cells

The motion in the cytosol of microorganisms such as bacteria and yeast has been observed to undergo a dramatic slowing down upon cell energy depletion. These observations have been interpreted as the motion being "glassy," but whether this notion is useful also for active, motor-protein-driven transport in eukaryotic cells is less clear. Here, we use fluorescence microscopy of beads in human (HeLa) cells to probe the motion of membrane-surrounded structures that are carried along the cytoskeleton by motor proteins. Evaluating several hallmarks of glassy dynamics, we show that at short length scales, the motion is heterogeneous, is nonergodic, is well described by a model for the displacement distribution in glassy systems, and exhibits a decoupling of the exchange and persistence times. Overall, these results suggest that the short length scale behavior of objects that can be transported actively by motor proteins in human cells shares features with the motion in glassy systems.

Random First Order Transition Theory for Glassy Dynamics in a Single Condensed Polymer

The number of compact structures of a single condensed polymer (SCP), with similar free energies, grows exponentially with the degree of polymerization. In analogy with structural glasses (SGs), we expect that at low temperatures chain relaxation should occur by activated transitions between the compact metastable states. By evolving the states of the SCP, linearly coupled to a reference state, we show that, below a dynamical transition temperature (T_{d}), the SCP is trapped in a metastable state leading to slow dynamics. At a lower temperature, T_{K}≠0, the configurational entropy vanishes, resulting in a thermodynamic random first order ideal glass transition. The relaxation time obeys the Vogel-Fulcher-Tamman law, diverging at T=T_{0}≈T_{K}. These findings, accord well with the random first order transition theory, establishing that SCP and SG exhibit similar universal characteristics.

Connecting real glasses to mean-field models

We propose a novel model for a glass-forming liquid, which allows us to switch in a continuous manner from a standard three-dimensional liquid to a fully connected mean-field model. This is achieved by introducing k additional particle-particle interactions, which thus augments the effective number of neighbors of each particle. Our computer simulations of this system show that the structure of the liquid does not change with the introduction of these pseudo-neighbors and by means of analytical calculations, and we determine the structural properties related to these additional neighbors. We show that the relaxation dynamics of the system slows down very quickly with the increase in k and that the onset and the mode-coupling temperatures increase. The systems with high values of k follow the mode-coupling theory power law behavior for a larger temperature range compared to the ones with lower values of k. The dynamic susceptibility indicates that the dynamic heterogeneity decreases with the increase in k, whereas the non-Gaussian parameter is independent of it. Thus, we conclude that with the increase in the number of pseudo-neighbors, the system becomes more mean-field-like. By comparing our results with previous studies on mean-field-like systems, we come to the conclusion that the details of how the mean-field limit is approached are important since they can lead to different dynamical behavior in this limit.

Connecting relaxation time to a dynamical length scale in athermal active glass formers

Supercooled liquids display dynamics that are inherently heterogeneous in space. This essentially means that at temperatures below the melting point, particle dynamics in certain regions of the liquid can be orders of magnitude faster than other regions. Often dubbed dynamical heterogeneity, this behavior has fascinated researchers involved in the study of glass transition for over two decades. A fundamentally important question in all glass transition studies is whether one can connect the growing relaxation time to a concomitantly growing length scale. In this paper, we go beyond the realm of ordinary glass forming liquids and study the origin of a growing dynamical length scale ξ in a self-propelled "active" glass former. This length scale, which is constructed using structural correlations, agrees well with the average size of the clusters of slow-moving particles that are formed as the liquid becomes spatially heterogeneous. We further report that the concomitantly growing α-relaxation time exhibits a simple scaling law, τ_{α}∼exp(μξ/T_{eff}), with μ as an effective chemical potential, T_{eff} as the effective temperature, and μξ as the growing free energy barrier for cluster rearrangements. The findings of our study are valid over four decades of persistence times, and hence they could be very useful in understanding the slow dynamics of a generic active liquid such as an active colloidal suspension, or a self-propelled granular medium.

Finite Dimensional Vestige of Spinodal Criticality above the Dynamical Glass Transition

Finite dimensional signatures of spinodal criticality are notoriously difficult to come by. The dynamical transition of glass-forming liquids, first described by mode-coupling theory, is a spinodal instability preempted by thermally activated processes that also limit how close the instability can be approached. We combine numerical tools to directly observe vestiges of the spinodal criticality in finite dimensional glass formers. We use the swap Monte Carlo algorithm to efficiently thermalize configurations beyond the mode-coupling crossover, and analyze their dynamics using a scheme to screen out activated processes, in spatial dimensions ranging from d=3 to d=10. We observe a strong softening of the mean-field square-root singularity in d=3 that is progressively restored as d increases above d=8, in surprisingly good agreement with perturbation theory.

Solvable Models of Supercooled Liquids in Three Dimensions

We introduce a supercooled liquid model and obtain parameter-free quantitative predictions that are in excellent agreement with numerical simulations, notably in the hard low-temperature region characterized by strong deviations from mode-coupling-theory behavior. The model is the Fredrickson-Andersen kinetically constrained model on the three-dimensional M-layer lattice. The agreement has implications beyond the specific model considered because the theory is potentially valid for many more systems, including realistic models and actual supercooled liquids.

An alternative approach to evidence the structural conditioning in the dynamic slowdown in a polymer glass-former

Dynamic slowdown of liquids, leading to a breakdown of Arrhenius behavior of relaxation and Stokes-Einstein relationship (SER), as the glass transition is approached, is still not fully understood despite decades of study. They are usually associated to the emergence of dynamic heterogeneity, that is, regions or clusters of particles that have high or low mobilities. But the physical origin of these dynamic heterogeneity, and in particular, the question whether they have a structural origin or they are a purely dynamical phenomenon, is still under debate. In this work we study through molecular dynamics simulations in a polymer model the dynamic slowdown and the breakdown of SER, in connection with dynamic susceptibility calculated for an isoconfigurational ensemble, such that the effects of structure on dynamics can be discriminated. The onset of structure effects on dynamical behavior is found to be coincident with the onset of slow dynamics and SER breakdown.

Common microscopic structural origin for water's thermodynamic and dynamic anomalies

Water displays a vast array of unique properties, known as water's anomalies, whose origin remains subject to hot debate. Our aim in this article is to provide a unified microscopic physical picture of water's anomalies in terms of locally favored structures, encompassing both thermodynamic and dynamic anomalies, which are often attributed to different origins. We first identify locally favored structures via a microscopic structural descriptor that measures local translational order and provide direct evidence that they have a hierarchical impact on the anomalies. At each state point, the strength of thermodynamic anomalies is directly proportional to the amount of locally favored structures, while the dynamic properties of each molecule depend on the local structure surrounding both itself and its nearest neighbors. To incorporate this, we develop a novel hierarchical two-state model. We show by extensive simulations of two popular water models that both thermodynamic and kinetic anomalies can be almost perfectly explained by the temperature and pressure dependence of these local and non-local versions of the same structural descriptor, respectively. Moreover, our scenario makes three unique predictions in supercooled water, setting it apart from other scenarios: (1) Presence of an "Arrhenius-to-Arrhenius" crossover upon cooling, as the origin of the apparent "fragile-to-strong" transition; (2) maximum of dynamic heterogeneity around 20 K below the Widom line and far above the glass transition; (3) violation of the Stokes-Einstein-Debye relation at ∼2T _g, rather than 1.2T _g typical of normal glass-formers. These predictions are verified by recent measurement of water's diffusion at very low temperatures (point 1) and discoveries from our extensive simulations (points 2-3). We suggest that the same scenario may generally apply to water-like anomalies in liquids tending to form locally favored structures, including not only other important tetrahedral liquids such as silicon, germanium, and silica, but also metallic and chalcogenide liquids.

Interplay between structure and relaxation in polyurea networks: the point of view from a novel method of cooperativity analysis of dielectric response

The influence of structural constraints on the relaxation dynamics of three polyurea networks with a varying degree of crosslinking, has been studied by means of a thorough analysis of broadband dielectric spectroscopy measurements. Two different relaxation processes are observed, namely, a fast process involving the soft poly(propylene oxide) chains, and a slower and much broader process associated with the immediate surroundings of the hard crosslinkers. Microphase separation in soft and hard domains characterizes the systems in the presence of hydrogen bonding. In this case, different confinement conditions are explored by varying the soft chain length; overall, so called "adsorption" effects dominate. With respect to both cooperativity and the rearrangement energy threshold in fast relaxation, it is found that the enhancement of configurational constraints is similar to cooling, but only on qualitative grounds. An upper bound of the hard domains' interface thickness, in which the slow relaxation is believed to take place, is estimated from the analysis of the fast relaxation in the system characterized by the highest degree of confinement, taking into account the results of the structural analysis. Dropping the hydrogen bonding mechanism, phase separation does not occur anymore and the configurational constraints at the ends of the soft chains are reduced, leaving just those imposed by the rigid crosslinkers. This leads to a significant increase in cooperativity on approaching the glass transition, and to a complex behavior that is thoroughly discussed in comparison with those observed in the micro-segregated systems.

The race to the bottom: approaching the ideal glass?

Key to resolving the scientific challenge of the glass transition is to understand the origin of the massive increase in viscosity of liquids cooled below their melting temperature (avoiding crystallisation). A number of competing and often mutually exclusive theoretical approaches have been advanced to describe this phenomenon. Some posit a bona fide thermodynamic phase to an 'ideal glass', an amorphous state with exceptionally low entropy. Other approaches are built around the concept of the glass transition as a primarily dynamic phenomenon. These fundamentally different interpretations give equally good descriptions of the data available, so it is hard to determine which-if any-is correct. Recently however this situation has begun to change. A consensus has emerged that one powerful means to resolve this longstanding question is to approach the putative thermodynamic transition sufficiently closely, and a number of techniques have emerged to meet this challenge. Here we review the results of some of these new techniques and discuss the implications for the existence-or otherwise-of the thermodynamic transition to an ideal glass.

Origin of the emergent fragile-to-strong transition in supercooled water

Liquids can be broadly classified into two categories, fragile and strong ones, depending on how their dynamical properties change with temperature. The dynamics of a strong liquid obey the Arrhenius law, whereas the fragile one displays a super-Arrhenius law, with a much steeper slowing down upon cooling. Recently, however, it was discovered that many materials such as water, oxides, and metals do not obey this simple classification, apparently exhibiting a fragile-to-strong transition far above [Formula: see text] Such a transition is particularly well known for water, and it is now regarded as one of water's most important anomalies. This phenomenon has been attributed to either an unusual glass transition behavior or the crossing of a Widom line emanating from a liquid-liquid critical point. Here by computer simulations of two popular water models and through analyses of experimental data, we show that the emergent fragile-to-strong transition is actually a crossover between two Arrhenius regimes with different activation energies, which can be naturally explained by a two-state description of the dynamics. Our finding provides insight into the fragile-to-strong transition observed in a wide class of materials.

Cooperativity in plastic crystals

A statistical mechanical model previously adopted for the analysis of the α-relaxation in structural glass formers is rederived within a general theoretical framework originally developed for systems approaching the ideal glassy state. The interplay between nonexponentiality and cooperativity is reconsidered in the light of energy landscape concepts. The method is used to estimate the cooperativity in orientationally disordered crystals, either from the analysis of literature data on linear dielectric response or from the enthalpy relaxation function obtained by temperature-modulated calorimetry. Knowledge of the specific heat step due to the freezing of the configurational or conformational modes at the glass transition is needed in order to properly account for the extent to which the relaxing system deviates from equilibrium during the rearrangement processes. A number of plastic crystals have been analyzed, and relatively higher cooperativities are found in the presence of hydrogen bonding interaction.

Dynamic and thermodynamic crossover scenarios in the Kob-Andersen mixture: Insights from multi-CPU and multi-GPU simulations

The physical behavior of glass-forming liquids presents complex features of both dynamic and thermodynamic nature. Some studies indicate the presence of thermodynamic anomalies and of crossovers in the dynamic properties, but their origin and degree of universality is difficult to assess. Moreover, conventional simulations are barely able to cover the range of temperatures at which these crossovers usually occur. To address these issues, we simulate the Kob-Andersen Lennard-Jones mixture using efficient protocols based on multi-CPU and multi-GPU parallel tempering. Our setup enables us to probe the thermodynamics and dynamics of the liquid at equilibrium well below the critical temperature of the mode-coupling theory, [Formula: see text]. We find that below [Formula: see text] the analysis is hampered by partial crystallization of the metastable liquid, which nucleates extended regions populated by large particles arranged in an fcc structure. By filtering out crystalline samples, we reveal that the specific heat grows in a regular manner down to [Formula: see text] . Possible thermodynamic anomalies suggested by previous studies can thus occur only in a region of the phase diagram where the system is highly metastable. Using the equilibrium configurations obtained from the parallel tempering simulations, we perform molecular dynamics and Monte Carlo simulations to probe the equilibrium dynamics down to [Formula: see text]. A temperature-derivative analysis of the relaxation time and diffusion data allows us to assess different dynamic scenarios around [Formula: see text]. Hints of a dynamic crossover come from analysis of the four-point dynamic susceptibility. Finally, we discuss possible future numerical strategies to clarify the nature of crossover phenomena in glass-forming liquids.

Growing length scale accompanying vitrification: A perspective based on nonsingular density fluctuations

A model for describing growing length scale accompanying the vitrification is introduced, in which we assume that in a subsystem whose density is above a certain threshold value, ρ_{c}, due to steric constraints, particle rearrangements are highly suppressed for a sufficiently long time period (∼structural relaxation time). We regard such a subsystem as a glassy cluster. With this assumption and without introducing any complicated thermodynamic arguments, we predict that with compression (increasing average density ρ) at a fixed temperature T in supercooled states, the characteristic length of the clusters, ξ, diverges as ξ∼(ρ_{c}-ρ)^{-2/d}, where d is the spatial dimensionality. This ξ measures the average persistence length of the steric constraints in blocking the rearrangement motions and is determined by the subsystem density. Additionally, with decreasing T at a fixed ρ, the length scale diverges in the same manner as ξ∼(T-T_{c})^{-2/d}, for which ρ is identical to ρ_{c} at T=T_{c}. The exponent describing the diverging length scale is the same as the one predicted by some theoretical models and indeed has been observed in some simulations and experiments. However, the basic mechanism for this divergence is different; that is, we do not invoke thermodynamic anomalies associated with the thermodynamic phase transition as the origin of the growing length scale. We further present arguements for the cooperative properties of the structural relaxation based on the clusters.

Unifying different interpretations of the nonlinear response in glass-forming liquids

This work aims at reconsidering several interpretations coexisting in the recent literature concerning nonlinear susceptibilities in supercooled liquids. We present experimental results on glycerol and propylene carbonate, showing that the three independent cubic susceptibilities have very similar frequency and temperature dependences, for both their amplitudes and phases. This strongly suggests a unique physical mechanism responsible for the growth of these nonlinear susceptibilities. We show that the framework proposed by two of us [J.-P. Bouchaud and G. Biroli, Phys. Rev. B 72, 064204 (2005)PRBMDO1098-012110.1103/PhysRevB.72.064204], where the growth of nonlinear susceptibilities is intimately related to the growth of glassy domains, accounts for all the salient experimental features. We then review several complementary and/or alternative models and show that the notion of cooperatively rearranging glassy domains is a key (implicit or explicit) ingredient to all of them. This paves the way for future experiments, which should deepen our understanding of glasses.

Study of dynamical heterogeneities in colloidal nanoclay suspensions approaching dynamical arrest

The dynamics of aqueous Laponite clay suspensions slow down with increasing sample waiting time (t w ). This behavior, and the material fragility that results, closely resemble the dynamical slowdown in fragile supercooled liquids with decreasing temperature, and are typically ascribed to the increasing sizes of distinct dynamical heterogeneities in the sample. In this article, we characterize the dynamical heterogeneities in Laponite suspensions by invoking the three-point dynamic susceptibility formalism. The average time-dependent two-point intensity autocorrelation and its sensitivity to t w are probed in dynamic light scattering experiments. Distributions of relaxation time scales, deduced from the Kohlrausch-Williams-Watts equation, are seen to widen with increasing t w . The calculated three-point dynamic susceptibility of Laponite suspensions exhibits a peak, with the peak height increasing with evolving t w at fixed volume fraction or with increasing volume fraction at fixed t w , thereby signifying the slowdown of the sample dynamics. The number of dynamically correlated particles, calculated from the peak-height, is seen to initially increase rapidly with increasing t w , before eventually slowing down close to the non-ergodic transition point. This observation is in agreement with published reports on supercooled liquids and hard sphere colloidal suspensions and offers a unique insight into the colloidal glass transition of Laponite suspensions.

Emergence of compositionally tunable nanoscale dynamical heterogeneity in model binary lipid biomembranes

While the existence of nanoscale dynamical heterogeneity in biological membranes has been suggested to act as an active functional platform for enabling various cellular processes like signal transduction and viral or bacterial entry, it has been extremely difficult to detect the existence of such domains. Model lipid bilayer membranes have been widely used to detect such dynamical heterogeneity in order to avoid complications arising from the compositional heterogeneity of cellular membranes. However, even in model biological membranes the issue of nanoscale lipid dynamics has remained controversial and unresolved due to the difficulty of detecting the existence of such dynamical heterogeneity on the scale of 10-300 nm. Here we report direct evidence of nanoscale lipid dynamical heterogeneity in model binary lipid bilayer membranes using a combination of super-resolution stimulated emission depletion (STED) microscopy and fluorescence correlation spectroscopy (FCS). We control the phase behavior of the lipid bilayers by varying their composition and discuss how this leads to the emergence of dynamical lipid domains on the scale of 80-150 nm, which is also dependent on the lipid phase in which such dynamics are observed. Notably, our work shows that the presence of cholesterol is not required for the existence of such domains even in fluid like bilayers, as has been widely believed, and specifies the minimal conditions required for the emergence of such dynamical heterogeneity in cellular membranes. Our work will thus not only be of great significance towards understanding the nanoscale dynamic organizing principles of cellular membranes but could also be useful in understanding the dynamics of related soft matter systems and nanoparticle-cell membrane interactions.

Disordered hyperuniform heterogeneous materials

Disordered hyperuniform many-body systems are distinguishable states of matter that lie between a crystal and liquid: they are like perfect crystals in the way they suppress large-scale density fluctuations and yet are like liquids or glasses in that they are statistically isotropic with no Bragg peaks. These systems play a vital role in a number of fundamental and applied problems: glass formation, jamming, rigidity, photonic and electronic band structure, localization of waves and excitations, self-organization, fluid dynamics, quantum systems, and pure mathematics. Much of what we know theoretically about disordered hyperuniform states of matter involves many-particle systems. In this paper, we derive new rigorous criteria that disordered hyperuniform two-phase heterogeneous materials must obey and explore their consequences. Two-phase heterogeneous media are ubiquitous; examples include composites and porous media, biological media, foams, polymer blends, granular media, cellular solids, and colloids. We begin by obtaining some results that apply to hyperuniform two-phase media in which one phase is a sphere packing in d-dimensional Euclidean space [Formula: see text]. Among other results, we rigorously establish the requirements for packings of spheres of different sizes to be 'multihyperuniform'. We then consider hyperuniformity for general two-phase media in [Formula: see text]. Here we apply realizability conditions for an autocovariance function and its associated spectral density of a two-phase medium, and then incorporate hyperuniformity as a constraint in order to derive new conditions. We show that some functional forms can immediately be eliminated from consideration and identify other forms that are allowable. Specific examples and counterexamples are described. Contact is made with well-known microstructural models (e.g. overlapping spheres and checkerboards) as well as irregular phase-separation and Turing-type patterns. We also ascertain a family of autocovariance functions (or spectral densities) that are realizable by disordered hyperuniform two-phase media in any space dimension, and present select explicit constructions of realizations. These studies provide insight into the nature of disordered hyperuniformity in the context of heterogeneous materials and have implications for the design of such novel amorphous materials.

Calculation of displacement correlation tensor indicating vortical cooperative motion in two-dimensional colloidal liquids

As an indicator of cooperative motion in a system of Brownian particles that models two-dimensional colloidal liquids, a displacement correlation tensor is calculated analytically and compared with numerical results. The key idea for the analytical calculation is to relate the displacement correlation tensor, which is a kind of four-point space-time correlation, to the Lagrangian two-time correlation of the deformation gradient tensor. Tensorial treatment of the statistical quantities, including the displacement correlation itself, allows capturing the vortical structure of the cooperative motion. The calculated displacement correlation also implies a negative long-time tail in the velocity autocorrelation, which is a manifestation of the cage effect. Both the longitudinal and transverse components of the displacement correlation are found to be expressible in terms of a similarity variable, suggesting that the cages are nested to form a self-similar structure in the space-time.

Hyperuniformity and its generalizations

Disordered many-particle hyperuniform systems are exotic amorphous states of matter that lie between crystal and liquid: They are like perfect crystals in the way they suppress large-scale density fluctuations and yet are like liquids or glasses in that they are statistically isotropic with no Bragg peaks. These exotic states of matter play a vital role in a number of problems across the physical, mathematical as well as biological sciences and, because they are endowed with novel physical properties, have technological importance. Given the fundamental as well as practical importance of disordered hyperuniform systems elucidated thus far, it is natural to explore the generalizations of the hyperuniformity notion and its consequences. In this paper, we substantially broaden the hyperuniformity concept along four different directions. This includes generalizations to treat fluctuations in the interfacial area (one of the Minkowski functionals) in heterogeneous media and surface-area driven evolving microstructures, random scalar fields, divergence-free random vector fields, and statistically anisotropic many-particle systems and two-phase media. In all cases, the relevant mathematical underpinnings are formulated and illustrative calculations are provided. Interfacial-area fluctuations play a major role in characterizing the microstructure of two-phase systems (e.g., fluid-saturated porous media), physical properties that intimately depend on the geometry of the interface, and evolving two-phase microstructures that depend on interfacial energies (e.g., spinodal decomposition). In the instances of random vector fields and statistically anisotropic structures, we show that the standard definition of hyperuniformity must be generalized such that it accounts for the dependence of the relevant spectral functions on the direction in which the origin in Fourier space is approached (nonanalyticities at the origin). Using this analysis, we place some well-known energy spectra from the theory of isotropic turbulence in the context of this generalization of hyperuniformity. Among other results, we show that there exist many-particle ground-state configurations in which directional hyperuniformity imparts exotic anisotropic physical properties (e.g., elastic, optical, and acoustic characteristics) to these states of matter. Such tunability could have technological relevance for manipulating light and sound waves in ways heretofore not thought possible. We show that disordered many-particle systems that respond to external fields (e.g., magnetic and electric fields) are a natural class of materials to look for directional hyperuniformity. The generalizations of hyperuniformity introduced here provide theoreticians and experimentalists new avenues to understand a very broad range of phenomena across a variety of fields through the hyperuniformity "lens."

Glass susceptibility: Growth kinetics and saturation under shear

We study the growth kinetics of glassy correlations in a structural glass by monitoring the evolution, within mode-coupling theory, of a suitably defined three-point function χ_{C}(t,t_{w}) with time t and waiting time t_{w}. From the complete wave-vector-dependent equations of motion for domain growth, we pass to a schematic limit to obtain a numerically tractable form. We find that the peak value χ_{C}^{P} of χ_{C}(t,t_{w}), which can be viewed as a correlation volume, grows as t_{w}^{0.5}, and the relaxation time as t_{w}^{0.8}, following a quench to a point deep in the glassy state. These results constitute a theoretical explanation of the simulation findings of Parisi [J. Phys. Chem. B 103, 4128 (1999)JPCBFK1520-610610.1021/jp983967m] and Kob and Barrat [Phys. Rev. Lett. 78, 4581 (1997)PRLTAO0031-900710.1103/PhysRevLett.78.4581], and they are also in qualitative agreement with Parsaeian and Castillo [Phys. Rev. E 78, 060105(R) (2008)PLEEE81539-375510.1103/PhysRevE.78.060105]. On the other hand, if the quench is to a point on the liquid side, the correlation volume grows to saturation. We present a similar calculation for the growth kinetics in a p-spin spin glass mean-field model where we find a slower growth, χ_{C}^{P}∼t_{w}^{0.13}. Further, we show that a shear rate γ[over ̇] cuts off the growth of glassy correlations when t_{w}∼1/γ[over ̇] for quench in the glassy regime and t_{w}=min(t_{r},1/γ[over ̇]) in the liquid, where t_{r} is the relaxation time of the unsheared liquid. The relaxation time of the steady-state fluid in this case is ∝γ[over ̇]^{-0.8}.

Universal scaling in the aging of the strong glass former SiO2

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

Relaxation and physical aging in network glasses: a review

Recent progress in the description of glassy relaxation and aging are reviewed for the wide class of network-forming materials such as GeO2, Ge x Se1-x , silicates (SiO2-Na2O) or borates (B2O3-Li2O), all of which have an important usefulness in domestic, geological or optoelectronic applications. A brief introduction of the glass transition phenomenology is given, together with the salient features that are revealed both from theory and experiments. Standard experimental methods used for the characterization of the slowing down of the dynamics are reviewed. We then discuss the important role played by aspects of network topology and rigidity for the understanding of the relaxation of the glass transition, while also permitting analytical predictions of glass properties from simple and insightful models based on the network structure. We also emphasize the great utility of computer simulations which probe the dynamics at the molecular level, and permit the calculation of various structure-related functions in connection with glassy relaxation and the physics of aging which reveal the non-equilibrium nature of glasses. We discuss the notion of spatial variations of structure which leads to the concept of 'dynamic heterogeneities', and recent results in relation to this important topic for network glasses are also reviewed.

Glassy dynamics of model colloidal polymers: The effect of "monomer" size

In recent years, attempts have been made to assemble colloidal particles into chains, which are termed "colloidal polymers." An apparent difference between molecular and colloidal polymers is the "monomer" size. Here, we propose a model to represent the variation from molecular polymer to colloidal polymer and study the quantitative differences in their glassy dynamics. For chains, two incompatible local length scales, i.e., monomer size and bond length, are manifested in the radial distribution function and intramolecular correlation function. The mean square displacement of monomers exhibits Rouse-like sub-diffusion at intermediate time/length scale and the corresponding exponent depends on the volume fraction and the monomer size. We find that the threshold volume fraction at which the caging regime emerges can be used as a rescaling unit so that the data of localization length versus volume fraction for different monomer sizes can gather close to an exponential curve. The increase of monomer size effectively increases the hardness of monomers and thus makes the colloidal polymers vitrify at lower volume fraction. Static and dynamic equivalences between colloidal polymers of different monomer sizes have been discussed. In the case of having the same peak time of the non-Gaussian parameter, the motion of monomers of larger size is much less non-Gaussian. The mode-coupling critical exponents for colloidal polymers are in agreement with that of flexible bead-spring chains.

Scaling and universality in glass transition

Kinetic facilitated models and the Mode Coupling Theory (MCT) model B are within those systems known to exhibit a discontinuous dynamical transition with a two step relaxation. We consider a general scaling approach, within mean field theory, for such systems by considering the behavior of the density correlator 〈q(t)〉 and the dynamical susceptibility 〈q(2)(t)〉 - 〈q(t)〉(2). Focusing on the Fredrickson and Andersen (FA) facilitated spin model on the Bethe lattice, we extend a cluster approach that was previously developed for continuous glass transitions by Arenzon et al. (Phys. Rev. E 90, 020301(R) (2014)) to describe the decay to the plateau, and consider a damage spreading mechanism to describe the departure from the plateau. We predict scaling laws, which relate dynamical exponents to the static exponents of mean field bootstrap percolation. The dynamical behavior and the scaling laws for both density correlator and dynamical susceptibility coincide with those predicted by MCT. These results explain the origin of scaling laws and the universal behavior associated with the glass transition in mean field, which is characterized by the divergence of the static length of the bootstrap percolation model with an upper critical dimension dc = 8.

Mean-field dynamic criticality and geometric transition in the Gaussian core model

We use molecular dynamics simulations to investigate dynamic heterogeneities and the potential energy landscape of the Gaussian core model (GCM). Despite the nearly Gaussian statistics of particles' displacements, the GCM exhibits giant dynamic heterogeneities close to the dynamic transition temperature. The divergence of the four-point susceptibility is quantitatively well described by the inhomogeneous version of the mode-coupling theory. Furthermore, the potential energy landscape of the GCM is characterized by large energy barriers, as expected from the lack of activated, hopping dynamics, and display features compatible with a geometric transition. These observations demonstrate that all major features of mean-field dynamic criticality can be observed in a physically sound, three-dimensional model.

Short-Time Beta Relaxation in Glass-Forming Liquids Is Cooperative in Nature

Temporal relaxation of density fluctuations in supercooled liquids near the glass transition occurs in multiple steps. Using molecular dynamics simulations for three model glass-forming liquids, we show that the short-time β relaxation is cooperative in nature. Using finite-size scaling analysis, we extract a growing length scale associated with beta relaxation from the observed dependence of the beta relaxation time on the system size. We find, in qualitative agreement with the prediction of the inhomogeneous mode coupling theory, that the temperature dependence of this length scale is the same as that of the length scale that describes the spatial heterogeneity of local dynamics in the long-time α-relaxation regime.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.

Slow dynamics of a tagged particle in a supercooled liquid

The ergodicity-nonergodicity (ENE) transition of the self-consistent mode-coupling theory (MCT) is marked by the point at which the time correlation of collective density fluctuations is not zero in the long-time limit. The nonergodic state, reaching beyond the ENE transition of simple MCT, is characterized by a finite shear modulus. The MCT, formulated in the current set of papers, predicts that the single-particle density correlation, unlike the collective density correlation, decays to zero at long times on either side of the ENE transition. The self-diffusion coefficient remains finite. This differs from the existing MCT results in which both collective and single-particle correlations are simultaniously frozen at the ENE transition. We discuss in this paper mechanisms by which a sharp fall in self-diffusion coefficient may occur within the present model. This overdamping or the so-called adiabatic approximation for the supercooled state does not maintain microscopic momentum conservation. Within this approximation, the self-diffusion constant approaches zero at the ENE transition point. This approximate result, which is similar to the prediction of the existing MCT models, further illustrates the process of cage formation with increase of density. At a qualitative level, our analysis shows that the self-diffusion process depends on the structure as well as short-time transport properties of the supercooled liquid. We solve the integral equations for the nonergodicity parameters to analyze the full implications of the adiabatic approximation.

Field induced anisotropic cooperativity in a magnetic colloidal glass

The translational dynamics of a repulsive colloidal glass-former is probed by time-resolved X-ray Photon Correlation Spectroscopy. In this dense dispersion of charge-stabilized and magnetic nanoparticles, the interaction potential can be tuned, from quasi-isotropic to anisotropic by applying an external magnetic field. This powerful control parameter finely tunes the anisotropy of the intricate energy landscape in the colloidal glass-former, which is seen here as a new tunable model-system to probe the dynamical heterogeneities at the approach of the glass transition. Both structural and dynamical anisotropies are reported on interparticle lengthscales associated with highly anisotropic cooperativity, almost two orders of magnitude larger in the field direction than in the perpendicular direction and in zero field.

How Local and Average Particle Diffusivities of Inhomogeneous Fluids Depend on Microscopic Dynamics

Computer simulations and a stochastic Fokker-Planck equation based approach are used to compare the single-particle diffusion coefficients of equilibrium hard-sphere fluids exhibiting identical inhomogeneous static structure and governed by either Brownian (i.e., overdamped Langevin) or Newtonian microscopic dynamics. The physics of inhomogeneity is explored via the imposition of one-dimensional sinusoidal density profiles of different wavelengths and amplitudes. When imposed density variations are small in magnitude for distances on the scale of a particle diameter, bulk-like average correlations between local structure and mobility are observed. In contrast, when density variations are significant on that length scale, qualitatively different structure-mobility correlations emerge that are sensitive to the governing microscopic dynamics. Correspondingly, a previously proposed scaling between long-time diffusivities for bulk isotropic fluids of particles exhibiting Brownian versus Newtonian dynamics [Pond et al. Soft Matter 2011, 7, 9859-9862] cannot be generalized to describe the position-dependent behaviors of strongly inhomogeneous fluids. While average diffusivities in the inhomogeneous and homogeneous directions are coupled, their qualitative dependencies on inhomogeneity wavelength are sensitive to the details of the microscopic dynamics. Nonetheless, average diffusivities of the inhomogeneous fluids can be approximately predicted for either type of dynamics based on knowledge of bulk isotropic fluid behavior and how inhomogeneity modifies the distribution of available volume. Analogous predictions for average diffusivities of experimental, inhomogeneous colloidal dispersions (based on known bulk behavior) suggest that they will exhibit qualitatively different trends than those predicted by models governed by overdamped Langevin dynamics that do not account for hydrodynamic interactions.

Evidence for a Disordered Critical Point in a Glass-Forming Liquid

Using computer simulations of an atomistic glass-forming liquid, we investigate the fluctuations of the overlap between a fluid configuration and a quenched reference system. We find that large fluctuations of the overlap develop as temperature decreases, consistent with the existence of the random critical point that is predicted by effective field theories. We discuss the scaling of fluctuations near the presumed critical point, comparing the observed behavior with that of the random-field Ising model. We argue that this critical point directly reveals the existence of an interfacial tension between amorphous metastable states, a quantity relevant both for equilibrium relaxation and for nonequilibrium melting of stable glass configurations.