Following the evolution of hard sphere glasses in infinite dimensions under external perturbations: compression and shear strain

Ductile and brittle yielding of athermal amorphous solids: A mean-field paradigm beyond the random-field Ising model

Amorphous solids can yield in either a ductile or brittle manner under strain: plastic deformation can set in gradually, or abruptly through a macroscopic stress drop. Developing a unified theory describing both ductile and brittle yielding constitutes a fundamental challenge of nonequilibrium statistical physics. Recently, it has been proposed that, in the absence of thermal effects, the nature of the yielding transition is controlled by physics akin to that of the quasistatically driven random field Ising model (RFIM), which has served as the paradigm for understanding the effect of quenched disorder in slowly driven systems with short-ranged interactions. However, this theoretical picture neglects both the dynamics of, and the elasticity-induced long-ranged interactions between, the mesoscopic material constituents. Here, we address these two aspects and provide a unified theory building on the Hébraud-Lequeux elastoplastic description. The first aspect is crucial to understanding the competition between the imposed deformation rate and the finite timescale of plastic rearrangements: We provide a dynamical description of the macroscopic stress drop, as well as predictions for the shifting of the brittle yield strain and the scaling of the peak susceptibility with inverse shear rate. The second is essential to capture properly the behavior in the limit of quasistatic driving, where avalanches of plasticity diverge with system size at any value of the strain. We fully characterise the avalanche behavior, which is radically different to that of the RFIM. In the quasistatic, infinite-size limit, we find that both models have mean-field Landau exponents, obscuring the effect of the interactions. We show, however, that the latter profoundly affect the behavior of finite systems approaching the spinodal-like brittle yield point, where we recover qualitatively the finite-size trends found in particle simulations. The interactions also modify the nature of the random critical point separating ductile and brittle yielding, where we predict critical behavior on top of the marginality present at any value of the strain. We finally discuss how all our predictions can be directly tested against particle simulations and eventually experiments, and make first steps in this direction.

The yielding of granular matter is marginally stable and critical

Amorphous materials undergo a transition from liquid-like to solid-like states through processes like rapid quenching or densification. Under external loads, they exhibit yielding, with minimal structural changes compared to crystals. However, these universal characteristics are rarely explored comprehensively in a single granular experiment due to the added complexity of inherent friction. The discernible differences between static configurations before and after yielding are largely unaddressed, and a comprehensive examination from both statistical physics and mechanical perspectives is lacking. To address these gaps, we conducted experiments using photoelastic disks, simultaneously tracking particles and measuring forces. Our findings reveal that the yielding transition demonstrates critical behavior from a statistical physics standpoint and marginal stability from a mechanical perspective, akin to the isotropic jamming transition. This criticality differs significantly from spinodal criticality in frictionless amorphous solids, highlighting unique characteristics of granular yielding. Furthermore, our analysis confirms the marginal stability of granular yielding by assessing the contact number and evaluating the balance between weak forces and small gaps. These factors serve as structural indicators for configurations before and after yielding. Our results not only contribute to advancing our understanding of the fundamental physics of granular materials but also bear significant implications for practical applications in various fields.

Brittle and ductile yielding in soft materials

Many soft materials yield under mechanical loading, but how this transition from solid-like behavior to liquid-like behavior occurs can vary significantly. Understanding the physics of yielding is of great interest for the behavior of biological, environmental, and industrial materials, including those used as inks in additive manufacturing and muds and soils. For some materials, the yielding transition is gradual, while others yield abruptly. We refer to these behaviors as being ductile and brittle. The key rheological signatures of brittle yielding include a stress overshoot in steady-shear-startup tests and a steep increase in the loss modulus during oscillatory amplitude sweeps. In this work, we show how this spectrum of yielding behaviors may be accounted for in a continuum model for yield stress materials by introducing a parameter we call the brittility factor. Physically, an increased brittility decreases the contribution of recoverable deformation to plastic deformation, which impacts the rate at which yielding occurs. The model predictions are successfully compared to results of different rheological protocols from a number of real yield stress fluids with different microstructures, indicating the general applicability of the phenomenon of brittility. Our study shows that the brittility of soft materials plays a critical role in determining the rate of the yielding transition and provides a simple tool for understanding its effects under various loading conditions.

Dynamical yield criterion for granular matter from first principles

We investigate, using a recently developed model of liquid state theory describing the rheology of dense granular flows, how a yield stress appears in granular matter at the yielding transition. Our model allows us to predict an analytical equation of the corresponding dynamical yield surface, which is compared to the usual models of solid fracture. In particular, this yield surface interpolates between the typical failure behaviors of soft and hard materials. This work also underlines the central role played by the effective friction coefficient at the yielding transition.

Glasslike caging with random planes

The richness of the mean-field solution of simple glasses leaves many of its features challenging to interpret. A minimal model that illuminates glass physics in the same way that the random energy model clarifies spin glass behavior would therefore be beneficial. Here we propose such a real-space model that is amenable to infinite-dimensional d→∞ analysis and is exactly solvable in finite d in some regimes. By joining analysis with numerical simulations, we uncover geometrical signatures of the dynamical and jamming transitions and obtain insight into the origin of activated processes. Translating these findings into the context of standard glass formers further reveals the role played by nonconvexity in the emergence of Gardner and jamming physics.

Experimental observations of fractal landscape dynamics in a dense emulsion

