Hypostatic jammed packings of frictionless nonspherical particles

We perform computational studies of static packings of a variety of nonspherical particles including circulo-lines, circulo-polygons, ellipses, asymmetric dimers, and dumbbells to determine which shapes form hypostatic versus isostatic packings and to understand why hypostatic packings of nonspherical particles can be mechanically stable despite having fewer contacts than that predicted from na\"ive constraint counting. To generate highly accurate force- and torque-balanced packings of circulo-lines and -polygons, we developed an interparticle potential that gives continuous forces and torques as a function of the particle coordinates. We show that the packing fraction and coordination number at jamming onset obey a master-like form for all of the nonspherical particle packings we studied when plotted versus the particle asphericity ${\cal A}$, which is proportional to the ratio of the squared perimeter to the area of the particle. Further, the eigenvalue spectra of the dynamical matrix for packings of different particle shapes collapse when plotted at the same ${\cal A}$. For hypostatic packings of nonspherical particles, we verify that the number of "quartic" modes along which the potential energy increases as the fourth power of the perturbation amplitude matches the number of missing contacts relative to the isostatic value. In addition, we calculate the principal curvatures of the inequality constraints for each contact in circulo-line packings and identify specific types of contacts with inequality constraints that possess convex curvature. These contacts can constrain multiple degrees of freedom and allow hypostatic packings of nonspherical particles to be mechanically stable.
Comment: 17 pages, 22 figures