Hopping and the Stokes-Einstein relation breakdown in simple glass formers

Unraveling relationship between complex lifetimes and microscopic diffusion in deep eutectic solvents

Aqueous mixtures of deep eutectic solvents (DESs) have emerged as a subject of interest in recent years for their tailored physicochemical properties. However, a comprehensive understanding of water's multifaceted influence on the microscopic dynamics, including its impact on improved transport properties of the DES, remains elusive. Additionally, the diffusion mechanisms within DESs manifest heterogeneous behavior, intricately tied to the formation and dissociation kinetics of complexes and hydrogen bonds. Therefore, it is imperative to explore the intricate interplay between bond kinetics, diffusion mechanism, and dynamical heterogeneity. This work employs water as an agent to explore their relationships by studying various relaxation phenomena in a DES based on acetamide and lithium perchlorate over a wide range of water concentrations. Notably, acetamide exhibits Fickian yet non-Gaussian diffusion across all water concentrations with Fickian (τf) and Gaussian (τg) timescales following a power-law relationship, τg∝τfγ, γ ∼ 1.4. The strength of coupling between bond kinetics and different diffusion timescales is estimated through various power-law relationships. Notably, acetamide-water hydrogen bond lifetime is linked to diffusive timescales through a single power-law over the entire water concentration studied. However, the relationship between diffusive timescales and the lifetime of acetamide-lithium complexes shows a sharp transition in behavior at 20 wt. % water, reflecting a change from vehicular diffusion below this concentration to structural diffusion above it. Our findings emphasize the critical importance of understanding bond dynamics within DESs, as they closely correlate with and regulate the molecular diffusion processes within these systems.

Local vs. cooperative: Unraveling glass transition mechanisms with SEER

Which phenomenon slows down the dynamics in supercooled liquids and turns them into glasses is a long-standing question of condensed matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or "excitations" govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, Systematic Excitation ExtRaction, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.

Dynamics and fluctuations of minimally structured glass formers

The mean-field theory (MFT) of simple structural glasses, which is exact in the limit of infinite spatial dimensions, d→∞, offers theoretical insight as well as quantitative predictions about certain features of d=3 systems. In order to more systematically relate the behavior of physical systems to MFT, however, various finite-d effects need to be accounted for. Although some efforts along this direction have already been undertaken, theoretical and technical challenges hinder progress. A general approach to sidestep many of these difficulties consists of simulating minimally structured models whose behavior smoothly converges to that described by the MFT as d increases, so as to permit a controlled dimensional extrapolation. Using this approach, we here extract the small fluctuations around the dynamical MFT captured by a standard liquid-state observable, the non-Gaussian parameter α_{2}. The results provide insight into the physical origin of these fluctuations as well as a quantitative reference with which to compare observations for more realistic glass formers.

Understanding the glassy dynamics from melting temperatures in binary glass-forming liquids

It is natural to expect that small particles in binary mixtures move faster than large ones. However, in binary glass-forming liquids with soft-core particle interactions, we observe the counterintuitive dynamic reversal between large and small particles along with the increase of pressure by performing molecular dynamics simulations. The structural relaxation (dynamic heterogeneity) of small particles is faster (weaker) than large ones at low pressures, but becomes slower (stronger) above a crossover pressure. In contrast, this dynamic reversal never happens in glass-forming liquids with hard-core interactions. We find that the difference of the effective melting temperatures felt by large and small particles can be used to understand the dynamic reversal. In binary mixtures, we derive effective melting temperatures of large and small particles simply from the conversion of units and find that particles with a higher effective melting temperature usually undergo a slower and more heterogeneous relaxation. The presence (absence) of the dynamic reversal in soft-core (hard-core) systems is simply due to the non-monotonic (monotonic) behavior of the melting temperature as a function of pressure. Interestingly, by manipulating the relative softness between large and small particles, we obtain a special case of soft-core systems, in which large particles always have higher effective melting temperatures than small ones. As a result, the dynamic reversal is totally eliminated. Our work provides another piece of evidence of the underlying connections between the properties of non-equilibrium glass-formers and equilibrium crystal-formers.

Multilayered noise model for transport in complex environments

Transport in complex fluidic environments often exhibits transient subdiffusive dynamics accompanied by non-Gaussian probability density profiles featuring a nonmonotonic non-Gaussian parameter. Such properties cannot be adequately explained by the original theory of Brownian motion. Based on an extension of kinetic theory, this study introduces a chain of hierarchically coupled random walks approach that effectively captures all these intriguing characteristics. If the environment consists of a series of independent white noise sources, then the problem can be expressed as a system of hierarchically coupled Ornstein-Uhlenbech equations. Due to the linearity of the system, the most essential transport properties have a closed analytical form.

