Simulation of shear strain at arbitrary angles as a probe of packing instabilities

This paper introduces a simulation tool for applying shear strain at arbitrary angles to reveal that packing instabilities in disordered solids form continuous lines in phase space, exhibit diverse interaction behaviors, and generate a proliferation of small hysterons as the strain angle approaches the point of instability disappearance.

The Gist

Imagine a crowded dance floor filled with people (the particles) packed so tightly they can barely move. This is what physicists call a "jammed packing" or a disordered solid, like sand, glass, or even a pile of oranges.

When you push this crowd gently, they shuffle a little and then snap back to their original spots when you stop. But if you push hard enough, the crowd gets chaotic. People bump into each other, trip, and rearrange themselves into a new formation. Once they do this, they can't go back to exactly where they started. This is failure or plastic deformation.

For a long time, scientists studied what happens when you push this crowd from just one specific direction (like pushing from the left). But in the real world, forces come from all angles. This paper introduces a new way to simulate pushing the crowd from any angle you want, not just straight on.

Here is a breakdown of their discoveries using simple analogies:

1. The "Magic Push" Tool

Usually, computer simulations of these crowds are stuck in a box with rigid walls. If you want to push the crowd diagonally, the box shape gets weird, and the simulation breaks.

• The Analogy: Imagine trying to stretch a rubber sheet with a pattern on it. If you pull it straight, it's easy. If you pull it diagonally, the pattern gets distorted, and the edges don't line up anymore.
• The Solution: The authors built a "magic tool" that reshapes the invisible box around the particles every time they change the angle of the push. This lets them push the crowd from any angle (0 to 180 degrees) without the simulation crashing, allowing them to see the "bulk" of the material without the walls interfering.

2. The "Instability Lines" (The Fault Lines)

When they pushed the crowd at different angles, they found that the "tripping points" (where the particles rearrange) aren't random. They form lines.

• The Analogy: Think of a map of a forest. If you walk in a straight line, you might hit a specific tree that falls over. If you walk at a slightly different angle, you might hit the same tree.
• The Discovery: They found that a single "weak spot" in the material can be triggered by pushing from a wide range of angles. It's like a specific tree that will fall whether you push it from the North, Northeast, or East. These "lines" of failure can stretch across huge angles (sometimes over 90 degrees!).

3. The "Ghost Dancers" (Crossing Lines)

Sometimes, two different "weak spots" exist in the crowd. When the authors changed the angle of the push, they saw these two weak spots interact.

• The Analogy: Imagine two dancers on a stage. As the music (the angle of the push) changes, they might cross paths. In a normal crowd, they would bump into each other and change their dance. But here, the scientists found that these two "instabilities" could pass right through each other like ghosts. They didn't interfere; they just kept dancing their own separate dances, even though they occupied the same space at the same time.

4. The "Fading Echo" (Closing Hysterons)

Some rearrangements are like a "memory" of the push. If you push the crowd forward, they move. If you push them back, they move back, but not quite to the start. The difference between the forward and backward move is called hysteresis (or a "hysteron").

• The Analogy: Think of a door that sticks. You have to push it hard to open it (forward), and it takes a different amount of force to close it (backward). The gap between the "open" and "close" forces is the hysteresis.
• The Discovery: As they changed the angle of the push toward a "magic angle" where the instability disappears, this "sticking" gap got smaller and smaller until it vanished.
  • The Result: Near this magic angle, the material creates thousands of tiny, almost invisible "stickiness" events. It's like the door is so well-oiled at that specific angle that it barely sticks at all. This explains why materials can store complex memories of how they were pushed.

5. The "Detour" (Circumventing Instability)

Finally, they found a way to trick the material.

• The Analogy: Imagine you are walking through a minefield. If you walk straight, you hit a mine. But if you walk in a big circle around the mine, then come back, you never trigger it.
• The Discovery: By changing the angle of the push, going around the "weak spot," and then returning, they could reach the exact same final arrangement of particles without ever triggering the rearrangement. This proves that the history of how you got to a state matters just as much as the state itself.

Why Does This Matter?

This research is like upgrading from a black-and-white TV to a 3D movie. Previously, scientists could only see how materials break when pushed from one direction. Now, they can see the full 360-degree landscape of failure.

This helps us understand:

• How materials remember: Why a glass or a pile of sand "remembers" how it was shaken in the past.
• Predicting failure: Understanding that weak spots aren't just single points, but long lines that can be triggered from many angles.
• Designing better materials: By knowing how these "ghost dancers" and "fading echoes" work, we might be able to design materials that are tougher or have specific memory properties.

In short, the authors built a new lens to look at how messy materials break, revealing that the story of failure is much more complex, interconnected, and fascinating than we previously thought.

Technical Summary

1. Problem Statement

Disordered solids (such as jammed packings and glasses) deform reversibly under small stresses but undergo catastrophic, irreversible failure when subjected to larger shear strains. This failure occurs via localized rearrangements at "soft spots," "shear-transformation zones," or "hysterons."

While previous research has focused on predicting where these failures occur, less attention has been paid to how the specific details of the forcing (specifically the direction of shear) affect the response. Traditional simulations often use fixed boundary conditions or uniaxial shear (either pure or simple shear along a single axis), which limits the ability to observe how instabilities evolve, interact, or disappear as the loading direction changes. The authors aim to develop a method to probe these instabilities across a continuous range of shear angles to understand the geometry and persistence of failure events in the phase space of strain directions.

