Maximally Random Jammed State Research Articles

Hyperuniformity of maximally random jammed packings of hyperspheres across spatial dimensions.

The maximally random jammed (MRJ) state is the most random (i.e., disordered) configuration of strictly jammed (mechanically rigid) nonoverlapping objects. MRJ packings are hyperuniform, meaning their long-wavelength density fluctuations are anomalously suppressed compared to typical disordered systems, i.e., their structure factors S(k) tend to zero as the wave number |k| tends to zero. Here we show that generating high-quality strictly jammed states for Euclidean space dimensions d=3,4, and 5 is of paramount importance in ensuring hyperuniformity and extracting precise values of the hyperuniformity exponent α>0 for MRJ states, defined by the power-law behavior of S(k)∼|k|^{α} in the limit |k|→0. Moreover, we show that for fixed d it is more difficult to ensure jamming as the particle number N increases, which results in packings that are nonhyperuniform. Free-volume theory arguments suggest that the ideal MRJ state does not contain rattlers, which act as defects in numerically generated packings. As d increases, we find that the fraction of rattlers decreases substantially. Our analysis of the largest truly jammed packings suggests that the ideal MRJ packings for all dimensions d≥3 are hyperuniform with α=d-2, implying the packings become more hyperuniform as d increases. The differences in α between MRJ packings and the recently proposed Manna-class random close packed (RCP) states, which were reported to have α=0.25 in d=3 and be nonhyperuniform (α=0) for d=4 and d=5, demonstrate the vivid distinctions between the large-scale structure of RCP and MRJ states in these dimensions. Our paper clarifies the importance of the link between true jamming and hyperuniformity and motivates the development of an algorithm to produce rattler-free three-dimensional MRJ packings.

Predicting maximally random jammed packing density of non-spherical hard particles via analytical continuation of fluid equation of state.

Dense packings of hard particles are useful models for condensed matters including crystalline and glassy state of solids, simple liquids, granular materials and composites. It is very challenging to devise predictive theories of random packings, due to the intrinsic non-equilibrium and non-local nature of the system. Here, we develop a formalism for accurately predicting the density (i.e., fraction of space covered by the particles) ηMRJ of the maximally random jammed (MRJ) packing state of a wide spectrum of congruent non-spherical hard particles in three-dimensional Euclidean space [Doublestruck R]3, via analytical continuation of the corresponding fluid equation of state (EOS). This formalism is based on the assumption that the fluid branch of the EOS can be analytically extended into the meta-stable region, which leads to a diverging pressure at the jamming point (i.e., the MRJ state). This allows us to estimate ηMRJ as the pole in the EOS, which can be expressed in terms of the virial coefficients encoding intrinsic local n-body packing information of the particles, and depending alone on particle shape. The accuracy of our formalism is verified using the hard sphere system and is subsequently applied to a wide spectrum of non-spherical shapes. The predictions are compared to numerical results whenever possible, and excellent agreements are found.

Hard convex lens-shaped particles: Characterization of dense disordered packings.

Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is "optimal" in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction ϕ_{MRJ}≃0.73 [G. Cinacchi and S. Torquato, Soft Matter 14, 8205 (2018)1744-683X10.1039/C8SM01519H]. This value is only a few percent lower than ϕ_{DKP}=0.76210⋯, the packing fraction of the corresponding densest-known crystalline (degenerate) packings [G. Cinacchi and S. Torquato, J. Chem. Phys. 143, 224506 (2015)JCPSA60021-960610.1063/1.4936938]. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method-based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows us to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated here can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of ten contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve.

Evolutions of packing properties of perfect cylinders under densification and crystallization.

Cylindrical particles are ubiquitous in nature and industry, and a cylinder is a representative shape of rod-like particles. However, the disordered packing results of cylinders in previous studies are quite inconsistent with each other. In this work, we obtain the MRJ (maximally random jammed) packings and the MDRPs (maximally dense random packings) of perfect cylinders with the aspect ratio (height/diameter) 0.2 ≤ w ≤ 6.0 using the ASC (adaptive shrinking cell) algorithm and the IMC (inverse Monte Carlo) method, respectively. The optimal aspect ratio corresponding to the maximal packing density is w = 0.9 in the MRJ state, while the value is w = 1.2 in the MDRP state. Then we investigate the evolutions of packing properties of perfect cylinders under densification and crystallization. We compare the different final packing states generated via the two methods with different compression rates and order constraints. In the densification procedure, we generate jammed and random packings of cylinders with various compression rates via the ASC and IMC method, respectively. When decreasing the compression rate, we find that the packing density increases but the optimal w remains the same in both methods. In the crystallization procedure, the order constraint in the IMC method is gradually released which means the degree of order in the packings is allowed to increase, and we find that the optimal w shifts from 1.2 to 0.9 while the packing density increases as well. Meanwhile, the random packings evolve into the jammed packings in the crystallization procedure which reflects the competition mechanism between the randomness and jamming. These results also indicate that the optimal w is solely related to the degree of order in the cylinder packings but not determined by the protocol or packing density. Furthermore, a uniform shape elongation effect on the random-packing densities of various shaped particles is found via a new proposed definition of the scaled aspect ratio. Finally, a rough linear relationship between the mean and standard deviation of the reduced Voronoi cell volumes is obtained only for the random packings. Our findings should lead to a better understanding toward the jammed and random packings and are helpful in guiding the granular material design.

The Perfect Glass Paradigm: Disordered Hyperuniform Glasses Down to Absolute Zero.

