An ellipsoid, the simplest nonspherical shape, has been extensively used as a model for elongated building blocks for a wide spectrum of molecular, colloidal, and granular systems. Yet the densest packing of congruent hard ellipsoids, which is intimately related to the high-density phase of many condensed matter systems, is still an open problem. We discover an unusual family of dense crystalline packings of self-dual ellipsoids (ratios of the semiaxes α:√α:1), containing 24 particles with a quasi-square-triangular (SQ-TR) tiling arrangement in the fundamental cell. The associated packing density ϕ exceeds that of the densest known SM2 crystal [ A. Donev et al., Phys. Rev. Lett. 92, 255506 (2004)] for aspect ratios α in (1.365, 1.5625), attaining a maximal ϕ≈0.75806... at α=93/64. We show that the SQ-TR phase derived from these dense packings is thermodynamically stable at high densities over the aforementioned α range and report a phase diagram for self-dual ellipsoids. The discovery of the SQ-TR crystal suggests organizing principles for nonspherical particles and self-assembly of colloidal systems.

The critical behaviors of a granular system at the jamming transition have been extensively studied from both mechanical and thermodynamic perspectives. In this work, we numerically investigate the jamming behaviors of a variety of frictionless non-spherical particles, including spherocylinder, ellipsoid, spherotetrahedron and spherocube. In particular, for a given particle shape, a series of random configurations at different fixed densities are generated and relaxed to minimize interparticle overlaps using the relaxation algorithm. We find that as the jamming point (i.e., point J">J) is approached, the number of iteration steps (defined as the “time-scale” for our systems) required to completely relax the interparticle overlaps exhibits a clear power-law divergence. The dependence of the detailed mathematical form of the power-law divergence on particle shapes is systematically investigated and elucidated, which suggests that the shape effects can be generally categorized as elongation and roundness. Importantly, we show the jamming transition density can be accurately determined from the analysis of time-scale divergence for different non-spherical shapes, and the obtained values agree very well with corresponding ones reported in literature. Moreover, we study the plastic behaviors of over-jammed packings of different particles under a compression–expansion procedure and find that the jamming of ellipsoid is much more robust than other non-spherical particles. This work offers an alternative approximate procedure besides conventional packing algorithms for studying athermal jamming transition in granular system of frictionless non-spherical particles.