Analytic Nonlinear Theory of Shear Banding in Amorphous Solids

This paper presents an analytic nonlinear theory of shear banding in athermally sheared amorphous solids, deriving equations that account for plastic-induced dipole screening to explain the instability mechanism, predict shear band width, and determine the critical stress threshold for failure.

The Gist

Imagine a block of glass or a pile of sand. In the world of physics, these are called "amorphous solids." Unlike a crystal (like a diamond), where atoms are lined up in perfect rows, the atoms in these materials are jumbled up randomly, like a crowd of people at a concert with no assigned seats.

For a long time, scientists tried to predict how these materials would break or deform using the same rules they use for perfect crystals. But those rules failed. When you push on glass or sand, it doesn't just bend; it suddenly snaps or forms a narrow, sharp line of damage called a shear band. Think of it like a crack forming in a windshield, but instead of a single line, it's a zone where the material slides past itself.

This paper by Avanish Kumar and Itamar Procaccia offers a new, mathematical "recipe" to predict exactly how and why these shear bands form, and what they look like. Here is the breakdown in simple terms:

1. The Problem: The "Hidden" Mess

When you push on a perfect crystal, it stretches smoothly. But when you push on amorphous solids, tiny, chaotic rearrangements happen inside. The authors call these "plastic events."

• The Analogy: Imagine a crowded room. If you push the crowd, people don't just move in a straight line; they bump into each other, shuffle sideways, and create little whirlpools of movement. In the paper, these whirlpools are called "quadrupoles" (four-pointed shapes of movement).
• The Old Theory: Previous theories treated these whirlpools as if they were evenly spread out, like sugar dissolved in tea. This worked for small pushes but failed to explain the sudden, violent formation of shear bands.
• The New Insight: The authors realized that when the material gets stressed, these whirlpools stop being evenly spread. They start clustering, creating "dipoles" (two-pointed forces) that act like screening charges.
  • Metaphor: Think of these dipoles like a crowd of people holding umbrellas. If they are spread out evenly, the rain (stress) hits everyone equally. But if they cluster together, they create a "shield" or a "screen" that blocks the rain in some spots and lets it pour through in others. This screening creates a specific "length scale"—a natural width for the damage zone.

2. The Big Breakthrough: Non-Linear Math

The paper argues that to understand shear bands, you can't use simple, straight-line math (linear equations). You need non-linear math.

• The Analogy: Imagine driving a car. At low speeds, if you turn the wheel a little, the car turns a little (linear). But at high speeds, a tiny turn of the wheel can send the car into a spin (non-linear).
• The authors derived a new set of equations that account for this "high-speed" behavior of the material. They included two main non-linear effects:
  1. How the shape of the material changes as it deforms (the strain-displacement relationship).
  2. How the "whirlpools" of movement interact with each other when they get crowded (the dipole interactions).

3. The Result: Predicting the "Crack"

By solving these complex equations, the authors found a way to predict the profile of the shear band.

• The "Ductile" (Soft) Case: In materials that are a bit more flexible, the shear band is wide and smooth.
  • Metaphor: Like a slow, gentle slope. The material slides gradually over a wide area. The math predicts this shape looks like a tangent hyperbolic (tanh) curve—a smooth S-shape.
• The "Brittle" (Hard) Case: In materials that are very stiff, the shear band is incredibly sharp and narrow.
  • Metaphor: Like a cliff edge. The material stays still on one side and slides instantly on the other. The math shows that in this case, the "core" of the band behaves differently than the edges, creating a very sharp transition.

4. The "Instability" Switch

The paper also explains when this happens.

• The Analogy: Imagine balancing a pencil on its tip. As long as the wind is light, it stands. But at a specific critical wind speed, it becomes unstable and falls.
• The authors calculated the exact "critical stress" (the wind speed) where the material loses its stability. They found that this happens when a specific mathematical value (an eigenvalue of the "Hessian," which is just a fancy way of saying a stability calculator) drops to zero.
• Once this point is reached, the material can no longer hold its shape uniformly, and the shear band "snaps" into existence.

5. Why This Matters (According to the Paper)

