Source author record

Edan Lerner

Edan Lerner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2023arXiv

Variability of mesoscopic mechanical disorder in disordered solids

Quantifying mechanical disorder in solids, either disordered crystals or glassy solids, and understanding its range of variability are of prime importance, e.g.~for discovering structure-properties relations. The bounds on the degree of mechanical fluctuations in disordered solids and how those depend on solids formation history remain unknown. Here, we study a broadly applicable quantifier of mesoscopic mechanical disorder $χ$, defined via the dimensionless fluctuations of the shear modulus, over a wide variety of disordered computer solids and upon varying different control parameters. $χ$ is intimately related to basic properties of disordered solids, such as elastic constants and plastic deformability, and can be experimentally extracted by wave-attenuation measurements. We find that a large variety of self-organized glassy solids, where disorder is an emergent property, appear to satisfy a generic lower bound on $χ$. On the other hand, we show that $χ$ is unbounded from above, and may diverge in systems driven towards the critical unjamming point. These results highlight basic properties of disordered solids and set the ground for systematically quantifying mechanical disorder across different systems.

preprint2022arXiv

Disordered crystals reveal soft quasilocalized glassy excitations

Structural glasses formed by quenching a melt are known to host a population of low-energy quasilocalized (nonphononic) excitations whose frequencies $ω$ follow a universal $\sim\!ω^4$ distribution as $ω\!\to\!0$, independently of the glass formation history, the interparticle interaction potential or spatial dimension. Here, we show that the universal quartic law of nonphononic excitations also holds in disordered crystals featuring finite long-range order, which is absent in their glassy counterparts. We thus establish that the degree of universality of the quartic law extends beyond structural glasses quenched from a melt. We further find that disordered crystals, whose level of disorder can be continuously controlled, host many more quasilocalized excitations than expected based on their degree of mechanical disorder -- quantified by the relative fluctuations of the shear modulus -- as compared to structural glasses featuring a similar degree of mechanical disorder. Finally, we show that the stability bound on nonlinear quasilocalized excitations -- previously established in structural glasses -- also holds in disordered crystals. Our results are related to glass-like anomalies experimentally observed in disordered crystals. More broadly, they constitute an important step towards tracing the essential ingredients necessary for the emergence of universal nonphononic excitations in disordered solids.

preprint2021arXiv

A simple and broadly-applicable definition of shear transformation zones

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

preprint2021arXiv

An energy-landscape-based crossover temperature in glass-forming liquids

The systematic identification of temperature scales in supercooled liquids that are key to understanding those liquids' underlying glass properties, and the latter's formation-history dependence, is a challenging task. Here we study the statistics of particles' squared displacements $δr^2$ between equilibrium liquid configurations at temperature $T$, and their underlying inherent states, using computer simulations of 11 different computer-glass-formers. We show that the relative fluctuations of $δr^2$ are nonmonotonic in $T$, exhibiting a maximum whose location defines the crossover temperature $T_{\small{\mathsf{X}}}$. Therefore, $T_{\small{\mathsf{X}}}$ marks the point of maximal heterogeneity during the process of tumbling down the energy landscape, starting from an equilibrium liquid state at temperature $T$, down to its underlying inherent state. We extract $T_{\small{\mathsf{X}}}$ for the 11 employed computer glasses, ranging from tetrahedral glasses to packings of soft elastic spheres, and demonstrate its usefulness in putting the elastic properties of different glasses on the same footing. Interestingly, we further show that $T_{\small{\mathsf{X}}}$ marks the crossover between two distinct regimes of the mean $\langleδr^2\rangle$: a high temperature regime in which $\langleδr^2\rangle$ scales approximately as $T^{0.5}$, and a deeply-supercooled regime in which $\langleδr^2\rangle$ scales approximately as $T^{1.3}$. Further research directions are discussed.

preprint2021arXiv

Brittle to ductile transitions in glasses: Roles of soft defects and loading geometry

Understanding the fracture toughness of glasses is of prime importance for science and technology. We study it here using extensive atomistic simulations in which the interaction potential, glass transition cooling rate and loading geometry are systematically varied, mimicking a broad range of experimentally accessible properties. Glasses' nonequilibrium mechanical disorder is quantified through $A_{\rm g}$, the dimensionless prefactor of the universal spectrum of nonphononic excitations, which measures the abundance of soft glassy defects that affect plastic deformability. We show that while a brittle-to-ductile transition might be induced by reducing the cooling rate, leading to a reduction in $A_{\rm g}$, iso-$\!A_{\rm g}$ glasses are either brittle or ductile depending on the degree of Poisson contraction under unconstrained uniaxial tension. Eliminating Poisson contraction using constrained tension reveals that iso-$\!A_{\rm g}$ glasses feature similar toughness, and that varying $A_{\rm g}$ under these conditions results in significant toughness variation. Our results highlight the roles played by both soft defects and loading geometry (which affects the activation of defects) in the toughness of glasses.

preprint2021arXiv

