Exact constructions of a family of dense periodic packings of tetrahedra

Kinetic Frustration Effects on Dense Two-Dimensional Packings of Convex Particles and Their Structural Characteristics

The study of hard-particle packings is of fundamental importance in physics, chemistry, cell biology, and discrete geometry. Much of the previous work on hard-particle packings concerns their densest possible arrangements. By contrast, we examine kinetic effects inevitably present in both numerical and experimental packing protocols. Specifically, we determine how changing the compression/shear rate of a two-dimensional packing of noncircular particles causes it to deviate from its densest possible configuration, which is always periodic. The adaptive shrinking cell (ASC) optimization scheme maximizes the packing fraction of a hard-particle packing by first applying random translations and rotations to the particles and then isotropically compressing and shearing the simulation box repeatedly until a possibly jammed state is reached. We use a stochastic implementation of the ASC optimization scheme to mimic different effective time scales by varying the number of particle moves between compressions/shears. We generate dense, effectively jammed, monodisperse, two-dimensional packings of obtuse scalene triangle, rhombus, curved triangle, lens, and "ice cream cone" (a semicircle grafted onto an isosceles triangle) shaped particles, with a wide range of packing fractions and degrees of order. To quantify these kinetic effects, we introduce the kinetic frustration index K, which measures the deviation of a packing from its maximum possible packing fraction. To investigate how kinetics affect short- and long-range ordering in these packings, we compute their spectral densities χ̃_V(k) and characterize their contact networks. We find that kinetic effects are most significant when the particles have greater asphericity, less curvature, and less rotational symmetry. This work may be relevant to the design of laboratory packing protocols.

Hard convex lens-shaped particles: Characterization of dense disordered packings

Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is "optimal" in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction ϕ_{MRJ}≃0.73 [G. Cinacchi and S. Torquato, Soft Matter 14, 8205 (2018)1744-683X10.1039/C8SM01519H]. This value is only a few percent lower than ϕ_{DKP}=0.76210⋯, the packing fraction of the corresponding densest-known crystalline (degenerate) packings [G. Cinacchi and S. Torquato, J. Chem. Phys. 143, 224506 (2015)JCPSA60021-960610.1063/1.4936938]. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method-based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows us to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated here can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of ten contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve.

Engineering entropy in soft matter: the bad, the ugly and the good

The role of entropic interactions, often subtle and sometimes crucial, on the structure and properties of soft matter has a well-recognized place in the classic and modern scientific literature. However, the lessons learned from many of those studies do not always form part of the standard arsenal of strategies that are taught or used for de novo studies relevant to the engineering of new materials. Fortunately, a growing number of examples exist where entropic effects have been designed a priori to achieve a desired or new outcome. This tutorial review describes some recent such examples, selected to illustrate the potential benefits of a more pro-active approach to harnessing the often overlooked power of entropy.

Packing density of rigid aggregates is independent of scale

Large planetary seedlings, comets, microscale pharmaceuticals, and nanoscale soot particles are made from rigid, aggregated subunits that are compacted under low compression into larger structures spanning over 10 orders of magnitude in dimensional space. Here, we demonstrate that the packing density (θf) of compacted rigid aggregates is independent of spatial scale for systems under weak compaction. The θf of rigid aggregated structures across six orders of magnitude were measured using nanoscale spherical soot aerosol composed of aggregates with ∼ 17-nm monomeric subunits and aggregates made from uniform monomeric 6-mm spherical subunits at the macroscale. We find θf = 0.36 ± 0.02 at both dimensions. These values are remarkably similar to θf observed for comet nuclei and measured values of other rigid aggregated systems across a wide variety of spatial and formative conditions. We present a packing model that incorporates the aggregate morphology and show that θf is independent of both monomer and aggregate size. These observations suggest that the θf of rigid aggregates subject to weak compaction forces is independent of spatial dimension across varied formative conditions.

Disordered strictly jammed binary sphere packings attain an anomalously large range of densities

