Local yield stress statistics in model amorphous solids

An elastoplastic model approach for the relaxation dynamics of active glasses

How activity affects the glassy dynamics is crucial for several biological processes. Furthermore, active glasses offer fascinating phenomenologies, extend the scope of equilibrium glass-forming liquids, and can provide novel insights into the original problem. We introduce a family of novel approaches to investigating the relaxation dynamics of active glasses via an active elastoplastic model (EPM). These approaches describe the relaxation dynamics via local plastic yielding and can provide improved insights as we can study various aspects of the system separately. Activity enters the model via three crucial features: activity-mediated plastic yielding, activated barrier crossing, and persistent rotational dynamics of the yielding direction. We first consider a minimal active EPM that adds the effect of active yielding to a thermal EPM. We show that this active EPM captures the known results of active glasses within a reasonable parameter space. The results also agree well with the analytical results for active glasses when activity is small. The minimal model breaks down at very low temperatures where other effects become important. Looking at the broader model class, we demonstrate that whereas active yielding primarily dominates the relaxation dynamics, the persistence of the yielding direction governs the dynamic heterogeneity in active glasses.

Anomalous Softness in Amorphous Matter in the Reversible Plastic Regime

We study an elastoplastic model of an amorphous solid subject to athermal quasistatic cyclic shear strain. We focus on cycling amplitudes in the so-called reversible-plastic regime where, after a transient, the system locks into a hysteretic limit cycle and returns to the same microscopic configuration after one or more strain cycles. We show that the ground state energy of the terminal limit cycle decreases with increasing cycling amplitude. In analogy to an annealed alloy or an aged colloidal glass, one would expect the states with lower energy to be mechanically harder and to require larger stresses and strains to trigger microscopic rearrangements. However, we show the opposite result: the systems with lower energy cycled at higher strain amplitude are mechanically softer and begin to exhibit plastic rearrangements at smaller stresses and strains within the cycle. We explain this anomaly quantitatively in terms of Eshelby inclusion theory where an inclusion is subjected to a particular negative stress value after it undergoes a yielding event. These results point the way toward measurements to be conducted in experiments and particle-based computer simulations on cyclically sheared amorphous solids.

Selecting relevant structural features for glassy dynamics by information imbalance

We numerically investigate the identification of relevant structural features that contribute to the dynamical heterogeneity in a model glass-forming liquid. By employing the recently proposed information imbalance technique, we select these features from a range of physically motivated descriptors. This selection process is performed in a supervised manner (using both dynamical and structural data) and an unsupervised manner (using only structural data). We then apply the selected features to predict future dynamics using a machine learning technique. One of the advantages of the information imbalance technique is that it does not assume any model a priori, i.e., it is a non-parametric method. Finally, we discuss the potential applications of this approach in identifying the dominant mechanisms governing the glassy slow dynamics.

Thermally activated intermittent flow in amorphous solids

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

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

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

Memory of shear flow in soft jammed materials

Cessation of flow in yield stress fluids results in a stress relaxation process that eventually leads to a finite residual stress. Both the rate of stress relaxation and the magnitude of the residual stresses systematically depend on the preceding flow conditions. To assess the microscopic origin of this memory effect, we combine experiments with large-scale computer simulations, exploring the behavior of jammed suspensions of soft repulsive particles. A spatiotemporal analysis of particle motion reveals that memory formation during flow is primarily governed by the emergence of domains of spatially correlated nonaffine displacements. These domains imprint the configuration of stress imbalances that drive dynamics upon flow cessation, as evidenced by a striking equivalence of the spatial correlation patterns in particle displacements observed during flow and upon flow cessation. Additional contributions to stress relaxation result from the particle packing that reorganizes to minimize the resistance to flow by decreasing the number of locally stiffer configurations. Regaining rigidity upon flow cessation drives further relaxation and effectively sets the magnitude of the residual stress. Our findings highlight that flow in yield stress fluids can be seen as a training process during which the material stores information of the flowing state through the development of domains of correlated particle displacements and the reorganization of particle packings optimized to sustain the flow. This encoded memory can then be retrieved in flow cessation experiments.

