Glassy Dynamics Research Articles

Selecting relevant structural features for glassy dynamics by information imbalance.

We numerically investigate the identification of relevant structural features that contribute to the dynamical heterogeneity in a model glass-forming liquid. By employing the recently proposed information imbalance technique, we select these features from a range of physically motivated descriptors. This selection process is performed in a supervised manner (using both dynamical and structural data) and an unsupervised manner (using only structural data). We then apply the selected features to predict future dynamics using a machine learning technique. One of the advantages of the information imbalance technique is that it does not assume any model a priori, i.e., it is a non-parametric method. Finally, we discuss the potential applications of this approach in identifying the dominant mechanisms governing the glassy slow dynamics.

Influence of Density and Pressure on Glass Formation in the Kob-Andersen Model.

We systematically study glass formation in the well-known Kob-Andersen model over a wide range of densities and pressures as a basis for judging the "universality" of glass formation through a comparison to a recent systematic study of model polymeric glass-forming liquids. Our purpose is to establish general characteristics of glass formation, to identify new relations, and to discern which properties of glass-forming liquids are material-specific and which are "universal." To this end, we analyze a number of characteristic properties of glass formation, such as the structural relaxation time, self-diffusion coefficient, viscosity, characteristic temperatures, and fragility. We also consider a suite of properties presumably related to dynamic heterogeneity in an attempt to better understand its relation to structural relaxation. We demonstrate that the glassy dynamics in the Kob-Andersen model exhibit many of the essential trends observed in polymeric glass-forming liquids, pointing to a remarkable "universality" of many aspects of glass formation.

Quantum Liquids: Emergent Higher-Rank Gauge Theory and Fractons

Fractons emerge from many-body systems, featuring subdimensional particles with restricted mobility. These particles have attracted interest for their roles across disciplines, including topological quantum codes, quantum field theory, emergent gravity, and quantum information. They display unique nonequilibrium behaviors such as nonergodicity and glassy dynamics. This review offers a structured overview of fracton phenomena, especially those of gapless fracton liquids, which enable collective modes similar to gauge fluctuations in Maxwell's electromagnetic framework, yet their phenomena are distinguished by a unique conservation law that restricts the mobility of individual charges and monopoles. We delve into the theoretical basis of three-dimensional (3D) fracton liquids, exploring emergent symmetric tensor gauge theories and their properties. We also discuss the material realization of fracton liquids in Yb-based pyrochlore lattices and other synthetic quantum matter platforms.

Observing the two-dimensional Bose glass in an optical quasicrystal

The presence of disorder substantially influences the behaviour of physical systems. It can give rise to slow or glassy dynamics, or to a complete suppression of transport as in Anderson insulators1, where normally extended wavefunctions such as light fields or electronic Bloch waves become exponentially localized. The combined effect of disorder and interactions is central to the richness of condensed-matter physics2. In bosonic systems, it can also lead to additional quantum states such as the Bose glass3,4—an insulating but compressible state without long-range phase coherence that emerges in disordered bosonic systems and is distinct from the well-known superfluid and Mott insulating ground states of interacting bosons. Here we report the experimental realization of the two-dimensional Bose glass using ultracold atoms in an eight-fold symmetric quasicrystalline optical lattice5. By probing the coherence properties of the system, we observe a Bose-glass-to-superfluid transition and map out the phase diagram in the weakly interacting regime. We furthermore demonstrate that it is not possible to adiabatically traverse the Bose glass on typical experimental timescales by examining the capability to restore coherence and discuss the connection to the expected non-ergodicity of the Bose glass. Our observations are in good agreement with recent quantum Monte Carlo predictions6 and pave the way for experimentally testing the connection between the Bose glass, many-body localization and glassy dynamics more generally7,8.

Diffusive Trends in Concentrated Oppositely-Charged Polyelectrolyte Solutions and Onset of Glassy Dynamics.

We utilize single particle tracking studies to investigate the diffusion of polylysine through concentrated matrices of cationic polylysine and anionic polyglutamic acid with no added salts. These studies show that diffusivity has a strong apparently exponential dependence on concentration in crowded systems that does not appear to be a function of the charge sign. These trends are consistent in both single-phase systems prepared at concentrated conditions and polymer-rich coacervate phases formed from dilute phase-separating systems. The likely origin of this behavior is the onset of glassy dynamics spurred by a decrease in plasticization by water and the large excluded volume associated with charge-bearing species. This effect can be contextualized through free volume-based theories such as the Vrentas-Duda model. Overall, we obtain dynamic behavior that is distinctly different from behavior observed in more dilute systems and warrants further investigation.

