Critical exponents of the yielding transition of amorphous solids

Thermally activated intermittent flow in amorphous solids

Using mean field theory and a mesoscale elastoplastic model, we analyze the steady state shear rheology of thermally activated amorphous solids. At sufficiently high temperature and driving rates, flow is continuous and described by well-established rheological flow laws such as Herschel-Bulkley and logarithmic rate dependence. However, we find that these flow laws change in the regime of intermittent flow, where collective events no longer overlap and serrated flow becomes pronounced. In this regime, we identify a thermal activation stress scale, xa(T,), that wholly captures the effect of driving rate  and temperature T on average flow stress, stress drop (avalanche) size and correlation lengths. Different rheological regimes are summarized in a dynamic phase diagram for the amorphous yielding transition. Theoretical predictions call for a need to re-examine the rheology of very slowly sheared amorphous matter much below the glass transition.

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.

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.

Quasistatic deformation of yield stress materials: Homogeneous or localized?

We analyze a mesoscopic model of a shear stress material with a three-dimensional slab geometry, under an external quasistatic deformation of a simple shear type. Relaxation is introduced in the model as a mechanism by which an unperturbed system achieves progressively mechanically more stable configurations. Although in all cases deformation occurs via localized plastic events (avalanches), we find qualitatively different behavior depending on the degree of relaxation in the model. For no or low relaxation, yielding is homogeneous in the sample, and even the largest avalanches become negligible in size compared with the system size (measured as the thickness of the slab L_{z}) when this is increased. On the contrary, for high relaxation, the deformation localizes in an almost two-dimensional region where all avalanches occur. Scaling analysis of the numerical results indicates that in this case, the linear size of the largest avalanches is comparable with L_{z}, even when this becomes very large. We correlate the two scenarios with a qualitative difference in the flow curve of the system in the two cases, which is monotonous in the first case and velocity weakening in the second case.

Mechanical excitation and marginal triggering during avalanches in sheared amorphous solids

We study plastic strain during individual avalanches in overdamped particle-scale molecular dynamics (MD) and mesoscale elastoplastic models (EPM) for amorphous solids sheared in the athermal quasistatic limit. We show that the spatial correlations in plastic activity exhibit a short length scale that grows as t^{3/4} in MD and ballistically in EPM, which is generated by mechanical excitation of nearby sites not necessarily close to their stability thresholds, and a longer lengthscale that grows diffusively for both models and is associated with remote marginally stable sites. These similarities in spatial correlations explain why simple EPMs accurately capture the size distribution of avalanches observed in MD, though the temporal profiles and dynamical critical exponents are quite different.

Dynamic phase diagram of plastically deformed amorphous solids at finite temperature

The yielding transition that occurs in amorphous solids under athermal quasistatic deformation has been the subject of many theoretical and computational studies. Here, we extend this analysis to include thermal effects at finite shear rate, focusing on how temperature alters avalanches. We derive a nonequilibrium phase diagram capturing how temperature and strain rate effects compete, when avalanches overlap, and whether finite-size effects dominate over temperature effects. The predictions are tested through simulations of an elastoplastic model in two dimensions and in a mean-field approximation. We find a scaling for temperature-dependent softening in the low-strain rate regime when avalanches do not overlap, and a temperature-dependent Herschel-Bulkley exponent in the high-strain rate regime when avalanches do overlap.

Volume-shear coupling in a mesoscopic model of amorphous materials

We present a two-dimensional mesoscopic model of a yield stress material that includes the possibility of local volume fluctuations coupled to shear in such a way that the shear strength of the material decreases as the local density decreases. The model reproduces a number of effects well known in the phenomenology of this kind of material. In particular, we find that the volume of the sample increases as the deformation rate increases; shear bands are no longer oriented at 45^{∘} with respect to the principal axis of the applied stress (as in the absence of volume-shear coupling); and homogeneous deformation becomes unstable at low enough deformation rates if volume-shear coupling is strong enough. We also discuss the effect of this coupling on some out-of-equilibrium configurations, which can be relevant to the study of the shear bands observed in metallic glasses.

Signatures of the spatial extent of plastic events in the yielding transition in amorphous solids

Amorphous solids are yield stress materials that flow when a sufficient load is applied. Their flow consists of periods of elastic loading interrupted by rapid stress drops, or avalanches, coming from microscopic rearrangements known as shear transformations (STs). Here we show that the spatial extent of avalanches in a steadily sheared amorphous solid has a profound effect on the distribution of local residual stresses that in turn determines the stress drop statistics. As reported earlier, the most unstable sites are located in a flat "plateau" region that decreases with system size. While the entrance into the plateau is set by the lower cutoff of the mechanical noise produced by individual STs, the departure from the usually assumed power-law (pseudogap) form of the residual stress distribution comes from far field effects related to spatially extended rearrangements. Interestingly, we observe that the average residual stress of the weakest sites is located in an intermediate power-law regime between the pseudogap and the plateau regimes, whose exponent decreases with system size. Our findings imply a new scaling relation linking the exponents characterizing the avalanche size and residual stress distributions.

Yielding in an Integer Automaton Model for Amorphous Solids under Cyclic Shear