Topological characterization of rearrangements in amorphous solids

In amorphous materials, plasticity is localized and occurs as shear transformations. It was recently shown by Wu et al. that these shear transformations can be predicted by applying topological defect concepts developed for liquid crystals to an analysis of vibrational eigenmodes [Z. W. Wu et al., Nat. Commun. 14, 2955 (2023)10.1038/s41467-023-38547-w]. This study relates the -1 topological defects to the displacement fields expected of an Eshelby inclusion, which are characterized by an orientation and the magnitude of the eigenstrain. A corresponding orientation and magnitude can be defined for each defect using the local displacement field around each defect. These parameters characterize the plastic stress relaxation associated with the local structural rearrangement and can be extracted using the fit to either the global displacement field or the local field. Both methods provide a reasonable estimation of the molecular-dynamics-measured stress drop, confirming the localized nature of the displacements that control both long-range deformation and stress relaxation.

Universal stress correlations in crystalline and amorphous packings

We present a universal characterization of stress correlations in athermal systems, across crystalline to amorphous packings. Via numerical analysis of static configurations of particles interacting through harmonic as well as Lennard-Jones potentials, for a variety of preparation protocols and ranges of microscopic disorder, we show that the properties of the stress correlations at large lengthscales are surprisingly universal across all situations, independent of structural correlations, or the correlations in orientational order. In the near-crystalline limit, we present exact results for the stress correlations for both models, which work surprisingly well at large lengthscales, even in the amorphous phase. Finally, we study the differences in stress fluctuations across the amorphization transition, where stress correlations reveal the loss of periodicity in the structure at short lengthscales with increasing disorder.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Thawed matrix method for computing local mechanical properties of amorphous solids

We present a method for computing locally varying nonlinear mechanical properties in particle simulations of amorphous solids. Plastic rearrangements outside a probed region are suppressed by introducing an external field that directly penalizes large nonaffine displacements. With increasing strength of the field, plastic deformation can be localized. We characterize the distribution of local plastic yield stresses (residual local stresses to instability) with our approach and assess the correlation of their spatial maps with plastic activity in a model two-dimensional amorphous solid. Our approach reduces artifacts inherent in a previous method known as the "frozen matrix" approach that enforces fully affine deformation and improves the prediction of plastic rearrangements from structural information.

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.

Topology of vibrational modes predicts plastic events in glasses

The plastic deformation of crystalline materials can be understood by considering their structural defects such as disclinations and dislocations. Although also glasses are solids, their structure resembles closely the one of a liquid and hence the concept of structural defects becomes ill-defined. As a consequence it is very challenging to rationalize on a microscopic level the mechanical properties of glasses close to the yielding point and to relate plastic events to structural properties. Here we investigate the topological characteristics of the eigenvector field of the vibrational excitations of a two-dimensional glass model, notably the geometric arrangement of the topological defects as a function of vibrational frequency. We find that if the system is subjected to a quasistatic shear, the location of the resulting plastic events correlate strongly with the topological defects that have a negative charge. Our results provide thus a direct link between the structure of glasses prior their deformation and the plastic events during deformation.

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.

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

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

Relevance of Shear Transformations in the Relaxation of Supercooled Liquids

While deeply supercooled liquids exhibit divergent viscosity and increasingly heterogeneous dynamics as the temperature drops, their structure shows only seemingly marginal changes. Understanding the nature of relaxation processes in this dramatic slowdown is key for understanding the glass transition. Here, we show by atomistic simulations that the heterogeneous dynamics of glass-forming liquids strongly correlate with the local residual plastic strengths along soft directions computed in the initial inherent structures. The correlation increases with decreasing temperature and is maximum in the vicinity of the relaxation time. For the lowest temperature investigated, this maximum is comparable with the best values from the literature dealing with the structure-property relationship. However, the nonlinear probe of the local shear resistance in soft directions provides here a real-space picture of relaxation processes. Our detection method of thermal rearrangements allows us to investigate the first passage time statistics and to study the scaling between the activation energy barriers and the residual plastic strengths. These results shed new light on the nature of relaxations of glassy systems by emphasizing the analogy between the thermal relaxations in viscous liquids and the plastic shear transformation in amorphous solids.

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.

