Dense Regular Packings of Irregular Non-Convex Particles Joost de Graaf,1 Ren´evan Roij,2 and Marjolein Dijkstra1 1Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands 2Institute for Theoretical Physics, Utrecht University, Leuvenlaan 4, 3584 CE Utrecht, The Netherlands (Dated: September 22, 2018) We present a new numerical scheme to study systems of non-convex, irregular, and punctured par- ticles in an efficient manner. We employ this method to analyze regular packings of odd-shaped bod- ies, not only from a nanoparticle but also both from a computational geometry perspective. Besides determining close-packed structures for many shapes, we also discover a new denser configuration for Truncated Tetrahedra. Moreover, we consider recently synthesized nanoparticles and colloids, where we focus on the excluded volume interactions, to show the applicability of our method in the investigation of their crystal structures and phase behavior. Extensions to the presented scheme include the incorporation of soft particle-particle interactions, the study of quasicrystalline systems, and random packings. PACS numbers: 61.66.-f, 61.46.Df, 82.70.Dd, 02.60.-x The synthesis of colloids and nanoparticles has ad- In this Letter, we present a novel numerical method vanced enormously over the last decade [1{6]. Currently to study systems of non-convex irregular polytopes, col- it is not only possible to synthesize spherical particles, loids, and nanoparticles. We employ this method to es- but also a wide variety of convex faceted shapes, such tablish a rigorous lower bound φlb to the volume fraction as tetrahedra, cubes and octahedra [1, 4]. Perhaps the of their densest regular packing and to predict candidate most remarkable advancement in synthesis techniques is crystal structures. Our numerical scheme elegantly and the capability to create with high precision and repro- systematically reduces simulation studies of these com- ducibility non-convex, irregular, and even punctured par- plex systems in a two-fold way. Firstly, the problem ticles, e.g., colloidal caps [5], tetrapods [2], octapods [6], of determining overlaps between irregular particles is re- and nanostars [3]. Along with the increased availabil- duced tremendously, especially when these particles are ity of complex shapes, there is the concurrent increase faceted [4], by approximating the particle shape with tri- in the study of their self-assembly into liquid [7], amor- angles [24], since triangle intersections can be straightfor- phous [8, 9], and ordered phases [10]. Especially the or- wardly determined. Secondly, predicting candidate crys- dered (quasi)crystalline structures [10{13] of odd-shaped tal structures is greatly simplified by using the Floppy particles, as well as their material properties, have re- Box Monte Carlo method [25] (FBMC), which efficiently ceived a lot of attention in the past years. Interestingly, analyzes the regular structure on a unit-cell level. We an- the research into dense packings of particles is not re- ticipate our investigation to be the starting point of an stricted to materials science, as it is also connected to exploration of (dis)ordered packings of arbitrarily shaped fields as diverse as discrete geometry and number the- particles, which rivals the current research effort on more ory [14{16], computer science and biophysics [17]. conventional bodies. Predictions obtained from computer simulations on To prove the accuracy of our method, we have repro- the phase behaviour and the self-assembled structures duced the lower bound φlb both for the Platonic and for of these particles have been essential in guiding experi- the Archimedean Solids, see Ref. [26]. We find excellent mental studies and in answering fundamental mathemat- agreement with Refs. [15, 19{21]. The deviation between ∗ −7 arXiv:1107.0603v2 [cond-mat.soft] 5 Jul 2011 ical questions on the packing of particles. Convex ob- φlb and the literature value φlb ranges from 2 · 10 to −3 ∗ jects such as spheres [7, 16] and ellipsoids [9], as well as 3·10 . For all results obtained by FBMC φlb < φlb, be- (semi)regular [13{15, 18{21] and space-filling [22] solids cause compression to the mathematically derived packing have been the subject of intense ongoing investigation. fraction [15, 21] is virtually impossible. However, FBMC However, ordered structures comprised of irregular non- performs remarkably well for these convex shapes: the convex particles have hardly been studied by simulation simulations we performed typically yielded a very nar- because of the numerical challenges in implementing ex- row distribution of crystal-structure candidates near the cluded volume interactions for such systems, due to the closest-packed configuration, the densest of which only complex particle shape and the additional rotational de- required minimal additional compression to achieve the grees of freedom. Only recently were the first attempts given value of φlb. Moreover, for Truncated Tetrahedra made to study such systems, namely for superdisks and we have discovered a new crystal structure, a dimer lat- superballs [23]. tice, with φlb = 0:988 ::: [26]. This is not only mathe- 2 matically interesting, it is also relevant to the study of PI DT DD TT TI ST nanoparticle systems, since Truncated Tetrahedral parti- cles have recently been synthesized [1]. 1 For all prisms and antiprisms with regular n-gonal bases (n > 3), for all Johnson Solids, and for all Cata- 0.8 lan Solids, we have verified Ulam's conjecture [27], which states that all convex objects pack denser than spheres. More importantly, for all of these sets we also find that 0.6 the obtained densest packing for centrally symmetric polyhedra is a Bravais lattice packing, in accordance φ with the conjecture of Ref. [19], and that for the con- 0.4 vex, congruent solids without central symmetry it is not φlb (FBMC) o a Bravais lattice, in accordance with the conjecture of φlb (OBBA) 0.2 o Ref. [20]. Here, we present only the data on the Cata- φlb (OSA) o lan Solids, which are dual to the Archimedean Solids. In φub (ISA) Fig. 1 several bounds to the densest packing fraction for φ 0 sph the Catalan Solids are shown. The outscribed-sphere ap- proximation (OSA) and oriented-bounding-box approx- PH PD DH DI RT TH RD imation (OBBA) yield lower bounds to the maximum packing fraction, the inscribed-sphere approximation [19] (ISA) provides an upper bound. Ref. [26] gives a more comprehensive account of the terminology and the em- ployed techniques. It also lists the values given in Fig. 1 PH PI PD DT DH with 6 digit precision. The ISA in combination with the FBMC result gives a much narrower region in which the densest packing can be found for this group of solids than was previously established. The narrowness of the region can be attributed to the relatively low asphericity of the DD DI TT RT Catalan Solids. Other particle shapes for which we have determined the densest-packing crystal structures include the regu- lar Penta- and Heptaprism: φlb = 0:921 ::: and φlb = TI TH ST RD 0:892 ::: , respectively [26]. The interest in Pentaprisms stems not only from the fact that such particles can nowa- FIG. 1. Upper and lower bounds to the densest pack- days be synthesized [28], but also from the fact that Pen- ing fraction for the 13 Catalan Solids. The graph shows the results for the outscribed-sphere approximation (OSA, tagons offer the exciting possibility of a structure with ◦ φlb, circles), oriented-bounding-box approximation (OBBA, quasicrystalline order in two dimensions [29]. Hence, ◦ φlb, squares), and FBMC method (φlb, connected crosses), we consider the prismatic equivalent of Pentagons, Pen- which give lower bounds to the maximum packing fraction, ◦ taprisms, in our three dimensional (3D) method. Hep- and the inscribed-sphere approximation [19] (ISA, φub, dia- taprisms are studied for the purpose of comparison. The monds), which gives an upper bound to this packing fraction. structures we obtain by FBMC for these Penta- and Hep- The value of the densest packing for spheres φsph is indicated taprisms are the 3D columnar continuation of the two di- by a red line. A visual representation of the Catalan Solids mensional (2D) double-lattices found for Pentagons and is given below the graph, the abbreviations are as follows: Pentagonal Hexecontrahedron (PH), Pentagonal Icositetra- Heptagons, respectively. The prisms in different columns, hedron (PI), Pentakis Dodecahedron (PD), Disdyakis Tria- which themselves form a 2D crystal, are out of register. contrahedron (DT), Deltoidal Hexecontrahedron (DH), Dis- The 2D results, to which we compare our 3D crystals, dyakis Dodecahedon (DD), Deltoidal Icositetrahedron (DI), were established by mathematical techniques [14] and Triakis Tetrahedron (TT) , Rhombic Triacontrahedron (RT), simulations [18]. Our results indicate that a 3D colum- Triakis Icosahedron (TI), Tetrakis Hexahedron (TH), Small nar continuation of a Pentaprism quasicrystalline con- Triakis Octahedron (ST), and Rhombic Dodecahedron (RD). figuration may also exist, since columnar continuations Note the improvement of FBMC with respect to the OSA and OBBA lower bound: the maximum packing fraction can appear to be common for non-space-filling prismatic sys- henceforth be located in the gray area above the FBMC line. tems. For space-filling Tri- and Hexaprisms (and Cubes) we observe both columnar continuations and randomly stacked layers, i.e., the prisms form crystalline planes, but the ordering in different planes is not in register. (irregular) shapes, which may even contain holes, thereby We have also applied our method to study non-convex going beyond existing studies. Fig. 2 shows representa- 3 a) b) c) d) e) f) g) h) FIG. 2. Four non-convex particle species and their densest-packed crystals. An Echidnahedron (a), a Great Stellated Dodecahedron (b), Jessen's Orthogonal Icosahedron (c), and a Stanford Bunny model (d). The crystal structures found by FBMC are given in (e) - (h) respectively, which show the unit cell and 7 of its periodic images.