Collective Relaxation Dynamics in a Three-Dimensional Lattice Glass Model

We numerically elucidate the microscopic mechanisms controlling the relaxation dynamics of a three-dimensional lattice glass model that has static properties compatible with the approach to a random first-order transition. At low temperatures, the relaxation is triggered by a small population of particles with low-energy barriers forming mobile clusters. These emerging quasiparticles act as facilitating defects responsible for the spatially heterogeneous dynamics of the system, whose characteristic length scales remain strongly coupled to thermodynamic fluctuations. We compare our findings both with existing theoretical models and atomistic simulations of glass formers.

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.

Method to probe the pronounced growth of correlation lengths in active glass-forming liquids using an elongated probe

The growth of correlation lengths in equilibrium glass-forming liquids near the glass transition is considered a critical finding in the quest to understand the physics of glass formation. These understandings helped us understand various dynamical phenomena observed in supercooled liquids. It is known that at least two different length scales exist; one is of thermodynamic origin, while the other is dynamical in nature. Recent observations of glassy dynamics in biological and synthetic systems where the external or internal driving source controls the dynamics, apart from the usual thermal noise, lead to the emergence of the field of active glassy matter. A question of whether the physics of glass formation in these active systems is also accompanied by growing dynamic and static lengths is indeed timely. In this article, we probe the growth of dynamic and static lengths in a model active glass system using rod-like elongated probe particles, an experimentally viable method. We show that the dynamic and static lengths in these nonequilibrium systems grow much more rapidly than their passive counterparts. We then offer an understanding of the violation of the Stokes-Einstein relation and Stokes-Einstein-Debye relation using these lengths via a scaling theory.

Benzothiadiazole-based rotation and possible antipolar order in carboxylate-based metal-organic frameworks

By modifying organic ligands of metal-organic framework with dipolar units, they turn suitable for various applications, e.g., in the field of sensor systems or switching of gas permeation. Dipolar linkers in the organic ligand are capable to rotate in certain temperature and frequency ranges. The copper-bearing paddlewheel shaped metal-organic frameworks ZJNU-40 and JLU-Liu30 possess such a polarizable dipole moment due to their benzothiadiazole moiety in the organic ligands. Here, we investigate the molecular rotor behavior of benzothiadiazole units of the two carboxylate-based MOFs by dielectric spectroscopy and computational simulation. Our dielectric results provide clear evidence for significant reorientational relaxation dynamics of these rotors, revealing various characteristics of glasslike freezing upon cooling. The calculated rotational energy barriers are consistent with experimentally determined barriers for single-dipole dynamics. Moreover, for JLU-Liu30 we find hints at antipolar ordering below about 300 K.

Coherent X-ray Scattering Reveals Nanoscale Fluctuations in Hydrated Proteins

Hydrated proteins undergo a transition in the deeply supercooled regime, which is attributed to rapid changes in hydration water and protein structural dynamics. Here, we investigate the nanoscale stress-relaxation in hydrated lysozyme proteins stimulated and probed by X-ray Photon Correlation Spectroscopy (XPCS). This approach allows us to access the nanoscale dynamics in the deeply supercooled regime (T = 180 K), which is typically not accessible through equilibrium methods. The observed stimulated dynamic response is attributed to collective stress-relaxation as the system transitions from a jammed granular state to an elastically driven regime. The relaxation time constants exhibit Arrhenius temperature dependence upon cooling with a minimum in the Kohlrausch-Williams-Watts exponent at T = 227 K. The observed minimum is attributed to an increase in dynamical heterogeneity, which coincides with enhanced fluctuations observed in the two-time correlation functions and a maximum in the dynamic susceptibility quantified by the normalized variance χ_T. The amplification of fluctuations is consistent with previous studies of hydrated proteins, which indicate the key role of density and enthalpy fluctuations in hydration water. Our study provides new insights into X-ray stimulated stress-relaxation and the underlying mechanisms behind spatiotemporal fluctuations in biological granular materials.

Nonlinear response theory for Markov processes. IV. The asymmetric double-well potential model revisited