Static self-induced heterogeneity in glass-forming liquids: Overlap as a microscope

We propose and numerically implement a local probe of the static self-induced heterogeneity characterizing glass-forming liquids. This method relies on the equilibrium statistics of the overlap between pairs of configurations measured in mesoscopic cavities with unconstrained boundaries. By systematically changing the location of the probed cavity, we directly detect spatial variations of the overlap fluctuations. We provide a detailed analysis of the statistics of a local estimate of the configurational entropy, and we infer an estimate of the surface tension between amorphous states, ingredients that are both at the basis of the random first-order transition theory of glass formation. Our results represent the first direct attempt to visualize and quantify the self-induced heterogeneity underpinning the thermodynamics of glass formation. They pave the way for the development of coarse-grained effective theories and for a direct assessment of the role of thermodynamics in the activated dynamics of deeply supercooled liquids.

Beyond the Average: Spatial and Temporal Fluctuations in Oxide Glass-Forming Systems

Atomic structure dictates the performance of all materials systems; the characteristic of disordered materials is the significance of spatial and temporal fluctuations on composition-structure-property-performance relationships. Glass has a disordered atomic arrangement, which induces localized distributions in physical properties that are conventionally defined by average values. Quantifying these statistical distributions (including variances, fluctuations, and heterogeneities) is necessary to describe the complexity of glass-forming systems. Only recently have rigorous theories been developed to predict heterogeneities to manipulate and optimize glass properties. This article provides a comprehensive review of experimental, computational, and theoretical approaches to characterize and demonstrate the effects of short-, medium-, and long-range statistical fluctuations on physical properties (e.g., thermodynamic, kinetic, mechanical, and optical) and processes (e.g., relaxation, crystallization, and phase separation), focusing primarily on commercially relevant oxide glasses. Rigorous investigations of fluctuations enable researchers to improve the fundamental understanding of the chemistry and physics governing glass-forming systems and optimize structure-property-performance relationships for next-generation technological applications of glass, including damage-resistant electronic displays, safer pharmaceutical vials to store and transport vaccines, and lower-attenuation fiber optics. We invite the reader to join us in exploring what can be discovered by going beyond the average.

Using delaunay triangularization to characterize non-affine displacement fields during athermal, quasistatic deformation of amorphous solids

We investigate the non-affine displacement fields that occur in two-dimensional Lennard-Jones models of metallic glasses subjected to athermal, quasistatic simple shear (AQS). During AQS, the shear stress versus strain displays continuous quasi-elastic segments punctuated by rapid drops in shear stress, which correspond to atomic rearrangement events. We capture all information concerning the atomic motion during the quasi-elastic segments and shear stress drops by performing Delaunay triangularizations and tracking the deformation gradient tensor F_α associated with each triangle α. To understand the spatio-temporal evolution of the displacement fields during shear stress drops, we calculate F_α along minimal energy paths from the mechanically stable configuration immediately before to that after the stress drop. We find that quadrupolar displacement fields form and dissipate both during the quasi-elastic segments and shear stress drops. We then perform local perturbations (rotation, dilation, simple and pure shear) to single triangles and measure the resulting displacement fields. We find that local pure shear deformations of single triangles give rise to mostly quadrupolar displacement fields, and thus pure shear strain is the primary type of local strain that is activated by bulk, athermal quasistatic simple shear. Other local perturbations, e.g. rotations, dilations, and simple shear of single triangles, give rise to vortex-like and dipolar displacement fields that are not frequently activated by bulk AQS. These results provide fundamental insights into the non-affine atomic motion that occurs in driven, glassy materials.

Atomic vibration as an indicator of the propensity for configurational rearrangements in metallic glasses

