Large-scale exponential correlations of nonaffine elastic response of strongly disordered materials

This paper demonstrates that while nonaffine displacement fields in strongly disordered materials exhibit power-law decay, their spatial derivatives reveal large-scale exponentially decaying correlations governed by a disorder-dependent heterogeneity length scale $\xi$, a phenomenon validated through both theoretical analysis of correlated random matrices and numerical simulations of various glassy and amorphous systems.

The Gist

The Big Picture: Why Do Glassy Materials Act Weird?

Imagine you have a perfect crystal, like a diamond or a salt crystal. The atoms are arranged in a neat, repeating grid, like soldiers standing in perfect rows. If you push on this crystal, every soldier moves in perfect unison. This is called an affine response. It's predictable, orderly, and easy to calculate.

Now, imagine a disordered material, like window glass, a polymer, or a metal glass. The atoms are jumbled up like a pile of marbles or a bowl of spaghetti. There is no neat grid. If you push on this material, the atoms don't move in unison. Some get pushed hard, some get squeezed, and some wiggle sideways. These messy, local movements are called nonaffine displacements.

For decades, scientists have been trying to figure out: How far does this "messiness" spread? If I push one atom, how far away does the effect of that push ripple through the material?

The Old Theory vs. The New Discovery

The Old Idea (The Ripple in a Pond):
Previously, scientists thought the "messiness" spread out like a ripple in a pond. The further you got from the push, the weaker the effect became, but it never truly stopped; it just faded away slowly. Mathematically, this is called a power-law decay. It's like shouting in a canyon; the echo gets quieter, but you can still hear a faint sound very far away.

The New Discovery (The Sponge and the Wall):
This paper, using advanced math (Random Matrix Theory) and computer simulations, found something surprising. While the displacement of the atoms does fade slowly, the strain (the stretching or twisting) behaves differently.

The authors found that the "messiness" actually hits a wall. It spreads out for a certain distance, but then it dies off exponentially.

• Analogy: Imagine trying to push a heavy couch across a room. If the floor is smooth (ordered), the whole room vibrates. But if the floor is covered in thick, sticky foam (disordered), the vibration stops after a few feet. The foam absorbs the energy.
• The "Heterogeneity Length Scale" ($\xi$): This is the size of that "sticky foam." It's a specific distance. Inside this distance, the material is chaotic and weird. Outside this distance, the material acts like a normal, smooth elastic solid.

The Two Types of "Messiness"

The paper makes a crucial distinction between two ways the material can deform, using a water analogy:

1. The Divergence (Squeezing/Expanding):

   • Imagine: Squeezing a sponge or blowing air into a balloon. The material gets denser in some spots and less dense in others.
   • The Finding: This type of messiness has a long reach. It follows the "sticky foam" rule. It spreads out to the distance $\xi$ and then drops off sharply. This is the "large-scale exponential correlation" mentioned in the title.

2. The Rotor (Twisting/Spinning):

   • Imagine: Twisting a towel or swirling water in a bucket. The material rotates locally without changing its density.
   • The Finding: This is where it gets really cool.
     • If you squeeze the material (volumetric deformation), the "twisting" messiness vanishes almost immediately. It doesn't have that long "sticky foam" reach. It's like a short, sharp snap.
     • However, if you shear (slide) the material, the twisting messiness does have a long reach, but it also has a tiny, lingering "tail" that fades very slowly (the power-law part).

How They Proved It (The Experiments)

The authors didn't just do math; they built digital worlds to test their theory. They simulated three very different types of messy materials:

1. A Broken Net (Rigidity Percolation): Imagine a net where some strings are cut. As they cut more strings, the net gets weaker. They found that as the net gets closer to falling apart, the "sticky foam" distance ($\xi$) gets huge. The material becomes sensitive to pushes over a massive area.
2. Plastic (Polystyrene): They simulated a block of plastic. Even though it's a complex polymer, they found the same "sticky foam" distance. It was about 1.4 nanometers (roughly 3 times the size of a single molecule).
3. Glass (Lennard-Jones Glass): They simulated a standard glass. They found the same pattern: a distinct distance where the chaos stops, and a tiny, lingering twist that fades slowly.