Different glassy characteristics are related to either caging or dynamical heterogeneity

Despite the enormous theoretical and application interests, a fundamental understanding of the glassy dynamics remains elusive. The static properties of glassy and ordinary liquids are similar, but their dynamics are dramatically different. What leads to this difference is the central puzzle of the field. Even the primary defining glassy characteristics, their implications, and if they are related to a single mechanism remain unclear. This lack of clarity is a severe hindrance to theoretical progress. Here, we combine analytical arguments and simulations of various systems in different dimensions and address these questions. Our results suggest that the myriad of glassy features are manifestations of two distinct mechanisms. Particle caging controls the mean, and coexisting slow- and fast-moving regions govern the distribution of particle displacements. All the other glassy characteristics are manifestations of these two mechanisms; thus, the Fickian yet non-Gaussian nature of glassy liquids is not surprising. We discover a crossover, from stretched exponential to a power law, in the behavior of the overlap function. This crossover is prominent in simulation data and forms the basis of our analyses. Our results have crucial implications on how the glassy dynamics data are analyzed, challenge some recent suggestions on the mechanisms governing glassy dynamics, and impose strict constraints that a correct theory of glasses must have.

Jamming, relaxation, and memory in a minimally structured glass former

Structural glasses form through various out-of-equilibrium processes, including temperature quenches, rapid compression (crunches), and shear. Although each of these processes should be formally understandable within the recently formulated dynamical mean-field theory (DMFT) of glasses, the numerical tools needed to solve the DMFT equations up to the relevant physical regime do not yet exist. In this context, numerical simulations of minimally structured (and therefore mean-field-like) model glass formers can aid the search for and understanding of such solutions, thanks to their ability to disentangle structural from dimensional effects. We study here the infinite-range Mari-Kurchan model under simple out-of-equilibrium processes, and we compare results with the random Lorentz gas [J. Phys. A 55, 334001 (2022)10.1088/1751-8121/ac7f06]. Because both models are mean-field-like and formally equivalent in the limit of infinite spatial dimensions, robust features are expected to appear in the DMFT as well. The comparison provides insight into temperature and density onsets, memory, as well as anomalous relaxation. This work also further enriches the algorithmic understanding of the jamming density.

Diffusion coefficient scaling of a free Brownian particle with velocity-dependent damping

An analytical expression for the diffusion coefficient D of a free Brownian particle with velocity-dependent damping γ(v) is derived from the Green-Kubo formula. A special case of damping that decreases monotonically with velocity is considered. At high temperature T, the diffusion coefficient is found to exhibit two scaling types: (i) for a power-law decrease of damping with the particle's kinetic energy, γ(v)∝1/v^{2α}, it scales as D∝T^{α+1}; (ii) for a Gaussian function γ(v), it diverges at temperatures above a critical value T_{c} and behaves as D∝1/sqrt[T_{c}-T] at T slightly below T_{c}. At T>T_{c}, the particle trajectory contains long flight events, which are not observed at T<T_{c} in case (ii) and at all temperatures in case (i).

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.

Collective dynamics in a glass-former with Mari-Kurchan interactions

We numerically study the equilibrium relaxation dynamics of a two-dimensional Mari-Kurchan glass model. The tree-like structure of particle interactions forbids both non-trivial structural motifs and the emergence of a complex free-energy landscape leading to a thermodynamic glass transition, while the finite-dimensional nature of the model prevents the existence of a mode-coupling singularity. Nevertheless, the equilibrium relaxation dynamics is shown to be in excellent agreement with simulations performed in conventional glass-formers. Averaged time-correlation functions display a phenomenology typical of supercooled liquids, including the emergence of an excess signal in relaxation spectra at intermediate frequencies. We show that this evolution is accompanied by strong signatures of collective and heterogeneous dynamics which cannot be interpreted in terms of single particle hopping and emerge from dynamic facilitation. Our study demonstrates that an off-lattice interacting particle model with extremely simple structural correlations displays quantitatively realistic glassy dynamics.

Scaling and universality in Brownian motion on a stochastic harmonic oscillator chain

Diffusion of a Brownian particle along a stochastic harmonic oscillator chain is investigated. In contrast to the usually discussed Brownian motion driven by Gaussian white noise, the particle at high temperatures performs long Lévy flights. At high temperatures T the diffusion coefficient scales as D∼T^{2+α}, where the parameter α determine the average damping force ∝1/(T^{α}P) on the particle at large momentum P and at high temperature. The exponent α depends on the particle-chain interaction and chain properties. It is shown that the mean time t[over ¯]_{f} necessary to perform a flight of l lattice constant scales with l as t[over ¯]_{f}∝l^{2/3} at high temperatures and flight lengths. Last, the flight length probability distribution is found to decay as 1/l^{β} with the exponent β=4/3 being universal, i.e., independent of the model parameters.

