Constraints and vibrations in static packings of ellipsoidal particles

TL;DR

This study numerically explores the mechanical stability and vibrational properties of static packings of ellipsoidal particles, revealing unique zero-energy modes and stabilization mechanisms that differ from spherical particle packings.

Contribution

It introduces a detailed decomposition of the dynamical matrix to understand the stability and vibrational modes of ellipsoidal particle packings, highlighting new zero-energy modes and stabilization effects.

Findings

Presence of zero eigenvalue modes in the stiffness matrix at finite compression.

Finite compression stabilizes packings through nearly eigenvector modes.

Potential energy scales as δ^4 along quartic modes at jamming onset.

Abstract

We numerically investigate the mechanical properties of static packings of ellipsoidal particles in 2D and 3D over a range of aspect ratio and compression $\Delta \phi$. While amorphous packings of spherical particles at jamming onset ($\Delta \phi=0$) are isostatic and possess the minimum contact number $z_{\rm iso}$ required for them to be collectively jammed, amorphous packings of ellipsoidal particles generally possess fewer contacts than expected for collective jamming ($z < z_{\rm iso}$) from naive counting arguments, which assume that all contacts give rise to linearly independent constraints on interparticle separations. To understand this behavior, we decompose the dynamical matrix $M=H-S$ for static packings of ellipsoidal particles into two important components: the stiffness $H$ and stress $S$ matrices. We find that the stiffness matrix possesses $N(z_{\rm iso} - z)$ eigenmodes ${\hat e}_0$ with zero eigenvalues even at finite compression, where $N$ is the number of particles. In addition, these modes ${\hat e}_0$ are nearly eigenvectors of the dynamical matrix with eigenvalues that scale as $\Delta \phi$, and thus finite compression stabilizes packings of ellipsoidal particles. At jamming onset, the harmonic response of static packings of ellipsoidal particles vanishes, and the total potential energy scales as $\delta^4$ for perturbations by amplitude $\delta$ along these `quartic' modes, ${\hat e}_0$. These findings illustrate the significant differences between static packings of spherical and ellipsoidal particles.