Search

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

The Eshelby problem refers to the response of a two-dimensional elastic sheet to cutting away a circle, deforming it into an ellipse, and pushing it back. The resulting response is dominated by the so-called Eshelby kernel, which was derived for purely elastic (infinite) material, but has been employed extensively to model the redistribution of stress after plastic events in amorphous solids with finite boundaries. Here, we discuss and solve the Eshelby problem directly for amorphous solids, taking into account possible screening effects and realistic boundary conditions. We find major modifications compared to the classical Eshelby solution. These modifications are needed for modeling correctly the spatial responses to plastic events in amorphous solids.

The response of amorphous solids to mechanical loads is accompanied by plasticity that is generically associated with "non-affine" quadrupolar events seen in the resulting displacement field. To develop a continuum theory, one needs to assess when these quadrupolar events have a finite density, allowing the development of a field theory. Is there a transition, as a function of the material parameters and the nature of the loads, from isolated plastic events whose density is zero to a regime governed by a finite density? And if so, what is the nature of this transition? The aim of the paper is to explore this issue. The motivation for the present study stems from recent research in which it was shown that gradients of the quadrupolar fields act as dipole charges that can screen elasticity. Analytically soluble examples of mechanical loading that lead to screening and emergent length scales (that are absent in classical elasticity) have been analyzed and tested. However, "gradients of quadrupolar fields" make sense only when the density of quadrupoles is finite, and hence, the issue is central to this article. The article introduces a notion of polarizability under the strain of Eshelby quadrupoles and concludes that the onset of a density of such quadrupoles with random orientations can only appear when the polarizability is finite., (© 2024 Author(s). Published under an exclusive license by AIP Publishing.)

In recent work it was shown that elasticity theory can break down in amorphous solids subjected to nonuniform static loads. The elastic fields are screened by geometric dipoles; these stem from gradients of the quadrupole field associated with plastic responses. Here we study the dynamical responses induced by oscillatory loads. The required modification to classical elasticity is described. Exact solutions for the displacement field in circular geometry are presented, demonstrating that dipole screening results in essential departures from the expected predictions of classical elasticity theory. Numerical simulations are conducted to validate the theoretical predictions and to delineate their range of validity.

A significant amount of attention was dedicated in recent years to the phenomenon of jamming of athermal amorphous solids by increasing the volume fraction of the microscopic constituents. At a critical value of the volume fraction, pressure shoots up from zero to finite values with a host of critical exponents discovered and discussed. In this paper, we advance evidence for the existence of a second transition, within the jammed state of two-dimensional granular systems, that separates two regimes of different mechanical responses. Explicitly, highly packed systems are quasielastic with quadrupole screening, and more loosely jammed systems exhibit anomalous mechanics with dipole screening. Evidence is given for a clear transition between these two regimes, reminiscent of the intermediate hexatic phase of crystal melting in two-dimensional crystals. Theoretical estimates of the screening parameters and the pressure where transition takes place are provided.

When amorphous solids are subjected to simple or pure strain, they exhibit elastic increase in stress, punctuated by plastic events that become denser (in strain) upon increasing the system size. It is customary to assume in theoretical models that the stress released in each plastic event is redistributed according to the linear Eshelby kernel, causing avalanches of additional stress release. Here we demonstrate that, contrary to the uniform affine strain resulting from simple or pure strain, each plastic event is associated with a nonuniform strain that gives rise to a displacement field that contains quadrupolar and dipolar charges that typically screen the linear elastic phenomenology and introduce anomalous length scales and influence the form of the stress redistribution. An important question that opens up is how to take this into account in elastoplastic models of shear induced phenomena like shear banding.

In recent work, we developed a screening theory for describing the effect of plastic events in amorphous solids on its emergent mechanics. The suggested theory uncovered an anomalous mechanical response of amorphous solids where plastic events collectively induce distributed dipoles that are analogous to dislocations in crystalline solids. The theory was tested against various models of amorphous solids in two dimensions, including frictional and frictionless granular media and numerical models of amorphous glass. Here we extend our theory to screening in three-dimensional amorphous solids and predict the existence of anomalous mechanics similar to the one observed in two-dimensional systems. We conclude by interpreting the mechanical response as the formation of nontopological distributed dipoles that have no analog in the crystalline defects literature. Having in mind that the onset of dipole screening is reminiscent of Kosterlitz-Thouless and hexatic transitions, the finding of dipole screening in three dimensions is surprising.

