Locally preferred structures and many-body static correlations in viscous liquids

Tuning higher order structure in colloidal fluids

Colloidal particles self assemble into a wide range of structures under external AC electric fields due to induced dipolar interactions [Yethiraj and Van Blaaderen, Nature, 2003, 421, 513]. As a result of these dipolar interactions, at low volume fraction the system is modulated between a hard-sphere like state (in the case of zero applied field) and a "string fluid" upon application of the field. Using both particle-resolved experiments and computer simulations, we investigate the emergence of the string fluid with a variety of structural measures including two-body and higher-order correlations. We probe the higher-order structure using three-body spatial correlation functions and a many-body approach based on minimum energy clusters of a dipolar-Lennard-Jones system. The latter constitutes a series of geometrically distinct minimum energy clusters upon increasing the strength of the dipolar interaction, which are echoed in the higher-order structure of the colloidal fluids we study here. We find good agreement between experiment and simulation at the two-body level. Higher-order correlations exhibit reasonable agreement between experiment and simulation, again with more discrepancy at higher field strength for three-body correlation functions. At higher field strength, the cluster population in our experiments and simulations is dominated by the minimum energy clusters for all sizes 8 ≤ m ≤ 12.

Collineations of particles in the Kob-Andersen system

Numerous indications suggest that subtle changes occurring in the structures of liquids on supercooling are connected to the phenomenon of the glass transition and that detailed understanding of these changes is crucial for the development of new glasses with desired properties. J. D. Bernal in his 1962 Bakerian lecture in particular reported about an observation of approximately linear chains of several particles, referred to as collineations. He found that in the studied hard sphere system, these collineations can contain up to eight particles. Since then the collineations of three particles have been discussed in many papers in the context of the splitting of the second peak in pair density functions of supercooled liquids and glasses. However, it appears that longer collineations involving more that three particles have not been systematically studied. Here we report on our study of such collineations for the Kob-Andersen system of particles on cooling for the parent and inherent structures. Contrary to intuition, our findings reveal that below the potential energy landscape crossover temperature, the number of collineations in the parent structures can exceed that of the corresponding inherent structures. We also introduce a model that connects long collineations with the pair density and angular density distribution functions and demonstrate that this model describes long collineations quite well. We also studied the diffusion of the particles forming collineations and demonstrated that these particles are slightly slower than the randomly chosen average particle. Preliminary investigation of the collineations' lifetimes suggests that there are vibrational and structural relaxation regimes in the decay of the considered correlation function. The second part of the paper explores potential connections between collineations and (1) the disclination lines associated with the geometric frustration approach, (2) low-energy clusters from the topological cluster classification approach, and (3) chainlike cooperative motion of particles in low-temperature supercooled liquids. For the studied system, according to the used methods, no clear connection was found between collineations and these phenomena.

Freezing, melting, and the onset of glassiness in binary mixtures

We clarify the relationship between freezing, melting, and the onset of glassy dynamics in a prototypical glass-forming mixture model. Our starting point is a precise operational definition of the onset of glassiness, as expressed by the emergence of inflections in time-dependent correlation functions. By scanning the temperature-composition phase diagram of the mixture, we find a disconnect between the onset of glassiness and freezing. Surprisingly, however, the onset temperature closely tracks the melting line, along which the excess entropy is approximately constant. At fixed composition, all characteristic temperatures display nonetheless similar pressure dependencies, which are very well predicted by the isomorph theory. While our results rule out a general connection between thermodynamic metastability and glassiness, they call for a reassessment of the role of crystalline precursors in glass-forming liquids.

Structural Origin of Dynamic Heterogeneity in Supercooled Liquids

