Critical test of the mode-coupling theory of the glass transition

Length-scale dependence of Stokes-Einstein breakdown in active glass-forming liquids

The Stokes-Einstein (SE) relation, which relates diffusion constants with the viscosity of a liquid at high temperatures in equilibrium, is violated in the supercooled temperature regime. Whether this relation is obeyed in nonequilibrium active liquids is a question of significant current interest to the statistical physics community trying to develop the theoretical framework of nonequilibrium statistical mechanics. Via extensive computer simulations of model active glass-forming liquids in three dimensions, we show that SE is obeyed at a high temperature similar to the equilibrium behavior, and it gets violated in the supercooled temperature regimes. The degree of violation increases systematically with the increasing activity which quantifies the amount the system is driven out of equilibrium. First passage-time (FPT) distributions helped us to gain insights into this enhanced breakdown from the increased short-time peak, depicting hoppers. Subsequently, we study the wave vector dependence of SE relation and show that it gets restored at a wave vector that decreases with increasing activity, and the crossover wave vector is found to be proportional to the inverse of the dynamical heterogeneity length scale in the system. Our work showed how SE violation in active supercooled liquids could be rationalized using the growth of dynamic length scale, which is found to grow enormously with increasing activity in these systems.

Microscopic theory of the elastic shear modulus and length-scale-dependent dynamic re-entrancy phenomena in very dense sticky particle fluids

We apply the hybrid projectionless dynamic theory (hybrid PDT) formulation of the elastically collective nonlinear Langevin equation (ECNLE) activated dynamics approach to study dense fluids of sticky spheres interacting with short range attractions. Of special interest is the problem of non-monotonic evolution with short range attraction strength of the elastic modulus ("re-entrancy") at very high packing fractions far beyond the ideal mode coupling theory (MCT) nonergodicity boundary. The dynamic force constraints explicitly treat the bare attractive forces that drive transient physical bond formation, while a projection approximation is employed for the singular hard-sphere potential. The resultant interference between repulsive and attractive forces contribution to the dynamic vertex results in the prediction of localization length and elastic modulus re-entrancy, qualitatively consistent with experiments. The non-monotonic evolution of the structural (alpha) relaxation time predicted by the ECNLE theory with the hybrid PDT approach is explored in depth as a function of packing fraction, attraction strength, and attraction range. Isochronal dynamic arrest boundaries based on activated relaxation display the classic non-monotonic glass melting form. Comparisons of these results with the corresponding predictions of ideal MCT, and also the ECNLE and NLE activated theories based on projection, reveal large qualitative differences. The consequences of stochastic trajectory fluctuations on intra-cage single particle dynamics with variable strength of attractions are also studied. Large dynamical heterogeneity effects in attractive glasses are properly captured. These include a rapidly increasing amplitude of the non-Gaussian parameter with packing fraction and a non-monotonic evolution with attraction strength, in qualitative accord with recent simulations. Extension of the microscopic theoretical approach to treat double yielding in attractive glass nonlinear rheology is possible.

Importance of Many Particle Correlations to the Collective Debye-Waller Factor in a Single-Particle Activated Dynamic Theory of the Glass Transition

We theoretically study the importance of many body correlations on the collective Debye-Waller (DW) factor in the context of the Nonlinear Langevin Equation (NLE) single-particle activated dynamics theory of glass transition and its extension to include collective elasticity (ECNLE theory). This microscopic force-based approach envisions structural alpha relaxation as a coupled local-nonlocal process involving correlated local cage and longer range collective barriers. The crucial question addressed here is the importance of the deGennes narrowing contribution versus a literal Vineyard approximation for the collective DW factor that enters the construction of the dynamic free energy in NLE theory. While the Vineyard-deGennes approach-based NLE theory and its ECNLE theory extension yields predictions that agree well with experimental and simulation results, use of a literal Vineyard approximation for the collective DW factor massively overpredicts the activated relaxation time. The current study suggests many particle correlations are crucial for a reliable description of activated dynamics theory of model hard sphere fluids.

Structural Signatures of Ultrastability in a Deposited Glassformer

Glasses obtained from vapor deposition on a cold substrate have superior thermodynamic and kinetic stability with respect to ordinary glasses. Here we perform molecular dynamics simulations of vapor deposition of a model glassformer and investigate the origin of its high stability compared to that of ordinary glasses. We find that the vapor deposited glass is characterized by locally favored structures (LFSs) whose occurrence correlates with its stability, reaching a maximum at the optimal deposition temperature. The formation of LFSs is enhanced near the free surface, hence supporting the idea that the stability of vapor deposited glasses is connected to the relaxation dynamics at the surface.

A fractal structural feature related to dynamic crossover in metallic glass-forming liquids

