Universal features of amorphous plasticity

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.

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.

Amorphous alloys surpass E/10 strength limit at extreme strain rates

Theoretical predictions of the ideal strength of materials range from E/30 to E/10 (E is Young's modulus). However, despite intense interest over the last decade, the value of the ideal strength achievable through experiments for metals remains a mystery. This study showcases the remarkable spall strength of Cu50Zr50 amorphous alloy that exceeds the E/10 limit at strain rates greater than 107 s-1 through laser-induced shock experiments. The material exhibits a spall strength of 11.5 GPa, approximately E/6 or 1/13 of its P-wave modulus, which sets a record for the elastic limit of metals. Electron microscopy and large-scale molecular dynamics simulations reveal that the primary failure mechanism at extreme strain rates is void nucleation and growth, rather than shear-banding. The rate dependence of material strength is explained by a void kinetic model controlled by surface energy. These findings help advance our understanding on the mechanical behavior of amorphous alloys under extreme strain rates.

Elasticity, Facilitation, and Dynamic Heterogeneity in Glass-Forming Liquids

We study the role of elasticity-induced facilitation on the dynamics of glass-forming liquids by a coarse-grained two-dimensional model in which local relaxation events, taking place by thermal activation, can trigger new relaxations by long-range elastically mediated interactions. By simulations and an analytical theory, we show that the model reproduces the main salient facts associated with dynamic heterogeneity and offers a mechanism to explain the emergence of dynamical correlations at the glass transition. We also discuss how it can be generalized and combined with current theories.

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.

Residual stress in athermal soft disordered solids: insights from microscopic and mesoscale models

In soft amorphous materials, shear cessation after large shear deformation leads to configurations having residual shear stress. The origin of these states and the distribution of the local shear stresses within the material is not well understood, despite its importance for the change in material properties and consequent applications. In this work, we use molecular dynamics simulations of a model dense non-Brownian soft amorphous material to probe the non-trivial relaxation process towards a residual stress state. We find that, similar to thermal glasses, an increase in shear rate prior to the shear cessation leads to lower residual stress states. We rationalise our findings using a mesoscopic elasto-plastic description that explicitly includes a long range elastic response to local shear transformations. We find that after flow cessation the initial stress relaxation indeed depends on the pre-sheared stress state, but the final residual stress is majorly determined by newly activated plastic events occurring during the relaxation process, a scenario consistent with the phenomenology of avalanche dynamics in the low shear rate limit of steadily sheared amorphous solids. Our simplified coarse grained description not only allows capturing the phenomenology of residual stress states but also rationalising the altered material properties that are probed using small and large deformation protocols applied to the relaxed material.

Microscopic dynamics underlying the stress relaxation of arrested soft materials

Arrested soft materials such as gels and glasses exhibit a slow stress relaxation with a broad distribution of relaxation times in response to linear mechanical perturbations. Although this macroscopic stress relaxation is an essential feature in the application of arrested systems as structural materials, consumer products, foods, and biological materials, the microscopic origins of this relaxation remain poorly understood. Here, we elucidate the microscopic dynamics underlying the stress relaxation of such arrested soft materials under both quiescent and mechanically perturbed conditions through X-ray photon correlation spectroscopy. By studying the dynamics of a model associative gel system that undergoes dynamical arrest in the absence of aging effects, we show that the mean stress relaxation time measured from linear rheometry is directly correlated to the quiescent superdiffusive dynamics of the microscopic clusters, which are governed by a buildup of internal stresses during arrest. We also show that perturbing the system via small mechanical deformations can result in large intermittent fluctuations in the form of avalanches, which give rise to a broad non-Gaussian spectrum of relaxation modes at short times that is observed in stress relaxation measurements. These findings suggest that the linear viscoelastic stress relaxation in arrested soft materials may be governed by nonlinear phenomena involving an interplay of internal stress relaxations and perturbation-induced intermittent avalanches.

Avalanches, Clusters, and Structural Change in Cyclically Sheared Silica Glass

We investigate avalanches and clusters associated with plastic rearrangements and the nature of structural change in the prototypical strong glass, silica, computationally. We perform a detailed analysis of avalanches, and of spatially disconnected clusters that constitute them, for a wide range of system sizes. Although qualitative aspects of yielding in silica are similar to other glasses, the statistics of clusters exhibits significant differences, which we associate with differences in local structure. Across the yielding transition, anomalous structural change and densification, associated with a suppression of tetrahedral order, is observed to accompany strain localization.

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.

Quasistatic kinetic avalanches and self-organized criticality in deviatorically loaded granular media

The behavior of granular media under quasistatic loading has recently been shown to attain a stable evolution state corresponding to a manifold in the space of micromechanical variables. This state is characterized by sudden transitions between metastable jammed states, involving the partial micromechanical rearrangement of the granular medium. Using numerical simulations of two-dimensional granular media under quasistatic biaxial compression, we show that the dynamics in the stable evolution state is characterized by scale-free avalanches well before the macromechanical stationary flow regime traditionally linked to a self-organized critical state. This, together with the nonuniqueness and the nonmonotony of macroscopic deformation curves, suggests that the statistical avalanche properties and the susceptibilities of the system cannot be reduced to a function of the macromechanical state. The associated scaling exponents are nonuniversal and depend on the interactions between particles. For stiffer particles (or samples at low confining pressure) we find distributions of avalanche properties compatible with the predictions of mean-field theory. The scaling exponents decrease below the mean-field values for softer interactions between particles. These lower exponents are consistent with observations for amorphous solids at their critical point. We specifically discuss the relationship between microscopic and macroscopic variables, including the relation between the external stress drop and the internal potential energy released during kinetic avalanches.

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.

Elastoplasticity Mediates Dynamical Heterogeneity Below the Mode Coupling Temperature

As liquids approach the glass transition temperature, dynamical heterogeneity emerges as a crucial universal feature of their behavior. Dynamic facilitation, where local motion triggers further motion nearby, plays a major role in this phenomenon. Here we show that long-ranged, elastically mediated facilitation appears below the mode coupling temperature, adding to the short-range component present at all temperatures. Our results suggest deep connections between the supercooled liquid and glass states, and pave the way for a deeper understanding of dynamical heterogeneity in glassy systems.

Effects of external noise on threshold-induced correlations in ferromagnetic systems

In the present paper we investigate the impact of the external noise and detection threshold level on the simulation data for the systems that evolve through metastable states. As a representative model of such systems we chose the nonequilibrium athermal random-field Ising model with two types of the external noise, uniform white noise and Gaussian white noise with various different standard deviations, imposed on the original response signal obtained in model simulations. We applied a wide range of detection threshold levels in analysis of the signal and show how these quantities affect the values of exponent γ_{S/T} (describing the scaling of the average avalanche size with duration), the shift of waiting time between the avalanches, and finally the collapses of the waiting time distributions. The results are obtained via extensive numerical simulations on the equilateral three-dimensional cubic lattices of various sizes and disorders.

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.

Microstructural behavior of magnetorheological elastomer undergoing durability evaluation by stress relaxation

The widespread use of magnetorheological elastomer (MRE) materials in various applications has yet to be limited due to the fact that there are substantial deficiencies in current experimental and theoretical research on its microstructural durability behavior. In this study, MRE composed of silicon rubber (SR) and 70 wt% of micron-sized carbonyl iron particles (CIP) was prepared and subjected to stress relaxation evaluation by torsional shear load. The microstructure and particle distribution of the obtained MRE was evaluated by a field emission scanning electron microscopy (FESEM). The influence of constant low strain at 0.01% is the continuing concern within the linear viscoelastic (LVE) region of MRE. Stress relaxation plays a significant role in the life cycle of MRE and revealed that storage modulus was reduced by 8.7%, normal force has weakened by 27%, and stress performance was reduced by 6.88% along approximately 84,000 s test duration time. This time scale was the longest ever reported being undertaken in the MRE stress relaxation study. Novel micro-mechanisms that responsible for the depleted performance of MRE was obtained by microstructurally observation using FESEM and in-phase mode of atomic force microscope (AFM). Attempts have been made to correlate strain localization produced by stress relaxation, with molecular deformation in MRE amorphous matrix. Exceptional attention was focused on the development of molecular slippage, disentanglement, microplasticity, microphase separation, and shear bands. The relation between these microstructural phenomena and the viscoelastic properties of MRE was diffusely defined and discussed. The presented MRE is homogeneous with uniform distribution of CIP. The most significant recent developments of systematic correlation between the effects of microstructural deformation and durability performance of MRE under stress relaxation has been observed and evaluated.

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 β/ν.

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.

Prediction of creep failure time using machine learning

A subcritical load on a disordered material can induce creep damage. The creep rate in this case exhibits three temporal regimes viz. an initial decelerating regime followed by a steady-state regime and a stage of accelerating creep that ultimately leads to catastrophic breakdown. Due to the statistical regularities in the creep rate, the time evolution of creep rate has often been used to predict residual lifetime until catastrophic breakdown. However, in disordered samples, these efforts met with limited success. Nevertheless, it is clear that as the failure is approached, the damage become increasingly spatially correlated, and the spatio-temporal patterns of acoustic emission, which serve as a proxy for damage accumulation activity, are likely to mirror such correlations. However, due to the high dimensionality of the data and the complex nature of the correlations it is not straightforward to identify the said correlations and thereby the precursory signals of failure. Here we use supervised machine learning to estimate the remaining time to failure of samples of disordered materials. The machine learning algorithm uses as input the temporal signal provided by a mesoscale elastoplastic model for the evolution of creep damage in disordered solids. Machine learning algorithms are well-suited for assessing the proximity to failure from the time series of the acoustic emissions of sheared samples. We show that materials are relatively more predictable for higher disorder while are relatively less predictable for larger system sizes. We find that machine learning predictions, in the vast majority of cases, perform substantially better than other prediction approaches proposed in the literature.

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.

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.

Bending Nanoindentation and Plasticity Noise in FCC Single and Polycrystals

We present a high-throughput nanoindentation study of in situ bending effects on incipient plastic deformation behavior of polycrystalline and single-crystalline pure aluminum and pure copper at ultranano depths (< 200 nm). We find that hardness displays a statistically inverse dependence on in-plane stress for indentation depths smaller than 10 nm, and the dependence disappears for larger indentation depths. In contrast, plastic noise in the nanoindentation force and displacement displays statistically robust noise features, independently of applied stresses. Our experimental results suggest the existence of a regime in Face Centered Cubic (FCC) crystals where ultranano hardness is sensitive to residual applied stresses, but plasticity pop-in noise is insensitive to it.

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.

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.

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.

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.