Many soft and biological materials display so-called 'soft glassy' dynamics; their constituents undergo anomalous random motions and complex cooperative rearrangements. A recent simulation model of one soft glassy material, a coarsening foam, suggested that the random motions of its bubbles are due to the system configuration moving over a fractal energy landscape in high-dimensional space. Here we show that the salient geometrical features of such high-dimensional fractal landscapes can be explored and reliably quantified, using empirical trajectory data from many degrees of freedom, in a model-free manner. For a mayonnaise-like dense emulsion, analysis of the observed trajectories of oil droplets quantitatively reproduces the high-dimensional fractal geometry of the configuration path and its associated local energy minima generated using a computational model. That geometry in turn drives the droplets' complex random motion observed in real space. Our results indicate that experimental studies can elucidate whether the similar dynamics in different soft and biological materials may also be due to fractal landscape dynamics.

Dynamic Gardner cross-over in a simple glass

The criticality of the jamming transition responsible for amorphous solidification has been theoretically linked to the marginal stability of a thermodynamic Gardner phase. While the critical exponents of jamming appear independent of the preparation history, the pertinence of Gardner physics far from equilibrium is an open question. To fill this gap, we numerically study the nonequilibrium dynamics of hard disks compressed toward the jamming transition using a broad variety of protocols. We show that dynamic signatures of Gardner physics can be disentangled from the aging relaxation dynamics. We thus define a generic dynamic Gardner cross-over regardless of the history. Our results show that the jamming transition is always accessed by exploring increasingly complex landscape, resulting in anomalous microscopic relaxation dynamics that remains to be understood theoretically.

How to count in hierarchical landscapes: A full solution to mean-field complexity

We derive the general solution for counting the stationary points of mean-field complex landscapes. It incorporates Parisi's solution for the ground state, as it should. Using this solution, we count the stationary points of two models: one with multistep replica symmetry breaking and one with full replica symmetry breaking.

Finite-Disorder Critical Point in the Yielding Transition of Elastoplastic Models

Upon loading, amorphous solids can exhibit brittle yielding, with the abrupt formation of macroscopic shear bands leading to fracture, or ductile yielding, with a multitude of plastic events leading to homogeneous flow. It has been recently proposed, and subsequently questioned, that the two regimes are separated by a sharp critical point, as a function of some control parameter characterizing the intrinsic disorder strength and the degree of stability of the solid. In order to resolve this issue, we have performed extensive numerical simulations of athermally driven elastoplastic models with long-range and anisotropic realistic interaction kernels in two and three dimensions. Our results provide clear evidence for a finite-disorder critical point separating brittle and ductile yielding, and we provide an estimate of the critical exponents in 2D and 3D.

Nonlinear elasticity, yielding, and entropy in amorphous solids

The holographic duality has proven successful in linking seemingly unrelated problems in physics. Recently, intriguing correspondences between the physics of soft matter and gravity are emerging, including strong similarities between the rheology of amorphous solids, effective field theories for elasticity, and the physics of black holes. However, direct comparisons between theoretical predictions and experimental/simulation observations remain limited. Here, we study the effects of nonlinear elasticity on the mechanical and thermodynamic properties of amorphous materials responding to shear, using effective field and gravitational theories. The predicted correlations among the nonlinear elastic exponent, the yielding strain/stress, and the entropy change due to shear are supported qualitatively by simulations of granular matter models. Our approach opens a path toward understanding the complex mechanical responses of amorphous solids, such as mixed effects of shear softening and shear hardening, and offers the possibility to study the rheology of solid states and black holes in a unified framework.

Determining the nonequilibrium criticality of a Gardner transition via a hybrid study of molecular simulations and machine learning

Apparent critical phenomena, typically indicated by growing correlation lengths and dynamical slowing down, are ubiquitous in nonequilibrium systems such as supercooled liquids, amorphous solids, active matter, and spin glasses. It is often challenging to determine if such observations are related to a true second-order phase transition as in the equilibrium case or simply a crossover and even more so to measure the associated critical exponents. Here we show that the simulation results of a hard-sphere glass in three dimensions are consistent with the recent theoretical prediction of a Gardner transition, a continuous nonequilibrium phase transition. Using a hybrid molecular simulation-machine learning approach, we obtain scaling laws for both finite-size and aging effects and determine the critical exponents that traditional methods fail to estimate. Our study provides an approach that is useful to understand the nature of glass transitions and can be generalized to analyze other nonequilibrium phase transitions.

Yielding in an amorphous solid subject to constant stress at finite temperatures

Understanding the nature of the yield transition is a long-standing problem in the physics of amorphous solids. Here we use molecular dynamics simulations to study the response of amorphous solids to constant stresses at finite temperatures. We compare amorphous solids that are prepared using fast and slow quenches and show that for thermal systems, the steady-state velocity exhibits a continuous transition from very slow creep to a finite strain rate as a function of the stress. This behavior is observed for both well-annealed and poorly annealed systems. However, the transient dynamics is different in the latter and involves overcoming an energy barrier. Due to the different simulation protocol, the strain rate as a function of stress and temperature follows a scaling relation that is different from the ones that are shown for systems where the strain is controlled. Collapsing the data using this scaling relation allows us to calculate critical exponents for the dynamics close to yield, including an exponent for thermal rounding. We also demonstrate that strain slips due to avalanche events above yield follow standard scaling relations and we extract critical exponents that are comparable to the ones obtained in previous studies that performed simulations of both molecular dynamics and elastoplastic models using strain-rate control.

Universality of plastic instability and mechanical yield in metallic glasses

