Probing a liquid to glass transition in equilibrium

Does the Naı̈ve Mode-Coupling Power Law Divergence Provide an Objective Determination of the Crossover Temperature in Glass Formation Behavior?

The glass formation temperature range is commonly divided into a weakly supercooled regime at higher temperatures and a deeply supercooled regime at lower temperatures, with a change in the physical mechanisms that govern dynamics often postulated to occur at the crossover between these regimes. This crossover temperature Tc is widely determined based on a fit of relaxation time vs temperature data to a power law divergence form predicted by the naı̈ve mode coupling theory (MCT). Here, we show, based on simulation data spanning polymeric, small molecule organic, metallic, and inorganic glass formers, that this approach does not yield an objective measure of a crossover temperature. Instead, the value of Tc is determined by the lowest temperature Tmin employed in the fit, and no regime of stationary or convergent Tc value is generally observed as Tmin is varied. Nor does the coefficient of determination R2 provide any robust means of selecting a fit range and thus a value of Tc. These results may require a re-evaluation of published results that have employed the fit MCT Tc value as a metric of temperature-dependent dynamics or a benchmark for depth of supercooling, and they highlight a need for the field to converge on a more objective determination of any posited crossover temperature between high and low temperature regimes of glass formation.

Collective dynamic length increases monotonically in pinned and unpinned glass forming systems

The Random First-Order Transition (RFOT) theory predicts that transport proceeds by the cooperative movement of particles in domains, whose sizes increase as a liquid is compressed above a characteristic volume fraction, ϕd. The rounded dynamical transition around ϕd, which signals a crossover to activated transport, is accompanied by a growing correlation length that is predicted to diverge at the thermodynamic glass transition density (>ϕd). Simulations and imaging experiments probed the single particle dynamics of mobile particles in response to pinning all the particles in a semi-infinite space or randomly pinning (RP) a fraction of particles in a liquid at equilibrium. The extracted dynamic length increases non-monotonically with a peak around ϕd, which not only depends on the pinning method but is also different from ϕd of the actual liquid. This finding is at variance with the results obtained using the small wavelength limit of a four-point structure factor for unpinned systems. To obtain a consistent picture of the growth of the dynamic length, one that is impervious to the use of RP, we introduce a multiparticle structure factor, Smpc(q,t), that probes collective dynamics. The collective dynamical length, calculated from the small wave vector limit of Smpc(q,t), increases monotonically as a function of the volume fraction in a glass-forming binary mixture of charged colloidal particles in both unpinned and pinned systems. This prediction, which also holds in the presence of added monovalent salt, may be validated using imaging experiments.

Effects of quenched disorder on the kinetics and pathways of phase transition in a soft colloidal system

Although impurities are unavoidable in real-world and experimental systems, most numerical studies on nucleation focus on pure (impurity-free) systems. As a result, the role of impurities in phase transitions remains poorly understood, especially for systems with complex free energy landscapes featuring one or more intermediate metastable phases. In this study, we employed Monte Carlo simulations to investigate the effects of static impurities (quenched disorder) of varying length scales and surface morphologies on the crystal nucleation mechanism and kinetics in the Gaussian core model system-a representative model for soft colloidal systems. We first explored how the nucleation free energy barrier and critical cluster size are influenced by the fraction of randomly pinned (or, static) particles (fp) and the size (np) of the pinned region or cluster. Both the nucleation free energy barrier and critical cluster size increase sharply with increasing fp but decrease as np grows for a given fraction of pinned particles, eventually approaching the homogeneous nucleation limit. On examining the impact of impurity's surface morphology on nucleation kinetics, we observed that the nucleation barrier significantly decreases with increasing the impurity (or, seed) size with crystalline surface morphologies with body-centered cubic showing the greatest facilitation. Interestingly, seeds with random surface roughness had little effect on nucleation kinetics. In addition, the polymorphic identity of particles in the final crystalline phase is influenced by both the seed's surface morphology and system size. This study further provides crucial insights into the intricate relationship between surface-induced local structural fluctuations and the selection of the polymorphic identity in the final crystalline phase, which is essential for understanding and controlling crystallization processes in experiments.

