Universal microstructure and mechanical stability of jammed packings

Critical slowing down and hyperuniformity on approach to jamming

Hyperuniformity characterizes a state of matter that is poised at a critical point at which density or volume-fraction fluctuations are anomalously suppressed at infinite wavelengths. Recently, much attention has been given to the link between strict jamming (mechanical rigidity) and (effective or exact) hyperuniformity in frictionless hard-particle packings. However, in doing so, one must necessarily study very large packings in order to access the long-ranged behavior and to ensure that the packings are truly jammed. We modify the rigorous linear programming method of Donev et al. [J. Comput. Phys. 197, 139 (2004)JCTPAH0021-999110.1016/j.jcp.2003.11.022] in order to test for jamming in putatively collectively and strictly jammed packings of hard disks in two dimensions. We show that this rigorous jamming test is superior to standard ways to ascertain jamming, including the so-called "pressure-leak" test. We find that various standard packing protocols struggle to reliably create packings that are jammed for even modest system sizes of N≈10^{3} bidisperse disks in two dimensions; importantly, these packings have a high reduced pressure that persists over extended amounts of time, meaning that they appear to be jammed by conventional tests, though rigorous jamming tests reveal that they are not. We present evidence that suggests that deviations from hyperuniformity in putative maximally random jammed (MRJ) packings can in part be explained by a shortcoming of the numerical protocols to generate exactly jammed configurations as a result of a type of "critical slowing down" as the packing's collective rearrangements in configuration space become locally confined by high-dimensional "bottlenecks" from which escape is a rare event. Additionally, various protocols are able to produce packings exhibiting hyperuniformity to different extents, but this is because certain protocols are better able to approach exactly jammed configurations. Nonetheless, while one should not generally expect exact hyperuniformity for disordered packings with rattlers, we find that when jamming is ensured, our packings are very nearly hyperuniform, and deviations from hyperuniformity correlate with an inability to ensure jamming, suggesting that strict jamming and hyperuniformity are indeed linked. This raises the possibility that the ideal MRJ packings have no rattlers. Our work provides the impetus for the development of packing algorithms that produce large disordered strictly jammed packings that are rattler free, which is an outstanding, challenging task.

Pinning Susceptibility: The Effect of Dilute, Quenched Disorder on Jamming

We study the effect of dilute pinning on the jamming transition. Pinning reduces the average contact number needed to jam unpinned particles and shifts the jamming threshold to lower densities, leading to a pinning susceptibility, χ_{p}. Our main results are that this susceptibility obeys scaling form and diverges in the thermodynamic limit as χ_{p}∝|ϕ-ϕ_{c}^{∞}|^{-γ_{p}} where ϕ_{c}^{∞} is the jamming threshold in the absence of pins. Finite-size scaling arguments yield these values with associated statistical (systematic) errors γ_{p}=1.018±0.026(0.291) in d=2 and γ_{p}=1.534±0.120(0.822) in d=3. Logarithmic corrections raise the exponent in d=2 to close to the d=3 value, although the systematic errors are very large.

The random packing density of nearly spherical particles

Obtaining general relations between macroscopic properties of random assemblies, such as density, and the microscopic properties of their constituent particles, such as shape, is a foundational challenge in the study of amorphous materials. By leveraging existing understanding of the random packing of spherical particles, we estimate the random packing density for all sufficiently spherical shapes. Our method uses the ensemble of random packing configurations of spheres as a reference point for a perturbative calculation, which we carry to linear order in the deformation. A fully analytic calculation shows that all sufficiently spherical shapes pack more densely than spheres. Additionally, we use simulation data for spheres to calculate numerical estimates for nonspherical particles and compare these estimates to simulations.

Geometric order parameters derived from the Voronoi tessellation show signatures of the jamming transition

A jammed packing of frictionless spheres at zero temperature is perfectly specified by the network of contact forces from which mechanical properties can be derived. However, we can alternatively consider a packing as a geometric structure, characterized by a Voronoi tessellation which encodes the local environment around each particle. We find that this local environment characterizes systems both above and below jamming and changes markedly at the transition. A variety of order parameters derived from this tessellation carry signatures of the jamming transition, complete with scaling exponents. Furthermore, we define a real space geometric correlation function which also displays a signature of jamming. Taken together, these results demonstrate the validity and usefulness of a purely geometric approach to jamming.

