Dynamical Heterogeneity in Glass-Forming Liquids
We review the phenomena of dynamical heterogeneity in glass-forming systems and its description within replica and mean-field theories of the glass transition.
Discover
Research tools
Network
Opportunities
Account
Source author record
Kunimasa Miyazaki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We review the phenomena of dynamical heterogeneity in glass-forming systems and its description within replica and mean-field theories of the glass transition.
Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalised 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 new framework allows for a novel and fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition.
We extend the replica liquid theory in order to describe the multiple glass transitions of binary mixtures with large size disparities, by taking into account the two-step replica symmetry breaking (2RSB). We determine the glass phase diagram of the mixture of large and small particles in the large-dimension limit where the mean-field theory becomes exact. When the size ratio of particles is beyond a critical value, the theory predicts three distinct glass phases; (i) the 1RSB double glass where both components vitrify simultaneously, (ii) the 1RSB single glass where only large particles are frozen while small particles remain mobile, and (iii) a new glass phase called the 2RSB double glass where both components vitrify simultaneously but with an energy landscape topography distinct from the 1RSB double glass.
We perform molecular dynamics simulations to investigate the effect of a glass preparation on its yielding transition under oscillatory shear. We use swap Monte Carlo to investigate a broad range of glass stabilities from poorly annealed to highly stable systems. We observe a qualitative change in the nature of yielding, which evolves from ductile to brittle as glass stability increases. Our results disentangle the relative role of mechanical and thermal annealing on the mechanical properties of amorphous solids, which is relevant for various experimental situations from the rheology of soft materials to fatigue failure in metallic glasses.
We investigate criticality near the jamming transition in both quiescent systems and those under shear by considering the effect of mechanical training on the jamming transition and nonlinear rheology. We simulate frictionless soft particles undergoing athermal quasi-static shear using initial configurations trained with athermal quasi-static cyclic volume deformations. The jamming transition density of the initial configuration $φ_{\rm J0}$ is systematically altered by tuning the ``depth'' of mechanical training. We exert a steady shear on these configurations and observe either shear jamming (gain of stiffness due to shear) or shear melting (loss of stiffness due to shear), depending on the depth of training and proximity to the jamming transition density. We also observe that the characteristic strains, at which shear jamming or melting occur, diverge at a unique density $φ_{\rm JS}$. This is due to the shift of the jamming transition density from $φ_{\rm J0}$ to $φ_{\rm JS}$ under shear, associated with loss of memory of the initial configuration. Finally, we thoroughly investigate nonlinear rheology near the jamming transition density, and contrary to previous works, we find a nonlinear ``softening'' takes place below as well as above the jamming transition density.
We study supercooled dynamics in quantum hard-sphere liquid using quantum mode-coupling formulation. In the moderate quantum regime, classical cage effects lead to slower dynamics compared to strongly quantum regime, where tunneling overcomes classical caging, leading to faster relaxation. As a result, the glass transition critical density can become significantly higher than for the classical liquids. Perturbative approach is used to solve time dependent quantum mode-coupling equations to study in detail the dynamics of the supercooled liquid in moderate quantum regime. Similar to the classical case, relaxation time shows power-law increase with increasing density in the supercooled regime. However, the power-law exponent is found to be dependent on the quantumness; it increases linearly as the quantumness is increased in the moderate quantum regime.
Using molecular dynamics simulation, we investigate the slow dynamics of a supercooled binary mixture of soft particles interacting with a generalized Hertzian potential. At low density, it displays typical slow dynamics near its glass transition temperature. At higher densities, particles bond together, forming clusters, and the clusters undergo the glass transition. The number of particles in a cluster increases one by one as the density increases. We demonstrate that there exist the multiple cluster-glass phases characterized by a different number of particles per cluster, each of which is separated by distinct minima. Surprisingly, a so-called higher order singularity of the mode-coupling theory signaled by a logarithmic relaxation is observed in the vicinity of the boundaries between monomer and cluster glass phases. The system also exhibits rich and anomalous dynamics in the cluster glass phases, such as the decoupling of the self- and collective dynamics.
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 a geometric transition and large energy barriers, as expected from the lack of activated, hopping dynamics. These observations demonstrate that all major features of mean-field dynamic criticality can be observed in a physically sound, three-dimensional model.
We perform molecular dynamics simulations for a silica glass former model proposed by Coslovich and Pastore (CP) over a wide range of densities. The density variation can be mapped onto the change of the potential depth between Si and O interactions of the CP model. By reducing the potential depth (or increasing the density), the anisotropic tetrahedral network structure observed in the original CP model transforms into the isotropic structure with the purely repulsive soft-sphere potential. Correspondingly, the temperature dependence of the relaxation time exhibits the crossover from the Arrhenius to the super-Arrhenius behavior. Being able to control the fragility over a wide range by tuning the potential of a single model system helps us to bridge the gap between the network and isotropic glass formers and to obtain the insight into the underlying mechanism of the fragility. We study the relationship between the fragility and dynamical properties such as the magnitude of the Stokes-Einstein violation and the stretch exponent in the density correlation function. We also demonstrate that the peak of the specific heat systematically shifts as the density increases, hinting that the fragility is correlated with the hidden thermodynamic anomalies of the system.
We study the effects of random pinning on the Fredrickson-Andersen model on the Bethe lattice. We find that the nonergodic transition temperature rises as the fraction of the pinned spins increases and the transition line terminates at a critical point. The freezing behavior of the spins is analogous to that of a randomly pinned p-spin mean-field spin glass model which has been recently reported. The diverging behavior of correlation lengths in the vicinity of the terminal critical point is found to be identical to the prediction of the inhomogeneous mode-coupling theory at the A3 singularity point for the glass transition.
Reply to the preceding letter by Chakrabarty et al. [Proc. Nat. Acad. Sci. USA, 112, no. 35, E4818-E4820 (2015)].
We analyse, using Inhomogenous Mode-Coupling Theory, the critical scaling behaviour of the dynamical susceptibility at a distance epsilon from continuous second-order glass transitions. We find that the dynamical correlation length xi behaves generically as epsilon^{-1/3} and that the upper critical dimension is equal to six. More surprisingly, we find activated dynamic scaling, where xi grows with time as [ln(t)]^2 exactly at criticality. All these results suggest a deep analogy between the glassy behaviour of attractive colloids or randomly pinned supercooled liquids and that of the Random Field Ising Model.
The Langevin equation with multiplicative noise and state-dependent transport coefficient has to be always complemented with the proper interpretation rule of the noise, such as the Ito and Stratonovich conventions. Although the mathematical relationship between the different rules and how to translate from one rule to another are well-established, it still remains controversial what is a more {\it physically} natural rule. In this communication, we derive the overdamped Langevin equation with multiplicative noise for Brownian particles, by systematically eliminating the fast degrees of freedom of the underdamped Langevin equation. The Langevin equations obtained here vary depending on the choice of the noise conventions but they are different representations for an identical phenomenon. The results apply to multi-variable, nonequilibrium, non-stationary systems, and other general settings.
In this work we numerically investigate a new method for the characterization of growing length scales associated with spatially heterogeneous dynamics of glass-forming liquids. This approach, motivated by the formulation of the inhomogeneous mode-coupling theory (IMCT) [Biroli G. et al. Phys. Rev. Lett. 2006 97, 195701], utilizes inhomogeneous molecular dynamics simulations in which the system is perturbed by a spatially modulated external potential. We show that the response of the two-point correlation function to the external field allows one to probe dynamic correlations. We examine the critical properties shown by this function, in particular the associated dynamic correlation length, that is found to be comparable to the one extracted from standardly-employed four-point correlation functions. Our numerical results are in qualitative agreement with IMCT predictions but suggest that one has to take into account fluctuations not included in this mean-field approach in order to reach quantitative agreement. Advantages of our approach over the more conventional one based on four-point correlation functions are discussed.
Recent studies show that volume fractions $\phiJ$ at the jamming transition of frictionless hard spheres and discs are not uniquely determined but exist over a continuous range. Motivated by this observation, we numerically investigate dependence of $\phiJ$ on the initial configurations of the parent fluids equilibrated at a fraction $\phiini$, before compressing to generate a jammed packing. We find that $\phiJ$ remains constant when $\phiini$ is small but sharply increases when $\phiini$ exceeds the dynamic transition point which the mode-coupling theory predicts. We carefully analyze configurational properties of both jammed packings and parent fluids and find that, while all jammed packings remain isostatic, the increase of $\phiJ$ is accompanied with subtle but distinct changes of (i) local orders, (ii) a static length scale, and (iii) an exponent of the finite size scaling. These results quantitatively support the scenario of the random first order transition theoryof the glass transition.
We numerically study dynamical properties of the one-component Gaussian Core Model in the supercooled states. We find that nucleation is suppressed as density increases. Concomitantly the system exhibits glassy slow dynamics characterized by the two-step and stretched exponential relaxation of the density correlation as well as drastic increase of the relaxation time. It is found that violation of the Stokes-Einstein relation is weaker and the non-Gaussian parameter is smaller than typical model glass formers, implying weaker dynamic heterogeneities. Besides, agreement of simulation data with the prediction of mode-coupling theory is exceptionally good, indicating that the nature of slow dynamics of this ultra-soft particle fluid is mean-field-like. This fact may be understood as the consequences of multiple overlaps of the constituent particles at high densities.
Reply to the preceding comment by Rolf Schilling and Bernhard Schmid [arXiv:1101.5577, Phys. Rev. Lett. 106, 049601 (2011)].
We numerically study crystal nucleation and glassy slow dynamics of the one-component Gaussian core model (GCM) at high densities. The nucleation rate at a fixed supersaturation is found to decrease as the density increases. At very high densities, the nucleation is not observed at all in the time window accessed by long molecular dynamics (MD) simulation. Concomitantly, the system exhibits typical slow dynamics of the supercooled fluids near the glass transition point. We compare the simulation results of the supercooled GCM with the predictions of mode-coupling theory (MCT) and find that the agreement between them is better than any other model glassformers studied numerically in the past. Furthermore, we find that a violation of the Stokes-Einstein relation is weaker and the non-Gaussian parameter is smaller than canonical glassformers. Analysis of the probability distribution of the particle displacement clearly reveals that the hopping effect is strongly suppressed in the high density GCM. We conclude from these observations that the GCM is more amenable to the mean-field picture of the glass transition than other models. This is attributed to the long-ranged nature of the interaction potential of the GCM in the high density regime. Finally, the intermediate scattering function at small wavevectors is found to decay much faster than its self part, indicating that dynamics of the large-scale density fluctuations decouples with the shorter-ranged caging motion.
Using molecular dynamics simulations, we study the slow dynamics of supercooled liquids confined in a random matrix of immobile obstacles. We study the dynamical crossover from glass-like to Lorentz-gas-like behavior in terms of the density correlation function, the mean square displacement, the nonlinear dynamic susceptibility, the non-Gaussian parameter, and the fragility. Cooperative and spatially heterogeneous dynamics are suppressed as the obstacle density increases, which lead to the more Arrhenius-like behavior in the temperature dependence of the relaxation time. Our findings are qualitatively consistent with the results of recent experimental and numerical studies for various classes of spatially heterogeneous systems. We also investigate the dependence of the dynamics of mobile particles on the protocol to generate the random matrix. A reentrant transition from the arrested phase to the liquid phase as the mobile particle density {\it increases} is observed for a class of protocols. This reentrance is explained in terms of the distribution of the volume of the voids that are available to the mobile particles.
A comprehensive microscopic dynamical theory is presented for the description of quantum fluids as they transform into glasses. The theory is based on a quantum extension of mode-coupling theory. Novel effects are predicted, such as reentrant behavior of dynamical relaxation times. These predictions are supported by path integral ring polymer molecular dynamics simulations. The simulations provide detailed insight into the factors that govern slow dynamics in glassy quantum fluids. Connection to other recent work on both quantum glasses as well as quantum optimization problems is presented.
We numerically study thermodynamic and structural properties of the one-component Gaussian core model (GCM) at very high densities. The solid-fluid phase boundary is carefully determined. We find that the density dependence of both the freezing and melting temperatures obey the asymptotic relation, $\log T_f$, $\log T_m \propto -ρ^{2/3}$, where $ρ$ is the number density, which is consistent with Stillinger's conjecture. Thermodynamic quantities such as the energy and pressure and the structural functions such as the static structure factor are also investigated in the fluid phase for a wide range of temperature above the phase boundary. We compare the numerical results with the prediction of the liquid theory with the random phase approximation (RPA). At high temperatures, the results are in almost perfect agreement with RPA for a wide range of density, as it has been already shown in the previous studies. In the low temperature regime close to the phase boundary line, although RPA fails to describe the structure factors and the radial distribution functions at the length scales of the interparticle distance, it successfully predicts their behaviors at shorter length scales. RPA also predicts thermodynamic quantities such as the energy, pressure, and the temperature at which the thermal expansion coefficient becomes negative, almost perfectly. Striking ability of RPA to predict thermodynamic quantities even at high densities and low temperatures is understood in terms of the decoupling of the length scales which dictate thermodynamic quantities from the interparticle distance which dominates the peak structures of the static structure factor due to the softness of the Gaussian core potential.
Mode-coupling theory (MCT) is conjectured to be a mean-field description of dynamics of the structural glass transition and the replica theory to be its thermodynamic counterpart. However, the relationship between the two theories remains controversial and quantitative comparison is lacking. In this Letter, we investigate MCT for monatomic hard sphere fluids at arbitrary dimensions above three and compare the results with replica theory. We find grave discrepancies between the predictions of two theories. While MCT describes the nonergodic parameter quantitatively better than the replica theory in three dimension, it predicts a completely different dimension dependence of the dynamical transition point. We find it to be due to the pathological behavior of the nonergodic parameters derived from MCT, which exhibit negative tails in real space at high dimensions.
The very nature of glass is somewhat mysterious: while relaxation times in glasses are of sufficient magnitude that large-scale motion on the atomic level is essentially as slow as it is in the crystalline state, the structure of glass appears barely different than that of the liquid that produced it. Quantum mechanical systems ranging from electron liquids to superfluid helium appear to form glasses, but as yet no unifying framework exists connecting classical and quantum regimes of vitrification. Here we develop new insights from theory and simulation into the quantum glass transition that surprisingly reveal distinct regions where quantum fluctuations can either promote or inhibit glass formation.