The generic response of a wide range of amorphous solids is the average increase of stress upon external loading until the yielding transition point, after which elasto-plastic steady state sets in. The stress-strain response comprises of a series of elastic branches interspersed with plastic drops. The ubiquitousness of these phenomena indicates universality, independent of material property, but the literature predominantly deals with specific materials. In pursuit of generality among different amorphous systems, we undertake a careful investigation in the mechanical response of metallic glasses using computer simulation. By comparing our results of multi-body metallic glass potentials to those obtained from pairwise Lennard-Jones glasses, we show that the mechanism of plastic instabilities is universal and independent of the details of the underlying potential. We also investigate the yielding transition in terms of the overlap parameterQ₁₂, which has been successfully used Lennard-Jones glasses. The yielding is unambiguously identified as a first-order phase transition. These observations conform the nature of plastic instabilities and mechanical yield as universal and independent of microscopic interactions.

A direct link between active matter and sheared granular systems

The similarity in mechanical properties of dense active matter and sheared amorphous solids has been noted in recent years without a rigorous examination of the underlying mechanism. We develop a mean-field model that predicts that their critical behavior-as measured by their avalanche statistics-should be equivalent in infinite dimensions up to a rescaling factor that depends on the correlation length of the applied field. We test these predictions in two dimensions using a numerical protocol, termed "athermal quasistatic random displacement," and find that these mean-field predictions are surprisingly accurate in low dimensions. We identify a general class of perturbations that smoothly interpolates between the uncorrelated localized forces that occur in the high-persistence limit of dense active matter and system-spanning correlated displacements that occur under applied shear. These results suggest a universal framework for predicting flow, deformation, and failure in active and sheared disordered materials.

The role of annealing in determining the yielding behavior of glasses under cyclic shear deformation

Yielding behavior in amorphous solids has been investigated in computer simulations using uniform and cyclic shear deformation. Recent results characterize yielding as a discontinuous transition, with the degree of annealing of glasses being a significant parameter. Under uniform shear, discontinuous changes in stresses at yielding occur in the high annealing regime, separated from the poor annealing regime in which yielding is gradual. In cyclic shear simulations, relatively poorly annealed glasses become progressively better annealed as the yielding point is approached, with a relatively modest but clear discontinuous change at yielding. To understand better the role of annealing on yielding characteristics, we perform athermal quasistatic cyclic shear simulations of glasses prepared with a wide range of annealing in two qualitatively different systems-a model of silica (a network glass) and an atomic binary mixture glass. Two strikingly different regimes of behavior emerge. Energies of poorly annealed samples evolve toward a unique threshold energy as the strain amplitude increases, before yielding takes place. Well-annealed samples, in contrast, show no significant energy change with strain amplitude until they yield, accompanied by discontinuous energy changes that increase with the degree of annealing. Significantly, the threshold energy for both systems corresponds to dynamical cross-over temperatures associated with changes in the character of the energy landscape sampled by glass-forming liquids.

A jamming plane of sphere packings

The concept of jamming has attracted great research interest due to its broad relevance in soft-matter, such as liquids, glasses, colloids, foams, and granular materials, and its deep connection to sphere packing and optimization problems. Here, we show that the domain of amorphous jammed states of frictionless spheres can be significantly extended, from the well-known jamming-point at a fixed density, to a jamming-plane that spans the density and shear strain axes. We explore the jamming-plane, via athermal and thermal simulations of compression and shear jamming, with initial equilibrium configurations prepared by an efficient swap algorithm. The jamming-plane can be divided into reversible-jamming and irreversible-jamming regimes, based on the reversibility of the route from the initial configuration to jamming. Our results suggest that the irreversible-jamming behavior reflects an escape from the metastable glass basin to which the initial configuration belongs to or the absence of such basins. All jammed states, either compression- or shear-jammed, are isostatic and exhibit jamming criticality of the same universality class. However, the anisotropy of contact networks nontrivially depends on the jamming density and strain. Among all state points on the jamming-plane, the jamming-point is a unique one with the minimum jamming density and the maximum randomness. For crystalline packings, the jamming-plane shrinks into a single shear jamming-line that is independent of initial configurations. Our study paves the way for solving the long-standing random close-packing problem and provides a more complete framework to understand jamming.

Dilatancy, shear jamming, and a generalized jamming phase diagram of frictionless sphere packings

Granular packings display the remarkable phenomenon of dilatancy, wherein their volume increases upon shear deformation. Conventional wisdom and previous results suggest that dilatancy, also being the related phenomenon of shear-induced jamming, requires frictional interactions. Here, we show that the occurrence of isotropic jamming densities φ_j above the minimal density (or the J-point density) φ_J leads both to the emergence of shear-induced jamming and dilatancy in frictionless packings. Under constant pressure shear, the system evolves into a steady-state at sufficiently large strains, whose density only depends on the pressure and is insensitive to the initial jamming density φ_j. In the limit of vanishing pressure, the steady-state exhibits critical behavior at φ_J. While packings with different φ_j values display equivalent scaling properties under compression, they exhibit striking differences in rheological behaviour under shear. The yield stress under constant volume shear increases discontinuously with density when φ_j > φ_J, contrary to the continuous behaviour in generic packings that jam at φ_J. Our results thus lead to a more coherent, generalised picture of jamming in frictionless packings, which also have important implications on how dilatancy is understood in the context of frictional granular matter.

Mechanical disorder of sticky-sphere glasses. II. Thermomechanical inannealability

