🔬 condensed matter

Analytic theory of shear localization in amorphous solids confined by Couette geometry

This paper presents an analytic theory explaining shear localization in two-dimensional amorphous solids under quasi-static Couette rotation, where stress accumulation triggers plastic avalanches that cause adjacent regions to rotate in opposite directions, enabling a-priori predictions of the resulting displacement fields.

Original authors: Yang Fu, Yuliang Jin, Itamar Procaccia

Published 2026-06-26

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Yang Fu, Yuliang Jin, Itamar Procaccia

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a block of glass. While it looks solid, it is technically an amorphous solid, meaning its internal building blocks are not arranged in a neat, repeating grid like in crystals. Instead, they are jumbled together, much like the grains in a pile of sand or the structure of metallic glasses.

Now, picture this material trapped between two giant, invisible rings: a small one in the middle and a large one on the outside.

The scientists in this paper decided to slowly twist the inner ring, like turning the handle of a jar, while holding the outer ring still. They wanted to see how the material inside would react.

Here is the story of what they found, explained simply:

The Problem: The "Surprise" Twist

If you push on a normal, perfectly elastic rubber band, it stretches smoothly and predictably. If you twist the inner ring of our amorphous material, you might expect the whole thing to twist smoothly, like a spiral staircase.

But that's not what happened. Instead, the material acted like a nervous system having a sudden "spasm." As the stress built up, the material would suddenly snap. When it snapped, something weird occurred: the part of the material near the inner ring started spinning one way, while a patch of material right next to it suddenly started spinning the opposite way.

It's as if you were turning a steering wheel, and suddenly, a section of the dashboard started spinning backward on its own. This is called shear localization. The "twist" (or shear) got stuck in a very specific, narrow zone, rather than spreading out evenly.

The Old Theory vs. The New Theory

Before this paper, scientists tried to explain this using standard "elasticity theory" (the math used for springs and rubber).

The Old Theory: Predicted a smooth, gentle curve where the material just twists gradually from the center out.
The Reality: The material did something chaotic and jagged. The old math failed miserably, like trying to predict a thunderstorm using only a calculator for sunny days.

The New Solution: The "Topological Charges"

The authors of this paper built a new mathematical model to explain the chaos. They realized that when the material snaps, it's not just a simple break. It's like the material is made of tiny, invisible magnets (which they call "topological charges" or "quadrupoles").

Think of these charges like tiny, invisible whirlpools inside the material.

The Setup: As you twist the inner ring, you are building up pressure, like winding up a spring.
The Snap: When the pressure gets too high, these invisible whirlpools suddenly activate. They push and pull on their neighbors.
The Result: These whirlpools create a "screening" effect. They cancel out the smooth twisting in some areas and amplify it in others, creating a sharp boundary where the direction of the twist flips.

The "Magic Numbers"

One of the coolest things the authors found is that the material doesn't just snap randomly. It seems to have a preference for specific "magic numbers" (mathematical values they call $\kappa_e$ ).

Imagine you are tuning a radio. You can't just pick any frequency; you have to hit the specific stations where the signal is clear.

If the material "tunes" to one specific number, the snapping happens right next to the inner ring.
If it tunes to a different number, the snapping happens in the middle of the material, far away from the rings.

The paper shows that the material naturally "chooses" these specific settings based on the size of the rings, much like a guitar string only vibrates at certain notes depending on its length.

Why This Matters (According to the Paper)

The authors created a formula that can predict exactly where this "spasm" will happen and what the movement of the particles will look like, right down to the smallest detail.

They didn't just guess: They compared their new math formula against computer simulations (virtual experiments) and found the two matched perfectly.
They explained the "Why": They showed that the "whirlpools" (plastic events) are the reason the material behaves this way, not just simple stretching.
The Limit: They noted that this theory works for slow, careful twisting. If you twist the material too fast or too hard (like breaking glass), the rules change, and that's a problem for a future study.

In a Nutshell

This paper is like finding the rulebook for a chaotic dance. The material (the dancers) usually moves smoothly, but when the music (stress) gets too loud, they suddenly break into a specific, localized spin. The authors figured out the mathematical steps that tell us exactly where that spin will happen and why the dancers choose to move in opposite directions. They replaced a failed "spring" model with a new "whirlpool" model that perfectly predicts the chaos.

Technical Summary: Analytic Theory of Shear Localization in Amorphous Solids

