Critical-like behaviour of glass-forming liquids

Power laws in pressure-induced structural change of glasses

Many glasses exhibit fractional power law (FPL) between the mean atomic volume v_a and the first diffraction peak position q₁, i.e. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_{\mathrm{a}} \propto q_1^{ - d}$$\end{document}va∝q1−d with d ≃ 2.5 deviating from the space dimension D = 3, under compression or composition change. What structural change causes such FPL and whether the FPL and d are universal remain controversial. Here our simulations show that the FPL holds in both two- and three-dimensional glasses under compression when the particle interaction has two length scales which can induce nonuniform local deformations. The exponent d is not universal, but varies linearly with the deformable part of soft particles. In particular, we reveal an unexpected crossover regime with d > D from crystal behavior (d = D) to glass behavior (d < D). The results are explained by two types of bond deformation. We further discover FPLs in real space from the radial distribution functions, which correspond to the FPLs in reciprocal space.

A puzzle in metallic glass research is the existence of the fractional power law in reciprocal space, whilst its origin remains controversial. Zhang et al. show that nonuniform local deformations under compression induce this phenomenon and quantify the power law exponent at both two and three dimensions.

Structural Ordering in Liquid Gallium under Extreme Conditions

The atomic-scale structure, melting curve, and equation of state of liquid gallium has been measured to high pressure (p) and high temperature (T) up to 26 GPa and 900 K by in situ synchrotron x-ray diffraction. Ab initio molecular dynamics simulations up to 33.4 GPa and 1000 K are in excellent agreement with the experimental measurements, providing detailed insight at the level of pair distribution functions. The results reveal an absence of dimeric bonding in the liquid state and a continuous increase in average coordination number n[over ¯]_{Ga}^{Ga} from 10.4(2) at 0.1 GPa approaching ∼12 by 25 GPa. Topological cluster analysis of the simulation trajectories finds increasing fractions of fivefold symmetric and crystalline motifs at high p-T. Although the liquid progressively resembles a hard-sphere structure towards the melting curve, the deviation from this simple description remains large (≥40%) across all p-T space, with specific motifs of different geometries strongly correlating with low local two-body excess entropy at high p-T.

An alternative approach to evidence the structural conditioning in the dynamic slowdown in a polymer glass-former

Dynamic slowdown of liquids, leading to a breakdown of Arrhenius behavior of relaxation and Stokes-Einstein relationship (SER), as the glass transition is approached, is still not fully understood despite decades of study. They are usually associated to the emergence of dynamic heterogeneity, that is, regions or clusters of particles that have high or low mobilities. But the physical origin of these dynamic heterogeneity, and in particular, the question whether they have a structural origin or they are a purely dynamical phenomenon, is still under debate. In this work we study through molecular dynamics simulations in a polymer model the dynamic slowdown and the breakdown of SER, in connection with dynamic susceptibility calculated for an isoconfigurational ensemble, such that the effects of structure on dynamics can be discriminated. The onset of structure effects on dynamical behavior is found to be coincident with the onset of slow dynamics and SER breakdown.

Curie-Weiss behavior of liquid structure and ideal glass state

We present the results of a structural study of metallic alloy liquids from high temperature through the glass transition. We use high energy X-ray scattering and electro-static levitation in combination with molecular dynamics simulation and show that the height of the first peak of the structure function, S(Q) - 1, follows the Curie-Weiss law. The structural coherence length is proportional to the height of the first peak, and we suggest that its increase with cooling may be related to the rapid increase in viscosity. The Curie temperature is negative, implying an analogy with spin-glass. The Curie-Weiss behavior provides a pathway to an ideal glass state, a state with long-range correlation without lattice periodicity, which is characterized by highly diverse local structures, reminiscent of spin-glass.

Structural order as a genuine control parameter of dynamics in simple glass formers

Glass transition is characterised by drastic dynamical slowing down upon cooling, accompanied by growing spatial heterogeneity. Its rationalisation by subtle changes in the liquid structure has been long debated but remains elusive, due to intrinsic difficulty in detecting the underlying complex structural ordering. Here we report that structural order parameter characterising local packing capability can well describe the glassy dynamics not only macroscopically but also microscopically, no matter whether it is driven by temperature or density. A Vogel-Fulcher-Tammann (VFT)-like relation is universally identified between the structural relaxation time and the order parameter for supercooled liquids with isotropic interactions. More importantly, we find such an intriguing VFT-like relation to be statistically valid even at a particle level, between spatially coarse-grained structural order and microscopic particle-level dynamics. Such a unified description of glassy dynamics based solely on structural order is expected to contribute to the ultimate understanding of the long-standing glass-transition problem.