We present results on an automaton model of an amorphous solid under cyclic shear. After a transient, the steady state falls into one of three cases in order of increasing strain amplitude: (i) pure elastic behavior with no plastic activity, (ii) limit cycles where the state recurs after an integer period of strain cycles, and (iii) irreversible plasticity with longtime diffusion. The number of cycles N required for the system to reach a periodic orbit diverges as the amplitude approaches the yielding transition between regimes (ii) and (iii) from below, while the effective diffusivity D of the plastic strain field vanishes on approach from above. Both of these divergences can be described by a power law. We further show that the average period T of the limit cycles increases on approach to yielding.

Properties of the density of shear transformations in driven amorphous solids

The strain load Δγthat triggers consecutive avalanches is a key observable in the slow deformation of amorphous solids. Its temporally averaged value ⟨Δγ⟩ displays a non-trivial system-size dependence that constitutes one of the distinguishing features of the yielding transition. Details of this dependence are not yet fully understood. We address this problem by means of theoretical analysis and simulations of elastoplastic models for amorphous solids. An accurate determination of the size dependence of ⟨Δγ⟩ leads to a precise evaluation of the steady-state distribution of local distances to instabilityx. We find that the usually assumed formP(x) ∼x^θ(withθbeing the so-called pseudo-gap exponent) is not accurate at lowxand that in generalP(x) tends to a system-size-dependentfinitelimit asx→ 0. We work out the consequences of this finite-size dependence standing on exact results for random-walks and disclosing an alternative interpretation of the mechanical noise felt by a reference site. We test our predictions in two- and three-dimensional elastoplastic models, showing the crucial influence of the saturation ofP(x) at smallxon the size dependence of ⟨Δγ⟩ and related scalings.

Residual stress distributions in amorphous solids from atomistic simulations

The distribution of local residual stresses (threshold to instability) that controls the statistical properties of plastic flow in athermal amorphous solids is examined with an atomistic simulation technique. For quiescent configurations, the distribution has a pseudogap (power-law) form with an exponent that agrees well with global yielding statistics. As soon as deformation sets in, the pseudogap region gives way to a system size dependent plateau at small residual stresses that can be understood from the statistics of local residual stress differences between plastic events. Results further suggest that the local yield stress in amorphous solids changes even if the given region does not participate in plastic activity.

Tensorial description of the plasticity of amorphous composites

We use a continuous mesoscopic model to address the yielding properties of plastic composites, formed by a host material and inclusions with different elastic and/or plastic properties. We investigate the flow properties of the composed material under a uniform externally applied deviatoric stress. We show that due to the heterogeneities induced by the inclusions, a scalar modeling in terms of a single deviatoric strain of the same symmetry as the externally applied deformation gives inaccurate results. A realistic modeling must include all possible shear deformations. Implementing this model in a two-dimensional system, we show that the effect of harder inclusions is very weak to relatively high concentrations. For softer inclusions instead, the effect is much stronger; even a small concentration of inclusions affecting the form of the flow curve and the critical stress. We also present the details of a full three-dimensional simulation scheme and obtain the corresponding results for harder and softer inclusions.

Criticality in elastoplastic models of amorphous solids with stress-dependent yielding rates

We analyze the behavior of different elastoplastic models approaching the yielding transition. We propose two kinds of rules for the local yielding events: yielding occurs above the local threshold either at a constant rate or with a rate that increases as the square root of the stress excess. We establish a family of "static" universal critical exponents which do not depend on this dynamic detail of the model rules: in particular, the exponents for the avalanche size distribution P(S) ∼S^-τ_Sf(S/L^d_f) and the exponents describing the density of sites at the verge of yielding, which we find to be of the form P(x) ≃P(0) + x^θ with P(0) ∼L^-a controlling the extremal statistics. On the other hand, we discuss "dynamical" exponents that are sensitive to the local yielding rule. We find that, apart form the dynamical exponent z controlling the duration of avalanches, also the flowcurve's (inverse) Herschel-Bulkley exponent β ([small gamma, Greek, dot above]∼ (σ-σ_c)^β) enters in this category, and is seen to differ in ½ between the two yielding rate cases. We give analytical support to this numerical observation by calculating the exponent variation in the Hébraud-Lequeux model and finding an identical shift. We further discuss an alternative mean-field approximation to yielding only based on the so-called Hurst exponent of the accumulated mechanical noise signal, which gives good predictions for the exponents extracted from simulations of fully spatial models.

Elastic Interfaces on Disordered Substrates: From Mean-Field Depinning to Yielding

We consider a model of an elastic manifold driven on a disordered energy landscape, with generalized long range elasticity. Varying the form of the elastic kernel by progressively allowing for the existence of zero modes, the model interpolates smoothly between mean-field depinning and finite dimensional yielding. We find that the critical exponents of the model change smoothly in this process. Also, we show that in all cases the Herschel-Buckley exponent of the flow curve depends on the analytical form of the microscopic pinning potential. Within the present elastoplastic description, all this suggests that yielding in finite dimensions is a mean-field transition.

Criticality at a Finite Strain Rate in Fluidized Soft Glassy Materials

We study the emergence of critical dynamics in the steady shear rheology of fluidized soft glassy materials. Within a mesoscale elastoplastic model accounting for a shear band instability, we show how additional noise can induce a transition from a phase separated to homogeneous flow, accompanied by critical-like fluctuations of the macroscopic shear rate. Both macroscopic quantities and fluctuations exhibit power law behaviors in the vicinity of this transition, consistent with previous experimental findings on vibrated granular media. Altogether, our results suggest a generic scenario for the emergence of criticality when shear weakening mechanisms compete with a fluidizing noise.