Applying constant tensile stress to a piece of amorphous solid results in a slow extension, followed by an eventual rapid mechanical collapse. This "creep" process is of paramount engineering concern, and as such was the subject of study in a variety of materials, for more than a century. Predictive theories for τ_{w}, the expected time of collapse, are incomplete, mainly due to its dependence on a bewildering variety of parameters, including temperature, system size, tensile force, but also the detailed microscopic interactions between constituents. The complex dependence of the collapse time on all the parameters is discussed below, using simulations of strip of amorphous material. Different scenarios are observed for ductile and brittle materials, resulting in serious difficulties in creating an all-encompassing theory that could offer safety measures for given conditions. A central aim of this paper is to employ scaling concepts, to achieve data collapse for the probability distribution function (pdf) of lnτ_{w}. The scaling ideas result in a universal function which provides a prediction of the pdf of lnτ_{w} for out-of-sample systems, from measurements at other values of these parameters. The predictive power of the scaling theory is demonstrated for both ductile and brittle systems. Finally, we present a derivation of universal scaling function for brittle materials. The ductile case appears to be due to a plastic necking instability and is left for future research.

Amorphous solids under mechanical strains are prone to plastic responses. Recent work showed that in amorphous granular systems these plastic events, that are typically quadrupolar in nature, can screen the elastic response. When the density of the quadrupoles is high, the gradients of the quadrupole field act as emergent dipole sources, leading to qualitative changes in the mechanical response, as seen for example in the displacement field. In this paper we examine the effect of screening in classical glass formers. These are made of point particles that interact via binary forces. Both inverse power law forces and Lennard-Jones interactions are examined, and it is shown that in both cases the elastic response can be strongly screened, in agreement with the novel theory. The degree of deviation from classical elasticity theory is parametrized by a proposed measure that is shown to have a functional dependence on the amount of energy lost to plastic responses.

The theoretical understanding of the low-frequency modes in amorphous solids at finite temperature is still incomplete. The study of the relevant modes is obscured by the dressing of interparticle forces by collision-induced momentum transfer that is unavoidable at finite temperatures. Recently, it was proposed that low-frequency modes of vibrations around the thermally averaged configurations deserve special attention. In simple model glasses with bare binary interactions, these included quasilocalized modes whose density of states appears to be universal, depending on the frequencies as D(ω)∼ω^{4}, in agreement with the similar law that is obtained with bare forces at zero temperature. In this paper, we report investigations of a model of silica glass at finite temperature; here the bare forces include binary and ternary interactions. Nevertheless, we can establish the validity of the universal law of the density of quasilocalized modes also in this richer and more realistic model glass.

Recent progress in studying the physics of amorphous solids has revealed that mechanical strains can be strongly screened by the formation of plastic events that are typically quadrupolar in nature. The theory stipulates that gradients in the density of the quadrupoles act as emergent dipole sources, leading to strong screening and to qualitative changes in the mechanical response, as seen, for example, in the displacement field. In this Letter we first offer direct measurements of the dipole field, independently of any theoretical assumptions, and second we demonstrate detailed agreement with the recently proposed theory. These two goals are achieved by using data from both simulations and experiments. Finally, we show how measurements of the dipole fields pinpoint the theory parameters that determine the profile of the displacement field.

In mechanical engineering Wöhler plots serve to measure the average number of load cycles before materials break, as a function of the maximal stress in each cycle. Although such plots have been prevalent in engineering for more than 150 years, their theoretical understanding is lacking. Recently a scaling theory of Wöhler plots in the context of cyclic bending was offered [Bhowmik, arXiv:2103.03040 (2021)]. Here we elaborate further on cyclic bending and extend the considerations to cyclic tensile loads on an amorphous strip of material; the scaling theory applies to both types of cyclic loading equally well. On the basis of atomistic simulations we conclude that the crucial quantities to focus on are the accumulated damage and the average damage per cycle. The dependence of these quantities on the loading determines the statistics of the number of cycles to failure. Finally, we consider the probability distribution functions of the number of cycles to failure and demonstrate that the scaling theory allows prediction of these distributions at one value of the forcing amplitude from measurements and another value.