Mechanical disorder of sticky-sphere glasses. I. Effect of attractive interactions

Recent literature indicates that attractive interactions between particles of a dense liquid play a secondary role in determining its bulk mechanical properties. Here we show that, in contrast with their apparent unimportance to the bulk mechanics of dense liquids, attractive interactions can have a major effect on macro- and microscopic elastic properties of glassy solids. We study several broadly-applicable dimensionless measures of stability and mechanical disorder in simple computer glasses, in which the relative strength of attractive interactions -- referred to as `glass stickiness' -- can be readily tuned. We show that increasing glass stickiness can result in the decrease of various quantifiers of mechanical disorder, on both macro- and microscopic scales, with a pair of intriguing exceptions to this rule. Interestingly, in some cases strong attractions can lead to a reduction of the number density of soft, quasilocalized modes, by up to an order of magnitude, and to a substantial decrease in their core size, similar to the effects of thermal annealing on elasticity observed in recent works. Contrary to the behavior of canonical glass models, we provide compelling evidence indicating that the stabilization mechanism in our sticky-sphere glasses stems predominantly from the self-organized depletion of interactions featuring large, \emph{negative} stiffnesses. Finally, we establish a fundamental link between macroscopic and microscopic quantifiers of mechanical disorder, which we motivate via scaling arguments. Future research directions are discussed.

preprint2021arXiv

Mechanical disorder of sticky-sphere glasses. II. Thermo-mechanical inannealability

Many structural glasses feature static and dynamic mechanical properties that can depend strongly on glass formation history. The degree of universality of this history-dependence, and what it is possibly affected by, are largely unexplored. Here we show that the variability of elastic properties of simple computer glasses under thermal annealing depends strongly on the strength of attractive interactions between the glasses' constituent particles -- referred to here as glass `stickiness'. We find that in stickier glasses the stiffening of the shear modulus with thermal annealing is strongly suppressed, while the thermal-annealing-induced softening of the bulk modulus is enhanced. Our key finding is that the characteristic frequency and density per frequency of soft quasilocalized modes becomes effectively invariant to annealing in very sticky glasses, the latter are therefore deemed `thermo-mechanically inannealable'. The implications of our findings and future research directions are discussed.

preprint2021arXiv

Note: Simple argument for emergent anisotropic stress correlations in disordered solids

It is now well-established that mechanical equilibrium in athermal disordered solids gives rise to anisotropic spatial correlations of the coarse-grained stress field that decay in space as $1/r^d$, where $r$ is the distance from the origin, and $d$ denotes the spatial dimension. In this note we present a simple, geometry based argument for the scaling form of the emergent spatial correlations of the stress field in disordered solids.

preprint2020arXiv

Finite-size effects in the nonphononic density of states in computer glasses

The universal form of the density of nonphononic, quasilocalized vibrational modes of frequency $ω$ in structural glasses, ${\cal D}(ω)$, was predicted theoretically decades ago, but only recently revealed in numerical simulations. In particular, it has been recently established that, in generic computer glasses, ${\cal D}(ω)$ increases from zero frequency as $ω^4$, independent of spatial dimension and of microscopic details. However, in [E. Lerner, and E. Bouchbinder, Phys. Rev. E 96, 020104(R) (2017)] it was shown that the preparation protocol employed to create glassy samples may affect the form of their resulting ${\cal D}(ω)$: glassy samples rapidly quenched from high temperature liquid states were shown to feature ${\cal D}(ω)\!\sim\!ω^β$ with $β\!<\!4$, presumably limiting the degree of universality of the $ω^4$ law. Here we show that exponents $β\!<\!4$ are only seen in small glassy samples quenched from high-temperatue liquid states --- whose sizes are comparable to or smaller than the size of the disordered core of soft quasilocalized vibrations --- while larger glassy samples made with the same protocol feature the universal $ω^4$ law. Our results demonstrate that observations of $β\!<\!4$ in the nonphononic density of states stem from finite-size effects, and we thus conclude that the $ω^4$ law should be featured by any sufficiently large glass quenched from a melt.

preprint2020arXiv

Nonlinear quasilocalized excitations in glasses. I. True representatives of soft spots

