From depinning transition to plastic yielding of amorphous media: A soft-modes perspective

Yielding Is an Absorbing Phase Transition with Vanishing Critical Fluctuations

The yielding transition in athermal complex fluids can be interpreted as an absorbing phase transition between an elastic, absorbing state with high mesoscopic degeneracy and a flowing, active state. We characterize quantitatively this phase transition in an elastoplastic model under fixed applied shear stress, using a finite-size scaling analysis. We find vanishing critical fluctuations of the order parameter (i.e., the shear rate), and relate this property to the convex character of the phase transition (β>1). We locate yielding within a family of models akin to fixed-energy sandpile (FES) models, only with long-range redistribution kernels with zero modes that result from mechanical equilibrium. For redistribution kernels with sufficiently fast decay, this family of models belongs to a short-range universality class distinct from the conserved directed percolation class of usual FES, which is induced by zero modes.

Loss of memory of an elastic line on its way to limit cycles

Oscillatory-driven amorphous materials forget their initial configuration and converge to limit cycles. Here we investigate this memory loss under a nonquasistatic drive in a minimal model system, with quenched disorder and memory encoded in a spatial pattern, where oscillating protocols are formally replaced by a positive-velocity drive. We consider an elastic line driven athermally in a quenched disorder with biperiodic boundary conditions and tunable system size, thus controlling the area swept by the line per cycle as would the oscillation amplitude. The convergence to disorder-dependent limit cycle is strongly coupled to the nature of its velocity dynamics depending on system size. Based on the corresponding phase diagram, we propose a generic scenario for memory formation in disordered systems under finite driving rate.

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.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

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.

Variety of scaling behaviors in nanocrystalline plasticity

We address the question of why larger, high-symmetry crystals are mostly weak, ductile, and statistically subcritical, while smaller crystals with the same symmetry are strong, brittle and supercritical. We link it to another question of why intermittent elasto-plastic deformation of submicron crystals features highly unusual size sensitivity of scaling exponents. We use a minimal integer-valued automaton model of crystal plasticity to show that with growing variance of quenched disorder, which can serve in this case as a proxy for increasing size, submicron crystals undergo a crossover from spin-glass marginality to criticality characterizing the second order brittle-to-ductile (BD) transition. We argue that this crossover is behind the nonuniversality of scaling exponents observed in physical and numerical experiments. The nonuniversality emerges only if the quenched disorder is elastically incompatible, and it disappears if the disorder is compatible.

Rejuvenation and shear banding in model amorphous solids

We measure the local yield stress, at the scale of small atomic regions, in a deeply quenched two-dimensional glass model undergoing shear banding in response to athermal quasistatic deformation. We find that the occurrence of essentially a single plastic event suffices to bring the local yield stress distribution to a well-defined value for all strain orientations, thus essentially erasing the memory of the initial structure. It follows that in a well-relaxed sample, plastic events cause the abrupt (nucleation-like) emergence of a local softness contrast and thus precipitate the formation of a band, which, in its early stages, is measurably softer than the steady-state flow. Moreover, this postevent yield stress ensemble presents a mean value comparable to that of the inherent states of a supercooled liquid around the mode-coupling temperature T_{MCT}. This, we argue, explains that the transition between brittle and ductile yielding in amorphous materials occurs around a comparable parent temperature. Our data also permit to capture quantitatively the contributions of pressure and density changes and demonstrate unambiguously that they are negligible compared with the changes of softness caused by structural rejuvenation.

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.

Avalanches, thresholds, and diffusion in mesoscale amorphous plasticity

We present results on a mesoscale model for amorphous matter in athermal, quasistatic (a-AQS), steady-state shear flow. In particular, we perform a careful analysis of the scaling with the lateral system size L of (i) statistics of individual relaxation events in terms of stress relaxation S, and individual event mean-squared displacement M, and the subsequent load increments Δγ, required to initiate the next event; (ii) static properties of the system encoded by x=σ_{y}-σ, the distance of local stress values from threshold; and (iii) long-time correlations and the emergence of diffusive behavior. For the event statistics, we find that the distribution of S is similar to, but distinct from, the distribution of M. We find a strong correlation between S and M for any particular event, with S∼M^{q} with q≈0.65. The exponent q completely determines the scaling exponents for P(M) given those for P(S). For the distribution of local thresholds, we find P(x) is analytic at x=0, and has a value P(x)|_{x=0}=p_{0} which scales with lateral system length as p_{0}∝L^{-0.6}. The size dependence of the average load increment 〈Δγ〉 appears to be asymptotically controlled by the plateau behavior of P(x) rather than by a subsequent apparent power-law behavior. Extreme value statistics arguments lead thus to a scaling relation between the exponents governing P(x) and those governing P(S). Finally, we study the long-time correlations via single-particle tracer statistics. The value of the diffusion coefficient is completely determined by 〈Δγ〉 and the scaling properties of P(M) (in particular from 〈M〉) rather than directly from P(S) as one might have naively guessed. Our results (i) further define the a-AQS universality class, (ii) clarify the relation between avalanches of stress relaxation and diffusive behavior, and (iii) clarify the relation between local threshold distributions and event statistics.