Many structural glasses feature static and dynamic mechanical properties that can depend strongly on glass formation history. The degree of universality of this history dependence and what it is possibly affected by are largely unexplored. Here we show that the variability of elastic properties of simple computer glasses under thermal annealing depends strongly on the strength of attractive interactions between the glasses' constituent particles-referred to here as glass "stickiness." We find that in stickier glasses the stiffening of the shear modulus with thermal annealing is strongly suppressed, while the thermal-annealing-induced softening of the bulk modulus is enhanced. Our key finding is that the characteristic frequency and density per frequency of soft quasilocalized modes becomes effectively invariant to annealing in very sticky glasses; the latter are therefore deemed "thermomechanically inannealable." The implications of our findings and future research directions are discussed.

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.

Ductile and Brittle Yielding in Thermal and Athermal Amorphous Materials

We study theoretically the yielding of sheared amorphous materials as a function of increasing levels of initial sample annealing prior to shear, in three widely used constitutive models and three widely studied annealing protocols. In thermal systems we find a gradual progression, with increasing annealing, from smoothly "ductile" yielding, in which the sample remains homogeneous, to abruptly "brittle" yielding, in which it becomes strongly shear banded. This progression arises from an increase with annealing in the size of an overshoot in the underlying stress-strain curve for homogeneous shear, which causes a shear banding instability that becomes more severe with increasing annealing. Ductile and brittle yielding thereby emerge as two limiting cases of a continuum of yielding transitions, from gradual to catastrophic. In contrast, athermal systems with a stress overshoot always show brittle yielding at low shear rates, however small the overshoot.

Jamming Energy Landscape is Hierarchical and Ultrametric

The free energy landscape of mean-field marginal glasses is ultrametric. We demonstrate that this feature persists in finite three-dimensional systems that are out of equilibrium by finding sets of minima, which are nearby in configuration space. By calculating the distance between these nearby minima, we produce a small region of the distance metric. This metric exhibits a clear hierarchical structure and shows the signature of an ultrametric space. That such a hierarchy exists for the jamming energy landscape provides direct evidence for the existence of a marginal phase along the zero temperature jamming line.

Effect of Pinning on the Yielding Transition of Amorphous Solids

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

Mean-field stability map of hard-sphere glasses

The response of amorphous solids to an applied shear deformation is an important problem, both in fundamental and in applied research. To tackle this problem, we focus on a system of hard spheres in infinite dimensions as a solvable model for colloidal systems and granular matter. The system is prepared above the dynamical glass transition density, and we discuss the phase diagram of the resulting glass under compression, decompression, and shear strain, expanding on previous results [Urbani and Zamponi, Phys. Rev. Lett. 118, 038001 (2017)PRLTAO0031-900710.1103/PhysRevLett.118.038001]. We show that the solid region is bounded by a "shear jamming" line, at which the solid reaches close packing, and a "shear yielding" line, at which the solid undergoes a spinodal instability towards a liquid flowing phase. Furthermore, we characterize the evolution of these lines upon varying the glass preparation density. This paper aims to provide a general overview on yielding and jamming phenomena in hard-sphere systems by a systematic exploration of the phase diagram.

Critical Jammed Phase of the Linear Perceptron

Criticality in statistical physics naturally emerges at isolated points in the phase diagram. Jamming of spheres is not an exception: varying density, it is the critical point that separates the unjammed phase where spheres do not overlap and the jammed phase where they cannot be arranged without overlaps. The same remains true in more general constraint satisfaction problems with continuous variables where jamming coincides with the (protocol dependent) satisfiability transition point. In this work we show that by carefully choosing the cost function to be minimized, the region of criticality extends to occupy a whole region of the jammed phase. As a working example, we consider the spherical perceptron with a linear cost function in the unsatisfiable jammed phase and we perform numerical simulations which show critical power laws emerging in the configurations obtained minimizing the linear cost function. We develop a scaling theory to compute the emerging critical exponents.

Gardner physics in amorphous solids and beyond

One of the most remarkable predictions to emerge out of the exact infinite-dimensional solution of the glass problem is the Gardner transition. Although this transition was first theoretically proposed a generation ago for certain mean-field spin glass models, its materials relevance was only realized when a systematic effort to relate glass formation and jamming was undertaken. A number of nontrivial physical signatures associated with the Gardner transition have since been considered in various areas, from models of structural glasses to constraint satisfaction problems. This perspective surveys these recent advances and discusses the novel research opportunities that arise from them.

Rejuvenation and Memory Effects in a Structural Glass

We show numerically that a three-dimensional model for structural glass displays aging, rejuvenation, and memory effects when subjected to a temperature cycle. These effects indicate that the free energy landscape of structural glasses may possess the complex hierarchical structure that characterizes materials such as spin and polymer glasses. We use the theoretical concept of marginal stability to interpret our results, and explain in which physical conditions a complex aging dynamics can emerge in dense supercooled liquids, paving the way for future experimental studies of complex aging dynamics in colloidal and granular glasses.

Morphology of renormalization-group flow for the de Almeida-Thouless-Gardner universality class