2. Methodology

The authors developed a novel simulation framework to apply shear strain at continuously variable angles ($\theta$) in systems with periodic boundary conditions (PBC).

• Formalism for Arbitrary Shear:

  • Standard PBCs require the simulation box to be a parallelogram where particles interacting across boundaries are replicas of those on the opposite side.
  • The authors derived transformation matrices ($\Gamma$) that define the box shape for pure shear ($\Gamma_p$) and simple shear ($\Gamma_s$) at any arbitrary angle $\theta$.
  • These matrices are derived by conjugating the standard shear matrices (at $\theta=0$) with a rotation matrix $R_\theta$. This ensures the box volume remains constant (for pure shear) and that the periodicity is maintained correctly.
  • The formalism is applicable in $d$-dimensions (demonstrated for $d=2$ and $d=3$).

• Simulation Protocol:

  • System: Jammed packings of soft particles (discs in 2D, spheres in 3D) interacting via a repulsive Hertzian potential.
  • Initialization: Random configurations are generated and quenched to mechanical equilibrium using the Fast Inertial Relaxation Engine (FIRE).
  • Shearing: Athermal quasistatic shear is applied in small strain steps ($5\times10^{-5}$ to $5\times10^{-3}$).
  • Instability Detection: Rearrangements are identified by testing for reversibility. The system is sheared forward one step, minimized, sheared back, and minimized again. If the total energy or particle positions differ, an irreversible rearrangement has occurred.
  • Correlation Analysis: To compare rearrangements at different angles, particle displacements are mapped back to a unit square cell. A normalized correlation coefficient ($C$) is calculated between displacement vectors of rearrangements at angles $\theta_1$ and $\theta_2$.

3. Key Contributions

1. Generalized Shear Formalism: The derivation of transformation matrices allowing for the application of pure and simple shear at arbitrary angles within periodic boundary conditions, enabling bulk material studies without surface effects.
2. Mapping Instability Lines: The ability to map "instability lines" in phase space, showing that a single rearrangement event can persist over a broad angular range (sometimes $>90^\circ$).
3. Characterization of Instability Interactions: A systematic classification of how instability lines behave as the shear angle varies:
   • Crossing: Lines can pass through one another with minimal interaction.
   • Merging/Transition: Lines can smoothly transition into one another or abruptly switch.
   • Termination: Lines can end abruptly or fade out.
4. Hysteron Scaling Laws: A detailed investigation of "hysterons" (rearrangements that undo themselves upon reversing shear) near their disappearance point, revealing specific scaling behaviors.

4. Key Results

• Persistence of Instabilities:

  • Instabilities are not isolated to specific angles; they form continuous "lines" in the $(\theta, \gamma)$ phase space.
  • A single rearrangement can persist over angular ranges exceeding $90^\circ$. This implies that the "soft spot" has an inherent $180^\circ$ symmetry, and the triggering condition depends on the projection of the strain matrix onto the feature's orientation, rather than the total strain magnitude.

• Interaction of Rearrangements:

  • Crossing: In small systems where spatial avoidance is impossible, two distinct rearrangement lines can cross. The particle displacements for the two events remain distinct, allowing them to "pass through" each other with minimal interaction.
  • Termination Modes:
    • Abrupt End: A line ends without a nearby alternative.
    • Smooth Transition: A line gradually changes into another (indicated by smooth color changes in correlation matrices).
    • Abrupt Switch: A line jumps to a different rearrangement without a large strain jump, detectable only via correlation analysis.

• Closing of Hysterons:

  • As the shear angle $\theta$ approaches the point where a hysteron disappears, the hysteresis loop width ($\gamma_h = \gamma_+ - \gamma_-$) shrinks to zero.
  • Scaling Behavior: As $\gamma_h \to 0$, the energy drop ($\Delta E$) and particle displacement norm ($\Delta r$) scale smoothly with the hysteresis width:
    • $\Delta E \propto \gamma_h^{1.55 \pm 0.04}$
    • $\Delta r \propto \gamma_h^{0.57 \pm 0.05}$
  • This scaling is intrinsic to individual rearrangements, distinct from ensemble averages.

• Path Dependence and Memory:

  • The authors demonstrated that one can circumvent an instability entirely by taking a path in phase space that goes "around" the rearrangement line (e.g., changing the angle to avoid the critical strain).
  • This confirms that identical final particle configurations can be reached via paths that include rearrangements and paths that do not, highlighting the complex topology of the energy landscape.

5. Significance

This work fundamentally advances the understanding of failure in amorphous solids by moving beyond uniaxial loading constraints.

• Geometric Signature of Failure: It suggests that failure events possess a geometric signature (orientation-dependent) that persists even when the local rearrangement does not look like a classic quadrupole.
• Memory Formation: By revealing how hysterons close and how rearrangements persist or merge, the study provides a deeper mechanism for how disordered solids "remember" cyclic deformations.
• Universality: The technique is applicable to various dimensions, system sizes, and material types (network solids, suspensions, fluids), offering a unified framework to study plasticity and jamming.
• Predictive Power: Understanding the angular dependence of instabilities helps in predicting how materials will respond to complex, multi-axial stress states, which is crucial for engineering applications involving disordered materials.