Mechanical failure in amorphous solids: Scale-free spinodal criticality

Brittle and ductile yielding in soft materials

Many soft materials yield under mechanical loading, but how this transition from solid-like behavior to liquid-like behavior occurs can vary significantly. Understanding the physics of yielding is of great interest for the behavior of biological, environmental, and industrial materials, including those used as inks in additive manufacturing and muds and soils. For some materials, the yielding transition is gradual, while others yield abruptly. We refer to these behaviors as being ductile and brittle. The key rheological signatures of brittle yielding include a stress overshoot in steady-shear-startup tests and a steep increase in the loss modulus during oscillatory amplitude sweeps. In this work, we show how this spectrum of yielding behaviors may be accounted for in a continuum model for yield stress materials by introducing a parameter we call the brittility factor. Physically, an increased brittility decreases the contribution of recoverable deformation to plastic deformation, which impacts the rate at which yielding occurs. The model predictions are successfully compared to results of different rheological protocols from a number of real yield stress fluids with different microstructures, indicating the general applicability of the phenomenon of brittility. Our study shows that the brittility of soft materials plays a critical role in determining the rate of the yielding transition and provides a simple tool for understanding its effects under various loading conditions.

Annealing glasses by cyclic shear deformation

A major challenge in simulating glassy systems is the ability to generate configurations that may be found in equilibrium at sufficiently low temperatures, in order to probe static and dynamic behavior close to the glass transition. A variety of approaches have recently explored ways of surmounting this obstacle. Here, we explore the possibility of employing mechanical agitation, in the form of cyclic shear deformation, to generate low energy configurations in a model glass former. We perform shear deformation simulations over a range of temperatures, shear rates, and strain amplitudes. We find that shear deformation induces faster relaxation toward low energy configurations, or overaging, in simulations at sufficiently low temperatures, consistently with previous results for athermal shear. However, for temperatures at which simulations can be run until a steady state is reached with or without shear deformation, we find that the inclusion of shear deformation does not result in any speed up of the relaxation toward low energy configurations. Although we find the configurations from shear simulations to have properties indistinguishable from an equilibrium ensemble, the cyclic shear procedure does not guarantee that we generate an equilibrium ensemble at a desired temperature. In order to ensure equilibrium sampling, we develop a hybrid Monte Carlo algorithm that employs cyclic shear as a trial generation step and has acceptance probabilities that depend not only on the change in internal energy but also on the heat dissipated (equivalently, work done). We show that such an algorithm, indeed, generates an equilibrium ensemble.

Brittle yielding in supercooled liquids below the critical temperature of mode coupling theory

Molecular dynamics computer simulations of a polydisperse soft-sphere model under shear are presented. The starting point for these simulations are deeply supercooled samples far below the critical temperature, T c, of mode coupling theory. These samples are fully equilibrated with the aid of the swap Monte Carlo technique. For states below T c, we identify a lifetime τlt that measures the time scale on which the system can be considered as an amorphous solid. The temperature dependence of τlt can be well described by an Arrhenius law. The existence of transient amorphous solid states below T c is associated with the possibility of brittle yielding, as manifested by a sharp stress drop in the stress–strain relation and shear banding. We show that brittle yielding requires, on the one hand, low shear rates and, on the other hand, the time scale corresponding to the inverse shear rate has to be smaller or of the order of τlt. Both conditions can only be met for a large lifetime τlt, i.e., for states far below T c.

Universality of plastic instability and mechanical yield in metallic glasses

The generic response of a wide range of amorphous solids is the average increase of stress upon external loading until the yielding transition point, after which elasto-plastic steady state sets in. The stress-strain response comprises of a series of elastic branches interspersed with plastic drops. The ubiquitousness of these phenomena indicates universality, independent of material property, but the literature predominantly deals with specific materials. In pursuit of generality among different amorphous systems, we undertake a careful investigation in the mechanical response of metallic glasses using computer simulation. By comparing our results of multi-body metallic glass potentials to those obtained from pairwise Lennard-Jones glasses, we show that the mechanism of plastic instabilities is universal and independent of the details of the underlying potential. We also investigate the yielding transition in terms of the overlap parameterQ₁₂, which has been successfully used Lennard-Jones glasses. The yielding is unambiguously identified as a first-order phase transition. These observations conform the nature of plastic instabilities and mechanical yield as universal and independent of microscopic interactions.