The dielectric response of noninteracting dipoles is discussed in the framework of the classical model of stochastic reorientations in an asymmetric double-well potential (ADWP). In the nonlinear regime, this model exhibits some pecularities in the static response. We find that the saturation behavior of the symmetric double-well potential model does not follow the Langevin function and only in the linear regime are the standard results recovered. If a finite asymmetry is assumed, then the nonlinear susceptibilities are found to change the sign at a number of characteristic temperatures that depend on the magnitude of the asymmetry, as has been observed earlier for the third-order and fifth-order responses. If the kinetics of the barrier crossing in the ADWP model is described as a two-state model, then we can give analytical expressions for the values of the characteristic temperatures. The results for the response obtained from a (numerical) solution of the Fokker-Planck equation for the Brownian motion in a model ADWP behaves very similarly to the two-state model for high barriers. For small barriers no clear-cut timescale separation between the barrier crossing process and the intrawell relaxation exists and the model exhibits a number of timescales. In this case, the frequency-dependent linear susceptibility at low temperatures is dominated by the fast intrawell transitions and at higher temperatures by the barrier crossing kinetics. We find that for nonlinear susceptibilities the latter process appears to be more important and the intrawell transitions play only a role at the lowest temperatures.

Dynamical susceptibilities near ideal glass transitions

Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalized mode-coupling theory of supercooled liquids to inhomogeneous environments. This provides a first-principles-based, systematic, and rigorous way of deriving high-point dynamical susceptibilities from variations of the many-body dynamic structure factors with respect to their conjugate field. This framework allows for a fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition. The behavior of these dynamical susceptibilities is then studied in the context of simplified self-consistent relaxation models.

Glass transition as a topological phase transition

The glass transition is described as a phase transition in the system of topologically protected excitations in matter structure. The critical behavior of the system is considered in both static and dynamic cases. It is shown that the proposed model reproduces most of the characteristic thermodynamic and kinetic properties of glass transition: the Vogel-Fulcher-Tammann law, the behavior of susceptibility and nonlinear susceptibilities, and heat capacity behavior as well as the appearance of a boson peak in the frequency dependence of the dynamic structure factor near the glass transition temperature.

Molecular Cross-correlations Govern Structural Rearrangements in a Nonassociating Polar Glass Former

Self- and cross-correlation dynamics of deeply supercooled liquids were recently identified using photon correlation spectroscopy on the one hand and dielectric investigations on the other. These results fueled a controversial discussion whether the "generic" response identified by photon correlation spectroscopy, or rather the nonuniversal dielectric response, reflect the liquid's structural relaxation. The present study employs physical aging and oscillatory shear rheology to directly access the structural relaxation of a nonassociating glass-forming liquid and reveals that collective equilibrium fluctuations of simple liquids and not single-particle dynamics govern their structural relaxation. The present results thus challenge recent views that the glassy response of polar supercooled liquids can generically be decomposed into a Debye-type, supramolecular response and a single-particle dynamics with the latter reflecting the "true" structural relaxation. Furthermore, the current findings underscore the pivotal role dielectric spectroscopy plays in glass science as one of the rare molecular-level reorientation techniques that senses dynamical cooperativity directly.

Strong increase of correlations in liquid glycerol observed by nonlinear dielectric techniques

Nonlinear dielectric measurements are an important tool to access material properties and dynamics concealed in their linear counterparts, but the available data are often intermittent and, on occasion, even contradictory. Employing and refining a recently developed technique for high ac field dielectric measurements in the static limit, we ascertain nonlinear effects in glycerol over a wide temperature range from 230 to 320 K. We find that the temperature dependence of the Piekara factor a, which quantifies the saturation effect, changes drastically around 290 K, from ∂a/∂T = +1.4 to -130 in units of 10^-18 V² m^-2 K^-1. These high values of |a| quantify not only elevated dielectric saturation effects but also indicate a temperature driven increase in higher-order orientational correlations and considerable correction terms with respect to the central limit theorem. No signature of this feature can be found in the corresponding low field data.

Numerical test of the replica-symmetric Hamiltonian for correlations of the critical state of spin glasses in a field

A growing body of evidence indicates that the sluggish low-temperature dynamics of glass formers (e.g., supercooled liquids, colloids, or spin glasses) is due to a growing correlation length. Which is the effective field theory that describes these correlations? The natural field theory was drastically simplified by Bray and Roberts in 1980. More than 40 years later, we confirm the tenets of Bray and Roberts's theory by studying the Ising spin glass in an externally applied magnetic field, both in four spatial dimensions (data obtained from the Janus collaboration) and on the Bethe lattice.

Q-dependent collective relaxation dynamics of glass-forming liquid Ca0.4K0.6(NO3)1.4 investigated by wide-angle neutron spin-echo

