Avalanche localization and crossover scaling in amorphous plasticity

Criticality in the fracture of silica glass: Insights from molecular mechanics

The universality of avalanches characterizing the inelastic response of disordered materials has the potential to bridge the gap from micro to macroscale. In this study, we explore the statistics and the scaling behavior of avalanches occurring during the fracture process in silica glass using molecular mechanics. We introduce a robust method for capturing and quantifying these avalanches, allowing us to perform rigorous statistical analyses, revealing universal power laws associated with critical phenomena. The influence of an initial crack is explored, observing deviations from mean-field predictions while maintaining the property of criticality. However, the avalanche exponents in the unnotched samples are predicted correctly by the mean-field depinning model. Furthermore, we investigate the strain-dependent probability density function, its cutoff function, and the interrelation between the critical exponents. Finally, we unveil distinct scaling behavior for small and large avalanches of the crack growth, shedding light on the underlying fracture mechanisms in silica glass.

Random Traction Yielding Transition in Epithelial Tissues

We investigate how randomly oriented cell traction forces lead to fluidization in a vertex model of epithelial tissues. We find that the fluidization occurs at a critical value of the traction force magnitude F_{c}. We show that this transition exhibits critical behavior, similar to the yielding transition of sheared amorphous solids. However, we find that it belongs to a different universality class, even though it satisfies the same scaling relations between critical exponents established in the yielding transition of sheared amorphous solids. Our work provides a fluidization mechanism through active force generation that could be relevant in biological tissues.

Identifying microscopic factors that influence ductility in disordered solids

There are empirical strategies for tuning the degree of strain localization in disordered solids, but they are system-specific and no theoretical framework explains their effectiveness or limitations. Here, we study three model disordered solids: a simulated atomic glass, an experimental granular packing, and a simulated polymer glass. We tune each system using a different strategy to exhibit two different degrees of strain localization. In tandem, we construct structuro-elastoplastic (StEP) models, which reduce descriptions of the systems to a few microscopic features that control strain localization, using a machine learning-based descriptor, softness, to represent the stability of the disordered local structure. The models are based on calculated correlations of softness and rearrangements. Without additional parameters, the models exhibit semiquantitative agreement with observed stress-strain curves and softness statistics for all systems studied. Moreover, the StEP models reveal that initial structure, the near-field effect of rearrangements on local structure, and rearrangement size, respectively, are responsible for the changes in ductility observed in the three systems. Thus, StEP models provide microscopic understanding of how strain localization depends on the interplay of structure, plasticity, and elasticity.

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.

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.

Critical yielding rheology: from externally deformed glasses to active systems

We use extensive computer simulations to study the yielding transition under two different loading schemes: standard simple shear dynamics and self-propelled dense active systems. In the active systems, a yielding transition toward an out-of-equilibrium flowing state known as the liquid phase is observed when self-propulsion is increased. The range of self-propulsions in which this pure liquid regime exists appears to vanish upon approaching the so-called 'jamming point' at which the solidity of soft-sphere packings is lost. Such an 'active yielding' transition shares similarities with the generic yielding transition for shear flows. A Herschel-Bulkley law is observed along the liquid regime in both loading scenarios, with a clear difference in the critical scaling exponents between the two, suggesting the existence of different universality classes for the yielding transition under different driving conditions. In addition, we present the direct measurements of growing length and time scales for both driving scenarios. A comparison with theoretical predictions from the recent literature reveals poor agreement with our numerical results.

Unified view of avalanche criticality in sheared glasses

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

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.

Criticality in sheared, disordered solids. I. Rate effects in stress and diffusion

Rate effects in sheared disordered solids are studied using molecular dynamics simulations of binary Lennard-Jones glasses in two and three dimensions. In the quasistatic (QS) regime, systems exhibit critical behavior: the magnitudes of avalanches are power-law distributed with a maximum cutoff that diverges with increasing system size L. With increasing rate, systems move away from the critical yielding point and the average flow stress rises as a power of the strain rate with exponent 1/β, the Herschel-Bulkley exponent. Finite-size scaling collapses of the stress are used to measure β as well as the exponent ν which characterizes the divergence of the correlation length. The stress and kinetic energy per particle experience fluctuations with strain that scale as L^{-d/2}. As the largest avalanche in a system scales as L^{α}, this implies α<d/2. The diffusion rate of particles diverges as a power of decreasing rate before saturating in the QS regime. A scaling theory for the diffusion is derived using the QS avalanche rate distribution and generalized to the finite strain rate regime. This theory is used to collapse curves for different system sizes and confirm β/ν.