The dynamic crossover in supercooled liquids initially predicted by model coupling theory has been widely accepted, but its underlying structural origin is still an open issue for glass-forming liquids. By molecular dynamics simulations of binary CuZr liquids, the present work verifies that high pressure could enhance this crossover, facilitating the studies on the structural features at the crossover temperature T_c. We discover that the topological connectivity of icosahedral clusters is responsible for this dynamic crossover, rather than all clusters. T_c is the temperature at which the connectivity degree between these clusters reaches a maximum and the dynamic heterogeneity begins to keep stable. Below T_c, the fractal topological structures appear in the medium-range order scale. The icosahedral clusters with a certain connectivity pattern can be regarded as a fractal structural unit. By employing the established fractal analysis method, the fractal dimension D of the icosahedral network is calculated. Our results indicate that the D value increases monotonically with increasing pressure and the fractal behavior of the icosahedral network is an inherent feature of metallic glasses. We also find similar fractal behavior in clusters with high local five-fold symmetry. Our findings shed light on the origin of a dynamic crossover in the deep supercooled region of metallic glasses and also demonstrate the important role of icosahedral clusters in uncovering the fractal behavior of metallic glass.

Many-Body Correlations Are Non-negligible in Both Fragile and Strong Glassformers

It is widely believed that the emergence of slow glassy dynamics is encoded in a material's microstructure. First-principles theory [mode-coupling theory (MCT)] is able to predict the dramatic slowdown of the dynamics from only static two-point correlations as input, yet it cannot capture all of the observed dynamical behavior. Here we go beyond two-point spatial correlation functions by extending MCT systematically to include higher-order static and dynamic correlations. We demonstrate that only adding the static triplet direct correlations already qualitatively changes the predicted glass-transition diagram of binary hard spheres and silica. Moreover, we find a nontrivial competition between static triplet correlations that work to stabilize the glass state and dynamic higher-order correlations that destabilize it for both materials. We conclude that the conventionally neglected static triplet direct correlations as well as higher-order dynamic correlations are, in fact, non-negligible in both fragile and strong glassformers.

Comparing Microscopic and Macroscopic Dynamics in a Paradigmatic Model of Glass-Forming Molecular Liquid

Glass transition is a most intriguing and long-standing open issue in the field of molecular liquids. From a macroscopic perspective, glass-forming systems display a dramatic slowing-down of the dynamics, with the inverse diffusion coefficient and the structural relaxation times increasing by orders of magnitude upon even modest supercooling. At the microscopic level, single-molecule motion becomes strongly intermittent, and can be conveniently described in terms of “cage-jump” events. In this work, we investigate a paradigmatic glass-forming liquid, the Kob–Andersen Lennard–Jones model, by means of Molecular Dynamics simulations, and compare the macroscopic and microscopic descriptions of its dynamics on approaching the glass-transition. We find that clear changes in the relations between macroscopic timescales and cage-jump quantities occur at the crossover temperature where Mode Coupling-like description starts failing. In fact, Continuous Time Random Walk and lattice model predictions based on cage-jump statistics are also violated below the crossover temperature, suggesting the onset of a qualitative change in cage-jump motion. Interestingly, we show that a fully microscopic relation linking cage-jump time- and length-scales instead holds throughout the investigated temperature range.

Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions

We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.

Generalized mode-coupling theory for mixtures of Brownian particles

Generalized mode-coupling theory (GMCT) has recently emerged as a promising first-principles theory to study the poorly understood dynamics of glass-forming materials. Formulated as a hierarchical extension of standard mode-coupling theory (MCT), it is able to systematically improve its predictions by including the exact dynamics of higher-order correlation functions into its hierarchy. However, in contrast to Newtonian dynamics, a fully generalized version of the theory based on Brownian dynamics is still lacking. To close this gap, we provide a detailed derivation of GMCT for colloidal mixtures obeying a many-body Smoluchowski equation. We demonstrate that a hierarchy of coupled equations can again be established and show that these, consistent with standard MCT, are identical to the ones obtained from Newtonian GMCT when taking the overdamped limit. Consequently, the nontrivial similarity between Brownian and Newtonian MCT is maintained for our multicomponent GMCT. As a proof of principle, we also solve the generalized mode-coupling equations for the binary Kob-Andersen Lennard-Jones mixture undergoing Brownian dynamics and confirm the improved predictive power of the theory upon using more levels of the GMCT hierarchy of equations.

Energy landscapes for a modified repulsive Weeks-Chandler-Andersen potential

The short-range nature of the repulsive Weeks-Chandler-Anderson (WCA) potential can create free particles/rattlers in a condensed system. The presence of rattlers complicates the analysis of the energy landscape due to extra zero-frequency normal modes. By employing a long-range Gaussian tail modification, we remove the rattlers without changing the structure and the dynamics of the system, and successfully describe the potential energy landscape in terms of minima and transition states. This coarse-grained description of the landscape and the dynamical properties of the modified potential exhibit characteristic signatures of glass-forming liquids. However, we show that despite having qualitatively similar behaviour, the modified WCA potential is less frustrated compared to its attractive counterpart.

Glassy dynamics of sticky hard spheres beyond the mode-coupling regime