In a metallic glass (MG), the propensity for atomic rearrangements varies spatially from location to location in the amorphous solid, making the prediction of their likelihood a major challenge. One can attack this problem from the "structure controls properties" standpoint. But all the current structure-centric parameters are mostly based on local atomic packing information limited to short-range order, hence falling short in reliably forecasting how the local region would respond to external stimuli (e.g., temperature and/or stress). Alternatively, one can use indicators informed by physical properties to bridge the static structure on the one hand, and the response of the local configuration on the other. A sub-group of such physics-informed quantities consists of atomic vibration parameters, which will be singled out as the focus of this article. Here we use the Cu₆₄Zr₃₆ alloy to systematically demonstrate the following two points, all using a single model MG. First, we show in a comprehensive manner the interrelation among common vibrational parameters characterizing the atomic vibrational amplitude and frequency, including the atomic mean square displacement, flexibility volume, participation fraction in the low-frequency vibrational modes and boson peak intensity. Second, we demonstrate that these vibrational parameters fare much better than purely static structural parameters based on local geometrical packing in providing correlation with the propensity for local configurational transitions. These vibrational parameters also share a correlation length similar to that in structural rearrangements induced by external stimuli. This success, however, also poses a challenge, as it remains to be elucidated as to why short-time dynamical (vibrational) behavior at the bottom of the energy basin can be exploited to project the height of the energy barrier for cross-basin activities and in turn the propensity for locally collective atomic rearrangements.

Elastoplastic Approach Based on Microscopic Insights for the Steady State and Transient Dynamics of Sheared Disordered Solids

We develop a framework to study the mechanical response of athermal amorphous solids via a coupling of mesoscale and microscopic models. Using measurements of coarse-grained quantities from simulations of dense disordered particulate systems, we present a coherent elastoplastic model approach for deformation and flow of yield stress materials. For a given set of parameters, this model allows us to match consistently transient and steady state features of driven disordered systems with diverse preparation histories under both applied shear-rate and creep protocols.

Predicting orientation-dependent plastic susceptibility from static structure in amorphous solids via deep learning

It has been a long-standing materials science challenge to establish structure-property relations in amorphous solids. Here we introduce a rotationally non-invariant local structure representation that enables different predictions for different loading orientations, which is found essential for high-fidelity prediction of the propensity for stress-driven shear transformations. This novel structure representation, when combined with convolutional neural network (CNN), a powerful deep learning algorithm, leads to unprecedented accuracy for identifying atoms with high propensity for shear transformations (i.e., plastic susceptibility), solely from the static structure in both two- and three-dimensional model glasses. The data-driven models trained on samples at one composition and a given processing history are found transferrable to glass samples with different processing histories or at different compositions in the same alloy system. Our analysis of the new structure representation also provides valuable insight into key atomic packing features that influence the local mechanical response and its anisotropy in glasses.

Inferring the particle-wise dynamics of amorphous solids from the local structure at the jamming point

Jamming is a phenomenon shared by a wide variety of systems, such as granular materials, foams, and glasses in their high density regime. This has motivated the development of a theoretical framework capable of explaining many of their static critical properties with a unified approach. However, the dynamics occurring in the vicinity of the jamming point has received little attention and the problem of finding a connection with the local structure of the configuration remains unexplored. Here we address this issue by constructing physically well defined structural variables using the information contained in the network of contacts of jammed configurations, and then showing that such variables yield a resilient statistical description of the particle-wise dynamics near this critical point. Our results are based on extensive numerical simulations of systems of spherical particles that allow us to statistically characterize the trajectories of individual particles in terms of their first two moments. We first demonstrate that, besides displaying a broad distribution of mobilities, particles may also have preferential directions of motion. Next, we associate each of these features with a structural variable computed uniquely in terms of the contact vectors at jamming, obtaining considerably high statistical correlations. The robustness of our approach is confirmed by testing two types of dynamical protocols, namely molecular dynamics and Monte Carlo, with different types of interaction. We also provide evidence that the dynamical regime we study here is dominated by anharmonic effects and therefore it cannot be described properly in terms of vibrational modes. Finally, we show that correlations decay slowly and in an interaction-independent fashion, suggesting a universal rate of information loss.

Simple and Broadly Applicable Definition of Shear Transformation Zones

