Critical fluctuations of elastic moduli in jammed solids

This study reveals that while the average shear modulus in jammed solids depends on the interparticle potential, its sample-to-sample fluctuations exhibit a universal critical exponent independent of both the potential and spatial dimension, providing a key foundation for a unified theoretical description of jamming and its link to sound scattering.

The Gist

Imagine you have a giant jar filled with thousands of tiny marbles. If you shake it gently, they flow like a liquid. But if you squeeze the jar tight enough, the marbles lock together, stop flowing, and suddenly the whole jar becomes hard as a rock. This sudden change from "squishy" to "solid" is called the Jamming Transition. It happens in everything from sandcastles and toothpaste to the cells in your body.

Scientists have long known that as you squeeze these marbles closer to the point where they jam, the material gets stiffer. But this new paper by Shiraishi and Mizuno asks a deeper question: What happens when you look at the differences between different jars of marbles?

Here is the story of their discovery, explained simply.

1. The "Average" vs. The "Wild Cards"

Imagine you have 1,000 different jars, each packed with marbles in a slightly different random pattern.

• The Average: If you measure the stiffness of all 1,000 jars and take the average, you get a predictable result. The paper confirms that this average stiffness depends heavily on the "personality" of the marbles. Are they soft rubber balls (Harmonic spheres) or hard steel balls (Hertzian spheres)? The average stiffness changes differently for each type.
• The Fluctuations (The Wild Cards): Now, look at how much the jars differ from that average. Some jars are surprisingly stiff; others are surprisingly squishy, even though they are packed at the same pressure. The authors call these differences "fluctuations."

The Big Surprise: While the average stiffness depends on the type of marble, the size of the fluctuations does not. Whether you use soft rubber or hard steel, the "wild card" behavior of the jars follows the exact same mathematical rule as they get closer to jamming. It's as if the chaos of the packing process has its own universal language that ignores the specific material.

2. The "Glassy" vs. The "Critical" Zone

The researchers found that the behavior changes depending on how tightly packed the marbles are.