"Remote triggering" refers to the inducement of earthquakes by weak perturbations that emanate from faraway sources, typically intense earthquakes that happen at much larger distances than their nearby aftershocks, sometimes even around the globe. Here, we propose a mechanism for this phenomenon; the proposed mechanism is generic, resulting from the breaking of Hamiltonian symmetry due to the existence of friction. We allow a transition from static to dynamic friction. Linearly stable stressed systems display giant sensitivity to small perturbations of arbitrary frequency (without a need for resonance), which trigger an instability with exponential oscillatory growth. Once nonlinear effects kick in, the blow up in mean-square displacements can reach 15-20 orders of magnitude. Analytic and numerical results of the proposed model are presented and discussed.

Amorphous solids appear to react elastically to small external strains, but in contrast to ideal elastic media, plastic responses abound immediately at any value of the strain. Such plastic responses are quasilocalized in nature, with the "cheapest" one being a quadrupolar source. The existence of such plastic responses results in screened elasticity in which strains and stresses can either quantitatively or qualitatively differ from the unscreened theory, depending on the specific screening mechanism. Here we offer a theory of such screening effects by plastic quadrupoles, dipoles, and monopoles, explain their natural appearance, and point out the analogy to electrostatic screening by electric charges and dipoles. For low density of quadrupoles the effect is to normalize the elastic moduli without a qualitative change compared to pure elasticity theory; for higher density of quadrupoles the screening effects result in qualitative changes. Predictions for the spatial dependence of displacement fields caused by local sources of strains are provided and compared to numerical simulations. We find that anomalous elasticity is richer than electrostatics in having a screening mode that does not appear in the electrostatic analog.

It has been established that the low frequency quasilocalized modes of amorphous solids at zero temperature exhibit universal density of states, depending on the frequencies as D(ω)∼ω^{4}. It remains an open question whether this universal law extends to finite temperatures. In this Letter we show that well quenched model glasses at temperatures as high as T_{g}/3 possess the same universal density of states. The only condition required is that average particle positions stabilize before thermal diffusion destroys the cage structure of the material. The universal density of quasilocalized low frequency modes refers then to vibrations around the thermally averaged configuration of the material.

The rupture of a polymer chain maintained at temperature T under fixed tension is prototypical to a wide array of systems failing under constant external stress and random perturbations. Past research focused on analytic and numerical studies of the mean rate of collapse of such a chain. Surprisingly, an analytic calculation of the probability distribution function (PDF) of collapse rates appears to be lacking. Since rare events of rapid collapse can be important and even catastrophic, we present here a theory of this distribution, with a stress on its tail of fast rates. We show that the tail of the PDF is a power law with a universal exponent that is theoretically determined. Extensive numerics validate the offered theory. Lessons pertaining to other problems of the same type are drawn.

Theoretical treatments of frictional granular matter often assume that it is legitimate to invoke classical elastic theory to describe its coarse-grained mechanical properties. Here, we show, based on experiments and numerical simulations, that this is generically not the case since stress autocorrelation functions decay more slowly than what is expected from elasticity theory. It was theoretically shown that standard elastic decay demands pressure and torque density fluctuations to be normal, with possibly one of them being hyperuniform. However, generic compressed frictional assemblies exhibit abnormal pressure fluctuations, failing to conform with the central limit theorem. The physics of this failure is linked to correlations built in the material during compression from a dilute configuration prior to jamming. By changing the protocol of compression, one can observe different pressure fluctuations, and stress autocorrelations decay at large scales.

Diffusion-limited aggregation (DLA) has served for 40 years as a paradigmatic example for the creation of fractal growth patterns. In spite of thousands of references, no exact result for the fractal dimension D of DLA is known. In this Letter we announce an exact result for off-lattice DLA grown on a line embedded in the plane D=3/2. The result relies on representing DLA with iterated conformal maps, allowing one to prove self-affinity, a proper scaling limit, and a well-defined fractal dimension. Mathematical proofs of the main results are available in N. Berger, E. B. Procaccia, and A. Turner, Growth of stationary Hastings-Levitov, arXiv:2008.05792.

