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 an γ=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.
When is the ball a local pessimum for covering?
Y. Kallus (show abstract)
We consider the problem of identifying the worst point-symmetric shape for covering n-dimensional
Euclidean space with lattice translates. Here we focus on the dimensions where
the thinnest lattice covering with balls is known and ask whether the ball is
a pessimum for covering in these dimensions compared to all point-symmetric convex shapes.
We find that the ball is a local pessimum in 3 dimensions,
but not so for 4 and 5 dimensions.
Discrete and Computational Geometry 54, 232–245 (2015)
1301.5895
10.1007/s00454-015-9690-1
Pessimal packing shapes
Y. Kallus (show abstract)
We address the question of which convex shapes, when packed as densely as possible under certain restrictions, fill the least space and leave the most empty space. In each different dimension and under each different set of restrictions, this question is expected to have a different answer or perhaps no answer at all. As the problem of identifying global minima in most cases appears to be beyond current reach, in this paper we focus on local minima. We review some known results and prove these new results: in two dimensions, the regular heptagon is a local minimum of the double-lattice packing density, and in three dimensions, the directional derivative (in the sense of Minkowski addition) of the double-lattice packing density at the point in the space of shapes corresponding to the ball is in every direction positive.
Geometry & Topology 19, 343 (2015)
1305.0289
10.2140/gt.2015.19.343
Marginal stability in jammed packings: quasicontacts and weak contacts
Y. Kallus and S. Torquato (show abstract)
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θ.
Physical Review E 90, 022114 (2014)
1405.2087
10.1103/PhysRevE.90.022114
The 3-ball is a local pessimum for packing
Y. Kallus (show abstract)
It was conjectured by Ulam that the ball has the lowest optimal packing fraction out of all convex, three-dimensional solids. Here we
prove that any origin-symmetric convex solid of sufficiently small asphericity can be packed at a higher efficiency than balls. We also
show that in dimensions 4, 5, 6, 7, 8, and 24 there are origin-symmetric convex bodies of arbitrarily small asphericity that cannot be
packed using a lattice as efficiently as balls can be.
Advances in Mathematics 264:355 (2014)
1212.2551
10.1016/j.aim.2014.07.015
Jammed lattice sphere packings
Y. Kallus, É. Marcotte, and S. Torquato (show abstract)
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.
Physical Review E 88, 062151 (2013)
1310.5971
10.1103/PhysRevE.88.062151
Inextensible domains
Y. Kallus (show abstract)
We develop a theory of planar, origin-symmetric, convex domains that are inextensible with respect to lattice covering, that is, domains such that augmenting them in any way allows fewer domains to cover the same area. We show that origin-symmetric inextensible domains are exactly the origin-symmetric convex domains with a circle of outer billiard triangles. We address a conjecture by Genin and Tabachnikov about convex domains, not necessarily symmetric, with a circle of outer billiard triangles, and show that it follows immediately from a result of Sas.
Geometriae Dedicata 173, 177 (2014)
1301.5880
10.1007/s10711-013-9935-1
Statistical mechanics of the lattice sphere packing problem
Y. Kallus (show abstract)
We present an efficient Monte Carlo method for the lattice sphere packing problem in d dimensions. We use this method to numerically discover de novo the densest lattice sphere packing in dimensions 9 through 20. Our method goes beyond previous methods, not only in exploring higher dimensions but also in shedding light on the statistical mechanics underlying the problem in question. We observe evidence of a phase transition in the thermodynamic limit d➝∞. In the dimensions explored in the present work, the results are consistent with a first-order crystallization transition but leave open the possibility that a glass transition is manifested in higher dimensions.
Physical Review E 87 063307 (2013).
1305.1310
10.1103/PhysRevE.87.063307
Dense-packing crystal structures of physical tetrahedra
Y. Kallus and V. Elser (show abstract)
We present a method for discovering dense packings of general convex hard particles and apply it to study the dense packing behavior of a one-parameter family of particles with tetrahedral symmetry representing a deformation of the ideal mathematical tetrahedron into a less ideal, physical, tetrahedron and all the way to the sphere. Thus, we also connect the two well-studied problems of sphere packing and tetrahedron packing on a single axis. Our numerical results uncover a rich optimal-packing behavior, compared to that of other continuous families of particles previously studied. We present four structures as candidates for the optimal packing at different values of asphericity, providing an atlas of crystal structures that might be observed in systems of nanoparticles with tetrahedral symmetry.
Physical Review E 83, 036703 (2011)
1011.4034
10.1103/PhysRevE.83.036703
Upper bound on the packing density of regular tetrahedra and octahedra
S. Gravel, V. Elser, and Y. Kallus (show abstract)
Aristotle contended that (regular) tetrahedra tile space, an opinion that remained widespread until it was observed that non-overlapping
tetrahedra cannot subtend a solid angle of 4π around a point if this point lies on a tetrahedron edge. From this 15th century argument,
we can deduce that tetrahedra do not tile space but, more than 500 years later, we are unaware of any known non-trivial upper bound to the
packing density of tetrahedra. In this article, we calculate such a bound. To this end, we show the existence, in any packing of regular
tetrahedra, of a set of disjoint spheres centered on tetrahedron edges, so that each sphere is not fully covered by the packing. The bound
on the amount of space that is not covered in each sphere is obtained in a recursive way by building on the solid angle argument. The argument
can be readily modified to apply to other polyhedra. The resulting lower bound on the fraction of empty space in a packing of regular tetrahedra
is 2.6…×10−25 and reaches 1.4…×10−12 for regular octahedra.
Discrete and Computational Geometry 46, 799 (2011)
1008.2830
10.1007/s00454-010-9304-x
Method for dense packing discovery
Y. Kallus, V. Elser, and S. Gravel (show abstract)
The problem of packing a system of particles as densely as possible is foundational in the field of discrete geometry and is a powerful model in the material and biological sciences. As packing problems retreat from the reach of solution by analytic constructions, the importance of an efficient numerical method for conducting de novo (from-scratch) searches for dense packings becomes crucial. In this paper, we use the divide and concur framework to develop a general search method for the solution of periodic constraint problems, and we apply it to the discovery of dense periodic packings. An important feature of the method is the integration of the unit-cell parameters with the other packing variables in the definition of the configuration space. The method we present led to previously reported improvements in the densest-known tetrahedron packing. Here, we use the method to reproduce the densest-known lattice sphere packings and the best-known lattice kissing arrangements in up to 14 and 11 dimensions, respectively, providing numerical evidence for their optimality. For nonspherical particles, we report a dense packing of regular four-dimensional simplices with density φ=128/219≈0.5845 and with a similar structure to the densest-known tetrahedron packing.
Physical Review E 82, 056707 (2010)
1003.3301
10.1103/PhysRevE.82.056707
Dense periodic packings of tetrahedra with small repeating units
Y. Kallus, V. Elser, and S. Gravel (show abstract)
We present a one-parameter family of periodic packings of regular tetrahedra, with the packing fraction 100/117≈0.8547, that are simple in the sense
that they are transitive and their repeating units involve only four tetrahedra. The construction of the packings was inspired from results of a numerical
search that yielded a similar packing. We present an analytic construction of the packings and a description of their properties. We also present a transitive
packing with a repeating unit of two tetrahedra and a packing fraction (139+40√10)/369≈0.7194.
Discrete and Computational Geometry 44, 245 (2010)
0910.5226
10.1007/s00454-010-9254-3