Vibrational modes identify soft spots in a sheared disordered packing

Long-range correlations in elastic moduli and local stresses at the unjamming transition

We explore the behaviour of spatially heterogeneous elastic moduli as well as the correlations between local moduli in model solids with short-range repulsive potentials. We show through numerical simulations that local elastic moduli exhibit long-range correlations, similar to correlations in the local stresses. Specifically, the correlations in local shear moduli exhibit anisotropic behavior at large lengthscales characterized by pinch-point singularities in Fourier space, displaying a structural pattern akin to shear stress correlations. Focussing on two-dimensional jammed solids approaching the unjamming transition, we show that stress correlations exhibit universal properties, characterized by a quadratic p2 dependence of the correlations as the pressure p approaches zero, independent of the details of the model. In contrast, the modulus correlations exhibit a power-law dependence with different exponents depending on the specific interaction potential. Furthermore, we illustrate that while affine responses lack long-range correlations, the total modulus, which encompasses non-affine behavior, exhibits long-range correlations.

Predicting the pathways of string-like motions in metallic glasses via path-featurizing graph neural networks

String-like motions (SLMs)-cooperative, "snake"-like movements of particles-are crucial for dynamics in diverse glass formers. Despite their ubiquity, questions persist: Do SLMs prefer specific paths? If so, can we predict these paths? Here, in Al-Sm glasses, our isoconfigurational ensemble simulations reveal that SLMs do follow certain paths. By designing a graph neural network (GNN) to featurize the environment around directional paths, we achieve a high-fidelity prediction of likely SLM pathways, solely based on the static structure. GNN gauges a structural measure to assess each path's propensity to engage in SLMs, akin to a "softness" metric, but for paths rather than for atoms. Our GNN interpretation reveals the critical role of the bottleneck zone along a path in steering SLMs. By monitoring "path softness," we elucidate that SLM-favored paths transit from fragmented to interconnected upon glass transition. Our findings reveal that, beyond atoms or clusters, glasses have another dimension of structural heterogeneity: "paths."

Information decomposition in complex systems via machine learning

One of the fundamental steps toward understanding a complex system is identifying variation at the scale of the system's components that is most relevant to behavior on a macroscopic scale. Mutual information provides a natural means of linking variation across scales of a system due to its independence of functional relationship between observables. However, characterizing the manner in which information is distributed across a set of observables is computationally challenging and generally infeasible beyond a handful of measurements. Here, we propose a practical and general methodology that uses machine learning to decompose the information contained in a set of measurements by jointly optimizing a lossy compression of each measurement. Guided by the distributed information bottleneck as a learning objective, the information decomposition identifies the variation in the measurements of the system state most relevant to specified macroscale behavior. We focus our analysis on two paradigmatic complex systems: a Boolean circuit and an amorphous material undergoing plastic deformation. In both examples, the large amount of entropy of the system state is decomposed, bit by bit, in terms of what is most related to macroscale behavior. The identification of meaningful variation in data, with the full generality brought by information theory, is made practical for studying the connection between micro- and macroscale structure in complex systems.

Rearrangements during slow compression of a jammed two-dimensional emulsion

As amorphous materials get jammed, both geometric and dynamic heterogeneity are observed. We investigate the correlation between the local geometric heterogeneity and local rearrangements in a slowly compressed bidisperse quasi-two-dimensional emulsion system. The compression is driven by evaporation of the continuous phase and causes the area packing fraction to increase from 0.88 to 0.99. We quantify the structural heterogeneity of the system using the radical Voronoi tessellation following the method of Rieser et al. [Phys. Rev. Lett. 116, 088001 (2016)]0031-900710.1103/PhysRevLett.116.088001. We define two structural quantities characterizing local structure, the first of which considers nearest neighbors and the second of which includes information from second-nearest neighbors. We find that droplets in heterogeneous local regions are more likely to have local rearrangements. These rearrangements are generally T1 events where two droplets converge toward a void, and two droplets move away from the void to make room for the converging droplets. Thus, the presence of the voids tends to orient the T1 events. The presence of a correlation between the structural quantities and the rearrangement dynamics remains qualitatively unchanged over the entire range of packing fractions observed.

Low-frequency vibrational density of states of ordinary and ultra-stable glasses

A remarkable feature of disordered solids distinct from crystals is the violation of the Debye scaling law of the low-frequency vibrational density of states. Because the low-frequency vibration is responsible for many properties of solids, it is crucial to elucidate it for disordered solids. Numerous recent studies have suggested power-law scalings of the low-frequency vibrational density of states, but the scaling exponent is currently under intensive debate. Here, by classifying disordered solids into stable and unstable ones, we find two distinct and robust scaling exponents for non-phononic modes at low frequencies. Using the competition of these two scalings, we clarify the variation of the scaling exponent and hence reconcile the debate. Via the study of both ordinary and ultra-stable glasses, our work reveals a comprehensive picture of the low-frequency vibration of disordered solids and sheds light on the low-frequency vibrational features of ultra-stable glasses on approaching the ideal glass.

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.

Solid-that-Flows Picture of Glass-Forming Liquids

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

Instantaneous normal modes of glass-forming liquids during the athermal relaxation process of the steepest descent algorithm