The relaxation behavior of glass formers exhibits spatial heterogeneity and dramatically changes upon cooling towards the glass transition. However, the underlying mechanisms of the dynamics at different microscopic length scales are not fully understood. Employing the recently developed wide-angle neutron spin-echo spectroscopy technique, we measured the Q-dependent coherent intermediate scattering function of a prototypical ionic glass former Ca_0.4K_0.6(NO₃)_1.4, in the highly viscous liquid state. In contrast to the structure modulated dynamics for Q < 2.4 Å⁻¹, i.e., at and below the structure factor main peak, for Q > 2.4 Å⁻¹, beyond the first minimum above the structure factor main peak, the stretching exponent exhibits no temperature dependence and concomitantly the relaxation time shows smaller deviations from Arrhenius behavior. This finding indicates a change in the dominant relaxation mechanisms around a characteristic length of 2π/(2.4 Å⁻¹) ≈ 2.6 Å, below which the relaxation process exhibits a temperature independent distribution and more Arrhenius-like behavior.

Length scale dependence is important for understanding the collective relaxation dynamics in glass-forming liquids. Here, the authors find in liquid Ca_0.4K_0.6(NO₃)_1.4 a change in the dominant relaxation mechanisms around 2.6 Å, below which the relaxation process exhibits a temperature independent distribution and more Arrhenius-like behavior.

Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time

The noncrystalline glassy state of matter plays a role in virtually all fields of materials science and offers complementary properties to those of the crystalline counterpart. The caveat of the glassy state is that it is out of equilibrium and therefore exhibits physical aging, i.e., material properties change over time. For half a century, the physical aging of glasses has been known to be described well by the material-time concept, although the existence of a material time has never been directly validated. We do this here by successfully predicting the aging of the molecular glass 4-vinyl-1,3-dioxolan-2-one from its linear relaxation behavior. This establishes the defining property of the material time. Via the fluctuation-dissipation theorem, our results imply that physical aging can be predicted from thermal-equilibrium fluctuation data, which is confirmed by computer simulations of a binary liquid mixture.

Kinetic Yvon-Born-Green theory of the linear dielectric constant and complex permittivity of isotropic polar fluids

The theory of the linear static dielectric constant and linear complex permittivity of isotropic polar fluids is formulated starting from the coupled Langevin equations describing the rototranslational dynamics of long-range interacting molecules with thermal agitation and subjected to external forces and torques. To this aim, adequate reduced densities are introduced and equations governing their dynamics derived. In the equilibrium zero frequency limit, integral expressions for the Kirkwood correlation factor g_{K} are given, transparently showing that the popular method consisting in comparing g_{K} with 1 in order to deduce pair dipolar ordering has no serious theoretical grounding. In the dynamical situation, the complex permittivity spectrum of a simple liquid is shown to exhibit an infinite discrete set of relaxation times, some of which may have thermally activated behavior. The theory is also shown to contain all previous results derived in the area provided molecular inertial effects are ignored, so restricting the range of validity of the theory to frequencies much below the far-infrared region. Finally, the theory can be adapted without much effort to relaxation of interacting magnetic nanoparticles for which macroscopic magnetic anisotropy arising from the assembly of nanoparticles is neglected.

Medium-range atomic correlation in simple liquids. I. Distinction from short-range order

Physical properties of liquids and glasses are controlled not only by the short-range order (SRO) in the nearest-neighbor atoms but also by the medium-range order (MRO) observed for atoms beyond the nearest neighbors. In this article the nature of the MRO as the descriptor of point-to-set atomic correlation is discussed focusing on simple liquids, such as metallic liquids. Through the results of x-ray diffraction and simulation with classical potentials we show that the third peak of the pair-distribution function, which describes the MRO, shows a distinct change in temperature dependence at the glass transition, whereas the first peak, which represents the SRO, changes smoothly through the glass transition. The result suggests that the glass transition is induced by the freezing of the MRO rather than that of the SRO, implying a major role of the MRO on the viscosity of supercooled liquid.

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

Synchrotron x-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol, and D-sorbitol) from above the glass transition temperatures T_g to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently, reflecting the different natures of bonding. Upon cooling from 400 K, the O⋯O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrests at T_g; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above T_g, glycerol loses the fewest HBs per OH for a given temperature rise scaled by T_g, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C⋯C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above T_g causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol being the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

Medium-range atomic correlation in simple liquids. II. Theory of temperature dependence