Superionic UO_{2} crystal: How to model its relaxation and diffusion via a microscopic theory of glass-forming liquids

UO_{2} is a crucial nuclear material but its behaviors are elusive due to the impact of superionic diffusion. Herein, we introduce a simple but effective theoretical model to describe the intrinsic superionicity of UO_{2} at the quantitative level. Our idea stems from a close similarity between superionic crystals and supercooled liquids. Namely, we view UO_{2} as a randomly pinned hard-sphere fluid in the framework of the elastically collective nonlinear Langevin equation theory. This treatment allows us to fully evaluate the contribution of local, collective, pinning, and screening effects to the molecular dynamics of UO_{2} without complex computational processes. Finite-temperature factors are considered via volumetric expansion during isobaric heating. On that basis, we satisfactorily explain recent large-scale atomistic simulations on UO_{2} under various thermodynamic conditions. Our calculations also reveal that UO_{2} is equivalent to an intermediate glass former. Its structural relaxation, self-diffusion, and shear deformation are strongly correlated near the onset of superionicity. These intimate correlations are reminiscent of the famed Dyre shoving model in the soft-matter community. Our results promise to facilitate the development of diverse energy applications of UO_{2}.

Effect of the presence of pinned particles on the structural parameters of a liquid and correlation between structure and dynamics at the local level

Pinning particles at the equilibrium configuration of the liquid is expected not to affect the structure and any property that depends on the structure while slowing down the dynamics. This leads to a breakdown of the structure dynamics correlation. Here, we calculate two structural quantities: the pair excess entropy, S2, and the mean field caging potential, the inverse of which is our structural order parameter (SOP). We show that when the pinned particles are treated the same way as the mobile particles, both S2 and SOP of the mobile particles remain the same as those of the unpinned system, and the structure dynamics correlation decreases with an increase in pinning density, "c." However, when we treat the pinned particles as a different species, even if we consider that the structure does not change, the expression of S2 and SOP changes. The microscopic expressions show that the interaction between a pinned particle and a mobile particle affects S2 and SOP more than the interaction between two mobile particles. We show that a similar effect is also present in the calculation of the excess entropy and is the primary reason for the well-known vanishing of the configurational entropy at high temperatures. We further show that, contrary to the common belief, the pinning process does change the structure. When these two effects are considered, both S2 and SOP decrease with an increase in "c," and the correlation between the structural parameters and the dynamics continues even for higher values of "c."

Anisotropic-Isotropic Transition of Cages at the Glass Transition

Characterizing the local structural evolution is an essential step in understanding the nature of glass transition. In this work, we probe the evolution of Voronoi cell geometry in simple glass models by simulations and colloid experiments, and find that the individual particle cages deform anisotropically in supercooled liquid and isotropically in glass. We introduce an anisotropy parameter k for each Voronoi cell, whose mean value exhibits a sharp change at the mode-coupling glass transition ϕ_{c}. Moreover, a power law of packing fraction ϕ∝q_{1}^{d} is discovered in the supercooled liquid regime with d>D, in contrast to d=D in the glass regime, where q_{1} is the first peak position of structure factor, and D is the space dimension. This power law is qualitatively explained by the change of k. The active motions in supercooled liquid are spatially correlated with long axes rather than short axes of Voronoi cells. In addition, the dynamic slowing down approaching the glass transition can be well characterized through a modified free-volume model based on k. These findings reveal that the structural parameter k is effective in identifying the structure-dynamics correlations and the glass transition in these systems.

Glassy phases of the Gaussian core model

We present results from molecular dynamics simulations exploring the supercooled dynamics of the Gaussian Core Model in the low- and intermediate-density regimes. In particular, we analyse the transition from the low-density hard-sphere-like glassy dynamics to the high-density one. The dynamics at low densities is well described by the caging mechanism, giving rise to intermittent dynamics. At high densities, the particles undergo a more continuous motion in which the concept of cage loses its meaning. We elaborate on the idea that these different supercooled dynamics are in fact the precursors of two different glass states.

Soft pinning: Experimental validation of static correlations in supercooled molecular glass-forming liquids