A Relation between the Formation of a Hydrogen-Bond Network and a Time-Scale Separation of Translation and Rotation in Molecular Liquids

We study the relation between the translational and rotational motions of liquids, which is anticipated in the framework of the Stokes-Einstein-Debye (SED) treatment. For this purpose, we exploit the fact that ¹H field-cycling nuclear magnetic resonance relaxometry and molecular dynamics simulations provide access to both modes of motion. The experimental and computational findings are fully consistent and show that the time-scale separation between translation and rotation increases from the van der Waals liquid o-terphenyl over ethylene glycol to the hydrogen-bonded liquid glycerol, indicating an increasing degree of breakdown of the SED relation. Furthermore, the simulation results for two ethylene glycol models with different molecular conformations indicate that the translation is more retarded than the rotation when the density of intermolecular hydrogen bonds increases. We conclude that an increasing connectivity of a hydrogen-bond network leads to an increasing time-scale separation and, thus, to a stronger SED violation.

Efficient d-dimensional molecular dynamics simulations for studies of the glass-jamming transition

We develop an algorithm suitable for parallel molecular dynamics simulations in d spatial dimensions and describe its implementation in C++. All routines work in arbitrary d; the maximum simulated d is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass-jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. The scalings of simulation runtimes with the number of particles N and number of simulation threads n_{threads} are comparable to popular molecular dynamics codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled d=6 liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.

Local Dynamical Heterogeneity in Simple Glass Formers

We study the local dynamical fluctuations in glass-forming models of particles embedded in d-dimensional space, in the mean-field limit of d→∞. Our analytical calculation reveals that single-particle observables, such as squared particle displacements, display divergent fluctuations around the dynamical (or mode-coupling) transition, due to the emergence of nontrivial correlations between displacements along different directions. This effect notably gives rise to a divergent non-Gaussian parameter, α_{2}. The d→∞ local dynamics therefore becomes quite rich upon approaching the glass transition. The finite-d remnant of this phenomenon further provides a long sought-after, first-principle explanation for the growth of α_{2} around the glass transition that is not based on multiparticle correlations.

Fickian Non-Gaussian Diffusion in Glass-Forming Liquids

Fickian yet non-Gaussian diffusion (FnGD), a most intriguing open issue in soft matter, is generically associated with some dynamical and/or structural heterogeneity of the environment. Here we investigate the features of FnGD in glass-forming liquids, the epitome of dynamical heterogeneity, drawing on experiments on hard-sphere colloidal suspensions and simulations of a simple model of molecular liquid. We demonstrate that FnGD strengthens on approaching the glass transition, by identifying distinct timescales for Fickianity, τ_{F}, and for restoring of Gaussianity, τ_{G}>τ_{F}, as well as their associated length scales, ξ_{F} and ξ_{G}. We find τ_{G}∝τ_{F}^{γ} with γ≃1.8 for both systems. In the deep FnGD regime, the displacement distributions display exponential tails. We show that, in simulations, the time-dependent decay lengths l(t) at different temperatures all collapse onto a power-law master curve [l(t)/(ξ_{G})]∝(t/τ_{G})^{α}, with α=0.33. A similar collapse, if less sharp, is also found in experiments, seemingly with the same exponent α. We further discuss the connections of the timescales and length scales characterizing FnGD with structural relaxation and dynamic heterogeneity.

The dimensional evolution of structure and dynamics in hard sphere liquids

The formulation of the mean-field infinite-dimensional solution of hard sphere glasses is a significant milestone for theoretical physics. How relevant this description might be for understanding low-dimensional glass-forming liquids, however, remains unclear. These liquids indeed exhibit a complex interplay between structure and dynamics, and the importance of this interplay might only slowly diminish as dimension d increases. A careful numerical assessment of the matter has long been hindered by the exponential increase in computational costs with d. By revisiting a once common simulation technique involving the use of periodic boundary conditions modeled on D_d lattices, we here partly sidestep this difficulty, thus allowing the study of hard sphere liquids up to d = 13. Parallel efforts by Mangeat and Zamponi [Phys. Rev. E 93, 012609 (2016)] have expanded the mean-field description of glasses to finite d by leveraging the standard liquid-state theory and, thus, help bridge the gap from the other direction. The relatively smooth evolution of both the structure and dynamics across the d gap allows us to relate the two approaches and to identify some of the missing features that a finite-d theory of glasses might hope to include to achieve near quantitative agreement.

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.

Searching for structural predictors of plasticity in dense active packings