The spatial atomic correlations in liquids and glasses extend often significantly beyond the nearest neighbors. Such correlations, called the medium-range order (MRO), affect many physical properties, but their nature is not well understood. In this article the variation of the MRO with temperature is calculated based upon the concept of the atomic-level pressure, focusing on simple liquids, such as metallic liquids. It is shown that the structural coherence length that characterizes MRO follows the Curie-Weiss law with a negative Curie temperature as observed by experiment and simulation. It is also shown that the glass transition is induced by freezing of the MRO, rather than the freezing of the nearest-neighbor shell. The implications of these results are discussed.

Data-driven criterion for the solid-liquid transition of two-dimensional self-propelled colloidal particles far from equilibrium

We establish an explicit data-driven criterion for identifying the solid-liquid transition of two-dimensional self-propelled colloidal particles in the far from equilibrium parameter regime, where the transition points predicted by different conventional empirical criteria for melting and freezing diverge. This is achieved by applying a hybrid machine learning approach that combines unsupervised learning with supervised learning to analyze a huge amount of the system's configurations in the nonequilibrium parameter regime on an equal footing. Furthermore, we establish a generic data-driven evaluation function, according to which the performance of different empirical criteria can be systematically evaluated and improved. In particular, by applying this evaluation function, we identify a new nonequilibrium threshold value for the long-time diffusion coefficient, based on which the predictions of the corresponding empirical criterion are greatly improved in the far from equilibrium parameter regime. These data-driven approaches provide a generic tool for investigating phase transitions in complex systems where conventional empirical ones face difficulties.

First- and third-order shear nonlinearities across the structural relaxation peak of the deeply supercooled pharmaceutical liquid indomethacin

Nonlinear rheological properties of viscous indomethacin are studied in the frequency range of its structural relaxation, that is, in a range so far inaccessible to standard techniques involving medium-amplitude oscillatory shear amplitudes. The first- and third-order nonlinearity parameters thus recorded using a sequence of small and large shear excitations in a time efficient manner are compared with predictions from rheological models. By properly phase cycling the shear amplitudes, build-up and decay transients are recorded. Analogous to electrical-field experiments, these transients yield direct access to the structural relaxation times under linear and nonlinear shearing conditions. To demonstrate the broader applicability of the present approach, transient analyses are also carried out for the glass formers glycerol, ortho-terphenyl, and acetaminophen.

Long-Range Static and Dynamic Previtreous Effects in Supercooled Squalene-Impact of Strong Electric Field

This article presents evidence for the long-range previtreous changes of two static properties: the dielectric constant (ε) and its strong electric field related counterpart, the nonlinear dielectric effect (NDE). Important evidence is provided for the functional characterizations of ε(T) temperature changes by the ‘Mossotti Catastrophe’ formula, as well as for the NDE vs. T evolution by the relations resembling those developed for critical liquids. The analysis of the dynamic properties, based on the activation energy index, excluded the Vogel–Fulcher–Tammann (VFT) relation as a validated tool for portraying the evolution of the primary relaxation time. This result questions the commonly applied ‘Stickel operator’ routine as the most reliable tool for determining the dynamic crossover temperature. In particular, the strong electric field radically affects the distribution of the relaxation times, the form of the evolution of the primary relaxation time, and the fragility. The results obtained in this paper support the concept of a possible semi-continuous phase transition hidden below T_g. The studies were carried out in supercooled squalene, a material with an extremely low electric conductivity, a strongly elongated molecule, and which is vitally important for biology and medicine related issues.

On Structural Rearrangements Near the Glass Transition Temperature in Amorphous Silica

The formation of clusters was analyzed in a topologically disordered network of bonds of amorphous silica (SiO₂) based on the Angell model of broken bonds termed configurons. It was shown that a fractal-dimensional configuron phase was formed in the amorphous silica above the glass transition temperature T_g. The glass transition was described in terms of the concepts of configuron percolation theory (CPT) using the Kantor-Webman theorem, which states that the rigidity threshold of an elastic percolating network is identical to the percolation threshold. The account of configuron phase formation above T_g showed that (i) the glass transition was similar in nature to the second-order phase transformations within the Ehrenfest classification and that (ii) although being reversible, it occurred differently when heating through the glass–liquid transition to that when cooling down in the liquid phase via vitrification. In contrast to typical second-order transformations, such as the formation of ferromagnetic or superconducting phases when the more ordered phase is located below the transition threshold, the configuron phase was located above it.