Structural glasses formed by quenching a melt possess a population of soft quasilocalized excitations --- often called `soft spots' --- that are believed to play a key role in various thermodynamic, transport and mechanical phenomena. Under a narrow set of circumstances, quasilocalized excitations assume the form of vibrational (normal) modes, that are readily obtained by a harmonic analysis of the multi-dimensional potential energy. In general, however, direct access to the population of quasilocalized modes via harmonic analysis is hindered by hybridizations with other low-energy excitations, e.g.~phonons. In this series of papers we re-introduce and investigate the statistical-mechanical properties of a class of low-energy quasilocalized modes --- coined here \emph{nonlinear quasilocalized excitations} (NQEs) --- that are defined via an anharmonic (nonlinear) analysis of the potential energy landscape of a glass, and do not hybridize with other low-energy excitations. In this first paper, we review the theoretical framework that embeds a micromechanical definition of NQEs. We demonstrate how harmonic quasilocalized modes hybridize with other soft excitations, whereas NQEs properly represent soft spots without hybridization. We show that NQEs' energies converge to the energies of the softest, non-hybridized harmonic quasilocalized modes, cementing their status as true representatives of soft spots in structural glasses. Finally, we perform a statistical analysis of the mechanical properties of NQEs, which results in a prediction for the distribution of potential energy barriers that surround typical inherent states of structural glasses, as well as a prediction for the distribution of local strain thresholds to plastic instability.

preprint2020arXiv

Robustness of Density of Low Frequency States in Amorphous Solids

Low frequency quasi-localized modes of amorphous glasses appear to exhibit universal density of states, depending on the frequencies as $D(ω) \sim ω^4$. To date various models of glass formers with short range binary interaction, and network glasses with both binary and ternary interactions, were shown to conform with this law. In this paper we examine granular amorphous solids with long-range electrostatic interactions, and find that they exhibit the same law. To rationalize this wide universality class we return to a model proposed by Gurevich, Parshin and Schober (GPS) and analyze its predictions for interaction laws with varying spatial decay, exploring this wider than expected universality class. Numerical and analytic results are provided for both the actual system with long range interaction and for the GPS model.

preprint2020arXiv

Statistical mechanics of local force dipole responses in computer glasses

Soft quasilocalized modes (QLMs) are universally featured by structural glasses quenched from a melt, and are supposedly involved in a number of glassy anomalies such as the low temperature scaling of their thermal conductivity and specific heat, and sound attenuation at intermediate frequencies. In computer glasses, QLMs may assume the form of harmonic vibrational modes under a narrow set of circumstances, however direct access to their full distribution over frequency is hindered by hybridizations of QLMs with other low-frequency modes (e.g.~phonons). Previous studies to overcome this issue have demonstrated that the response of a glass to local force dipoles serves as a good proxy for its QLMs; we therefore study here the statistical-mechanical properties of the responses to local force dipoles in computer glasses, over a large range of glass stabilities and in various spatial dimensions, with the goal of revealing properties of the yet-inaccessible full distribution of QLMs' frequencies. We find that, as opposed to the spatial-dimension-independent universal distribution of QLMs' frequencies $ω$ (and, consequently, also of their stiffness $κ\!=\!ω^2$), the distribution of stiffnesses associated with responses to local force dipoles features a (weak) dependence on spatial dimension. We rationalize this dependence by introducing a lattice model that incorporates both the real-space profiles of QLMs --- associated with dimension-dependent long-range elastic fields --- and the universal statistical properties of their frequencies. Finally, we discuss possible connections between our findings and basic aspects of the glass transition problem, and to finite-size effects in plastic activity of ultrastable glasses.

preprint2020arXiv

Universality of the nonphononic vibrational spectrum across different classes of computer glasses

It has been recently established that the low-frequency spectrum of simple computer glass models is populated by soft, quasilocalized nonphononic vibrational modes whose frequencies $ω$ follow a gapless, universal distribution ${\cal D}(ω)\!\sim\!ω^4$. While this universal nonphononic spectrum has been shown to be robust to varying the glass history and spatial dimension, it has so far only been observed in simple computer glasses featuring radially-symmetric, pairwise interaction potentials. Consequently, the relevance of the universality of nonphononic spectra seen in simple computer glasses to realistic laboratory glasses remains unclear. Here we demonstrate the emergence of the universal $ω^4$ nonphononic spectrum in a broad variety of realistic computer glass models, ranging from tetrahedral network glasses with three-body interactions, through molecular glasses and glassy polymers, to bulk metallic glasses (BMGs). Taken together with previous observations, our results indicate that the low-frequency nonphononic vibrational spectrum of any glassy solid quenched from a melt features the universal $ω^4$ law, independently of the nature of its microscopic interactions.

preprint2016arXiv

Effect of particle collisions in dense suspension flows

We study non-local effects associated with particle collisions in dense suspension flows, in the context of the affine solvent model known to capture various aspects of the jamming transition. We show that an individual collision changes significantly the velocity field on a characteristic volume $Ω_c\sim 1/δz$ that diverges as jamming is approached, where $δz$ is the deficit in coordination number required to jam the system. Such an event also affects the contact forces between particles on that same volume $Ω_c$, but this change is modest in relative terms, of order $f_{coll}\sim \bar{f}^{0.8}$, where $\bar{f}$ is the typical contact force scale. We then show that the requirement that coordination is stationary (such that a collision has a finite probability to open one contact elsewhere in the system) yields the scaling of the viscosity (or equivalently the viscous number) with coordination deficit $δz$. The same scaling result was derived in [E.~DeGiuli, G.~Düring, E.~Lerner, and M.~Wyart, Phys.~Rev.~E {\bf 91}, 062206 (2015)] via different arguments making an additional assumption. The present approach gives a mechanistic justification as to why the correct finite size scaling volume behaves as $1/δz$, and can be used to recover a marginality condition known to characterize the distributions of contact forces and gaps in jammed packings.

preprint2016arXiv

The micromechanics of nonlinear plastic modes

Nonlinear plastic modes (NPMs) are collective displacements that are indicative of imminent plastic instabilities in elastic solids. In this work we formulate the atomistic theory that describes the reversible evolution of NPMs and their associated stiffnesses under external deformations. The deformation-dynamics of NPMs is compared to those of the analogous observables derived from atomistic linear elastic theory, namely destabilizing eigenmodes of the dynamical matrix and their associated eigenvalues. The key result we present and explain is that the dynamics of NPMs and of destabilizing eigenmodes under external deformations follow different scaling laws with respect to the proximity to imminent instabilities. In particular, destabilizing modes vary with a singular rate, whereas NPMs' exhibit no such singularity. As a result, NPMs converge much earlier than destabilizing eigenmodes to their common final form at plastic instabilities. This dynamical difference between NPMs and linear destabilizing eigenmodes underlines the usefulness of NPMs for predicting the locus and geometry of plastic instabilities, compared to their linear-elastic counterparts.

preprint2015arXiv

Nonlinear plastic modes in disordered solids

We propose a framework within which a robust mechanical definition of precursors to plastic instabilities, often termed `soft-spots', naturally emerges. They are shown to be collective displacements (modes) $\hat{z}_0$ that correspond to local minima of the `barrier function' $b(\hat{z})$. The latter is derived from the cubic approximation of the variation $δU_{\hat{z}}(s)$ of the potential energy upon displacing particles a distance $s$ along $\hat{z}$. We show that modes $\hat{z}_0$ corresponding to low-lying minima of $b(\hat{z})$ lead to transitions over energy barriers in the glass, and are therefore associated with highly asymmetric variations $δU_{\hat{z}}(s)$ with $s$. We further demonstrate how a heuristic search for local minima of $b(\hat{z})$ can a-priori detect the locus and geometry of imminent plastic instabilities with remarkable accuracy, at strains as large as $γ_c-γ\sim 10^{-2}$ away from the instability strain $γ_c$, where the non-affine displacements under shear are still largely delocalized. Our findings suggest that the a-priori detection of plastic instabilities can be effectively carried out by the investigation of the landscape of $b(\hat{z})$.

