Numerical Simulation of Random Close Packings in Particle Deformation from Spheres to Cubes

Variation of packing density in particle deforming from spheres to cubes is studied. A new model is presented to describe particle deformation between different particle shapes. Deformation is simulated by relative motion of component spheres in the sphere assembly model of a particle. Random close packings of particles in deformation form spheres to cubes are simulated with an improved relaxation algorithm. Packings in both 2D and 3D cases are simulated. With the simulations, we find that the packing density increases while the particle sphericity decreases in the deformation. Spheres and cubes give the minimum (0.6404) and maximum (0.7755) of packing density in the deformation respectively. In each deforming step, packings starting from a random configuration and from the final packing of last deforming step are both simulated. The packing density in the latter case is larger than the former in two dimensions, but is smaller in three dimensions. The deformation model can be applied to other particle shapes as well.     

Polytetrahedral nature of the dense disordered packings of hard spheres

We study the structure of numerically simulated hard sphere packings at different densities by investigating local tetrahedral configurations of the spheres. Clusters of tetrahedra adjacent by faces present relatively dense aggregates of spheres atypical for crystals. The number of spheres participating in such polytetrahedral configurations increases with densification of the packing, and at the Bernal's limiting density (the packing fraction around 0.64) all spheres of the packing become involved in such tetrahedra. Thus the polytetrahedral packing cannot provide further increase in the density, and alternative structural change (formation of crystalline nuclei) begins henceforth.     

直径任意分布球填充的数值模拟

提出球填充数值算法的新分类方法.改进原有的松弛算法,使其能够模拟直径任意分布的球填充问题,采用可变循环周期使不同球数情形下的填充率基本保持不变.算例数据表明,该算法的填充率和配位数均高于原算法.由于采用背景网格搜索和双向链表组数据结构,使得邻接球搜索效率有相当大的提高,算法的时间复杂度为O(N)(N为球数).在一台AMD Athlon 3200+PC上,对于10000个等径球的随机密排列,只需217s,填充率即可达到0.64.     

Boerdijk-Coxeter helix and biological helices

Helices and dense packing of spherical objects are two closely related problems. For instance, the Boerdijk-Coxeter helix, which is obtained as a linear packing of regular tetrahedra, is a very efficient solution to some close-packing problems. The shapes of biological helices result from various kinds of interaction forces, including steric repulsion. Thus, the search for a maximum density can lead to structures related to the Boerdijk-Coxeter helix. Examples are presented for the -helix structure in proteins and for the structure of the protein collagen, but there are other examples of helical packings at different scales in biology. Models based on packing efficiency related to the Boerdijk-Coxeter helix, explain, mainly from topological arguments, why the number of amino acids per turn is close to 3.6 in -helices and 2.7 in collagen. Received 26 November 1998 and Received in final form 12 April 1999     

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.     

Degenerate quasicrystal of hard triangular bipyramids

We report a degenerate quasicrystal in Monte Carlo simulations of hard triangular bipyramids each composed of two regular tetrahedra sharing a single face. The dodecagonal quasicrystal is similar to that recently reported for hard tetrahedra [Haji-Akbari et al., Nature (London) 462, 773 (2009)] but degenerate in the pairing of tetrahedra, and self-assembles at packing fractions above 54%. This notion of degeneracy differs from the degeneracy of a quasiperiodic random tiling arising through phason flips. Free energy calculations show that a triclinic crystal is preferred at high packing fractions.     

Dense packing of binary and polydisperse hard spheres

The packing of binary and polydisperse unimodal and bimodal ensembles of hard spheres in the limit of high pressure is studied using a sequential addition algorithm. Upon fixing the number of particles, and their size distribution, the average (maximum) packing fraction is determined for systems of up to 20 000 particles. The structures obtained correspond to amorphous states close to the dense random close packing density. Binary distributions obtained are denser than the equivalent monodisperse distribution and agree with the theoretical prediction for an infinite size ratio limit. Unimodal normal and lognormal polydisperse distributions obtained compare favourably with available simulation and experimental data. Results for bimodal lognormal distributions are presented. In all cases it is seen how an increase in polydispersity increases the packing fraction of the system. The results can be employed to gain insight into optimal formulations for dense emulsions.     

Micromechanical simulation and analysis of one-dimensional vibratory sphere packing

We present a numerical method capable of reproducing the densification process from the so-called random loose to dense packing of uniform spheres under vertical vibration. The effects of vibration amplitude and frequency are quantified, and the random close packing is shown to be achieved only if both parameters are properly controlled. Two densification mechanisms are identified: pushing filling by which the contact between spheres is maintained and jumping filling by which the contact between particles is periodically broken. In general, pushing filling occurs when the vibration intensity is low and jumping filling becomes dominant when the vibration intensity is high.     

受振颗粒“毛细”系统中的对流与有序化