A replica-symmetry-breaking phase transition is predicted in a host of disordered media. The criticality of the transition has, however, long been questioned below its upper critical dimension, six, due to the absence of a critical fixed point in the renormalization-group flows at one-loop order. A recent two-loop analysis revealed a possible strong-coupling fixed point, but given the uncontrolled nature of perturbative analysis in the strong-coupling regime, debate persists. Here we examine the nature of the transition as a function of spatial dimension and show that the strong-coupling fixed point can go through a Hopf bifurcation, resulting in a critical limit cycle and a concomitant discrete scale invariance. We further investigate a different renormalization scheme and argue that the basin of attraction of the strong-coupling fixed point (or limit cycle) may stay finite for all dimensions.

Marginally stable phases in mean-field structural glasses

A novel form of amorphous matter characterized by marginal stability was recently discovered in the mean-field theory of structural glasses. Using this approach, we provide complete phase diagrams delimiting the location of the marginally stable glass phase for a large variety of pair interactions and physical conditions, extensively exploring physical regimes relevant to granular matter, foams, emulsions, hard and soft colloids, and molecular glasses. We find that all types of glasses may become marginally stable, but the extent of the marginally stable phase highly depends on the preparation protocol. Our results suggest that marginal phases should be observable for colloidal and non-Brownian particles near jamming and for poorly annealed glasses. For well-annealed glasses, two distinct marginal phases are predicted. Our study unifies previous results on marginal stability in mean-field models and will be useful to guide numerical simulations and experiments aimed at detecting marginal stability in finite-dimensional amorphous materials.

The physics of jamming for granular materials: a review

Granular materials consist of macroscopic grains, interacting via contact forces, and unaffected by thermal fluctuations. They are one of a class systems that undergo jamming, i.e. a transition between fluid-like and disordered solid-like states. Roughly twenty years ago, proposals by Cates et al for the shear response of colloidal systems and by Liu and Nagel, for a universal jamming diagram in a parameter space of packing fraction, ϕ, shear stress, τ, and temperature, T raised key questions. Contemporaneously, experiments by Howell et al and numerical simulations by Radjai et al and by Luding et al helped provide a starting point to explore key insights into jamming for dry, cohesionless, granular materials. A recent experimental observation by Bi et al is that frictional granular materials have a a re-entrant region in their jamming diagram. In a range of ϕ, applying shear strain, γ, from an initially force/stress free state leads to fragile (in the sense of Cates et al), then anisotropic shear jammed states. Shear jamming at fixed ϕ is presumably conjugate to Reynolds dilatancy, involving dilation under shear against deformable boundaries. Numerical studies by Radjai and Roux showed that Reynolds dilatancy does not occur for frictionless systems. Recent numerical studies by several groups show that shear jamming occurs for finite, but not infinite, systems of frictionless grains. Shear jamming does not lead to known ordering in position space, but Sarkar et al showed that ordering occurs in a space of force tiles. Experimental studies seeking to understand random loose and random close packings (rlp and rcp) and dating back to Bernal have probed granular packings and their response to shear and intruder motion. These studies suggest that rlp's are anisotropic and shear-jammed-like, whereas rcp's are likely isotropically jammed states. Jammed states are inherently static, but the jamming diagram may provide a context for understanding rheology, i.e. dynamic shear in a variety of systems that include granular materials and suspensions.

A stability-reversibility map unifies elasticity, plasticity, yielding, and jamming in hard sphere glasses

Amorphous solids, such as glasses, have complex responses to deformations, with substantial consequences in material design and applications. In this respect, two intertwined aspects are important: stability and reversibility. It is crucial to understand, on the one hand, how a glass may become unstable due to increased plasticity under shear deformations, and, on the other hand, to what extent the response is reversible, meaning how much a system is able to recover the original configuration once the perturbation is released. Here, we focus on assemblies of hard spheres as the simplest model of amorphous solids such as colloidal glasses and granular matter. We prepare glass states quenched from equilibrium supercooled liquid states, which are obtained by using the swap Monte Carlo algorithm and correspond to a wide range of structural relaxation time scales. We exhaustively map out their stability and reversibility under volume and shear strains using extensive numerical simulations. The region on the volume-shear strain phase diagram where the original glass state remains solid is bounded by the shear yielding and the shear jamming lines that meet at a yielding-jamming crossover point. This solid phase can be further divided into two subphases: the stable glass phase, where the system deforms purely elastically and is totally reversible, and the marginal glass phase, where it experiences stochastic plastic deformations at mesoscopic scales and is partially irreversible. The details of the stability-reversibility map depend strongly on the quality of annealing of the glass. This study provides a unified framework for understanding elasticity, plasticity, yielding, and jamming in amorphous solids.

Microscopic Theory of Two-Step Yielding in Attractive Colloids

Attractive colloids display two distinct amorphous solid phases: the attractive glass due to particle bonding and the repulsive glass due to the hard-core repulsion. By means of a microscopic mean field approach, we analyze their response to a quasistatic shear strain. We find that the presence of two distinct interaction length scales may result in a sharp two-step yielding process, which can be associated with a hysteretic stress response or with a reversible but nonmonotonic stress-strain curve. We derive a generic phase diagram characterized by two distinct yielding lines, an inverse yielding and a critical point separating the hysteretic and reversible regimes. Our results should be applicable to a large class of glassy materials characterized by two distinct interaction length scales.

Spin-glass-like aging in colloidal and granular glasses