preprint2014arXiv

Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Athermal plastic flows were simulated for the Kob-Andersen binary Lennard-Jones system and its repulsive version in which the sign of the attractive terms is changed to a plus. Properties evaluated from simulations at different densities include the distributions of energy drops, stress drops, and strain intervals between the flow events. By reference to hidden scale invariance we show that simulations at a single density in conjunction with an equilibrium-liquid simulation at the same density allows one to predict the plastic flow-event properties at other densities. We furthermore demonstrate quasiuniversality of the flow-event statistics.

preprint2014arXiv

Effects of coordination and pressure on sound attenuation, boson peak and elasticity in amorphous solids

Connectedness and applied stress strongly affect elasticity in solids. In various amorphous solids, mechanical stability can be lost either by reducing connectedness or by increasing pressure. We present an effective medium theory of elasticity that extends previous approaches by incorporating the effect of compression, of amplitude $e$, allowing one to describe quantitative features of sound propagation, transport, the boson peak, and elastic moduli near the elastic instability occurring at a compression $e_c$. The theory disentangles several frequencies characterizing the vibrational spectrum: the onset frequency $ω_0\sim \sqrt{e_c-e}$ where strongly-scattered modes appear in the vibrational spectrum, the pressure-independent frequency $ω_*$ where the density of states displays a plateau, the boson peak frequency $ω_{BP}$, and the Ioffe-Regel frequency $ω_{IR}$ where scattering length and wavelength become equal. We predict that sound attenuation crosses over from $ω^4$ to $ω^2$ behaviour at $ω_0$. We predict that a frequency-dependent length scale $l_s(ω)$ and speed of sound $ν(ω)$ characterize vibrational modes, and could be extracted from scattering data. One key result is the prediction of a flat diffusivity above $ω_0$, in agreement with previously unexplained observations. We find that the shear modulus does not vanish at the elastic instability, but drops by a factor of 2. We check our predictions in packings of soft particles and study the case of covalent networks and silica. Overall, our approach unifies sound attenuation, transport and length scales entering elasticity in a single framework where disorder is not the main parameter controlling the boson peak, in agreement with observations. This framework leads to a phase diagram where various glasses can be placed, connecting microscopic structure to vibrational properties.

preprint2014arXiv

Length scales and self-organization in dense suspension flows