Enormous enhancement in the viscosity of a liquid near its glass transition is a hallmark of glass transition. Within a class of theoretical frameworks, it is connected to growing many-body static correlations near the transition, often called "amorphous ordering." At the same time, some theories do not invoke the existence of such a static length scale in the problem. Thus, proving the existence and possible estimation of the static length scales of amorphous order in different glass-forming liquids is very important to validate or falsify the predictions of these theories and unravel the true physics of glass formation. Experiments on molecular glass-forming liquids become pivotal in this scenario as the viscosity grows several folds (∼ 10 14 ), and simulations or colloidal glass experiments fail to access these required long-time scales. Here we design an experiment to extract the static length scales in molecular liquids using dilute amounts of another large molecule as a pinning site. Results from dielectric relaxation experiments on supercooled Glycerol with different pinning concentrations of Sorbitol and Glucose, as well as the simulations on a few model glass-forming liquids with pinning sites, indicate the versatility of the proposed method, opening possible new avenues to study the physics of glass transition in other molecular liquids.

Non-phononic density of states of two-dimensional glasses revealed by random pinning

The vibrational density of states of glasses is considerably different from that of crystals. In particular, there exist spatially localized vibrational modes in glasses. The density of states of these non-phononic modes has been observed to follow g(ω) ∝ ω4, where ω is the frequency. However, in two-dimensional systems, the abundance of phonons makes it difficult to accurately determine this non-phononic density of states because they are strongly coupled to non-phononic modes and yield strong system-size and preparation-protocol dependencies. In this article, we utilize the random pinning method to suppress phonons and disentangle their coupling with non-phononic modes and successfully calculate their density of states as g(ω) ∝ ω4. We also study their localization properties and confirm that low-frequency non-phononic modes in pinned systems are truly localized without far-field contributions. We finally discuss the excess density of states over the Debye value that results from the hybridization of phonons and non-phononic modes.

World beyond the nearest neighbors

The structure beyond the nearest neighbor atoms in liquid and glass is characterized by the medium-range order (MRO). In the conventional approach, the MRO is considered to result directly from the short-range order (SRO) in the nearest neighbors. To this bottom-up approach starting with the SRO, we propose to add a top-down approach in which global collective forces drive liquid to form density waves. The two approaches are in conflict with each other, and the compromise produces the structure with the MRO. The driving force to produce density waves provides the stability and stiffness to the MRO, and controls various mechanical properties. This dual framework provides a novel perspective for description of the structure and dynamics of liquid and glass.

Creating bulk ultrastable glasses by random particle bonding

A recent breakthrough in glass science has been the synthesis of ultrastable glasses via physical vapor deposition techniques. These samples display enhanced thermodynamic, kinetic and mechanical stability, with important implications for fundamental science and technological applications. However, the vapor deposition technique is limited to atomic, polymer and organic glass-formers and is only able to produce thin film samples. Here, we propose a novel approach to generate ultrastable glassy configurations in the bulk, via random particle bonding, and using computer simulations we show that this method does indeed allow for the production of ultrastable glasses. Our technique is in principle applicable to any molecular or soft matter system, such as colloidal particles with tunable bonding interactions, thus opening the way to the design of a large class of ultrastable glasses.

Low-frequency vibrational states in ideal glasses with random pinning

Glasses exhibit spatially localized vibrations in the low-frequency regime. These localized modes emerge below the boson peak frequency ω_{BP}, and their vibrational densities of state follow g(ω)∝ω^{4} (ω is frequency). Here, we attempt to address how the localized vibrations behave through the ideal glass transition. To do this, we employ a random pinning method, which enables us to study the thermodynamic glass transition. We find that the localized vibrations survive even in equilibrium glass states. Remarkably, the localized vibrations still maintain the properties of appearance below ω_{BP} and g(ω)∝ω^{4}. Our results provide important insight into the material properties of ideal glasses.

Evidence for growing structural correlation length in colloidal supercooled liquids