Understanding glass formation by quenching remains a challenge in soft condensed matter physics. Recent numerical studies on steepest descent dynamics, which is one of the simplest models of quenching, revealed that quenched liquids undergo slow relaxation with a power law towards mechanical equilibrium and that the late stage of this process is governed by local rearrangements of particles. These advances motivate the detailed study of instantaneous normal modes during the relaxation process because the glassy dynamics is considered to be governed by stationary points of the potential energy landscape. Here, we performed a normal mode analysis of configurations during the steepest descent dynamics and found that the dynamics is driven by almost flat directions of the potential energy landscape at long times. These directions correspond to localized modes and we characterized them in terms of their statistics and structure using methods developed in the study of local minima of the potential energy landscape.

Development of broad modulus profile upon polymer-polymer interface formation between immiscible glassy-rubbery domains

Interfaces of glassy materials such as thin films, blends, and composites create strong unidirectional gradients to the local heterogeneous dynamics that can be used to elucidate the length scales and mechanisms associated with the dynamic heterogeneity of glasses. We focus on bilayer films of two different polymers with very different glass transition temperatures ([Formula: see text]) where previous work has demonstrated a long-range (∼200 nm) profile in local [Formula: see text] is established between immiscible glassy and rubbery polymer domains when the polymer-polymer interface is formed to equilibrium. Here, we demonstrate that an equally long-ranged gradient in local modulus [Formula: see text] is established when the polymer-polymer interface ([Formula: see text]5 nm) is formed between domains of glassy polystyrene (PS) and rubbery poly(butadiene) (PB), consistent with previous reports of a broad [Formula: see text] profile in this system. A continuum physics model for the shear wave propagation caused by a quartz crystal microbalance across a PB/PS bilayer film is used to measure the viscoelastic properties of the bilayer during the evolution of the PB/PS interface showing the development of a broad gradient in local modulus [Formula: see text] spanning [Formula: see text]180 nm between the glassy and rubbery domains of PS and PB. We suggest these broad profiles in [Formula: see text] and [Formula: see text] arise from a coupling of the spectrum of vibrational modes across the polymer-polymer interface as a result of acoustic impedance matching of sound waves with [Formula: see text] nm during interface broadening that can then trigger density fluctuations in the neighboring domain.

Thawed matrix method for computing local mechanical properties of amorphous solids

We present a method for computing locally varying nonlinear mechanical properties in particle simulations of amorphous solids. Plastic rearrangements outside a probed region are suppressed by introducing an external field that directly penalizes large nonaffine displacements. With increasing strength of the field, plastic deformation can be localized. We characterize the distribution of local plastic yield stresses (residual local stresses to instability) with our approach and assess the correlation of their spatial maps with plastic activity in a model two-dimensional amorphous solid. Our approach reduces artifacts inherent in a previous method known as the "frozen matrix" approach that enforces fully affine deformation and improves the prediction of plastic rearrangements from structural information.

Athermal quasistatic cavitation in amorphous solids: Effect of random pinning

Amorphous solids are known to fail catastrophically via fracture, and cavitation at nano-metric scales is known to play a significant role in such a failure process. Micro-alloying via inclusions is often used as a means to increase the fracture toughness of amorphous solids. Modeling such inclusions as randomly pinned particles that only move affinely and do not participate in plastic relaxations, we study how the pinning influences the process of cavitation-driven fracture in an amorphous solid. Using extensive numerical simulations and probing in the athermal quasistatic limit, we show that just by pinning a very small fraction of particles, the tensile strength is increased, and also the cavitation is delayed. Furthermore, the cavitation that is expected to be spatially heterogeneous becomes spatially homogeneous by forming a large number of small cavities instead of a dominant cavity. The observed behavior is rationalized in terms of screening of plastic activity via the pinning centers, characterized by a screening length extracted from the plastic-eigenmodes.

Effects of polydispersity on the plastic behaviors of dense two-dimensional granular systems under shear

We study particle-scale motion in sheared highly polydisperse amorphous materials, in which the largest particles are as much as ten times the size of the smallest. We find strikingly different behavior from the more commonly studied amorphous systems with low polydispersity. In particular, an analysis of the nonaffine motion of particles reveals qualitative differences between large and small particles: The smaller particles have dramatically more nonaffine motion, which is induced by the presence of the large particles. We characterize how the nonaffine motion changes from the low- to high-polydispersity regimes. We further demonstrate a quantitative way to distinguish between "large" and "small" particles in systems with broad distributions of particle sizes. A macroscopic consequence of the nonaffine motion is a decrease in the energy dissipation rate for highly polydisperse samples, which is due both to a geometric consequence of the changing jamming conditions for higher polydispersity and to the changing character of nonaffine motion.

Detecting low-energy quasilocalized excitations in computer glasses

Soft, quasilocalized excitations (QLEs) are known to generically emerge in a broad class of disordered solids and to govern many facets of the physics of glasses, from wave attenuation to plastic instabilities. In view of this key role of QLEs, shedding light upon several open questions in glass physics depends on the availability of computational tools that allow one to study QLEs' statistical mechanics. The latter is a formidable task since harmonic analyses are typically contaminated by hybridizations of QLEs with phononic excitations at low frequencies, obscuring a clear picture of QLEs' abundance, typical frequencies, and other important micromechanical properties. Here we present an efficient algorithm to detect the field of quasilocalized excitations in structural computer glasses. The algorithm introduced takes a computer-glass sample as input and outputs a library of QLEs embedded in that sample. We demonstrate the power of the algorithm by reporting the spectrum of glassy excitations in two-dimensional computer glasses featuring a huge range of mechanical stability, which is inaccessible using conventional harmonic analyses due to phonon hybridizations. Future applications are discussed.