Relaxation and vibrational properties in metal alloys and other disordered systems

The relaxation dynamics and the vibrational spectra of amorphous solids, such as metal alloys, have been intensely investigated as well separated topics in the past. The aim of this review is to summarize recent results in both these areas in an attempt to establish, or unveil, deeper connections between the two phenomena of relaxation and vibration. Theoretical progress in the area of slow relaxation dynamics of liquid and glassy systems and in the area of vibrational spectra of glasses and liquids is reviewed. After laying down a generic modelling framework to connect vibration and relaxation, the physics of metal alloys is considered where the emergence of power-law exponents has been identified both in the vibrational density of states (VDOS) as well as in density correlations. Also, theoretical frameworks which connect the VDOS to the relaxation behaviour and mechanical viscoelastic response in metallic glasses are reviewed. The same generic interpretative framework is then applied to the case of molecular glass formers where the emergence of stretched-exponential relaxation in dielectric relaxation can be put in quantitative relation with the VDOS by means of memory-function approaches. Further connections between relaxation and vibration are provided by the study of phonon linewidths in liquids and glasses, where a natural starting point is given by hydrodynamic theories. Finally, an agenda of outstanding issues including the appearance of compressed exponential relaxation in the intermediate scattering function of experimental and simulated systems (metal alloys, colloidal gels, jammed packings) is presented in light of available (or yet to be developed) mathematical models, and compared to non-exponential behaviour measured with macroscopic means such as mechanical spectroscopy/rheology.

On interevent time distributions of avalanche dynamics

Physical systems characterized by stick-slip dynamics often display avalanches. Regardless of the diversity of their microscopic structure, these systems are governed by a power-law distribution of avalanche size and duration. Here we focus on the interevent times between avalanches and show that, unlike their distributions of size and duration, the interevent time distributions are able to distinguish different mechanical states of the system. We use experiments on granular systems and numerical simulations of emulsions to show that systems having the same probability distribution for avalanche size and duration can have different interevent time distributions. Remarkably, these interevent time distributions look similar to those for earthquakes and, if different from an exponential, are indirect evidence of non trivial space-time correlations among avalanches. Our results therefore indicate that interevent time statistics are essential to characterise the dynamics of avalanches.

Elementary plastic events in amorphous silica

Plastic instabilities in amorphous materials are often studied using idealized models of binary mixtures that do not capture accurately molecular interactions and bonding present in real glasses. Here we study atomic-scale plastic instabilities in a three-dimensional molecular dynamics model of silica glass under quasistatic shear. We identify two distinct types of elementary plastic events, one is a standard quasilocalized atomic rearrangement while the second is a bond-breaking event that is absent in simplified models of fragile glass formers. Our results show that both plastic events can be predicted by a drop of the lowest nonzero eigenvalue of the Hessian matrix that vanishes at a critical strain. Remarkably, we find very high correlation between the associated eigenvectors and the nonaffine displacement fields accompanying the bond-breaking event, predicting the locus of structural failure. Both eigenvectors and nonaffine displacement fields display an Eshelby-like quadrupolar structure for both failure modes, rearrangement, and bond breaking. Our results thus clarify the nature of atomic-scale plastic instabilities in silica glasses, providing useful information for the development of mesoscale models of amorphous plasticity.

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.

Avalanche statistics during coarsening dynamics

We study the coarsening dynamics of a two-dimensional system via numerical simulations. The system under consideration is a biphasic system consisting of domains of a dispersed phase closely packed together in a continuous phase and separated by thin interfaces. Such a system is elastic and typically out of equilibrium. The equilibrium state is attained via the coarsening dynamics, wherein the dispersed phase slowly diffuses through the interfaces, causing the domains to change in size and eventually rearrange abruptly. The effect of rearrangements is propagated throughout the system via the intrinsic elastic interactions and may cause rearrangements elsewhere, resulting in intermittent bursts of activity and avalanche behaviour. Here we aim at quantitatively characterizing the corresponding avalanche statistics (i.e. size, duration, and inter-avalanche time). Despite the coarsening dynamics is triggered by an internal driving mechanism, we find quantitative indications that such avalanche statistics displays scaling-laws very similar to those observed in the response of disordered materials to external loads.

