Robust Algorithm to Generate a Diverse Class of Dense Disordered and Ordered Sphere Packings via Linear Programming
arXiv:1008.2747 · doi:10.1103/PhysRevE.82.061302
Abstract
We have formulated the problem of generating periodic dense paritcle packings as an optimization problem called the Adaptive Shrinking Cell (ASC) formulation [S. Torquato and Y. Jiao, Phys. Rev. E {\bf 80}, 041104 (2009)]. Because the objective function and impenetrability constraints can be exactly linearized for sphere packings with a size distribution in -dimensional Euclidean space , it is most suitable and natural to solve the corresponding ASC optimization problem using sequential linear programming (SLP) techniques. We implement an SLP solution to produce robustly a wide spectrum of jammed sphere packings in for and with a diversity of disorder and densities up to the maximally densities. This deterministic algorithm can produce a broad range of inherent structures besides the usual disordered ones with very small computational cost by tuning the radius of the {\it influence sphere}. In three dimensions, we show that it can produce with high probability a variety of strictly jammed packings with a packing density anywhere in the wide range . We also apply the algorithm to generate various disordered packings as well as the maximally dense packings for and 6. Compared to the LS procedure, our SLP protocol is able to ensure that the final packings are truly jammed, produces disordered jammed packings with anomalously low densities, and is appreciably more robust and computationally faster at generating maximally dense packings, especially as the space dimension increases.
34 pages, 6 figures
References in corpus (7)
- Jamming at Zero Temperature and Zero Applied Stress: the Epitome of Disorder
- Packing Hyperspheres in High-Dimensional Euclidean Spaces
- Hypoconstrained Jammed Packings of Nonspherical Hard Particles: Ellipses and Ellipsoids
- Dense Packings of Superdisks and the Role of Symmetry
- Exactly Solvable Disordered Sphere-Packing Model in Arbitrary-Dimension Euclidean Spaces
- Toward the Jamming Threshold of Sphere Packings: Tunneled Crystals
- Geometrical Frustration: A Study of 4d Hard Spheres