Sticky hard spheres, i.e., hard particles decorated with a short-ranged attractive interaction potential, constitute a relatively simple model with highly non-trivial glassy dynamics. The mode-coupling theory of the glass transition (MCT) offers a qualitative account of the complex reentrant dynamics of sticky hard spheres, but the predicted glass transition point is notoriously underestimated. Here we apply an improved first-principles-based theory, referred to as generalized mode-coupling theory (GMCT), to sticky hard spheres. This theoretical framework seeks to go beyond MCT by hierarchically expanding the dynamics in higher-order density correlation functions. We predict the phase diagrams from the first few levels of the GMCT hierarchy and the dynamics-related critical exponents, all of which are much closer to the empirical observations than MCT. Notably, the prominent reentrant glassy dynamics, the glass-glass transition, and the higher-order bifurcation singularity classes (A3 and A4) of sticky hard spheres are found to be preserved within GMCT at arbitrary order. Moreover, we demonstrate that when the hierarchical order of GMCT increases, the effect of the short-ranged attractive interactions becomes more evident in the dynamics. This implies that GMCT is more sensitive to subtle microstructural differences than MCT, and that the framework provides a promising first-principles approach to systematically go beyond the MCT regime.

Multi-component generalized mode-coupling theory: predicting dynamics from structure in glassy mixtures

The emergence of glassy dynamics and the glass transition in dense disordered systems is still not fully understood theoretically. Mode-coupling theory (MCT) has shown to be effective in describing some of the non-trivial features of glass formation, but it cannot explain the full glassy phenomenology due to the strong approximations on which it is based. Generalized mode-coupling theory (GMCT) is a hierarchical extension of the theory, which is able to outclass MCT by carefully describing the dynamics of higher-order correlations in its generalized framework. Unfortunately, the theory has so far only been developed for single-component systems and as a result works poorly for highly polydisperse materials. In this paper, we solve this problem by developing GMCT for multi-component systems. We use it to predict the glassy dynamics of the binary Kob-Andersen Lennard-Jones mixture, as well as its purely repulsive Weeks-Chandler-Andersen analogue. Our results show that each additional level of the GMCT hierarchy gradually improves the predictive power of GMCT beyond its previous limit. This implies that our theory is able to harvest more information from the static correlations, thus being able to better understand the role of attraction in supercooled liquids from a first-principles perspective.

How attractive and repulsive interactions affect structure ordering and dynamics of glass-forming liquids

The theory developed in our previous papers [Phys. Rev. E 99, 030101(R) (2019)10.1103/PhysRevE.99.030101; Phys. Rev. E 103, 032611 (2021)10.1103/PhysRevE.103.032611] is applied in this paper to investigate the dependence of slowing down of dynamics of glass-forming liquids on the attractive and repulsive parts of intermolecular interactions. Through an extensive comparison of the behavior of a Lennard-Jones glass-forming liquid and that of its WCA reduction to a model with truncated pair potential without attractive tail, we demonstrate why the two systems exhibit very different dynamics despite having nearly identical pair correlation functions. In particular, we show that local structures characterized by the number of mobile and immobile particles around a central particle markedly differ in the two systems at densities and temperatures where their dynamics show large difference and nearly identical where dynamics nearly overlap. We also show how the parameter ψ(T) that measures the role of fluctuations embedded in the system on size of the cooperatively reorganizing cluster (CRC) and the crossover temperature T_{a} depend on the intermolecular interactions. These parameters stemming from the intermolecular interactions characterize the temperature and density dependence of structural relaxation time τ_{α}. The quantitative and qualitative agreements found with simulation results for the two systems suggest that our theory brings out the underlying features that determine the dynamics of glass-forming liquids.

Microscopic Theory of Softness in Supercooled Liquids

We introduce a new measure of the structure of a liquid which is the softness of the mean-field potential developed by us earlier. We find that this softness is sensitive to small changes in the structure. Then, we study its correlation with the supercooled liquid dynamics. The study involves a wide range of liquids (fragile, strong, attractive, repulsive, and active) and predicts some universal behaviors such as the softness being linearly proportional to the temperature and inversely proportional to the activation barrier of the dynamics with system dependent proportionality constants. We establish a master equation between the dynamics and the softness parameter and show that, indeed, the dynamics, when scaled by the temperature and system dependent parameters, show a data collapse when plotted against softness. The dynamics of fragile liquids show a strong softness dependence, whereas that of strong liquids show a much weaker softness dependence. We also connect the present study with the earlier studies of softness involving machine learning (ML), thus, providing a theoretical framework for understanding the ML results.

MD simulations of charged binary mixtures reveal a generic relation between high- and low-temperature behavior

Experimental studies of the glassy slowdown in molecular liquids indicate that the high-temperature activation energy E_∞ of glass-forming liquids is directly related to their glass transition temperature T_g. To further investigate such a possible relation between high- and low-temperature dynamics in glass-forming liquids, we analyze the glassy dynamics of binary mixtures using molecular dynamics simulations. We consider a binary mixture of charged Lennard-Jones particles and vary the partial charges of the particles and, thus, the high-temperature activation energy and the glass transition temperature of the system. Based on previous results, we introduce a phenomenological model describing relaxation times over the whole temperature regime from high temperatures to temperatures well inside the supercooled regime. By investigating the dynamics of both particle species on molecular and diffusive length scales along isochoric and isobaric pathways, we find a quadratic charge dependence of both E_∞ and T_g, resulting in an approximately constant ratio of both quantities independent of the underlying observable, the thermodynamic ensemble, and the particle species, and this result is robust against the actual definition of T_g. This generic relation between the activation energy and the glass transition temperature indicates that high-temperature dynamics and the glassy slowdown are related phenomena, and the knowledge of E_∞ may allow us to approximately predict T_g.