Seismicity in sheared granular matter

We report on experiments investigating the dynamics of a slider that is pulled by a spring across a granular medium consisting of a vertical layer of photoelastic disks. The motion proceeds through a sequence of discrete events, analogous to seismic shocks, in which elastic energy stored in the spring is rapidly released. We measure the statistics of several properties of the individual events: the energy loss in the spring, the duration of the movement, and the temporal profile of the slider motion. We also study certain conditional probabilities and the statistics of mainshock-aftershock sequences. At low driving rates, we observe crackling with Omori-Utsu, Båth, and waiting time laws similar to those observed in seismic dynamics. At higher driving rates, where the sequence of events shows strong periodicity, we observe scaling laws and asymmetrical event shapes that are clearly distinguishable from those in the crackling regime.

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.

Avalanche Behavior in Creep Failure of Disordered Materials

We present a mesoscale elastoplastic model of creep in disordered materials, which considers temperature-dependent stochastic activation of localized deformation events that are coupled by internal stresses, leading to collective avalanche dynamics. We generalize this stochastic plasticity model by introducing damage in terms of a local strength that decreases, on statistical average, with increasing local plastic strain. The model captures failure in terms of strain localization in a catastrophic shear band concomitant with a finite-time singularity of the creep rate. The statistics of avalanches in the run-up to failure is characterized by a decreasing avalanche exponent τ that, at failure, approaches the value τ=1.5 typical of a critical branching process. The average avalanche rate exhibits an inverse Omori law as a function of time to failure. The distribution of interavalanche times turns out to be consistent with the epidemic-type aftershock sequences (ETAS) model of earthquake statistics.

Orientation of plastic rearrangements in two-dimensional model glasses under shear

The plastic deformation of amorphous solids is mediated by localized shear transformations involving small groups of particles rearranging irreversibly in an elastic background. We introduce and compare three different computational methods to extract the size and orientation of these shear transformations in simulations of a two-dimensional athermal model glass under simple shear. We find that the shear angles are broadly distributed around the macroscopic shear direction, with a more or less Gaussian distribution with a standard deviation of around 20^{∘}. The distributions of sizes and orientations of shear transformations display no substantial sensitivity to the shear rate. These results can notably be used to refine the description of rearrangements in elastoplastic models.

Soft modes and strain redistribution in continuous models of amorphous plasticity: the Eshelby paradigm, and beyond?

The deformation of disordered solids relies on swift and localised rearrangements of particles. The inspection of soft vibrational modes can help predict the locations of these rearrangements, while the strain that they actually redistribute mediates collective effects. Here, we study soft modes and strain redistribution in a two-dimensional continuous mesoscopic model based on a Ginzburg-Landau free energy for perfect solids, supplemented with a plastic disorder potential that accounts for shear softening and rearrangements. Regardless of the disorder strength, our numerical simulations show soft modes that are always sharply peaked at the softest point of the material (unlike what happens for the depinning of an elastic interface). Contrary to widespread views, the deformation halo around this peak does not always have a quadrupolar (Eshelby-like) shape. Instead, for finite and narrowly-distributed disorder, it looks like a fracture, with a strain field that concentrates along some easy directions. These findings are rationalised with analytical calculations in the case where the plastic disorder is confined to a point-like 'impurity'. In this case, we unveil a continuous family of elastic propagators, which are identical for the soft modes and for the equilibrium configurations. This family interpolates between the standard quadrupolar propagator and the fracture-like one as the anisotropy of the elastic medium is increased. Therefore, we expect to see a fracture-like propagator when extended regions on the brink of failure have already softened along the shear direction and thus rendered the material anisotropic, but not failed yet. We speculate that this might be the case in carefully aged glasses just before macroscopic failure.

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.

Universal features of amorphous plasticity