Dense non-Brownian suspension flows of hard particles display mystifying properties: as the jamming threshold is approached, the viscosity diverges, as well as a length scale that can be identified from velocity correlations. To unravel the microscopic mechanism governing dissipation and its connection to the observed long-range correlations, we develop an analogy between suspension flows and the rigidity transition occurring when floppy networks are pulled -- a transition believed to be associated to the stress-stiffening of certain gels. After deriving the critical properties near the rigidity transition, we show numerically that suspensions flows lie close to it. We find that this proximity causes a decoupling between viscosity and the correlation length of velocities ξ, which scales as the length l_c characterizing the response of the velocity in flow to a local perturbation, previously predicted to follow l_c\sim 1/\sqrt{z_c-z}\sim p^{0.18} where p is the dimensionless particle pressure, z the coordination of the contact network made by the particles and z_c is twice the spatial dimension. We confirm these predictions numerically, predict the existence of a larger length scale l_r\sim 1/\sqrt{p} with mild effects on velocity correlation and the existence of a vanishing strain δγ\sim 1/p that characterizes de-correlation in flow.

preprint2014arXiv

Phonon gap and localization lengths in floppy materials

Gels of semi-flexible polymers, network glasses made of low valence elements, softly compressed ellipses and dense suspensions under flow are examples of floppy materials. These systems present collective motions with almost no restoring force. We study theoretically and numerically the frequency-dependence of the response of these materials, and the length scales that characterize their elasticity. We show that isotropic floppy elastic networks present a phonon gap for frequencies smaller than a frequency $ω^*$ governed by coordination, and that the elastic response is localized on a length scale $l_c\sim 1/\sqrt{ω^*}$ that diverges as the phonon gap vanishes (with a logarithmic correction in the two dimensional case). $l_c$ also characterizes velocity correlations under shear, whereas another length scale $l^*\sim 1/ω^*$ characterizes the effect of pinning boundaries on elasticity. We discuss the implications of our findings for suspensions flows, and the correspondence between floppy materials and amorphous solids near unjamming, where $l_c$ and $l^*$ have also been identified but where their roles are not fully understood.

preprint2014arXiv

Scaling description of the yielding transition in soft amorphous solids at zero temperature

Yield stress materials flow if a sufficiently large shear stress is ap- plied. Although such materials are ubiquitous and relevant for indus- try, there is no accepted microscopic description of how they yield, even in the simplest situations where temperature is negligible and where flow inhomogeneities such as shear bands or fractures are ab- sent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the Herschel-Bulkley expo- nent characterizing the singularity of the flow curve near the yield stress Σc, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x be- yond which they yield. We argue that the critical exponents of the yielding transition can be expressed in terms of three independent exponents θ, df and z, characterizing respectively the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elasto-plastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions both in two and three dimensions.

preprint2013arXiv

Breakdown of continuum elasticity in amorphous solids

We show numerically that the response of simple amorphous solids (elastic networks and particle packings) to a local force dipole is characterized by a lengthscale $\ell_c$ that diverges as unjamming is approached as $\ell_c \sim (z - 2d)^{-1/2}$, where $z \ge 2d$ is the mean coordination, and $d$ is the spatial dimension, at odds with previous numerical claims. We also show how the magnitude of the lengthscale $\ell_c$ is amplified by the presence of internal stresses in the disordered solid. Our data suggests a divergence of $\ell_c\sim (p_c-p)^{-1/4}$ with proximity to a critical internal stress $p_c$ at which soft elastic modes become unstable.

preprint2013arXiv

Low-energy non-linear excitations in sphere packings