The glass-forming materials exhibit dynamical slowing down together with spatial heterogeneity at microscales, but their origin remains debated. Tong and Tanaka show that this phenomenon can be unified based on a structural order parameter capable of detecting subtle ordering in instantaneous liquid states.

Nonmonotonic behavior in dense assemblies of active colloids

We study experimentally a sediment of self-propelled Brownian particles with densities ranging from dilute to ergodic supercooled to nonergodic glass to nonergodic polycrystal. In a companion paper, we observe a nonmonotonic response to activity of relaxation of the nonergodic glass state: a dramatic slowdown when particles become weakly self-propelled, followed by a speedup at higher activities. Here we map ergodic supercooled states to standard passive glassy physics, provided a monotonic shift of the glass packing fraction and the replacement of the ambient temperature by the effective temperature. However, we show that this mapping fails beyond glass transition. This failure is responsible for the nonmonotonic response. Furthermore, we generalize our finding by examining the dynamical response of another class of nonergodic systems: polycrystals. We observe the same nonmonotonic response to activity. To explain this phenomenon, we measure the size of domains where particles move in the same direction. This size also shows a nonmonotonic response, with small lengths corresponding to slow relaxation. This suggests that the failure of the mapping of nonergodic active states to a passive situation is general and is linked to anisotropic relaxation mechanisms specific to active matter.

What does the instantaneous normal mode spectrum tell us about dynamical heterogeneity in glass-forming fluids?

We examine the instantaneous normal mode spectrum of model metallic and polymeric glass-forming liquids. We focus on the localized modes in the unstable part of the spectrum [unstable localized (UL) modes] and find that the particles making the dominant contribution to the participation ratio form clusters that grow upon cooling in a fashion similar to the dynamical heterogeneity in glass-forming fluids, i.e., highly mobile (or immobile) particles form clusters that grow upon cooling; however, a comparison of the UL mode clusters to the mobile and immobile particle clusters indicates that they are distinct entities. We also show that the cluster size provides an alternate method to distinguish localized and delocalized modes, offering a significant practical advantage over the finite-size scaling approach. We examine the trajectories of particles contributing most to the UL modes and find that they have a slightly enhanced mobility compared to the average, and we determine a characteristic time quantifying the persistence time of this excess mobility. This time scale is proportional to the structural relaxation time τ_α of the fluid, consistent with a prediction by Zwanzig [Phys. Rev. 156, 190 (1967)] for the lifetime of collective excitations in cooled liquids. Evidently, these collective excitations serve to facilitate relaxation but do not actually participate in the motion associated with barrier crossing events governing activated transport. They also serve as a possible concrete realization of the "facilitation" clusters postulated in previous modeling of glass-forming liquids.

Orientational correlations in fluids with quenched disorder

Snapshots of colloidal particles moving on disordered two-dimensional substrates can be used to extract equal-time many-body correlations in their positions. To understand the systematics of these correlations, we perform Monte Carlo simulations of a two-dimensional model fluid placed in a quenched disordered background. We use configurations generated from these simulations to compute translational and orientational two-point correlations at equal time, concentrating on correlations in local orientational order as a function of density and disorder strength. We calculate both the disorder averaged version of conventional two-point correlation functions for orientational order, as well as the disorder averaged version of a novel correlation function of time-averaged disorder-induced inhomogeneities in local orientation analogous to the Edwards-Anderson correlation function in spin systems. We demonstrate that these correlations can exhibit interesting nonmonotonic behavior in proximity to the underlying fluid-solid transition and suggest that this prediction should be experimentally accessible.

Static and dynamic correlation lengths in supercooled polymers

A key point to understand the glass transition is the relationship between structural and dynamic behavior experienced by a glass former when it approaches T_g. In this work, the relaxation in a simple bead-spring polymer system in the supercooled regime near its glass transition temperature was investigated with molecular dynamic simulations. We develop a new manner to look at the dynamic length scales in a supercooled polymeric system, focusing on correlated motion of particles in an isoconfigurational ensemble (that is, associated with the structure), as measured by Pearson's correlation coefficient. We found that while the usual dynamic four-point correlation length deviates from the structural (mosaic or point-to-set) length scale at low temperatures, Pearson's length behaves similarly to the static length in the whole temperature range. The results lead to a consensus of similar scaling of structural and dynamical length scales, reinforcing the idea of the theories of Adam-Gibbs and random first order transition.

