Compression-driven jamming of athermal frictionless spherocylinders in two dimensions.

Abstract

We simulate numerically the compression-driven jamming of athermal, frictionless, soft-core spherocylinders in two dimensions, for a range of particle aspect ratios α. We find the critical packing fraction ϕ_{J}(α) for the jamming transition and the average number of contacts per particle z_{J}(α) at jamming. We find that both are nonmonotonic, with a peak at α≈1. We find that configurations at the compression-driven jamming point are always hypostatic for all α, with z_{J}<z_{iso}=2d_{f}=6 the isostatic value. We show that, for moderately elongated spherocylinders, there is no orientational ordering upon athermal compression through jamming. We analyze in detail the eigenmodes of the dynamical matrix close to the jamming point for a few different values of the aspect ratio, from nearly circular to moderately elongated. We find that there are low frequency bands containing N(z_{iso}-z_{J})/2 modes, such that the frequencies of these modes vanish as ϕ→ϕ_{J}. We consider the extended versus localized nature of these low frequency modes, and the extent to which they involve translational or rotational motion, and find many low frequency sliding modes where particles can move with little rotation. We highlight the importance of treating side-to-side contacts, along flat sides of the spherocylinder, properly for the correct determination of z_{J}. We note the singular nature of taking the α→0 limit. We discuss the similarities and differences with previous work on jammed ellipses and ellipsoids, to illustrate the effects that different particle shapes have on configurations at jamming.