Motivated by the mean field prediction of a Gardner phase transition between a "normal glass" and a "marginally stable glass", we investigate the off-equilibrium dynamics of three-dimensional polydisperse hard spheres, used as a model for colloidal or granular glasses. Deep inside the glass phase, we find that a sharp crossover pressure PG separates two distinct dynamical regimes. For pressure P < PG, the glass behaves as a normal solid, displaying fast dynamics that quickly equilibrate within the glass free energy basin. For P > PG, instead, the dynamics become strongly anomalous, displaying very large equilibration timescales, aging, and a constantly increasing dynamical susceptibility. The crossover at PG is strongly reminiscent of the one observed in three-dimensional spin-glasses in an external field, suggesting that the two systems could be in the same universality class, consistent with theoretical expectations.

Random critical point separates brittle and ductile yielding transitions in amorphous materials

We combine an analytically solvable mean-field elasto-plastic model with molecular dynamics simulations of a generic glass former to demonstrate that, depending on their preparation protocol, amorphous materials can yield in two qualitatively distinct ways. We show that well-annealed systems yield in a discontinuous brittle way, as metallic and molecular glasses do. Yielding corresponds in this case to a first-order nonequilibrium phase transition. As the degree of annealing decreases, the first-order character becomes weaker and the transition terminates in a second-order critical point in the universality class of an Ising model in a random field. For even more poorly annealed systems, yielding becomes a smooth crossover, representative of the ductile rheological behavior generically observed in foams, emulsions, and colloidal glasses. Our results show that the variety of yielding behaviors found in amorphous materials does not necessarily result from the diversity of particle interactions or microscopic dynamics but is instead unified by carefully considering the role of the initial stability of the system.

Aging, Jamming, and the Limits of Stability of Amorphous Solids

Apart from not having crystallized, supercooled liquids can be considered as being properly equilibrated and thus can be described by a few thermodynamic control variables. In contrast, glasses and other amorphous solids can be arbitrarily far away from equilibrium and require a description of the history of the conditions under which they formed. In this paper we describe how the locality of interactions intrinsic to finite-dimensional systems affects the stability of amorphous solids far off equilibrium. Our analysis encompasses both structural glasses formed by cooling and colloidal assemblies formed by compression. A diagram outlining regions of marginal stability can be adduced which bears some resemblance to the quasi-equilibrium replica meanfield theory phase diagram of hard sphere glasses in high dimensions but is distinct from that construct in that the diagram describes not true phase transitions but kinetic transitions that depend on the preparation protocol. The diagram exhibits two distinct sectors. One sector corresponds to amorphous states with relatively open structures, the other to high density, more closely packed ones. The former transform rapidly owing to there being motions with no free energy barriers; these motions are string-like locally. In the dense region, amorphous systems age via compact activated reconfigurations. The two regimes correspond, in equilibrium, to the collisional or uniform liquid and the so-called landscape regime, respectively. These are separated by a spinodal line of dynamical crossovers. Owing to the rigidity of the surrounding matrix in the landscape, high-density part of the diagram, a sufficiently rapid pressure quench adds compressive energy which also leads to an instability toward string-like motions with near vanishing barriers. Conversely, a dilute collection of rigid particles, such as a colloidal suspension leads, when compressed, to a spatially heterogeneous structure with percolated mechanically stable regions. This jamming corresponds to the onset of activation when the spinodal line is traversed from the low density side. We argue that a stable glass made of sufficiently rigid particles can also be viewed as exhibiting sporadic and localized buckling instabilities that result in local jammed structures. The lines of instability we discuss resemble the Gardner transition of meanfield systems but, in contrast, do not result in true criticality owing to being short-circuited by activated events. The locally marginally stable modes of motion in amorphous solids correspond to secondary relaxation processes in structural glasses. Their relevance to the low temperature anomalies in glasses is also discussed.

Isolating long-wavelength fluctuation from structural relaxation in two-dimensional glass: cage-relative displacement

It has recently been revealed that long-wavelength fluctuation exists in two-dimensional (2D) glassy systems, having the same origin as that given by the Mermin-Wagner theorem for 2D crystalline solids. In this paper, we discuss how to characterise quantitatively the long-wavelength fluctuation in a molecular dynamics simulation of a lightly supercooled liquid. We employ the cage-relative mean-square displacement (MSD), defined on relative displacement to its cage, to quantitatively separate the long-wavelength fluctuation from the original MSD. For increasing system size the amplitude of acoustic long wavelength fluctuations not only increases but shifts to later times causing a crossover with structural relaxation of caging particles. We further analyse the dynamic correlation length using the cage-relative quantities. It grows as the structural relaxation becomes slower with decreasing temperature, uncovering an overestimation by the four-point correlation function due to the long-wavelength fluctuation. These findings motivate the usage of cage-relative MSD as a starting point for analysis of 2D glassy dynamics.

Absence of Marginal Stability in a Structural Glass

Marginally stable solids have peculiar physical properties that were first analyzed in the context of the jamming transition. We theoretically investigate the existence of marginal stability in a prototypical model for structural glass formers, combining analytical calculations in infinite dimensions to computer simulations in three dimensions. While mean-field theory predicts the existence of a Gardner phase transition towards a marginally stable glass phase at low temperatures, simulations show no hint of diverging time scales or length scales, but reveal instead the presence of sparse localized defects. Our results suggest that the Gardner transition is deeply affected by finite dimensional fluctuations, and raise issues about the relevance of marginal stability in structural glasses far away from jamming.

Mechanical failure in amorphous solids: Scale-free spinodal criticality