Previous theories could say "a shear band will form," but they couldn't tell you how wide it would be or what its shape would look like without just guessing or using computer simulations.

• This paper provides an analytic theory, meaning it gives a direct formula.
• It explains that the width of the shear band is determined by a competition between the material's stiffness and the "screening" effect of those internal whirlpools.
• It distinguishes between brittle (sharp, sudden breaks) and ductile (slow, wide sliding) materials based on the math of these equations.

Summary

In short, the authors built a new mathematical model that treats amorphous solids (like glass or sand) not as simple springs, but as complex crowds of moving particles. By accounting for how these particles "screen" each other's movements and how they behave non-linearly under stress, they derived a formula that predicts exactly when a material will break and what the resulting "crack" (shear band) will look like, from a smooth slide to a sharp snap.

Technical Summary

Technical Summary: Analytic Nonlinear Theory of Shear Banding in Amorphous Solids

Problem Statement
The mechanical response of amorphous solids (e.g., glasses, granular media) under athermal quasi-static shear differs fundamentally from that of crystalline solids. While crystals exhibit purely elastic responses until yielding, amorphous solids undergo plastic deformations at any strain level, typically manifesting as quadrupolar "Eshelby inclusions." A critical instability in these materials is shear banding, where deformation localizes into narrow bands. This phenomenon ranges from sharp, brittle-like jumps in displacement across a plane to smoother, ductile localization.

Existing theories face significant limitations:

1. Linear Theories: Previous analytic approaches relied on linear approximations between strain and displacement and quadratic expansions of Lagrangians. While these successfully described screening effects and renormalized elastic moduli, they could not predict the specific profile of a shear band or the critical stress for instability onset.
2. Elastoplastic Models: Standard mesoscale models capture avalanche statistics and stress redistribution but are often local in space. They lack an intrinsic continuum length scale, meaning the width of a shear band is determined by numerical regularization or block size rather than physical material properties. Consequently, they cannot uniquely determine the internal strain profile.

The authors aim to develop a nonlinear analytic theory that derives the displacement field equations, predicts the critical accumulated stress for instability, provides an analytic solution for the shear band profile, and distinguishes between ductile and brittle characteristics based on material properties.

Methodology
The authors construct a variational theory based on a Lagrangian formulation that incorporates two essential nonlinearities previously neglected:

1. Nonlinear Strain-Displacement Relation: Moving beyond the linear relation $u_{\alpha\beta} = \frac{1}{2}(\partial_\alpha d_\beta + \partial_\beta d_\alpha)$, the theory utilizes the full geometric definition:
   $$u_{\alpha\beta} = \frac{1}{2} [\partial_\alpha d_\beta + \partial_\beta d_\alpha] + \frac{1}{2}\partial_\beta d_\gamma \partial_\alpha d_\gamma$$
2. Higher-Order Dipole Expansion: While plastic events induce quadrupoles, gradients in the quadrupole density create effective dipoles. The authors expand the Lagrangian to include quartic terms in the dipole field ($P^\alpha$), which are necessary to capture the nonlinear saturation of the instability.

The derivation proceeds through the following steps:

• Derivation of Nonlinear Equations: Starting from an energy functional including background stress ($\Sigma_{\alpha\beta}$), elastic energy, and quadrupole/dipole interactions, the authors minimize the energy with respect to the displacement field ($d$) and quadrupole field ($Q$). This yields a nonlinear partial differential equation (Eq. 32) for the displacement field.
• Projection to a Scalar Amplitude Equation: To solve the complex vector equation, the authors employ a strategy based on the soft mode of the system's Hessian. They argue that near the instability, the displacement field can be projected onto a scalar function $f(\xi)$ varying along the band normal $\xi = \mathbf{m} \cdot \mathbf{x}$. This reduces the problem to a one-dimensional "amplitude equation" (Eq. 61).
• Energy Functional and Hessian Analysis: An energy functional is constructed whose Euler-Lagrange equation matches the derived amplitude equation. The Hessian of this functional is analyzed to identify the critical point where the lowest eigenvalue vanishes, signaling the onset of instability.
• Analytic Solutions: The nonlinear amplitude equation is solved in two limits:
  • Ductile Limit: Where slopes are small, neglecting higher-order gradient terms.
  • Brittle Limit: Where the core of the band is dominated by quartic gradient terms, requiring a matched asymptotic expansion (inner and outer solutions).