Mutual Information in Molecular and Macromolecular Systems

The relaxation properties of viscous liquids close to their glass transition (GT) have been widely characterised by the statistical tool of time correlation functions. However, the strong influence of ubiquitous non-linearities calls for new, alternative tools of analysis. In this respect, information theory-based observables and, more specifically, mutual information (MI) are gaining increasing interest. Here, we report on novel, deeper insight provided by MI-based analysis of molecular dynamics simulations of molecular and macromolecular glass-formers on two distinct aspects of transport and relaxation close to GT, namely dynamical heterogeneity (DH) and secondary Johari–Goldstein (JG) relaxation processes. In a model molecular liquid with significant DH, MI reveals two populations of particles organised in clusters having either filamentous or compact globular structures that exhibit different mobility and relaxation properties. In a model polymer melt, MI provides clearer evidence of JG secondary relaxation and sharper insight into its DH. It is found that both DH and MI between the orientation and the displacement of the bonds reach (local) maxima at the time scales of the primary and JG secondary relaxation. This suggests that, in (macro)molecular systems, the mechanistic explanation of both DH and relaxation must involve rotation/translation coupling.

Glass Transition in Rice Pasta as Observed by Combined Neutron Scattering and Time-Domain NMR

Experimental protocols aiming at the characterisation of glass transition often suffer from ambiguity. The ambition of the present study is to describe the glass transition in a complex, micro heterogeneous system, the dry rice pasta, in a most unambiguous manner, minimising the influence of technique-specific bias. To this end, we apply an unprecedented combination of experimental techniques. Apart from the usually used NMR and DSC, we employ, in a concurrent manner, neutron transmission, diffraction, and Compton scattering. This enables us to investigate the glass transition over a range of spatio-temporal scales that stretches over seven orders of magnitude. The results obtained by neutron diffraction and DSC reveal that dry rice pasta is almost entirely amorphous. Moreover, the glass transition is evidenced by neutron transmission and diffraction data and manifested as a significant decrease of the average sample number density in the temperature range between 40 and 60 °C. At the microscopic level, our NMR, neutron transmission and Compton scattering results provide evidence of changes in the secondary structure of the starch within the dry rice pasta accompanying the glass transition, whereby the long-range order provided by the polymer structure within the starch present in the dry rice pasta is partially lost.

Amorphous Order and Nonlinear Susceptibilities in Glassy Materials

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

Quantifying dielectric permittivities in the nonlinear regime

In contrast to the static dielectric permittivity,ε, associated with linear response, its high-field counterpart,ε_E, is not a material specific quantity, but rather depends on the experimental method used to determine the nonlinear dielectric effect (NDE). Here, we defineε_Ein a manner consistent with how high field permittivities are typically derived from a capacitance measurement using high voltages. Based upon characterizing the materials nonlinear behavior via its third order susceptibility,χ₃, the relations between a givenχ₃and the observableε_Eis calculated for six different experimental or theoretical approaches to NDEs in the static limit. It is argued that the quantityχ₃is superior overε_Eor the Piekara factor, (ε_E-ε)/E², because it facilitates an unambiguous comparison among different experimental techniques and it provides a more robust connection between experiment and theory.

Modeling non-linear dielectric susceptibilities of supercooled molecular liquids

Advances in high-precision dielectric spectroscopy have enabled access to non-linear susceptibilities of polar molecular liquids. The observed non-monotonic behavior has been claimed to provide strong support for theories of dynamic arrest based on the thermodynamic amorphous order. Here, we approach this question from the perspective of dynamic facilitation, an alternative view focusing on emergent kinetic constraints underlying the dynamic arrest of a liquid approaching its glass transition. We derive explicit expressions for the frequency-dependent higher-order dielectric susceptibilities exhibiting a non-monotonic shape, the height of which increases as temperature is lowered. We demonstrate excellent agreement with the experimental data for glycerol, challenging the idea that non-linear response functions reveal correlated relaxation in supercooled liquids.

Translational dynamics of a rod-like probe in supercooled liquids: an experimentally realizable method to study Stokes-Einstein breakdown, dynamic heterogeneity, and amorphous order