In amorphous solids subject to shear or thermal excitation, so-called structural indicators have been developed that predict locations of future plasticity or particle rearrangements. An open question is whether similar tools can be used in dense active materials, but a challenge is that under most circumstances, active systems do not possess well-defined solid reference configurations. We develop a computational model for a dense active crowd attracted to a point of interest, which does permit a mechanically stable reference state in the limit of infinitely persistent motion. Previous work on a similar system suggested that the collective motion of crowds could be predicted by inverting a matrix of time-averaged two-particle correlation functions. Seeking a first-principles understanding of this result, we demonstrate that this active matter system maps directly onto a granular packing in the presence of an external potential, and extend an existing structural indicator based on linear response to predict plasticity in the presence of noisy dynamics. We find that the strong pressure gradient necessitated by the directed activity, as well as a self-generated free boundary, strongly impact the linear response of the system. In low-pressure regions the linear-response-based indicator is predictive, but it does not work well in the high-pressure interior of our active packings. Our findings motivate and inform future work that could better formulate structure-dynamics predictions in systems with strong pressure gradients.

Configurational entropy from a replica approach: A density-functional model

We study a field-theoretic model for the metastable liquid using a nonlocal free-energy functional with density ρ(x) is the order parameter and three-point correlation effects included in the formulation. We assume fragmentation of the free-energy landscape into distinct basins of local minima and evaluate the partition function for the many-particle system through mapping into a composite system of m identical replicas. Static correlations and configurational entropy S_{c} are calculated in the m=1 limit. The Kauzman packing fraction η_{K} obtained are in agreement with other works.

Breakdown of the Stokes-Einstein relation in supercooled liquids: A cage-jump perspective

The breakdown of the Stokes-Einstein relation in supercooled liquids, which is the increase in the ratio τ_ατ_D between the two macroscopic times for structural relaxation and diffusion on decreasing the temperature, is commonly ascribed to dynamic heterogeneities, but a clear-cut microscopic interpretation is still lacking. Here, we tackle this issue exploiting the single-particle cage-jump framework to analyze molecular dynamics simulations of soft disk assemblies and supercooled water. We find that τ_ατ_D∝⟨t_p⟩⟨t_c⟩, where ⟨t_p⟩ and ⟨t_c⟩ are the cage-jump times characterizing slow and fast particles, respectively. We further clarify that this scaling does not arise from a simple term-by-term proportionality; rather, the relations τ_α∝⟨t_p⟩⟨Δr_J ²⟩ and τ_D∝⟨t_c⟩⟨Δr_J ²⟩ effectively connect the macroscopic and microscopic timescales, with the mean square jump length ⟨Δr_J ²⟩ shrinking on cooling. Our work provides a microscopic perspective on the Stokes-Einstein breakdown and generalizes previous results on lattice models to the case of more realistic glass-formers.

High-dimensional percolation criticality and hints of mean-field-like caging of the random Lorentz gas

The random Lorentz gas (RLG) is a minimal model for transport in disordered media. Despite the broad relevance of the model, theoretical grasp over its properties remains weak. For instance, the scaling with dimension d of its localization transition at the void percolation threshold is not well controlled analytically nor computationally. A recent study [Biroli et al., Phys. Rev. E 103, L030104 (2021)2470-004510.1103/PhysRevE.103.L030104] of the caging behavior of the RLG motivated by the mean-field theory of glasses has uncovered physical inconsistencies in that scaling that heighten the need for guidance. Here we first extend analytical expectations for asymptotic high-d bounds on the void percolation threshold and then computationally evaluate both the threshold and its criticality in various d. In high-d systems, we observe that the standard percolation physics is complemented by a dynamical slowdown of the tracer dynamics reminiscent of mean-field caging. A simple modification of the RLG is found to bring the interplay between percolation and mean-field-like caging down to d=3.

Spatial Dimensionality Dependence of Heterogeneity, Breakdown of the Stokes-Einstein Relation, and Fragility of a Model Glass-Forming Liquid

We investigate the heterogeneity of dynamics, the breakdown of the Stokes-Einstein relation and fragility in a model glass forming liquid, a binary mixture of soft spheres with a harmonic interaction potential for spatial dimensions from 3 to 8. The dynamical heterogeneity is quantified through the dynamical susceptibility χ₄ and the non-Gaussian parameter α₂. We find that the fragility, the degree of breakdown of the Stokes-Einstein relation, and the heterogeneity of the dynamics decrease with increasing spatial dimensionality. We briefly describe the dependence of fragility on the density and use it to resolve an apparent inconsistency with previous results.

Thermodynamic stability of hard sphere crystals in dimensions 3 through 10

Although much is known about the metastable liquid branch of hard spheres-from low dimension d up to [Formula: see text]-its crystal counterpart remains largely unexplored for [Formula: see text]. In particular, it is unclear whether the crystal phase is thermodynamically stable in high dimensions and thus whether a mean-field theory of crystals can ever be exact. In order to determine the stability range of hard sphere crystals, their equation of state is here estimated from numerical simulations, and fluid-crystal coexistence conditions are determined using a generalized Frenkel-Ladd scheme to compute absolute crystal free energies. The results show that the crystal phase is stable at least up to [Formula: see text], and the dimensional trends suggest that crystal stability likely persists well beyond that point.