As a liquid is supercooled toward the glass transition point, its dynamics slow significantly, provided that crystallization is avoided. With increased supercooling, the particle dynamics become more spatially heterogeneous, a phenomenon known as dynamic heterogeneity. Since its discovery, this characteristic of metastable supercooled liquids has garnered considerable attention in glass science. However, the precise physical origins of dynamic heterogeneity remain elusive and widely debated. In this perspective, we examine the relationship between dynamic heterogeneity and structural order, based on numerical simulations of fragile liquids with isotropic potentials and strong liquids with directional interactions. We demonstrate that angular ordering, arising from many-body steric interactions, plays a crucial role in the slow dynamics and dynamic cooperativity of fragile liquids. Additionally, we explore how the growth of static order correlates with slower dynamics. In fragile liquids exhibiting super-Arrhenius behavior, the spatial extent of regions with high angular order grows upon cooling, and the sequential propagation of particle rearrangements within these ordered regions increases the activation energy for particle motion. In contrast, strong liquids with spatially constrained local ordering display a distinct "two-state" dynamic characteristic, marked by a transition between two Arrhenius-type behaviors. We argue that dynamic heterogeneity, irrespective of a liquid's fragility, arises from underlying structural order, with its spatial extent determined by static ordering. This perspective aims to deepen our understanding of the interplay between structural and dynamic properties in metastable supercooled liquids.

Liquidus curve of a Lennard-Jones mixture

The temperature below which the homogeneous liquid state of a mixture is not thermodynamically stable is called the liquidus temperature. This temperature varies with composition, and its composition dependence is represented by the liquidus curve. This curve provides fundamental reference points on the composition-temperature plane for characterizing the behavior of liquids, glasses, and crystals. In this paper, using molecular dynamics simulations, we determine the liquidus curve of the Wahnström mixture, which consists of large and small atoms interacting via the Lennard-Jones potential. Since this system is one of the standard models used to study the behavior of liquids and glasses, the liquidus curve presented in this work will contribute to a deeper understanding of disordered materials in general.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Dynamics and Time Scales of Higher-Order Correlations in Supercooled Colloidal Systems

The dynamics and time scales of higher-order correlations are studied in supercooled colloidal systems. A combination of X-ray photon correlation spectroscopy (XPCS) and X-ray cross-correlation analysis (XCCA) shows the typical slowing of the dynamics of a hard sphere system when approaching the glass transition. The time scales of higher-order correlations are probed using a novel time correlation function g_C, tracking the time evolution of cross-correlation function C. With an increasing volume fraction, the ratio of relaxation times of g_C to the standard individual particle relaxation time obtained by XPCS increases from ∼0.4 to ∼0.9. While a value of ∼0.5 is expected for free diffusion, the increasing values suggest that the local orders within the sample are becoming more long-lived for larger volume fractions. Furthermore, the dynamics of local order is more heterogeneous than the individual particle dynamics. These results indicate that not only the presence but also the lifetime of locally favored structures increases close to the glass transition.

Static self-induced heterogeneity in glass-forming liquids: Overlap as a microscope

We propose and numerically implement a local probe of the static self-induced heterogeneity characterizing glass-forming liquids. This method relies on the equilibrium statistics of the overlap between pairs of configurations measured in mesoscopic cavities with unconstrained boundaries. By systematically changing the location of the probed cavity, we directly detect spatial variations of the overlap fluctuations. We provide a detailed analysis of the statistics of a local estimate of the configurational entropy, and we infer an estimate of the surface tension between amorphous states, ingredients that are both at the basis of the random first-order transition theory of glass formation. Our results represent the first direct attempt to visualize and quantify the self-induced heterogeneity underpinning the thermodynamics of glass formation. They pave the way for the development of coarse-grained effective theories and for a direct assessment of the role of thermodynamics in the activated dynamics of deeply supercooled liquids.

Amorphous Order and Nonlinear Susceptibilities in Glassy Materials

We review 15 years of theoretical and experimental work on the nonlinear response of glassy systems. We argue that an anomalous growth of the peak value of nonlinear susceptibilities is a signature of growing "amorphous order" in the system, with spin-glasses as a case in point. Experimental results on supercooled liquids are fully compatible with the random first-order transition (RFOT) prediction of compact "glassites" of increasing volume as temperature is decreased, or as the system ages. We clarify why such a behavior is hard to explain within purely kinetic theories of glass formation, despite recent claims to the contrary.

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.

Neglecting polydispersity degrades propensity measurements in supercooled liquids