Generalized mode-coupling theory of the glass transition. II. Analytical scaling laws

Generalized mode-coupling theory (GMCT) constitutes a systematically correctable, first-principles theory to study the dynamics of supercooled liquids and the glass transition. It is a hierarchical framework that, through the incorporation of increasingly many particle density correlations, can remedy some of the inherent limitations of the ideal mode-coupling theory (MCT). However, despite MCT's limitations, the ideal theory also enjoys several remarkable successes, notably including the analytical scaling laws for the α- and β-relaxation dynamics. Here, we mathematically derive similar scaling laws for arbitrary-order multi-point density correlation functions obtained from GMCT under arbitrary mean-field closure levels. More specifically, we analytically derive the asymptotic and preasymptotic solutions for the long-time limits of multi-point density correlators, the critical dynamics with two power-law decays, the factorization scaling laws in the β-relaxation regime, and the time-density superposition principle in the α-relaxation regime. The two characteristic power-law-divergent relaxation times for the two-step decay and the non-trivial relation between their exponents are also obtained. The validity ranges of the leading-order scaling laws are also provided by considering the leading preasymptotic corrections. Furthermore, we test these solutions for the Percus-Yevick hard-sphere system. We demonstrate that GMCT preserves all the celebrated scaling laws of MCT while quantitatively improving the exponents, rendering the theory a promising candidate for an ultimately quantitative first-principles theory of glassy dynamics.

Role of Attractive Interactions in Structure Ordering and Dynamics of Glass-Forming Liquids

A key question in glass physics is what the origin of slow glassy dynamics is. The liquid structure is a natural candidate; however, an apparently severe counterexample has been known. Two model glass-forming liquids, with the standard Lennard-Jones interaction potential and its Weeks-Chandler-Andersen variation without the attractive tail, exhibit very similar structures at the two-body level but drastically different dynamical behaviors in the supercooled states. Here we look at the liquid structure through a (many-body) structural order parameter Θ characterizing the packing capability of local particle arrangements. We show that the structures of these two systems seen by Θ are actually very different at a many-body level, but, quite surprisingly, the macroscopic structure (Θ)-dynamics (τ_{α}) relationships commonly follow a Vogel-Fulcher-Tammann-like function. Furthermore, the mutual information analysis reveals strong local structure-dynamics correlations. Therefore, we conclude that attractive interactions affect the liquid structure in a nonperturbative manner, but a general structural origin of slow dynamics holds for these systems.

Attractive versus truncated repulsive supercooled liquids: The dynamics is encoded in the pair correlation function

We compare glassy dynamics in two liquids that differ in the form of their interaction potentials. Both systems have the same repulsive interactions but one has also an attractive part in the potential. These two systems exhibit very different dynamics despite having nearly identical pair correlation functions. We demonstrate that a properly weighted integral of the pair correlation function, which amplifies the subtle differences between the two systems, correctly captures their dynamical differences. The weights are obtained from a standard machine learning algorithm.

Continuous time random walk concepts applied to extended mode coupling theory: a study of the Stokes-Einstein breakdown

In an attempt to extend the mode coupling theory (MCT) to lower temperatures, some years back an Unified theory was proposed which within the MCT framework incorporated the activated dynamics via the random first order transition theory (RFOT). The theory successfully showed that there is hopping induced diffusive dynamics and the modified MCT coupled to the activated motion continues till low temperatures. Here we show that the theory although successful in describing other properties of supercooled liquids is unable to capture the Stokes-Einstein breakdown. We then show using continuous time random work (CTRW) formalism that the Unified theory is equivalent to a CTRW dynamics in presence of two waiting time distributions. It is known from earlier work on CTRW that in such cases the total dynamics is dominated by the fast motion. This explains the failure of the Unified theory in predicting the SE breakdown as both the structural relaxation and the diffusion process are described by the comparatively fast MCT like dynamics. The study also predicts that other forms of extended MCT with Markovian hopping kernel will face a similar issue. We next modify the Unified theory by applying the concept of renewal theory, usually used in CTRW models where the distribution has a long tail. According to this theory the first jump given by the persistent time is slower than the subsequent jumps given by the exchange time. We first show that for systems with two waiting time distributions even when both the distributions are exponential the persistent time is larger than the exchange time. We also identify the persistent time with the slower activated process. The extended Unified theory can now explain the SE breakdown. In this extended theory at low temperatures the structural relaxation is described by the activated dynamics whereas the diffusion is primarily determined by the MCT like dynamics leading to a decoupling between them. We also calculate a dynamic lengthscale from the wavenumber dependence of the relaxation time. We find that this dynamic length scale grows faster than the static length scale.

Effective one-component model of binary mixture: molecular arrest induced by the spatially correlated stochastic dynamics