We study theoretically and numerically how hard frictionless particles in random packings can rearrange. We demonstrate the existence of two distinct unstable non-linear modes of rearrangement, both associated with the opening and the closing of contacts. Mode one, whose density is characterized by some exponent θ', corresponds to motions of particles extending throughout the entire system. Mode two, whose density is characterized by an exponent θ != θ', corresponds to the local buckling of a few particles. Mode one is shown to yield at a much higher rate than mode two when a stress is applied. We show that the distribution of contact forces follows P(f) f^{min(θ',θ)}, and that imposing that the packing cannot be densified further leads to the bounds γ >= 1/(2+θ') and γ >= (1-θ)/2, where γ characterizes the singularity of the pair distribution function g(r) at contact. These results extend the theoretical analysis of [M. Wyart, Phys. Rev. Lett 109, 125502 (2012)] where the existence of mode two was not considered. We perform numerics that support that these bounds are saturated with γ \approx 0.38, θ \approx 0.17 and θ' \approx 0.44. We measure systematically the stability of all such modes in packings, and confirm their marginal stability. The principle of marginal stability thus allows to make clearcut predictions on the ensemble of configurations visited in these out-of-equilibrium systems, and on the contact forces and pair distribution functions. It also reveals the excitations that need to be included in a description of plasticity or flow near jamming, and suggests a new path to study two-level systems and soft spots in simple amorphous solids of repulsive particles.

preprint2013arXiv

On the density of shear transformation zones in amorphous solids

We study the stability of amorphous solids, focusing on the distribution P(x) of the local stress increase x that would lead to an instability. We argue that this distribution is singular P(x)x^θ, where the exponent θ is non-zero if the elastic interaction between rearranging regions is non-monotonic, and increases with the interaction range. For a class of finite dimensional models we show that stability implies a lower bound on θ, which is found to lie near saturation. For quadrupolar interactions these models yield θ ~ 0.6 for d=2 and θ~ 0.4 in d=3 where d is the spatial dimension, accurately capturing previously unresolved observations in atomistic models, both in quasi-static flow and after a fast quench.

preprint2013arXiv

Simulations of driven overdamped frictionless hard spheres

We introduce an event-driven simulation scheme for overdamped dynamics of frictionless hard spheres subjected to external forces, neglecting hydrodynamic interactions. Our event-driven approach is based on an exact equation of motion which relates the driving force to the resulting velocities through the geometric information characterizing the underlying network of contacts between the hard spheres. Our method allows for a robust extraction of the instantaneous coordination of the particles as well as contact force statistics and dynamics, under any chosen driving force, in addition to shear flow and compression. It can also be used for generating high-precision jammed packings under shear, compression, or both. We present a number of additional applications of our method.

preprint2012arXiv

A unified framework for non-Brownian suspension flows and soft amorphous solids

While the rheology of non-Brownian suspensions in the dilute regime is well-understood, their behavior in the dense limit remains mystifying. As the packing fraction of particles increases, particle motion becomes more collective, leading to a growing length scale and scaling properties in the rheology as the material approaches the jamming transition. There is no accepted microscopic description of this phenomenon. However, in recent years it has been understood that the elasticity of simple amorphous solids is governed by a critical point, the unjamming transition where the pressure vanishes, and where elastic properties display scaling and a diverging length scale. The correspondence between these two transitions is at present unclear. Here we show that for a simple model of dense flow, which we argue captures the essential physics near the jamming threshold, a formal analogy can be made between the rheology of the flow and the elasticity of simple networks. This analogy leads to a new conceptual framework to relate microscopic structure to rheology. It enables us to define and compute numerically normal modes and a density of states. We find striking similarities between the density of states in flow, and that of amorphous solids near unjamming: both display a plateau above some frequency scale ω* ~ |z_c-z|, where z is the coordination of the network of particles in contact, z_c = 2D where D is the spatial dimension. However, a spectacular difference appears: the density of states in flow displays a single mode at another frequency scale ω_{min} << ω* governing the divergence of the viscosity.

preprint2012arXiv

Toward a microscopic description of flow near the jamming threshold

We study the relationship between microscopic structure and viscosity in non-Brownian suspensions. We argue that the formation and opening of contacts between particles in flow effectively leads to a negative selection of the contacts carrying weak forces. We show that an analytically tractable model capturing this negative selection correctly reproduces scaling properties of flows near the jamming transition. In particular, we predict that (i) the viscosity η diverges with the coordination z as η ~ (z_c-z)^{-(3+θ)/(1+θ)}, (ii) the operator that governs flow displays a low-frequency mode that controls the divergence of viscosity, at a frequency ω_min\sim(z_c-z)^{(3+θ)/(2+2θ)}, and (iii) the distribution of forces displays a scale f* that vanishes near jamming as f*/<f>\sim(z_c-z)^{1/(1+θ)} where θ characterizes the distribution of contact forces P(f)\simf^θ at jamming, and where z_c is the Maxwell threshold for rigidity.

preprint2011arXiv

Direct Estimate of the Static Length-Scale Accompanying the Glass Transition

Characterizing the glass state remains elusive since its distinction from a liquid state is not obvious. Glasses are liquids whose viscosity has increased so much that they cannot flow. Accordingly there have been many attempts to define a static length-scale associated with the dramatic slowing down of supercooled liquid with decreasing temperature. Here we present a simple method to extract the desired length-scale which is highly accessible both for experiments and for numerical simulations. The fundamental idea is that low lying vibrational frequencies come in two types, those related to elastic response and those determined by plastic instabilities. The minimal observed frequency is determined by one or the other, crossing at a typical length-scale which is growing with the approach of the glass transition. This length-scale characterizes the correlated disorder in the system, where on longer length-scales the details of the disorder become irrelevant, dominated by the Debye model of elastic modes.

preprint2011arXiv

Do Athermal Amorphous Solids Exist?

We study the elastic theory of amorphous solids made of particles with finite range interactions in the thermodynamic limit. For the elastic theory to exist one requires all the elastic coefficients, linear and nonlinear, to attain a finite thermodynamic limit. We show that for such systems the existence of non-affine mechanical responses results in anomalous fluctuations of all the nonlinear coefficients of the elastic theory. While the shear modulus exists, the first nonlinear coefficient B_2 has anomalous fluctuations and the second nonlinear coefficient B_3 and all the higher order coefficients (which are non-zero by symmetry) diverge in the thermodynamic limit. These results put a question mark on the existence of elasticity (or solidity) of amorphous solids at finite strains, even at zero temperature. We discuss the physical meaning of these results and propose that in these systems elasticity can never be decoupled from plasticity: the nonlinear response must be very substantially plastic.

preprint2010arXiv

Athermal Nonlinear Elastic Constants of Amorphous Solids

We derive expressions for the lowest nonlinear elastic constants of amorphous solids in athermal conditions (up to third order), in terms of the interaction potential between the constituent particles. The effect of these constants cannot be disregarded when amorphous solids undergo instabilities like plastic flow or fracture in the athermal limit; in such situations the elastic response increases enormously, bringing the system much beyond the linear regime. We demonstrate that the existing theory of thermal nonlinear elastic constants converges to our expressions in the limit of zero temperature. We motivate the calculation by discussing two examples in which these nonlinear elastic constants play a crucial role in the context of elasto-plasticity of amorphous solids. The first example is the plasticity-induced memory that is typical to amorphous solids (giving rise to the Bauschinger effect). The second example is how to predict the next plastic event from knowledge of the nonlinear elastic constants. Using the results of this paper we derive a simple differential equation for the lowest eigenvalue of the Hessian matrix in the external strain near mechanical instabilities; this equation predicts how the eigenvalue vanishes at the mechanical instability and the value of the strain where the mechanical instability takes place.

preprint2010arXiv

Density Scaling of Avalanche Statistics in Amorphous Solids

Stress vs. strain fluctuations in athermal amorphous solids are an example of `crackling noise' of the type studied extensively in the context of elastic membranes moving through random potentials. Contrary to the latter, we do not have a stochastic equation whose statistics agree with the measured ones. On the other hand we show in this Letter that the statistics of the former exhibit 'density scaling' when the interparticle potential can be well approximated by a power law. The distributions of sizes of dissipative events for a wide range of densities and system sizes follow a single scaling function. We find that both the system-size scaling of energy drops and the entire strain interval statistics are invariant to changes in density. We use this to determine accurately the exponents in the scaling laws, establishing that the present crackling noise belongs to a different universality class.

preprint2010arXiv

Plasticity-Induced Anisotropy in Amorphous Solids: the Bauschinger Effect

Amorphous solids that underwent a strain in one direction such that they responded in a plastic manner `remember' that direction also when relaxed back to a state with zero mean stress. We address the question `what is the order parameter that is responsible for this memory?' and is therefore the reason for the different subsequent responses of the material to strains in different directions. We identify such an order parameter which is readily measurable, we discuss its trajectory along the stress-strain curve, and propose that it and its probability distribution function must form a necessary component of a theory of elasto-plasticity.

preprint2010arXiv

Predicting plastic flow events in athermal shear-strained amorphous solids

We propose a method to predict the value of the external strain where a generic amorphous solid will fail by a plastic response (i.e. an irreversible deformation), solely on the basis of measurements of the nonlinear elastic moduli. While usually considered fundamentally different, with the elastic properties describing reversible phenomena and plastic failure epitomizing irreversible behavior, we show that the knowledge of some nonlinear elastic moduli is enough to predict where plasticity sets in.

preprint2010arXiv

Statistical Physics of Elasto-Plastic Steady States in Amorphous Solids: Finite Temperatures and Strain Rates

The effect of finite temperature $T$ and finite strain rate $\dotγ$ on the statistical physics of plastic deformations in amorphous solids made of $N$ particles is investigated. We recognize three regimes of temperature where the statistics are qualitatively different. In the first regime the temperature is very low, $T<T_{\rm cross}(N)$, and the strain is quasi-static. In this regime the elasto-plastic steady state exhibits highly correlated plastic events whose statistics are characterized by anomalous exponents. In the second regime $T_{\rm cross}(N)<T<T_{\rm max}(\dotγ)$ the system-size dependence of the stress fluctuations becomes normal, but the variance depends on the strain rate. The physical mechanism of the cross-over is different for increasing temperature and increasing strain rate, since the plastic events are still dominated by the mechanical instabilities (seen as an eigenvalue of the Hessian matrix going to zero), and the effect of temperature is only to facilitate the transition. A third regime occurs above the second cross-over temperature $T_{\rm max}(\dotγ)$ where stress fluctuations become dominated by thermal noise. Throughout the paper we demonstrate that scaling concepts are highly relevant for the problem at hand, and finally we present a scaling theory that is able to collapse the data for all the values of temperatures and strain rates, providing us with a high degree of predictability.

preprint2010arXiv

Statistical Physics of the Yielding Transition in Amorphous Solids

The art of making structural, polymeric and metallic glasses is rapidly developing with many applications. A limitation to their use is their mechanical stability: under increasing external strain all amorphous solids respond elastically to small strains but have a finite yield stress which cannot be exceeded without effecting a plastic response which typically leads to mechanical failure. Understanding this is crucial for assessing the risk of failure of glassy materials under mechanical loads. Here we show that the statistics of the energy barriers ΔE that need to be surmounted changes from a probability distribution function (pdf) that goes smoothly to zero to a pdf which is finite at ΔE=0. This fundamental change implies a dramatic transition in the mechanical stability properties with respect to external strain. We derive exact results for the scaling exponents that characterize the magnitudes of average energy and stress drops in plastic events as a function of system size.

preprint2010arXiv

What Determines the Yield Stress in Amorphous Solids?

A crucially important material parameter for all amorphous solids is the yield stress, which is the value of the stress for which the material yields to plastic flow when it is strained quasi-statically at zero temperature. It is difficult in laboratory experiments to determine what parameters of the inter-particle potential effect the value of the yield stress. Here we use the versatility of numerical simulations to study the dependence of the yield stress on the parameters of the inter-particle potential. We find a very simple dependence on the fundamental scales which characterize the repulsive and attractive parts of the potential respectively, and offer a scaling theory that collapses the data for widely different potentials and in different space dimensions.

preprint2009arXiv

Locality and Non-locality in Elasto-plastic Responses of Amorphous Solids

A number of current theories of plasticity in amorphous solids assume at their basis that plastic deformations are spatially localized. We present in this paper a series of numerical experiments to test the degree of locality of plastic deformation. These experiments increase in terms of the stringency of the removal of elastic contributions to the observed elasto-plastic deformations. It is concluded that for all our simulational protocols the plastic deformations are not localized, and their scaling is sub-extensive. We offer a number of measures of the magnitude of the plastic deformation, all of which display sub-extensive scaling characterized by non-trivial exponents. We provide some evidence that the scaling exponents governing the sub-extensive scaling laws are non-universal, depending on the degree of disorder and on the parameters of the systems. Nevertheless understanding what determines these exponents should shed considerable light on the physics of amorphous solids.

preprint2009arXiv

Predictive Statistical Mechanics for Glass Forming Systems

Using two extremely different models of glass formers in two and three dimensions we demonstrate how to encode the subtle changes in the geometric rearrangement of particles during the scenario of the glass transition. We construct a statistical mechanical description that is able to explain and predict the geometric rearrangement, the temperature dependent thermodynamic functions and the $α$-relaxation time within the measured temperature range and beyond. The theory is based on an up-scaling to proper variables (quasi-species) which is validated using a simple criterion. Once constructed, the theory provides an accurate predictive tool for quantities like the specific heat or the entropy at temperatures that cannot be reached by measurements. In addition, the theory identifies a rapidly increasing typical length scale $ξ$ as the temperature decreases. This growing spatial length scale determines the $α$-relaxation time as $τ_α\sim \exp(μξ/T)$ where $μ$ is a typical chemical potential per unit length.

preprint2009arXiv

Relations between the Material Mechanical Parameters and the Inter-particle Potential in Amorphous Solids

The shear-modulus and yield-stress of amorphous solids are important material parameters, with the former determining the rate of increase of stress under external strain and the latter being the stress value at which the material flows in a plastic manner. It is therefore important to understand how these parameters can be related to the inter-particle potential. Here a scaling theory is presented such that given the inter-particle potential, the dependence of the yield stress and the shear modulus on the density of the solid can be predicted in the athermal limit. It is explained when such prediction is possible at all densities and when it is only applicable at high densities. These results open up exciting possibilities for designing in principle new materials with desirable mechanical properties.

preprint2009arXiv

Scaling Theory for Steady State Plastic Flows in Amorphous Solids

Strongly correlated amorphous solids are a class of glass-formers whose inter-particle potential admits an approximate inverse power-law form in a relevant range of inter-particle distances. We study the steady-state plastic flow of such systems, firstly in the athermal, quasi-static limit, and secondly at finite temperatures and strain rates. In all cases we demonstrate the usefulness of scaling concepts to reduce the data to universal scaling functions where the scaling exponents are determined a-priori from the inter-particle potential. In particular we show that the steady plastic flow at finite temperatures with efficient heat extraction is uniquely characterized by two scaled variables; equivalently, the steady state displays an equation of state that relates one scaled variable to the other two. We discuss the range of applicability of the scaling theory, and the connection to density scaling in supercooled liquid dynamics. We explain that the description of transient states calls for additional state variables whose identity is still far from obvious.

preprint2009arXiv

Size of Plastic Events in Strained Amorphous Solids at Finite Temperatures

We address the system-size dependence of typical plastic flow events when an amorphous solid is put under a fixed external strain rate at a finite temperature. For system sizes that are accessible to numerical simulations at reasonable strain rates and at low temperatures the magnitude of plastic events grows with the system size. We explain however that this must be a finite size effect; for larger systems there exist two cross-over length-scales $ξ_1$ and $ξ_2$, the first determined by the elastic time-scale and the second by the thermal energy-scale. For system of linear size $L$ larger than any of these scales the magnitude of plastic events must be bounded. For systems of size $L\gg ξ$ there must exist $(L/ξ)^d$ uncorrelated plastic events which occur simultaneously. We present a scaling theory that culminates with the dependence of the cross-over scales on temperature and strain rate. Finally we relate these findings to the temperature and size dependence of the stress fluctuations. We comment on the importance of these considerations for theories of elasto-plasticity.