Why Does This Matter?

Think of this discovery as finding the "Rule of Thumb" for how disorder works in materials.

• For Engineers: If you are designing a material that needs to be stiff (like a car part) or soft (like a medical implant), you need to know this "sticky foam" distance. If you put a tiny nanoparticle inside a polymer, the material around it becomes stiffer, but only up to this distance $\xi$. Knowing $\xi$ helps engineers design better nanocomposites.
• For Physicists: It solves a long-standing debate. For years, people argued: "Does disorder fade slowly (power-law) or quickly (exponential)?" The answer is: Both. The displacement fades slowly, but the strain (the actual stress) fades quickly, except for a tiny twist that lingers.

Summary in One Sentence

This paper reveals that in messy, glassy materials, the "chaos" caused by a push travels a specific distance (the heterogeneity length) before hitting a wall and stopping, acting like a sponge that absorbs the disorder, while a tiny bit of "twisting" chaos lingers on forever.

Technical Summary

1. Problem Statement

Amorphous solids (glasses, polymers, gels) exhibit nonaffine elastic deformations, where local atomic displacements deviate from the macroscopic uniform strain due to structural disorder. While the existence of these displacements is well-established, the spatial correlation properties of the nonaffine field remain a subject of debate.

• The Conflict: Previous theoretical work (e.g., DiDonna & Lubensky) suggested that nonaffine correlations decay purely via power laws ($r^{-1}$ in 3D). However, some numerical studies have observed exponential decay at certain length scales.
• The Gap: There is no unified theory explaining when and why exponential decay occurs, how it relates to the disorder strength, and how it manifests in different components of the deformation field (divergence vs. rotor). Furthermore, the role of mechanical stability (positive-definiteness of the force constant matrix) in generating these correlations has not been fully integrated into a random matrix framework.

2. Methodology

The authors employ a combination of Random Matrix Theory (RMT) and numerical simulations to resolve the correlation behavior.

A. Theoretical Framework: Correlated Random Matrices

• Mechanical Stability: The authors model the force-constant matrix $\hat{\Phi}$ of a stable amorphous solid as a correlated Wishart ensemble: $\hat{\Phi} = \hat{A} \cdot \hat{A}^T$. Here, $\hat{A}$ is a rectangular matrix of random variables representing bond interactions. This construction ensures $\hat{\Phi}$ is positive-semidefinite, a prerequisite for mechanical stability.
• Disorder Parameter: They introduce a control parameter $\kappa = 1 - N'_{dof}/N_b$, representing the distance from the isostatic point (where the number of bonds equals the number of nontrivial degrees of freedom).
  • $\kappa \to 0$ corresponds to the critical, strongly disordered limit.
  • $\tilde{\kappa} = 1 - \theta_0$ is used as a generalized measure of disorder, where $\theta_0$ is the largest eigenvalue of the correlation operator.