Spatially correlated noise (SCN), i.e. the thermal noise that affects neighbouring particles in a similar manner, is ubiquitous in soft matter systems. In this work, we apply the over-damped SCN-driven Langevin equations as an effective, one-component model of the dynamics in dense binary mixtures. We derive the thermodynamically consistent fluctuation-dissipation relation for SCN to show that it predicts the molecular arrest resembling the glass transition, i.e. the critical slow-down of dynamics in the disordered phases. We show that the mechanism of singular dissipation is embedded in the dissipation matrix, accompanying SCN. We are also able to identify the characteristic length of collective dissipation, which diverges at critical packing. This novel physical quantity conveniently describes the difference between the ergodic and non-ergodic dynamics. The model is fully analytically solvable, one-dimensional and admits arbitrary interactions between the particles. It qualitatively reproduces several different modes of arrested disorder encountered in binary mixtures, including e.g. the re-entrant arrest. The model can be effectively compared to the mode coupling theory.

Dynamic Heterogeneities in Colloidal Supercooled Liquids: Experimental Tests of Inhomogeneous Mode Coupling Theory

The dynamics in supercooled liquids slow enormously upon approaching the glass transition, albeit without significant change of liquid structure. This empirical observation has stimulated development of many theoretical models which attempt to elucidate microscopic mechanisms in glasses and glass precursors. Here, quasi-two-dimensional colloidal supercooled liquids and glasses are employed to experimentally test predictions of widely used models: mode coupling theory (MCT) and its important extension, inhomogeneous MCT (IMCT). We measure two-point dynamic correlation functions in the glass forming liquids to determine structural relaxation times, τ_α, and mode coupling exponents, a, b, and γ; these parameters are then used to extract the mode coupling dynamic crossover packing area-fraction, ϕ _c. This information, along with our measurements of supercooled liquid spatiotemporal dynamics, permits characterization of dynamic heterogeneities in the samples and facilitates direct experimental tests of the scaling predictions of IMCT. The time scales at which dynamic heterogeneities are largest, and their spatial sizes, exhibit power law growth on approaching ϕ _c. Within experimental error, the exponents of the measured power laws are close to the predictions of IMCT.

Role of the Pair Correlation Function in the Dynamical Transition Predicted by Mode Coupling Theory

In a recent study, we have found that for a large number of systems the configurational entropy at the pair level S_{c2}, which is primarily determined by the pair correlation function, vanishes at the dynamical transition temperature T_{c}. Thus, it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this, we describe the dynamics of the system at the mean field level and, using the concepts of the dynamical density functional theory, show that the dynamical transition temperature depends only on the pair correlation function. Thus, this theory is similar in spirit to the microscopic mode coupling theory (MCT). However, unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature that is similar to T_{c}. This implies that the information of the dynamical transition temperature is embedded in the pair correlation function.

Glass transition in hard-core fluids and beyond, using an effective static structure in the mode coupling theory

The dynamical arrest in classical fluids is studied using a simple modification of the mode coupling theory (MCT) aimed at correcting its overestimation of the tendency to glass formation while preserving its overall structure. As in previous attempts, the modification is based on the idea of tempering the static pair correlations used as input. It is implemented in this work by computing the static structure at a different state point than the one used to solve the MCT equation for the intermediate scattering function, using the pure hard-sphere glass for calibration. The location of the glass transition predicted from this modification is found to agree with simulations data for a variety of systems --pure fluids and mixtures with either purely repulsive interaction potentials or ones with attractive contributions. Besides improving the predictions in the long-time limit, and so reducing the non-ergodicity domain, the same modification works as well for the time-dependent correlators.

Probing heterogeneous dynamics from spatial density correlation in glass-forming liquids

We numerically investigate the connection between spatial density correlation and dynamical heterogeneity in glass-forming liquids. We demonstrate that the cluster size defined by the spatial aggregation of densely packed particles (DPPs) can better capture the difference between the dynamics of the Lennard-Jones glass model and the Weeks-Chandler-Andersen truncation model than the commonly used pair correlation functions. More interestingly, we compare the mobility of DPPs and loosely packed particles, and we find that high local density correlates well with slow dynamics in systems with relatively hard repulsive interactions but links to mobile ones in the system with soft repulsive interactions at one relaxation time scale. Our results show clear evidence that the above model dependence behavior stems from the hopping motion of DPPs at the end of the caging stage due to the compressive nature of soft repulsive spheres, which activates the dynamics of DPPs in the α relaxation stage.

Shear-induced solidification of athermal systems with weak attraction

We find that unjammed packings of frictionless particles with rather weak attraction can always be driven into solidlike states by shear. The structure of shear-driven solids evolves continuously with packing fraction from gel-like to jamminglike, but is almost independent of the shear stress. In contrast, both the density of vibrational states (DOVS) and force network evolve progressively with the shear stress. There exists a packing fraction independent shear stress σ_{c}, at which the shear-driven solids are isostatic and have a flattened DOVS. Solidlike states induced by a shear stress greater than σ_{c} possess properties of marginally jammed solids and are thus strictly defined shear jammed states. Below σ_{c}, shear-driven solids with rather different structures are all under isostaticity and share common features in the DOVS and force network. Our study leads to a jamming phase diagram for weakly attractive particles, which reveals the significance of the shear stress in determining properties of shear-driven solids and the connections and distinctions between jamminglike and gel-like states.