Using video microscopy, we measure the long-time diffusion coefficients of colloidal particles at different concentrations. The measured diffusion coefficients start to deviate from theoretical predictions based on random collision models upon entering the supercooled regime. The theoretical diffusion relation is recovered by assigning an effective mass proportional to the size of structurally correlated clusters to the diffusing particles, providing an indirect method to probe the growth of static correlation length scales approaching the glass transition. This method is tested and validated in the crystallization of mono-disperse colloids in quasi-two-dimensional experiments. The correlation length obtained for a binary colloidal liquid increases by a power law toward a critical packing fraction of ∼0.79. The system relaxation time exhibits a power-law dependence on the correlation length in agreement with dynamical facilitation theories.

Screening and collective effects in randomly pinned fluids: a new theoretical framework

We propose a theoretical framework for the dynamics of bulk isotropic hard-sphere systems in the presence of randomly pinned particles and apply this theory to supercooled water to validate it. Structural relaxation is mainly governed by local and non-local activated process. As the pinned fraction grows, a local caging constraint becomes stronger and the long range collective aspect of relaxation is screened by immobile obstacles. Different responses of the local and cooperative motions results in subtle predictions for how the alpha relaxation time varies with pinning and density. Our theoretical analysis for the relaxation time of water with pinned molecules quantitatively well describe previous simulations. In addition, the thermal dependence of relaxation for unpinned bulk water is also consistent with prior computational and experimental data.

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.

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.

Revealing the characteristic length of random close packing via critical-like random pinning

By randomly pinning particles in fluidized states and finding the local energy minima, we form static packings of mono-disperse disks that resemble random close packing, when only n_c = 2.6% of the particles are pinned. The packings are isostatic and exhibit typical critical scalings of the jamming transition. The non-triviality of n_c is manifested mainly in two aspects. First, n_c acts as a critical point, leading to bifurcated critical scalings in its vicinity. The criticality of n_c is also demonstrated in the packings of weakly polydisperse disks. Second, n_c sets a length scale in agreement with the characteristic length of random close packing. With robust evidence, we show that this agreement is generally true for both mono- and poly-disperse particles and in both two and three dimensions. The randomness inherited from fluidized states by random pinning thus interprets the randomness of random close packing from a unique perspective.

Lattice Glass Model in Three Spatial Dimensions

The understanding of thermodynamic glass transition has been hindered by the lack of proper models beyond mean-field theories. Here, we propose a three-dimensional lattice glass model on a simple cubic lattice that exhibits the typical dynamics observed in fragile supercooled liquids such as two-step relaxation, super-Arrhenius growth in the relaxation time, and dynamical heterogeneity. Using advanced Monte Carlo methods, we compute the thermodynamic properties deep inside the glassy temperature regime, well below the onset temperature of the slow dynamics. The specific heat has a finite jump towards the thermodynamic limit with critical exponents close to those expected from the hyperscaling and the random first-order transition theory for the glass transition. We also study an effective free energy of glasses, the Franz-Parisi potential, as a function of the overlap between equilibrium and quenched configurations. The effective free energy indicates the existence of a first-order phase transition, consistent with the random first-order transition theory. These findings strongly suggest that the glassy dynamics of the model has its origin in thermodynamics.

Structured randomness: jamming of soft discs and pins

Simulations are used to find the zero temperature jamming threshold, φj, for soft, bidisperse disks in the presence of small fixed particles, or "pins", arranged in a lattice. The presence of pins leads, as one expects, to a decrease in φj. Structural properties of the system near the jamming threshold are calculated as a function of the pin density. While the correlation length exponent remains ν = 1/2 at low pin densities, the system is mechanically stable with more bonds, yet fewer contacts than the Maxwell criterion implies in the absence of pins. In addition, as pin density increases, novel bond orientational order and long-range spatial order appear, which are correlated with the square symmetry of the pin lattice.

Glass forming liquids in a quenched random potential

The response of a model two-dimensional colloidal glass former to an externally imposed spatially random potential, which acts as a quenched disorder, is investigated using numerical simulations, motivated by recent experiments and also mean field predictions. The external potential induces the onset of the glassy dynamics at increasingly smaller field roughness, with increasing packing fraction of the particulate assembly, and the existence of aging processes within the glassy regime is also observed. Furthermore, along the axis of increasing field roughness, the dynamical slowdown is not correlated to the hexatic order within the supercooled regime.