Previous attempts to simulate disordered binary sphere packings have been limited in producing mechanically stable, isostatic packings across a broad spectrum of packing fractions. Here we report that disordered strictly jammed binary packings (packings that remain mechanically stable under general shear deformations and compressions) can be produced with an anomalously large range of average packing fractions 0.634≤φ≤0.829 for small to large sphere radius ratios α restricted to α≥0.100. Surprisingly, this range of average packing fractions is obtained for packings containing a subset of spheres (called the backbone) that are exactly strictly jammed, exactly isostatic, and also generated from random initial conditions. Additionally, the average packing fractions of these packings at certain α and small sphere relative number concentrations x approach those of the corresponding densest known ordered packings. These findings suggest for entropic reasons that these high-density disordered packings should be good glass formers and that they may be easy to prepare experimentally. We also identify an unusual feature of the packing fraction of jammed backbones (packings with rattlers excluded). The backbone packing fraction is about 0.624 over the majority of the α-x plane, even when large numbers of small spheres are present in the backbone. Over the (relatively small) area of the α-x plane where the backbone is not roughly constant, we find that backbone packing fractions range from about 0.606 to 0.829, with the volume of rattler spheres comprising between 1.6% and 26.9% of total sphere volume. To generate isostatic strictly jammed packings, we use an implementation of the Torquato-Jiao sequential linear programming algorithm [Phys. Rev. E 82, 061302 (2010)], which is an efficient producer of inherent structures (mechanically stable configurations at the local maxima in the density landscape). The identification and explicit construction of binary packings with such high packing fractions could have important practical implications for granular composites where density is critical both to material properties and fabrication cost, including for solid propellants, concrete, and ceramics. The densities and structures of jammed binary packings at various α and x are also relevant to the formation of a glass phase in multicomponent metallic systems.

Jammed frictional tetrahedra are hyperstatic

We prepare packings of frictional tetrahedra with volume fractions ϕ ranging from 0.469 to 0.622 using three different experimental protocols under isobaric conditions. Analysis via x-ray microtomography reveals that the contact number Z grows with ϕ, but does depend on the preparation protocol. While there exist four different types of contacts in tetrahedra packings, our analysis shows that the edge-to-face contacts contribute about 50% of the total increase in Z. The number of constraints per particle C increases also with ϕ and even the loosest packings are strongly hyperstatic, i.e., mechanically overdetermined with C approximately twice the degrees of freedom each particle possesses.

Shape effects on the random-packing density of tetrahedral particles

Regular tetrahedra have been demonstrated recently giving high packing density in random configurations. However, it is unknown whether the random-packing density of tetrahedral particles with other shapes can reach an even higher value. A numerical investigation on the random packing of regular and irregular tetrahedral particles is carried out. Shape effects of rounded corner, eccentricity, and height on the packing density of tetrahedral particles are studied. Results show that altering the shape of tetrahedral particles by rounding corners and edges, by altering the height of one vertex, or by lateral displacement of one vertex above its opposite face, all individually have the effect of reducing the random-packing density. In general, the random-packing densities of irregular tetrahedral particles are lower than that of regular tetrahedra. The ideal regular tetrahedron should be the shape which has the highest random-packing density in the family of tetrahedra, or even among convex bodies. An empirical formula is proposed to describe the rounded corner effect on the packing density, and well explains the density deviation of tetrahedral particles with different roundness ratios. The particles in the simulations are verified to be randomly packed by studying the pair correlation functions, which are consistent with previous results. The spherotetrahedral particle model with the relaxation algorithm is effectively applied in the simulations.

Organizing principles for dense packings of nonspherical hard particles: not all shapes are created equal

We have recently devised organizing principles to obtain maximally dense packings of the Platonic and Archimedean solids and certain smoothly shaped convex nonspherical particles [Torquato and Jiao, Phys. Rev. E 81, 041310 (2010)]. Here we generalize them in order to guide one to ascertain the densest packings of other convex nonspherical particles as well as concave shapes. Our generalized organizing principles are explicitly stated as four distinct propositions. All of our organizing principles are applied to and tested against the most comprehensive set of both convex and concave particle shapes examined to date, including Catalan solids, prisms, antiprisms, cylinders, dimers of spheres, and various concave polyhedra. We demonstrate that all of the densest known packings associated with this wide spectrum of nonspherical particles are consistent with our propositions. Among other applications, our general organizing principles enable us to construct analytically the densest known packings of certain convex nonspherical particles, including spherocylinders, "lens-shaped" particles, square pyramids, and rhombic pyramids. Moreover, we show how to apply these principles to infer the high-density equilibrium crystalline phases of hard convex and concave particles. We also discuss the unique packing attributes of maximally random jammed packings of nonspherical particles.