The nonequilibrium glassy dynamics of self-propelled particles

We study the glassy dynamics taking place in dense assemblies of athermal active particles that are driven solely by a nonequilibrium self-propulsion mechanism. Active forces are modeled as an Ornstein-Uhlenbeck stochastic process, characterized by a persistence time and an effective temperature, and particles interact via a Lennard-Jones potential that yields well-studied glassy behavior in the Brownian limit, which is obtained as the persistence time vanishes. By increasing the persistence time, the system departs more strongly from thermal equilibrium and we provide a comprehensive numerical analysis of the structure and dynamics of the resulting active fluid. Finite persistence times profoundly affect the static structure of the fluid and give rise to nonequilibrium velocity correlations that are absent in thermal systems. Despite these nonequilibrium features, for any value of the persistence time we observe a nonequilibrium glass transition as the effective temperature is decreased. Surprisingly, increasing departure from thermal equilibrium is found to promote (rather than suppress) the glassy dynamics. Overall, our results suggest that with increasing persistence time, microscopic properties of the active fluid change quantitatively, but the general features of the nonequilibrium glassy dynamics observed with decreasing the effective temperature remain qualitatively similar to those of thermal glass-formers.

Role of density modulation in the spatially resolved dynamics of strongly confined liquids

Confinement by walls usually produces a strong modulation in the density of dense liquids near the walls. Using molecular dynamics simulations, we examine the effects of the density modulation on the spatially resolved dynamics of a liquid confined between two parallel walls, using a resolution of a fraction of the interparticle distance in the liquid. The local dynamics is quantified by the relaxation time associated with the temporal autocorrelation function of the local density. We find that this local relaxation time varies in phase with the density modulation. The amplitude of the spatial modulation of the relaxation time can be quite large, depending on the characteristics of the wall and thermodynamic parameters of the liquid. To disentangle the effects of confinement and density modulation on the spatially resolved dynamics, we compare the dynamics of a confined liquid with that of an unconfined one in which a similar density modulation is induced by an external potential. We find several differences indicating that density modulation alone cannot account for all the features seen in the spatially resolved dynamics of confined liquids. We also examine how the dynamics near a wall depends on the separation between the two walls and show that the features seen in our simulations persist in the limit of large wall separation.

Unraveling the success and failure of mode coupling theory from consideration of entropy

We analyze the dynamics of model supercooled liquids in a temperature regime where predictions of mode coupling theory (MCT) are known to be valid qualitatively. In this regime, the Adam-Gibbs (AG) relation, based on an activation picture of dynamics, also describes the dynamics satisfactorily, and we explore the mutual consistency and interrelation of these descriptions. Although entropy and dynamics are related via phenomenological theories, the connection between MCT and entropy has not been argued for. In this work, we explore this connection and provide a microscopic derivation of the phenomenological Rosenfeld theory. At low temperatures, the overlap between the MCT power law regime and AG relation implies that the AG relation predicts an avoided divergence at Tc, the origin of which can be related to the vanishing of pair configurational entropy, which we find occurring at the same temperature. We also show that the residual multiparticle entropy plays an important role in describing the relaxation time.

Quantitative approximation schemes for glasses

By means of a systematic expansion around the infinite-dimensional solution, we obtain an approximation scheme to compute properties of glasses in low dimensions. The resulting equations take as input the thermodynamic and structural properties of the equilibrium liquid, and from this they allow one to compute properties of the glass. They are therefore similar in spirit to the Mode Coupling approximation scheme. Our scheme becomes exact, by construction, in dimension d→∞, and it can be improved systematically by adding more terms in the expansion.

Microscopic Theory for the Role of Attractive Forces in the Dynamics of Supercooled Liquids

We formulate a microscopic, no adjustable parameter, theory of activated relaxation in supercooled liquids directly in terms of the repulsive and attractive forces within the framework of pair correlations. Under isochoric conditions, attractive forces can nonperturbatively modify slow dynamics, but at high enough density their influence vanishes. Under isobaric conditions, attractive forces play a minor role. High temperature apparent Arrhenius behavior and density-temperature scaling are predicted. Our results are consistent with recent isochoric simulations and isobaric experiments on a deeply supercooled molecular liquid. The approach can be generalized to treat colloidal gelation and glass melting, and other soft matter slow dynamics problems.

Mode-coupling approach for the slow dynamics of a liquid on a spherical substrate

We study the dynamics of a one-component liquid constrained on a spherical substrate, a 2-sphere, and investigate how the mode-coupling theory (MCT) can describe the new features brought by the presence of curvature. To this end we have derived the MCT equations in a spherical geometry. We find that, as seen from the MCT, the slow dynamics of liquids in curved space at low temperature does not qualitatively differ from that of glass-forming liquids in Euclidean space. The MCT predicts the right trend for the evolution of the relaxation slowdown with curvature but is dramatically off at a quantitative level.