将球形颗粒倒入内径较窄的管状容器时,管壁的曲率会对颗粒的堆积结构产生影响,存在壁效应. 实验表明通过连续的竖直方向的振动,壁效应可以被强化,颗粒可以经由对流由无序排列转变为稳定的同轴筒形“壳层”结构. 每一壳层内,颗粒是二维的六角密堆积, 在这一转变过程中,颗粒堆积率的径向分布由初始的衰减振荡转变为等幅振荡. 分析了堆积率的不均匀性及空气在对流中的作用,以及形成“壳层”结构的动力学过程,对“壳层”结构的稳定性亦进行了讨论. 关键词： 颗粒物质 对流 有序化 毛细     

Is random close packing of spheres well defined?

Despite its long history, there are many fundamental issues concerning random packings of spheres that remain elusive, including a precise definition of random close packing (RCP). We argue that the current picture of RCP cannot be made mathematically precise and support this conclusion via a molecular dynamics study of hard spheres using the Lubachevsky-Stillinger compression algorithm. We suggest that this impasse can be broken by introducing the new concept of a maximally random jammed state, which can be made precise.     

Jamming II: Edwards’ statistical mechanics of random packings of hard spheres

The problem of finding the most efficient way to pack spheres has an illustrious history, dating back to the crystalline arrays conjectured by Kepler and the random geometries explored by Bernal in the 1960s. This problem finds applications spanning from the mathematician’s pencil, the processing of granular materials, the jamming and glass transitions, all the way to fruit packing in every grocery. There are presently numerous experiments showing that the loosest way to pack spheres gives a density of ∼55% (named random loose packing, RLP) while filling all the loose voids results in a maximum density of ∼63%-64% (named random close packing, RCP). While those values seem robustly true, to this date there is no well-accepted physical explanation or theoretical prediction for them. Here we develop a common framework for understanding the random packings of monodisperse hard spheres whose limits can be interpreted as the experimentally observed RLP and RCP. The reason for these limits arises from a statistical picture of jammed states in which the RCP can be interpreted as the ground state of the ensemble of jammed matter with zero compactivity, while the RLP arises in the infinite compactivity limit. We combine an extended statistical mechanics approach ‘a la Edwards’ (where the role traditionally played by the energy and temperature in thermal systems is substituted by the volume and compactivity) with a constraint on mechanical stability imposed by the isostatic condition. We show how such approaches can bring results that can be compared to experiments and allow for an exploitation of the statistical mechanics framework. The key result is the use of a relation between the local Voronoi volumes of the constituent grains (denoted the volume function) and the number of neighbors in contact that permits us to simply combine the two approaches to develop a theory of volume fluctuations in jammed matter. Ultimately, our results lead to a phase diagram that provides a unifying view of the disordered hard sphere packing problem and further sheds light on a diverse spectrum of data, including the RLP state. Theoretical results are well reproduced by numerical simulations that confirm the essential role played by friction in determining both the RLP and RCP limits. The RLP values depend on friction, explaining why varied experimental results can be obtained.     

Mean-field cage theory for the random close packed state of a metastable hard-sphere glass

In 2005, we developed a mean-field cage theory for the freezing of a stable hard-sphere fluid using the character of a stable hard-sphere fluid, some observations and the mean-configuration approximation [X.Z. Wang, J. Chem. Phys. 122 (2005) 044515]. It was found that near the freezing point, a thermal fluctuation of a cage causes the hard sphere in this cage to exchange positions with one of its nearest neighbors. In this paper, we extend the theory to the random close packed state of a metastable hard-sphere glass. It is found that near the random close packing point, a thermal fluctuation of a cage sets three of the hard spheres in this cage and its nearest cages into the local circulatory motion, resulting in indirect position exchanges among these three hard spheres. We obtain an analytical formula for the random close packing density. The predicted values are in good agreement with the experimental and simulation results for spatial dimensions d=2–7$d=2\text{–}7$.     

Dense packing in the monodisperse hard-sphere system: A numerical study

We report a numerical study of the close packing of monodisperse hard spheres. The close packings of hard spheres are produced by the Lubachesky-Stillinger (LS) compression algorithm and span the range from the disordered states to the ordered states. We provide quantitative evidence for the claim that the density and structural order of the arrested close packing can be determined by the compression rate, i.e., with slower rates producing denser and more ordered structures. Through deeply analyzing the structure of the resulting arrested close packings, a transition region has been identified in the plane of density and reciprocal compression rate, in between what have been historically thought of as amorphous and crystalline packings. We also find clear system size dependences in studying the structural properties of the packings from the disordered ones to the ordered ones. These detailed investigations, on the structure of the arrested close packings, may provide a link between the glassy states and the crystalline states in the hard spheres.     

Spectral properties of contact matrix: Application to proteins