Finite-size effects in the microscopic critical properties of jammed configurations: A comprehensive study of the effects of different types of disorder

Jamming criticality defines a universality class that includes systems as diverse as glasses, colloids, foams, amorphous solids, constraint satisfaction problems, neural networks, etc. A particularly interesting feature of this class is that small interparticle forces (f) and gaps (h) are distributed according to nontrivial power laws. A recently developed mean-field (MF) theory predicts the characteristic exponents of these distributions in the limit of very high spatial dimension, d→∞ and, remarkably, their values seemingly agree with numerical estimates in physically relevant dimensions, d=2 and 3. These exponents are further connected through a pair of inequalities derived from stability conditions, and both theoretical predictions and previous numerical investigations suggest that these inequalities are saturated. Systems at the jamming point are thus only marginally stable. Despite the key physical role played by these exponents, their systematic evaluation has yet to be attempted. Here, we carefully test their value by analyzing the finite-size scaling of the distributions of f and h for various particle-based models for jamming. Both dimension and the direction of approach to the jamming point are also considered. We show that, in all models, finite-size effects are much more pronounced in the distribution of h than in that of f. We thus conclude that gaps are correlated over considerably longer scales than forces. Additionally, remarkable agreement with MF predictions is obtained in all but one model, namely near-crystalline packings. Our results thus help to better delineate the domain of the jamming universality class. We furthermore uncover a secondary linear regime in the distribution tails of both f and h. This surprisingly robust feature is understood to follow from the (near) isostaticity of our configurations.

Mean-Field Caging in a Random Lorentz Gas

The random Lorentz gas (RLG) is a minimal model of both percolation and glassiness, which leads to a paradox in the infinite-dimensional, d → ∞ limit: the localization transition is then expected to be continuous for the former and discontinuous for the latter. As a putative resolution, we have recently suggested that, as d increases, the behavior of the RLG converges to the glassy description and that percolation physics is recovered thanks to finite-d perturbative and nonperturbative (instantonic) corrections [Biroli et al. Phys. Rev. E 2021, 103, L030104]. Here, we expand on the d → ∞ physics by considering a simpler static solution as well as the dynamical solution of the RLG. Comparing the 1/d correction of this solution with numerical results reveals that even perturbative corrections fall out of reach of existing theoretical descriptions. Comparing the dynamical solution with the mode-coupling theory (MCT) results further reveals that, although key quantitative features of MCT are far off the mark, it does properly capture the discontinuous nature of the d → ∞ RLG. These insights help chart a path toward a complete description of finite-dimensional glasses.

Interplay between percolation and glassiness in the random Lorentz gas

The random Lorentz gas (RLG) is a minimal model of transport in heterogeneous media that exhibits a continuous localization transition controlled by void space percolation. The RLG also provides a toy model of particle caging, which is known to be relevant for describing the discontinuous dynamical transition of glasses. In order to clarify the interplay between the seemingly incompatible percolation and caging descriptions of the RLG, we consider its exact mean-field solution in the infinite-dimensional d→∞ limit and perform numerics in d=2...20. We find that for sufficiently high d the mean-field caging transition precedes and prevents the percolation transition, which only happens on timescales diverging with d. We further show that activated processes related to rare cage escapes destroy the glass transition in finite dimensions, leading to a rich interplay between glassiness and percolation physics. This advance suggests that the RLG can be used as a toy model to develop a first-principle description of particle hopping in structural glasses.

Memory Formation in Jammed Hard Spheres

Liquids equilibrated below an onset condition share similar inherent states, while those above that onset have inherent states that markedly differ. Although this type of materials memory was first reported in simulations over 20 years ago, its physical origin remains controversial. Its absence from mean-field descriptions, in particular, has long cast doubt on its thermodynamic relevance. Motivated by a recent theoretical proposal, we reassess the onset phenomenology in simulations using a fast hard sphere jamming algorithm and find it to be both thermodynamically and dimensionally robust. Remarkably, we also uncover a second type of memory associated with a Gardner-like regime of the jamming algorithm.

Glassy dynamics of a binary Voronoi fluid: a mode-coupling analysis

The binary Voronoi mixture is a fluid model whose interactions are derived from the Voronoi-Laguerre tessellation of the configurations of the system. The resulting interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent coherent and incoherent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter λ ≈ 0.746 comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good qualitative agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob-Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.

Violation of Stokes-Einstein and Stokes-Einstein-Debye relations in polymers at the gas-supercooled liquid coexistence

