pith. sign in

arxiv: 2602.10677 · v2 · submitted 2026-02-11 · ❄️ cond-mat.stat-mech

Analytic Nonlinear Theory of Shear Banding in Amorphous Solids

Pith reviewed 2026-05-16 03:39 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords shear bandingamorphous solidsnonlinear elasticitydisplacement fieldplastic dipolesenergy functionalinstabilityelastic moduli
0
0 comments X

The pith

A derived energy functional has a vanishing Hessian eigenvalue that predicts the critical stress for shear band formation in amorphous solids.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops an analytic nonlinear theory for the shear banding instability in athermal quasi-statically sheared amorphous solids. It derives nonlinear equations for the displacement field that incorporate the effects of plastic events, modeled as distributed dipoles causing stress screening and new length scales. Solving the weakly nonlinear amplitude equation yields analytic expressions for displacement fields in shear bands and shows how elastic moduli control band width from ductile to brittle responses. The central result is an energy functional whose Hessian develops an eigenvalue that reaches zero at the instability, giving a prediction for the critical accumulated stress. A sympathetic reader would care because this offers a first-principles analytic route to the onset of failure in amorphous materials without full numerical simulation.

Core claim

We derive nonlinear equations for the displacement field, including the consequences of plastic deformation on the mechanical response of amorphous solids. The plastic events collectively induce distributed dipoles that are responsible for screening effects and the creation of typical length-scales that are absent in classical elasticity theory. The nonlinear theory exposes an instability that results in the creation of shear bands. By solving the weakly nonlinear amplitude equation we present analytic expressions for the displacement fields associated with shear bands, explaining the role of the elastic moduli that determine the width of a shear band from ductile to brittle characteristics.

What carries the argument

The energy functional derived from the nonlinear displacement equations closed by dipole screening, whose Hessian eigenvalue vanishes at the shear-banding instability.

If this is right

  • Analytic displacement fields become available for shear bands and depend on elastic moduli.
  • Shear band width varies continuously between ductile and brittle limits set by material properties.
  • The critical accumulated stress is given by the vanishing of the Hessian eigenvalue.
  • Instability arises specifically from screening and length scales introduced by plastic dipoles.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The framework could be extended to describe band propagation after onset.
  • Material properties might be tuned to select desired ductility by shifting the critical stress.
  • Dipole screening ideas may apply to instabilities in other disordered systems.
  • Direct comparison to molecular dynamics of sheared glasses would test the predicted critical stress.

Load-bearing premise

Collective plastic events are modeled as distributed dipoles that induce screening and new length scales to close the nonlinear equations for the displacement field.

What would settle it

Measurement of the accumulated stress at shear band onset in a sheared amorphous solid, checked against the value where the Hessian eigenvalue of the derived energy functional reaches zero.

Figures

Figures reproduced from arXiv: 2602.10677 by Avanish Kumar, Itamar Procaccia.