Densest binary sphere packings

The densest binary sphere packings in the α-x plane of small to large sphere radius ratio α and small sphere relative concentration x have historically been very difficult to determine. Previous research had led to the prediction that these packings were composed of a few known "alloy" phases including, for example, the AlB(2) (hexagonal ω), HgBr(2), and AuTe(2) structures, and to XY(n) structures composed of close-packed large spheres with small spheres (in a number ratio of n to 1) in the interstices, e.g., the NaCl packing for n=1. However, utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [Torquato and Jiao, Phys. Rev. E 82, 061302 (2010)], we have discovered that many more structures appear in the densest packings. For example, while all previously known densest structures were composed of spheres in small to large number ratios of one to one, two to one, and very recently three to one, we have identified densest structures with number ratios of seven to three and five to two. In a recent work [Hopkins et al., Phys. Rev. Lett. 107, 125501 (2011)], we summarized these findings. In this work, we present the structures of the densest-known packings and provide details about their characteristics. Our findings demonstrate that a broad array of different densest mechanically stable structures consisting of only two types of components can form without any consideration of attractive or anisotropic interactions. In addition, the structures that we have identified may correspond to currently unidentified stable phases of certain binary atomic and molecular systems, particularly at high temperatures and pressures.

Self-assembly of uniform polyhedral silver nanocrystals into densest packings and exotic superlattices

Understanding how polyhedra pack into extended arrangements is integral to the design and discovery of crystalline materials at all length scales. Much progress has been made in enumerating and characterizing the packing of polyhedral shapes, and the self-assembly of polyhedral nanocrystals into ordered superstructures. However, directing the self-assembly of polyhedral nanocrystals into densest packings requires precise control of particle shape, polydispersity, interactions and driving forces. Here we show with experiment and computer simulation that a range of nanoscale Ag polyhedra can self-assemble into their conjectured densest packings. When passivated with adsorbing polymer, the polyhedra behave as quasi-hard particles and assemble into millimetre-sized three-dimensional supercrystals by sedimentation. We also show, by inducing depletion attraction through excess polymer in solution, that octahedra form an exotic superstructure with complex helical motifs rather than the densest Minkowski lattice. Such large-scale Ag supercrystals may facilitate the design of scalable three-dimensional plasmonic metamaterials for sensing, nanophotonics and photocatalysis.

Dense regular packings of irregular nonconvex particles

We present a new numerical scheme to study systems of nonconvex, irregular, and punctured particles in an efficient manner. We employ this method to analyze regular packings of odd-shaped bodies, both from a nanoparticle and from a computational geometry perspective. Besides determining close-packed structures for 17 irregular shapes, we confirm several conjectures for the packings of a large set of 142 convex polyhedra and extend upon these. We also prove that we have obtained the densest packing for both rhombicuboctahedra and rhombic enneacontrahedra and we have improved upon the packing of enneagons and truncated tetrahedra.

Isostaticity of constraints in amorphous jammed systems of soft frictionless Platonic solids

The average number of constraints per particle <C(total)> in mechanically stable amorphous systems of Platonic solids approaches the isostatic limit at the jamming point (<C(total)>→12), though average number of contacts are hypostatic. By introducing angular alignment metrics to classify the degree of constraint imposed by each contact, constraints are shown to arise as a direct result of local orientational order reflected in edge-face and face-face alignment angle distributions. With approximately one face-face contact per particle at jamming, chainlike face-face clusters form with finite extent--a signature of amorphous jammed systems.

Experimental evidence of icosahedral and decahedral packing in one-dimensional nanostructures

The packing of spheres is a subject that has drawn the attention of mathematicians and philosophers for centuries and that currently attracts the interest of the scientific community in several fields. At the nanoscale, the packing of atoms affects the chemical and structural properties of the material and, hence, its potential applications. This report describes the experimental formation of 5-fold nanostructures by the packing of interpenetrated icosahedral and decahedral units. These nanowires, formed by the reaction of a mixture of metal salts (Au and Ag) in the presence of oleylamine, are obtained when the chemical composition is specifically Ag/Au = 3:1. The experimental images of the icosahedral nanowires have a high likelihood with simulated electron micrographs of structures formed by two or three Boerdijk-Coxeter-Bernal helices roped on a single structure, whereas for the decahedral wires, simulations using a model of adjacent decahedra match the experimental structures. To our knowledge, this is the first report of the synthesis of nanowires formed by the packing of structures with 5-fold symmetry. These icosahedral nanowire structures are similar to those of quasicrystals, which can only be formed if at least two atomic species are present and in which icosahedral and decahedral packing has been found for bulk crystals.