Molecular dynamics simulations are performed on a system of model linear polymers to look at the violations of Stokes-Einstein (SE) and Stokes-Einstein-Debye (SED) relations near the mode coupling theory transition temperatureT_cat three (one higher and two lower) densities. At low temperatures, both lower density systems show stable gas-supercooled-liquid coexistence whereas the higher density system is homogeneous. We show that monomer density relaxation exhibits SE violation for all three densities, whereas molecular density relaxation shows a weak violation of the SE relation nearT_cin both lower density systems. This study identifies disparity in monomer mobility and observation of jumplike motion in the typical monomer trajectories resulting in the SE violations. In addition to the SE violation, a weak SED violation is observed in the gas-supercooled-liquid coexisting domains of the lower densities. Both lower density systems also show a decoupling of translational and rotational dynamics in this polymer system.

Breakdown of the Stokes-Einstein Equation for Solutions of Water in Oil Reverse Micelles

An experimental study is presented for the reverse micellar system of 15% by mass polydisperse hexaethylene glycol monodecylether (C₁₀E₆) in cyclohexane with varying amounts of added water up to 4% by mass. Measurements of viscosity and self-diffusion coefficients were taken as a function of temperature between 10 and 45 °C at varying sample water loads but fixed C₁₀E₆/cyclohexane composition. The results were used to inspect the validity of the Stokes-Einstein equation for this system. Unreasonably small reverse average micelle radii and aggregation numbers were obtained with the Stokes-Einstein equation, but reasonable values for these quantities were obtained using the ratio of surfactant-to-cyclohexane self-diffusion coefficients. While bulk viscosity increased with increasing water load, a concurrent expected decrease of self-diffusion coefficient was only observed for the surfactant and water but not for cyclohexane, which showed independence of water load. Moreover, a spread of self-diffusion coefficients was observed for the protons associated with the ethylene oxide repeat unit in samples with polydisperse C₁₀E₆ but not in a sample with monodisperse C₁₀E₆. These findings were interpreted by the presence of reverse micelle to reverse micelle hopping motions that with higher water load become increasingly selective toward C₁₀E₆ molecules with short ethylene oxide repeat units, while those with long ethylene oxide repeat units remain trapped within the reverse micelle because of the increased hydrogen bonding interactions with the water inside the growing core of the reverse micelle. Despite the observed breakdown of the Stokes-Einstein equation, the temperature dependence of the viscosities and self-diffusion coefficients was found to follow Arrhenius behavior over the investigated range of temperatures.

Deciphering the Dynamics of Organic Nanoaggregates with AIEE Effect and Excited States: Lipophilic Benzothiadiazole Derivatives as Selective Cell Imaging Probes

An aggregation-induced emission enhancement (AIEE) effect in fluorescent lipophilic 2,1,3-benzothiadiazole (BTD) derivatives and their organic nanoaggregates were studied. A set of techniques such as single-crystal X-ray, dynamic light scattering (DLS), electron paramagnetic resonance (EPR), UV-vis, fluorescence, and density functional theory (DFT) calculations have been used to decipher the formation/break (kinetics), properties, and dynamics of the organic nanoaggregates of three BTD small organic molecules. An in-depth study of the excited-state also revealed the preferential relaxation emissive pathways for the BTD derivatives and the dynamics associated with it. The results described herein, for the first time, explain the formation of fluorescent BTD nanoaggregate derivatives and allow for the understanding of their dynamics in solution as well as the ruling forces of both aggregation and break processes along with the involved equilibrium. One of the developed dyes could be used at a nanomolar concentration to selectively stain lipid droplets emitting an intense and bright fluorescence at the red channel. The other two BTDs could also stain lipid droplets at very low concentrations and were visualized preferentially at the blue channel.

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.

Excess-entropy scaling in supercooled binary mixtures

Transport coefficients, such as viscosity or diffusion coefficient, show significant dependence on density or temperature near the glass transition. Although several theories have been proposed for explaining this dynamical slowdown, the origin remains to date elusive. We apply here an excess-entropy scaling strategy using molecular dynamics computer simulations and find a quasiuniversal, almost composition-independent, relation for binary mixtures, extending eight orders of magnitude in viscosity or diffusion coefficient. Metallic alloys are also well captured by this relation. The excess-entropy scaling predicts a quasiuniversal breakdown of the Stokes-Einstein relation between viscosity and diffusion coefficient in the supercooled regime. Additionally, we find evidence that quasiuniversality extends beyond binary mixtures, and that the origin is difficult to explain using existing arguments for single-component quasiuniversality.

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.

Numerical solution of the dynamical mean field theory of infinite-dimensional equilibrium liquids