The mechanical failure of amorphous media is a ubiquitous phenomenon from material engineering to geology. It has been noticed for a long time that the phenomenon is "scale-free," indicating some type of criticality. In spite of attempts to invoke "Self-Organized Criticality," the physical origin of this criticality, and also its universal nature, being quite insensitive to the nature of microscopic interactions, remained elusive. Recently we proposed that the precise nature of this critical behavior is manifested by a spinodal point of a thermodynamic phase transition. Demonstrating this requires the introduction of an "order parameter" that is suitable for distinguishing between disordered amorphous systems. At the spinodal point there exists a divergent correlation length which is associated with the system-spanning instabilities (known also as shear bands) which are typical to the mechanical yield. The theory, the order parameter used and the correlation functions which exhibit the divergent correlation length are universal in nature and can be applied to any amorphous solid that undergoes mechanical yield. The phenomenon is seen at its sharpest in athermal systems, as is explained below; in this paper we extend the discussion also to thermal systems, showing that at sufficiently high temperatures the spinodal phenomenon is destroyed by thermal fluctuations.

Shear bands as manifestation of a criticality in yielding amorphous solids

Amorphous solids increase their stress as a function of an applied strain until a mechanical yield point whereupon the stress cannot increase anymore, afterward exhibiting a steady state with a constant mean stress. In stress-controlled experiments, the system simply breaks when pushed beyond this mean stress. The ubiquity of this phenomenon over a huge variety of amorphous solids calls for a generic theory that is free of microscopic details. Here, we offer such a theory: The mechanical yield is a thermodynamic phase transition, where yield occurs as a spinodal phenomenon. At the spinodal point, there exists a divergent correlation length that is associated with the system-spanning instabilities (also known as shear bands), which are typical to the mechanical yield. The theory, the order parameter used, and the correlation functions that exhibit the divergent correlation length are universal in nature and can be applied to any amorphous solids that undergo mechanical yield.

Nontrivial Critical Fixed Point for Replica-Symmetry-Breaking Transitions

The transformation of the free-energy landscape from smooth to hierarchical is one of the richest features of mean-field disordered systems. A well-studied example is the de Almeida-Thouless transition for spin glasses in a magnetic field, and a similar phenomenon-the Gardner transition-has recently been predicted for structural glasses. The existence of these replica-symmetry-breaking phase transitions has, however, long been questioned below their upper critical dimension, d_{u}=6. Here, we obtain evidence for the existence of these transitions in d<d_{u} using a two-loop calculation. Because the critical fixed point is found in the strong-coupling regime, we corroborate the result by resumming the perturbative series with inputs from a three-loop calculation and an analysis of its large-order behavior. Our study offers a resolution of the long-lasting controversy surrounding phase transitions in finite-dimensional disordered systems.

Exploring the complex free-energy landscape of the simplest glass by rheology

For amorphous solids, it has been intensely debated whether the traditional view on solids, in terms of the ground state and harmonic low energy excitations on top of it, such as phonons, is still valid. Recent theoretical developments of amorphous solids revealed the possibility of unexpectedly complex free-energy landscapes where the simple harmonic picture breaks down. Here we demonstrate that standard rheological techniques can be used as powerful tools to examine nontrivial consequences of such complex free-energy landscapes. By extensive numerical simulations on a hard sphere glass under quasistatic shear at finite temperatures, we show that above the so-called Gardner transition density, the elasticity breaks down, the stress relaxation exhibits slow, and ageing dynamics and the apparent shear modulus becomes protocol-dependent. Being designed to be reproducible in laboratories, our approach may trigger explorations of the complex free-energy landscapes of a large variety of amorphous materials.

Breakdown of nonlinear elasticity in stress-controlled thermal amorphous solids

In recent work it was clarified that amorphous solids under strain control do not possess nonlinear elastic theory in the sense that the shear modulus exists but nonlinear moduli exhibit sample-to-sample fluctuations that grow without bound with the system size. More relevant, however, for experiments are the conditions of stress control. In the present Rapid Communication we show that also under stress control the shear modulus exists, but higher-order moduli show unbounded sample-to-sample fluctuation. The unavoidable consequence is that the characterization of stress-strain curves in experiments should be done with a stress-dependent shear modulus rather than with nonlinear expansions.

Shear Yielding and Shear Jamming of Dense Hard Sphere Glasses

We investigate the response of dense hard sphere glasses to a shear strain in a wide range of pressures ranging from the glass transition to the infinite-pressure jamming point. The phase diagram in the density-shear strain plane is calculated analytically using the mean-field infinite-dimensional solution. We find that just above the glass transition, the glass generically yields at a finite shear strain. The yielding transition in the mean-field picture is a spinodal point in presence of disorder. At higher densities, instead, we find that the glass generically jams at a finite shear strain: the jamming transition prevents yielding. The shear yielding and shear jamming lines merge in a critical point, close to which the system yields at extremely large shear stress. Around this point, highly nontrivial yielding dynamics, characterized by system-spanning disordered fractures, is expected.

Understanding soft glassy materials using an energy landscape approach

Many seemingly different soft materials-such as soap foams, mayonnaise, toothpaste and living cells-display strikingly similar viscoelastic behaviour. A fundamental physical understanding of such soft glassy rheology and how it can manifest in such diverse materials, however, remains unknown. Here, by using a model soap foam consisting of compressible spherical bubbles, whose sizes slowly evolve and whose collective motion is simply dictated by energy minimization, we study the foam's dynamics as it corresponds to downhill motion on an energy landscape function spanning a high-dimensional configuration space. We find that these downhill paths, when viewed in this configuration space, are, surprisingly, fractal. The complex behaviour of our model, including power-law rheology and non-diffusive bubble motion and avalanches, stems directly from the fractal dimension and energy function of these paths. Our results suggest that ubiquitous soft glassy rheology may be a consequence of emergent fractal geometry in the energy landscapes of many complex fluids.