The use of probe molecules to extract the local dynamical and structural properties of complex dynamical systems is an age-old technique both in simulations and in experiments. A lot of important information which is not immediately accessible from bulk measurements can be accessed via these local measurements. Still, a detailed understanding of how a probe particle dynamics is affected by the surrounding liquid medium is lacking, especially in the supercooled temperature regime. This work shows how the translational dynamics of a rod-like particle immersed in a supercooled liquid can give us information on the growth of the correlation length scales associated with dynamical heterogeneity and the multi-body static correlations in the medium. This work also provides an understanding of the breakdown of Stokes-Einstein and Stokes-Einstein-Debye relations in supercooled liquids along with a unified scaling theory that rationalizes all the observed results. Finally, this work proposes a novel yet simple method accessible in experiments to measure the growth of these important length scales in molecular glass-forming liquids.

Understanding Slow and Heterogeneous Dynamics in Model Supercooled Glass-Forming Liquids

Glasses are ubiquitous in nature. Many common items such as ketchups, cosmetic products, toothpaste, etc. and metallic glasses are examples of such glassy materials whose dynamical and rheological properties matter in our daily life. The dynamics of these glass-forming systems are known to be very sluggish and heterogeneous, but a detailed understanding of the origin of such slowing down is still lacking. Slow heterogeneous dynamics occur in a wide variety of systems at scales ranging from microscopic to macroscopic. Polymeric liquids, granular material, such as powder and sand, gels, and foams and also metallic alloys show such complex glassy dynamics at appropriate conditions. Recently, the existence of dynamical heterogeneity has also been found in biological systems starting from collective cell migration in a monolayer of cells to embryonic morphogenesis, cancer invasion, and wound healing. Extensive research in the past decade or so lead to the understanding that there are growing dynamic and static correlation lengths associated with the observed dynamical heterogeneity and rapid rise in viscosity. In this review, we have highlighted the recent developments on measuring these correlation lengths in glass-forming liquids and their possible implications in the physics of the glass transition.

Searching for the Gardner Transition in Glassy Glycerol

We search for a Gardner transition in glassy glycerol, a standard molecular glass, measuring the third harmonics cubic susceptibility χ_{3}^{(3)} from slightly below the usual glass transition temperature down to 10 K. According to the mean-field picture, if local motion within the glass were becoming highly correlated due to the emergence of a Gardner phase then χ_{3}^{(3)}, which is analogous to the dynamical spin-glass susceptibility, should increase and diverge at the Gardner transition temperature T_{G}. We find instead that upon cooling |χ_{3}^{(3)}| decreases by several orders of magnitude and becomes roughly constant in the regime 100-10  K. We rationalize our findings by assuming that the low temperature physics is described by localized excitations weakly interacting via a spin-glass dipolar pairwise interaction in a random magnetic field. Our quantitative estimations show that the spin-glass interaction is twenty to fifty times smaller than the local random field contribution, thus rationalizing the absence of the spin-glass Gardner phase. This hints at the fact that a Gardner phase may be suppressed in standard molecular glasses, but it also suggests ways to favor its existence in other amorphous solids and by changing the preparation protocol.

Scaling Law Describes the Spin-Glass Response in Theory, Experiments, and Simulations

The correlation length ξ, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer.

Mechanical and dielectric response within and beyond the linear regime

In this work a comparison of dielectric and mechanical data is presented based on experiments within the linear response limit and beyond that limit. The linear dynamic and shear-mechanical response is discussed in terms of the molecular supercooled liquid tetramethyl-tetraphenyl-trisiloxane. As the dynamics measured by the two methods depict the same temperature-dependence, the underlying cause for the observed responses is assumed to be identical for both methods, namely structural relaxation. The comparison of dielectric and mechanical measurements under high excitation amplitudes reveals that this cannot be assumed for the nonlinear response: Mechanical experiments on metallic glasses suggest that involved energies are clearly beyond k _B T, with observed nonlinear effects based on the activation of microstructural plastic rearrangements. In contrast, nonlinear dielectric measurements on another molecular glass-former involve energies clearly below k _B T, so that nonlinear dielectric effects occur due to energy uptake from the electric field or entropy-based changes in the dynamics, but are very unlikely connected to the triggering of plastic rearrangements by the applied electric field.

Universal correlations between the fragility and interparticle repulsion of glass-forming liquids