Plastic deformation in amorphous solids is known to be carried by stress-induced localized rearrangements of a few tens of particles, accompanied by the conversion of elastic energy to heat. Despite their central role in determining how glasses yield and break, the search for a simple and generally applicable definition of the precursors of those plastic rearrangements-the so-called shear transformation zones (STZs)-is still ongoing. Here we present a simple definition of STZs-based solely on the harmonic approximation of a glass's energy. We explain why and demonstrate directly that our proposed definition of plasticity carriers in amorphous solids is more broadly applicable compared to anharmonic definitions put forward previously. Finally, we offer an open-source library that analyzes low-lying STZs in computer glasses and in laboratory materials such as dense colloidal suspensions for which the harmonic approximation is accessible. Our results constitute a physically motivated methodological advancement towards characterizing mechanical disorder in glasses, and understanding how they yield.

Time-resolved shear transformations in the transient plastic regime of sheared amorphous silicon

The accumulation of shear transformations (STs) in space and time is responsible for plastic deformation in amorphous solids. Here we study the effect of finite strain rates on STs during simulations of athermal shear deformation in an atomistic model of amorphous silicon. We present a time-resolved analysis of STs by mapping the plastic events identified in the atomistic simulations on a collection of Eshelby inclusions, which are characterized in terms of number, effective volume, lifetime, and orientation. Our analysis led us to distinguish between small and large events. We find that the main effect of a lower strain rate is to allow for a larger number of small events, roughly identified by an effective volume γ_{0}V_{0}<20 Å^{3}, while the number and characteristics of larger events are surprisingly independent of the strain rate. We show that at low strains, the decrease of the stress observed at lower strain rates is mainly due to the excess of small events, while at larger strains, when the glass approaches the yield point where a shear band forms, larger events start to play a role and organize due to their elastic interactions. This phenomenology is compared with the predictions of mesoscale elastoplastic models. The technique developed here can be used as a systematic tool to analyze plasticity during molecular dynamics simulations. It can also give valuable information to develop physically grounded mesoscale models of plasticity, providing quantitative predictions of the mechanical properties of amorphous materials.

Residual stress distributions in amorphous solids from atomistic simulations

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

Origin of reversible and irreversible atomic-scale rearrangements in a model two-dimensional network glass

In this contribution, we investigate the fundamental mechanism of plasticity in a model two-dimensional network glass. The glass is generated by using a Monte Carlo bond-switching algorithm and subjected to athermal simple shear deformation, followed by subsequent unloading at selected deformation states. This enables us to investigate the topological origin of reversible and irreversible atomic-scale rearrangements. It is shown that some events that are triggered during loading recover during unloading, while some do not. Thus, two kinds of elementary plastic events are observed, which can be linked to the network topology of the model glass.

Computational study of transient shear banding in soft jammed solids

We have designed three-dimensional numerical simulations of a soft spheres model, with size polidispersity and in athermal conditions, to study the transient shear banding that occurs during yielding of jammed soft solids. We analyze the effects of different types of drag coefficients used in the simulations and compare the results obtained using Lees-Edwards periodic boundary conditions with the case in which the same model solid is confined between two walls. The specific damping mechanism and the different boundary conditions indeed modify the load curves and the velocity profiles in the transient regime. Nevertheless, we find that the presence of a stress overshoot and of a related transient banding phenomenon, for large enough samples, is a robust feature for overdamped systems, where their presence do not depend on the specific drag used and on the different boundary conditions.

Emergence and persistence of flow inhomogeneities in the yielding and fluidization of dense soft solids

In three-dimensional computer simulations of model non-Brownian jammed suspensions, we compute the time required to reach homogeneous flow upon yielding, by analyzing stresses and particle packing at different shear rates, with and without confinement. We show that the stress overshoot and persistent shear banding preceding the complete fluidization are controlled by the presence of overconstrained microscopic domains in the initial solids. Such domains, identifiable with icosahedrally packed regions in the model used, allow for stress accumulation during the shear startup. Their structural reorganization under deformation controls the emergence and the persistence of the shear banding.

Origin of the Bauschinger Effect in Amorphous Solids