Particle pinning suppresses spinodal criticality in the shear-banding instability

Strained amorphous solids often fail mechanically by creating a shear band. It had been understood that the shear-banding instability is usefully described as crossing a spinodal point (with disorder) in an appropriate thermodynamic description. It remained contested, however, whether the spinodal is critical (with divergent correlation length) or not. Here we offer evidence for critical spinodal by using particle pinning. For a finite concentration of pinned particles the correlation length is bounded by the average distance between pinned particles, but without pinning it is bounded by the system size.

Effect of Pinning on the Yielding Transition of Amorphous Solids

Using numerical simulations, we have studied the yielding response, in the athermal quasistatic limit, of a model amorphous material having inclusions in the form of randomly pinned particles. We show that, with increasing pinning concentration, the plastic activity becomes more spatially localized, resulting in smaller stress drops, and a corresponding increase in the magnitude of strain where yielding occurs. We demonstrate that, unlike the spatially heterogeneous and avalanche led yielding in the case of the unpinned glass, for the case of large pinning concentration, yielding takes place via a spatially homogeneous proliferation of localized events.

Relation between Heterogeneous Frozen Regions in Supercooled Liquids and Non-Debye Spectrum in the Corresponding Glasses

Recent numerical studies on glassy systems provide evidence for a population of non-Goldstone modes (NGMs) in the low-frequency spectrum of the vibrational density of states D(ω). Similarly to Goldstone modes (GMs), i.e., phonons in solids, NGMs are soft low-energy excitations. However, differently from GMs, NGMs are localized excitations. Here we first show that the parental temperature T^{*} modifies the GM/NGM ratio in D(ω). In particular, the phonon attenuation is reflected in a parental temperature dependency of the exponent s(T^{*}) in the low-frequency power law D(ω)∼ω^{s(T^{*})}, with 2≤s(T^{*})≤4. Second, by comparing s(T^{*}) with s(p), i.e., the same quantity obtained by pinning a p particle fraction, we suggest that s(T^{*}) reflects the presence of dynamical heterogeneous regions of size ξ^{3}∝p. Finally, we provide an estimate of ξ as a function of T^{*}, finding a mild power law divergence, ξ∼(T^{*}-T_{d})^{-α/3}, with T_{d} the dynamical crossover temperature and α falling in the range α∈[0.8,1.0].

Vapor-Deposited Ethylbenzene Glasses Approach "Ideal Glass" Density

Spectroscopic ellipsometry was used to characterize vapor-deposited glasses of ethylbenzene (T_g = 115.7 K). For this system, previous calorimetric experiments have established that a transition to the ideal glass state is expected to occur near 101 K (the Kauzmann temperature, T_K) if the low-temperature supercooled liquid has the properties expected based upon extrapolation from above T_g. Ethylbenzene glasses were vapor-deposited at substrate temperatures between 100 (∼0.86 T_g) and 116 K (∼T_g), using deposition rates of 0.02-2.1 nm/s. Down to 103 K, glasses prepared in the limit of low deposition rate have densities consistent with the extrapolated supercooled liquid. The highest density glass is within 0.15% of the density expected for the ideal glass. These results support the hypothesis that the extrapolated properties of supercooled ethylbenzene are correct to within just a few Kelvin of T_K, consistent with the existence of a phase transition to an ideal glass state at T_K.

Configurational entropy of glass-forming liquids

The configurational entropy is one of the most important thermodynamic quantities characterizing supercooled liquids approaching the glass transition. Despite decades of experimental, theoretical, and computational investigation, a widely accepted definition of the configurational entropy is missing, its quantitative characterization remains fraught with difficulties, misconceptions, and paradoxes, and its physical relevance is vividly debated. Motivated by recent computational progress, we offer a pedagogical perspective on the configurational entropy in glass-forming liquids. We first explain why the configurational entropy has become a key quantity to describe glassy materials, from early empirical observations to modern theoretical treatments. We explain why practical measurements necessarily require approximations that make its physical interpretation delicate. We then demonstrate that computer simulations have become an invaluable tool to obtain precise, nonambiguous, and experimentally relevant measurements of the configurational entropy. We describe a panel of available computational tools, offering for each method a critical discussion. This perspective should be useful to both experimentalists and theoreticians interested in glassy materials and complex systems.