Emergence and percolation of rigid domains during the colloidal glass transition

Using video microscopy, we measure local spatial constraints in disordered binary colloidal samples, ranging from dilute fluids to jammed glasses, and probe their spatial and temporal correlations to local dynamics during the glass transition. We observe the emergence of significant correlations between constraints and local dynamics within the Lindemann criterion, which coincides with the onset of glassy dynamics in supercooled liquids. Rigid domains in fluids are identified based on local constraints and demonstrate a percolation transition near the glass transition, accompanied by the emergence of dynamical heterogeneities. Our results show that spatial constraint instead of the geometry of amorphous structures is the key that connects the complex spatial-temporal correlations in disordered materials.

A universal state and its relaxation mechanisms of long-range interacting polygons

Using polygonal magnetic particles, we conduct experiments to explore the space-filling properties of anisotropic blocks with long-range interactions. In contrast to previous studies, we obtain the surprising finding that our systems' structures do not depend on the shape of building blocks: a single state, the hexagonal plastic crystal, appears as a universal attractor for a wide range of different polygons. This robust particle-shape independency appears as the interactions go beyond nearest neighbors. Particle shape plays an essential role in system relaxation, and determines the basic relaxation dynamics through a microscopic control parameter, internal roughness, produced by particle vertices. Thus our study reveals a new pattern-forming paradigm, in which particle shape plays little role in the static structure but determines the essential relaxation dynamics. Due to the ubiquity of long-range interactions and anisotropic building blocks, our discovery may shed new light on diverse problems involving structure formation, self-assembly, and packing.

Identity ordering and metastable clusters in fluids with random interactions

We use Langevin dynamics simulations to study dense two-dimensional systems of particles where all binary interactions are different in the sense that each interaction parameter is characterized by a randomly chosen number. We compare two systems that differ by the probability distributions from which the interaction parameters are drawn: uniform (U) and exponential (E). Both systems undergo neighborhood identity ordering and form metastable clusters in the fluid phase near the liquid-solid transition, but the effects are much stronger in E than in U systems. Possible implications of our results for the control of the structure of multicomponent alloys are discussed.

Dynamic Adsorption of Ions into Like-Charged Nanospace: A Dynamic Density Functional Theory Study

The adsorption processes of ions into charged nanospace are associated with many practical applications. Whereas a large number of microporous materials have been prepared toward efficient adsorption of ions from solutions, theoretical models that allow for capturing the characteristics of ion dynamic adsorption into like-charged nanopores are still few. The difficulty originates from the overlapping of electric potentials inside the pores. Herein, a theoretical model is proposed by incorporating dynamic density functional theory with modified Poisson equation for investigating the dynamic adsorption of ions into like-charged nanoslits. This model is rationalized by comparing the theoretical predictions with corresponding simulation results. Afterward, by analyzing the adsorption dynamics, we show that the overlapping effect is associated with the pore size, ion bulk concentration, and surface charge density, and it plays a dominant role in the coupling between the total adsorption amount of ions and total adsorption time. Specifically, with weak overlapping effect, the total adsorption amount is intuitively proportional to the total adsorption time; however, when the overlapping effect is strong, the total adsorption amount may be inversely proportional to the total adsorption time, indicating that both high adsorption amount and short adsorption time can be achieved simultaneously. This work provides a meaningful insight toward the rational design and optimization of microporous materials for efficient ion adsorption.

Evidence of a two-dimensional glass transition in graphene: Insights from molecular simulations

Liquids exhibit a sudden increase in viscosity when cooled fast enough, avoiding thermodynamically predicted route of crystallization. This phenomenon, known as glass transition, leads to the formation of non-periodic structures known as glasses. Extensive studies have been conducted on model materials to understand glass transition in two dimensions. However, despite the synthesis of disordered/amorphous single-atom thick structures of carbon, little attention has been given to glass transition in realistic two-dimensional materials such as graphene. Herein, using molecular dynamics simulation, we demonstrate the existence of glass transition in graphene leading to a realistic two-dimensional glassy structure, namely glassy graphene. We show that the resulting glassy structure exhibits excellent agreement with experimentally realized disordered graphene. Interestingly, this glassy graphene exhibits a wrinkled but stable structure, with reduced thermal vibration in comparison to its crystalline counterpart. We suggest that the topological disorder induced by glass transition governs the unique properties of this structure.

