Strong dynamical heterogeneity and universal scaling in driven granular fluids

Method to probe the pronounced growth of correlation lengths in active glass-forming liquids using an elongated probe

The growth of correlation lengths in equilibrium glass-forming liquids near the glass transition is considered a critical finding in the quest to understand the physics of glass formation. These understandings helped us understand various dynamical phenomena observed in supercooled liquids. It is known that at least two different length scales exist; one is of thermodynamic origin, while the other is dynamical in nature. Recent observations of glassy dynamics in biological and synthetic systems where the external or internal driving source controls the dynamics, apart from the usual thermal noise, lead to the emergence of the field of active glassy matter. A question of whether the physics of glass formation in these active systems is also accompanied by growing dynamic and static lengths is indeed timely. In this article, we probe the growth of dynamic and static lengths in a model active glass system using rod-like elongated probe particles, an experimentally viable method. We show that the dynamic and static lengths in these nonequilibrium systems grow much more rapidly than their passive counterparts. We then offer an understanding of the violation of the Stokes-Einstein relation and Stokes-Einstein-Debye relation using these lengths via a scaling theory.

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.

Homogeneous Crystallization in Cyclically Sheared Frictionless Grains

Many experiments over the past half century have shown that, for a range of protocols, granular materials compact under pressure and repeated small disturbances. A recent experiment on cyclically sheared spherical grains showed significant compaction via homogeneous crystallization (Rietz et al., 2018). Here we present numerical simulations of frictionless, purely repulsive spheres undergoing cyclic simple shear via Newtonian dynamics with linear viscous drag at fixed vertical load. We show that for sufficiently small strain amplitudes, cyclic shear gives rise to homogeneous crystallization at a volume fraction ϕ=0.646±0.001. This result indicates that neither friction nor gravity is essential for homogeneous crystallization in driven granular media. Understanding how crystal formation is initiated within a homogeneous disordered state gives key insights into the old open problem of glass formation in fluids.

Brownian motion of ellipsoidal particles on a granular magnetic bath

We study the Brownian motion of ellipsoidal particles lying on an agitated granular bath composed of magnetic particles. We quantify the mobility of different floating ellipsoidal particles using the mean-square displacement and the mean-square angular displacement, and relate the diffusion coefficients to the bath particle motion. In terms of the particle major radius R, we find the translational diffusion coefficient scales roughly as 1/R^{2} and the rotational diffusion coefficient scales as roughly 1/R^{4}; this is consistent with the assumption that diffusion arises from random kicks of the bath particles underneath the floating particle. By varying the magnetic forcing, the bath particles' diffusivity changes by a factor of ten; over this range, the translational and rotational diffusion of the floating particles change by a factor of 50. However, the ratio of the two diffusion constants for the floating particles is forcing-independent. Unusual aspects of the floating particle motion include non-Gaussian statistics for their displacements.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Persistent Anti-Correlations in Brownian Dynamics Simulations of Dense Colloidal Suspensions Revealed by Noise Suppression

Transport properties of a hard-sphere colloidal fluid are investigated by Brownian dynamics simulations. We implement a novel algorithm for the time-dependent velocity-autocorrelation function (VACF) essentially eliminating the noise of the bare random motion. The measured VACF reveals persistent anti-correlations manifested by a negative algebraic power-law tail t^{-5/2} at all densities. At small packing fractions the simulations fully agree with the analytic low-density prediction, yet the amplitude of the tail becomes dramatically suppressed as the packing fraction is increased. The mode-coupling theory of the glass transition provides a qualitative explanation for the strong variation in terms of the static compressibility as well as the slowing down of the structural relaxation.

Non-Gaussianity of the van Hove function and dynamic-heterogeneity length scale

Non-Gaussian nature of the probability distribution of particles' displacements in the supercooled temperature regime in glass-forming liquids are believed to be one of the major hallmarks of glass transition. It has already been established that this probability distribution, which is also known as the van Hove function, shows universal exponential tail. The origin of such an exponential tail in the distribution function is attributed to the hopping motion of particles observed in the supercooled regime. The non-Gaussian nature can also be explained if one assumes that the system has heterogeneous dynamics in space and time. Thus exponential tail is the manifestation of dynamic heterogeneity. In this work we directly show that non-Gaussianity of the distribution of particles' displacements occur over the dynamic heterogeneity length scale and the dynamical behavior course grained over this length scale becomes homogeneous. We study the non-Gaussianity of the van Hove function by systematically coarse graining at different length scales and extract the length scale of dynamic heterogeneity at which the shape of the van Hove function crosses over from non-Gaussian to Gaussian. The obtained dynamic heterogeneity scale is found to be in very good agreement with the scale obtained from other conventional methods.