We conduct molecular dynamics simulations of a bidisperse Kob-Andersen (KA) glass former, modified to add in additional polydispersity. The original KA system is known to exhibit dynamical heterogeneity. Prior work defined propensity, the mean motion of a particle averaged over simulations reconstructing the initial positions of all particles but with randomized velocities. The existence of propensity shows that structure and dynamics are connected. In this paper, we study systems which mimic problems that would be encountered in measuring propensity in a colloidal glass former, where particles are polydisperse (they have slight size variations). We mimic polydispersity by altering the bidisperse KA system into a quartet consisting of particles both slightly larger and slightly smaller than the parent particles in the original bidisperse system. We then introduce errors into the reconstruction of the initial positions that mimic mistakes one might make in a colloidal experiment. The mistakes degrade the propensity measurement, in some cases nearly completely; one no longer has an iso-configurational ensemble in any useful sense. Our results show that a polydisperse sample is suitable for propensity measurements provided one avoids reconstruction mistakes.

Isomorph Invariance of Higher-Order Structural Measures in Four Lennard-Jones Systems

In the condensed liquid phase, both single- and multicomponent Lennard-Jones (LJ) systems obey the "hidden-scale-invariance" symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank-Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob-Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.

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.

Structure of the simple harmonic-repulsive system in liquid and glassy states studied by the triple correlation function

An efficient description of the structures of liquids and, in particular, the structural changes that happen with liquids on supercooling remains to be a challenge. The systems composed of soft particles are especially interesting in this context because they often demonstrate non-trivial local orders that do not allow to introduce the concept of the nearest-neighbor shell. For this reason, the use of some methods, developed for the structure analysis of atomic liquids, is questionable for the soft-particle systems. Here we report about our investigations of the structure of the simple harmonic-repulsive liquid in 3D using the triple correlation function (TCF), i.e., the method that does not rely on the nearest neighbor concept. The liquid is considered at reduced pressure (P = 1.8) at which it exhibits remarkable stability against crystallization on cooling. It is demonstrated that the TCF allows addressing the development of the orientational correlations in the structures that do not allow drawing definite conclusions from the studies of the bond-orientational order parameters. Our results demonstrate that the orientational correlations, if measured by the heights of the peaks in the TCF, significantly increase on cooling. This rise in the orientational ordering is not captured properly by the Kirkwood's superposition approximation. Detailed considerations of the peaks' shapes in the TCF suggest the existence of a link between the orientational ordering and the slowdown of the system's dynamics. Our findings support the view that the development of the orientational correlations in liquids may play a significant role in the liquids' dynamics and that the considerations of the pair distribution function may not be sufficient to understand intuitively all the structural changes that happen with liquids on supercooling. In general, our results demonstrate that the considerations of the TCF are useful in the discussions of the liquid's structures beyond the pair density function and interpreting the results obtained with the bond-orientational order parameters.

Stable glassy configurations of the Kob-Andersen model using swap Monte Carlo

The swap Monte Carlo algorithm allows the preparation of highly stable glassy configurations for a number of glass-formers but is inefficient for some models, such as the much studied binary Kob-Andersen (KA) mixture. We have recently developed generalizations to the KA model where swap can be very effective. Here, we show that these models can, in turn, be used to considerably enhance the stability of glassy configurations in the original KA model at no computational cost. We successfully develop several numerical strategies both in and out of equilibrium to achieve this goal and show how to optimize them. We provide several physical measurements indicating that the proposed algorithms considerably enhance mechanical and thermodynamic stability in the KA model, including a transition toward brittle yielding behavior. Our results thus pave the way for future studies of stable glasses using the KA model.

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.

Elucidation of the Nature of Structural Relaxation in Glassy d-Sorbitol

The nature and origin of the glass transition is one of the great unsolved problems of condensed-matter science. With the rapid increase of viscosity upon cooling the liquid near the glass-transition temperature, a range of dynamical motifs are observed, revealing the sheer complexity of interactions between the amorphous units. Yet, the causal link between those motifs and the solidification process remains unclear. Here, we apply a novel approach for exploring nontrivial interactions between structural units in d-sorbitol, a canonical example of a hydrogen-bonded organic glass, by introducing a dihedral-rearrangement-indicator analysis to shed light on relaxation processes and dynamical heterogeneity, which are known for their association with the stability of a glass. We find that both α- and β-relaxation processes are governed by cooperative and heterogeneous changes in hydrogen-bond dynamics that can be described by spatial and dihedral-angle-rearrangement indicators. The methodology and findings are of general applicability to other glass-forming systems.