Experimental Test of the Cooperative Free Volume Rate Model under 1D Confinement: The Interplay of Free Volume, Temperature, and Polymer Film Thickness in Driving Segmental Mobility

We show that the Cooperative Free Volume (CFV) rate model, successful at modeling pressure-dependent dynamics, can be employed to describe the temperature and thickness dependence of the segmental time of polymers confined in thin films (1D confinement). The CFV model is based on an activation free energy that increases with the number of cooperating segments, which is determined by the system's free volume. Here, we apply the CFV model to new experimental results on the segmental relaxation of 1D confined poly(4-chlorostyrene), P4ClS, and find remarkable agreement over the whole temperature and thickness ranges investigated. This work further validates the robustness of the CFV model, which relates the effects of confinement on dynamics to pressure changes in the bulk, and supports the idea that confinement effects originate from local perturbations in density.

Vibrational Modes and Dynamic Heterogeneity in a Near-Equilibrium 2D Glass of Colloidal Kites

Using video microscopy and particle-tracking techniques developed for dense Brownian systems of polygons, we study the structure-dynamics relationship in a near-equilibrium 2D glass consisting of anisotropic Penrose kite-shaped colloids. Detailed vibrational properties of kite glasses, both translational and rotational, are obtained using covariance matrix techniques. Different from other colloidal glasses of spheres and ellipsoids, the vibrational modes of kite glasses at low frequencies show a strong translational character with spatially localized rotational modes and extended translational modes. Low-frequency quasilocalized soft modes commonly found in sphere glasses are absent in the translational phonon modes of kite glasses. Soft modes are observed predominantly in the rotational vibrations and correlate well with the spatial distribution of Debye-Waller factors. The local structural entropy field shows a strong correlation with the observed dynamic heterogeneity.

Dynamic properties of different liquid states in systems with competing interactions studied with lysozyme solutions

Recent studies of colloidal systems with a short-range attraction and long-range repulsion (SALR) have been demonstrated to have a generalized phase diagram with multiple liquid states defined by their structures. In this paper, we identify the different liquid states of previous experimentally studied lysozyme samples within this proposed generalized state diagram and explore the dynamic properties of each liquid state. We show that most lysozyme samples studied here and previously at low and intermediate concentrations are dispersed fluids while a few high concentration samples are randomly percolated liquids. In the dispersed fluid region, the short-time diffusion coefficient measured by neutron spin echo agrees well with the long time diffusion coefficient estimated with the solution viscosity. This dynamic feature is maintained even for some samples in the random percolated region. However, the short-time and long-time diffusion coefficients of random percolated fluids deviate at larger concentration and attraction strength. At high enough concentrations, the mean square displacement can be as slow as those of many glassy colloidal systems at time scales near the characteristic diffusion time even though these lysozyme samples remain in liquid states at the long-time limit. We thus identify the region in the generalized phase diagram where these equilibrium states with extremely slow local dynamics exist relative to bulk percolation and kinetic arrest (gel and glassy) transitions.

Surface Fluctuations Dominate the Slow Glassy Dynamics of Polymer-Grafted Colloid Assemblies

Colloids grafted with a corona layer of polymers show glassy behavior that covers a wide range of fragilities, with this behavior being tunable through variations in grafting density and grafting chain length. We find that the corona roughness, which is maximized for long chain lengths and sparse grafting, is directly correlated to the concentration-dependence of the system relaxation time (fragility). Relatively rougher colloids result in stronger liquids because their rotational motions become orientationally correlated across the whole system even at low particle loadings leading to an essentially Arrhenius-like concentration-dependence of the relaxation times near the glass transition. The smoother colloids do not show as much orientational correlation except at higher densities leading to fragile behavior. We therefore propose that these materials are an ideal model to study the physical properties of the glass transition.

Direct determination of structural heterogeneity in metallic glasses using four-dimensional scanning transmission electron microscopy