Hyperuniformity, quasi-long-range correlations, and void-space constraints in maximally random jammed particle packings. II. Anisotropy in particle shape

We extend the results from the first part of this series of two papers by examining hyperuniformity in heterogeneous media composed of impenetrable anisotropic inclusions. Specifically, we consider maximally random jammed (MRJ) packings of hard ellipses and superdisks and show that these systems both possess vanishing infinite-wavelength local-volume-fraction fluctuations and quasi-long-range pair correlations scaling as r(-(d+1)) in d Euclidean dimensions. Our results suggest a strong generalization of a conjecture by Torquato and Stillinger [Phys. Rev. E 68, 041113 (2003)], namely, that all strictly jammed saturated packings of hard particles, including those with size and shape distributions, are hyperuniform with signature quasi-long-range correlations. We show that our arguments concerning the constrained distribution of the void space in MRJ packings directly extend to hard-ellipse and superdisk packings, thereby providing a direct structural explanation for the appearance of hyperuniformity and quasi-long-range correlations in these systems. Additionally, we examine general heterogeneous media with anisotropic inclusions and show unexpectedly that one can decorate a periodic point pattern to obtain a hard-particle system that is not hyperuniform with respect to local-volume-fraction fluctuations. This apparent discrepancy can also be rationalized by appealing to the irregular distribution of the void space arising from the anisotropic shapes of the particles. Our work suggests the intriguing possibility that the MRJ states of hard particles share certain universal features independent of the local properties of the packings, including the packing fraction and average contact number per particle.

Mesophase behaviour of polyhedral particles

Translational and orientational excluded-volume fields encoded in particles with anisotropic shapes can lead to purely entropy-driven assembly of morphologies with specific order and symmetry. To elucidate this complex correlation, we carried out detailed Monte Carlo simulations of six convex space-filling polyhedrons, namely, truncated octahedrons, rhombic dodecahedrons, hexagonal prisms, cubes, gyrobifastigiums and triangular prisms. Simulations predict the formation of various new liquid-crystalline and plastic-crystalline phases at intermediate volume fractions. By correlating these findings with particle anisotropy and rotational symmetry, simple guidelines for predicting phase behaviour of polyhedral particles are proposed: high rotational symmetry is in general conducive to mesophase formation, with low anisotropy favouring plastic-solid behaviour and intermediate anisotropy (or high uniaxial anisotropy) favouring liquid-crystalline behaviour. It is also found that dynamical disorder is crucial in defining mesophase behaviour, and that the apparent kinetic barrier for the liquid-mesophase transition is much lower for liquid crystals (orientational order) than for plastic solids (translational order).

Maximum and minimum stable random packings of Platonic solids

Motivated by the relation between particle shape and packing, we measure the volume fraction ϕ occupied by the Platonic solids which are a class of polyhedrons with congruent sides, vertices, and dihedral angles. Tetrahedron-, cube-, octahedron-, dodecahedron-, and icosahedron-shaped plastic dice were fluidized or mechanically vibrated to find stable random loose packing ϕ(rlp)=0.51,0.54,0.52,0.51,0.50 and densest packing ϕ(rcp)=0.64,0.67,0.64,0.63,0.59, respectively, with standard deviation of ≃ ± 0.01. We find that ϕ obtained by all protocols peak at the cube, which is the only Platonic solid that can tessellate space, and then monotonically decrease with number of sides. This overall trend is similar but systematically lower than the maximum ϕ reported for frictionless Platonic solids and below ϕ(rlp) of spheres for the loose packings. Experiments with ceramic tetrahedron were also conducted, and higher friction was observed to lead to lower ϕ.

Athermal jamming of soft frictionless Platonic solids