Intermittent dynamics and logarithmic domain growth during the spinodal decomposition of a glass-forming liquid

We use large-scale molecular dynamics simulations of a simple glass-forming system to investigate how its liquid-gas phase separation kinetics depends on temperature. A shallow quench leads to a fully demixed liquid-gas system whereas a deep quench makes the dense phase undergo a glass transition and become an amorphous solid. This glass has a gel-like bicontinuous structure that evolves very slowly with time and becomes fully arrested in the limit where thermal fluctuations become negligible. We show that the phase separation kinetics changes qualitatively with temperature, the microscopic dynamics evolving from a surface tension-driven diffusive motion at high temperature to a strongly intermittent, heterogeneous, and thermally activated dynamics at low temperature, with a logarithmically slow growth of the typical domain size. These results elucidate the microscopic mechanisms underlying a specific class of viscoelastic phase separation.

Simulated glass-forming polymer melts: dynamic scattering functions, chain length effects, and mode-coupling theory analysis

We present molecular-dynamics simulations for a fully flexible model of polymer melts with different chain length N ranging from short oligomers (N = 4) to values near the entanglement length (N = 64). For these systems we explore the structural relaxation of the supercooled melt near the critical temperature T c of mode-coupling theory (MCT). Coherent and incoherent scattering functions are analyzed in terms of the idealized MCT. For temperatures T > T c we provide evidence for the space-time factorization property of the β relaxation and for the time-temperature superposition principle (TTSP) of the α relaxation, and we also discuss deviations from these predictions for T ≈ T c. For T larger than the smallest temperature where the TTSP holds we perform a quantitative analysis of the dynamics with the asymptotic MCT predictions for the late β regime. Within MCT a key quantity, in addition to T c, is the exponent parameter λ. For the fully flexible polymer models studied we find that λ is independent of N and has a value (λ = 0.735 ) typical of simple glass-forming liquids. On the other hand, the critical temperature increases with chain length toward an asymptotic value T c (∞) . This increase can be described by T c (∞) - T c(N) ∼ 1/N and may be interpreted in terms of the N dependence of the monomer density ρ, if we assume that the MCT glass transition is ruled by a soft-sphere-like constant coupling parameter Γ c = ρ c T c (-1/4), where ρ c is the monomer density at T c. In addition, we also estimate T c from a Hansen-Verlet-like criterion and MCT calculations based on structural input from the simulation. For our polymer model both the Hansen-Verlet criterion and the MCT calculations suggest T c to decrease with increasing chain length, in contrast to the direct analysis of the simulation data.

Role of Structure and Entropy in Determining Differences in Dynamics for Glass Formers with Different Interaction Potentials

We present a study of two model liquids with different interaction potentials, exhibiting similar structure but significantly different dynamics at low temperatures. By evaluating the configurational entropy, we show that the differences in the dynamics of these systems can be understood in terms of their thermodynamic differences. Analyzing their structure, we demonstrate that differences in pair correlation functions between the two systems, through their contribution to the entropy, dominate the differences in their dynamics, and indeed overestimate the differences. Including the contribution of higher order structural correlations to the entropy leads to smaller estimates for the relaxation times, as well as smaller differences between the two studied systems.

Crossovers in the dynamics of supercooled liquids probed by an amorphous wall

We study the relaxation dynamics of a binary Lennard-Jones liquid in the presence of an amorphous wall generated from equilibrium particle configurations. In qualitative agreement with the results presented by Kob et al. [Nat. Phys. 8, 164 (2012).] for a liquid of harmonic spheres, we find that our binary mixture shows a saturation of the dynamical length scale close to the mode-coupling temperature T(c). Furthermore we show that, due to the broken symmetry imposed by the wall, signatures of an additional change in dynamics become apparent at a temperature well above T(c). We provide evidence that this modification in the relaxation dynamics occurs at a recently proposed dynamical crossover temperature T(s) > T(c), which is related to the breakdown of the Stokes-Einstein relation. We find that this dynamical crossover at T(s) is also observed for the harmonic spheres as well as a WCA liquid, showing that it may be a general feature of glass-forming systems.

Structural signatures of (two) characteristic dynamical temperatures in lithium metasilicate

We report on the dynamic and structural characterization of lithium metasilicate Li2SiO3, a network-forming ionic glass, by means of molecular dynamics simulations. The system is characterized by a network of SiO4 tetrahedra disrupted by Li ions which diffuse through the network. Measures of mean square displacement and the diffusion constant of Si and O atoms allow us to identify the mode-coupling temperature, Tc ≈ 1500 K. At a much lower temperature, a change in the slope of the specific volume versus temperature singles out the glass transition at Tg ≈ 1000 K, the temperature below which the system goes out of equilibrium. We find signatures of both dynamical temperatures in structural order parameters related to the orientation of the tetrahedra. At lower temperatures we find that a set of order parameters which measure the relative orientation of neighbouring tetrahedra cease to increase and stay constant below Tc. Nevertheless, the bond orientational order parameter, which in this system measures local tetrahedral order, is found to continue growing below Tc until Tg, below which it remains constant. Although these structural signatures of the two dynamical temperatures do not imply any real thermodynamic transition in terms of the order parameters, they do give insight into the relaxation processes that occur between Tc and Tg, in particular they allow us to characterize the nature of the crossover happening around Tc.