Quantitative approximation schemes for glasses

By means of a systematic expansion around the infinite-dimensional solution, we obtain an approximation scheme to compute properties of glasses in low dimensions. The resulting equations take as input the thermodynamic and structural properties of the equilibrium liquid, and from this they allow one to compute properties of the glass. They are therefore similar in spirit to the Mode Coupling approximation scheme. Our scheme becomes exact, by construction, in dimension d→∞, and it can be improved systematically by adding more terms in the expansion.

Scaling collapse at the jamming transition

The jamming transition of particles with finite-range interactions is characterized by a variety of critical phenomena, including power-law distributions of marginal contacts. We numerically study a recently proposed simple model of jamming, which is conjectured to lie in the same universality class as the jamming of spheres in all dimensions. We extract numerical estimates of the critical exponents, θ=0.451±0.006 and γ=0.404±0.004, that match the exponents observed in sphere packing systems. We analyze finite-size scaling effects that manifest in a subcritical cutoff regime and size-independent but protocol-dependent scaling curves. Our results support the conjectured link with sphere jamming, provide more precise measurements of the critical exponents than previously reported, and shed light on the finite-size scaling behavior of continuous constraint satisfiability transitions.

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.

Structural and mechanical features of the order-disorder transition in experimental hard-sphere packings

Here we present an experimental and numerical investigation on the grain-scale geometrical and mechanical properties of partially crystallized structures made of macroscopic frictional grains. Crystallization is inevitable in arrangements of monosized hard spheres with packing densities exceeding Bernal's limiting density ϕ(Bernal)≈0.64. We study packings of monosized hard spheres whose density spans over a wide range (0.59<ϕ<0.72). These experiments harness x-ray computed tomography, three-dimensional image analysis, and numerical simulations to access precisely the geometry and the 3D structure of internal forces within the sphere packings. We show that clear geometrical transitions coincide with modifications of the mechanical backbone of the packing both at the grain and global scale. Notably, two transitions are identified at ϕ(Bernal)≈0.64 and ϕ(c)≈0.68. These results provide insights on how geometrical and mechanical features at the grain scale conspire to yield partially crystallized structures that are mechanically stable.

Jamming criticality revealed by removing localized buckling excitations

Recent theoretical advances offer an exact, first-principles theory of jamming criticality in infinite dimension as well as universal scaling relations between critical exponents in all dimensions. For packings of frictionless spheres near the jamming transition, these advances predict that nontrivial power-law exponents characterize the critical distribution of (i) small interparticle gaps and (ii) weak contact forces, both of which are crucial for mechanical stability. The scaling of the interparticle gaps is known to be constant in all spatial dimensions d-including the physically relevant d=2 and 3, but the value of the weak force exponent remains the object of debate and confusion. Here, we resolve this ambiguity by numerical simulations. We construct isostatic jammed packings with extremely high accuracy, and introduce a simple criterion to separate the contribution of particles that give rise to localized buckling excitations, i.e., bucklers, from the others. This analysis reveals the remarkable dimensional robustness of mean-field marginality and its associated criticality.

A master equation for the probability distribution functions of forces in soft particle packings

We study the microscopic response of force-chain networks in jammed soft particles to quasi-static isotropic (de)compressions by molecular dynamics simulations. We show that not only contacts but also interparticle gaps between the nearest neighbors must be considered for the stochastic evolution of the probability distribution functions (PDFs) of forces, where the mutual exchange of contacts and interparticle gaps, i.e. opening and closing contacts, are also crucial to the incremental system behavior. By numerically determining the transition rates for all changes of contacts and gaps, we formulate a Master equation for the PDFs of forces, where the insight one gets from the transition rates is striking: the mean change of forces reflects non-affine system responses, while their fluctuations obey uncorrelated Gaussian statistics. In contrast, interparticle gaps react mostly affine in average, but imply multi-scale correlations according to a much wider stable distribution function.

Devising a protocol-related statistical mechanics framework for granular materials

