Time-dependent elastic response to a local shear transformation in amorphous solids

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.

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.

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.

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.

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.

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypoelastoplastic mechanical model, and develop a projection method to enforce quasistatic equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and which imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

Compaction of elastic granular materials: inter-particles friction effects and plastic events

The uni-axial compaction of granular materials made of elastic neo-Hookean particles is investigated in the quasi-static regime. Two-dimensional disk assemblies are simulated using the Finite Element model coupled with Contact Dynamics method for dealing both with finite deformations of the particles and contact interactions. Due to large deformations of the particles, the packing fraction of the system increases continuously during the compaction process, reaching values close to 1. The influence of the coefficient of friction between the particles on the macroscopic and micro-structural behaviors of the system is thoroughly discussed.

Glassy dynamics in models of confluent tissue with mitosis and apoptosis

Recent work on particle-based models of tissues has suggested that any finite rate of cell division and cell death is sufficient to fluidize an epithelial tissue. At the same time, experimental evidence has indicated the existence of glassy dynamics in some epithelial layers despite continued cell cycling. To address this discrepancy, we quantify the role of cell birth and death on glassy states in confluent tissues using simulations of an active vertex model that includes cell motility, cell division, and cell death. Our simulation data is consistent with a simple ansatz in which the rate of cell-life cycling and the rate of relaxation of the tissue in the absence of cell cycling contribute independently and additively to the overall rate of cell motion. Specifically, we find that a glass-like regime with caging behavior indicated by subdiffusive cell displacements can be achieved in systems with sufficiently low rates of cell cycling.

Probing the Degree of Heterogeneity within a Shear Band of a Model Glass

Recent experiments provide evidence for density variations along shear bands in metallic glasses with a length scale of a few hundred nanometers. Via molecular dynamics simulations of a generic binary glass model, here we show that this is strongly correlated with variations of composition, coordination number, viscosity, and heat generation. Individual shear events along the shear band path show a mean distance of a few nanometers, comparable to recent experimental findings on medium range order. The aforementioned variations result from these localized perturbations, mediated by elasticity.

Creep dynamics of athermal amorphous materials: a mesoscopic approach

Yield stress fluids display complex dynamics, in particular when driven into the transient regime between the solid and the flowing state. Inspired by creep experiments on dense amorphous materials, we implement mesoscale elasto-plastic descriptions to analyze such transient dynamics in athermal systems. Both our mean-field and space-dependent approaches consistently reproduce the typical experimental strain rate responses to different applied steps in stress. Moreover, they allow us to understand basic processes involved in the strain rate slowing down (creep) and the strain rate acceleration (fluidization) phases. The fluidization time increases in a power-law fashion as the applied external stress approaches a static yield stress. This stress value is related to the stress over-shoot in shear start-up experiments, and it is known to depend on sample preparation and age. By calculating correlations of the accumulated plasticity in the spatially resolved model, we reveal different modes of cooperative motion during the creep dynamics.

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.

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.

Correlation and shear bands in a plastically deformed granular medium

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

Mean-Field Scenario for the Athermal Creep Dynamics of Yield-Stress Fluids

We develop a theoretical description based on an existent mean-field model for the transient dynamics prior to the steady flow of yielding materials. The mean-field model not only reproduces the experimentally observed nonlinear time dependence of the shear-rate response to an external stress, but also allows for the determination of the different physical processes involved in the onset of the reacceleration phase after the initial slowing down and a distinct fluidization phase. The fluidization time displays a power-law dependence on the distance of the applied stress to an age-dependent yield stress, which is not universal but strongly dependent on initial conditions.

Effect of composition and pressure on the shear strength of sodium silicate glasses: An atomic scale simulation study

The elastoplastic behavior of sodium silicate glasses is studied at different scales as a function of composition and pressure, with the help of quasistatic atomistic simulations. The samples are first compressed and then sheared at constant pressure to calculate yield strength and permanent plastic deformations. Changes occurring in the global response are then compared to the analysis of local plastic rearrangements and strain heterogeneities. It is shown that the plastic response results from the succession of well-identified localized irreversible deformations occurring in a nanometer-size area. The size and the number of these local rearrangements, as well as the amount of internal deviatoric and volumetric plastic deformation, are sensitive to the composition and to the pressure. In the early stages of the deformation, plastic rearrangements are driven by sodium mobility. Consequently, the elastic yield strength decreases when the sodium content increases, and the same when pressure increases. Finally, good correlation was found between global and local stress-strain relationships, reinforcing again the role of sodium ions as local initiators of the plastic behavior observed at larger scales.

Plastic response and correlations in athermally sheared amorphous solids

The onset of irreversible deformation in low-temperature amorphous solids is due to the accumulation of elementary events, consisting of spatially and temporally localized atomic rearrangements involving only a few tens of atoms. Recently, numerical and experimental work addressed the issue of spatiotemporal correlations between these plastic events. Here, we provide further insight into these correlations by investigating, via molecular dynamics (MD) simulations, the plastic response of a two-dimensional amorphous solid to artificially triggered local shear transformations. We show that while the plastic response is virtually absent in as-quenched configurations, it becomes apparent if a shear strain was previously imposed on the system. Plastic response has a fourfold symmetry, which is characteristic of the shear stress redistribution following the local transformation. At high shear rate we report evidence for a fluctuation-dissipation relation, connecting plastic response and correlation, which seems to break down if lower shear rates are considered.