We study the structural origin of the Bauschinger effect by accessing numerically the local plastic thresholds in the steady state flow of a two-dimensional model glass under athermal quasistatic deformation. More specifically, we compute the local residual strength, Δτ^{c}, for arbitrary loading orientations and find that plastic deformation generically induces material polarization, i.e., a forward-backward asymmetry in the Δτ^{c} distribution. In steady plastic flow, local packings are on average closer to forward (rather than backward) instabilities, due to the stress-induced bias of barriers. However, presumably due to mechanical noise, a significant fraction of zones lie close to reverse (backward) yielding, as the distribution of Δτ^{c} for reverse shearing extends quasilinearly down to zero local residual strength. By constructing an elementary model of the early plastic response, we then show that unloading causes reverse plasticity of a growing amplitude, i.e., reverse softening, while it shifts away forward-yielding barriers. This result in an inversion of polarization in the low-Δτ^{c} region and, consequently, in the Bauschinger effect. This scenario is quite generic, which explains the pervasiveness of the effect.

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.

Elastic avalanches reveal marginal behavior in amorphous solids

Mechanical deformation of amorphous solids can be described as consisting of an elastic part in which the stress increases linearly with strain, up to a yield point at which the solid either fractures or starts deforming plastically. It is well established, however, that the apparent linearity of stress with strain is actually a proxy for a much more complex behavior, with a microscopic plasticity that is reflected in diverging nonlinear elastic coefficients. Very generally, the complex structure of the energy landscape is expected to induce a singular response to small perturbations. In the athermal quasistatic regime, this response manifests itself in the form of a scale-free plastic activity. The distribution of the corresponding avalanches should reflect, according to theoretical mean-field calculations [S. Franz and S. Spigler, Phys. Rev. E 95, 022139 (2017)], the geometry of phase space in the vicinity of a typical local minimum. In this work, we characterize this distribution for simple models of glass-forming systems, and we find that its scaling is compatible with the mean-field predictions for systems above the jamming transition. These systems exhibit marginal stability, and scaling relations that hold in the stationary state are examined and confirmed in the elastic regime. By studying the respective influence of system size and age, we suggest that marginal stability is systematic in the thermodynamic limit.

Anisotropic structural predictor in glassy materials

There is growing evidence that relaxation in glassy materials, both spontaneous and externally driven, is mediated by localized soft spots. Recent progress made it possible to identify the soft spots inside glassy structures and to quantify their degree of softness. These softness measures, however, are typically scalars, not taking into account the tensorial, anisotropic nature of soft spots, which implies orientation-dependent coupling to external deformation. Here, we derive from first principles the linear response coupling between the local heat capacity of glasses, previously shown to provide a measure of glassy softness, and external deformation in different directions. We first show that this linear response quantity follows an anomalous, fat-tailed distribution related to the universal ω^{4} density of states of quasilocalized, nonphononic excitations in glasses. We then construct a structural predictor as the product of the local heat capacity and its linear response to external deformation, and show that it offers an enhanced predictability of plastic rearrangements under deformation in different directions, compared to the purely scalar predictor.

Machine learning characterization of structural defects in amorphous packings of dimers and ellipses

Structural defects within amorphous packings of symmetric particles can be characterized using a machine learning approach that incorporates structure functions of radial distances and angular arrangement. This yields a scalar field, softness, that correlates with the probability that a particle is about to be rearranged. However, when particle shapes are elongated, as in the case of dimers and ellipses, we find that the standard structure functions produce imprecise softness measurements. Moreover, ellipses exhibit deformation profiles in stark contrast to circular particles. In order to account for the effects of orientation and alignment, we introduce structure functions to recover the predictive performance of softness, as well as provide physical insight into local and extended dynamics. We study a model disordered solid, a bidisperse two-dimensional granular pillar, driven by uniaxial compression and composed entirely of monomers, dimers, or ellipses. We demonstrate how the computation of softness via a support vector machine extends to dimers and ellipses with the introduction of orientational structure functions. Then we highlight the spatial extent of rearrangements and defects, as well as their cross correlation, for each particle shape. Finally, we demonstrate how an additional machine learning algorithm, recursive feature elimination, provides an avenue to better understand how softness arises from particular structural aspects. We identify the most crucial structure functions in determining softness and discuss their physical implications.

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.