We present a numerical solution of the dynamical mean field theory of infinite-dimensional equilibrium liquids established by Maimbourg et al. [Phys. Rev. Lett. 116, 015902 (2016)]. For soft sphere interactions, we obtain the numerical solution by an iterative algorithm and a straightforward discretization of time. We also discuss the case of hard spheres for which we first derive analytically the dynamical mean field theory as a non-trivial limit of that of soft spheres. We present numerical results for the memory function and the mean square displacement. Our results reproduce and extend kinetic theory in the dilute or short-time limit, while they also describe dynamical arrest toward the glass phase in the dense strongly interacting regime.

Nonempirical Free Volume Viscosity Model for Alkane Lubricants under Severe Pressures

Viscosities η and diffusion coefficients D_{s} of linear and branched alkanes at pressure 0<P<0.7  GPa and temperature T=500-600  K are calculated from molecular dynamics simulations. Combining Stokes-Einstein, free volume, and random walk concepts results in an accurate viscosity model for the considered P and T. All model parameters (hydrodynamic radius, random walk step size, and step frequency) are extracted from equilibrium molecular dynamics via microscopic ensemble averages rendering η(P,T) a parameter-free predictor for lubrication simulations.

Insights into the Emerging Networks of Voids in Simulated Supercooled Water

The structural evolution of supercooled liquid water as we approach the glass transition temperature continues to be an active area of research. Here, we use molecular dynamics simulations of TIP4P/ice water to study the changes in the connected regions of empty space within the liquid, which we investigate using the Voronoi-voids network. We observe two important features: supercooling enhances the fraction of nonspherical voids and different sizes of voids tend to cluster forming a percolating network. By examining order parameters such as the local structure index (LSI), tetrahedrality and topological defects, we show that water molecules near large void clusters tend to be slightly more tetrahedral than those near small voids, with a lower population of under- and overcoordinated defects. We show further that the distribution of closed rings of water molecules around small and large void clusters maintain a balance between 6 and 7 membered rings. Our results highlight the changes of the dual voids and water network as a structural hallmark of supercooling and provide insights into the molecular origins of cooperative effects underlying density fluctuations on the subnanometer and nanometer length scale. In addition, the percolation of the voids and the hydrogen bond network around the voids may serve as useful order parameters to investigate density fluctuations in supercooled water.

Two-step relaxation and the breakdown of the Stokes-Einstein relation in glass-forming liquids

It is well known that glass-forming liquids exhibit a number of anomalous dynamical phenomena, most notably a two-step relaxation in the self-intermediate scattering function and the breakdown of the Stokes-Einstein (SE) relation, as they are cooled toward the glass transition temperature. While these phenomena are generally ascribed to dynamic heterogeneity, specifically to the presence of slow- and fast-moving particles, a quantitative elucidation of the two-step relaxation and the violation of the SE relation in terms of these concepts has not been successful. In this work, we propose a classification of particles according to the rank order of their displacements (from an arbitrarily defined origin of time), and we divide the particles into long-distance (LD), medium-distance, and short-distance (SD) traveling particle groups. Using molecular-dynamics simulation data of the Kob-Andersen model, we show quantitatively that the LD group is responsible for the fast relaxation in the two-step relaxation process in the intermediate scattering function, while the SD group gives rise to the slow (α) relaxation. Furthermore, our analysis reveals that τ_{α} is controlled by the SD group, while the ensemble-averaged diffusion coefficient D is controlled by both the LD and SD groups. The combination of these two features provides a natural explanation for the breakdown in the SE relation at low temperature. In addition, we find that the α-relaxation time, τ_{α}, of the overall system is related to the relaxation time of the LD particles, τ_{LD}, as τ_{α}=τ_{0}exp(Ωτ_{LD}/k_{B}T).

Active glasses

Active glassy matter has recently emerged as a novel class of non-equilibrium soft matter, combining energy-driven, active particle movement with dense and disordered glass-like behavior. Here we review the state-of-the-art in this field from an experimental, numerical, and theoretical perspective. We consider both non-living and living active glassy systems, and discuss how several hallmarks of glassy dynamics (dynamical slowdown, fragility, dynamical heterogeneity, violation of the Stokes-Einstein relation, and aging) are manifested in such materials. We start by reviewing the recent experimental evidence in this area of research, followed by an overview of the main numerical simulation studies and physical theories of active glassy matter. We conclude by outlining several open questions and possible directions for future work.

Understanding, predicting, and tuning the fragility of vitrimeric polymers

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

Gardner physics in amorphous solids and beyond

One of the most remarkable predictions to emerge out of the exact infinite-dimensional solution of the glass problem is the Gardner transition. Although this transition was first theoretically proposed a generation ago for certain mean-field spin glass models, its materials relevance was only realized when a systematic effort to relate glass formation and jamming was undertaken. A number of nontrivial physical signatures associated with the Gardner transition have since been considered in various areas, from models of structural glasses to constraint satisfaction problems. This perspective surveys these recent advances and discusses the novel research opportunities that arise from them.