On the temperature dependence of the density of states of liquids at low energies

We report neutron-scattering measurements of the density of states (DOS) of water and liquid Fomblin in a wide range of temperatures. In the liquid phase, we confirm the presence of a universal low-energy linear scaling of the experimental DOS as a function of the frequency, g(ω)=a(T)ω\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g(\\omega )= a(T) \\omega $$\\end{document}, which persists at all temperatures. The low-frequency scaling of the DOS exhibits a sharp jump at the melting point of water, below which the standard Debye’s law, g(ω)∝ω2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$g(\\omega ) \\propto \\omega ^2$$\\end{document}, is recovered. On the contrary, in Fomblin, we observe a continuous transition between the two exponents reflecting its glassy dynamics, which is confirmed by structure measurements. More importantly, in both systems, we find that the slope a(T) grows with temperature following an exponential Arrhenius-like form, a(T)∝exp(-⟨E⟩/T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$a(T) \\propto \\exp (-\\langle E \\rangle /T)$$\\end{document}. We confirm this experimental trend using molecular dynamics simulations and show that the prediction of instantaneous normal mode (INM) theory for a(T) is in qualitative agreement with the experimental data.

Machine learning that predicts well may not learn the correct physical descriptions of glassy systems

The complexity of glasses makes it challenging to explain their dynamics. Machine learning (ML) has emerged as a promising pathway for understanding glassy dynamics by linking their structural features to rearrangement dynamics. Support vector machine (SVM) was one of the first methods used to detect such correlations. Specifically, a certain output of SVMs trained to predict dynamics from structure, the distance from the separating hyperplane, was interpreted as being linearly related to the activation energy for the rearrangement. By numerical analysis of toy models, we explore under which conditions it is possible to infer the energy barrier to rearrangements from the distance to the separating hyperplane. We observe that such successful inference is possible only under very restricted conditions. Typical tests, such as the apparent Arrhenius dependence of the probability of rearrangement on the inferred energy and the temperature, or high cross-validation accuracy do not guarantee success. Since even in such relatively simple toy models, prediction success of ML models does not necessarily translate into success of learning the underlying physics, we suggest that more careful investigations are needed when such claims are made. For this, we propose practical approaches for measuring the quality of the energy inference and for modifying the inferred model to improve the inference, which should be usable in the context of realistic datasets. Published by the American Physical Society 2024

Advancements in Novel Mechano-Rheological Probes for Studying Glassy Dynamics in Nanoconfined Thin Polymer Films.

The nanoconfinement effects of glassy polymer thin films on their thermal and mechanical properties have been investigated thoroughly, especially with an emphasis on its altered glass transition behavior compared to bulk polymer, which has been known for almost three decades. While research in this direction is still evolving, reaching new heights to unravel the underlying physics of phenomena observed in confined thin polymer films, we have a much clearer picture now. This, in turn, has promoted their application in miniaturized and functional applications. To extract the full potential of such confined films, starting from their fabrication, function, and various applications, we must realize the necessity to have an understanding and availability of robust characterization protocols that specifically target thin film thermo-mechanical stability. Being nanometer-sized in thickness, often atop a solid substrate, direct mechanical testing on such films becomes extremely challenging and often encounters serious complexity from the dominating effect of the substrate. In this review, we have compiled together a few important novel and promising techniques for mechano-rheological characterization of glassy polymer thin films. The conceptual background involved in each technique, constitutive equations, methodology, and current status of research are touched upon following a pedagogical tutorial approach. Further, we discussed each technique's success and limitations, carefully covering the puzzling or contradicting observations reported within the broad nexus of glass transition temperature-viscosity-modulus-molecular mobility (including diffusion and relaxation).

Is the glassy dynamics same in 2D as in 3D? The Adam Gibbs relation test.

It has been recognized of late that even amorphous, glass-forming materials in two dimensions (2D) are affected by Mermin-Wagner-type long wavelength thermal fluctuation, which is inconsequential in three dimensions (3D). We consider the question of whether the effect of spatial dimension on dynamics is only limited to such fluctuations or if the nature of glassy dynamics is intrinsically different in 2D. To address it, we study the relationship between dynamics and thermodynamics using the Adam-Gibbs (AG) relation and the random first order transition (RFOT) theory. Using two model glass-forming liquids, we find that even after removing the effect of long wavelength fluctuations, the AG relation breaks down in two dimensions. Next, we consider the effect of anharmonicity of vibrational entropy-a second factor that affects the thermodynamics but not dynamics. Using the potential energy landscape formalism, we explicitly compute the configurational entropy, both with and without the anharmonic correction. We show that even with both the corrections, the AG relation still breaks down in 2D. The extent of deviation from the AG relation crucially depends on the attractive vs repulsive nature of interparticle interactions, choice of representative timescale (diffusion coefficient vs α-relaxation time), and implies that the RFOT scaling exponents also depend on these factors. Thus, our results suggest that some differences in the nature of glassy dynamics between 2D and 3D remain that are not explained by long wavelength fluctuations.

Measurement of the depth-dependent local dynamics in thin polymer films through rejuvenation of ultrastable glasses

The depth dependence of structural relaxation dynamics is a key part of understanding thin glassy films. Despite this importance and decades of research, a method to provide this information has proved elusive. We measure the isothermal rejuvenation of stable glass films of poly(styrene), and demonstrate that the propagation of the front responsible for the transformation to a supercooled-liquid state serves as a highly localized probe of the local dynamics of the supercooled liquid. We use this connection to probe the depth-dependent relaxation rate with nanometric precision for a series of polystyrene films over a range of temperatures near the bulk glass transition temperature. The analysis shows the spatial extent of enhanced surface mobility and reveals the existence of an unexpected large dynamical length scale in the system. The results are compared with the cooperative-string model for glassy dynamics. The data reveals that the film-thickness dependence of whole film properties arises mainly from the volume fraction of the near-surface region. While the dynamics farthest from the free surface shows the expected bulk-like temperature dependence, the dynamics in the near-surface region shows very little dependence on temperature. This technique can be used in a broad range of thin film materials to gain previously unattainable information about localized structural relaxation.

Extreme dependence of dynamics on concentration in highly crowded polyelectrolyte solutions.

Charge-carrying species, such as polyelectrolytes, are vital to natural and synthetic processes that rely on their dynamic behavior. Through single-particle tracking techniques, the diffusivity of individual polyelectrolyte chains and overall system viscosity are determined for concentrated polylysine solutions. These studies show scaling dependences of D ~ c-6.1 and η ~ c7.2, much stronger than theoretical predictions, drawing the applicability of power law fits into question. Similar trends are observed in concentrated solutions prepared at various pH and counterion conditions. These hindered system dynamics appear universal to polyelectrolyte systems and are attributed to the large effective excluded volumes of polyelectrolyte chains inducing glassy dynamics. The framework of the Vrentas-Duda free-volume theory is used to compare polyelectrolyte and neutral systems. Supported by this theory, excluding counterion mass from total polymer mass results in all environmental conditions collapsing onto a common trendline. These results are applicable to crowded biological systems, such as intracellular environments where protein mobility is strongly inhibited.

Mode coupling behavior and fragile to strong transition of trehalose in a binary mixture with water upon supercooling.

We perform molecular dynamics simulations of a binary mixture of water and trehalose with the TIP4P/Ice water model. We analyze the slow dynamics of trehalose molecules in the mildly supercooled region for concentrations of 3.66 and 18.57 wt. %. We previously studied the dynamics of water in the same mixtures. Supercooled TIP4P/Ice water solvating trehalose molecules was found to follow the Mode Coupling Theory (MCT) and to undergo a transition from a fragile to a strong behavior for both concentrations. Here, we show that also the dynamics of trehalose molecules follows the MCT and displays a fragile to strong crossover (FSC). The results show that trehalose in binary mixtures with water shares with it the dynamical behavior typical of glass forming liquids. Moreover, the FSC for trehalose structural relaxation times is found to occur at temperatures close to those previously obtained for water in the same solutions, showing that the dynamics of the solute is strongly coupled to that of the solvent. We also perform a MCT test showing that the trehalose dynamics obeys the MCT time-temperature superposition principle and that the exponents derived from the theory and the ones obtained from fitting procedure of the relaxation times are comparable, confirming that trehalose molecules in supercooled water solutions follow the MCT of glassy dynamics. Moreover, as predicted by the theory, trehalose particles have MCT parameters comparable to those of water in the same mixtures. This is an important result, given that MCT was originally formulated for monoatomic particles.

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel–Fulcher–Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

A perspective on active glassy dynamics in biological systems

AbstractDynamics is central to living systems. Many experiments in the last two decades have revealed glassy dynamics in diverse biological systems, showing a transition between a solid-like and a fluid-like state. The biological systems have nontrivial characteristics: they are active with novel control parameters and immense complexity. Moreover, glassiness in these systems has many nontrivial features, such as the behavior of dynamical heterogeneity and readily found sub-Arrhenius relaxation dynamics. Theoretical treatments of these systems are generally challenging due to their nonequilibrium nature and large number of control parameters. We first discuss the primary characteristics of a glassy system and then review the experiments that started this field and simulations that have led to a deeper understanding. We also show that despite many challenges in these systems, it has been possible to develop theories that have played a significant role in unifying diverse phenomena and bringing insights. The field is at the interface of physics and biology, freely borrowing tools from both disciplines. We first discuss the known equilibrium scenario and then present the primary changes under activity.

Rotational inertia-induced glassy transition in chiral particle systems.

The dense active matter exhibits characteristics reminiscent of traditional glassy phenomena, yet the role of rotational inertia in glass dynamics remains elusive. In this study, we investigate the glass dynamics of chiral active particles influenced by rotational inertia. Rotational inertia endows exponential memory to particle orientation, restricting its alteration and amplifying the effective persistence time. At lower spinning frequencies, the diffusion coefficient exhibits a peak function relative to rotational inertia for shorter persistence times, while it steadily increases with rotational inertia for longer persistence times. In the realm of high-frequency spinning, the impact of rotational inertia on diffusion behavior becomes more pronounced, resulting in a nonmonotonic and intricate relationship between the diffusion coefficient and rotational inertia. Consequently, the introduction of rotational inertia significantly alters the glassy dynamics of chiral active particles, allowing for the control over transitions between fluid and glassy states by modulating rotational inertia. Moreover, our findings indicate that at a specific spinning temperature, there exists an optimal spinning frequency at which the diffusion coefficient attains its maximum value.

Emergent facilitation and glassy dynamics in supercooled liquids

In supercooled liquids, dynamical facilitation refers to a phenomenon where microscopic motion begets further motion nearby, resulting in spatially heterogeneous dynamics. This is central to the glassy relaxation dynamics of such liquids, which show super-Arrhenius growth of relaxation timescales with decreasing temperature. Despite the importance of dynamical facilitation, there is no theoretical understanding of how facilitation emerges and impacts relaxation dynamics. Here, we present a theory that explains the microscopic origins of dynamical facilitation. We show that dynamics proceeds by localized bond-exchange events, also known as excitations, resulting in the accumulation of elastic stresses with which new excitations can interact. At low temperatures, these elastic interactions dominate and facilitate the creation of new excitations near prior excitations. Using the theory of linear elasticity and Markov processes, we simulate a model, which reproduces multiple aspects of glassy dynamics observed in experiments and molecular simulations, including the stretched exponential decay of relaxation functions, the super-Arrhenius behavior of relaxation timescales as well as their two-dimensional finite-size effects. The model also predicts the subdiffusive behavior of the mean squared displacement (MSD) on short, intermediate timescales. Furthermore, we derive the phonon contributions to diffusion and relaxation, which when combined with the excitation contributions produce the two-step relaxation processes, and the ballistic-subdiffusive-diffusive crossover MSD behaviors commonly found in supercooled liquids.

Quantum fluctuations lead to glassy electron dynamics in the good metal regime of electron doped KTaO3

One of the central challenges in condensed matter physics is to comprehend systems that have strong disorder and strong interactions. In the strongly localized regime, their subtle competition leads to glassy electron dynamics which ceases to exist well before the insulator-to-metal transition is approached as a function of doping. Here, we report on the discovery of glassy electron dynamics deep inside the good metal regime of an electron-doped quantum paraelectric system: KTaO3. We reveal that upon excitation of electrons from defect states to the conduction band, the excess injected carriers in the conduction band relax in a stretched exponential manner with a large relaxation time, and the system evinces simple aging phenomena—a telltale sign of glassy dynamics. Most significantly, we observe a critical slowing down of carrier dynamics below 35 K, concomitant with the onset of quantum paraelectricity in the undoped KTaO3. Our combined investigation using second harmonic generation technique, density functional theory and phenomenological modeling demonstrates quantum fluctuation-stabilized soft polar modes as the impetus for the glassy behavior. This study addresses one of the most fundamental questions regarding the potential promotion of glassiness by quantum fluctuations and opens a route for exploring glassy dynamics of electrons in a well-delocalized regime.

Mean-field analysis of the glassy dynamics of an elastoplastic model of super-cooled liquids

We present a mean-field theory of a coarse-grained model of a super-cooled liquid in which relaxation occurs via local plastic rearrangements. Local relaxation can be induced by thermal fluctuations or by the long-range elastic consequences of other rearrangements. We extract the temperature dependence of both the relaxation time and the length scale of dynamical correlations. We find two dynamical regimes. First, a regime in which the characteristic time and length scales diverge as a power law at a critical temperature T c . This regime is found by an approximation that neglects activated relaxation channels, which can be interpreted as akin to the one found by the mode-coupling transition of glasses. In reality, only a crossover takes place at T c . The residual plastic activity leads to a second regime characterised by an Arrhenius law below T c . In this case, we show that the length scale governing dynamical correlations diverges as a power law as , and is logarithmically related to the relaxation time.

Combined experimental and theoretical studies on glasslike transitions in the frustrated molecular conductors θ−(BEDT−TTF)2MM′(SCN)4

We present results of the coefficient of thermal expansion for the frustrated quasi-two-dimensional molecular conductor θ−(BEDT−TTF)2RbZn(SCN)4 for temperatures 1.5≤T≤290 K. A pronounced first-order phase transition anomaly is observed at the combined charge-order/structural transition at 215 K. Furthermore, clear evidence is found for two separate glasslike transitions at Tg=90–100 K and Tg†=120–130 K, similar to previous findings for θ−(BEDT−TTF)2CsZn(SCN)4 and θ−(BEDT−TTF)2CsCo(SCN)4 reported by Thomas [], both of which lack the charge-order/structural transition. Our findings indicate that these glasslike transitions are common features for the θ−(BEDT−TTF)2MM′(SCN)4 family, with M=(Rb, Cs) and M′=(Co, Zn), irrespective of the presence or absence of charge order. These results are consistent with our model calculations on the glasslike dynamics associated with the flexible ethylene end groups of the BEDT-TTF molecules for various θ−(BEDT−TTF)2MM′(SCN)4 salts, predicting two different conformational glass transitions. Moreover, calculations of the hopping integrals show a substantial degree of dependence on the end groups' conformation, suggesting a significant coupling to the electronic degrees of freedom. Our findings support the possibility that the glassy freezing of the ethylene end groups could drive or enhance glassy charge dynamics. Published by the American Physical Society 2024

Keldysh crossover in one-dimensional Mott insulators

Recent advancements in pulse laser technology have facilitated the exploration of nonequilibrium spectroscopy of electronic states in the presence of strong electric fields across a broad range of photon energies. The Keldysh crossover serves as an indicator that distinguishes between excitations resulting from photon absorption triggered by near-infrared multicycle pulses and those arising from quantum tunneling induced by terahertz pulses. Using a time-dependent density-matrix renormalization group, we investigate the emergence of the Keldysh crossover in a one-dimensional (1D) Mott insulator. We find that the Drude weight is proportional to photo-doped doublon density when a pump pulse induces photon absorption. In contrast, the Drude weight is suppressed when a terahertz pulse introduces doublons and holons via quantum tunneling. The suppressed Drude weight accompanies glassy dynamics with suppressed diffusion, which is a consequence of strong correlations and exhibits finite polarization decaying slowly after pulse irradiation. In the quantum tunneling region, entanglement entropy slowly grows logarithmically. These contrasting behaviors between the photon-absorption and quantum tunneling regions are a manifestation of the Keldysh crossover in 1D Mott insulators and provide a novel methodology for designing the localization and symmetry of electronic states called subcycle-pulse engineering.