It was recently shown that different simple models of glass formers with binary interactions define a universality class in terms of the density of states of their quasilocalized low-frequency modes. Explicitly, once the hybridization with standard Debye (extended) modes is avoided, a number of such models exhibit a universal density of states, depending on the mode frequencies as D(ω)∼ω^{4}. It is unknown, however, how wide this universality class is, and whether it also pertains to more realistic models of glass formers. To address this issue we present analysis of the quasilocalized modes in silica, a network glass that has both binary and ternary interactions. We conclude that in three dimensions silica exhibits the very same frequency dependence at low frequencies, suggesting that this universal form is a generic consequence of amorphous glassiness.

The extreme slowing down associated with glass formation in experiments and in simulations results in serious difficulties to prepare deeply quenched, well annealed, glassy material. Recently, methods to achieve such deep quenching were proposed, including vapor deposition on the experimental side and "swap Monte Carlo" and oscillatory shearing on the simulation side. The relation between the resulting glasses under different protocols remains unclear. Here we show that oscillatory shear and swap Monte Carlo result in thermodynamically equivalent glasses sharing the same statistical mechanics and similar mechanical responses under external strain.

Catastrophic events in nature can be often triggered by small perturbations, with "remote triggering" of earthquakes being an important example. Here we present a mechanism for the giant amplification of small perturbations that is expected to be generic in systems whose dynamics is not derivable from a Hamiltonian. We offer a general discussion of the typical instabilities involved (being oscillatory with an exponential increase of noise) and examine in detail the normal forms that determine the relevant dynamics. The high sensitivity to external perturbations is explained for systems with and without dissipation. Numerical examples are provided using the dynamics of frictional granular matter. Finally, we point out the relationship of the presently discussed phenomenon to the highly topical issue of "exceptional points" in quantum models with non-Hermitian Hamiltonians.

The dynamics of amorphous granular matter with frictional interactions cannot be derived in general from a Hamiltonian and therefore displays oscillatory instabilities stemming from the onset of complex eigenvalues in the stability matrix. These instabilities were discovered in the context of one- and two-dimensional systems, while the three-dimensional case was never studied in detail. Here we fill this gap by deriving and demonstrating the presence of oscillatory instabilities in a three-dimensional granular packing. We study binary assemblies of spheres of two sizes interacting via classical Hertz and Mindlin force laws for the longitudinal and tangent interactions, respectively. We formulate analytically the stability matrix in three dimensions and observe that a couple of complex eigenvalues emerge at the onset of the instability as in the case of frictional disks in two dimensions. The dynamics then shows oscillatory exponential growth in the mean-square displacement, followed by a catastrophic event in which macroscopic portions of mechanical stress and energy are lost. The generality of these results for any choice of forces that break the symplectic Hamiltonian symmetry is discussed.

Electrostatic theory preserves charges, but allows dipolar excitations. Elasticity theory preserves dipoles, but allows quadrupolar (Eshelby-like) plastic events. Charged amorphous granular systems are interesting in their own right; here we focus on their plastic instabilities and examine their mechanical response to external strain and to an external electric field, to expose the competition between elasticity and electrostatics. In this paper a generic model is offered, its mechanical instabilities are examined, and a theoretical analysis is presented. Plastic instabilities are discussed as saddle-node bifurcations that can be fully understood in terms of eigenvalues and eigenfunctions of the relevant Hessian matrix. This system exhibits moduli that describe how electric polarization and stress are influenced by strain and the electric field. Theoretical expression for these moduli are offered and compared to the measurements in numerical simulations.

The nature of an instability that controls the transition from static to dynamical friction is studied in the context of an array of frictional disks that are pressed from above on a substrate. In this case the forces are all explicit and Newtonian dynamics can be employed without any phenomenological assumptions. We show that an oscillatory instability that had been discovered recently is responsible for the transition, allowing individual disks to spontaneously reach the Coulomb limit and slide with dynamic friction. The transparency of the model allows a full understanding of the phenomenon, including the speeds of the waves that travel from the trailing to the leading edge and vice versa.

In amorphous solids at finite temperatures the particles follow chaotic trajectories which, at temperatures sufficiently lower than the glass transition, are trapped in "cages." Averaging their positions for times shorter than the diffusion time, one can define a time-averaged configuration. Under strain or stress, these average configurations undergo sharp plastic instabilities. In athermal glasses the understanding of plastic instabilities is furnished by the Hessian matrix and its eigenvalues and eigenfunctions. Here we propose an uplifting of Hessian methods to thermal glasses, with the aim of understanding the plastic responses in the time-averaged configuration. We discuss a number of potential definitions of Hessians and identify which of these can provide eigenvalues and eigenfunctions which can explain and predict the instabilities of the time-averaged configurations. The conclusion is that the nonaffine changes in the average configurations during an instability is accurately predicted by the eigenfunctions of the low-lying eigenvalues of the chosen Hessian.

