(a) p = 0 . 85 (b) p = 1 . 5 FIG. 5: (color online). Typical configurations of MRJ packings of superballs with different values of the deformation parameter p . 9 superdisks. We do not use monodispersed superdisk systems here because they are easily crystallized into ordered packings [25]. For p = 1, one obtains the binary circular-disk sys- tem which has been intensively studied as a prototypical glass former [41]. Typical jammed packing configurations are shown in Fig. 3. The density φ and the average contact number per particle Z as a function of p are shown in Fig. 4, which reveals that the initial rapid in- creases of φ and Z are linear in | p − 1 | [42]. The density φ increases monotonically as p moves away from unity and shows no signs of a plateau, even for relatively large p . In addition, φ quickly surpasses the density of the optimal binary circular-disk packing associated with the size and molar ratios employed here, which contains phase-separate regions of triangular lattice packings of different sized circular disks [41]. The quantity Z quickly reaches its plateau value Z ∗ ≈ 4 . 7 at p ≈ 1 . 3, which is smaller than 2 f = 6, indicating the packings are hypostatic. The cubatic order parameter P$_{4}$ is defined as P$_{4}$ = 〈 (35 cos 4 θ − 30 cos 2 θ +3) / 8 〉 , where θ is the angle between the particle axis and the director, along which the principle axes of the particles have maximum mutual alignment [37]. The measured cubatic order parameter is P$_{4}$ ≈ 0 . 06 to 0 . 32 with the tendency to increase as | p − 1 | grows.