We report the first direct quantification of the structural heterogeneity in metallic glasses using intensity variance and angular correlation analyses of the 4-dimensional (4-D) scanning transmission electron microscopy (STEM) data. We demonstrate that the real-space reconstruction and analyses of the 4-D nanodiffraction data acquired using a pixelated fast STEM detector enables quantitative determination of the details of local structural heterogeneity, including the type, size, volume fraction and spatial distribution of local ordering at the nano- to meso-scale, beyond the limits of the previous measurements using conventional detectors. We show that different types of local ordering are present in Zr₅₅Co₂₅Al₂₀ glass, leading to a high degree of structural heterogeneity, with the total volume of locally ordered regions making up to ∼14% of the entire volume. These findings are significant, as the structure-property relationship in metallic glasses and other amorphous materials has been difficult to establish because of the lack of detailed structural information from experiments.

The race to the bottom: approaching the ideal glass?

Key to resolving the scientific challenge of the glass transition is to understand the origin of the massive increase in viscosity of liquids cooled below their melting temperature (avoiding crystallisation). A number of competing and often mutually exclusive theoretical approaches have been advanced to describe this phenomenon. Some posit a bona fide thermodynamic phase to an 'ideal glass', an amorphous state with exceptionally low entropy. Other approaches are built around the concept of the glass transition as a primarily dynamic phenomenon. These fundamentally different interpretations give equally good descriptions of the data available, so it is hard to determine which-if any-is correct. Recently however this situation has begun to change. A consensus has emerged that one powerful means to resolve this longstanding question is to approach the putative thermodynamic transition sufficiently closely, and a number of techniques have emerged to meet this challenge. Here we review the results of some of these new techniques and discuss the implications for the existence-or otherwise-of the thermodynamic transition to an ideal glass.

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.

Long-range stress correlations in viscoelastic and glass-forming fluids

A simple and rigorous approach to obtain stress correlations in viscoelastic liquids (including supercooled liquid and equilibrium amorphous systems) is proposed. The long-range dynamical correlations of local shear stress are calculated and analyzed in 2-dimensional space. It is established how the long-range character of the stress correlations gradually emerges as the relevant dynamical correlation length l grows in time. The correlation range l is defined by momentum propagation due to acoustic waves and vorticity diffusion which are the basic mechanisms for transmission of shear stress perturbations. We obtain the general expression defining the time- and distance-dependent stress correlation tensor in terms of material functions (generalized relaxation moduli). The effect of liquid compressibility is quantitatively analyzed; it is shown to be important at large distances and/or short times. The revealed long-range stress correlation effect is shown to be dynamical in nature and unconnected with static structural correlations in liquids (correlation length ξ_s). Our approach is based on the assumption that ξ_s is small enough as reflected in weak wave-number dependencies of the generalized relaxation moduli. We provide a simple physical picture connecting the elucidated long-range fluctuation effect with anisotropic correlations of the (transient) inherent stress field, and discuss its implications.

Glass Transition in Supercooled Liquids with Medium-Range Crystalline Order

The origin of the rapid dynamical slowdown in glass forming liquids in the growth of static length scales, possibly associated with identifiable structural ordering, is a much debated issue. Growth of medium range crystalline order (MRCO) has been observed in various model systems to be associated with glassy behavior. Such observations raise the question of whether molecular mechanisms for the glass transition in liquids with and without MRCO are the same. In this study we perform extensive molecular dynamics simulations of a number of glass forming liquids and show that the static and dynamic properties of glasses with MRCO are different from those of other glass forming liquids with no predominant local order. We also resolve an important issue regarding the so-called point-to-set method for determining static length scales, and demonstrate it to be a robust method for determining static correlation lengths in glass formers.

Local structure in deeply supercooled liquids exhibits growing lengthscales and dynamical correlations

Glasses are among the most widely used of everyday materials, yet the process by which a liquid's viscosity increases by 14 decades to become a glass remains unclear, as often contradictory theories provide equally good descriptions of the available data. Knowledge of emergent lengthscales and higher-order structure could help resolve this, but this requires time-resolved measurements of dense particle coordinates-previously only obtained over a limited time interval. Here we present an experimental study of a model colloidal system over a dynamic window significantly larger than previous measurements, revealing structural ordering more strongly linked to dynamics than previously found. Furthermore we find that immobile regions and domains of local structure grow concurrently with density, and that these regions have low configurational entropy. We thus show that local structure plays an important role at deep supercooling, consistent with a thermodynamic interpretation of the glass transition rather than a principally dynamic description.

Electric Field Assembly of Colloidal Superstructures