Compressive Failure as a Critical Transition: Experimental Evidence and Mapping onto the Universality Class of Depinning

Acoustic emission (AE) measurements performed during the compressive loading of concrete samples with three different microstructures (aggregate sizes and porosity) and four sample sizes revealed that failure is preceded by an acceleration of the rate of fracturing events, power law distributions of AE energies and durations near failure, and a divergence of the fracturing correlation length and time towards failure. This argues for an interpretation of compressive failure of disordered materials as a critical transition between an intact and a failed state. The associated critical exponents were found to be independent of sample size and microstructural disorder and close to mean-field depinning values. Although compressive failure differs from classical depinning in several respects, including the nature of the elastic redistribution kernel, an analogy between the two processes allows deriving (finite-) sizing effects on strength that match our extensive data set. This critical interpretation of failure may have also important consequences in terms of natural hazards forecasting, such as volcanic eruptions, landslides, or cliff collapses.

Diffusion in Mesoscopic Lattice Models of Amorphous Plasticity

We present results on tagged particle diffusion in a mesoscale lattice model for sheared amorphous material in athermal quasistatic conditions. We find a short time diffusive regime and a long time diffusive regime whose diffusion coefficients depend on system size in dramatically different ways. At short time, we find that the diffusion coefficient, D, scales roughly linearly with system length, D∼L^{1.05}. This short time behavior is consistent with particle-based simulations. The long-time diffusion coefficient scales like D∼L^{1.6}, close to previous studies which found D∼L^{1.5}. Furthermore, we show that the near-field details of the interaction kernel do not affect the short time behavior but qualitatively and dramatically affect the long time behavior, potentially causing a saturation of the mean-squared displacement at long times. Our finding of a D∼L^{1.05} short time scaling resolves a long standing puzzle about the disagreement between the diffusion coefficient measured in particle-based models and mesoscale lattice models of amorphous plasticity.

Crackling Dynamics in the Mechanical Response of Knitted Fabrics

Crackling noise, which occurs in a wide range of situations, is characterized by discrete events of various sizes, often correlated in the form of avalanches. We report experimental evidence that the mechanical response of a knitted fabric displays such broadly distributed events both in the force signal and in the deformation field, with statistics analogous to that of earthquakes or soft amorphous materials. A knit consists of a regular network of frictional contacts, linked by the elasticity of the yarn. When deformed, the fabric displays spatially extended avalanchelike yielding events resulting from collective interyarn contact slips. We measure the size distribution of these avalanches, at the stitch level from the analysis of nonelastic displacement fields and externally from force fluctuations. The two measurements yield consistent power law distributions reminiscent of those found in other avalanching systems. Our study shows that a knitted fabric is not only a thread-based metamaterial with highly sought after mechanical properties, but also an original, model system, with topologically protected structural order, where an intermittent, scale-invariant response emerges from minimal ingredients, and thus a significant landmark in the study of out-of-equilibrium universality.

Correlation and shear bands in a plastically deformed granular medium

Recent experiments (Le Bouil et al., Phys. Rev. Lett., 2014, 112, 246001) have analyzed the statistics of local deformation in a granular solid undergoing plastic deformation. Experiments report strongly anisotropic correlation between events, with a characteristic angle that was interpreted using elasticity theory and the concept of Eshelby transformations with dilation; interestingly, the shear bands that characterize macroscopic failure occur at an angle that is different from the one observed in microscopic correlations. Here, we interpret this behavior using a mesoscale elastoplastic model of solid flow that incorporates a local Mohr-Coulomb failure criterion. This differs from the interpretation of Le Bouil et al., which is based on purely elastic considerations ignoring the potential role of local friction on deformation patterns. We show that the angle observed in the microscopic correlations can be understood by combining the elastic interactions associated with Eshelby transformation with the local failure criterion. At large strains, we also induce permanent shear bands at an angle that is different from the one observed in the correlation pattern. We interpret this angle as the one that leads to the maximal instability of slip lines.

Local yield stress statistics in model amorphous solids