Devising a statistical mechanics framework for jammed granular materials is a challenging task as those systems do not share some important properties required to characterize them with statistical thermodynamics tools. In a recent paper [Asenjo et al. Phys. Rev. Lett. 112, 098002 (2014)], a new definition of a granular entropy, which puts the protocol used to generate the packings at its roots, has been proposed. Following up these results, it is shown that the protocol used in Asenjo et al. can be recast as a canonical ensemble with a particular value of the temperature. Signature of gaussianity for large system sizes strongly suggests an asymptotic equivalence with a corresponding microcanonical ensemble where jammed states with certain basin volumes are sampled uniformly. We argue that this microcanonical ensemble is not Edwards's microcanonical ensemble and generalize this argument to other protocols.

Existence of isostatic, maximally random jammed monodisperse hard-disk packings

We generate jammed packings of monodisperse circular hard-disks in two dimensions using the Torquato-Jiao sequential linear programming algorithm. The packings display a wide diversity of packing fractions, average coordination numbers, and order as measured by standard scalar order metrics. This geometric-structure approach enables us to show the existence of relatively large maximally random jammed (MRJ) packings with exactly isostatic jammed backbones and a packing fraction (including rattlers) of [Formula: see text]. By contrast, the concept of random close packing (RCP) that identifies the most probable packings as the most disordered misleadingly identifies highly ordered disk packings as RCP in 2D. Fundamental structural descriptors such as the pair correlation function, structure factor, and Voronoi statistics show a strong contrast between the MRJ state and the typical hyperstatic, polycrystalline packings with [Formula: see text] that are more commonly obtained using standard packing protocols. Establishing that the MRJ state for monodisperse hard disks is isostatic and qualitatively distinct from commonly observed polycrystalline packings contradicts conventional wisdom that such a disordered, isostatic packing does not exist due to a lack of geometrical frustration and sheds light on the nature of disorder. This prompts the question of whether an algorithm may be designed that is strongly biased toward generating the monodisperse disk MRJ state.

Force distribution affects vibrational properties in hard-sphere glasses

We theoretically and numerically study the elastic properties of hard-sphere glasses and provide a real-space description of their mechanical stability. In contrast to repulsive particles at zero temperature, we argue that the presence of certain pairs of particles interacting with a small force f soften elastic properties. This softening affects the exponents characterizing elasticity at high pressure, leading to experimentally testable predictions. Denoting P(f) ~ f(θ(e)), the force distribution of such pairs and ϕ(c) the packing fraction at which pressure diverges, we predict that (i) the density of states has a low-frequency peak at a scale ω*, rising up to it as D(ω) ~ ω(2+a), and decaying above ω* as D(ω) ~ ω(-a) where a = (1 - θ(e))/(3 + θ(e)) and ω is the frequency, (ii) shear modulus and mean-squared displacement are inversely proportional with ⟨δR²⟩ ~ 1/μ ~ (ϕ(c) - ϕ)(κ), where κ = 2 - 2/(3 + θ(e)), and (iii) continuum elasticity breaks down on a scale ℓ(c) ~ 1/√(δz) ~ (ϕ(c) - ϕ)(-b), where b = (1 + θ(e))/(6 + 2θ(e)) and δz = z - 2d, where z is the coordination and d the spatial dimension. We numerically test (i) and provide data supporting that θ(e) ≈ 0.41 in our bidisperse system, independently of system preparation in two and three dimensions, leading to κ ≈ 1.41, a ≈ 0.17, and b ≈ 0.21. Our results for the mean-square displacement are consistent with a recent exact replica computation for d = ∞, whereas some observations differ, as rationalized by the present approach.

Cavity method for force transmission in jammed disordered packings of hard particles