A protein can be modelled by a set of points representing its amino acids. Topologically, this set of points is entirely defined by its contact matrix (adjacency matrix in graph theory). The contact matrix characterizing the relation between neighboring amino acids is deduced from Voronoi or Laguerre decomposition. This method allows contact matrices to be defined without any arbitrary cut-off that could induce arbitrary effects. Eigenvalues of these matrices are related with elementary excitations in proteins. We present some spectral properties of these matrices that reflect global properties of proteins. The eigenvectors indicate participation of each amino acids to the excitation modes of the proteins. It is interesting to compare the protein modelled as a close packing of amino acids, with a random close packing of spheres. The main features of the protein are those of a packing, a result that confirms the importance of the dense packing model for proteins. Nevertheless there are some properties, specific to the hierarchical organization of the protein: the primary chain order, the secondary structures and the domain structures.     

Diffusive Wave Spectroscopy of a random close packing of spheres

We are interested in the propagation of light in a random packing of dielectric spheres within the geometrical optics approximation. Numerical simulations are performed using a ray tracing algorithm. The effective refractive indexes and the transport mean free path are computed for different refractive indexes of spheres and intersticial media. The variations of the optical path length under small deformations of the spheres assembly are also computed and compared to the results of Diffusive Wave Spectroscopy experiments. Finally, we propose a measure of the transport mean free path and a Diffusive Wave Spectroscopy experiment on a packing of glass spheres. The results of those experiments agree with the predictions of this ray tracing approach.     

Jamming III: Characterizing randomness via the entropy of jammed matter

The nature of randomness in disordered packings of frictional and frictionless spheres is investigated using theory and simulations of identical spherical grains. The entropy of the packings is defined through the force and volume ensemble of jammed matter and this is shown to be difficult to calculate analytically. A mesoscopic ensemble of isostatic states is then utilized in an effort to predict the entropy through the definition of a volume function that is dependent on the coordination number. Equations of state are obtained relating entropy, volume fraction and compactivity characterizing the different states of jammed matter, and elucidating the phase diagram for jammed granular matter. Analytical calculations are compared to numerical simulations using volume fluctuation analysis and graph theoretical methods, with reasonable agreement. The entropy of the jammed system reveals that random loose packings are more disordered than random close packings, allowing for an unambiguous interpretation of both limits. Ensemble calculations show that the entropy vanishes at random close packing (RCP), while numerical simulations show that a finite entropy remains in the microscopic states at RCP. The notion of a negative compactivity, which explores states with volume fractions below those achievable by existing simulation protocols, is also explored, expanding the equations of state. The mesoscopic theory reproduces the simulations results in shape well, though a difference in magnitude implies that the entire entropy of the packing may not be captured by the methods presented herein. We discuss possible extensions to the present mesoscopic approach describing packings from random loose packing (RLP) to RCP to the ordered branch of the equation of state in an effort to understand the entropy of jammed matter in the full range of densities from RLP to face-centered cubic (FCC) packing.     

Definition and properties of ideal amorphous solids

It is proposed that two ideal amorphous structures, type I and type II, based on maximally random jammed packing of spheres of equal size, form a distinct class of ideal amorphous solids. The ideal amorphous structures contain wide variations in local density, limited by the condition of solidity. Four distinct characteristics, based on statistical geometry and topology, are shown to define this class. Voronoi tessellations carried out on simulated cells of random packed spheres and amorphous polymers give a broad distribution of individual volumes, skewed, with a tail at the high volume end.     

Entropy of jammed matter

We investigate the nature of randomness in disordered packings of frictional spheres. We calculate the entropy of 3D packings through the force and volume ensemble of jammed matter, a mesoscopic ensemble and numerical simulations using volume fluctuation analysis and graph theoretical methods. Equations of state are obtained relating entropy, volume fraction and compactivity characterizing the different states of jammed matter. At the mesoscopic level the entropy vanishes at random close packing, while the microscopic states contribute to a finite entropy. The entropy of the jammed system reveals that the random loose packings are more disordered than random close packings, allowing for an unambiguous interpretation of both limits.     

A first-order phase transition defines the random close packing of hard spheres

Randomly packing spheres of equal size into a container consistently results in a static configuration with a density of ∼64%. The ubiquity of random close packing (RCP) rather than the optimal crystalline array at 74% begs the question of the physical law behind this empirically deduced state. Indeed, there is no signature of any macroscopic quantity with a discontinuity associated with the observed packing limit. Here we show that RCP can be interpreted as a manifestation of a thermodynamic singularity, which defines it as the “freezing point” in a first-order phase transition between ordered and disordered packing phases. Despite the athermal nature of granular matter, we show the thermodynamic character of the transition in that it is accompanied by sharp discontinuities in volume and entropy. This occurs at a critical compactivity, which is the intensive variable that plays the role of temperature in granular matter. Our results predict the experimental conditions necessary for the formation of a jammed crystal by calculating an analogue of the “entropy of fusion”. This approach is useful since it maps out-of-equilibrium problems in complex systems onto simpler established frameworks in statistical mechanics.