Application of Edwards' statistical mechanics to high-dimensional jammed sphere packings

Jammed disks of two sizes and weights in a channel: Alternating sequences

Disks of two sizes and weights in alternating sequence are confined to a long and narrow channel. The axis of the channel is horizontal and its plane vertical. The channel is closed off by pistons that freeze jammed microstates out of loose disk configurations subject to moderate pressure, gravity, and random agitations. Disk sizes and channel width are such that under jamming no disk remains loose and all disks touch one wall. We present exact results for the characterization of jammed macrostates including volume and entropy. The rigorous analysis divides the disk sequences of jammed microstates into overlapping tiles from which we construct a small number of species of statistically interacting particles. Jammed macrostates depend on dimensionless control parameters inferred from ratios between measures of expansion work against the pistons, gravitational potential energy, and intensity of random agitations. These control parameters enter the configurational statistics via the activation energies prior to jamming of the particles. The range of disk weights naturally divides into regimes where qualitatively different features come into play. We sketch a path toward generalizations that include random sequences under a modified jamming protocol.

Adhesive loose packings of small dry particles

We explore adhesive loose packings of small dry spherical particles of micrometer size using 3D discrete-element simulations with adhesive contact mechanics and statistical ensemble theory. A dimensionless adhesion parameter (Ad) successfully combines the effects of particle velocities, sizes and the work of adhesion, identifying a universal regime of adhesive packings for Ad > 1. The structural properties of the packings in this regime are well described by an ensemble approach based on a coarse-grained volume function that includes the correlation between bulk and contact spheres. Our theoretical and numerical results predict: (i) an equation of state for adhesive loose packings that appear as a continuation from the frictionless random close packing (RCP) point in the jamming phase diagram and (ii) the existence of an asymptotic adhesive loose packing point at a coordination number Z = 2 and a packing fraction ϕ = 1/2(3). Our results highlight that adhesion leads to a universal packing regime at packing fractions much smaller than the random loose packing (RLP), which can be described within a statistical mechanical framework. We present a general phase diagram of jammed matter comprising frictionless, frictional, adhesive as well as non-spherical particles, providing a classification of packings in terms of their continuation from the spherical frictionless RCP.

Marginal stability in jammed packings: quasicontacts and weak contacts

Maximally random jammed (MRJ) sphere packing is a prototypical example of a system naturally poised at the margin between underconstraint and overconstraint. This marginal stability has traditionally been understood in terms of isostaticity, the equality of the number of mechanical contacts and the number of degrees of freedom. Quasicontacts, pairs of spheres on the verge of coming in contact, are irrelevant for static stability, but they come into play when considering dynamic stability, as does the distribution of contact forces. We show that the effects of marginal dynamic stability, as manifested in the distributions of quasicontacts and weak contacts, are consequential and nontrivial. We study these ideas first in the context of MRJ packing of d-dimensional spheres, where we show that the abundance of quasicontacts grows at a faster rate than that of contacts. We reexamine a calculation of Jin et al. [Phys. Rev. E 82, 051126 (2010)], where quasicontacts were originally neglected, and we explore the effect of their inclusion in the calculation. This analysis yields an estimate of the asymptotic behavior of the packing density in high dimensions. We argue that this estimate should be reinterpreted as a lower bound. The latter part of the paper is devoted to Bravais lattice packings that possess the minimum number of contacts to maintain mechanical stability. We show that quasicontacts play an even more important role in these packings. We also show that jammed lattices are a useful setting for studying the Edwards ensemble, which weights each mechanically stable configuration equally and does not account for dynamics. This ansatz fails to predict the power-law distribution of near-zero contact forces, P(f)∼f(θ).

Fundamental challenges in packing problems: from spherical to non-spherical particles

Random packings of objects of a particular shape are ubiquitous in science and engineering. However, such jammed matter states have eluded any systematic theoretical treatment due to the strong positional and orientational correlations involved. In recent years progress on a fundamental description of jammed matter could be made by starting from a constant volume ensemble in the spirit of conventional statistical mechanics. Recent work has shown that this approach, first introduced by S. F. Edwards more than two decades ago, can be cast into a predictive framework to calculate the packing fractions of both spherical and non-spherical particles.

Statistical theory of correlations in random packings of hard particles