Zero-temperature glass transition in two dimensions

Liquids cooled towards the glass transition temperature transform into amorphous solids that have a wide range of applications. While the nature of this transformation is understood rigorously in the mean-field limit of infinite spatial dimensions, the problem remains wide open in physical dimensions. Nontrivial finite-dimensional fluctuations are hard to control analytically, and experiments fail to provide conclusive evidence regarding the nature of the glass transition. Here, we develop Monte Carlo methods for two-dimensional glass-forming liquids that allow us to access equilibrium states at sufficiently low temperatures to directly probe the glass transition in a regime inaccessible to experiments. We find that the liquid state terminates at a thermodynamic glass transition which occurs at zero temperature and is associated with an entropy crisis and a diverging static correlation length. Our results thus demonstrate that a thermodynamic glass transition can occur in finite dimensional glass-formers.

Ideal Glass States Are Not Purely Vibrational: Insight from Randomly Pinned Glasses

We use computer simulations to probe the thermodynamic and dynamic properties of a glass former that undergoes an ideal glass transition because of the presence of randomly pinned particles. We find that even deep in the equilibrium glass state, the system relaxes to some extent because of the presence of localized excitations that allow the system to access different inherent structures, thus giving rise to a nontrivial contribution to the entropy. By calculating with high accuracy the vibrational part of the entropy, we show that also in the equilibrium glass state thermodynamics and dynamics give a coherent picture, and that glasses should not be seen as a disordered solid in which the particles undergo just vibrational motion but instead as a system with a highly nonlinear internal dynamics.

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.

A geometrical criterion for glass transition in soft-sphere fluids

As glass-forming fluids become colder and denser, structural rearrangements become slow and eventually cease. For hard-sphere fluids, percolation of particles unable to change neighbors (T1-inactive particles) signals the glass transition. To investigate this geometrical criterion for mobility in soft-sphere systems, we simulate monodisperse fluids interacting with a generalized Weeks-Chandler-Andersen (WCA) potential in metastable equilibrium, using our previously developed crystal-avoiding method. We find that the vanishing diffusivity as the glass transition is approached can be described by a power law below the onset temperature of super-Arrhenius behavior. By mapping the soft spheres to hard spheres based on mean collision energy, we find that the diffusivity versus effective volume fraction curves collapse onto the hard-sphere curve for all systems studied. We find that the onset of super-Arrhenius behavior and the MCT dynamic glass transition correlate well with temperature when the T1-inactive particles form clusters of two particles on average and when the T1-inactive clusters percolate the entire system, respectively. Our findings provide new insight into the structural origin of glassy dynamics.

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.

Probing the non-Debye low-frequency excitations in glasses through random pinning

We investigate the properties of the low-frequency spectrum in the density of states [Formula: see text] of a 3D model glass former. To magnify the non-Debye sector of the spectrum, we introduce a random pinning field that freezes a finite particle fraction to break the translational invariance and shifts all of the vibrational frequencies of the extended modes toward higher frequencies. We show that non-Debye soft localized modes progressively emerge as the fraction p of pinned particles increases. Moreover, the low-frequency tail of [Formula: see text] goes to zero as a power law [Formula: see text], with [Formula: see text] and [Formula: see text] above a threshold fraction [Formula: see text].

Atomic clusters with addressable complexity

