🔬 condensed matter

Exceptionally Dense and Resilient Polydisperse Disk Packings

This paper demonstrates that a novel generalized simultaneous particle swap algorithm can systematically construct exceptionally dense and resilient polydisperse disk packings with predictable critical packing fractions by linking the energy landscape of overjammed states to the jamming transition.

Original authors: Sangwoo Kim, Sascha Hilgenfeldt

Published 2026-04-14

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Sangwoo Kim, Sascha Hilgenfeldt

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 are trying to pack a suitcase for a trip. You have a bag full of oddly shaped items: some are huge beach balls, some are tiny marbles, and some are medium-sized shoes.

The Problem:
Usually, when you just throw these items in randomly, you end up with a lot of empty space. If you try to squeeze them in tighter, they might shift around, creating a mess, or they might get stuck in a way that leaves gaps. In the world of physics, this is called "jamming." Scientists want to know: What is the absolute tightest way to pack these specific items so they don't move, but also so they don't break?

For a long time, scientists thought there was a "standard" answer for how full a bag could get before it jammed (about 84% full). But they also knew that how you packed it mattered. If you shook the bag gently, you got one result. If you slammed it, you got another.

The Breakthrough:
This paper introduces a clever new trick to find the absolute densest, most unbreakable packing possible for any mix of sizes, without needing a perfect crystal structure (like a honeycomb).

Here is how they did it, using a few simple analogies:

1. The "Over-Packed" Hotel (The Starting Point)

Instead of trying to pack the suitcase from empty to full, the researchers started with the suitcase completely stuffed, even more than it should be. Imagine the particles (the items) are soft balloons. They inflated them so much that they were squishing against each other hard. This is called an "overjammed" state.

In this squished state, the balloons have a lot of "stress energy" (like a spring being squeezed).

2. The "Magic Swap" (The Algorithm)

Here is the genius part. The researchers invented a computer program that acts like a magical shape-shifter.

The Old Way: You try to move items around one by one to find a better spot. This is slow and often gets stuck in a "good enough" spot.
The New Way: The program looks at the whole crowd of balloons at once. It asks: "If I swap the size of this big balloon with that small one, and rearrange the whole group slightly, can we make the whole system more relaxed?"

They do this swap all at the same time for the whole group. It's like if everyone in a crowded elevator suddenly changed their height and position in a coordinated dance to make everyone more comfortable.

3. Finding the "Goldilocks" State

By controlling how much they "swap" and "relax" the system, they can dial the "stress energy" up or down.

High Energy: The balloons are still very squished. When they slowly let the air out (decompress), the packing falls apart easily at a low density (the standard 84%).
Low Energy: The balloons are in a state of perfect, relaxed harmony. When they let the air out, the packing stays tight and rigid all the way down to a density of 90%.

4. Why This Matters (The "Super-Resilient" Suitcase)

The most exciting discovery is about resilience.

Imagine two suitcases packed to 90% capacity:

Suitcase A (The Old Way): If you push on it or shake it, the items inside shift, slide, and rearrange. It's wobbly.
Suitcase B (The New Way): Because the researchers found a "low-energy" arrangement, the items are locked in a perfect, non-moving dance. If you push on this suitcase, the whole thing moves together like a solid rock. The items don't slide past each other; they move in perfect unison.

The "Secret Sauce" (The Angle Order)

How do they know they found the perfect packing? They didn't look at how many neighbors each item had (a common measure). Instead, they looked at the angles.

Think of three friends standing in a circle holding hands.

In a messy crowd, their arms are stretched at weird, uncomfortable angles.
In the "perfect" packing, their arms form the most natural, comfortable triangles possible given their sizes. The researchers found that if the angles between neighbors are "just right," the whole system becomes incredibly strong and dense.

The Big Picture

This paper tells us that if you have a specific set of particles (like sand, glass beads, or even cells in your body), you don't have to settle for a "random" mess. By using this "simultaneous swap" trick, you can find a hidden, ultra-dense, ultra-strong arrangement that nature usually misses.

It's like finding a way to fold a shirt so perfectly that it takes up half the space of a normal fold, and if you drop it, it doesn't unfold. This could help engineers design stronger materials, better batteries, or understand how our own cells pack together without getting crushed.

1. Problem Statement

The jamming transition—the point at which disordered particle systems transition from a fluid-like state to a rigid, jammed state—is a fundamental problem in soft matter physics. A critical parameter in this transition is the critical packing fraction ( $\phi_c$ ), the density at which the system becomes rigid.

The Challenge: While $\phi_c$ was historically thought to be a universal constant for specific particle types (e.g., $\phi_{rcp} \approx 0.84$ for random close packing of spheres), recent research shows it is highly protocol-dependent. Different preparation methods yield different $\phi_c$ values.
The Gap: Existing protocols (like the Lubachevsky-Stillinger algorithm) typically approach jamming from below (unjamming states) or rely on random compression, often failing to access the highest possible density states for a given set of particle sizes. Furthermore, traditional structural measures (like hexagonal order parameters or radial distribution functions) fail to distinguish between jammed states of different densities or predict their mechanical resilience.
Goal: To develop a method to systematically construct polydisperse disk packings with exceptionally high critical packing fractions (approaching the theoretical limit of hexagonal crystals, $\phi_{hex} \approx 0.907$ ) without crystallization, and to understand the structural and energetic origins of this density.

2. Methodology

The authors employ a novel computational approach based on soft-particle models and a generalized simultaneous particle swap algorithm.

