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

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.

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.

Low-energy quasilocalized excitations in structural glasses

Glassy solids exhibit a wide variety of generic thermomechanical properties, ranging from universal anomalous specific heat at cryogenic temperatures to nonlinear plastic yielding and failure under external driving forces, which qualitatively differ from their crystalline counterparts. For a long time, it has been believed that many of these properties are intimately related to nonphononic, low-energy quasilocalized excitations (QLEs) in glasses. Indeed, recent computer simulations have conclusively revealed that the self-organization of glasses during vitrification upon cooling from a melt leads to the emergence of such QLEs. In this Perspective, we review developments over the past three decades toward understanding the emergence of QLEs in structural glasses and the degree of universality in their statistical and structural properties. We discuss the challenges and difficulties that hindered progress in achieving these goals and review the frameworks put forward to overcome them. We conclude with an outlook on future research directions and open questions.

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.

Extracting the properties of quasilocalized modes in computer glasses: Long-range continuum fields, contour integrals, and boundary effects

Low-frequency nonphononic modes and plastic rearrangements in glasses are spatially quasilocalized, i.e., they feature a disorder-induced short-range core and known long-range decaying elastic fields. Extracting the unknown short-range core properties, potentially accessible in computer glasses, is of prime importance. Here we consider a class of contour integrals, performed over the known long-range fields, which are especially designed for extracting the core properties. We first show that, in computer glasses of typical sizes used in current studies, the long-range fields of quasilocalized modes experience boundary effects related to the simulation box shape and the widely employed periodic boundary conditions. In particular, image interactions mediated by the box shape and the periodic boundary conditions induce the fields' rotation and orientation-dependent suppression of their long-range decay. We then develop a continuum theory that quantitatively predicts these finite-size boundary effects and support it by extensive computer simulations. The theory accounts for the finite-size boundary effects and at the same time allows the extraction of the short-range core properties, such as their typical strain ratios and orientation. The theory is extensively validated in both two and three dimensions. Overall, our results offer a useful tool for extracting the intrinsic core properties of nonphononic modes and plastic rearrangements in computer glasses.

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.

Transition rates for slip-avalanches in soft athermal disks under quasi-static simple shear deformations

We study slip-avalanches in two-dimensional soft athermal disks by quasi-static simulations of simple shear deformations. Sharp drops in shear stress, or slip-avalanches, are observed intermittently during steady state. Such stress drops are caused by restructuring of the contact networks, accompanied by drastic changes of the interaction forces, Δf. The changes of the forces happen heterogeneously in space, indicating that collective non-affine motions of the disks are most pronounced when slip-avalanches occur. We analyze and predict the statistics for the force changes, Δf, by transition rates of the contact forces and angles, where slip-avalanches are characterized by wide power-law tails. We find that the transition rates are described by a q-Gaussian distribution regardless of the area fraction of the disks. Because the transition rates quantify structural changes of the force-chains, our findings are an important step towards linking macroscopic observations to a microscopic theory of slip-avalanches in the experimentally accessible quasi-static regime.