Plastic yielding of amorphous solids occurs by power-law distributed deformation avalanches whose universality is still debated. Experiments and molecular dynamics simulations are hampered by limited statistical samples, and although existing stochastic models give precise exponents, they require strong assumptions about fixed deformation directions, at odds with the statistical isotropy of amorphous materials. Here, we introduce a fully tensorial, stochastic mesoscale model for amorphous plasticity that links the statistical physics of plastic yielding to engineering mechanics. It captures the complex shear patterning observed for a wide variety of deformation modes, as well as the avalanche dynamics of plastic flow. Avalanches are described by universal size exponents and scaling functions, avalanche shapes, and local stability distributions, independent of system dimensionality, boundary and loading conditions, and stress state. Our predictions consistently differ from those of mean-field depinning models, providing evidence that plastic yielding is a distinct type of critical phenomenon.

The universality class for plastic yield in amorphous materials remains controversial. Here authors present a tensorial mesoscale model that captures both complex shear patterns and avalanche scaling behaviour, which differs from mean-field models and suggests a distinct type of critical phenomenon.

Scaling of slip avalanches in sheared amorphous materials based on large-scale atomistic simulations

Atomistic simulations of binary amorphous systems with over 4 million atoms are performed. Systems of two interatomic potentials of the Lennard-Jones type, LJ12-6 and LJ9-6, are simulated. The athermal quasistatic shearing protocol is adopted, where the shear strain is applied in a stepwise fashion with each step followed by energy minimization. For each avalanche event, the shear stress drop (Δσ), the hydrostatic pressure drop (Δσ_{h}), and the potential energy drop (ΔE) are computed. It is found that, with the avalanche size increasing, the three become proportional to each other asymptotically. The probability distributions of avalanche sizes are obtained and values of scaling exponents fitted. In particular, the distributions follow a power law, P(ΔU)∼ΔU^{-τ}, where ΔU is a measure of avalanche sizes defined based on shear stress drops. The exponent τ is 1.25±0.1 for the LJ12-6 systems, and 1.15±0.1 for the LJ9-6 systems. The value of τ for the LJ12-6 systems is consistent with that from an earlier atomistic simulation study by Robbins et al. [Phys. Rev. Lett. 109, 105703 (2012)]PRLTAO0031-900710.1103/PhysRevLett.109.105703, but the fitted values of other scaling exponents differ, which may be because the shearing protocol used here differs from that in their study.

Direct Observation of Percolation in the Yielding Transition of Colloidal Glasses

When strained beyond the linear regime, soft colloidal glasses yield to steady-state plastic flow in a way that is similar to the deformation of conventional amorphous solids. Because of the much larger size of the colloidal particles with respect to the atoms comprising an amorphous solid, colloidal glasses allow us to obtain microscopic insight into the nature of the yielding transition, as we illustrate here combining experiments, atomistic simulations, and mesoscopic modeling. Our results unanimously show growing clusters of nonaffine deformation percolating at yielding. In agreement with percolation theory, the spanning cluster is fractal with a fractal dimension d_{f}≃2, and the correlation length diverges upon approaching the critical yield strain. These results indicate that percolation of highly nonaffine particles is the hallmark of the yielding transition in disordered glassy systems.

Stress and temperature dependence of the avalanche dynamics during creep deformation of metallic glasses

The understanding of the mesoscopic origin of plasticity in metallic glasses remains still an open issue. At the microscopic level, Shear Transformation Zones (STZ), composed by dozens of atoms, have been identified as the basic unit of the deformation process. Macroscopically, metallic glasses perform either homogeneous or inhomogeneous flow depending on the experimental conditions. However, the emergence of macroscopic behavior resulting from STZ interactions is still an open issue and is of great interest. In the current work we present an approach to analyze the different interaction mechanisms of STZ's by studying the statistics of the avalanches produced by a metallic glass during tensile creep deformation. We identified a crossover between different regimes of avalanches, and we analyzed the dependence of such crossover on the experimental conditions, namely stress and temperature. We interpret such crossover as a transition from 3D random STZ activity to localized 2D nano-shear bands. The experimental time at which the crossover takes place seems to depend on the overall strain and strain rate in the sample.

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

A mesoscopic model of amorphous plasticity is discussed in the context of depinning models. After embedding in a d+1-dimensional space, where the accumulated plastic strain lives along the additional dimension, the gradual plastic deformation of amorphous media can be regarded as the motion of an elastic manifold in a disordered landscape. While the associated depinning transition leads to scaling properties, the quadrupolar Eshelby interactions at play in amorphous plasticity induce specific additional features like shear-banding and weak ergodicity breakdown. The latters are shown to be controlled by the existence of soft modes of the elastic interaction, the consequence of which is discussed in the context of depinning.