Problem Statement
Shear localization and shear banding are prevalent phenomena in amorphous solids under shear strain, often leading to catastrophic failure in materials like metallic glasses. While it is established that these phenomena arise from plastic events (avalanches) rather than pure elasticity, current theoretical frameworks lack the ability to predict or describe the detailed displacement fields associated with stress drops during these events. Specifically, there is no existing theory that can provide a-priori predictions for the displacement field (particle movement) during plastic avalanches in a quasi-statically loaded amorphous solid confined by Couette geometry. Furthermore, while shear banding (system-spanning events at yielding) is distinct, the mechanisms governing smaller-scale shear localization occurring both before and after yielding require a dedicated theoretical approach.

Methodology
The authors combine numerical simulations with a newly derived analytic theory to address this gap.

Numerical Simulations: The study utilizes a two-dimensional model of an amorphous solid consisting of $N=3534$ bidisperse repulsive disks (radii 1 and 1.4) confined between two concentric rings ( $R_{in}=28$ , $R_{out}=80.8$ ). The system is prepared at an area fraction $\phi=0.91$ and brought to mechanical equilibrium. A quasi-static loading protocol is applied where the inner ring is rotated anti-clockwise in small steps ( $\delta\theta = 0.024^\circ$ ), with the system re-equilibrated via damped molecular dynamics after each step. The outer boundary is fixed. The study focuses on the displacement fields associated with sharp stress drops (plastic events) rather than the elastic response.
Analytic Theory: The authors extend previous theoretical frameworks that treat plastic events as topological quadrupoles and dipoles in the displacement field.
- They formulate a Lagrangian ( $U$ ) that includes the cost of quadrupoles, quadrupole-quadrupole interactions, and the coupling between the strain field and quadrupole fields.
- Crucially, the theory accounts for a substantial background stress ( $\Sigma_{r\theta}$ ) accumulated prior to the plastic event, a factor not fully addressed in previous single-step load studies.
- By minimizing the Lagrangian with respect to the quadrupolar field ( $Q$ ) and the displacement field ( $d$ ), the authors derive a renormalized elastic tensor.
- The inclusion of dipolar contributions (arising from gradients in the quadrupolar field) leads to a governing differential equation for the displacement field: $\Delta d + (1 + \tilde{\lambda})\nabla(\nabla \cdot d) + \Gamma d = 0$ .
- This equation is solved analytically using Fourier series expansions for the radial ( $d_r$ ) and tangential ( $d_\theta$ ) components, subject to the specific boundary conditions of the Couette geometry.

Key Contributions and Results

Analytic Prediction of Displacement Fields: The primary contribution is the derivation of an analytic solution (Eq. 14) for the angle-averaged tangential displacement field $d_\theta(r)$ . This solution, expressed in terms of Bessel functions ( $J_1$ and $Y_1$ ) and a screening parameter $\kappa_e$ , successfully reproduces the complex displacement profiles observed in simulations, including the reversal of particle rotation direction (clockwise vs. anti-clockwise) that classical elasticity theory fails to predict.
Mechanism of Shear Localization: The theory demonstrates that shear localization is determined by discrete values of the screening parameter $\kappa_e$ $κ_{e}$ . These values are dictated solely by the geometric parameters ( $R_{in}$ $R_{in}$ and $R_{out}$ $R_{o u t}$ ) of the system, not the specific state of the material.
- When $\kappa_e \approx 0.067$ , shear localization occurs adjacent to the inner boundary.
- When $\kappa_e \approx 0.123$ , shear localization occurs in the bulk of the material.
- The system preferentially selects $\kappa_e$ values where the denominator of the analytic solution approaches zero (singularities), corresponding to maximal system response.
Validation: The analytic predictions show excellent agreement with numerical simulations. The theory correctly predicts the location and shape of the shear localization for the majority of the 200 largest stress drops analyzed. The paper notes that large stress drops tend to localize near the boundary, while smaller ones localize in the bulk, consistent with the selection of different $\kappa_e$ modes.
Distinction from Shear Banding: The authors clarify that their linear theory (expanded to quadratic order in the Lagrangian) applies to quasi-static shear steps where the background stress is divergence-less. It is explicitly stated that this theory does not describe shear banding associated with yielding, which requires a nonlinear extension (e.g., adding quartic dipole terms) beyond the scope of this work.

Significance
The paper claims to provide the first analytic theory capable of predicting the angle-averaged displacement field associated with plastic drops in sheared amorphous solids. By incorporating topological charges (quadrupoles and dipoles) and accounting for background stress, the theory bridges the gap between pure elasticity (which fails to describe these phenomena) and complex numerical simulations. The ability to predict a priori where shear localization will occur based on geometric constraints offers a fundamental understanding of ductile response in amorphous materials. The authors emphasize that while the theory explains localization in ductile solids, the mechanism for brittle shear banding remains a subject for future nonlinear extensions.