Rapid cooling of liquids below a certain temperature range can result in a transition to glassy states. The traditional understanding of glasses includes their thermodynamic metastability with respect to crystals. However, here we present specific examples of interactions that eliminate the possibilities of crystalline and quasicrystalline phases, while creating mechanically stable amorphous glasses down to absolute zero temperature. We show that this can be accomplished by introducing a new ideal state of matter called a “perfect glass”. A perfect glass represents a soft-interaction analog of the maximally random jammed (MRJ) packings of hard particles. These latter states can be regarded as the epitome of a glass since they are out of equilibrium, maximally disordered, hyperuniform, mechanically rigid with infinite bulk and shear moduli, and can never crystallize due to configuration-space trapping. Our model perfect glass utilizes two-, three-, and four-body soft interactions while simultaneously retaining the salient attributes of the MRJ state. These models constitute a theoretical proof of concept for perfect glasses and broaden our fundamental understanding of glass physics. A novel feature of equilibrium systems of identical particles interacting with the perfect-glass potential at positive temperature is that they have a non-relativistic speed of sound that is infinite.

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.

Jamming and crystallization in athermal polymer packings

Dense packings of chains of hard spheres possess characteristic features that do not have a counterpart in corresponding packings of monomeric spheres especially near the maximally random jammed (MRJ) state. From the modelling perspective the additional requirement that spheres keep their connectivity while maximizing the occupied volume fraction imposes severe constraints on generation algorithms of dense chain configurations. The extremely sluggish dynamics imposed by the uncrossability of chains precludes the use of deterministic or stochastic dynamics to generate all but dilute polymer packings. As a viable alternative, especially tailored chain-connectivity-altering Monte Carlo (MC) algorithms have been developed that bypass this kinetic hindrance and have actually been able to produce packings of hard-sphere chains in a volume fraction range spanning from infinite dilution up to the MRJ state. Such very dense athermal polymer packings share a number of structural features with packings of monomeric hard spheres, but also display unique characteristics due to the constraints imposed by connectivity. We give an overview of the most relevant results of our recent modeling work on packings of freely-jointed chains of tangent hard spheres about the MRJ state, local structure, chain dimensions and their scaling with density, topological constraints in the form of entanglements and knots, contact network at jamming, and entropically driven crystallization.

Contact network in nearly jammed disordered packings of hard-sphere chains

We present salient results of the analysis of the geometrical structure of a large fully equilibrated ensemble of nearly jammed packings of linear freely jointed chains of tangent hard spheres generated via extensive Monte Carlo simulations. In spite the expected differences due to chain connectivity, both the pair-correlation function and the contact network for chain packings are found to strongly resemble those in packings of monomeric hard spheres at the maximally random jammed (MRJ) state. A remarkable finding of the present work is the tendency of chains to form closed loops at the MRJ state as a consequence of chain collapse. Our simulations on disordered nearly jammed chain packings yield an average coordination number of 6, which fulfills the isostaticity condition and is in excellent agreement with the corresponding simulation [A. Donev, S. Torquato, and F. H. Stillinger, Phys. Rev. E 71, 011105 (2005)] and experimental [T. Aste, M. Saadatfar, and T. J. Senden, Phys. Rev. E 71, 061302 (2005)] findings for jammed packings of monatomic hard spheres. An exact correspondence between the statistical-mechanical ensembles of monomeric spheres and of hard-sphere chains offers insights regarding the structure and topology of the contact network of hard-sphere systems at the MRJ state.

The structure of random packings of freely jointed chains of tangent hard spheres

We analyze the structure of dense random packings of freely jointed chains of tangent hard spheres as a function of concentration (packing density) with particular emphasis placed on the behavior in the vicinity of their maximally random jammed (MRJ) state. Representative configurations over the whole density range are generated through extensive off-lattice Monte Carlo simulations on systems of average chain lengths ranging from N=12 to 1000 hard spheres. Several measures of order are used to quantitatively describe either local structure (sphere arrangements and bonded geometry) or global behavior (chain conformations and statistics). In addition, the employed measures are used to elucidate the effect of connectivity on structure, by comparing monatomic and chain assemblies of hard spheres at the MRJ state.

Dense and Nearly Jammed Random Packings of Freely Jointed Chains of Tangent Hard Spheres

Dense packings of freely jointed chains of tangent hard spheres are produced by a novel Monte Carlo method. Within statistical uncertainty, chains reach a maximally random jammed (MRJ) state at the same volume fraction as packings of single hard spheres. A structural analysis shows that as the MRJ state is approached (i) the radial distribution function for chains remains distinct from but approaches that of single hard sphere packings quite closely, (ii) chains undergo progressive collapse, and (iii) a small but increasing fraction of sites possess highly ordered first coordination shells.

Diversity of order and densities in jammed hard-particle packings.

Recently the conventional notion of random close packing has been supplanted by the more appropriate concept of the maximally random jammed (MRJ) state. This inevitably leads to the necessity of distinguishing the MRJ state among the entire collection of jammed packings. While the ideal method of addressing this question would be to enumerate and classify all possible jammed hard-sphere configurations, practical limitations prevent such a method from being employed. Instead, we generate numerically a large number of representative jammed hard-sphere configurations (primarily relying on a slight modification of the Lubachevsky-Stillinger algorithm to do so) and evaluate several commonly employed order metrics for each of these packings. Our investigation shows that, even in the large-system limit, jammed systems of hard spheres can be generated with a wide range of packing fractions from phi approximately 0.52 to the fcc limit (phi approximately 0.74). Moreover, at a fixed packing fraction, the variation in the order can be substantial, indicating that the density alone does not uniquely characterize a packing. Interestingly, each order metric evaluated yielded a relatively consistent estimate for the packing fraction of the maximally random jammed state (phi(MRJ) approximately 0.63). This estimate, however, is compromised by the weaknesses in the order metrics available, and we propose several guiding principles for future efforts to define more broadly applicable metrics.