Shearing a glass and the role of pinning delay in models of interface depinning

When a disordered solid is sheared, yielding is followed by the onset of intermittent response that is characterized by slip in local regions usually labeled shear-transformation zones. Such intermittent response resembles the behavior of earthquakes or contact depinning, where a well-defined landscape of pinning disorder prohibits the deformation of an elastic medium. Nevertheless, a disordered solid is evidently different in that pinning barriers of particles are due to neighbors that are also subject to motion. Microscopic yielding leads to destruction of the local microstructure and local heating. It is natural to assume that locally a liquid emerges for a finite timescale before cooling down to a transformed configuration. For including this characteristic transient in glass depinning models, we propose a general mechanism that involves a "pinning delay" time T(pd), during which each region that slipped evolves as a fluid. The new timescale can be as small as a single avalanche time step. This is a local, effective, and dynamical in nature mechanism that may be thought as dynamical softening. We demonstrate that the inclusion of this mechanism causes a drift of the critical exponents toward higher values for the slip sizes τ, until a transition to permanent shear-banding behavior happens causing almost oscillatory, stick-slip response. Moreover, it leads to a proliferation of large events that are highly inhomogeneous and resemble sharp slip band formation.

Driving Rate Dependence of Avalanche Statistics and Shapes at the Yielding Transition

We study stress time series caused by plastic avalanches in athermally sheared disordered materials. Using particle-based simulations and a mesoscopic elastoplastic model, we analyze system size and shear-rate dependence of the stress-drop duration and size distributions together with their average temporal shape. We find critical exponents different from mean-field predictions, and a clear asymmetry for individual avalanches. We probe scaling relations for the rate dependency of the dynamics and we report a crossover towards mean-field results for strong driving.

Finite-size effects in a model for plasticity of amorphous composites

We discuss the plastic behavior of an amorphous matrix reinforced by hard particles. A mesoscopic depinning-like model accounting for Eshelby elastic interactions is implemented. Only the effect of a plastic disorder is considered. Numerical results show a complex size dependence of the effective flow stress of the amorphous composite. In particular, the departure from the mixing law shows opposite trends associated to the competing effects of the matrix and the reinforcing particles, respectively. The reinforcing mechanisms and their effects on localization are discussed. Plastic strain is shown to gradually concentrate on the weakest band of the system. This correlation of the plastic behavior with the material structure is used to design a simple analytical model. The latter nicely captures reinforcement size effects in (logN/N)(1/2), where N is the linear size of the system, observed numerically. Predictions of the effective flow stress accounting for further logarithmic corrections show a very good agreement with numerical results.

Cooperativity flows and shear-bandings: a statistical field theory approach

Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin and A. Ajdari, Phys. Rev. Lett., 2009, 103, 036001, we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity) are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. The compacton coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material.

Probing relevant ingredients in mean-field approaches for the athermal rheology of yield stress materials

Although the notion of mechanical noise is expected to play a key role in the non-linear rheology of athermally sheared amorphous systems, its characterization has so far remained elusive. Here, we show using molecular dynamic simulations that in spite of the presence of strong spatio-temporal correlations in the system, the local stress exhibits normal diffusion under the effect of the mechanical noise in the finite driving regime. The diffusion constant appears to be proportional to the mean plastic activity. Our data suggests that the corresponding proportionality constant is density independent, and can be directly related to the specific form of the rheological flow curve, pointing the way to a generic way of modeling mechanical noise in mean-field equations.

Moving line model and avalanche statistics of Bingham fluid flow in porous media

In this article, we propose a simple model to understand the critical behavior of path opening during flow of a yield stress fluid in porous media as numerically observed by Chevalier and Talon (2015). This model can be mapped to the problem of a contact line moving in an heterogeneous field. Close to the critical point, this line presents an avalanche dynamic where the front advances by a succession of waiting time and large burst events. These burst events are then related to the non-flowing (i.e. unyielded) areas. Remarkably, the statistics of these areas reproduce the same properties as in the direct numerical simulations. Furthermore, even if our exponents seem to be close to the mean field universal exponents, we report an unusual bump in the distribution which depends on the disorder. Finally, we identify a scaling invariance of the cluster spatial shape that is well fit, to first order, by a self-affine parabola.