Role of Attractive Forces in the Relaxation Dynamics of Supercooled Liquids

The attractive tail of the intermolecular interaction affects very weakly the structural properties of liquids, while it affects dramatically their dynamical ones. Via the numerical simulations of model systems not prone to crystallization, both in three and in two spatial dimensions, here we demonstrate that the nonperturbative dynamical effects of the attractive forces are tantamount to a rescaling of the activation energy by the glass transition temperature T_{g}: systems only differing in their attractive interaction have the same structural and dynamical properties if compared at the same value of T/T_{g}.

Assessing the utility of structure in amorphous materials

This paper presents a set of general strategies for the analysis of structure in amorphous materials and a general approach to assessing the utility of any selected structural description. Two measures of structure are defined, "diversity" and "utility," and applied to two model glass forming binary atomic alloys, Cu₅₀Zr₅₀ and a Lennard-Jones A₈₀B₂₀ mixture. We show that the change in diversity associated with selecting Voronoi structures with high localization or low energy, while real, is too weak to support claims that specific structures are the prime cause of these local physical properties. In addition, a new structure-free measure of incipient crystal-like organization in mixtures is introduced, suitable for cases where the stable crystal is a compound structure.

Glass forming phase diagram and local structure of Kob-Andersen binary Lennard-Jones nanoparticles

Molecular dynamics simulation is used to study glass formation in Kob-Andersen binary Lennard-Jones nanoparticles and determine the glass forming phase diagram for the system as a function of composition. The radial distribution function, a Steinhardt bond-orientational order parameter, and favored local structure analysis are used to distinguish between glassy and ordered systems. We find that surface enrichment of the large atoms alters the nanoparticle core composition, leading to an overall shift of the glass forming region to lower small atom mole fractions, relative to the bulk system. At small atom mole fraction, x_B = 0.1, the nanoparticles form a solid with an amorphous core, enriched with the small atoms, surrounded by a partially ordered surface region, enriched with the large atom component. The most disordered glass nanoparticles occur at x_B ≈ 0.3, but the surface-core enrichment leads to the crystallization of the nanoparticle to the CsCl crystal above x_B ≈ 0.35, which is lower than observed in the bulk. The glass transition temperatures of the nanoparticles are also significantly reduced. This allows the liquid to remain dynamic to low temperatures and sample the low energy inherent structure minima on the potential energy surface containing a high abundance of favoured local structures.

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.

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

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

Role of structure in the α and β dynamics of a simple glass-forming liquid

The elusive connection between dynamics and local structure in supercooled liquids is an important piece of the puzzle in the unsolved problem of the glass transition. The Johari-Goldstein β relaxation, ubiquitous in glass-forming liquids, exhibits mean properties that are strongly correlated to the long-time α dynamics. However, the former comprises simpler, more localized motion, and thus has perhaps a more straightforward connection to structure. Molecular dynamics simulations were carried out on a two-dimensional, rigid diatomic molecule (the simplest structure exhibiting a distinct β process) to assess the role of the local liquid structure on both the Johari-Goldstein β and the α relaxation. Although the average properties for these two relaxations are correlated, there is no connection between the β and α properties of a given (single) molecule. The propensity for motion at long times is independent of the rate or strength of a molecule's β relaxation. The mobility of a molecule averaged over many initial energies, a measure of the influence of structure, was found to be heterogeneous, with clustering at both the β and α time scales. This heterogeneity is less extended spatially for the β than for the α dynamics, as expected; however, the local structure is the more dominant control parameter for the β process. In the glassy state, the arrangement of neighboring molecules determines entirely the relaxation properties, with no discernible effect from the particle momenta.

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.

Weak links between fast mobility and local structure in molecular and atomic liquids