A recently published analytical model describing and predicting elasticity, viscosity, and fragility of metallic melts is applied for the analysis of about 30 nonmetallic glassy systems, ranging from oxide network glasses to alcohols, low-molecular-weight liquids, polymers, plastic crystals, and even ionic glass formers. The model is based on the power-law exponent λ representing the steepness parameter of the repulsive part of the inter-atomic or inter-molecular potential and the thermal-expansion parameter α_T determined by the attractive anharmonic part of the effective interaction. It allows fitting the typical super-Arrhenius temperature variation of the viscosity or dielectric relaxation time for various classes of glass-forming matter, over many decades. We discuss the relation of the model parameters found for all these different glass-forming systems to the fragility parameter m and detect a correlation of λ and m for the non-metallic glass formers, in accord with the model predictions. Within the framework of this model, the fragility of glass formers can be traced back to microscopic model parameters characterizing the intermolecular interactions.

Dynamical phase transitions and their relation to structural and thermodynamic aspects of glass physics

We review recent developments in structural-dynamical phase transitions in trajectory space based on dynamic facilitation theory. An open question is how the dynamic facilitation perspective on the glass transition may be reconciled with thermodynamic theories that posit collective reorganization accompanied by a growing static length scale and, eventually, a vanishing configurational entropy. In contrast, dynamic facilitation theory invokes a dynamical phase transition between an active phase (close to the normal liquid) and an inactive phase, which is glassy and whose order parameter is either a time-averaged dynamic or structural quantity. In particular, the dynamical phase transition in systems with non-trivial thermodynamics manifests signatures of a lower critical point that lies between the mode-coupling crossover and the putative Kauzmann temperature, at which a thermodynamic phase transition to an ideal glass state would occur. We review these findings and discuss such criticality in the context of the low-temperature decrease in configurational entropy predicted by thermodynamic theories of the glass transition.

On Viscous Flow in Glass-Forming Organic Liquids

The two-exponential Sheffield equation of viscosity η(T) = A1·T·[1 + A2·exp(Hm/RT)]·[1 + C·exp(Hd/RT)], where A1, A2, Hm, C, and Hm are material-specific constants, is used to analyze the viscous flows of two glass-forming organic materials-salol and α-phenyl-o-cresol. It is demonstrated that the viscosity equation can be simplified to a four-parameter version: η(T) = A·T·exp(Hm/RT)]·[1 + C·exp(Hd/RT)]. The Sheffield model gives a correct description of viscosity, with two exact Arrhenius-type asymptotes below and above the glass transition temperature, whereas near the Tg it gives practically the same results as well-known and widely used viscosity equations. It is revealed that the constants of the Sheffield equation are not universal for all temperature ranges and may need to be updated for very high temperatures, where changes occur in melt properties leading to modifications of A and Hm for both salol and α-phenyl-o-cresol.

The structural α-relaxation times of prilocaine confined in 1 nm pores of molecular sieves: quantitative explanation by the coupling model

The molecular glass-former and pharmaceutical, prilocaine, distinguishes itself by exhibiting seven general and fundamental dynamic and thermodynamic properties [Z. Wojnarowska, et al., J. Phys. Chem. B, 2015, 39, 12699.], all of which have been explained using the coupling model. What has not been studied before are the changes in properties of the structural α-relaxation of prilocaine when subjected to extreme nano-confinement in spaces with a size of about 1 nm. Recently, Ruis et al. [G. N. Ruiz, et al., Phys. Chem. Chem. Phys., 2019, 21, 15576.] measured the α-relaxation times, τ_α,conf(T), of prilocaine confined in 1 nm pores of molecular sieves. They found that τ_α,conf(T) are significantly reduced from those of bulk prilocaine, τ_α,bulk(T), and assume a weaker temperature dependence. The data in toto pose a challenge for any theory of glass transition to explain quantitatively. The coupling model (CM) was applied to this problem to predict the α-relaxation times of prilocaine when cooperativity is removed, which is expected because only a few prilocaine molecules can fit into the 1 nm pores. The results from the CM are in quantitative agreement with the experimental values of τ_α,conf(T) and the temperature dependence. The success is nontrivial because no other extant theory can do the same to the best of our knowledge.

Fragile-to-strong crossover, growing length scales, and dynamic heterogeneity in Wigner glasses