Relevance of structural defects to the mechanism of mechanical deformation in metallic glasses

It is known that deformation in disordered materials such as metallic glasses and supercooled liquids occurs via the cooperative rearrangement of atoms or constituent particles at dynamical heterogeneities, commonly regarded as point-like defects. We show via molecular-dynamics simulations that there is no apparent relationship between atomic rearrangements and the local atomic environment as measured by the atomic-level stresses, kinetic and potential energies, and the per-atom Voronoi volume. In addition, there is only a weak correlation between atomic rearrangements and the largest and smallest eigenvalues of the dynamical matrix. Our results confirm the transient nature of dynamical heterogeneities and suggest that the notion of defects may be less relevant than that of a propensity for rearrangement.

Enhanced vibrational stability in glass droplets

We show through simulations of amorphous solids prepared in open-boundary conditions that they possess significantly fewer low-frequency vibrational modes compared to their periodic boundary counterparts. Specifically, using measurements of the vibrational density of states, we find that the D ( ω ) ∼ ω 4 law changes to D ( ω ) ∼ ω δ with δ ≈ 5 in two dimensions and δ ≈ 4.5 in three dimensions. Crucially, this enhanced stability is achieved when utilizing slow annealing protocols to generate solid configurations. We perform an anharmonic analysis of the minima corresponding to the lowest frequency modes in such open-boundary systems and discuss their correlation with the density of states. A study of various system sizes further reveals that small systems display a higher degree of localization in vibrations. Lastly, we confine open-boundary solids in order to introduce macroscopic stresses in the system, which are absent in the unconfined system and find that the D ( ω ) ∼ ω 4 behavior is recovered.

Designing the pressure-dependent shear modulus using tessellated granular metamaterials

Jammed packings of granular materials display complex mechanical response. For example, the ensemble-averaged shear modulus 〈G〉 increases as a power law in pressure p for static packings of soft spherical particles that can rearrange during compression. We seek to design granular materials with shear moduli that can either increase or decrease with pressure without particle rearrangements even in the large-system limit. To do this, we construct tessellated granular metamaterials by joining multiple particle-filled cells together. We focus on cells that contain a small number of bidisperse disks in two dimensions. We first study the mechanical properties of individual disk-filled cells with three types of boundaries: periodic boundary conditions (PBC), fixed-length walls (FXW), and flexible walls (FLW). Hypostatic jammed packings are found for cells with FLW, but not in cells with PBC and FXW, and they are stabilized by quartic modes of the dynamical matrix. The shear modulus of a single cell depends linearly on p. We find that the slope of the shear modulus with pressure λ_{c}<0 for all packings in single cells with PBC where the number of particles per cell N≥6. In contrast, single cells with FXW and FLW can possess λ_{c}>0, as well as λ_{c}<0, for N≤16. We show that we can force the mechanical properties of multicell granular metamaterials to possess those of single cells by constraining the end points of the outer walls and enforcing an affine shear response. These studies demonstrate that tessellated granular metamaterials provide a platform for the design of soft materials with specified mechanical properties.

Topology of vibrational modes predicts plastic events in glasses

The plastic deformation of crystalline materials can be understood by considering their structural defects such as disclinations and dislocations. Although also glasses are solids, their structure resembles closely the one of a liquid and hence the concept of structural defects becomes ill-defined. As a consequence it is very challenging to rationalize on a microscopic level the mechanical properties of glasses close to the yielding point and to relate plastic events to structural properties. Here we investigate the topological characteristics of the eigenvector field of the vibrational excitations of a two-dimensional glass model, notably the geometric arrangement of the topological defects as a function of vibrational frequency. We find that if the system is subjected to a quasistatic shear, the location of the resulting plastic events correlate strongly with the topological defects that have a negative charge. Our results provide thus a direct link between the structure of glasses prior their deformation and the plastic events during deformation.

Theory of rigidity and numerical analysis of density of states of two-dimensional amorphous solids with dispersed frictional grains in the linear response regime

Using the Jacobian matrix, we obtain a theoretical expression of rigidity and the density of states of two-dimensional amorphous solids consisting of frictional grains in the linear response to an infinitesimal strain, in which we ignore the dynamical friction caused by the slip processes of contact points. The theoretical rigidity agrees with that obtained by molecular dynamics simulations. We confirm that the rigidity is smoothly connected to the value in the frictionless limit. We find that there are two modes in the density of states for sufficiently small k_{T}/k_{N}, which is the ratio of the tangential to normal stiffness. Rotational modes exist at low frequencies or small eigenvalues, whereas translational modes exist at high frequencies or large eigenvalues. The location of the rotational band shifts to the high-frequency region with an increase in k_{T}/k_{N} and becomes indistinguishable from the translational band for large k_{T}/k_{N}.

Local and global measures of the shear moduli of jammed disk packings