The assembly of materials from building blocks has been in the core of a wide range of applications from catalysis to photonics and electronics. External electric fields enable the interactions between building blocks to be controlled via induced dipoles. Dipolar interactions were used so far to obtain one-dimensional chains or three-dimensional non-close-packed lattices. However, complex colloidal assemblies and clusters of simple spherical particles are rare. Here we demonstrate a novel self-assembly approach enabling the formation of regular axially symmetric clusters, an array of colloidal assemblies as per design of posts, and hierarchical complex assemblies by using posts and dipolar interactions or combining them. Regulating the polarization of the particles from positive to negative allows us to control the interparticle interactions from attractive to repulsive at the poles or equator of the particles. Therefore, such particle-particle interactions enable switching between Saturn ring-like and candle-flame-like axially symmetric assemblies, which could potentially be exploited for display applications.

Nanoscale Structural Evolution and Anomalous Mechanical Response of Nanoglasses by Cryogenic Thermal Cycling

One of the central themes in the amorphous materials research is to understand the nanoscale structural responses to mechanical and thermal agitations, the decoding of which is expected to provide new insights into the complex amorphous structural-property relationship. For common metallic glasses, their inherent atomic structural inhomogeneities can be rejuvenated and amplified by cryogenic thermal cycling, thus can be decoded from their responses to mechanical and thermal agitations. Here, we reported an anomalous mechanical response of a new kind of metallic glass (nanoglass) with nanoscale interface structures to cryogenic thermal cycling. As compared to those metallic glasses by liquid quenching, the Sc₇₅Fe₂₅ (at. %) nanoglass exhibits a decrease in the Young's modulus but a significant increase in the yield strength after cryogenic cycling treatments. The abnormal mechanical property change can be attributed to the complex atomic rearrangements at the short- and medium- range orders due to the intrinsic nonuniformity of the nanoglass architecture. The present work gives a new route for designing high-performance metallic glassy materials by manipulating their atomic structures and helps for understanding the complex atomic structure-property relationship in amorphous materials.

Weak temperature dependence of ageing of structural properties in atomistic model glassformers

Ageing phenomena are investigated from a structural perspective in two binary Lennard-Jones glassformers, the Kob-Andersen and Wahnström mixtures. In both, the geometric motif assumed by the glassformer upon supercooling, the locally favoured structure (LFS), has been established. The Kob-Andersen mixture forms bicapped square antiprisms; the Wahnström model forms icosahedra. Upon ageing, we find that the structural relaxation time has a time-dependence consistent with a power law. However, the LFS population and potential energy increase and decrease, respectively, in a logarithmic fashion. Remarkably, over the time scales investigated, which correspond to a factor of 10⁴ change in relaxation times, the rate at which these quantities age appears almost independent of temperature. Only at temperatures far below the Vogel-Fulcher-Tamman temperature do the ageing dynamics slow.

Configuration correlation governs slow dynamics of supercooled metallic liquids

The origin of dramatic slowing down of dynamics in metallic glass-forming liquids toward their glass transition temperatures is a fundamental but unresolved issue. Through extensive molecular dynamics simulations, here we show that, contrary to the previous beliefs, it is not local geometrical orderings extracted from instantaneous configurations but the intrinsic correlation between configurations that captures the structural origin governing slow dynamics. More significantly, it is demonstrated by scaling analyses that it is the correlation length extracted from configuration correlation rather than dynamic correlation lengths that is the key to determine the drastic slowdown of supercooled metallic liquids. The key role of the configuration correlation established here sheds important light on the structural origin of the mysterious glass transition and provides an essential piece of the puzzle for the development of a universal theoretical understanding of glass transition in glasses.

Yielding of a model glass former: An interpretation with an effective system of icosahedra

We consider the yielding under simple shear of a binary Lennard-Jones glass former whose super-Arrhenius dynamics are correlated with the formation of icosahedral structures. We recast this glass former as an effective system of icosahedra [Pinney et al., J. Chem. Phys. 143, 244507 (2015)JCPSA60021-960610.1063/1.4938424]. Looking at the small-strain region of sheared simulations, we observe that shear rates affect the shear localization behavior particularly at temperatures below the glass transition as defined with a fit to the Vogel-Fulcher-Tamman equation. At higher temperature, shear localization starts immediately on shearing for all shear rates. At lower temperatures, faster shear rates can result in a delayed start in shear localization, which begins close to the yield stress. Building from a previous work which considered steady-state shear [Pinney et al., J. Chem. Phys. 143, 244507 (2015)JCPSA60021-960610.1063/1.4938424], we interpret the response to shear and the shear localization in terms of a local effective temperature with our system of icosahedra. We find that the effective temperatures of the regions undergoing shear localization increase significantly with increasing strain (before reaching a steady-state plateau).