Plastic instabilities in amorphous materials are often studied using idealized models of binary mixtures that do not capture accurately molecular interactions and bonding present in real glasses. Here we study atomic-scale plastic instabilities in a three-dimensional molecular dynamics model of silica glass under quasistatic shear. We identify two distinct types of elementary plastic events, one is a standard quasilocalized atomic rearrangement while the second is a bond-breaking event that is absent in simplified models of fragile glass formers. Our results show that both plastic events can be predicted by a drop of the lowest nonzero eigenvalue of the Hessian matrix that vanishes at a critical strain. Remarkably, we find very high correlation between the associated eigenvectors and the nonaffine displacement fields accompanying the bond-breaking event, predicting the locus of structural failure. Both eigenvectors and nonaffine displacement fields display an Eshelby-like quadrupolar structure for both failure modes, rearrangement, and bond breaking. Our results thus clarify the nature of atomic-scale plastic instabilities in silica glasses, providing useful information for the development of mesoscale models of amorphous plasticity.

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.

It was discovered recently that frictional granular materials can exhibit an important mechanism for instabilities, i.e., the appearance of pairs of complex eigenvalues in their stability matrix. The consequence is an oscillatory exponential growth of small perturbations which are tamed by dynamical nonlinearities. The amplification can be giant, many orders of magnitude, and it ends with a catastrophic system-spanning plastic event. Here we follow up on this discovery, explore the scaling laws characterizing the onset of the instability, the scenarios of the development of the instability with and without damping, and the nature of the eventual system-spanning events. The possible relevance to earthquake physics and to the transition from static to dynamic friction is discussed.

We study agitated frictional disks in two dimensions with the aim of developing a scaling theory for their diffusion over time. As a function of the area fraction ϕ and mean-square velocity fluctuations 〈v^{2}〉 the mean-square displacement of the disks 〈d^{2}〉 spans four to five orders of magnitude. The motion evolves from a subdiffusive form to a complex diffusive behavior at long times. The statistics of 〈d^{n}〉 at all times are multiscaling, since the probability distribution function (PDF) of displacements has very broad wings. Even where a diffusion constant can be identified it is a complex function of ϕ and 〈v^{2}〉. By identifying the relevant length and time scales and their interdependence one can rescale the data for the mean-square displacement and the PDF of displacements into collapsed scaling functions for all ϕ and 〈v^{2}〉. These scaling functions provide a predictive tool, allowing one to infer from one set of measurements (at a given ϕ and 〈v^{2}〉) what are the expected results at any value of ϕ and 〈v^{2}〉 within the scaling range.

Frictional granular matter is shown to be fundamentally different in its plastic responses to external strains from generic glasses and amorphous solids without friction. While regular glasses exhibit plastic instabilities due to the vanishing of a real eigenvalue of the Hessian matrix, frictional granular materials can exhibit a previously unnoticed additional mechanism for instabilities, i.e., the appearance of a pair of complex eigenvalues leading to oscillatory exponential growth of perturbations that are tamed by dynamical nonlinearities. This fundamental difference appears crucial for the understanding of plasticity and failure in frictional granular materials. The possible relevance to earthquake physics is discussed.

Shearing with a finite shear rate a compressed granular system results in a region of grains flowing over a compact, static assembly. Perforce this region is dilated to a degree that depends on the shear rate, the loading pressure, gravity, various material parameters, and the preparation protocol. In spite of numerous studies of granular flows a predictive theory of the amount of dilation is still lacking. Here, we offer a scaling theory that is focused on such a prediction as a function of shear rate and the dissipative parameters of the granular assembly. The resulting scaling laws are universal with respect to changing the interparticle force laws.

Three-dimensional anisotropic turbulence in classical fluids tends towards isotropy and homogeneity with decreasing scales, allowing-eventually-the abstract model of homogeneous and isotropic turbulence to be relevant. We show here that the opposite is true for superfluid ^{4}He turbulence in three-dimensional counterflow channel geometry. This flow becomes less isotropic upon decreasing scales, becoming eventually quasi-two-dimensional. The physical reason for this unusual phenomenon is elucidated and supported by theory and simulations.