Strain-controlled isotropic compression gives rise to jammed packings of repulsive, frictionless disks with either positive or negative global shear moduli. We carry out computational studies to understand the contributions of the negative shear moduli to the mechanical response of jammed disk packings. We first decompose the ensemble-averaged, global shear modulus as 〈G〉=(1-F_{-})〈G_{+}〉+F_{-}〈G_{-}〉, where F_{-} is the fraction of jammed packings with negative shear moduli and 〈G_{+}〉 and 〈G_{-}〉 are the average values from packings with positive and negative moduli, respectively. We show that 〈G_{+}〉 and 〈|G_{-}|〉 obey different power-law scaling relations above and below pN^{2}∼1. For pN^{2}>1, both 〈G_{+}〉N and 〈|G_{-}|〉N∼(pN^{2})^{β}, where β∼0.5 for repulsive linear spring interactions. Despite this, 〈G〉N∼(pN^{2})^{β^{'}} with β^{'}≳0.5 due to the contributions from packings with negative shear moduli. We show further that the probability distribution of global shear moduli P(G) collapses at fixed pN^{2} and different values of p and N. We calculate analytically that P(G) is a Γ distribution in the pN^{2}≪1 limit. As pN^{2} increases, the skewness of P(G) decreases and P(G) becomes a skew-normal distribution with negative skewness in the pN^{2}≫1 limit. We also partition jammed disk packings into subsystems using Delaunay triangulation of the disk centers to calculate local shear moduli. We show that the local shear moduli defined from groups of adjacent triangles can be negative even when G>0. The spatial correlation function of local shear moduli C(r[over ⃗]) displays weak correlations for pn_{sub}^{2}<10^{-2}, where n_{sub} is the number of particles within each subsystem. However, C(r[over ⃗]) begins to develop long-ranged spatial correlations with fourfold angular symmetry for pn_{sub}^{2}≳10^{-2}.

Eigenvalue analysis of stress-strain curve of two-dimensional amorphous solids of dispersed frictional grains with finite shear strain

The stress-strain curve of two-dimensional frictional dispersed grains interacting with a harmonic potential without considering the dynamical slip under a finite strain is determined by using eigenvalue analysis of the Hessian matrix. After the configuration of grains is obtained, the stress-strain curve based on the eigenvalue analysis is in almost perfect agreement with that obtained by the simulation, even if there are plastic deformations caused by stress avalanches. Unlike the naive expectation, the eigenvalues in our model do not indicate any precursors to the stress-drop events.

Gas-liquid phase separation at zero temperature: mechanical interpretation and implications for gelation

The relationship between glasses and gels has been intensely debated for decades; however, the transition between these two phases remains elusive. To investigate a gel formation process in the zero-temperature limit and its relation to the glass phase, we conducted numerical experiments on athermal quasistatic decompression. During decompression, the system experiences a cavitation event similar to phase separation and this is a gelation process at zero temperature. A normal mode analysis revealed that the phase separation is signaled by the vanishing of the lowest eigenenergy, similar to plastic events of glasses under shear. One primary difference from the shear-induced plasticity is that the vanishing mode experiences a qualitative change in its spatial energy distribution at the phase separation point. These findings enable us to define the glass-gel phase boundary based on mechanics.

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.

Predicting Defects in Soft Sphere Packings near Jamming Using the Force Network Ensemble

Amorphous systems of soft particles above jamming have more contacts than are needed to achieve mechanical equilibrium. The force network of a granular system with a fixed contact network is thus underdetermined and can be characterized as a random instantiation within the space of the force network ensemble. In this Letter, we show that defect contacts that are not necessary for stability of the system can be uniquely identified by examining the boundaries of this space of allowed force networks. We further show that, for simulations in the near jamming limit, this identification is nearly always correct and that defect contacts are broken under decompression of the system.

Mechanical annealing and memories in a disordered solid

Shearing a disordered or amorphous solid for many cycles with a constant strain amplitude can anneal it, relaxing a sample to a steady state that encodes a memory of that amplitude. This steady state also features a remarkable stability to amplitude variations that allows one to read the memory. Here, we shed light on both annealing and memory by considering how to mechanically anneal a sample to have as little memory content as possible. In experiments, we show that a "ring-down" protocol reaches a comparable steady state but with no discernible memories and minimal structural anisotropy. We introduce a method to characterize the population of rearrangements within a sample and show how it connects with the response to amplitude variation and the size of annealing steps. These techniques can be generalized to other forms of glassy matter and a wide array of disordered solids, especially those that yield by flowing homogeneously.

Density of states below the first sound mode in 3D glasses

Glasses feature universally low-frequency excess vibrational modes beyond Debye prediction, which could help rationalize, e.g., the glasses’ unusual temperature dependence of thermal properties compared to crystalline solids. The way the density of states of these low-frequency excess modes D( ω) depends on the frequency ω has been debated for decades. Recent simulation studies of 3D glasses suggest that D( ω) scales universally with ω4 in a low-frequency regime below the first sound mode. However, no simulation study has ever probed as low frequencies as possible to test directly whether this quartic law could work all the way to extremely low frequencies. Here, we calculated D( ω) below the first sound mode in 3D glasses over a wide range of frequencies. We find D( ω) scales with ω β with β < 4 at very low frequencies examined, while the ω4 law works only in a limited intermediate-frequency regime in some glasses. Moreover, our further analysis suggests our observation does not depend on glass models or glass stabilities examined. The ω4 law of D( ω) below the first sound mode is dominant in current simulation studies of 3D glasses, and our direct observation of the breakdown of the quartic law at very low frequencies thus leaves an open but important question that may attract more future numerical and theoretical studies.

Correlation of plastic events with local structure in jammed packings across spatial dimensions

SignificanceMean-field theories, exact in the limit of infinite spatial dimensions, succeed in describing many features of glasses and amorphous solids in low dimensions, leading to considerable effort to understand how behavior evolves with dimension. Until now, all evidence has supported a picture in which "localized physics," responsible for deviations from mean-field behavior in low dimensions, fades away with rising dimension. Our work shows that rearrangements, in which particles change relative positions leading to fluid-like response, reveal a different picture of dimensional cross-over. We find that rearrangements, which are localized in two- and three-dimensional systems and correlated with local structure, remain just as correlated with local structure up to five dimensions, suggesting that local structure is important even in high dimensions.

Cooperative effects driving the multi-periodic dynamics of cyclically sheared amorphous solids

When subject to cyclic forcing, amorphous solids can reach periodic, repetitive states, where the system behaves plastically, but the particles return to their initial positions after one or more forcing cycles, where the latter response is called multi-periodic. It is known that plasticity in amorphous materials is mediated by local rearrangements called ``soft spots' or ``shear transformation zones'.Experiments and simulations indicate that soft spots can be modeled as hysteretic two-state entities interacting via quadrupolar displacement fields generated when they switch states and that these interactions can give rise to multi-periodic behavior. However, how interactions facilitate multi-periodicity is unknown. Here we show, using a model of random interacting two-state systems and molecular dynamics simulations, that multi-periodicity arises from oscillations in the magnitudes of the switching field of soft spots which cause soft spots to be active during some forcing cycles and idle during others. We demonstrate that these oscillations result from cooperative effects facilitated by the frustrated interactions between the soft spots. The presence of such mechanisms has implications for manipulating memory in frustrated hysteretic systems.

Identifying structural signature of dynamical heterogeneity via the local softness parameter

In this work we study the relationship between the softness of a mean-field caging potential and dynamics at the local level. We first describe the local softness, which shows a distribution, thus identifying structural heterogeneity. We show that the lifetime of the softness parameter is connected to the lifetime of the well-known cage structure in supercooled liquids. Finally, our theory predicts that the local softness and the local dynamics is causal below the onset temperature where there is a decoupling between the short and long time dynamics, thus allowing a static description of the cage. With the decrease in temperature, the correlation between structure and dynamics increases. The study shows that at lower temperatures, the structural heterogeneity increases, and since the structure becomes a better predictor of the dynamics, it leads to an increase in the dynamical heterogeneity. We also find that the softness of a hard, immobile region evolves with time and becomes soft and eventually mobile due to the rearrangements in the neighborhood, confirming the well-known facilitation effect.

Avalanche dynamics in sheared athermal particle packings occurs via localized bursts predicted by unstable linear response

Under applied shear strain, granular and amorphous materials deform via particle rearrangements, which can be small and localized or organized into system-spanning avalanches. While the statistical properties of avalanches under quasi-static shear are well-studied, the dynamics during avalanches is not. In numerical simulations of sheared soft spheres, we find that avalanches can be decomposed into bursts of localized deformations, which we identify using an extension of persistent homology methods. We also study the linear response of unstable systems during an avalanche, demonstrating that eigenvalue dynamics are highly complex during such events, and that the most unstable eigenvector is a poor predictor of avalanche dynamics. Instead, we modify existing tools that identify localized excitations in stable systems, and apply them to these unstable systems with non-positive definite Hessians, quantifying the evolution of such excitations during avalanches. We find that bursts of localized deformations in the avalanche almost always occur at localized excitations identified using the linear spectrum. These new tools will provide an improved framework for validating and extending mesoscale elastoplastic models that are commonly used to explain avalanche statistics in glasses and granular matter.

Role of Local Structure in the Enhanced Dynamics of Deformed Glasses

External stress can accelerate molecular mobility of amorphous solids by several orders of magnitude. The changes in mobility are commonly interpreted through the Eyring model, which invokes an empirical activation volume. Here, we analyze constant-stress molecular dynamics simulations and propose a structure-dependent Eyring model, connecting activation volume to a machine-learned field, softness. We show that stress has a heterogeneous effect on the mobility that depends on local structure through softness. The barrier impeding relaxation reduces more for well-packed particles, which explains the narrower distribution of relaxation time observed under stress.

Understanding slow compression and decompression of frictionless soft granular matter by network analysis

We consider dense granular systems in three spatial dimensions exposed to slow compression and decompression, below, during, above and well above jamming. The evolution of granular systems under slow deformation is non-trivial and involves smooth, continuous, reversible (de)compression periods, interrupted by fast, discontinuous, irreversible transition events. These events are often, but not always, associated with rearrangements of particles and of the contact network. How many particles are involved in these transitions between two states can range from few to almost all in the system. An analysis of the force network that is built on top of the contact network is carried out using the tools of persistent homology. Results involve the observation that kinetic energy is correlated with the intensity of rearrangements, while the evolution of global mechanical measures, such as pressure, is strongly correlated with the evolution of the topological measures quantifying loops in the force network. Surprisingly, some transitions are clearly detected by persistent homology even though motion/rearrangement of particles is much weaker, i.e., much harder to detect or, in some cases, not observed at all.

Searching for structural predictors of plasticity in dense active packings

In amorphous solids subject to shear or thermal excitation, so-called structural indicators have been developed that predict locations of future plasticity or particle rearrangements. An open question is whether similar tools can be used in dense active materials, but a challenge is that under most circumstances, active systems do not possess well-defined solid reference configurations. We develop a computational model for a dense active crowd attracted to a point of interest, which does permit a mechanically stable reference state in the limit of infinitely persistent motion. Previous work on a similar system suggested that the collective motion of crowds could be predicted by inverting a matrix of time-averaged two-particle correlation functions. Seeking a first-principles understanding of this result, we demonstrate that this active matter system maps directly onto a granular packing in the presence of an external potential, and extend an existing structural indicator based on linear response to predict plasticity in the presence of noisy dynamics. We find that the strong pressure gradient necessitated by the directed activity, as well as a self-generated free boundary, strongly impact the linear response of the system. In low-pressure regions the linear-response-based indicator is predictive, but it does not work well in the high-pressure interior of our active packings. Our findings motivate and inform future work that could better formulate structure-dynamics predictions in systems with strong pressure gradients.

Quantifying local rearrangements in three-dimensional granular materials: Rearrangement measures, correlations, and relationship to stresses

Quantifying the ways in which local particle rearrangements contribute to macroscopic plasticity is one of the fundamental pursuits of granular mechanics and soft matter physics. Here we examine local rearrangements that occur naturally during the deformation of three samples of 3D granular materials subjected to distinct boundary conditions by employing in situ x-ray measurements of particle-resolved structure and stress. We focus on five distinct rearrangement measures, their statistics, interrelationships, contributions to macroscopic deformation, repeatability, and dependence on local structure and stress. Our most significant findings are that local rearrangements (1) are correlated on a scale of three to four particle diameters, (2) exhibit volumetric strain-shear strain and nonaffine displacement-rotation coupling, (3) exhibit correlations that suggest either rearrangement repeatability or that rearrangements span multiple steps of incremental sample strain, and (4) show little dependence on local stress but correlate with quantities describing local structure, such as porosity. Our results are presented in the context of relevant plasticity theories and are consistent with recent findings suggesting that local structure may play at least as important of a role as local stress in determining the nature of local rearrangements.

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of LJ glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale ξ_s of the structures grows as the density is lowered. Second, the LJ gels are extremely soft, with both shear G and bulk K moduli being orders of magnitude smaller than those of LJ glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency ω_* in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by ξ_s, G, K, and ω_*, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout this work, we also reveal that LJ gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

Low-Frequency Excess Vibrational Modes in Two-Dimensional Glasses

Glasses possess more low-frequency vibrational modes than predicted by Debye theory. These excess modes are crucial for the understanding of the low temperature thermal and mechanical properties of glasses, which differ from those of crystalline solids. Recent simulational studies suggest that the density of the excess modes scales with their frequency ω as ω^{4} in two and higher dimensions. Here, we present extensive numerical studies of two-dimensional model glass formers over a large range of glass stabilities. We find that the density of the excess modes follows D_{exc}(ω)∼ω^{2} up to around the boson peak, regardless of the glass stability. The stability dependence of the overall scale of D_{exc}(ω) correlates with the stability dependence of low-frequency sound attenuation. However, we also find that, in small systems, where the first sound mode is pushed to higher frequencies, at frequencies below the first sound mode, there are excess modes with a system size independent density of states that scales as ω^{3}.

The initiation of shear band formation in deformed metallic glasses from soft localized domains

It has long been thought that shear band (SB) formation in amorphous solids initiates from relatively "soft" regions in the material in which large-scale non-affine deformations become localized. The test of this hypothesis requires an effective means of identifying "soft" regions and their evolution as the material is deformed to varying degrees, where the metric of "softness" must also account for the effect of temperature on local material stiffness. We show that the mean square atomic displacement on a caging timescale ⟨u²⟩, the "Debye-Waller factor," provides a useful method for estimating the shear modulus of the entire material and, by extension, the material stiffness at an atomic scale. Based on this "softness" metrology, we observe that SB formation indeed occurs through the strain-induced formation of localized soft regions in our deformed metallic glass free-standing films. Unexpectedly, the critical strain condition for SB formation occurs when the softness (⟨u²⟩) distribution within the emerging soft regions approaches that of the interfacial region in its undeformed state, initiating an instability with similarities to the transition to turbulence. Correspondingly, no SBs arise when the material is so thin that the entire material can be approximately described as being "interfacial" in nature. We also quantify relaxation in the glass and the nature and origin of highly non-Gaussian particle displacements in the dynamically heterogeneous SB regions at times longer than the caging time.

Low-energy quasilocalized excitations in structural glasses

Glassy solids exhibit a wide variety of generic thermomechanical properties, ranging from universal anomalous specific heat at cryogenic temperatures to nonlinear plastic yielding and failure under external driving forces, which qualitatively differ from their crystalline counterparts. For a long time, it has been believed that many of these properties are intimately related to nonphononic, low-energy quasilocalized excitations (QLEs) in glasses. Indeed, recent computer simulations have conclusively revealed that the self-organization of glasses during vitrification upon cooling from a melt leads to the emergence of such QLEs. In this Perspective, we review developments over the past three decades toward understanding the emergence of QLEs in structural glasses and the degree of universality in their statistical and structural properties. We discuss the challenges and difficulties that hindered progress in achieving these goals and review the frameworks put forward to overcome them. We conclude with an outlook on future research directions and open questions.

Unifying Description of the Vibrational Anomalies of Amorphous Materials

The vibrational density of states D(ω) of solids controls their thermal and transport properties. In crystals, the low-frequency modes are extended phonons distributed in frequency according to Debye's law, D(ω)∝ω^{2}. In amorphous solids, phonons are damped, and at low frequency D(ω) comprises extended modes in excess over Debye's prediction, leading to the so-called boson peak in D(ω)/ω^{2} at ω_{bp}, and quasilocalized ones. Here we show that boson peak and phonon attenuation in the Rayleigh scattering regime are related, as suggested by correlated fluctuating elasticity theory, and that amorphous materials can be described as homogeneous isotropic elastic media punctuated by quasilocalized modes acting as elastic heterogeneities. Our numerical results resolve the conflict between theoretical approaches attributing amorphous solids' vibrational anomalies to elastic disorder and localized defects.

Bond-space operator disentangles quasilocalized and phononic modes in structural glasses

The origin of several emergent mechanical and dynamical properties of structural glasses is often attributed to populations of localized structural instabilities, coined quasilocalized modes (QLMs). Under a restricted set of circumstances, glassy QLMs can be revealed by analyzing computer glasses' vibrational spectra in the harmonic approximation. However, this analysis has limitations due to system-size effects and hybridization processes with low-energy phononic excitations (plane waves) that are omnipresent in elastic solids. Here we overcome these limitations by exploring the spectrum of a linear operator defined on the space of particle interactions (bonds) in a disordered material. We find that this bond-force-response operator offers a different interpretation of QLMs in glasses and cleanly recovers some of their important statistical and structural features. The analysis presented here reveals the dependence of the number density (per frequency) and spatial extent of QLMs on material preparation protocol (annealing). Finally, we discuss future research directions and possible extensions of this work.

Marginal stability in memory training of jammed solids

Memory encoding by cyclic shear is a reliable process to store information in jammed solids, yet its underlying mechanism and its connection to the amorphous structure are not fully understood. When a jammed sphere packing is repeatedly sheared with cycles of the same strain amplitude, it optimizes its mechanical response to the cyclic driving and stores a memory of it. We study memory by cyclic shear training as a function of the underlying stability of the amorphous structure in marginally stable and highly stable packings, the latter produced by minimizing the potential energy using both positional and radial degrees of freedom. We find that jammed solids need to be marginally stable in order to store a memory by cyclic shear. In particular, highly stable packings store memories only after overcoming brittle yielding and the cyclic shear training takes place in the shear band, a region which we show to be marginally stable.

Using delaunay triangularization to characterize non-affine displacement fields during athermal, quasistatic deformation of amorphous solids

We investigate the non-affine displacement fields that occur in two-dimensional Lennard-Jones models of metallic glasses subjected to athermal, quasistatic simple shear (AQS). During AQS, the shear stress versus strain displays continuous quasi-elastic segments punctuated by rapid drops in shear stress, which correspond to atomic rearrangement events. We capture all information concerning the atomic motion during the quasi-elastic segments and shear stress drops by performing Delaunay triangularizations and tracking the deformation gradient tensor F_α associated with each triangle α. To understand the spatio-temporal evolution of the displacement fields during shear stress drops, we calculate F_α along minimal energy paths from the mechanically stable configuration immediately before to that after the stress drop. We find that quadrupolar displacement fields form and dissipate both during the quasi-elastic segments and shear stress drops. We then perform local perturbations (rotation, dilation, simple and pure shear) to single triangles and measure the resulting displacement fields. We find that local pure shear deformations of single triangles give rise to mostly quadrupolar displacement fields, and thus pure shear strain is the primary type of local strain that is activated by bulk, athermal quasistatic simple shear. Other local perturbations, e.g. rotations, dilations, and simple shear of single triangles, give rise to vortex-like and dipolar displacement fields that are not frequently activated by bulk AQS. These results provide fundamental insights into the non-affine atomic motion that occurs in driven, glassy materials.

Atomic vibration as an indicator of the propensity for configurational rearrangements in metallic glasses

In a metallic glass (MG), the propensity for atomic rearrangements varies spatially from location to location in the amorphous solid, making the prediction of their likelihood a major challenge. One can attack this problem from the "structure controls properties" standpoint. But all the current structure-centric parameters are mostly based on local atomic packing information limited to short-range order, hence falling short in reliably forecasting how the local region would respond to external stimuli (e.g., temperature and/or stress). Alternatively, one can use indicators informed by physical properties to bridge the static structure on the one hand, and the response of the local configuration on the other. A sub-group of such physics-informed quantities consists of atomic vibration parameters, which will be singled out as the focus of this article. Here we use the Cu₆₄Zr₃₆ alloy to systematically demonstrate the following two points, all using a single model MG. First, we show in a comprehensive manner the interrelation among common vibrational parameters characterizing the atomic vibrational amplitude and frequency, including the atomic mean square displacement, flexibility volume, participation fraction in the low-frequency vibrational modes and boson peak intensity. Second, we demonstrate that these vibrational parameters fare much better than purely static structural parameters based on local geometrical packing in providing correlation with the propensity for local configurational transitions. These vibrational parameters also share a correlation length similar to that in structural rearrangements induced by external stimuli. This success, however, also poses a challenge, as it remains to be elucidated as to why short-time dynamical (vibrational) behavior at the bottom of the energy basin can be exploited to project the height of the energy barrier for cross-basin activities and in turn the propensity for locally collective atomic rearrangements.

Does mesoscopic elasticity control viscous slowing down in glassforming liquids?

The dramatic slowing down of relaxation dynamics of liquids approaching the glass transition remains a highly debated problem, where the crux of the puzzle resides in the elusive increase in the activation barrier ΔE(T) with decreasing temperature T. A class of theoretical frameworks-known as elastic models-attribute this temperature dependence to the variations of the liquid's macroscopic elasticity, quantified by the high-frequency shear modulus G_∞(T). While elastic models find some support in a number of experimental studies, these models do not take into account the spatial structures, length scales, and heterogeneity associated with structural relaxation in supercooled liquids. Here, we propose and test the possibility that viscous slowing down is controlled by a mesoscopic elastic stiffness κ(T), defined as the characteristic stiffness of response fields to local dipole forces in the liquid's underlying inherent structures. First, we show that κ(T)-which is intimately related to the energy and length scales characterizing quasilocalized, nonphononic excitations in glasses-increases more strongly with decreasing T than the macroscopic inherent structure shear modulus G(T) [the glass counterpart of liquids' G_∞(T)] in several computer liquids. Second, we show that the simple relation ΔE(T) ∝ κ(T) holds remarkably well for some computer liquids, suggesting a direct connection between the liquid's underlying mesoscopic elasticity and enthalpic energy barriers. On the other hand, we show that for other computer liquids, the above relation fails. Finally, we provide strong evidence that what distinguishes computer liquids in which the ΔE(T) ∝ κ(T) relation holds from those in which it does not is that the latter feature highly fragmented/granular potential energy landscapes, where many sub-basins separated by low activation barriers exist. Under such conditions, it appears that the sub-basins do not properly represent the landscape properties relevant for structural relaxation.

Structural and short-time vibrational properties of colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition

We investigate the short-time vibrational properties and structure of two-dimensional, bidisperse, colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition, as a function of interparticle depletion attraction strength. The long-time spatiotemporal dynamics of the samples are measured to be non-monotonic, confirming that the suspensions evolve from repulsive glass to supercooled liquid to attractive glass with increasing depletion attraction. Here, we search for vibrational signatures of the re-entrant behavior in the short-time spatiotemporal dynamics, i.e., dynamics associated with particle motion inside its nearest-neighbor cage. Interestingly, we observe that the anharmonicity of these in-cage vibrations varies non-monotonically with increasing attraction strength, consistent with the non-monotonic long-time structural relaxation dynamics of the re-entrant glass. We also extract effective spring constants between neighboring particles; we find that spring stiffness involving small particles also varies non-monotonically with increasing attraction strength, while stiffness between large particles increases monotonically. Last, from study of depletion-dependent local structure and vibration participation fractions, we gain microscopic insight into the particle-size-dependent contributions to short-time vibrational modes in the glass and supercooled liquid states.

Unified view of avalanche criticality in sheared glasses

Plastic events in sheared glasses are considered an example of so-called avalanches, whose sizes obey a power-law probability distribution with the avalanche critical exponent τ. Although the so-called mean-field depinning (MFD) theory predicts a universal value of this exponent, τ_{MFD}=1.5, such a simplification is now known to connote qualitative disagreement with realistic systems. Numerically and experimentally, different values of τ have been reported depending on the literature. Moreover, in the elastic regime, it has been noted that the critical exponent can be different from that in the steady state, and even criticality itself is a matter of debate. Because these confusingly varying results have been reported under different setups, our knowledge of avalanche criticality in sheared glasses is greatly limited. To gain a unified understanding, in this work, we conduct a comprehensive numerical investigation of avalanches in Lennard-Jones glasses under athermal quasistatic shear. In particular, by excluding the ambiguity and arbitrariness that has crept into the conventional measurement schemes, we achieve high-precision measurement and demonstrate that the exponent τ in the steady state follows the prediction of MFD theory, τ_{MFD}=1.5. Our results also suggest that there are two qualitatively different avalanche events. This binariness leads to the nonuniversal behavior of the avalanche size distribution and is likely to be the cause of the varying values of τ reported thus far. To investigate the dependence of criticality and universality on applied shear, we further study the statistics of avalanches in the elastic regime and the ensemble of the first avalanche event in different samples, which provide information about the unperturbed system. We show that while the unperturbed system is indeed off-critical, criticality gradually develops as shear is applied. The degree of criticality is encoded in the fractal dimension of the avalanches, which starts from zero in the off-critical unperturbed state and saturates in the steady state. Moreover, the critical exponent τ is consistent with the prediction of the MFD τ_{MFD} universally, regardless of the amount of applied shear, once the system becomes critical.

Multiperiodic orbits from interacting soft spots in cyclically sheared amorphous solids

When an amorphous solid is deformed cyclically, it may reach a steady state in which the paths of constituent particles trace out closed loops that repeat in each driving cycle. A remarkable variant has been noticed in simulations where the period of particle motions is a multiple of the period of driving, but the reasons for this behavior have remained unclear. Motivated by mesoscopic features of displacement fields in experiments on jammed solids, we propose and analyze a simple model of interacting soft spots-locations where particles rearrange under stress and that resemble two-level systems with hysteresis. We show that multiperiodic behavior can arise among just three or more soft spots that interact with each other, but in all cases it requires frustrated interactions, illuminating this otherwise elusive type of interaction. We suggest directions for seeking this signature of frustration in experiments and for achieving it in designed systems.