A random packing of hard particles represents a fundamental model for granular matter. Despite its importance, analytical modeling of random packings remains difficult due to the existence of strong correlations which preclude the development of a simple theory. Here, we take inspiration from liquid theories for the n-particle angular correlation function to develop a formalism of random packings of hard particles from the bottom up. A progressive expansion into a shell of particles converges in the large layer limit under a Kirkwood-like approximation of higher-order correlations. We apply the formalism to hard disks and predict the density of two-dimensional random close packing (RCP), ϕ(rcp) = 0.85 ± 0.01, and random loose packing (RLP), ϕ(rlp) = 0.67 ± 0.01. Our theory also predicts a phase diagram and angular correlation functions that are in good agreement with experimental and numerical data.

Jammed lattice sphere packings

We generate and study an ensemble of isostatic jammed hard-sphere lattices. These lattices are obtained by compression of a periodic system with an adaptive unit cell containing a single sphere until the point of mechanical stability. We present detailed numerical data about the densities, pair correlations, force distributions, and structure factors of such lattices. We show that this model retains many of the crucial structural features of the classical hard-sphere model and propose it as a model for the jamming and glass transitions that enables exploration of much higher dimensions than are usually accessible.

Random perfect lattices and the sphere packing problem

Motivated by the search for best lattice sphere packings in Euclidean spaces of large dimensions we study randomly generated perfect lattices in moderately large dimensions (up to d=19 included). Perfect lattices are relevant in the solution of the problem of lattice sphere packing, because the best lattice packing is a perfect lattice and because they can be generated easily. Their number, however, grows superexponentially with the dimension, so to get an idea of their properties we propose to study a randomized version of the generating algorithm and to define a random ensemble with an effective temperature in a way reminiscent of a Monte Carlo simulation. We therefore study the distribution of packing fractions and kissing numbers of these ensembles and show how as the temperature is decreased the best known packers are easily recovered. We find that, even at infinite temperature, the typical perfect lattices are considerably denser than known families (like A(d) and D(d)), and we propose two hypotheses between which we cannot distinguish in this paper: one in which they improve the Minkowsky bound φ~2(-(0.84±0.06)d), and a competitor in which their packing fraction decreases superexponentially, namely, φ~d(-ad) but with a very small coefficient a=0.06±0.04. We also find properties of the random walk which are suggestive of a glassy system already for moderately small dimensions. We also analyze local structure of network of perfect lattices conjecturing that this is a scale-free network in all dimensions with constant scaling exponent 2.6±0.1.

Edwards thermodynamics of the jamming transition for frictionless packings: ergodicity test and role of angoricity and compactivity

This paper illustrates how the tools of equilibrium statistical mechanics can help to describe a far-from-equilibrium problem: the jamming transition in frictionless granular materials. Edwards ideas consist of proposing a statistical ensemble of volume and stress fluctuations through the thermodynamic notion of entropy, compactivity, X, and angoricity, A (two temperature-like variables). We find that Edwards thermodynamics is able to describe the jamming transition (J point) in frictionless packings. Using the ensemble formalism we elucidate the following: (i) We test the combined volume-stress ensemble by comparing the statistical properties of jammed configurations obtained by dynamics with those averaged over the ensemble of minima in the potential energy landscape as a test of ergodicity. Agreement between both methods supports the idea of ergodicity and "thermalization" at a given angoricity and compactivity. (ii) A microcanonical ensemble analysis supports the maximum entropy principle for grains. (iii) The intensive variables A and X describe the approach to jamming through a series of scaling relations as A → 0+ and X → 0-. Due to the force-strain coupling in the interparticle forces, the jamming transition is probed thermodynamically by a "jamming temperature" T(J) composed of contributions from A and X. (iv) The thermodynamic framework reveals the order of the jamming phase transition by showing the absence of critical fluctuations at jamming in static observables like pressure and volume, and we discuss other critical scenarios for the jamming transition. (v) Finally, we elaborate on a comparison with relevant studies by Gao, Blawzdziewicz, and O'Hern [Phys. Rev. E 74, 061304 (2006)], showing a breakdown of equiprobability of microstates obtained via fast quenches. A network analysis of the energy landscape reveals the origin of the inhomogeneities in the uneven distribution of the areas of the basins. Such inhomogeneities are also found in other out-of-equilibrium systems like Lennard-Jones glasses and their existence does not preclude the use of statistical mechanics for jammed systems.

Glass transition and random close packing above three dimensions

Motivated by a recently identified severe discrepancy between a static and a dynamic theory of glasses, we numerically investigate the behavior of dense hard spheres in spatial dimensions 3 to 12. Our results are consistent with the static replica theory, but disagree with the dynamic mode-coupling theory, indicating that key ingredients of high-dimensional physics are missing from the latter. We also obtain numerical estimates of the random close packing density, which provides new insights into the mathematical problem of packing spheres in large dimensions.