The role of annealing in determining the yielding behavior of glasses under cyclic shear deformation

Yielding behavior in amorphous solids has been investigated in computer simulations using uniform and cyclic shear deformation. Recent results characterize yielding as a discontinuous transition, with the degree of annealing of glasses being a significant parameter. Under uniform shear, discontinuous changes in stresses at yielding occur in the high annealing regime, separated from the poor annealing regime in which yielding is gradual. In cyclic shear simulations, relatively poorly annealed glasses become progressively better annealed as the yielding point is approached, with a relatively modest but clear discontinuous change at yielding. To understand better the role of annealing on yielding characteristics, we perform athermal quasistatic cyclic shear simulations of glasses prepared with a wide range of annealing in two qualitatively different systems-a model of silica (a network glass) and an atomic binary mixture glass. Two strikingly different regimes of behavior emerge. Energies of poorly annealed samples evolve toward a unique threshold energy as the strain amplitude increases, before yielding takes place. Well-annealed samples, in contrast, show no significant energy change with strain amplitude until they yield, accompanied by discontinuous energy changes that increase with the degree of annealing. Significantly, the threshold energy for both systems corresponds to dynamical cross-over temperatures associated with changes in the character of the energy landscape sampled by glass-forming liquids.

Ductile and Brittle Yielding in Thermal and Athermal Amorphous Materials

We study theoretically the yielding of sheared amorphous materials as a function of increasing levels of initial sample annealing prior to shear, in three widely used constitutive models and three widely studied annealing protocols. In thermal systems we find a gradual progression, with increasing annealing, from smoothly "ductile" yielding, in which the sample remains homogeneous, to abruptly "brittle" yielding, in which it becomes strongly shear banded. This progression arises from an increase with annealing in the size of an overshoot in the underlying stress-strain curve for homogeneous shear, which causes a shear banding instability that becomes more severe with increasing annealing. Ductile and brittle yielding thereby emerge as two limiting cases of a continuum of yielding transitions, from gradual to catastrophic. In contrast, athermal systems with a stress overshoot always show brittle yielding at low shear rates, however small the overshoot.

Flow heterogeneities in supercooled liquids and glasses under shear

Using extensive nonequilibrium molecular dynamics simulations, we investigate a glass-forming binary Lennard-Jones mixture under shear. Both supercooled liquids and glasses are considered. Our focus is on the characterization of inhomogeneous flow patterns such as shear bands that appear as a transient response to the external shear. For the supercooled liquids, we analyze the crossover from Newtonian to non-Newtonian behavior with increasing shear rate γ[over ̇]. Above a critical shear rate γ[over ̇]_{c} where a non-Newtonian response sets in, the transient dynamics are associated with the occurrence of short-lived vertical shear bands, i.e., bands of high mobility that form perpendicular to the flow direction. In the glass states, long-lived horizontal shear bands, i.e., bands of high mobility parallel to the flow direction, are observed in addition to vertical ones. The systems with shear bands are characterized in terms of mobility maps, stress-strain relations, mean-squared displacements, and (local) potential energies. The initial formation of a horizontal shear band provides an efficient stress release, corresponds to a local minimum of the potential energy, and is followed by a slow broadening of the band towards the homogeneously flowing fluid in the steady state. Whether a horizontal or a vertical shear band forms cannot be predicted from the initial undeformed sample. Furthermore, we show that with increasing system size, the probability for the occurrence of horizontal shear bands increases.

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.

Renormalization of elastic quadrupoles in amorphous solids

Plasticity in amorphous solids is mediated by localized quadrupolar instabilities, but the mechanism by which an amorphous solid eventually fails or melts is debated. In this work we argue that these phenomena can be investigated in the model problem of an elastic continuum with quadrupolar defects, at finite temperature. This problem is posed and the collective behavior of the defects is analytically investigated. Using both renormalization group and field-theoretic techniques, it is found that the model has a yielding/melting transition of spinodal type.

Critical Slowing Down at the Abrupt Mott Transition: When the First-Order Phase Transition Becomes Zeroth Order and Looks Like Second Order

We report that the thermally induced Mott transition in vanadium sesquioxide shows critical slowing down and enhanced variance ("critical opalescence") of the order parameter fluctuations measured through low-frequency resistance-noise spectroscopy. Coupled with the observed increase of the phase-ordering time, these features suggest that the strong abrupt transition is controlled by a critical-like singularity in the hysteretic metastable phase. The singularity is identified with the spinodal point and is a likely consequence of the strain-induced long-range interaction.