The force distribution of jammed disordered packings has always been considered a central object in the physics of granular materials. However, many of its features are poorly understood. In particular, analytic relations to other key macroscopic properties of jammed matter, such as the contact network and its coordination number, are still lacking. Here we develop a mean-field theory for this problem, based on the consideration of the contact network as a random graph where the force transmission becomes a constraint satisfaction problem. We can thus use the cavity method developed in the past few decades within the statistical physics of spin glasses and hard computer science problems. This method allows us to compute the force distribution P(f) for random packings of hard particles of any shape, with or without friction. We find a new signature of jamming in the small force behavior P(f) ∼ f(θ), whose exponent has attracted recent active interest: we find a finite value for P(f = 0), along with θ = 0. Furthermore, we relate the force distribution to a lower bound of the average coordination number z[combining macron](μ) of jammed packings of frictional spheres with coefficient μ. This bridges the gap between the two known isostatic limits z[combining macron]c (μ = 0) = 2D (in dimension D) and z[combining macron]c(μ → ∞) = D + 1 by extending the naive Maxwell's counting argument to frictional spheres. The theoretical framework describes different types of systems, such as non-spherical objects in arbitrary dimensions, providing a common mean-field scenario to investigate force transmission, contact networks and coordination numbers of jammed disordered packings.

Mechanical characterization of partially crystallized sphere packings

We study grain-scale mechanical and geometrical features of partially crystallized packings of frictional spheres, produced experimentally by a vibrational protocol. By combining x-ray computed tomography, 3D image analysis, and discrete element method simulations, we have access to the 3D structure of internal forces. We investigate how the network of mechanical contacts and intergranular forces change when the packing structure evolves from amorphous to near perfect crystalline arrangements. We compare the behavior of the geometrical neighbors (quasicontracts) of a grain to the evolution of the mechanical contacts. The mechanical coordination number Z(m) is a key parameter characterizing the crystallization onset. The high fluctuation level of Z(m) and of the force distribution in highly crystallized packings reveals that a geometrically ordered structure still possesses a highly random mechanical backbone similar to that of amorphous packings.

Jamming in finite systems: stability, anisotropy, fluctuations, and scaling

Athermal packings of soft repulsive spheres exhibit a sharp jamming transition in the thermodynamic limit. Upon further compression, various structural and mechanical properties display clean power-law behavior over many decades in pressure. As with any phase transition, the rounding of such behavior in finite systems close to the transition plays an important role in understanding the nature of the transition itself. The situation for jamming is surprisingly rich: the assumption that jammed packings are isotropic is only strictly true in the large-size limit, and finite-size has a profound effect on the very meaning of jamming. Here, we provide a comprehensive numerical study of finite-size effects in sphere packings above the jamming transition, focusing on stability as well as the scaling of the contact number and the elastic response.

Shear modulus of glasses: results from the full replica-symmetry-breaking solution

We compute the shear modulus of amorphous hard and soft spheres, using the exact solution in infinite spatial dimensions that has been developed recently. We characterize the behavior of this observable in the whole phase diagram, and in particular around the glass and jamming transitions. Our results are consistent with other theoretical approaches, which are unified within this general picture, and they are also consistent with numerical and experimental results. Furthermore, we discuss some properties of the out-of-equilibrium dynamics after a deep quench close to the jamming transition, and we show that a combined measure of the shear modulus and of the mean square displacement allows one to probe experimentally the complex structure of phase space predicted by the full replica-symmetry-breaking solution.

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(θ).

Physics in ordered and disordered colloidal matter composed of poly(N-isopropylacrylamide) microgel particles

This review collects and describes experiments that employ colloidal suspensions to probe physics in ordered and disordered solids and related complex fluids. The unifying feature of this body of work is its clever usage of poly(N-isopropylacrylamide) (PNIPAM) microgel particles. These temperature-sensitive colloidal particles provide experimenters with a 'knob' for in situ control of particle size, particle interaction and particle packing fraction that, in turn, influence the structural and dynamical behavior of the complex fluids and solids. A brief summary of PNIPAM particle synthesis and properties is given, followed by a synopsis of current activity in the field. The latter discussion describes a variety of soft matter investigations including those that explore formation and melting of crystals and clusters, and those that probe structure, rearrangement and rheology of disordered (jammed/glassy) and partially ordered matter. The review, therefore, provides a snapshot of a broad range of physics phenomenology which benefits from the unique properties of responsive microgel particles.

Fractal free energy landscapes in structural glasses