We investigate by molecular-dynamics simulations, the fast mobility-the rattling amplitude of the particles temporarily trapped by the cage of the neighbors-in mildly supercooled states of dense molecular (linear trimers) and atomic (binary mixtures) liquids. The mixture particles interact by the Lennard-Jones potential. The non-bonded particles of the molecular system are coupled by the more general Mie potential with variable repulsive and attractive exponents in a range which is a characteristic of small n-alkanes and n-alcohols. Possible links between the fast mobility and the geometry of the cage (size and shape) are searched. The correlations on a per-particle basis are rather weak. Instead, if one groups either the particles in fast-mobility subsets or the cages in geometric subsets, the increase of the fast mobility with both the size and the asphericity of the cage is revealed. The observed correlations are weak and differ in states with equal relaxation time. Local forces between a tagged particle and the first-neighbour shell do not correlate with the fast mobility in the molecular liquid. It is concluded that the cage geometry alone is unable to provide a microscopic interpretation of the known, universal link between the fast mobility and the slow structural relaxation. We suggest that the particle fast dynamics is affected by regions beyond the first neighbours, thus supporting the presence of collective, extended fast modes.

Brownian dynamics: from glassy to trivial

We endow a system of interacting particles with two distinct, local, Markovian, and reversible microscopic dynamics that both converge to the Boltzmann-Gibbs equilibrium of standard liquids. While the first, standard, one leads to glassy dynamics, we use field-theoretical techniques to show that the latter displays no sign of glassiness. The approximations we use, akin to the mode-coupling approximation, are famous for magnifying glassy aspects of the dynamics, supposedly through the neglect of activated events. Despite this, the modified dynamics seem to stick to standard liquid relaxation. This finding singles out as applying to a realistic system of interacting particles in low dimensions and questions the role of the dynamical rules used to explore a given static free-energy landscape. Moreover, our peculiar choice of dynamical rules offers the possibility of a direct connection with replica theory, and our findings therefore call for a clarification of the interplay between replica theory and the underlying dynamics of the system.

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.

Correlation of local order with particle mobility in supercooled liquids is highly system dependent

We investigate the connection between local structure and dynamical heterogeneity in supercooled liquids. Through the study of four different models, we show that the correlation between a particle's mobility and the degree of local order in nearby regions is highly system dependent. Our results suggest that the correlation between local structure and dynamics is weak or absent in systems that conform well to the mean-field picture of glassy dynamics and strong in those that deviate from this paradigm. Finally, we investigate the role of order-agnostic point-to-set correlations and reveal that they provide similar information content to local structure measures, at least in the system where local order is most pronounced.

Investigating isomorphs with the topological cluster classification

Isomorphs are lines in the density-temperature plane of certain "strongly correlating" or "Roskilde simple" liquids where two-point structure and dynamics have been shown to be close to identical up to a scale transformation. Here we consider such a liquid, a Lennard-Jones glass former, and investigate the behavior along isomorphs of higher-order structural and dynamical correlations. We then consider an inverse power law reference system mapped to the Lennard-Jones system [Pedersen et al., Phys. Rev. Lett. 105, 157801 (2010)]. Using the topological cluster classification to identify higher-order structures, in both systems we find bicapped square antiprisms, which are known to be a locally favored structure in the Lennard-Jones glass former. The population of these locally favored structures is up to 80% higher in the Lennard-Jones system than the equivalent inverse power law system. The structural relaxation time of the two systems, on the other hand, is almost identical, and the four-point dynamical susceptibility is marginally higher in the inverse power law system. Upon cooling, the lifetime of the locally favored structures in the Lennard-Jones system is up to 40% higher relative to the reference system.

Structural evolution in the aging process of supercooled colloidal liquids

When a liquid is rapidly quenched to a temperature below the glass-transition point, it is driven out of equilibrium; it then slowly relaxes to a (quasi)equilibrium state. This slow relaxation process is called aging. By definition, any glasses are inevitably in the process of aging and actually slowly evolving with time. Thus the study of aging phenomena is of fundamental importance for understanding not only the nonequilibrium nature of the glass transition, but also the stability of glassy materials. Here we consider aging after a rather shallow quench, for which a system is still able to reach (metastable) equilibrium. By using polydisperse colloidal liquids as a model, we show the validity of dynamical scaling that there is only one relevant length scale not only for a quasiequilibrium supercooled state but also for a nonequilibrium process of aging, which is reminiscent of dynamical critical phenomena. Our finding indicates that the aging toward (metastable) equilibrium may be regarded as the growth process of critical-like fluctuations of static order associated with low-free-energy configurations, further suggesting that this ordering is the origin of cooperative slow dynamics in the systems studied. The generality of this statement for other glass-forming systems remains for a future study.

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.