Quenched pinning and collective dislocation dynamics

Several experiments show that crystalline solids deform in a bursty and intermittent fashion. Power-law distributed strain bursts in compression experiments of micron-sized samples, and acoustic emission energies from larger-scale specimens, are the key signatures of the underlying critical-like collective dislocation dynamics - a phenomenon that has also been seen in discrete dislocation dynamics (DDD) simulations. Here we show, by performing large-scale two-dimensional DDD simulations, that the character of the dislocation avalanche dynamics changes upon addition of sufficiently strong randomly distributed quenched pinning centres, present e.g. in many alloys as immobile solute atoms. For intermediate pinning strength, our results adhere to the scaling picture of depinning transitions, in contrast to pure systems where dislocation jamming dominates the avalanche dynamics. Still stronger disorder quenches the critical behaviour entirely.

Atomistic mechanisms of intermittent plasticity in metals: dislocation avalanches and defect cluster pinning

Intermittent plastic deformation in crystals with power-law behaviors has been reported in previous experimental studies. The power-law behavior is reminiscent of self-organized criticality, and mesoscopic models have been proposed that describe this behavior in crystals. In this paper, we show that intermittent plasticity in metals under tensile deformation can be observed in molecular dynamics models, using embedded atom method potentials for Ni, Cu, and Al. Power-law behaviors of stress drop and waiting time of plastic deformation events are observed. It is shown that power-law behavior is due to dislocation avalanche motions in Cu and Ni. A different mechanism of dislocation pinning is found in Al. These different stress relaxation mechanisms give different power-law exponents. We propose a probabilistic model to describe the novel dislocation motion in Al and analytically deduce the power-law behavior.

Morphological study of elastic-plastic-brittle transitions in disordered media

We use a spring lattice model with springs following a bilinear elastoplastic-brittle constitutive behavior with spatial disorder in the yield and failure thresholds to study patterns of plasticity and damage evolution. The elastic-perfectly plastic transition is observed to follow percolation scaling with the correlation length critical exponent ν≈1.59, implying the universality class corresponding to the long-range correlated percolation. A quantitative analysis of the plastic strain accumulation reveals a dipolar anisotropy (for antiplane loading) which vanishes with increasing hardening modulus. A parametric study with hardening modulus and ductility controlled through the spring level constitutive response demonstrates a wide spectrum of behaviors with varying degree of coupling between plasticity and damage evolution.

Scaling description of the yielding transition in soft amorphous solids at zero temperature

Yield stress materials flow if a sufficiently large shear stress is applied. Although such materials are ubiquitous and relevant for industry, there is no accepted microscopic description of how they yield, even in the simplest situations in which temperature is negligible and in which flow inhomogeneities such as shear bands or fractures are absent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the Herschel-Bulkley exponent characterizing the singularity of the flow curve near the yield stress Σc, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x beyond which they yield. We argue that the critical exponents of the yielding transition may be expressed in terms of three independent exponents, θ, df, and z, characterizing, respectively, the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elasto-plastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions in both two and three dimensions.

Universal and non-universal features in coarse-grained models of flow in disordered solids

We study the two-dimensional (2D) shear flow of amorphous solids within variants of an elastoplastic model, paying particular attention to spatial correlations and time fluctuations of, e.g., local stresses. The model is based on the local alternation between an elastic regime and plastic events during which the local stress is redistributed. The importance of a fully tensorial description of the stress and of the inclusion of (coarse-grained) convection in the model is investigated; scalar and tensorial models yield similar results, while convection enhances fluctuations and breaks the spurious symmetry between the flow and velocity gradient directions, for instance when shear localisation is observed. Besides, correlation lengths measured with diverse protocols are discussed. One class of such correlation lengths simply scale with the spacing between homogeneously distributed, simultaneous plastic events. This leads to a scaling of the correlation length with the shear rate as γ̇(-1/2) in 2D in the athermal regime, regardless of the details of the model. The radius of the cooperative disk, defined as the near-field region in which plastic events induce a stress redistribution that is not amenable to a mean-field treatment, notably follows this scaling. On the other hand, the cooperative volume measured from the four-point stress susceptibility and its dependence on the system size and the shear rate are model-dependent.