Disordered Contact Networks in Jammed Packings of Frictionless Disks

Abstract

We analyse properties of contact networks formed in packings of soft frictionless disks near the unjamming transition. We construct polygonal tilings and triangulations of the contact network that partitions space into convex regions which are either covered or uncovered. This allows us to characterize the local spatial structure of the packing near the transition using well-defined geometric objects. We construct bounds on the number of polygons and triangulation vectors that appear in such packings. We study these networks using simulations of bidispersed disks interacting via a one-sided linear spring potential. We find that several underlying geometric distributions are reproducible and display self averaging properties. We find that the total covered area is a reliable real space parameter that can serve as a substitute for the packing fraction. We find that the unjamming transition occurs at a fraction of covered area \(A_G^{*} = 0.446(1)\). We determine scaling exponents of the excess covered area as the energy of the system approaches zero \(E_G \to 0^+\), and the coordination number \(\langle z_g \rangle\) approaches its isostatic value \(\Delta Z = \langle z_g \rangle - \langle z_g \rangle_{\rm iso} \to 0^{+}\). We find \(\Delta A_G \sim \Delta {E_G}^{0.28(1)}\) and \(\Delta A_G \sim \Delta Z^{1.00(1)}\), representing new structural critical exponents. We use the distribution functions of local areas to study the underlying geometric disorder in the packings. We find that a finite fraction of order \(\Psi_O^* = 0.369(1)\) persists as the transition is approached.