Generating ultradense jammed ellipse packings using biased SWAP

Abstract

Using a Lubachevsky-Stillinger-like growth algorithm combined with biased SWAP Monte Carlo and transient degrees of freedom, we generate ultradense disordered jammed ellipse packings. For all aspect ratios $\alpha$, these packings exhibit significantly smaller intermediate-wavelength density fluctuations and greater local nematic order than their less-dense counterparts. The densest packings are disordered despite having packing fractions ${\varphi }_{\mathrm{J}}(\alpha )$ that are within less than 0.5% of that of the monodisperse-ellipse crystal [${\varphi }_{\mathrm{xtal}}=\pi /(2\sqrt{3})\simeq .9069$] over the range $1.25\lesssim \alpha \lesssim 1.4$ and coordination numbers $Z_{\mathrm{J}}(\alpha )$ that are within less than 0.5% of isostaticity [$Z_{\mathrm{iso}}=6$] over the range $1.3\lesssim \alpha \lesssim 2.0$. Lower-$\alpha$ packings are strongly fractionated and consist of polycrystals of intermediate-size particles, with the largest and smallest particles isolated at the grain boundaries. Higher-$\alpha$ packings are also fractionated, but in a qualitatively-different fashion; they are composed of increasingly-large locally-nematic domains reminiscent of liquid glasses.