---
title: Maximally Random Jammed Packings
url: https://www.emergentmind.com/topics/maximally-random-jammed-packings
type: topic
---

# Maximally Random Jammed Packings

Maximally random jammed (MRJ) packings represent the most disordered configurations of nonoverlapping particles that are strictly jammed, i.e., mechanically stable against all collective and global deformations. MRJ packings form a paradigmatic model for disordered, rigid, athermal materials—glassy states poised at the margins of mechanical and dynamic stability. Central features include isostaticity, hyperuniform suppression of large-scale density fluctuations, universal quasi-long-range correlations, and dimension-dependent exponents quantifying their structural singularities. The MRJ paradigm generalizes beyond monodisperse spheres to arbitrary convex shapes, polydisperse mixtures, and high-dimensional Euclidean and confined geometries.

## 1. Variational Definition, Mechanical Stability, and Isostaticity

The MRJ state is defined as the strictly jammed configuration that minimizes a specified scalar order metric ψ, quantifying disorder. Strict jamming implies incompatibility with any collective displacement or volume-nonincreasing boundary deformation; no strictly jammed packing can be continuously unjammed without overlaps. Maxwell counting for frictionless spheres in dimension $d$ yields an isostatic contact number per particle, $z=2d$, such that the total number of mechanical contacts $N_c$ satisfies $2N_c = Nd$ in the thermodynamic limit. For monodisperse spheres in $d=3$, the MRJ state occurs at $z\approx6.00$ and a packing fraction $\phi_{MRJ}\approx0.64$; almost every backbone sphere participates in exactly the minimal number of contacts required for mechanical stability [1101.1327, 1008.2982].

MRJ states are uniquely distinguished among jammed packings by being simultaneously mechanically rigid and maximally disordered with respect to ψ. The variational principle $\min_\Omega \psi(\Omega)$ subject to jamming, and (if desired) fixed $\phi$, singles out the MRJ configuration from the generically nonunique set of jammed packings at each density [1101.1327].

## 2. Universality of Hyperuniformity and Quasi-Long-Range Correlations

A defining feature of MRJ packings is disordered hyperuniformity: the anomalous suppression of infinite-wavelength density (number or local-volume-fraction) fluctuations. Hyperuniformity is verified by the structure factor $S(\mathbf{k})$ vanishing as $|\mathbf{k}|\to0$, typically with a power law $S(k)\sim A k^\alpha$ where $\alpha>0$; the corresponding real-space pair correlations decay as $g_2(r)-1\sim -C r^{-(d+\alpha)}$ for large $r$. For monodisperse spheres, $\alpha=1$ in $d=3$, leading to $g_2(r)-1\sim -C r^{-4}$, and in general, $\alpha=d-2$ for $d\geq3$ [2311.06970]. These quasi-long-range (QLR) correlations are universal across particle shapes (spheres, ellipsoids, superdisks, superballs, truncated tetrahedra), sizes (monodisperse and binary), and even in confined geometries [1008.2548, 1104.4093, 1001.0423, 2111.14954, 1406.1535, 1512.08740].

Hyperuniformity is consistently observed in the two-phase medium (local-volume-fraction fluctuations), even when the point pattern of centers is not hyperuniform due to polydispersity or anisotropy [1104.4091]. The suppression of large-scale fluctuations is tied to the rigidity of the void space—loops of contacting particles enclosing all voids in finite, highly regular shapes [1104.4091].

## 3. Marginal Stability, Quasicontacts, and Weak-Force Distributions

MRJ packings are poised at a double marginal: statically at the isostatic contact number, and dynamically at the onset of unstable collective rearrangements. This marginality is manifest in two signature power-laws:

- **Quasicontacts**: The distribution of near contacts—pairs at gap $\xi\ll1$—exhibits a singularity in the pair correlation function, $g(r)\sim A_d(r-1)^{-\gamma}$ for $r\to1^+$, with $\gamma\approx0.4$ in low dimensions. The count of quasicontacts within gap $\xi$ scales as $Z_{nc}(\xi)\sim A_d \xi^{1-\gamma}$.
- **Weak contact forces**: The distribution of normal contact forces follows $P(f)\sim f^{\theta}/\langle f\rangle^{1+\theta}$ for $f\ll\langle f\rangle$, with $\theta\approx0.3$ in $2\leq d\leq 10$.

Dynamic marginality relates these exponents by the Wyart stability bound: $\gamma\geq1/(2+\theta)$; in MRJ packings, this is nearly saturated [1405.2087].

Quasicontacts, negligible for static stability, become dynamically relevant by providing the shortest paths for contact closure upon infinitesimal perturbations. In high dimensions, the abundance parameter $A_d$ of quasicontacts grows faster than that of contacts, altering prefactors in theoretical density estimates and dominating in certain Bravais lattice cases [1405.2087].

## 4. Structural Characteristics and Statistical Descriptors

MRJ packings display a suite of structural features:

- **Isostatic backbone and rattlers**: The jammed network consists of a mechanically rigid backbone interspersed by rattlers (locally unjammed particles), with the rattler fraction decreasing with dimension and anisotropy. Ideal MRJ states in the infinite system limit may be rattler-free [2311.06970, 1008.2982].
- **Order metrics**: The MRJ state minimizes translational and/or orientational order metrics—translational order parameter $\mathcal{T}$, bond-orientational order $Q_6$, or general ψ [1008.2982, 1101.1327].
- **Voronoi anticorrelations**: Two-point correlation functions of Minkowski functionals (volume, surface area, mean curvature) of Voronoi cells exhibit strong anticorrelations, consistent with the suppression of large-scale density fluctuations. Large Voronoi cells are statistically paired with small ones, sharply distinguishing MRJ from unjammed liquids or uncorrelated Poisson processes [1501.00593].
- **Packing spectrum**: MRJ packings anchor the lower boundary of the φ–order parameter plane among all jammed packings, with the continuous spectrum extending upwards to fully ordered crystals; the full set of jammed packings at varying density displays protocol dependence, except for the uniquely defined MRJ point [1101.1327, 1008.2982].

## 5. Generalization to Shape, Mixtures, and Geometries

The MRJ paradigm—and associated structural features—extends robustly:

- **Nonspherical shapes**: MRJ packing fractions and contact statistics depend drastically on shape deformation parameters (e.g., superdisk/superball exponent $p$). For superballs, $\phi_{MRJ}$ increases monotonically and nonanalytically with $|p-1|$, and jammed networks become hypostatic, supported by nongeneric local arrangements of contacts [1001.0423, 2111.14954].
- **Binary and polydisperse mixtures**: Binary hard-disk MRJ states in 2D display a hyperuniformity exponent $\alpha$ maximized (up to $\approx0.45$) for certain size ratios $\beta\approx1.4$ [2412.10883]. In 3D, polydispersity preserves hyperuniformity in the medium even if number density fluctuations are not hyperuniform [1104.4091].
- **High dimensions and Bravais lattices**: In $d\gg1$, quasicontacts significantly increase the theoretical prefactor for packing density, and in certain isostatic Bravais lattices, dominate the asymptotic scaling exponent [1405.2087].
- **Confinement**: Under plane or wall confinement, MRJ packings retain global hyperuniformity while exhibiting layering, heightened rattler fractions, and discontinuous transitions in response to geometric frustration [1512.08740].

## 6. Transport, Electromagnetic, and Physical Property Implications

The disordered hyperuniformity of MRJ packings engenders exceptional transport, diffusive, and electromagnetic properties compared to conventional amorphous or crystalline media:

- **Transport**: Rigorous upper/lower bounds and universal scaling laws for permeability and mean survival time in porous MRJ media show dramatically reduced values, reflecting the uniformity and smallness of the voids [1710.04268, 1608.01984]. The principal relaxation time and trapping constants approach those of periodic lattices as the local void geometry is regularized by jamming.
- **Wave propagation**: MRJ packings are nearly dissipationless for long-wavelength electromagnetic waves due to hyperuniformity—Im$\epsilon_{e}\to 0$ in the quasistatic limit. Stealthy hyperuniform variants could even yield finite-wavelength complete photonic gaps [1710.04268].
- **Topology of the void space**: MRJ packings possess a narrow, compact-support pore-size distribution set by local loops in the contact network, with void regions sharply bounded by the jamming and saturation constraints [1104.4091].

## 7. Algorithms, Protocols, and Open Challenges

Protocols for generating MRJ packings include event-driven compression (Lubachevsky–Stillinger), sequential (or adaptive) linear programming (Torquato–Jiao), and Monte Carlo adaptive-shrinking-cell methods for non-spherical particles. The protocol-dependence of jammed states outside the MRJ point has been demonstrated; only the MRJ state specified by variational minimization over order metrics is universally defined [1101.1327, 1008.2982, 1406.1535].

Critical slowing down complicates generation of large, truly strictly jammed, rattler-free MRJ packings. Ensuring exact jamming and hyperuniformity for $N\gg10^4$ remains a fundamental challenge, intimately connected to the saturation of collective rearrangements near the jamming threshold [1606.05227, 2311.06970].

Further, capturing the full spectral density—including the amplitude of small-$k$ linear scaling and the direct prediction of void-space statistics from first principles—remains an open theoretical frontier [1104.4093].

---

### References (arXiv ID)

- Marginal stability, quasicontacts, and power-law force/gap distributions: [1405.2087]
- Non-universality and variational definition of MRJ: [1101.1327]
- Hyperuniformity and universal QLR correlations: [2311.06970], [1008.2548], [1104.4093], [1001.0423], [2111.14954], [1406.1535], [1512.08740]
- Protocol-dependence and the φ–order map: [1101.1327], [1008.2982]
- Voronoi statistics and Minkowski correlations: [1501.00593]
- Transport and EM property implications: [1608.01984], [1710.04268]
- Rattlers, algorithmic obstacles, and scaling: [1606.05227], [2412.10883]
- Binary and polydisperse MRJ packings: [2412.10883], [1104.4091]
- Geometric-structure theory for nonspherical shapes: [1605.06716]
- Analytic RCP/MRJ density via isostaticity: [2201.04541]

Source: https://www.emergentmind.com/topics/maximally-random-jammed-packings