Figure 1
Figure 1. Figure 1: FIG. 1. Examples of shear-band profiles in the ductile case. [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The shear associated with the solution plotted in [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Examples of shear-band profiles in the brittle case. [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Examples of shear profiles in the brittle case. The [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
read the original abstract

The aim of this paper is to offer an analytic theory of the shear banding instability in amorphous solids that are subjected to athermal quasi-static shear. To this aim we derive nonlinear equations for the displacement field, including the consequences of plastic deformation on the mechanical response of amorphous solids. The plastic events collectively induce distributed dipoles that are responsible for screening effects and the creation of typical length-scales that are absent in classical elasticity theory. The nonlinear theory exposes an instability that results in the creation of shear bands. By solving the weakly nonlinear amplitude equation we present analytic expressions for the displacement fields that is associated with shear bands, explaining the role of the elastic moduli that determine the width of a shear band from ductile to brittle characteristics. We derive an energy functional whose Hessian possesses an eigenvalue that goes to zero at the shear-banding instability, providing a prediction for the critical value of the accumulated stress that results in an instability.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper derives nonlinear equations for the displacement field in athermal quasi-static shear of amorphous solids, closing them by modeling collective plastic events as distributed dipoles that generate screening and new length scales absent from classical elasticity. An energy functional is constructed whose Hessian possesses an eigenvalue that vanishes at the shear-banding instability, yielding an analytic prediction for the critical accumulated stress; the weakly nonlinear amplitude equation is solved to obtain explicit displacement fields inside shear bands and to relate elastic moduli to band width and the ductile-to-brittle crossover.

Significance. If the central construction holds, the work supplies an analytic route to the onset and morphology of shear bands in amorphous solids, with explicit dependence on elastic moduli that could unify ductile and brittle responses. This would be a notable advance for the field provided the dipole closure can be placed on firmer footing.

major comments (2)
  1. [Abstract] Abstract: the prediction that the Hessian eigenvalue reaches zero at a specific critical stress is presented as following from the closed nonlinear system, yet the dipole representation used to introduce screening lengths is introduced by postulate rather than derived; because this closure directly determines the functional whose eigenvalue supplies the instability threshold, the critical-stress claim rests on an unverified modeling assumption whose validity in the nonlinear regime is not independently checked.
  2. [Abstract] Abstract (and the derivation of the energy functional): it is unclear whether the dipole strengths or the resulting screening length are fixed by independent microscopic considerations or are effectively calibrated to the same instability data used to test the zero-eigenvalue condition; if the latter, the reported critical stress is not a genuine prediction but a consistency check within the same closure.
minor comments (1)
  1. [Abstract] Abstract: the phrase 'weakly nonlinear amplitude equation' is used without stating its explicit form or the order of the expansion employed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful and detailed report. We address each major comment below. The dipole closure is a modeling assumption motivated by the Eshelby picture, but we agree that its status and the independence of parameters require clearer exposition; we will revise accordingly while maintaining that the critical stress follows directly from the closed equations.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the prediction that the Hessian eigenvalue reaches zero at a specific critical stress is presented as following from the closed nonlinear system, yet the dipole representation used to introduce screening lengths is introduced by postulate rather than derived; because this closure directly determines the functional whose eigenvalue supplies the instability threshold, the critical-stress claim rests on an unverified modeling assumption whose validity in the nonlinear regime is not independently checked.

    Authors: We acknowledge that the distributed-dipole representation is introduced as a closure to capture collective screening by plastic events rather than being derived from the microscopic dynamics within this manuscript. This choice is standard in the amorphous-plasticity literature because each localized rearrangement is mechanically equivalent to an Eshelby inclusion (force dipole). The resulting screening length then enters the energy functional, so the vanishing of the Hessian eigenvalue is indeed a property of the closed model. We will revise the abstract and the derivation section to state the modeling assumption explicitly and to note that its validity in the fully nonlinear regime remains to be verified by direct comparison with particle simulations; no such independent check is provided in the present work. revision: partial

  2. Referee: [Abstract] Abstract (and the derivation of the energy functional): it is unclear whether the dipole strengths or the resulting screening length are fixed by independent microscopic considerations or are effectively calibrated to the same instability data used to test the zero-eigenvalue condition; if the latter, the reported critical stress is not a genuine prediction but a consistency check within the same closure.

    Authors: The dipole strength and the emergent screening length are fixed by independent microscopic considerations: they are expressed in terms of the typical plastic-event size and the local yield strain, both taken from prior athermal-quasistatic simulations that characterize individual rearrangements (references to be added). These parameters are not adjusted to reproduce the instability threshold; once inserted, the nonlinear equations yield the critical accumulated stress as the value at which the lowest Hessian eigenvalue crosses zero. We will add an explicit paragraph in the derivation of the energy functional that traces the screening length back to these microscopic inputs, thereby clarifying that the critical stress is a prediction of the closed theory rather than a consistency check. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation proceeds from explicit modeling assumptions to eigenvalue prediction.

full rationale

The paper postulates collective plastic events as distributed dipoles to close the nonlinear displacement equations and introduce screening lengths, then derives the energy functional from that closed system. The zero-eigenvalue condition on the Hessian is obtained directly from the resulting functional and supplies the critical-stress prediction. No step reduces by construction to a fitted parameter renamed as prediction, no self-citation chain bears the load of the central claim, and the dipole closure is introduced as an ansatz rather than smuggled via prior self-citation. The derivation is therefore self-contained relative to its stated premises and does not exhibit any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The central claim rests on treating plastic events as distributed dipoles (an invented modeling device) and on standard continuum assumptions for the displacement field; no free parameters are explicitly named in the abstract, but elastic moduli appear as inputs that control band width.

free parameters (1)
  • elastic moduli
    Enter the analytic solution for band width and control the ductile-to-brittle transition.
axioms (1)
  • domain assumption Plastic events collectively induce distributed dipoles responsible for screening and new length scales
    Invoked to close the nonlinear equations and to generate the instability absent in classical elasticity.
invented entities (1)
  • distributed dipoles from plastic events no independent evidence
    purpose: To produce screening effects and characteristic length scales in the displacement field
    Postulated without independent falsifiable evidence in the abstract; the entity is introduced to explain the instability.

pith-pipeline@v0.9.0 · 5451 in / 1346 out tokens · 116064 ms · 2026-05-16T03:39:40.641575+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Topological Defects in Amorphous Solids

    cond-mat.mtrl-sci 2026-04 unverdicted novelty 2.0

    Topological defects can be identified in glasses and may provide a first-principles framework for their mechanical response and spatiotemporal dynamics.

Reference graph

Works this paper leans on

60 extracted references · 60 canonical work pages · cited by 1 Pith paper

  1. [1]

    Analysis of the LHS Let us label the first four terms on the LHS of Eq. (35) as Z(0) k , Z(1) k , Z(2) k , Z(3) k , and the final long cubic group as Z(4) k : LHSk = Σ αβ dk,αβ   Z(0) k + ˜Aαkγδ dγ,δα   Z(1) k + 1 2 ˜Aαkγδ ( dµ,γα dµ,δ +dµ,γ dµ,δα )    Z(2) k (43) + ˜Aαβγδ ( dγ,δα dk,β +dγ,δ dk,βα )    Z(3) k + 1 2 ˜Aαβγδ ( dµ,γα dµ,δ dk,β +...

  2. [2]

    (55) Using this in the definitions (48), (50), (52) and (54) we find µ (1) =µ, µ (2) =µ (3) = 0, µ (4) =λ + 2µ

    Simplifying for homogeneous isotropic systems Presently we are ready to simplify the result of the previous subsection, taking into account the explicit form of the elastic tensor in homogeneous isotropic systems: ˜Aαβγδ =λ δαβ δγδ +µ (δαγ δβδ +δαδδβγ ). (55) Using this in the definitions (48), (50), (52) and (54) we find µ (1) =µ, µ (2) =µ (3) = 0, µ (4) =...

  3. [3]

    To this aim we choose the background stress in the form Σαβ = ( 0 Σ Σ 0 )

    Simplifying for simple shear Next we will specialize our equations for the case of simple shear. To this aim we choose the background stress in the form Σαβ = ( 0 Σ Σ 0 ) . (57) At present this tensor is off diagonal in an unspecified co- ordinate frame that we must make explicit. To this aim we introduce a pair n, t of orthonormal unit vectors as- sociated...

  4. [4]

    In that case we can estimate α 4 (f ′)4 ≈ B 4 (f 2 − f 2 0 )2

    Inner solution In the core of the shear band we assume that the quar- tic term in f ′ is dominant over the quadratic term, i.e ǫ2(f ′)2/ 2 ≫ ǫ1. In that case we can estimate α 4 (f ′)4 ≈ B 4 (f 2 − f 2 0 )2. (104) 11 Since shear band solution f (ξ) satisfied |f (ξ)| ≤ f0, we can take twice a square root of Eq. (104) to find f ′ ≈ ( B ǫ2 ) 1/ 4 √ f 2 0 − f 2...

  5. [5]

    The outer solution Finding the outer solution is easy, since f ′ → 0 far from the core, and therefore we face again the condition ǫ1 ≫ ǫ2f ′2, which is identical to the limit considered in the ductile limit Subsect. VIII B. Therefore we can write without further ado an outer solution fout(ξ) = f0 tanh (ξ − ξ0 ℓ0 ) . (110)

  6. [6]

    The matched solution Using the results of the last two subsections we can now offer a matched solution, f (ξ) ≈ f0 sin (ξ − ξ0 Lc ) W (ξ) + f0 tanh (ξ − ξ0 ℓ0 ) [1 − W (ξ)] , (111) where W (ξ) is a smooth function that is ≈ 1 in the core and ≈ 0 in the tails. It is interesting to note that in the brittle limit the width Lc of the shear band is determined b...

  7. [7]

    Guided by Eq

    The unprojected energy functional U [d] We note that the governing equilibrium equation is an Euler-Lagrange equation so it should be derived form the variation of anenergy functional U [d], that is, δU[d] = 0 for all admissible variations η satisfying η|∂ Ω = 0, such that after integration by parts we get a PDE of the form (61). Guided by Eq. 61, we writ...

  8. [8]

    Full unprojected Hessian from the Energy Functional U [d] Let us take two arbitrary perturbations η,ζ . We write the Hessian which is defined as the bilinear second vari- ation of U , δ2U [d;η,ζ ] = ∂ 2 ∂s∂t ⏐ ⏐ ⏐ ⏐ s=t=0 U [d +sη +tζ], (B12) together with periodic boundary conditions, or d fixed on ∂Ω so that the admissible variations satisfy ηk|∂ Ω = 0, ζ...

  9. [9]

    local stiffness

    Second variation: Full unprojected Hessian in bilinear form a. Prestress-gradient term From (B15): 16 δU∞ [d;η] = ∫ Σαβ dk,β ηk,α . We vary this with respect to d in direction ζ: dk,β ↦→ dk,β +tζk,β . Thus δ2U∞ [d;η,ζ ] = ∫ Ω Σαβ ζk,β ηk,α d2x. (B20) This is already symmetric in ( η,ζ ) if Σ αβ is symmetric (interchanging dummy indices α ↔ β ). b. Local t...

  10. [10]

    This is the full unprojected Hessian

    Full unprojected Hessian in bilinear form Collect (B20), (B23), and (B28): δ2U [d;η,ζ ] = ∫ Ω Σαβ ηk,α ζk,β d2x +δ2Uel[d;η,ζ ] + ∫ Ω ηkMkζ (d)ζζd2x, (B29) with Mkζ (d) = Lkζ − 3Tkζqrdqdr and δ2Uel given by (B28). This is the full unprojected Hessian. 17

  11. [11]

    (B30) We must integrate by parts in each term to move derivatives off η

    W riting to operator form HU [d] explicitly Now we define HU [d] such that, δ2U [d;η,ζ ] = ∫ Ω ηk(HU [d]ζ)kd2x for all η. (B30) We must integrate by parts in each term to move derivatives off η. Now let us compute the operator term- by-term. a. Prestress operator part Start from (B20): ∫ Σαβ ηk,α ζk,β = − ∫ ηk(Σαβ ζk,β ),α d2x. So (H∞ ζ)k = − (Σαβ ζk,β ),α ...

  12. [12]

    Material

    Elastic operator part Take the elastic bilinear form (B28) and rewrite it as∫ ηk(· · ·)d2x. (i) “Material” elastic part I1 := ∫ ˜Aαβγδ δuαβ [d;η]δuγδ [d;ζ]d2x. (B33) Insert explicit δuαβ [d;η] from (B17): δuαβ [d;η] = 1 2 ( ηα,β +ηβ,α +ηµ,α dµ,β +dµ,α ηµ,β ) . ThereforeI1 is a sum of terms each linear in ηk,α . After collecting each them, we obtain: I1 = ...

  13. [13]

    (B45) Self-adjointness of HU [d] means: ⟨η, HU [d]ζ⟩ = ⟨ζ, HU [d]η⟩ for all admissible η,ζ

    Self-adjointness We use the L2 inner product ⟨a,b ⟩ = ∫ Ω akbkd2x. (B45) Self-adjointness of HU [d] means: ⟨η, HU [d]ζ⟩ = ⟨ζ, HU [d]η⟩ for all admissible η,ζ. (B46) a. Symmetry of the second variation Because U [d] is a scalar functional in d, its mixed sec- ond derivatives commute: δ2U [d;η,ζ ] = ∂ 2 ∂s∂t ⏐ ⏐ ⏐ ⏐ 0 U [d +sη +tζ] = ∂ 2 ∂t∂s ⏐ ⏐ ⏐ ⏐ 0 U [d...

  14. [14]

    Which deformation pattern wants to appear first?

    Linearization and Eigenvalue Problem In the linear stability analysis, our primary aim is to outline: “Which deformation pattern wants to appear first?” Therefore, it is a selection problem. Selection is always governed by linear stability. Nonlinear terms com- pete after instability starts, but which mode appears first is determined by linear terms

  15. [15]

    Also δuγδ [0;ζ] = 1 2 (ζγ,δ +ζδ,γ )

    Linearization at d = 0 At d = 0: dk,β = 0, σαβ (0) = 0, and Mkζ (0) = Lkζ . Also δuγδ [0;ζ] = 1 2 (ζγ,δ +ζδ,γ ). Thus (HU [0]ζ)k = − (Σαβ ζk,β ),α − (δσαk [0;ζ]),α +Lkζζζ, (B51) where δσαk [0;ζ] = ˜Aαkγδ 1 2 (ζγ,δ +ζδ,γ ). (B52)

  16. [16]

    soft mode

    The instability operator The equilibrium condition is, Ek[d] = − δU δd. (B53) Therefore its linearization is δEk[0] = −H U [0]. (B54) This is the operator used to study linear instability of the equilibrium equation (the “soft mode” of the PDE). So define the instability (Jacobian) operator J[0] := δEk[0] = −H U [0]. (B55) Applying minus to (B54): (J[0]ζ)k...

  17. [17]

    Plane waves are used because they diagonalize the linear operator as shown below

    Fourier modes and the eigenvalue problem We note that in the linearized stability analysis, the equilibrium operator (the equilibrium PDE) and the Hes- sian have constant coefficients, so its eigenfunctions are plane waves (eiq·x) labeled by a continuous wavevector q. Plane waves are used because they diagonalize the linear operator as shown below. They ide...

  18. [18]

    Shear polarization

    Computing the full 2 × 2 eigenvalue matrix We consider the isotropic systems such that the “local” tensor is isotropic: Lkζ = Aδkζ. (B71) Now combine (B66), (B70), (B71) into (B63): M(κ,m ) = κ 2 [ Σ(m)I+µI +(λ +µ )m⊗m ] −A I. (B72) Here,m ⊗ m is the dyadic product. It is a second-rank tensor with components, ( m ⊗ m)ij = mimj. It acts as a projection ope...

  19. [19]

    applied shear di- rection

    Applied Macroscopic Shear and the Shear Polarization Selected by Hessian To avoid any confusion between the “applied shear di- rection” and the “shear polarization selected by the Hes- sian” here we outline the basic difference between these two quantities. The applied shear is a macroscopic quantity. When we apply a pure shear stress Σ αβ = Σ 12 = Σ. This...

  20. [20]

    displacement has a component along a Hessian mode

    Reduction of Displacement Field Along the Soft-Mode Note that, as shown earlier, the linearized Hessian Hlin is a self-adjoint operator acting on vector fields. It has an eigenbasis Hlinψ (n) =λnψ (n), ⟨ψ (m), ψ (n)⟩ =δmn, (B91) where ⟨a, b⟩ := ∫ Ω ak(x)bk(x)d2x. Then any admissible displacement can be expanded as: dk(x) = ∑ n anψ (n) k (x), a n = ⟨ψ (n), ...

  21. [21]

    soft-subspace

    Structure of the soft eigenmode ψ (1) k (x) The linear Hessian Hlin has constant coefficients since we still have d = 0. Therefore the eigenmodes can be chosen as plane waves ψ (1) k (x) = ˆψk eiq·x (B105) Sinceq = |q|n,q ·x = |q|(m ·x). We define a new scalar coordinate ξ, such that m ·x = ξ. Therefore the soft-mode is simplified to ψ (1) k (x) = ˆψk ei|q|ξ....

  22. [22]

    (ansatz) Then, we have dζ =eζf (ξ), d pdqdr =epeqerf 3(ξ)

    Projection along a soft-mode We have dk(x) = ekf (ξ), with ξ =m ·x . (ansatz) Then, we have dζ =eζf (ξ), d pdqdr =epeqerf 3(ξ). (products) Substitute this into RHS k and separate powers of u, −Lkζdζ = −Lkζ (eζf (ξ)) = − (Lkζeζ )f (ξ) , (linear) Tkpqrdpdqdr = (Tkpqrepeqer)f 3(ξ) . (cubic) So altogether, RHSk = − (Lkζeζ )f (ξ) + ( Tkpqrepeqer)f 3(ξ). (C4) N...

  23. [23]

    derivative of (( f ′)2) times derivative operator

    Isotropic and homogenous systems In isotropic and homogeneous systems we can use the following forms Γk α =gδk α, Λαβ = Λδαβ , Λαβ = 1 Λδαβ . (C9) Regarding the isotropic form of the tensor Gαβγδ one could think that it involves two parameters like the elastic tensor Eq. (55). In fact this tensor is fully symmetric, since it originates from a fourth order...

  24. [24]

    J. D. Eshelby. The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 241, 376 (1957)

  25. [25]

    D. L. Malandro and D. J. Lacks. Relationships of shear- induced changes in the potential energy landscape to the mechanical properties of ductile glasses. J. Chem. Phys. , 4593 (1999)

  26. [26]

    C. E. Maloney and A. Lemaitre. Amorphous systems in athermal, quasistatic shear. Phys. Rev. E 74, 016118 (2006)

  27. [27]

    M. E. Gurtin. An introduction to continuum mechanics, Vol. 158. Academic press (1982)

  28. [28]

    J. A. Hudson and J. P. Harrison. Engineering rock 27 mechanics: an introduction to the principles. Elsevier (2000)

  29. [29]

    Steif, F

    P. Steif, F. Spaepen, and J. Hutchinson. Strain local- ization in amorphous metals. Acta Metallurgica 30, 447 (1982)

  30. [30]

    Y. Shi, M. B. Katz, H. Li, and M. L. Falk. Evaluation of the disorder temperature and free-volume formalisms via simulations of shear banding in amorphous solids. Phys. Rev. Lett. 98, 185505 (2007)

  31. [31]

    Y. Fu, Y. Jin, and I. Procaccia. Analytic theory of shear localization in amorphous solids confined by couette ge- ometry. arXiv preprint arXiv:2511.14522 (2025)

  32. [32]

    Barbot, M

    A. Barbot, M. Lerbinger, A. Lemaitre, D. Vandem- broucq, and S. Patinet. Rejuvenation and shear banding in model amorphous solids. Phys. Rev. E 101, 033001 (2020)

  33. [33]

    Ozawa, L

    M. Ozawa, L. Berthier, G. Biroli, A. Rosso, and G. Tar- jus. Random critical point separates brittle and ductile yielding transitions in amorphous materials. Proceedings of the National Academy of Sciences 115, 6656 (2018)

  34. [34]

    Lemaitre, C

    A. Lemaitre, C. Mondal, M. Moshe, I. Procaccia, S. Roy, and K. Screiber-Re’em. Anomalous elasticity and plastic screening in amorphous solids. Phys. Rev. E 104, 024904 (2021)

  35. [35]

    Mondal, M

    C. Mondal, M. Moshe, I. Procaccia, S. Roy, J. Shang, and J. Zhang. Experimental and numerical verifica- tion of anomalous screening theory in granular matter. Chaos, Solitons and Fractals , in press, arXiv preprint arXiv:2108.13334 (2022)

  36. [36]

    B. P. Bhowmik, M. Moshe, and I. Procaccia. Direct mea- surement of dipoles in anomalous elasticity of amorphous solids. Phys. Rev. E 105, L043001 (2022)

  37. [37]

    Kumar, M

    A. Kumar, M. Moshe, I. Procaccia, and M. Singh. Anomalous elasticity in classical glass formers. Phys. Rev. E 106, 015001 (2022)

  38. [38]

    Kupferman, M

    R. Kupferman, M. Moshe, and J. P. Solomon. Met- ric description of singular defects in isotropic materials . Archive for Rational Mechanics and Analysis 216, 1009 (2015)

  39. [39]

    Moshe, I

    M. Moshe, I. Levin, H. Aharoni, R. Kupferman, and E. Sharon. Geometry and mechanics of two-dimensional de- fects in amorphous materials. Proceedings of the National Academy of Sciences 112, 10873 (2015)

  40. [40]

    J. W. Rudnicki and J. Rice. Conditions for the localiza- tion of deformation in pressure-sensitive dilatant materi - als. Journal of the Mechanics and Physics of Solids 23, 371 (1975)

  41. [41]

    J. R. Rice and J. Rudnicki. A note on some features of the theory of localization of deformation. International Journal of solids and structures 16, 597 (1980)

  42. [42]

    Charan, M

    H. Charan, M. Moshe, and I. Procaccia. Anomalous elas- ticity and emergent dipole screening in three-dimensional amorphous solids. Phys. Rev. E 107, 055005 (2023)

  43. [43]

    L. D. Landau and E. M. Lifshitz. Theory of Elasticity, Course of Theoretical Physics (1970)

  44. [44]

    Karmakar, E

    S. Karmakar, E. Lerner, and I. Procaccia. Athermal non- linear elastic constants of amorphous solids. Phys. Rev. E 82, 026105 (2010)

  45. [45]

    H. G. E. Hentschel, S. Karmakar, E. Lerner, and I. Pro- caccia. Do athermal amorphous solids exist? Phys. Rev. E 83, 061101 (2011)

  46. [46]

    Dasgupta, H

    R. Dasgupta, H. G. E. Hentschel, and I. Procaccia. Yield strain in shear banding amorphous solids. Phys. Rev. E 87, 022810 (2013)

  47. [47]

    Fr¨ ohlich

    H. Fr¨ ohlich. Theory of dielectrics: dielectric constant and dielectric loss. Clarendon Press (1949)

  48. [48]

    J. N. Nampoothiri, Y. Wang, K. Ramola, J. Zhang, S. Bhattacharjee, and B. Chakraborty. Emergent elasticity in amorphous solids. Phys. Rev. Lett. 125, 118002 (2020)

  49. [49]

    Y. Fu, H. G. E. Hentschel, P. Kaur, A. Kumar, and I. Procaccia. Odd dipole screening in radial inflation. Phys. Rev. E 110, 065003 (2024)

  50. [50]

    P. Kaur, I. Procaccia, and T. Samanta. Selection prin- ciple for the screening parameters in the mechanical re- sponse of amorphous solids. Phys. Rev. E 111, 015506 (2025)

  51. [51]

    Cohen, A

    Y. Cohen, A. Schiller, D. Wang, J. A. Dijksman, and M. Moshe. Odd-dipole screening in disordered matter. Phys. Rev. Mater. 9, 065602 (2025)

  52. [52]

    Y. Jin, I. Procaccia, and T. Samanta. Intermediate phase between jammed and unjammed amorphous solids. Phys. Rev. E 109, 014902 (2024)

  53. [53]

    Procaccia and T

    I. Procaccia and T. Samanta. Dipole-induced transition in 3 dimensions. Proceedings of the National Academy of Sciences 122, e2427273122 (2025)

  54. [54]

    L. A. Segel. The structure of non-linear cellular solutions to the boussinesq equations. Journal of Fluid Mechanics 21, 345 (1965)

  55. [55]

    M. C. Cross and P. C. Hohenberg. Pattern formation outside of equilibrium. Reviews of modern physics 65, 851 (1993)

  56. [56]

    Vandembroucq and S

    D. Vandembroucq and S. Roux. Mechanical noise de- pendent aging and shear banding behavior of a meso- scopic model of amorphous plasticity. Physical Review B—Condensed Matter and Materials Physics 84, 134210 (2011)

  57. [57]

    Benzi, M

    R. Benzi, M. Sbragaglia, M. Bernaschi, S. Succi, and F. Toschi. Cooperativity flows and shear-bandings: a statis- tical field theory approach. Soft Matter 12, 514 (2016)

  58. [58]

    Nicolas, E

    A. Nicolas, E. E. Ferrero, K. Martens, and J.-L. Barrat. Deformation and flow of amorphous solids: Insights from elastoplastic models. Rev. Mod. Phys. 90, 045006 (2018)

  59. [59]

    Rossi, G

    S. Rossi, G. Biroli, M. Ozawa, G. Tarjus, and F. Zam- poni. Finite-disorder critical point in the yielding trans i- tion of elastoplastic models. Phys. Rev. Lett. 129, 228002 (2022)

  60. [60]

    Rossi and G

    S. Rossi and G. Tarjus. 2022, 093301 (2022)