Glasses are amorphous solids whose constituent particles are caged by their neighbours and thus cannot flow. This sluggishness is often ascribed to the free energy landscape containing multiple minima (basins) separated by high barriers. Here we show, using theory and numerical simulation, that the landscape is much rougher than is classically assumed. Deep in the glass, it undergoes a 'roughness transition' to fractal basins, which brings about isostaticity and marginal stability on approaching jamming. Critical exponents for the basin width, the weak force distribution and the spatial spread of quasi-contacts near jamming can be analytically determined. Their value is found to be compatible with numerical observations. This advance incorporates the jamming transition of granular materials into the framework of glass theory. Because temperature and pressure control what features of the landscape are experienced, glass mechanics and transport are expected to reflect the features of the topology we discuss here.

Geometric signatures of jamming in the mechanical vacuum

Jamming has traditionally been studied as a mechanical phenomenon and characterized with mechanical order parameters. However, this approach is not meaningful in the "mechanical vacuum" of systems below jamming in which all mechanical properties are precisely zero. We find that the network of nearest neighbors and the geometric structure of the Voronoi cell contain well-defined and meaningful order parameters for jamming, which exist on both sides of the transition. We observe critical exponents in these order parameters and an upper critical dimension of 3. Further, we present evidence for a new incipient-jamming phase below the jamming transition marked by additional symmetry in the Voronoi tessellation.

Jamming graphs: a local approach to global mechanical rigidity

We revisit the concept of minimal rigidity as applied to frictionless, repulsive soft sphere packings in two dimensions with the introduction of the jamming graph. Minimal rigidity is a purely combinatorial property encoded via Laman's theorem in two dimensions. It constrains the global, average coordination number of the graph, for example. However, minimal rigidity does not address the geometry of local mechanical stability. The jamming graph contains both properties of global mechanical stability at the onset of jamming and local mechanical stability. We demonstrate how jamming graphs can be constructed using local moves via the Henneberg construction such that these graphs fall under the jurisdiction of correlated percolation. We then probe how jamming graphs destabilize, or become unjammed, by deleting a bond and computing the resulting rigid cluster distribution. We also study how the system restabilizes with the addition of new contacts and how a jamming graph with extra (redundant) contacts destabilizes. The latter endeavor allows us to probe a disk packing in the rigid phase and uncover a potentially new diverging length scale associated with the random deletion of contacts as compared to the study of cut-out (or frozen-in) subsystems.

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.

Detailed characterization of rattlers in exactly isostatic, strictly jammed sphere packings

We generate jammed disordered packings of 100≤N≤2000 monodisperse hard spheres in three dimensions whose strictly jammed backbones are demonstrated to be exactly isostatic with unprecedented numerical accuracy. This is accomplished by using the Torquato-Jiao (TJ) packing algorithm as a means of studying the maximally random jammed (MRJ) state. The rattler fraction of these packings converges towards 0.015 in the infinite-system limit, which is markedly lower than previous estimates for the MRJ state using the Lubachevsky-Stillinger protocol. This is because the packings that the TJ algorithm creates are closer to the true MRJ state, as shown using bond-orientational and translational order metrics. The rattler pair correlation statistics exhibit strongly correlated behavior contrary to the conventional understanding that they be randomly (Poisson) distributed. Dynamically interacting "polyrattlers" may be found imprisoned in shared cages as well as interacting through "bottlenecks" in the backbone and these clusters are mainly responsible for the sharp increase in the rattler pair correlation function near contact. We discover the surprising existence of polyrattlers with cluster sizes of up to five rattlers (which is expected to increase with system size) and present a distribution of polyrattler occurrence as a function of cluster size and system size. We also enumerate all of the rattler interaction topologies we observe and present images of several examples, showing that MRJ packings of monodisperse spheres can contain large rattler cages while still obeying the strict jamming criterion. The backbone spheres that encage the rattlers are significantly hypostatic, implying that correspondingly hyperstatic regions must exist elsewhere in these isostatic packings. We also observe that rattlers in hard-sphere packings share an apparent connection with the low-temperature two-level system anomalies that appear in real amorphous insulators and semiconductors.

Bond percolation in higher dimensions

We collect results for bond percolation on various lattices from two to fourteen dimensions that, in the limit of large dimension d or number of neighbors z, smoothly approach a randomly diluted Erdős-Rényi graph. We include results on bond-diluted hypersphere packs in up to nine dimensions, which show the mean coordination, excess kurtosis, and skewness evolving smoothly with dimension towards the Erdős-Rényi limit.