Transport dynamics of complex fluids

Thermal motion in complex fluids is a complicated stochastic process but ubiquitously exhibits initial ballistic, intermediate subdiffusive, and long-time diffusive motion, unless interrupted. Despite its relevance to numerous dynamical processes of interest in modern science, a unified, quantitative understanding of thermal motion in complex fluids remains a challenging problem. Here, we present a transport equation and its solutions, which yield a unified quantitative explanation of the mean-square displacement (MSD), the non-Gaussian parameter (NGP), and the displacement distribution of complex fluids. In our approach, the environment-coupled diffusion kernel and its time correlation function (TCF) are the essential quantities that determine transport dynamics and characterize mobility fluctuation of complex fluids; their time profiles are directly extractable from a model-free analysis of the MSD and NGP or, with greater computational expense, from the two-point and four-point velocity autocorrelation functions. We construct a general, explicit model of the diffusion kernel, comprising one unbound-mode and multiple bound-mode components, which provides an excellent approximate description of transport dynamics of various complex fluidic systems such as supercooled water, colloidal beads diffusing on lipid tubes, and dense hard disk fluid. We also introduce the concepts of intrinsic disorder and extrinsic disorder that have distinct effects on transport dynamics and different dependencies on temperature and density. This work presents an unexplored direction for quantitative understanding of transport and transport-coupled processes in complex disordered media.

Coupling and decoupling between translational and rotational dynamics in supercooled monodisperse soft Janus particles

We perform dynamics simulations to investigate the translational and rotational glassy dynamics in a glass-forming liquid of monodisperse soft Janus particles. We find that, with decreasing temperature, the mean-square angular displacement shows no clear plateau in the caging region, in contrast with the apparent caging behavior of translational motion. By defining a reorientational mean-square angular displacement, the caging behavior of rotational motion can be recognized. On approaching the glass transition (decreasing temperature), the coupling between translational and rotational relaxation increases, while the coupling between translational and rotational diffusion decreases, whereas the coupling between translational and reorientational diffusion increases. The strong decoupling between translational and rotational diffusion is due to the suppressed translational mobility but promoted rotational mobility of soft Janus particles. We think that the low-T SE and SED decoupling is mainly attributed to hopping motion of soft Janus particles, whereas the high-T SE and SED decoupling is mainly attributed to collective cage motion of soft Janus particles. Our results demonstrate that interaction anisotropy has a critical effect on the translational and rotational dynamics of soft Janus particles.

Dynamics around the site percolation threshold on high-dimensional hypercubic lattices

Recent advances on the glass problem motivate reexamining classical models of percolation. Here we consider the displacement of an ant in a labyrinth near the percolation threshold on cubic lattices both below and above the upper critical dimension of simple percolation, d_{u}=6. Using theory and simulations, we consider the scaling regime and obtain that both caging and subdiffusion scale logarithmically for d≥d_{u}. The theoretical derivation, which considers Bethe lattices with generalized connectivity and a random graph model, confirms that logarithmic scalings should persist in the limit d→∞. The computational validation employs accelerated random walk simulations with a transfer-matrix description of diffusion to evaluate directly the dynamical critical exponents below d_{u} as well as their logarithmic scaling above d_{u}. Our numerical results improve various earlier estimates and are fully consistent with our theoretical predictions.

Bypassing sluggishness: SWAP algorithm and glassiness in high dimensions

The recent implementation of a swap Monte Carlo algorithm (SWAP) for polydisperse glass forming mixtures bypasses computational sluggishness and closes the gap between experimental and simulation timescales in physical dimensions d=2 and 3. Here, we consider suitably optimized systems in d=2,3,⋯,8 to obtain insights into the performance and underlying physics of SWAP. We show that the speedup obtained decays rapidly with increasing the dimension. SWAP nonetheless delays systematically the onset of the activated dynamics by an amount that remains seemingly finite in the limit d→∞. This shows that the glassy dynamics in high dimensions d>3 is now computationally accessible using SWAP, thus opening the door for the systematic consideration of finite-dimensional deviations from the mean-field description.

Determination of onset temperature from the entropy for fragile to strong liquids

In this paper, we establish a connection between the onset temperature of glassy dynamics with the change in the entropy for a wide range of model systems. We identify the crossing temperature of pair and excess entropies as the onset temperature. Below the onset temperature, the residual multiparticle entropy, the difference between excess and pair entropies, becomes positive. The positive entropy can be viewed as equivalent to the larger phase space exploration of the system. The new method of onset temperature prediction from entropy is less ambiguous, as it does not depend on any fitting parameter like the existing methods.