Theories of glass formation and the glass transition

This key-issues review is a plea for a new focus on simpler and more realistic models of glass-forming fluids. It seems to me that we have too often been led astray by sophisticated mathematical models that beautifully capture some of the most intriguing features of glassy behavior, but are too unrealistic to provide bases for predictive theories. As illustrations of what I mean, the first part of this article is devoted to brief summaries of imaginative, sensible, but disparate and often contradictory ideas for solving glass problems. Almost all of these ideas remain alive today, with their own enthusiastic advocates. I then describe numerical simulations, mostly by H Tanaka and coworkers, in which it appears that very simple, polydisperse systems of hard disks and spheres develop long range, Ising-like, bond-orientational order as they approach glass transitions. Finally, I summarize my recent proposal that topologically ordered clusters of particles, in disordered environments, tend to become aligned with each other as if they were two-state systems, and thus produce the observed Ising-like behavior. Neither Tanaka's results nor my proposed interpretation of them fit comfortably within any of the currently popular glass theories.

Probing the glass transition from structural and vibrational properties of zero-temperature glasses

We find that the density dependence of the glass transition temperature of Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) systems can be predicted from properties of the zero-temperature (T=0) glasses. Below a crossover density ρs, LJ and WCA glasses show different structures, leading to different vibrational properties and consequently making LJ glasses more stable with higher glass transition temperatures than WCA ones. Above ρs, structural and vibrational quantities of the T=0 glasses show scaling collapse. From scaling relations and dimensional analysis, we predict a density scaling of the glass transition temperature, in excellent agreement with simulation results. We also propose an empirical expression of the glass transition temperature using structural and vibrational properties of the T=0 glasses, which works well over a wide range of densities.

Yield stress in amorphous solids: a mode-coupling-theory analysis

The yield stress is a defining feature of amorphous materials which is difficult to analyze theoretically, because it stems from the strongly nonlinear response of an arrested solid to an applied deformation. Mode-coupling theory predicts the flow curves of materials undergoing a glass transition and thus offers predictions for the yield stress of amorphous solids. We use this approach to analyze several classes of disordered solids, using simple models of hard-sphere glasses, soft glasses, and metallic glasses for which the mode-coupling predictions can be directly compared to the outcome of numerical measurements. The theory correctly describes the emergence of a yield stress of entropic nature in hard-sphere glasses, and its rapid growth as density approaches random close packing at qualitative level. By contrast, the emergence of solid behavior in soft and metallic glasses, which originates from direct particle interactions is not well described by the theory. We show that similar shortcomings arise in the description of the caging dynamics of the glass phase at rest. We discuss the range of applicability of mode-coupling theory to understand the yield stress and nonlinear rheology of amorphous materials.

Dynamic equivalences in the hard-sphere dynamic universality class

We perform systematic simulation experiments on model systems with soft-sphere repulsive interactions to test the predicted dynamic equivalence between soft-sphere liquids with similar static structure. For this we compare the simulated dynamics (mean squared displacement, intermediate scattering function, α-relaxation time, etc.) of different soft-sphere systems, between them and with the hard-sphere liquid. We then show that the referred dynamic equivalence does not depend on the (Newtonian or Brownian) nature of the microscopic laws of motion of the constituent particles, and hence, applies independently to colloidal and to atomic simple liquids. Finally, we verify another more recently proposed dynamic equivalence, this time between the long-time dynamics of an atomic liquid and its corresponding Brownian fluid (i.e., the Brownian system with the same interaction potential).

Observation and characterization of the vestige of the jamming transition in a thermal three-dimensional system

We study the dependence on the packing fraction of the pair-correlation function g(r) and particle mobility in a dense three-dimensional packing of soft colloids made of poly N-isopropyl acrylamide (pNIPAM), a thermosensitive gel. We find that g(r) for our samples is qualitatively like that of a liquid at all packing fractions. There is a peak in g(1), the height of the first peak of g(r), as a function of the packing fraction. This peak is identified as a vestige, which remains at finite temperature, of the divergence found at the jamming transition in simulations of soft frictionless spheres at zero temperature. As the density is increased, the particle dynamics slow down and near the packing fraction where there is a peak in g(1) the particles become arrested on the time scale of the experiment.

Quantitative field theory of the glass transition

We develop a full microscopic replica field theory of the dynamical transition in glasses. By studying the soft modes that appear at the dynamical temperature, we obtain an effective theory for the critical fluctuations. This analysis leads to several results: we give expressions for the mean field critical exponents, and we analytically study the critical behavior of a set of four-points correlation functions, from which we can extract the dynamical correlation length. Finally, we can obtain a Ginzburg criterion that states the range of validity of our analysis. We compute all these quantities within the hypernetted chain approximation for the Gibbs free energy, and we find results that are consistent with numerical simulations.