Spatiotemporal intermittency and localized dynamic fluctuations upon approaching the glass transition

We introduce a robust approach for characterizing spatially and temporally heterogeneous behavior within a system based on the evolution of dynamic fluctuations averaged over different space lengths and timescales. We apply it to investigate the dynamics in two canonical systems as the glass transition is approached: simulated Lennard-Jones liquids and experimental dense colloidal suspensions. In both cases the onset of glassiness is marked by spatially localized dynamic fluctuations originating in regions of correlated mobile particles. By removing the trivial system size dependence we show that the spatial heterogeneity of the dynamics extends to large length scales containing tens to hundreds of particles, corresponding to the timescale of maximally non-Gaussian dynamics.

Block Analysis for the Calculation of Dynamic and Static Length Scales in Glass-Forming Liquids

We present block analysis, an efficient method of performing finite-size scaling for obtaining the length scale of dynamic heterogeneity and the point-to-set length scale for generic glass-forming liquids. This method involves considering blocks of varying sizes embedded in a system of a fixed (large) size. The length scale associated with dynamic heterogeneity is obtained from a finite-size scaling analysis of the dependence of the four-point dynamic susceptibility on the block size. The block size dependence of the variance of the α relaxation time yields the static point-to-set length scale. The values of the obtained length scales agree quantitatively with those obtained from other conventional methods. This method provides an efficient experimental tool for studying the growth of length scales in systems such as colloidal glasses for which performing finite-size scaling by carrying out experiments for varying system sizes may not be feasible.

Slow and long-ranged dynamical heterogeneities in dissipative fluids

A two-dimensional bidisperse granular fluid is shown to exhibit pronounced long-ranged dynamical heterogeneities as dynamical arrest is approached. Here we focus on the most direct approach to study these heterogeneities: we identify clusters of slow particles and determine their size, Nc, and their radius of gyration, RG. We show that , providing direct evidence that the most immobile particles arrange in fractal objects with a fractal dimension, df, that is observed to increase with packing fraction ϕ. The cluster size distribution obeys scaling, approaching an algebraic decay in the limit of structural arrest, i.e., ϕ→ϕc. Alternatively, dynamical heterogeneities are analyzed via the four-point structure factor S4(q,t) and the dynamical susceptibility χ4(t). S4(q,t) is shown to obey scaling in the full range of packing fractions, 0.6 ≤ϕ≤ 0.805, and to become increasingly long-ranged as ϕ→ϕc. Finite size scaling of χ4(t) provides a consistency check for the previously analyzed divergences of χ4(t) ∝ (ϕ-ϕc)(-γχ) and the correlation length ξ∝ (ϕ-ϕc)(-γξ). We check the robustness of our results with respect to our definition of mobility. The divergences and the scaling for ϕ→ϕc suggest a non-equilibrium glass transition which seems qualitatively independent of the coefficient of restitution.

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.

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.

Fundamental differences between glassy dynamics in two and three dimensions

The two-dimensional freezing transition is very different from its three-dimensional counterpart. In contrast, the glass transition is usually assumed to have similar characteristics in two and three dimensions. Using computer simulations, here we show that glassy dynamics in supercooled two- and three-dimensional fluids are fundamentally different. Specifically, transient localization of particles on approaching the glass transition is absent in two dimensions, whereas it is very pronounced in three dimensions. Moreover, the temperature dependence of the relaxation time of orientational correlations is decoupled from that of the translational relaxation time in two dimensions but not in three dimensions. Last, the relationships between the characteristic size of dynamically heterogeneous regions and the relaxation time are very different in two and three dimensions. These results strongly suggest that the glass transition in two dimensions is different than in three dimensions.

The existing glass transition theories show little dimensional dependence. Here Flenner et al. disprove this general consensus by finding that, for instance, the bond-orientational relaxation is decoupled from the translational relaxation in two dimensions, but not in three dimensions.