A general formulation for constructing addressable atomic clusters is introduced, based on one or more reference structures. By modifying the well depths in a given interatomic potential in favour of nearest-neighbour interactions that are defined in the reference(s), the potential energy landscape can be biased to make a particular permutational isomer the global minimum. The magnitude of the bias changes the resulting potential energy landscape systematically, providing a framework to produce clusters that should self-organise efficiently into the target structure. These features are illustrated for small systems, where all the relevant local minima and transition states can be identified, and for the low-energy regions of the landscape for larger clusters. For a 55-particle cluster, it is possible to design a target structure from a transition state of the original potential and to retain this structure in a doubly addressable landscape. Disconnectivity graphs based on local minima that have no direct connections to a lower minimum provide a helpful way to visualise the larger databases. These minima correspond to the termini of monotonic sequences, which always proceed downhill in terms of potential energy, and we identify them as a class of biminimum. Multiple copies of the target cluster are treated by adding a repulsive term between particles with the same address to maintain distinguishable targets upon aggregation. By tuning the magnitude of this term, it is possible to create assemblies of the target cluster corresponding to a variety of structures, including rings and chains.

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.

Structural-dynamical transition in the Wahnström mixture

In trajectory space, dynamical heterogeneities in glass-forming liquids correspond to the emergence of a dynamical phase transition between an active phase poor in local structure and an inactive phase which is rich in local structure. We support this scenario with the study of a model additive mixture of Lennard-Jones particles, quantifying how the choice of the relevant structural and dynamical observable affects the transition in trajectory space. We find that the low mobility, structure-rich phase is dominated by icosahedral order. Applying a non-equilibrium rheological protocol, we connect local order to the emergence of mechanical rigidity.

Pinning susceptibility: a novel method to study growth of amorphous order in glass-forming liquids

The existence and growth of amorphous order in supercooled liquids approaching glass transition is a subject of intense research. Even after decades of work, there is still no clear consensus on the molecular mechanisms that lead to a rapid slowing down of liquid dynamics approaching this putative transition. The existence of a correlation length associated with amorphous order has recently been postulated and has also been estimated using multi-point correlation functions which cannot be calculated easily in experiments. Thus the study of growing amorphous order remains mostly restricted to systems like colloidal glasses and simulations of model glass-forming liquids. In this work, we propose an experimentally realizable yet simple susceptibility to study the growth of amorphous order. We then demonstrate the validity of this approach for a few well-studied model supercooled liquids and obtain results which are consistent with other conventional methods.

Study on the Glass Transition Process of Polymer System Using Differential Scanning Calorimetry and Fourier Transform Infrared Spectroscopy

The change in the infrared spectrum of polymer samples with temperature and their differential scanning calorimetry (DSC) experimental results are analyzed. According to the van't Hoff equation at constant pressure, the changes in the absorbance ratio corresponding to high and low vibrational states are calculated, and the apparent enthalpy differences of the vibration energy states transformation of the characteristic group can be obtained. From the experimental results, we can find that characteristic vibration modes of a chemical group in a polymer are under the influence of the glass transition process of the polymer with a different extent. The characteristic vibration modes of the same chemical group behave differently due to the influence of the polymer system at which the chemical moiety is situated.

Experimental Evidence for a Structural-Dynamical Transition in Trajectory Space

Among the key insights into the glass transition has been the identification of a nonequilibrium phase transition in trajectory space which reveals phase coexistence between the normal supercooled liquid (active phase) and a glassy state (inactive phase). Here, we present evidence that such a transition occurs in experiments. In colloidal hard spheres, we find a non-Gaussian distribution of trajectories leaning towards those rich in locally favored structures (LFSs), associated with the emergence of slow dynamics. This we interpret as evidence for a nonequilibrium transition to an inactive LFS-rich phase. Reweighting trajectories reveals a first-order phase transition in trajectory space between a normal liquid and a LFS-rich phase. We also find evidence for a purely dynamical transition in trajectory space.

Nonmonotonic dynamic correlations in quasi-two-dimensional confined glass-forming liquids

It has been broadly accepted that the behavior of glass-forming liquids, where their relaxation dynamics exhibit a pronounced slowdown as they are cooled toward the glass transition temperature, is caused by the increase in one or more correlation lengths. However, the role of length scales in the dynamics of glass-forming liquids is not clearly established, and past simulation work that suggests a surprising nonmonotonic temperature evolution of spatial dynamical correlations near the mode-coupling crossover temperature has been both questioned and supported by subsequent work. Here, using molecular dynamics simulation, we also show a striking maximum in the dynamic length scale ξ_{c}^{dyn} at a given temperature, but the temperature of this maximum is found to shift as the size of the confined system increases. Furthermore, we find that such a maximum disappears for all geometry sizes considered when a rough wall is replaced with a smooth, hard wall, suggesting that the nature of the nonmonotonic temperature dependence of ξ_{c}^{dyn} does not reflect an intrinsic property of bulk liquids, but originates from wall effects. Our results provide new insights into the dynamics of glass-forming liquids, particularly for quasi-two-dimensional systems.