Micro-structure and motion of two-dimensional dense short spherocylinder liquids

We numerically investigate the micro-structure and motion of 2D liquids composed of dense short spherocylinders, by reducing the shape aspect ratio from 3. It is found that reducing shape aspect ratio from 3 causes a smooth transition from heterogeneous structures composed of crystalline ordered domains with good tetratic alignment order to those with good hexagonal bond-orientational order at an aspect ratio equaling 1.35. In the intermediate regime, both structural orders are strongly deteriorated, and the translational hopping rate reaches a maximum due to the poor particle interlocking of the disordered structure. Shortening rod length allows easier rotation, induces monotonic increase of rotational hopping rates, and resumes the separation of rotational and translational hopping time scales at the small aspect ratio end, after the crossover of their rates in the intermediate regime. At the large shape aspect ratio end, the poor local tetratic order has the same positive effects on facilitating local rotational and translational hopping. In contrast, at the small shape aspect ratio end, the poor local bond orientational order has the opposite effects on facilitating local rotational and translational hopping.

Revealing the Link between Structural Relaxation and Dynamic Heterogeneity in Glass-Forming Liquids

Despite the use of glasses for thousands of years, the nature of the glass transition is still mysterious. On approaching the glass transition, the growth of dynamic heterogeneity has long been thought to play a key role in explaining the abrupt slowdown of structural relaxation. However, it still remains elusive whether there is an underlying link between structural relaxation and dynamic heterogeneity. Here, we unravel the link by introducing a characteristic time scale hiding behind an identical dynamic heterogeneity for various model glass-forming liquids. We find that the time scale corresponds to the kinetic fragility of liquids. Moreover, it leads to scaling collapse of both the structural relaxation time and dynamic heterogeneity for all liquids studied, together with a characteristic temperature associated with the same dynamic heterogeneity. Our findings imply that studying the glass transition from the viewpoint of dynamic heterogeneity is more informative than expected.

Growing length scale accompanying vitrification: A perspective based on nonsingular density fluctuations

A model for describing growing length scale accompanying the vitrification is introduced, in which we assume that in a subsystem whose density is above a certain threshold value, ρ_{c}, due to steric constraints, particle rearrangements are highly suppressed for a sufficiently long time period (∼structural relaxation time). We regard such a subsystem as a glassy cluster. With this assumption and without introducing any complicated thermodynamic arguments, we predict that with compression (increasing average density ρ) at a fixed temperature T in supercooled states, the characteristic length of the clusters, ξ, diverges as ξ∼(ρ_{c}-ρ)^{-2/d}, where d is the spatial dimensionality. This ξ measures the average persistence length of the steric constraints in blocking the rearrangement motions and is determined by the subsystem density. Additionally, with decreasing T at a fixed ρ, the length scale diverges in the same manner as ξ∼(T-T_{c})^{-2/d}, for which ρ is identical to ρ_{c} at T=T_{c}. The exponent describing the diverging length scale is the same as the one predicted by some theoretical models and indeed has been observed in some simulations and experiments. However, the basic mechanism for this divergence is different; that is, we do not invoke thermodynamic anomalies associated with the thermodynamic phase transition as the origin of the growing length scale. We further present arguements for the cooperative properties of the structural relaxation based on the clusters.

Molecular dynamic heterogeneity in relation to free volume and relaxation dynamics in organic glass-formers: oligomeric cis-1,4-poly(isoprene)

Herein, a combined study of the molecular rotation dynamics and free volume in cis-1,4-poly(isoprene) using two external probing techniques via ESR and PALS together with relaxation dynamics of the host medium via BDS is presented. The spectral evolution of the spin probe TEMPO from simulations over a wide range from 100 K up to 300 K exhibits three different regions of its correlation time consisting of a slow regime at low temperatures followed by the molecular dynamic heterogeneity zone from T = T = 155 K = 0.82 × T up to T_c ≅ 236 K = 1.26 × T and ending with a fast regime at high temperatures with the further characteristic ESR temperatures, T = 186 K ≅ T and T = 260 K. These are in close coincidence with four characteristic PALS temperatures: T = 160 K, T = 190 K, T = 227 K, and T = 263 K. Finally, using BDS, we revealed that the high-frequency features of the structural relaxation of 1,4-PIP 0.8k were related to the observed effects in the ESR and PALS response of the liquid state.