In thermal glasses at temperatures sufficiently lower than the glass transition, the constituent particles are trapped in their cages for a sufficiently long time such that their time-averaged positions can be determined before diffusion and structural relaxation takes place. The effective forces are those that hold these average positions in place. In numerical simulations the effective forces F_{ij} between any pair of particles can be measured as a time average of the bare forces f_{ij}(r_{ij}(t)). In general, even if the bare forces come from two-body interactions, thermal dynamics dress the effective forces to contain many-body interactions. Here, we develop the effective theory for systems with generic interactions, where the effective forces are derivable from an effective potential and in turn they give rise to an effective Hessian whose eigenvalues are all positive when the system is stable. In this Rapid Communication, we offer analytic expressions for the effective theory, and demonstrate the usefulness and the predictive power of the approach.

We report a joint experimental and theoretical investigation of the probability distribution functions (PDFs) of the normal and tangential (frictional) forces in amorphous frictional media. We consider both the joint PDF of normal and tangential forces together, and the marginal PDFs of normal forces separately and tangential forces separately. A maximum entropy formalism is utilized for all these cases after identifying the appropriate constraints. Excellent agreements with both experimental and simulation data are reported. The proposed joint PDF predicts giant slip events at low pressures, again in agreement with observations.

In very recent work the mean field theory of the jamming transition in infinite-dimensional hard sphere models was presented. Surprisingly, this theory predicts quantitatively the numerically determined characteristics of jamming in two and three dimensions. This is a rare and unusual finding. Here we argue that this agreement is nongeneric: only for hard sphere models does it happen that sufficiently close to the jamming density (at any temperature) the effective interactions are binary, in agreement with mean field theory, justifying the truncation of many-body interactions (which is the exact protocol in infinite dimensions). Any softening of the bare hard sphere interactions results in many-body effective interactions that are not mean field at any density, making the d=∞ results not applicable.

The mechanical failure of amorphous media is a ubiquitous phenomenon from material engineering to geology. It has been noticed for a long time that the phenomenon is "scale-free," indicating some type of criticality. In spite of attempts to invoke "Self-Organized Criticality," the physical origin of this criticality, and also its universal nature, being quite insensitive to the nature of microscopic interactions, remained elusive. Recently we proposed that the precise nature of this critical behavior is manifested by a spinodal point of a thermodynamic phase transition. Demonstrating this requires the introduction of an "order parameter" that is suitable for distinguishing between disordered amorphous systems. At the spinodal point there exists a divergent correlation length which is associated with the system-spanning instabilities (known also as shear bands) which are typical to the mechanical yield. The theory, the order parameter used and the correlation functions which exhibit the divergent correlation length are universal in nature and can be applied to any amorphous solid that undergoes mechanical yield. The phenomenon is seen at its sharpest in athermal systems, as is explained below; in this paper we extend the discussion also to thermal systems, showing that at sufficiently high temperatures the spinodal phenomenon is destroyed by thermal fluctuations.

We revisit the problem of the stress distribution in a frictional sandpile with both normal and tangential (frictional) inter-granular forces, under gravity, equipped with a new numerical method of generating such assemblies. Numerical simulations allow a determination of the spatial dependence of all the components of the stress field, principle stress axis, angle of repose, as a function of systems size, the coefficient of static friction and the frictional interaction with the bottom surface. We compare these results with the predictions of a theory based on continuum equilibrium mechanics. Basic to the theory of sandpiles are assumptions about the form of scaling solutions and constitutive relations for cohesive-less hard grains for which no typical scale is available. We find that these constitutive relations must be modified; moreover for smaller friction coefficients and smaller piles these scaling assumptions break down in the bulk of the sandpile due to the presence of length scales that must be carefully identified. Fortunately, for larger friction coefficient and for larger piles the breaking of scaling is weak in the bulk, allowing an approximate analytic theory which agrees well with the observations. After identifying the crucial scale, triggering the breaking of scaling, we provide a predictive theory to when scaling solutions are expected to break down. At the bottom of the pile the scaling assumption breaks always, due to the different interactions with the bottom surface. The consequences for measurable quantities like the pressure distribution and shear stress at the bottom of the pile are discussed. For example one can have a transition from no dip in the base-pressure to a dip at the center of the pile as friction increases.