We develop and extend a method presented by Patinet, Vandembroucq, and Falk [Phys. Rev. Lett. 117, 045501 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.045501] to compute the local yield stresses at the atomic scale in model two-dimensional Lennard-Jones glasses produced via differing quench protocols. This technique allows us to sample the plastic rearrangements in a nonperturbative manner for different loading directions on a well-controlled length scale. Plastic activity upon shearing correlates strongly with the locations of low yield stresses in the quenched states. This correlation is higher in more structurally relaxed systems. The distribution of local yield stresses is also shown to strongly depend on the quench protocol: the more relaxed the glass, the higher the local plastic thresholds. Analysis of the magnitude of local plastic relaxations reveals that stress drops follow exponential distributions, justifying the hypothesis of an average characteristic amplitude often conjectured in mesoscopic or continuum models. The amplitude of the local plastic rearrangements increases on average with the yield stress, regardless of the system preparation. The local yield stress varies with the shear orientation tested and strongly correlates with the plastic rearrangement locations when the system is sheared correspondingly. It is thus argued that plastic rearrangements are the consequence of shear transformation zones encoded in the glass structure that possess weak slip planes along different orientations. Finally, we justify the length scale employed in this work and extract the yield threshold statistics as a function of the size of the probing zones. This method makes it possible to derive physically grounded models of plasticity for amorphous materials by directly revealing the relevant details of the shear transformation zones that mediate this process.

Local and global avalanches in a two-dimensional sheared granular medium

We present the experimental and numerical studies of a two-dimensional sheared amorphous material composed of bidisperse photoelastic disks. We analyze the statistics of avalanches during shear including the local and global fluctuations in energy and changes in particle positions and orientations. We find scale-free distributions for these global and local avalanches denoted by power laws whose cutoffs vary with interparticle friction and packing fraction. Different exponents are found for these power laws depending on the quantity from which variations are extracted. An asymmetry in time of the avalanche shapes is evidenced along with the fact that avalanches are mainly triggered by the shear bands. A simple relation independent of the intensity is found between the number of local avalanches and the global avalanches they form. We also compare these experimental and numerical results for both local and global fluctuations to predictions from mean-field and depinning theories.

Different universality classes at the yielding transition of amorphous systems

We study the yielding transition of a two-dimensional amorphous system under shear by using a mesoscopic elasto-plastic model. The model combines a full (tensorial) description of the elastic interactions in the system and the possibility of structural reaccommodations that are responsible for the plastic behavior. The possible structural reaccommodations are encoded in the form of a "plastic disorder" potential, which is chosen independently at each position of the sample to account for local heterogeneities. We observe that the stress must exceed a critical value σ_{c} in order for the system to yield. In addition, when the system yields a flow curve (relating stress σ and strain rate γ[over ̇]) of the form γ[over ̇]∼(σ-σ_{c})^{β} is obtained. Remarkably, we observe the value of β to depend on some details of the plastic disorder potential. For smooth potentials a value of β≃2.0 is obtained, whereas for potentials obtained as a concatenation of smooth pieces a value β≃1.5 is observed in the simulations. This indicates a dependence of critical behavior on details of the plastic behavior. In addition, by integrating out nonessential, harmonic degrees of freedom, we derive a simplified scalar version of the model that represents a collection of interacting Prandtl-Tomlinson particles. A mean-field treatment of this interaction reproduces the difference of β exponents for the two classes of plastic disorder potentials and provides values of β that compare favorably with those found in the full simulations.

Dynamic phases of active matter systems with quenched disorder

Depinning and nonequilibrium transitions within sliding states in systems driven over quenched disorder arise across a wide spectrum of size scales ranging from atomic friction at the nanoscale, flux motion in type II superconductors at the mesoscale, colloidal motion in disordered media at the microscale, and plate tectonics at geological length scales. Here we show that active matter or self-propelled particles interacting with quenched disorder under an external drive represents a class of system that can also exhibit pinning-depinning phenomena, plastic flow phases, and nonequilibrium sliding transitions that are correlated with distinct morphologies and velocity-force curve signatures. When interactions with the substrate are strong, a homogeneous pinned liquid phase forms that depins plastically into a uniform disordered phase and then dynamically transitions first into a moving stripe coexisting with a pinned liquid and then into a moving phase-separated state at higher drives. We numerically map the resulting dynamical phase diagrams as a function of external drive, substrate interaction strength, and self-propulsion correlation length. These phases can be observed for active matter moving through random disorder. Our results indicate that intrinsically nonequilibrium systems can exhibit additional nonequilibrium transitions when subjected to an external drive.