Rigidity-Controlled Crossover: From Spinodal to Critical Failure

Failure in disordered solids is accompanied by intermittent fluctuations extending over a broad range of scales. The implied scaling has been previously associated with either spinodal or critical points. We use an analytically transparent mean-field model to show that both analogies are relevant near the brittle-to-ductile transition. Our study indicates that in addition to the strength of quenched disorder, an appropriately chosen global measure of rigidity (connectivity) can be also used to tune the system to criticality. By interpreting rigidity as a timelike variable we reveal an intriguing parallel between earthquake-type critical failure and Burgers turbulence.

Particle pinning suppresses spinodal criticality in the shear-banding instability

Strained amorphous solids often fail mechanically by creating a shear band. It had been understood that the shear-banding instability is usefully described as crossing a spinodal point (with disorder) in an appropriate thermodynamic description. It remained contested, however, whether the spinodal is critical (with divergent correlation length) or not. Here we offer evidence for critical spinodal by using particle pinning. For a finite concentration of pinned particles the correlation length is bounded by the average distance between pinned particles, but without pinning it is bounded by the system size.

Effect of Pinning on the Yielding Transition of Amorphous Solids

Using numerical simulations, we have studied the yielding response, in the athermal quasistatic limit, of a model amorphous material having inclusions in the form of randomly pinned particles. We show that, with increasing pinning concentration, the plastic activity becomes more spatially localized, resulting in smaller stress drops, and a corresponding increase in the magnitude of strain where yielding occurs. We demonstrate that, unlike the spatially heterogeneous and avalanche led yielding in the case of the unpinned glass, for the case of large pinning concentration, yielding takes place via a spatially homogeneous proliferation of localized events.

Mean-field stability map of hard-sphere glasses

The response of amorphous solids to an applied shear deformation is an important problem, both in fundamental and in applied research. To tackle this problem, we focus on a system of hard spheres in infinite dimensions as a solvable model for colloidal systems and granular matter. The system is prepared above the dynamical glass transition density, and we discuss the phase diagram of the resulting glass under compression, decompression, and shear strain, expanding on previous results [Urbani and Zamponi, Phys. Rev. Lett. 118, 038001 (2017)PRLTAO0031-900710.1103/PhysRevLett.118.038001]. We show that the solid region is bounded by a "shear jamming" line, at which the solid reaches close packing, and a "shear yielding" line, at which the solid undergoes a spinodal instability towards a liquid flowing phase. Furthermore, we characterize the evolution of these lines upon varying the glass preparation density. This paper aims to provide a general overview on yielding and jamming phenomena in hard-sphere systems by a systematic exploration of the phase diagram.

Microscopic Theory of Two-Step Yielding in Attractive Colloids

Attractive colloids display two distinct amorphous solid phases: the attractive glass due to particle bonding and the repulsive glass due to the hard-core repulsion. By means of a microscopic mean field approach, we analyze their response to a quasistatic shear strain. We find that the presence of two distinct interaction length scales may result in a sharp two-step yielding process, which can be associated with a hysteretic stress response or with a reversible but nonmonotonic stress-strain curve. We derive a generic phase diagram characterized by two distinct yielding lines, an inverse yielding and a critical point separating the hysteretic and reversible regimes. Our results should be applicable to a large class of glassy materials characterized by two distinct interaction length scales.

Random critical point separates brittle and ductile yielding transitions in amorphous materials

We combine an analytically solvable mean-field elasto-plastic model with molecular dynamics simulations of a generic glass former to demonstrate that, depending on their preparation protocol, amorphous materials can yield in two qualitatively distinct ways. We show that well-annealed systems yield in a discontinuous brittle way, as metallic and molecular glasses do. Yielding corresponds in this case to a first-order nonequilibrium phase transition. As the degree of annealing decreases, the first-order character becomes weaker and the transition terminates in a second-order critical point in the universality class of an Ising model in a random field. For even more poorly annealed systems, yielding becomes a smooth crossover, representative of the ductile rheological behavior generically observed in foams, emulsions, and colloidal glasses. Our results show that the variety of yielding behaviors found in amorphous materials does not necessarily result from the diversity of particle interactions or microscopic dynamics but is instead unified by carefully considering the role of the initial stability of the system.

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.