• Dyson-Schwinger Equations: Using diagrammatic techniques, they derive coupled equations for the resolvent (Green's function) $\hat{G}$ and a secondary resolvent $\hat{G}_b$.
• Eigenvalue Decomposition: The four-point resolvent (covariance of displacements) is analyzed via eigenvalue decomposition of the superoperators. The spectrum consists of:
  • An upper branch ($\theta_0$) close to 1, which dominates large-scale behavior.
  • Lower branches ($\theta_n, n>0$) which are further from 1.
• Key Insight: The proximity of $\theta_0$ to 1 generates a heterogeneity length scale $\xi \sim \tilde{\kappa}^{-1/2}$. As disorder increases ($\tilde{\kappa} \to 0$), $\xi$ diverges, potentially exceeding structural correlation lengths.

B. Numerical Validation

The theory is tested against three distinct models:

1. Rigidity Percolation Model: A distorted FCC lattice with randomly "cut" bonds near the percolation threshold.
2. Molecular Dynamics (MD) of Amorphous Polystyrene: Coarse-grained simulations using the MARTINI force field, quenched to zero temperature.
3. Lennard-Jones (LJ) Glass: A binary Kob-Andersen mixture, quenched at a slower rate to achieve a more equilibrated structure.

In all cases, small strains ($\epsilon = 10^{-5}$) are applied, and the nonaffine displacement fields are projected onto a reference lattice to compute correlation functions.

3. Key Contributions and Results

A. Dual Nature of Correlations

The paper demonstrates that the spatial covariance of nonaffine displacements, $K_{\alpha\beta}(r)$, is a superposition of two distinct terms:

1. Ladder Term ($K^{ld}$): Exhibits power-law decay ($r^{2-d}$). This recovers the classical results of DiDonna and Lubensky and corresponds to the Green's function of a continuum elastic medium.
2. Twisted Term ($K^{tw}$): Exhibits exponential decay ($e^{-r/\xi}$) modulated by power laws. This term arises specifically from the strong disorder and the near-critical nature of the system ($\theta_0 \approx 1$).

B. Divergence vs. Rotor Correlations

A major theoretical advance is the derivation of correlation functions for the divergence (density fluctuations) and rotor (local rotations) of the nonaffine field:

• Divergence Correlation ($K_{div}$):

  • Contains a delta-function (white noise) term from the power-law part.
  • Contains a large-scale exponential tail $e^{-r/\xi}$ governed by the heterogeneity length scale $\xi$.
  • This exponential tail exists for all deformation types (volumetric and shear).

• Rotor Correlation ($K_{rot}$):

  • Volumetric Deformation: The exponential tail vanishes. The correlation is dominated by a delta function and a short-range atomic scale. This is a unique prediction: under pure volume change, the long-range nonaffine correlations in the rotational field are suppressed.
  • Shear Deformation: The exponential tail exists with length scale $\xi$.
  • Power-Law Tail: The theory predicts a small, additional power-law tail ($\sim 1/r$) in the rotor correlation for shear deformation, arising from lower eigenvalue branches (odd branches) that were previously neglected.

C. Numerical Verification

• Rigidity Percolation: Confirmed the exponential decay of $K_{div}$ and $K_{rot}$ (shear) with $\xi$ diverging as $(p-p_c)^{-\nu}$ near the percolation threshold. The critical exponent $\nu \approx 0.37$ was found (slightly deviating from the theoretical 0.5 due to fractal heterogeneity).
• Polystyrene: Observed $\xi \approx 1.4$ nm for divergence (shear and volumetric), which is significantly larger than the monomer size. Crucially, the rotor correlation under volumetric strain showed a much faster decay ($\xi_0 \approx 0.4$ nm), confirming the theoretical vanishing of the long-range tail.
• LJ Glass: Confirmed the existence of the large length scale $\xi \approx 2.5$ (in LJ units) and, importantly, observed the predicted power-law tail ($1/r$) in the rotor correlation function at large distances, validating the contribution of lower eigenvalue branches.

4. Significance and Implications

1. Resolution of the Decay Debate: The paper resolves the conflict between power-law and exponential decay theories by showing that both exist simultaneously. The power law dominates at intermediate scales, while the exponential tail (governed by $\xi$) dominates the large-scale behavior of specific field components (divergence and shear-rotor).
2. Mechanical Stability as a Driver: The work establishes that the mechanical stability constraint (positive-definiteness of $\hat{\Phi}$) is the fundamental origin of the long-range exponential correlations. The disorder is not just random noise; it is correlated to ensure stability, leading to the emergence of the scale $\xi$.
3. New Probe for Criticality: The heterogeneity length scale $\xi$ serves as a new, direct probe for criticality in disordered systems (e.g., rigidity percolation and jamming), offering a way to measure the distance from the isostatic point.
4. Application to Nanocomposites: The findings explain the "stiffening" effect observed around nanoparticles in amorphous matrices. The suppression of nonaffine displacements creates an effective elastic shell with a thickness determined by $\xi$.
5. Methodological Advance: The application of correlated Wishart ensembles to the static elastic response provides a robust analytical tool for studying amorphous solids, extending beyond vibrational properties (like the Boson peak) to static deformation correlations.

In summary, this paper provides a comprehensive theoretical and numerical framework showing that strongly disordered materials possess a large-scale heterogeneity length scale $\xi$ that governs the exponential decay of nonaffine correlations, a feature that is distinct from the structural correlation length and highly sensitive to the type of deformation and the mechanical stability of the system.