Amorphous solids increase their stress as a function of an applied strain until a mechanical yield point whereupon the stress cannot increase anymore, afterward exhibiting a steady state with a constant mean stress. In stress-controlled experiments, the system simply breaks when pushed beyond this mean stress. The ubiquity of this phenomenon over a huge variety of amorphous solids calls for a generic theory that is free of microscopic details. Here, we offer such a theory: The mechanical yield is a thermodynamic phase transition, where yield occurs as a spinodal phenomenon. At the spinodal point, there exists a divergent correlation length that is associated with the system-spanning instabilities (also known as shear bands), which are typical to the mechanical yield. The theory, the order parameter used, and the correlation functions that exhibit the divergent correlation length are universal in nature and can be applied to any amorphous solids that undergo mechanical yield., Competing Interests: The authors declare no conflict of interest.

In recent work it was clarified that amorphous solids under strain control do not possess nonlinear elastic theory in the sense that the shear modulus exists but nonlinear moduli exhibit sample-to-sample fluctuations that grow without bound with the system size. More relevant, however, for experiments are the conditions of stress control. In the present Rapid Communication we show that also under stress control the shear modulus exists, but higher-order moduli show unbounded sample-to-sample fluctuation. The unavoidable consequence is that the characterization of stress-strain curves in experiments should be done with a stress-dependent shear modulus rather than with nonlinear expansions.

Amorphous media at finite temperatures, be them liquids, colloids, or glasses, are made of interacting particles that move chaotically due to thermal energy, continuously colliding and scattering off each other. When the average configuration in these systems relaxes only at long times, one can introduce effective interactions that keep the mean positions in mechanical equilibrium. We introduce a framework to determine the effective force laws that define an effective Hessian that can be employed to discuss stability properties and the density of states of the amorphous system. We exemplify the approach with a thermal glass of hard spheres; these experience zero forces when not in contact and infinite forces when they touch. Close to jamming we recapture the effective interactions that at temperature T depend on the gap h between spheres as T/h [C. Brito and M. Wyart, Europhys. Lett. 76, 149 (2006)EULEEJ0295-507510.1209/epl/i2006-10238-x]. For hard spheres at lower densities or for systems whose binary bare interactions are longer ranged (at any density), the emergent force laws include ternary, quaternary, and generally higher-order many-body terms, leading to a temperature-dependent effective Hessian.

It is known [H. G. E. Hentschel et al., Phys. Rev. E 83, 061101 (2011)PLEEE81539-375510.1103/PhysRevE.83.061101] that amorphous solids at zero temperature do not possess a nonlinear elasticity theory: besides the shear modulus, which exists, none of the higher order coefficients exist in the thermodynamic limit. Here we show that the same phenomenon persists up to temperatures comparable to that of the glass transition. The zero-temperature mechanism due to the prevalence of dangerous plastic modes of the Hessian matrix is replaced by anomalous stress fluctuations that lead to the divergence of the variances of the higher order elastic coefficients. The conclusion is that in amorphous solids elasticity can never be decoupled from plasticity: the nonlinear response is very substantially plastic.

We consider the problem of how to determine the force laws in an amorphous system of interacting particles. Given the positions of the centers of mass of the constituent particles we propose an algorithm to determine the interparticle force laws. Having n different types of constituents we determine the coefficients in the Laurent polynomials for the n(n+1)/2 possibly different force laws. A visual providing the particle positions in addition to a measurement of the pressure is all that is required. The algorithm proposed includes a part that can correct for experimental errors in the positions of the particles. Such a correction of unavoidable measurement errors is expected to benefit many experiments in the field.

Plastic events in amorphous solids can be much more than just "shear transformation zones" when the positional degrees of freedom are coupled nontrivially to other degrees of freedom. Here we consider magnetic amorphous solids where mechanical and magnetic degrees of freedom interact, leading to rather complex plastic events whose nature must be disentangled. In this paper we uncover the anatomy of the various contributions to some typical plastic events. These plastic events are seen as Barkhausen noise or other "serrated noises." Using theoretical considerations we explain the observed statistics of the various contributions to the considered plastic events. The richness of contributions and their different characteristics imply that in general the statistics of these serrated noises cannot be universal, but rather highly dependent on the state of the system and on its microscopic interactions.