Multiple ordering transitions in a liquid stabilized by low symmetry structures

We present a numerical study of a lattice model of a liquid characterized by a low-symmetry favored local structure. We find that the freezing point is depressed far enough to reveal an exotic liquid-liquid transition characterized by the appearance of an extended chirality, prior to freezing. The ordered liquid can be readily supercooled to zero temperature, as the combination of critical slowing down and competing crystal polymorphs results in a dramatically slow crystallization process. These results provide an explicit scenario by which the ordering of a liquid can proceed via an intermediate liquid-liquid transition, a scenario that may prove helpful in the analysis of low temperature liquids interacting by more realistic interactions.

Spatiotemporal heterogeneity of local free volumes in highly supercooled liquid

We discuss the spatiotemporal behavior of local density and its relation to dynamical heterogeneity in a highly supercooled liquid by using molecular dynamics simulations of a binary mixture with different particle sizes in two dimensions. To trace voids heterogeneously existing with lower local densities, which move along with the structural relaxation, we employ the minimum local density for each particle in a time window whose width is set along with the structural relaxation time. Particles subject to free volumes correspond well to the configuration rearranging region of dynamical heterogeneity. While the correlation length for dynamical heterogeneity grows with temperature decrease, no growth in the correlation length of heterogeneity in the minimum local density distribution takes place. A comparison of these results with those of normal mode analysis reveals that superpositions of lower-frequency soft modes extending over the free volumes exhibit spatial correlation with the broken bonds. This observation suggests a possibility that long-ranged vibration modes facilitate the interactions between fragile regions represented by free volumes, to induce dynamical correlations at a large scale.

Decorrelation of the static and dynamic length scales in hard-sphere glass formers

We show that, in the equilibrium phase of glass-forming hard-sphere fluids in three dimensions, the static length scales tentatively associated with the dynamical slowdown and the dynamical length characterizing spatial heterogeneities in the dynamics unambiguously decorrelate. The former grow at a much slower rate than the latter when density increases. This observation is valid for the dynamical range that is accessible to computer simulations, which roughly corresponds to that accessible in colloidal experiments. We also find that, in this same range, no one-to-one correspondence between relaxation time and point-to-set correlation length exists. These results point to the coexistence of several relaxation mechanisms in the dynamically accessible regime of three-dimensional hard-sphere glass formers.

Ultrastable glasses from in silico vapour deposition

Glasses are generally prepared by cooling from the liquid phase, and their properties depend on their thermal history. Recent experiments indicate that glasses prepared by vapour deposition onto a substrate can exhibit remarkable stability, and might correspond to equilibrium states that could hitherto be reached only by glasses aged for thousands of years. Here we create ultrastable glasses by means of a computer-simulation process that mimics physical vapour deposition. These stable glasses have, far below the conventional glass-transition temperature, the properties expected for the equilibrium supercooled liquid state, and optimal stability is attained when deposition occurs at the Kauzmann temperature. We also show that the glasses' extraordinary stability is associated with distinct structural motifs, in particular the abundance of regular Voronoi polyhedra and the relative lack of irregular polyhedra.

Using mutual information to measure order in model glass formers

Whether or not there is growing static order accompanying the dynamical heterogeneity and increasing relaxation times seen in glassy systems is a matter of dispute. An obstacle to resolving this issue is that the order is expected to be amorphous and so not amenable to simple order parameters. We use mutual information to provide a general measurement of order that is sensitive to multiparticle correlations. We apply this to two glass-forming systems (two-dimensional binary mixtures of hard disks with different size ratios to give varying amounts of hexatic order) and show that there is little growth of amorphous order in the system without crystalline order. In both cases we measure the dynamical length with a four-point correlation function and find that it increases significantly faster than the static lengths in the system as density is increased. We further show that we can recover the known scaling of the dynamic correlation length in a kinetically constrained model, the two-vacancy-assisted-hopping triangular lattice gas.