Breakdown of nonlinear elasticity in amorphous solids at finite temperatures

It is known [H. G. E. Hentschel et al., Phys. Rev. E 83, 061101 (2011)PLEEE81539-375510.1103/PhysRevE.83.061101] that amorphous solids at zero temperature do not possess a nonlinear elasticity theory: besides the shear modulus, which exists, none of the higher order coefficients exist in the thermodynamic limit. Here we show that the same phenomenon persists up to temperatures comparable to that of the glass transition. The zero-temperature mechanism due to the prevalence of dangerous plastic modes of the Hessian matrix is replaced by anomalous stress fluctuations that lead to the divergence of the variances of the higher order elastic coefficients. The conclusion is that in amorphous solids elasticity can never be decoupled from plasticity: the nonlinear response is very substantially plastic.

Failure of the Volume Function in Granular Statistical Mechanics and an Alternative Formulation

We first show that the currently accepted statistical mechanics for granular matter is flawed. The reason is that it is based on the volume function, which depends only on a minute fraction of all the structural degrees of freedom and is unaffected by most of the configurational microstates. Consequently, the commonly used partition function underestimates the entropy severely. We then propose a new formulation, replacing the volume function with a connectivity function that depends on all the structural degrees of freedom and accounts correctly for the entire entropy. We discuss the advantages of the new formalism and derive explicit results for two- and three-dimensional systems. We test the formalism by calculating the entropy of an experimental two-dimensional system, as a function of system size, and showing that it is an extensive variable.

Mechanical Yield in Amorphous Solids: A First-Order Phase Transition

Amorphous solids yield at a critical value of the strain (in strain-controlled experiments); for larger strains, the average stress can no longer increase-the system displays an elastoplastic steady state. A long-standing riddle in the materials community is what the difference is between the microscopic states of the material before and after yield. Explanations in the literature are material specific, but the universality of the phenomenon begs a universal answer. We argue here that there is no fundamental difference in the states of matter before and after yield, but the yield is a bona fide first-order phase transition between a highly restricted set of possible configurations residing in a small region of phase space to a vastly rich set of configurations which include many marginally stable ones. To show this, we employ an order parameter of universal applicability, independent of the microscopic interactions, that is successful in quantifying the transition in an unambiguous manner.

Quantitative approximation schemes for glasses

By means of a systematic expansion around the infinite-dimensional solution, we obtain an approximation scheme to compute properties of glasses in low dimensions. The resulting equations take as input the thermodynamic and structural properties of the equilibrium liquid, and from this they allow one to compute properties of the glass. They are therefore similar in spirit to the Mode Coupling approximation scheme. Our scheme becomes exact, by construction, in dimension d→∞, and it can be improved systematically by adding more terms in the expansion.

The structural origin of the hard-sphere glass transition in granular packing

Glass transition is accompanied by a rapid growth of the structural relaxation time and a concomitant decrease of configurational entropy. It remains unclear whether the transition has a thermodynamic origin, and whether the dynamic arrest is associated with the growth of a certain static order. Using granular packing as a model hard-sphere glass, we show the glass transition as a thermodynamic phase transition with a ‘hidden' polytetrahedral order. This polytetrahedral order is spatially correlated with the slow dynamics. It is geometrically frustrated and has a peculiar fractal dimension. Additionally, as the packing fraction increases, its growth follows an entropy-driven nucleation process, similar to that of the random first-order transition theory. Our study essentially identifies a long-sought-after structural glass order in hard-sphere glasses.

Glass transition shows dramatic dynamic slowdown, but its origin remains unclear. Here, Xia et al. observe in granular systems the rapid growth of a geometrically frustrated polytetrahedral order with packing fraction, which is spatially correlated with the slow dynamics.

Numerical detection of the Gardner transition in a mean-field glass former

Recent theoretical advances predict the existence, deep into the glass phase, of a novel phase transition, the so-called Gardner transition. This transition is associated with the emergence of a complex free energy landscape composed of many marginally stable sub-basins within a glass metabasin. In this study, we explore several methods to detect numerically the Gardner transition in a simple structural glass former, the infinite-range Mari-Kurchan model. The transition point is robustly located from three independent approaches: (i) the divergence of the characteristic relaxation time, (ii) the divergence of the caging susceptibility, and (iii) the abnormal tail in the probability distribution function of cage order parameters. We show that the numerical results are fully consistent with the theoretical expectation. The methods we propose may also be generalized to more realistic numerical models as well as to experimental systems.

Calorimetric glass transition in a mean-field theory approach

The study of the properties of glass-forming liquids is difficult for many reasons. Analytic solutions of mean-field models are usually available only for systems embedded in a space with an unphysically high number of spatial dimensions; on the experimental and numerical side, the study of the properties of metastable glassy states requires thermalizing the system in the supercooled liquid phase, where the thermalization time may be extremely large. We consider here a hard-sphere mean-field model that is solvable in any number of spatial dimensions; moreover, we easily obtain thermalized configurations even in the glass phase. We study the 3D version of this model and we perform Monte Carlo simulations that mimic heating and cooling experiments performed on ultrastable glasses. The numerical findings are in good agreement with the analytical results and qualitatively capture the features of ultrastable glasses observed in experiments.