Quasistatic strain-controlled measurements of stress versus strain curves in macroscopic amorphous solids result in a nonlinear-looking curve that ends up either in mechanical collapse or in a steady state with fluctuations around a mean stress that remains constant with increasing strain. It is therefore very tempting to fit a nonlinear expansion of the stress in powers of the strain. We argue here that at low temperatures the meaning of such an expansion needs to be reconsidered. We point out the enormous difference between quenched and annealed averages of the stress versus strain curves and propose that a useful description of the mechanical response is given by a stress (or strain) -dependent shear modulus for which a theoretical evaluation exists. The elastic response is piecewise linear rather than nonlinear.

Amorphous solids yield at a critical value of the strain (in strain-controlled experiments); for larger strains, the average stress can no longer increase-the system displays an elastoplastic steady state. A long-standing riddle in the materials community is what the difference is between the microscopic states of the material before and after yield. Explanations in the literature are material specific, but the universality of the phenomenon begs a universal answer. We argue here that there is no fundamental difference in the states of matter before and after yield, but the yield is a bona fide first-order phase transition between a highly restricted set of possible configurations residing in a small region of phase space to a vastly rich set of configurations which include many marginally stable ones. To show this, we employ an order parameter of universal applicability, independent of the microscopic interactions, that is successful in quantifying the transition in an unambiguous manner.

The determination of the normal and transverse (frictional) interparticle forces within a granular medium is a long-standing, daunting, and yet unresolved problem. We present a new formalism that employs the knowledge of the external forces and the orientations of contacts between particles (of any given size), to compute all the interparticle forces. Having solved this problem, we exemplify the efficacy of the formalism showing that the force chains in such systems are determined by an expansion in the eigenfunctions of a newly defined operator.

We focus on the probability distribution function (PDF) P(Δγ;γ) where Δγ are the measured strain intervals between plastic events in a athermal strained amorphous solids, and γ measures the accumulated strain. The tail of this distribution as Δγ→0 (in the thermodynamic limit) scales like Δγ(η). The exponent η is related via scaling relations to the tail of the PDF of the eigenvalues of the plastic modes of the Hessian matrix P(λ) which scales like λ(θ), η=(θ-1)/2. The numerical values of η or θ can be determined easily in the unstrained material and in the yielded state of plastic flow. Special care is called for in the determination of these exponents between these states as γ increases. Determining the γ dependence of the PDF P(Δγ;γ) can shed important light on plasticity and yield. We conclude that the PDF's of both Δγ and λ are not continuous functions of γ. In slowly quenched amorphous solids they undergo two discontinuous transitions, first at γ=0(+) and then at the yield point γ=γ(Y) to plastic flow. In quickly quenched amorphous solids the second transition is smeared out due to the nonexisting stress peak before yield. The nature of these transitions and scaling relations with the system size dependence of 〈Δγ〉 are discussed.

Much of the progress achieved in understanding plasticity and failure in amorphous solids had been achieved using experiments and simulations in which the materials were loaded using strain control. There is paucity of results under stress control. Here we present a method that was carefully geared to allow loading under stress control either at T=0 or at any other temperature, using Monte Carlo techniques. The method is applied to a model-perfect crystalline solid, to a crystalline solid contaminated with topological defects, and to a generic glass. The highest yield stress belongs to the crystal, the lowest to the crystal with a few defects, with the glass in between. Although the glass is more disordered than the crystal with a few defects, it yields stress much higher than that of the latter. We explain this fact by considering the actual microscopic interactions that are typical of glass-forming materials, pointing out the reasons for the higher cohesive nature of the glass. The main conclusion of this paper is that the instabilities encountered in stress-control condition are the identical saddle-node bifurcation seen in strain control. Accordingly one can use the latter condition to infer the former. Finally we discuss temperature effects and comment on the time needed to see a stress-controlled material failure.

In this paper we focus on the mechanical properties of oligomeric glasses (waxes), employing a microscopic model that provides, via numerical simulations, information about the shear modulus of such materials, the failure mechanism via plastic instabilities, and the geometric responses of the oligomers themselves to a mechanical load. We present a microscopic theory that explains the numerically observed phenomena, including an exact theory of the shear modulus and of the plastic instabilities, both local and system spanning. In addition we present a model to explain the geometric changes in the oligomeric chains under increasing strains.