Structure and dynamics in network-forming materials

The study of the structure and dynamics of network-forming materials is reviewed. Experimental techniques used to extract key structural information are briefly considered. Strategies for building simulation models, based on both targeting key (experimentally-accessible) materials and on systematically controlling key model parameters, are discussed. As an example of the first class of materials, a key target system, SiO₂, is used to highlight how the changing structure with applied pressure can be effectively modelled (in three dimensions) and used to link to both experimental results and simple structural models. As an example of the second class the topology of networks of tetrahedra in the MX₂ stoichiometry are controlled using a single model parameter linked to the M-X-M bond angles. The evolution of ordering on multiple length-scales is observed as are the links between the static structure and key dynamical properties. The isomorphous relationship between the structures of amorphous Si and SiO₂ is discussed as are the similarities and differences in the phase diagrams, the latter linked to potential polyamorphic and 'anomalous' (e.g. density maxima) behaviour. Links to both two-dimensional structures for C, Si and Ge and near-two-dimensional bilayers of SiO₂ are discussed. Emerging low-dimensional structures in low temperature molten carbonates are also uncovered.

The role of fivefold symmetry in suppressing crystallization

Although long assumed to have an important role in the suppression of crystallization and the development of glassformers, the effect of local fivefold symmetry has never been directly tested. Here we consider whether such suppression of crystallization has a kinetic or thermodynamic nature and investigate its mechanism. We introduce a model in which the degree of fivefold symmetry can be tuned by favouring arrangements of particles in pentagonal bipyramids. We thus show that fivefold symmetry has both kinetic and thermodynamic effects on the mechanism of crystallization to a face-centred cubic crystal. Our results suggest that the mechanism of crystallization suppression is related to the surface tension between fluid and crystal. Interestingly, the degree of fivefold symmetry has little effect on crystal growth rate, suggesting that growth may be only weakly coupled to fluid structure in hard sphere like systems. Upon increasing the fivefold symmetry, we find a first-order transition to an alternative icosahedra-rich phase. At intermediate bias strengths we find a one-component glassformer.

The suppression of crystallization due to the appearance of structures with fivefold symmetry is widely adopted, but its kinetic and thermodynamic origin remains elusive. Taffs et al. show that fivefold symmetry substantially slows down the nucleation rate but not the crystal growth rate as expected.

Point-to-set lengths, local structure, and glassiness

The growing sluggishness of glass-forming liquids is thought to be accompanied by growing structural order. The nature of such order, however, remains hotly debated. A decade ago, point-to-set (PTS) correlation lengths were proposed as measures of amorphous order in glass formers, but recent results raise doubts as to their generality. Here, we extend the definition of PTS correlations to agnostically capture any type of growing order in liquids, be it local or amorphous. This advance enables the formulation of a clear distinction between slowing down due to conventional critical ordering and that due to glassiness, and provides a unified framework to assess the relative importance of specific local order and generic amorphous order in glass formation.

Phase Transition for Quenched Coupled Replicas in a Plaquette Spin Model of Glasses

We study a three-dimensional plaquette spin model whose low temperature dynamics is glassy, due to localized defects and effective kinetic constraints. The thermodynamics of this system is smooth at all temperatures. We show that coupling it to a second system with a fixed (quenched) configuration leads to a phase transition, at finite coupling. The order parameter is the overlap between the copies, and the transition is between phases of low and high overlap. We find critical points whose properties are consistent with random-field Ising universality. We analyze the interfacial free energy cost between the high- and low-overlap states that coexist at (and below) the critical point, and we use this cost as the basis for a finite-size scaling analysis. We discuss these results in the context of mean-field and dynamical facilitation theories of the glass transition.