Key Contributions and Results

1. Derivation of the Nonlinear Amplitude Equation:
   The central result is Eq. (61), a nonlinear differential equation for the shear band profile $f(\xi)$:
   $$(\mu + \Sigma(m)) f''(\xi) + \frac{3}{2}(\lambda + 2\mu) (f'(\xi))^2 f''(\xi) = -A f(\xi) + B f^3(\xi)$$
   Here, $\Sigma(m)$ represents the projection of the background stress, and $A$ and $B$ are parameters derived from the coupling tensors and screening parameters.

2. Prediction of Critical Stress:
   By analyzing the linearized Hessian of the energy functional, the authors derive a criterion for the onset of instability. The critical condition occurs when the lowest eigenvalue of the Hessian vanishes:
   $$(\mu + \Sigma(m)) \left(\frac{2\pi}{L}\right)^2 = A$$
   This provides a prediction for the critical accumulated stress required to trigger shear banding, linking it directly to the screening parameter $A$ (related to $\tilde{\kappa}_e^2$).

3. Analytic Profiles for Ductile and Brittle Bands:

   • Ductile Solution: In the limit of small slopes, the solution is a tanh profile:
     $$f(\xi) = f_0 \tanh\left(\frac{\xi - \xi_0}{\ell}\right)$$
     The width $\ell$ is controlled by the screening parameter and elastic moduli: $\ell \sim \sqrt{(\mu + \Sigma)/A}$.
   • Brittle Solution: In the limit of large slopes, the core of the band is governed by the quartic term, leading to a sinusoidal core matched to a tanh tail. The core width $L_c$ is determined by the coefficients of the cubic dipole term ($G$), distinct from the screening length. This explains the sharper, more localized nature of brittle shear bands.

4. Intrinsic Length Scales:
   The theory demonstrates that shear banding requires an intrinsic length scale generated by the coupling tensors in the Lagrangian. This scale arises from the competition between destabilizing stress terms and regularizing gradient terms (screening). Unlike elastoplastic models where width is arbitrary, here the width is a physical property determined by material parameters ($\mu, \lambda, A, B, G$).

5. Zero Mode of the Nonlinear Hessian:
   The authors prove that the derivative of the shear band solution, $f'(\xi)$, is the zero mode (eigenfunction with zero eigenvalue) of the nonlinear Hessian. This connects the translational invariance of the kink solution directly to the stability analysis.

Significance and Claims
The paper claims to offer the first analytic theory capable of predicting both the critical stress for shear banding and the detailed spatial profile of the band, distinguishing between ductile and brittle behaviors.

• Resolution of Previous Failures: The authors argue that previous theories (e.g., Rudnicki and Rice) failed to provide a solution because they lacked the necessary nonlinearities to saturate the instability. The cubic term ($B f^3$) in the amplitude equation tames the linear instability, while the nonlinear gradient term controls the core shape.
• Physical Selection of Localization: The theory posits that shear localization is a "true screened continuum object." The screening of long-range elastic interactions (Eshelby fields) by plastic dipoles is essential; without this screening, deformation would spread globally. The theory naturally generates a finite-width band without ad-hoc regularization.
• Distinction from Elastoplastic Models: The paper emphasizes that standard elastoplastic models are incomplete for predicting shear band profiles because they lack intrinsic length scales. In contrast, this theory builds the length scale directly into the Lagrangian via coupling parameters, making the band width a predictive, material-dependent quantity.
• Modesty on Parameters: The authors acknowledge that while the functional form of the solution is derived analytically, specific numerical values for parameters like $B$ (which determines the inner length scale in the brittle limit) are not yet measured. They propose that future numerical simulations on classical glass-forming models are required to determine these parameters and fully validate the theory.

In summary, the paper provides a rigorous variational framework that unifies the concepts of plastic screening, nonlinear elasticity, and instability theory to explain the emergence, criticality, and morphology of shear bands in amorphous solids.