Eshelby inclusions in granular matter: Theory and simulations

We present a numerical implementation of an active inclusion in a granular material submitted to a biaxial test. We discuss the dependence of the response to this perturbation on two parameters: the intragranular friction coefficient on one hand, and the degree of the loading on the other hand. We compare the numerical results to theoretical predictions taking into account the change of volume of the inclusion as well as the anisotropy of the elastic matrix.

Analysis of structural correlations in a model binary 3D liquid through the eigenvalues and eigenvectors of the atomic stress tensors

It is possible to associate with every atom or molecule in a liquid its own atomic stress tensor. These atomic stress tensors can be used to describe liquids' structures and to investigate the connection between structural and dynamic properties. In particular, atomic stresses allow to address atomic scale correlations relevant to the Green-Kubo expression for viscosity. Previously correlations between the atomic stresses of different atoms were studied using the Cartesian representation of the stress tensors or the representation based on spherical harmonics. In this paper we address structural correlations in a 3D model binary liquid using the eigenvalues and eigenvectors of the atomic stress tensors. This approach allows to interpret correlations relevant to the Green-Kubo expression for viscosity in a simple geometric way. On decrease of temperature the changes in the relevant stress correlation function between different atoms are significantly more pronounced than the changes in the pair density function. We demonstrate that this behaviour originates from the orientational correlations between the eigenvectors of the atomic stress tensors. We also found correlations between the eigenvalues of the same atomic stress tensor. For the studied system, with purely repulsive interactions between the particles, the eigenvalues of every atomic stress tensor are positive and they can be ordered: λ1 ≥ λ2 ≥ λ3 ≥ 0. We found that, for the particles of a given type, the probability distributions of the ratios (λ2/λ1) and (λ3/λ2) are essentially identical to each other in the liquids state. We also found that λ2 tends to be equal to the geometric average of λ1 and λ3. In our view, correlations between the eigenvalues may represent "the Poisson ratio effect" at the atomic scale.

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.

Effects of Inertia on the Steady-Shear Rheology of Disordered Solids

We study the finite-shear-rate rheology of disordered solids by means of molecular dynamics simulations in two dimensions. By systematically varying the damping strength ζ in the low-temperature limit, we identify two well-defined flow regimes, separated by a thin (temperature-dependent) crossover region. In the overdamped regime, the athermal rheology is governed by the competition between elastic forces and viscous forces, whose ratio gives the Weissenberg number Wi∝ζγ[over ˙]; the macroscopic stress Σ follows the frequently encountered Herschel-Bulkley law Σ=Σ_{0}+ksqrt[Wi], with yield stress Σ_{0}>0. In the underdamped (inertial) regime, dramatic changes in the rheology are observed for low damping: the flow curve becomes nonmonotonic. This change is not caused by longer-lived correlations in the particle dynamics at lower damping; instead, for weak dissipation, the sample heats up considerably due to, and in proportion to, the driving. By thermostating more or less underdamped systems, we are able to link quantitatively the rheology to the kinetic temperature and the shear rate, rescaled with Einstein's vibration frequency.

Glass transition of two-dimensional 80-20 Kob-Andersen model at constant pressure

We reconsider numerically the two-dimensional version of the Kob-Andersen model (KA2d) with a fraction of 80% of large spheres. A constant moderate pressure is imposed while the temperature T is systematically quenched from the liquid limit through the glass transition at [Formula: see text] down to very low temperatures. Monodisperse Lennard-Jones (mdLJ) bead systems, forming a crystal phase at low temperatures, are used to highlight several features of the KA2d model. As can be seen, e.g. from the elastic shear modulus G(T), determined using the stress-fluctuation formalism, our KA2d model is a good glass-former. A continuous cusp-singularity, [Formula: see text] with [Formula: see text], is observed in qualitative agreement with other recent numerical and theoretical work, however in striking conflict with the additive jump discontinuity predicted by mode-coupling theory.

Role of inertia in the rheology of amorphous systems: A finite-element-based elastoplastic model

A simple finite-element analysis with varying damping strength is used to model the athermal shear rheology of densely packed glassy systems at a continuum level. We focus on the influence of dissipation on bulk rheological properties. Our numerical studies, done over a wide range of damping coefficients, identify two well-separated rheological regimes along with a crossover region controlled by a critical damping. In the overdamped limit, inertial effects are negligible and the rheological response is well described by the commonly observed Herschel-Bulkley equation. In stark contrast, inertial vibrations in the underdamped regime prompt a significant drop in the mean-stress level, leading to a nonmonotonic constitutive relation. The observed negative slope in the flow curve, which is a signature of mechanical instability and thus permanent shear banding, arises from the sole influence of inertia, in qualitative agreement with the recent molecular dynamics study of Nicolas et al., Phys. Rev. Lett. 116, 058303 (2016).