System Model:
- Particles: Frictionless, circular disks with a continuous size distribution (polydisperse).
- Interaction: Repulsive harmonic potential $V(r_{ij})$ acting only when particles overlap ( $r_{ij} < \sigma_{ij}$ ).
- Size Distributions: Gamma, log-normal, and uniform distributions with a coefficient of variation $c_A$ between 0.2 and 0.5.
The Generalized Swap Algorithm:
Unlike traditional methods that compress from a fluid state, this method starts with overjammed states ( $\phi > \phi_c$ ) at a high density (specifically $\phi=1$ ) and relaxes them.
1. Initialization: Generate an initial configuration of $N$ particles with random positions and sizes within a periodic box at $\phi=1$ .
2. Relaxation: Use the FIRE algorithm to relax the system to a metastable state (MS).
3. Network Construction: Identify neighboring particles via radical tessellation to form a contact network.
4. Generalized Angle Swap:
  - Treat the contact network as a system with fixed topology but variable node sizes and positions.
  - Assign harmonic potentials to network bonds.
  - Optimize the network to find a new set of particle sizes $\{A'_i\}$ that minimize energy for a target average network coordination number ( $z_n$ ).
  - Crucial Step: Reassign the original particle sizes $\{A_i\}$ to the new particle positions based on the sorted order of the optimized sizes $\{A'_i\}$ . This effects a simultaneous size swap for the entire system.
5. Final Relaxation: Anneal the new configuration to the nearest metastable state.
6. Decompression: Quasi-statically decompress the system (shrinking particles) until the energy drops to zero, identifying the critical jamming point $\phi_c$ .
Control Parameter: The average network coordination number ( $z_n$ ) serves as the control parameter to tune the energy of the metastable state and, consequently, the resulting $\phi_c$ .

3. Key Contributions

Algorithmic Innovation: The introduction of a generalized simultaneous particle swap algorithm that allows for the systematic construction of overjammed metastable states with specific energy levels and coordination numbers.
Energy-Density Correlation: Establishing a strong, predictable negative correlation between the energy of an overjammed metastable state and the critical packing fraction ( $\phi_c$ ) of the resulting jammed state. Lower energy states yield higher $\phi_c$ .
Discovery of Ultra-Dense States: Demonstrating that for a wide range of polydisperse size distributions, it is possible to find jammed states with $\phi_c$ significantly higher than the standard random close packing ( $\phi_{rcp} \approx 0.84$ ), reaching values $\phi_c > 0.89$ (approaching the monodisperse hexagonal limit of 0.907) without long-range crystalline order.
Identification of a Structural Predictor: Showing that traditional disorder measures (hexagonal order parameter $\Psi_6$ , radial distribution $g(r)$ ) fail to predict $\phi_c$ . Instead, the angle order parameter ( $\Theta$ ), which quantifies the deviation of local triplet angles from sterically optimal packing, is an excellent predictor of both metastable energy and $\phi_c$ .

4. Key Results

Tunability of $\phi_c$ : By varying the target coordination number $z_n$ $z_{n}$ during the swap algorithm, the authors can generate a continuous spectrum of critical packing fractions.
- High-energy states ( $z_n \approx 0$ ) yield $\phi_c \approx 0.84$ (standard random packing).
- Low-energy states ( $z_n \approx 6$ ) yield $\phi_c > 0.89$ .
Structural Resilience:
- Compression/Decompression: High- $\phi_c$ states exhibit affine motion under volumetric changes. Particles move almost perfectly in proportion to the strain, maintaining their local structure. In contrast, low- $\phi_c$ states undergo significant non-affine rearrangements.
- Shear Resistance: High- $\phi_c$ states are exceptionally resilient to shear. They can withstand large cyclic shear strains ( $\gamma$ ) with negligible particle rearrangement, whereas low- $\phi_c$ states rearrange immediately. This behavior mimics "ultrastable glasses."
Universality: The relationship between energy, the angle order parameter $\Theta$ , and $\phi_c$ holds across different size distribution shapes (Gamma, Log-normal, Uniform) and system sizes ( $N=144$ to $1024$).
Structural Indicator: The angle order parameter $\Theta$ (measuring how close local triplets are to a sterically optimal, touching configuration) correlates strongly with $\phi_c$ . Low $\Theta$ (highly optimal local packing) corresponds to high $\phi_c$ .

5. Significance

Fundamental Physics: The work bridges the gap between the energy landscape of overjammed soft particles and the jamming transition. It suggests that the jamming transition is not a single point but a continuous range of states determined by the preparation protocol's ability to access deep minima in the energy landscape.
Glassy Dynamics: The high- $\phi_c$ states resemble ideal glass states (minimal entropy states). The findings suggest that simultaneous swap algorithms are a powerful tool for equilibrating glassy systems to high densities and stability, potentially offering insights into the nature of the glass transition.
Material Design: The ability to construct disordered materials with extraordinary mechanical properties (high density, high shear resistance, structural robustness) without crystallization has potential applications in designing amorphous materials with tailored mechanical, electrical, or thermal behaviors.
Predictive Power: The identification of the angle order parameter $\Theta$ as a robust predictor allows researchers to screen for high-density, stable packings without needing to perform full mechanical simulations.

In summary, the paper demonstrates that by manipulating the energy landscape of overjammed soft particles through a novel swap algorithm, one can efficiently construct disordered, polydisperse packings that are both denser and mechanically more robust than previously thought possible, challenging the conventional understanding of random close packing limits.