A mechanically based structural optimization method is utilized to explore the phenomena of jamming for assemblies of frictionless Platonic solids. Systems of these regular convex polyhedra exhibit mechanically stable phases with density substantially less than optimal for a given shape, revealing that thermal motion is necessary to access high-density phases. We confirm that the large system jamming threshold of 0.623 ± 0.003 for tetrahedra is consistent with experiments on tetrahedral dice. Also, the extremely short-ranged translational correlations of packed tetrahedra observed in experiments are confirmed here, in contrast with those of thermally simulated glasses. Although highly ordered phases are observed to form for small numbers of cubes and dodecahedra, the short correlation length scale suppresses ordering in large systems, resulting in packings that are mechanically consistent with "orientationally disordered" contacts (point-face and edge-edge contacts). Mild nematic ordering is observed for large systems of cubes, whereas angular correlations for the remaining shapes are ultrashort ranged. In particular the angular correlation function of tetrahedra agrees with that recently observed experimentally for tetrahedral dice. Power-law scaling exponents for energy with respect to distance from the jamming threshold exhibit a clear dependence on the "highest-order" percolating contact topology. These nominal exponents are 6, 4, and 2 for configurations having percolating point-face (or edge-edge), edge-face, and face-face contacts, respectively. Jamming contact number is approximated for small systems of tetrahedra, icosahedra, dodecahedra, and octahedra with order and packing representative of larger systems. These Platonic solids exhibit hypostatic behavior, with average jamming contact number between the isostatic value for spheres and that of asymmetric particles. These shapes violate the isostatic conjecture, displaying contact number that decreases monotonically with sphericity. The common symmetry of dual polyhedra results in local translational structural similarity. Systems of highly spherical particles possessing icosahedral symmetry, such as icosahedra or dodecahedra, exhibit structural behavior similar to spheres, including jamming contact number and radial distribution function. These results suggest that although continuous rotational symmetry is broken by icosahedra and dodecahedra, the structural features of disordered packings of these particles are well replicated by spheres. Octahedra and cubes, which possess octahedral symmetry, exhibit similar local translational ordering, despite exhibiting strong differences in nematic ordering. In general, the structural features of systems with tetrahedra, octahedra, and cubes differ significantly from those of sphere packings.

Method for dense packing discovery

The problem of packing a system of particles as densely as possible is foundational in the field of discrete geometry and is a powerful model in the material and biological sciences. As packing problems retreat from the reach of solution by analytic constructions, the importance of an efficient numerical method for conducting de novo (from-scratch) searches for dense packings becomes crucial. In this paper, we use the divide and concur framework to develop a general search method for the solution of periodic constraint problems, and we apply it to the discovery of dense periodic packings. An important feature of the method is the integration of the unit-cell parameters with the other packing variables in the definition of the configuration space. The method we present led to previously reported improvements in the densest-known tetrahedron packing. Here, we use the method to reproduce the densest-known lattice sphere packings and the best-known lattice kissing arrangements in up to 14 and 11 dimensions, respectively, providing numerical evidence for their optimality. For nonspherical particles, we report a dense packing of regular four-dimensional simplices with density ϕ=128/219≈0.5845 and with a similar structure to the densest-known tetrahedron packing.

Phase behavior of colloidal superballs: shape interpolation from spheres to cubes

The phase behavior of hard superballs is examined using molecular dynamics within a deformable periodic simulation box. A superball's interior is defined by the inequality |x|(2q)+|y|(2q)+|z|(2q)≤1 , which provides a versatile family of convex particles (q≥0.5) with cubelike and octahedronlike shapes as well as concave particles (q<0.5) with octahedronlike shapes. Here, we consider the convex case with a deformation parameter q between the sphere point (q=1) and the cube (q=∞). We find that the asphericity plays a significant role in the extent of cubatic ordering of both the liquid and crystal phases. Calculation of the first few virial coefficients shows that superballs that are visually similar to cubes can have low-density equations of state closer to spheres than to cubes. Dense liquids of superballs display cubatic orientational order that extends over several particle lengths only for large q. Along the ordered, high-density equation of state, superballs with 1<q<3 exhibit clear evidence of a phase transition from a crystal state to a state with reduced long-ranged orientational order upon the reduction of density. For q≥3 , long-ranged orientational order persists until the melting transition. The width of the apparent coexistence region between the liquid and ordered, high-density phase decreases with q up to q=4.0. The structures of the high-density phases are examined using certain order parameters, distribution functions, and orientational correlation functions. We also find that a fixed simulation cell induces artificial phase transitions that are out of equilibrium. Current fabrication techniques allow for the synthesis of colloidal superballs and thus the phase behavior of such systems can be investigated experimentally.