Analysis of spatial correlations in a model two-dimensional liquid through eigenvalues and eigenvectors of atomic-level stress matrices

Considerations of local atomic-level stresses associated with each atom represent a particular approach to address structures of disordered materials at the atomic level. We studied structural correlations in a two-dimensional model liquid using molecular dynamics simulations in the following way. We diagonalized the atomic-level stress tensor of every atom and investigated correlations between the eigenvalues and orientations of the eigenvectors of different atoms as a function of distance between them. It is demonstrated that the suggested approach can be used to characterize structural correlations in disordered materials. In particular, we found that changes in the stress correlation functions on decrease of temperature are the most pronounced for the pairs of atoms with separation distance that corresponds to the first minimum in the pair density function. We also show that the angular dependencies of the stress correlation functions previously reported by Wu et al. [Phys. Rev. E 91, 032301 (2015)10.1103/PhysRevE.91.032301] do not represent the anisotropic Eshelby's stress fields, as it is suggested, but originate in the rotational properties of the stress tensors.

Nonlinear plastic modes in disordered solids

We propose a theoretical framework within which a robust micromechanical definition of precursors to plastic instabilities, often termed soft spots, naturally emerges. They are shown to be collective displacements (modes) z[over ̂] that correspond to local minima of a barrier function b(z[over ̂]), which depends solely on inherent structure information. We demonstrate how some heuristic searches for local minima of b(z[over ̂]) can a priori detect the locus and geometry of imminent plastic instabilities with remarkable accuracy, at strains as large as γ_{c}-γ∼10^{-2} away from the instability strain γ_{c}. Our findings suggest that the a priori detection of the entire field of soft spots can be effectively carried out by a systematic investigation of the landscape of b(z[over ̂]).

The effect of a reversible shear transformation on plastic deformation of an amorphous solid

Langevin dynamics simulations are performed to investigate the plastic response of a model glass to a local shear transformation in a quiescent system. The deformation of the material is induced by a spherical inclusion that is gradually strained into an ellipsoid of the same volume and then reverted back into the sphere. We show that the number of cage-breaking events increases with increasing strain amplitude of the shear transformation. The results of numerical simulations indicate that the density of cage jumps is larger in the cases of weak damping or slow shear transformation. Remarkably, we also found that, for a given strain amplitude, the peak value of the density profiles is a function of the ratio of the damping coefficient and the time scale of the shear transformation.

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.

Plastic deformation of a model glass induced by a local shear transformation

The effect of a local shear transformation on plastic deformation of a three-dimensional amorphous solid is studied using molecular dynamics simulations. We consider a spherical inclusion, which is gradually transformed into an ellipsoid of the same volume and then converted back into the sphere. It is shown that at sufficiently large strain amplitudes, the deformation of the material involves localized plastic events that are identified based on the relative displacement of atoms before and after the shear transformation. We find that the density profiles of cage jumps decay away from the inclusion, which correlates well with the radial dependence of the local deformation of the material. At the same strain amplitude, the plastic deformation becomes more pronounced in the cases of weakly damped dynamics or large time scales of the shear transformation. We show that the density profiles can be characterized by the universal function of the radial distance multiplied by a dimensionless factor that depends on the friction coefficient and the time scale of the shear event.

Anisotropic stress correlations in two-dimensional liquids

In this paper we demonstrate the presence of anisotropic stress correlations in the simulated two-dimensional liquids. Whereas the temporal correlation of macroscopic shear stress is known to contribute to viscosity via the Green-Kubo formula, the general question regarding angular dependence of the spatial correlation among atomic-level stresses in liquids without external shear has not been explored. We observed the apparent anisotropicity with well-defined symmetry which can be explained in terms of the elastic continuum theory by Eshelby. In addition, we found that the shear stress correlation is screened compared to the prediction by the elastic continuum theory, and the screening length depends on temperature and follows the power law, suggesting divergence around the glass transition temperature. The success of the Eshelby theory to explain the anisotropy of the stress correlations justifies the idea that the mismatch between the atom and its nearest neighbor cage produces the atomic-level stress as well as the long-range stress fields.

Spatiotemporal correlations between plastic events in the shear flow of athermal amorphous solids

The slow flow of amorphous solids exhibits striking heterogeneities: swift localised particle rearrangements take place in the midst of a more or less homogeneously deforming medium. Recently, experimental as well as numerical work has revealed spatial correlations between these flow heterogeneities. Here, we use molecular dynamics (MD) simulations to characterise the rearrangements and systematically probe their correlations both in time and in space. In particular, these correlations display a four-fold azimuthal symmetry characteristic of shear stress redistribution in an elastic medium and we unambiguously detect their increase in range with time. With increasing shear rate, correlations become shorter-ranged. In addition, we study a coarse-grained model motivated by the observed flow characteristics and challenge its predictions directly with the MD simulations. While the model captures both macroscopic and local properties rather satisfactorily, the agreement with respect to the spatiotemporal correlations is at most qualitative. The discrepancies provide important insight into relevant physics that is missing in all related coarse-grained models that have been developed for the flow of amorphous materials so far, namely the finite shear wave velocity and the impact of elastic heterogeneities on stress redistribution.