• The Glassy Zone (Tight Packing): When the jar is very full, the marbles are stuck in a messy, frozen state (like a glass of water that's been supercooled). Here, the fluctuations behave one way, and the rules get complicated.
• The Critical Zone (Just About to Jam): As you loosen the pressure slightly, bringing the system to the very edge of jamming, a "critical" state emerges. In this zone, the fluctuations explode in size. The paper shows that this explosion follows a perfect, simple rule that works for both 3D jars (like our marbles) and 2D jars (like coins on a table).

3. The "Ripple Effect" Analogy

Why does this matter? The paper connects these fluctuations to how sound travels through the material.

Think of the jammed marbles as a crowd of people holding hands.

• If everyone holds hands with the exact same tension, a sound wave (a ripple) travels smoothly.
• But in a jammed solid, the "tension" in the connections varies wildly from spot to spot. Some spots are tight; others are loose. This is elastic heterogeneity.

The paper suggests that the "wild card" fluctuations they measured are the source of this uneven tension. When sound waves try to travel through this uneven crowd, they scatter and get absorbed, much like light scattering in fog. This is called Rayleigh scattering.

4. The "Universal Translator"

The most exciting part of the paper is that it helps bridge two different theories that scientists use to describe these materials:

1. Heterogeneous Elasticity Theory (HET): A theory that says "The material is messy, and that messiness causes sound to scatter."
2. Effective Medium Theory (EMT): A theory that tries to simplify the mess into an "average" material.

The authors found that the "messiness" (the fluctuations) they measured acts as a universal translator. It proves that the way sound gets scattered in these materials is directly linked to how much the stiffness varies from sample to sample. Crucially, this link works the same way whether you are looking at a 2D layer of coins or a 3D jar of marbles, and whether the marbles are soft or hard.

The Takeaway

In the world of jammed materials (like sand, foam, or biological tissue), the average behavior tells you about the specific ingredients (the type of particle). But the fluctuations (the differences between samples) tell you about the fundamental laws of nature that govern how disorder turns into rigidity.

This paper shows that these fluctuations are not random noise; they are a critical signal. They reveal a hidden, universal order that dictates how sound travels and how these materials break, regardless of what they are made of or how big they are. It's like discovering that while every crowd of people is unique, the way a rumor spreads through a crowd follows the exact same mathematical pattern everywhere in the world.

Technical Summary

Here is a detailed technical summary of the paper "Critical fluctuations of elastic moduli in jammed solids" by Shiraishi and Mizuno.

1. Problem Statement

The jamming transition is a critical phenomenon where disordered particle systems (foams, granular materials, colloids) acquire rigidity upon compression. While the scaling behavior of average physical quantities (like the shear modulus) near the jamming point is well-characterized, the behavior of sample-to-sample fluctuations remains less understood.

A central challenge in the field is bridging the gap between two major theoretical frameworks for amorphous solids:

• Heterogeneous-Elasticity Theory (HET): Describes anomalous properties (e.g., the Boson peak, Rayleigh scattering of sound) based on the spatial fluctuations of the local shear modulus, quantified by a disorder parameter $\gamma$.
• Effective Medium Theory (EMT): Relies on microscopic parameters like coordination number and pre-stress.

Currently, there is a lack of precise numerical characterization of elastic modulus fluctuations to validate the connection between these theories, particularly regarding how these fluctuations scale with pressure and system dimension, and whether they depend on the specific interparticle potential.

2. Methodology

The authors employed large-scale numerical simulations to investigate ensembles of randomly jammed packings.

• System Models:
  • Particle Types: Binary mixtures of spheres (diameter ratio 1:1.4) with equal mass.
  • Interaction Potentials: Two prototypical models were compared:
    • Harmonic spheres: Potential exponent $\alpha = 2$.
    • Hertzian spheres: Potential exponent $\alpha = 5/2$.
  • Dimensions: Simulations were conducted in both 3D and 2D.
• Simulation Protocol:
  • Packings were generated using the FIRE algorithm to minimize energy at high density, followed by global compression/decompression to target specific pressures ($p$).
  • "Rattler" particles (those with insufficient contacts) were removed to ensure mechanical stability.
  • Samples with negative shear moduli (unstable to shear) were discarded.
• Shear Modulus Calculation:
  • The shear modulus ($G$) was calculated via linear response theory, decomposed into an affine component ($G_A$) and a non-affine component ($G_N$).
  • Two types of moduli were analyzed:
    • Stressed modulus ($G_1$): Includes the effect of contact forces (pre-stress).
    • Unstressed modulus ($G_0$): Calculated by removing the first-derivative force terms in the Hessian matrix, effectively simulating a system with relaxed springs.
• Fluctuation Quantification:
  • To handle non-Gaussian distributions and outliers (common in stressed systems), the authors used a robust, median-based fluctuation quantifier:
    $$\chi_{X,med} = \sqrt{N} \frac{\text{median}[(X - \langle X \rangle)^2]}{\langle X \rangle}$$
  • This was applied to both the excess contact number ($\delta z = z - z_{iso}$) and the shear modulus ($G$).

3. Key Contributions

1. Decoupling of Average vs. Fluctuation Scaling: The study demonstrates that while the average shear modulus depends on the interaction potential ($\alpha$), the fluctuations of the shear modulus follow a universal scaling law independent of the potential.
2. Dimensional Independence of Elastic Fluctuations: The critical exponent governing shear modulus fluctuations is found to be identical in both 2D and 3D systems, contrasting with contact number fluctuations which are dimension-dependent.
3. Validation of HET Parameters: The results provide a direct numerical link between sample-to-sample fluctuations ($\chi_G$) and the disorder parameter ($\gamma$) used in Heterogeneous-Elasticity Theory, offering a pathway to predict acoustic attenuation in amorphous solids.

4. Key Results

A. Average Shear Modulus ($\langle G \rangle$)

• The average modulus exhibits potential-dependent scaling:
  • Harmonic spheres ($\alpha=2$): $\langle G \rangle \sim \delta z^1$
  • Hertzian spheres ($\alpha=5/2$): $\langle G \rangle \sim \delta z^2$
• The unstressed modulus ($G_0$) scales similarly to the stressed modulus ($G_1$) but is approximately half the magnitude near the jamming point, consistent with EMT predictions.

B. Fluctuations of Excess Contact Number ($\chi_{\delta z}$)

• Fluctuations diverge as the system approaches jamming ($\delta z \to 0$).
• 3D: $\chi_{\delta z} \sim \delta z^{-2/5}$
• 2D: $\chi_{\delta z} \sim \delta z^{-3/5}$
• Conclusion: Contact number fluctuations are dimension-dependent.

C. Fluctuations of Shear Modulus ($\chi_G$)

• Universal Scaling: Both stressed ($G_1$) and unstressed ($G_0$) moduli exhibit the same critical exponent regardless of the interaction potential ($\alpha$) or spatial dimension ($d$).
  • Scaling Law: $\chi_G \sim \delta z^{-1/2}$
• Distribution Characteristics:
  • The distribution of the unstressed modulus ($G_0$) remains symmetric and Gaussian-like.
  • The distribution of the stressed modulus ($G_1$) becomes skewed and non-Gaussian at low pressures, featuring significant outliers (samples with anomalously low stiffness) due to localized low-frequency vibrational modes.
• Dimensional Independence: The exponent $-1/2$ holds for both 2D and 3D packings.

5. Significance and Theoretical Implications

• Connection to Correlation Lengths: The divergence of shear modulus fluctuations ($\chi_G \sim \delta z^{-1/2}$) coincides with the divergence of the characteristic length scale ($l_c \sim \delta z^{-1/2}$) associated with low-frequency vibrational modes and non-affine displacements. This suggests that elastic criticality is controlled by a mesoscopic length scale that is robust against microscopic details (potential type) and spatial dimension.
• Implications for Heterogeneous-Elasticity Theory (HET):
  • HET predicts Rayleigh scattering of sound waves with an attenuation coefficient $\Gamma \propto \gamma \Omega^{d+1}$, where $\gamma$ is the disorder parameter derived from elastic fluctuations.
  • The authors show that if the disorder parameter is identified with the sample-to-sample fluctuation ($\gamma \sim \chi_G^2$), then $\gamma \sim \delta z^{-1}$.
  • This scaling aligns with numerical simulations of acoustic attenuation in 2D systems when the characteristic frequency is chosen as $\omega_G \sim \delta z^{1/2}$ (related to the shear modulus).
• Unified Description: The findings suggest that while the average mechanical response is sensitive to the specific interaction potential, the fluctuations and the associated critical length scales are universal. This provides a solid foundation for developing a unified theoretical description of the jamming critical phenomenon that bridges microscopic models (EMT) and mesoscopic theories (HET).

In summary, this work establishes that the critical fluctuations of elastic moduli in jammed solids are governed by a universal, dimension-independent exponent, providing a crucial numerical benchmark for theories describing the anomalous vibrational and acoustic properties of amorphous materials.