Colloidal particles, which are ubiquitous, have become ideal testing grounds for the structural glass transition theories. In these systems glassy behavior arises as the density of the particles is increased. Thus, soft colloidal particles with varying degree of softness capture diverse glass-forming properties, observed normally in molecular glasses. Brownian dynamics simulations for a binary mixture of micron-sized charged colloidal suspensions show that tuning the softness of the interaction potential, achievable by changing the monovalent salt concentration results in a continuous transition from fragile to strong behavior. Remarkably, this is found in a system where the well characterized interaction potential between the colloidal particles is isotropic. We also show that the predictions of the random first-order transition (RFOT) theory quantitatively describes the universal features such as the growing correlation length, ξ∼(ϕ_{K}/ϕ-1)^{-ν} with ν=2/3 where ϕ_{K}, the analog of the Kauzmann temperature, depends on the salt concentration. As anticipated by the RFOT predictions, we establish a causal relationship between the growing correlation length and a steep increase in the relaxation time and dynamic heterogeneity as the system is compressed. The broad range of fragility observed in Wigner glasses is used to draw analogies with molecular and polymer glasses. The large variations in the fragility are normally found only when the temperature dependence of the viscosity is examined for a large class of diverse glass-forming materials. In sharp contrast, this is vividly illustrated in a single system that can be experimentally probed. Our work also shows that the RFOT predictions are accurate in describing the dynamics over the entire density range, regardless of the fragility of the glasses.

Predicting the α-relaxation time of glycerol confined in 1.16 nm pores of zeolitic imidazolate frameworks

Uhl et al. [J. Chem. Phys., 2019, 150, 024504] studied the molecular dynamics of glycerol confined in a microporous zeolitic imidazolate framework (ZIF-8) with well-defined pore diameters of 1.16 nm by broadband dielectric spectroscopy. Of interest is a fast process in the central part of the pores identified as the α-relaxation of the confined supercooled glycerol with relaxation times τα,conf(T) reduced from τα,bulk(T) of bulk glycerol and having a temperature dependence different from the super-Arrhenius temperature of the latter. The focus of Uhl et al. was relating the confined molecular dynamics to the cooperativity length scales Lcorr(T) of molecular motion above the glass transition, and deducing the limiting high-temperature value of the correlation length of about 1.22 nm. Not yet considered by anyone are the observed values of τα,conf(T) and temperature dependence. Since the cooperativity length scales Lcorr(T) were found to be larger than the pore size of ZIF-8 over the temperature range studied and the density of the glycerol in the pore is possibly lower than the bulk, the cooperativity of the α-relaxation of glycerol confined in ZIF-8 is drastically reduced. Thus, within the framework of the Coupling Model (CM), τα,conf(T) should be nearly the same as the primitive relaxation time τ0(T) for glycerol when devoid of intermolecular coupling and cooperativity. Consistent with the absence of cooperativity of the glycerol confined in ZIF-8, we find the calculated τα,conf(T) are either the same or slightly longer than the calculated values of τ0(T). The quantitative prediction of the CM is verified. At this time we know of no other theory that can make such a quantitative prediction.

Active glasses

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

Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids

Many systems, including polymers and molecular liquids, when adequately cooled and/or compressed, solidify into a disordered solid, i.e., a glass. The transition is not abrupt, featuring progressive decrease of the microscopic mobility and huge slowing down of the relaxation. A distinctive aspect of glass-forming materials is the microscopic dynamical heterogeneity (DH), i.e., the presence of regions with almost immobile particles coexisting with others where highly mobile ones are located. Following the first compelling evidence of a strong correlation between vibrational dynamics and ultraslow relaxation, we posed the question if the vibrational dynamics encodes predictive information on DH. Here, we review our results, drawn from molecular-dynamics numerical simulation of polymeric and molecular glass-formers, with a special focus on both the breakdown of the Stokes-Einstein relation between diffusion and viscosity, and the size of the regions with correlated displacements.

Mutual information does not detect growing correlations in the propensity of a model molecular liquid

The dynamical spatial correlations detected by the mutual information (MI) in the isoconfigurational particle displacements of a monodisperse molecular viscous liquid are studied via molecular-dynamics simulations by changing considerably both the molecular mobility and the degree of dynamical heterogeneity. Different from atomic liquids, the MI correlation length does not grow on approaching the glass transition by considering the liquid both in full detail as a collection of monomers and as a coarse-grained ensemble of molecular centers of mass. In the detailed picture, it is found that: (i) the MI correlations between monomers are largely due to inter-molecular correlations, (ii) the MI length scale is numerically identical, within the errors, to the correlation length scale of the displacement direction, as drawn by conventional correlation functions. The time evolution of the MI spatial correlations complies with the scaling between the fast vibrational dynamics and the long-time relaxation. Our findings suggest that the characteristics of the MI length scale are markedly system-dependent and not obviously related to dynamical heterogeneity.

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.