Linking slow dynamics and microscopic connectivity in dense suspensions of charged colloids

The quest to unravel the nature of the glass transition, where the viscosity of a liquid increases by many orders of magnitude, while its static structure remains largely unaffected, remains unresolved. While various structural and dynamical precursors to vitrification have been identified, a predictive and quantitative description of how subtle changes at the microscopic scale give rise to the steep growth in macroscopic viscosity is missing. It was recently proposed that the presence of long-lived bonded structures within the liquid may provide the long-sought connection between local structure and global dynamics. Here we directly observe and quantify the connectivity dynamics in liquids of charged colloids en route to vitrification using three-dimensional confocal microscopy. We determine the dynamic structure from the real-space van Hove correlation function and from the particle trajectories, providing upper and lower bounds on connectivity dynamics. Based on these data, we extend Dyre's model for the glass transition to account for particle-level structural dynamics; this results in a microscopic expression for the slowing down of relaxations in the liquid that is in quantitative agreement with our experiments. These results indicate how vitrification may be understood as a dynamical connectivity transition with features that are strongly reminiscent of rigidity percolation scenarios.

The role of local-geometrical-orders on the growth of dynamic-length-scales in glass-forming liquids

The precise nature of complex structural relaxation as well as an explanation for the precipitous growth of relaxation time in cooling glass-forming liquids are essential to the understanding of vitrification of liquids. The dramatic increase of relaxation time is believed to be caused by the growth of one or more correlation lengths, which has received much attention recently. Here, we report a direct link between the growth of a specific local-geometrical-order and an increase of dynamic-length-scale as the atomic dynamics in metallic glass-forming liquids slow down. Although several types of local geometrical-orders are present in these metallic liquids, the growth of icosahedral ordering is found to be directly related to the increase of the dynamic-length-scale. This finding suggests an intriguing scenario that the transient icosahedral connectivity could be the origin of the dynamic-length-scale in metallic glass-forming liquids.

Measurements of growing surface tension of amorphous-amorphous interfaces on approaching the colloidal glass transition

There is mounting evidence indicating that relaxation dynamics in liquids approaching their glass transition not only become increasingly cooperative, but the relaxing regions also become more compact in shape. Of the many theories of the glass transition, only the random first-order theory—a thermodynamic framework—anticipates the surface tension of relaxing regions to play a role in deciding both their size and morphology. However, owing to the amorphous nature of the relaxing regions, even the identification of their interfaces has not been possible in experiments hitherto. Here, we devise a method to directly quantify the dynamics of amorphous–amorphous interfaces in bulk supercooled colloidal liquids. Our procedure also helped unveil a non-monotonic evolution in dynamical correlations with supercooling in bulk liquids. We measure the surface tension of the interfaces and show that it increases rapidly across the mode-coupling area fraction. Our experiments support a thermodynamic origin of the glass transition.

The existence of interfaces, separating distinct relaxing regions, has been predicted in glass  theory, but a direct proof remains challenging due to the amorphous nature of glasses. Ganapathi et al. identify and measure the surface tension of these interfaces in bulk supercooled colloidal liquids.

Structural predictor for nonlinear sheared dynamics in simple glass-forming liquids

Glass-forming liquids subjected to sufficiently strong shear universally exhibit striking nonlinear behavior; for example, a power-law decrease of the viscosity with increasing shear rate. This phenomenon has attracted considerable attention over the years from both fundamental and applicational viewpoints. However, the out-of-equilibrium and nonlinear nature of sheared fluids have made theoretical understanding of this phenomenon very challenging and thus slower to progress. We find here that the structural relaxation time as a function of the two-body excess entropy, calculated for the extensional axis of the shear flow, collapses onto the corresponding equilibrium curve for a wide range of pair potentials ranging from harsh repulsive to soft and finite. This two-body excess entropy collapse provides a powerful approach to predicting the dynamics of nonequilibrium liquids from their equilibrium counterparts. Furthermore, the two-body excess entropy scaling suggests that sheared dynamics is controlled purely by the liquid structure captured in the form of the two-body excess entropy along the extensional direction, shedding light on the perplexing mechanism behind shear thinning.