Optimal packings of superballs

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.

Packing and self-assembly of truncated triangular bipyramids

Motivated by breakthroughs in the synthesis of faceted nano- and colloidal particles, as well as theoretical and computational studies of their packings, we investigate a family of truncated triangular bipyramids. We report dense periodic packings with small unit cells that were obtained via numerical and analytical optimization. The maximal packing fraction φ(max) changes continuously with the truncation parameter t. Eight distinct packings are identified based on discontinuities in the first and second derivatives of φ(max)(t). These packings differ in the number of particles in the fundamental domain (unit cell) and the type of contacts between the particles. In particular, we report two packings with four particles in the unit cell for which both φ(max)(t) and φ(max)'(t) are continuous and the discontinuity occurs in the second derivative only. In the self-assembly simulations that we perform for larger boxes with 2048 particles, only one out of eight packings is found to assemble. In addition, the degenerate quasicrystal reported previously for triangular bipyramids without truncation [Haji-Akbari et al., Phys. Rev. Lett. 107, 215702 (2011)] assembles for truncations as high as 0.45. The self-assembly propensities for the structures formed in the thermodynamic limit are explained using the isoperimetric quotient of the particles and the coordination number in the disordered fluid and in the assembled structure.

Shape-dependent plasmonic response and directed self-assembly in a new semiconductor building block, indium-doped cadmium oxide (ICO)

The influence of particle shape on plasmonic response and local electric field strength is well-documented in metallic nanoparticles. Morphologies such as rods, plates, and octahedra are readily synthesized and exhibit drastically different extinction spectra than spherical particles. Despite this fact, the influence of composition and shape on the optical properties of plasmonic semiconductor nanocrystals, in which free electrons result from heavy doping, has not been well-studied. Here, we report the first observation of plasmonic resonance in indium-doped cadmium oxide (ICO) nanocrystals, which exhibit the highest quality factors reported for semiconductor nanocrystals. Furthermore, we are able to independently control the shape and free electron concentration in ICO nanocrystals, allowing for the influence of shape on the optical response of a plasmonic semiconductor to be conclusively demonstrated. The highly uniform particles may be self-assembled into ordered single component and binary nanocrystal superlattices, and in thin films, exhibit negative permittivity in the near infrared (NIR) region, validating their use as a new class of tunable low-loss plasmonic building blocks for 3-D optical metamaterials.

Efficient linear programming algorithm to generate the densest lattice sphere packings

Finding the densest sphere packing in d-dimensional Euclidean space R(d) is an outstanding fundamental problem with relevance in many fields, including the ground states of molecular systems, colloidal crystal structures, coding theory, discrete geometry, number theory, and biological systems. Numerically generating the densest sphere packings becomes very challenging in high dimensions due to an exponentially increasing number of possible sphere contacts and sphere configurations, even for the restricted problem of finding the densest lattice sphere packings. In this paper we apply the Torquato-Jiao packing algorithm, which is a method based on solving a sequence of linear programs, to robustly reproduce the densest known lattice sphere packings for dimensions 2 through 19. We show that the TJ algorithm is appreciably more efficient at solving these problems than previously published methods. Indeed, in some dimensions, the former procedure can be as much as three orders of magnitude faster at finding the optimal solutions than earlier ones. We also study the suboptimal local density-maxima solutions (inherent structures or "extreme" lattices) to gain insight about the nature of the topography of the "density" landscape.

Effect of dimensionality on the percolation threshold of overlapping nonspherical hyperparticles

We study the effect of dimensionality on the percolation threshold η(c) of identical overlapping nonspherical convex hyperparticles in d-dimensional Euclidean space R(d). This is done by formulating a scaling relation for η(c) that is based on a rigorous lower bound [Torquato, J. Chem. Phys. 136, 054106 (2012)] and a conjecture that hyperspheres provide the highest threshold, for any d, among all convex hyperparticle shapes (that are not a trivial affine transformation of a hypersphere). This scaling relation also exploits the recently discovered principle that low-dimensional continuum percolation behavior encodes high-dimensional information. We derive an explicit formula for the exclusion volume v(ex) of a hyperparticle of arbitrary shape in terms of its d-dimensional volume v, surface area s, and radius of mean curvature R[over ¯] (or, equivalently, mean width). These basic geometrical properties are computed for a wide variety of nonspherical hyperparticle shapes with random orientations across all dimensions, including, among other shapes, various polygons for d=2, Platonic solids, spherocylinders, parallepipeds, and zero-volume plates for d=3 and their appropriate generalizations for d≥4. Using this information, we compute the lower bound and scaling relation for η(c) for this comprehensive set of continuum percolation models across dimensions. We demonstrate that the scaling relation provides accurate upper-bound estimates of the threshold η(c) across dimensions and becomes increasingly accurate as the space dimension increases.

Bending and elongation effects on the random packing of curved spherocylinders

Studies on the macroscopic and microscopic packing properties of nonconvex particles are scarce. As a common concave form, the curved spherocylinder is used in the simulations, and its bending and elongation effects on the random packings are investigated numerically with sphere assembly models and a relaxation algorithm. The aspect ratio is demonstrated to be the main factor regarding the packing density. However, at certain aspect ratios of low densities around 0.3-0.4, the density of curved spherocylinders may increase by 15% more than that of the straight ones, indicating that bending is also a contributor to the packing density. The excluded volume of the curved spherocylinder decreases with the increase of the bending angle, indicating that the excluded volume is applicable in explaining the bending effect on the packing density variation of nonconvex particles. The packings are verified to be randomly distributed in orientation with no significant layering or in-plane order. The local arrangements are further analyzed from the radial distribution function and contact results. The results show that the random packings of nonconvex particles have significant differences and richer characteristics on both the macroscopic and microscopic properties compared with convex objects.

Phase transitions and thermodynamic properties of dense assemblies of truncated nanocubes and cuboctahedra

Inspired by recent advances on the self-assembly of non-spherical nanoparticles, Monte Carlo simulations of the packing and thermodynamic properties of truncated nanocubes and cuboctahedra have been performed. The ergodicity problem was overcome by a modified Wang-Landau entropic sampling algorithm and equilibrium structural and thermodynamic properties were computed over a wide density range for both non-interacting and interacting particles. We found a structural transition from a simple cubic to a rhombohedral order when the degree of truncation exceeds a value of 0.9.

Self-assembly of colloidal cubes via vertical deposition

The vertical deposition technique for creating crystalline microstructures is applied for the first time to nonspherical colloids in the form of hollow silica cubes. Controlled deposition of the cubes results in large crystalline films with variable symmetry. The microstructures are characterized in detail with scanning electron microscopy and small-angle X-ray scattering. In single layers of cubes, distorted square to hexagonal ordered arrays are formed. For multilayered crystals, the intralayer ordering is predominantly hexagonal with a hollow site stacking, similar to that of the face centered cubic lattice for spheres. Additionally, a distorted square arrangement in the layers is also found to form under certain conditions. These crystalline films are promising for various applications such as photonic materials.

Entangled granular media

We study the geometrically induced cohesion of ensembles of granular "u particles" that mechanically entangle through particle interpenetration. We vary the length-to-width ratio l/w of the u particles and form them into freestanding vertical columns. In a laboratory experiment, we monitor the response of the columns to sinusoidal vibration (with peak acceleration Γ). Column collapse occurs in a characteristic time τ which follows the relation τ∝exp(Γ/Δ). Δ resembles an activation energy and is maximal at intermediate l/w. A simulation reveals that optimal strength results from competition between packing and entanglement.

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.

Shaping phases by phasing shapes

Incorporation of shape-shifting building blocks into self-assembled systems has emerged as a promising concept for dynamic structural control. The computational work by Nguyen et al. reported in this issue of ACS Nano examines the phase reconfigurations and kinetic pathways for systems built from shape-shifting building blocks. The studies illustrate several unique properties of such systems, including more efficient packings, novel structures that are distinctive from those obtained through conventional self-assembly, and reversible multistep shape-shifting pathways. The proposed assembly strategy is potentially applicable to a diverse range of systems because it relies on a change of geometrical constraints, which are common across all length scales. Recent developments in the areas of responsive materials and self-assembly methods provide feasible platforms for experimental realizations of shape-shifting reconfigurations; such systems might enable the next generation of dynamically switchable materials and reconfigurable devices.

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.

Continuous phase transformation in nanocube assemblies

The phase behavior of 3D assemblies of nanocubes in a ligand-rich solution upon solvent evaporation was experimentally investigated using small-angle x-ray scattering and electron microscopy. We observed a continuous transformation of assemblies between simple cubic and rhombohedral phases, where a variable angle of rhombohedral structure is determined by ligand thickness. We established a quantitative relationship between the particle shape evolution from cubes to quasispheres and the lattice distortion during the transformation, with a pathway exhibiting the highest known packing.

Sediments of soft spheres arranged by effective density

Colloidal sedimentation has been studied for decades from both thermodynamic and dynamic perspectives. In the present work, binary mixtures of colloidal spheres are observed to separate spontaneously into two distinct layers on sedimentation. Both layers have a high volume fraction and contain distinct compositions of particles. Although predicting these compositions using settling dynamics is challenging, here we show that the compositions are readily predicted thermodynamically by minimizing the gravitational energy of the system. As the random packing fraction of a mixture of spheres exceeds that of monodisperse spheres of either type, the mixture produces a denser suspension that forms the bottom phase. Experimentally, the use of charged particles and low-ionic-strength solutions provides interparticle repulsions that keep the packed particles mobile, avoiding a glassy state that would prevent particles from reaching their equilibrium configuration. We extend this work beyond binary systems, showing similar separated layers for a five-component mixture of particles.

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).

Dense-packing crystal structures of physical tetrahedra

We present a method for discovering dense packings of general convex hard particles and apply it to study the dense packing behavior of a one-parameter family of particles with tetrahedral symmetry representing a deformation of the ideal mathematical tetrahedron into a less ideal, physical, tetrahedron and all the way to the sphere. Thus, we also connect the two well-studied problems of sphere packing and tetrahedron packing on a single axis. Our numerical results uncover a rich optimal-packing behavior, compared to that of other continuous families of particles previously studied. We present four structures as candidates for the optimal packing at different values of asphericity, providing an atlas of crystal structures that might be observed in systems of nanoparticles with tetrahedral symmetry.

Robust algorithm to generate a diverse class of dense disordered and ordered sphere packings via linear programming

We have formulated the problem of generating dense packings of nonoverlapping, nontiling nonspherical particles within an adaptive fundamental cell subject to periodic boundary conditions as an optimization problem called the adaptive-shrinking cell (ASC) formulation [S. Torquato and Y. Jiao, Phys. Rev. E 80, 041104 (2009)]. Because the objective function and impenetrability constraints can be exactly linearized for sphere packings with a size distribution in d-dimensional Euclidean space R(d), it is most suitable and natural to solve the corresponding ASC optimization problem using sequential-linear-programming (SLP) techniques. We implement an SLP solution to produce robustly a wide spectrum of jammed sphere packings in R(d) for d=2, 3, 4, 5, and 6 with a diversity of disorder and densities up to the respective maximal densities. A novel feature of this deterministic algorithm is that it can produce a broad range of inherent structures (locally maximally dense and mechanically stable packings), besides the usual disordered ones (such as the maximally random jammed state), with very small computational cost compared to that of the best known packing algorithms by tuning the radius of the influence sphere. For example, in three dimensions, we show that it can produce with high probability a variety of strictly jammed packings with a packing density anywhere in the wide range [0.6, 0.7408...], where π/√18 = 0.7408... corresponds to the density of the densest packing. We also apply the algorithm to generate various disordered packings as well as the maximally dense packings for d=2, 4, 5, and 6. Our jammed sphere packings are characterized and compared to the corresponding packings generated by the well-known Lubachevsky-Stillinger (LS) molecular-dynamics packing algorithm. Compared to the LS procedure, our SLP protocol is able to ensure that the final packings are truly jammed, produces disordered jammed packings with anomalously low densities, and is appreciably more robust and computationally faster at generating maximally dense packings, especially as the space dimension increases.

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.

Distinctive features arising in maximally random jammed packings of superballs

Dense random packings of hard particles are useful models of granular media and are closely related to the structure of nonequilibrium low-temperature amorphous phases of matter. Most work has been done for random jammed packings of spheres and it is only recently that corresponding packings of nonspherical particles (e.g., ellipsoids) have received attention. Here we report a study of the maximally random jammed (MRJ) packings of binary superdisks and monodispersed superballs whose shapes are defined by |x1|2p+...+|xd|2p<or=1 with d=2 and 3, respectively, where p is the deformation parameter with values in the interval (0,infinity). As p increases from zero, one can get a family of both concave (0<p<0.5) and convex (p>or=0.5) particles with square symmetry (d=2), or octahedral and cubic symmetry (d=3). In particular, for p=1 the particle is a perfect sphere (circular disk) and for p-->infinity the particle is a perfect cube (square). We find that the MRJ densities of such packings increase dramatically and nonanalytically as one moves away from the circular-disk and sphere point (p=1). Moreover, the disordered packings are hypostatic, i.e., the average number of contacting neighbors is less than twice the total number of degrees of freedom per particle, and yet the packings are mechanically stable. As a result, the local arrangements of particles are necessarily nontrivially correlated to achieve jamming. We term such correlated structures "nongeneric." The degree of "nongenericity" of the packings is quantitatively characterized by determining the fraction of local coordination structures in which the central particles have fewer contacting neighbors than average. We also show that such seemingly "special" packing configurations are counterintuitively not rare. As the anisotropy of the particles increases, the fraction of rattlers decreases while the minimal orientational order as measured by the tetratic and cubatic order parameters increases. These characteristics result from the unique manner in which superballs break their rotational symmetry, which also makes the superdisk and superball packings distinctly different from other known nonspherical hard-particle packings.

Exact constructions of a family of dense periodic packings of tetrahedra

The determination of the densest packings of regular tetrahedra (one of the five Platonic solids) is attracting great attention as evidenced by the rapid pace at which packing records are being broken and the fascinating packing structures that have emerged. Here we provide the most general analytical formulation to date to construct dense periodic packings of tetrahedra with four particles per fundamental cell. This analysis results in six-parameter family of dense tetrahedron packings that includes as special cases recently discovered "dimer" packings of tetrahedra, including the densest known packings with density phi=4000/4671=0.856347... . This study strongly suggests that the latter set of packings are the densest among all packings with a four-particle basis. Whether they are the densest packings of tetrahedra among all packings is an open question, but we offer remarks about this issue. Moreover, we describe a procedure that provides estimates of upper bounds on the maximal density of tetrahedron packings, which could aid in assessing the packing efficiency of candidate dense packings.

Self-assembled structures of Gaussian nematic particles

We investigate the stable crystalline configurations of a nematic liquid crystal made of soft parallel ellipsoidal particles interacting via a repulsive, anisotropic Gaussian potential. For this purpose, we use genetic algorithms (GA) in order to predict all relevant and possible solid phase candidates into which this fluid can freeze. Subsequently we present and discuss the emerging novel structures and the resulting zero-temperature phase diagram of this system. The latter features a variety of crystalline arrangements, in which the elongated Gaussian particles in general do not align with any one of the high-symmetry crystallographic directions, a compromise arising from the interplay and competition between anisotropic repulsions and crystal ordering. Only at very strong degrees of elongation does a tendency of the Gaussian nematics to align with the longest axis of the elementary unit cell emerge.

Three-dimensional reconstruction of statistically optimal unit cells of polydisperse particulate composites from microtomography

In this paper, we present a systematic approach for characterization and reconstruction of statistically optimal representative unit cells of polydisperse particulate composites. Microtomography is used to gather rich three-dimensional data of a packed glass bead system. First-, second-, and third-order probability functions are used to characterize the morphology of the material, and the parallel augmented simulated annealing algorithm is employed for reconstruction of the statistically equivalent medium. Both the fully resolved probability spectrum and the geometrically exact particle shapes are considered in this study, rendering the optimization problem multidimensional with a highly complex objective function. A ten-phase particulate composite composed of packed glass beads in a cylindrical specimen is investigated, and a unit cell is reconstructed on massively parallel computers. Further, rigorous error analysis of the statistical descriptors (probability functions) is presented and a detailed comparison between statistics of the voxel-derived pack and the representative cell is made.

Dense packings of polyhedra: Platonic and Archimedean solids

Understanding the nature of dense particle packings is a subject of intense research in the physical, mathematical, and biological sciences. The preponderance of previous work has focused on spherical particles and very little is known about dense polyhedral packings. We formulate the problem of generating dense packings of nonoverlapping, nontiling polyhedra within an adaptive fundamental cell subject to periodic boundary conditions as an optimization problem, which we call the adaptive shrinking cell (ASC) scheme. This optimization problem is solved here (using a variety of multiparticle initial configurations) to find the dense packings of each of the Platonic solids in three-dimensional Euclidean space R3 , except for the cube, which is the only Platonic solid that tiles space. We find the densest known packings of tetrahedra, icosahedra, dodecahedra, and octahedra with densities 0.823..., 0.836..., 0.904..., and 0.947..., respectively. It is noteworthy that the densest tetrahedral packing possesses no long-range order. Unlike the densest tetrahedral packing, which must not be a Bravais lattice packing, the densest packings of the other nontiling Platonic solids that we obtain are their previously known optimal (Bravais) lattice packings. We also derive a simple upper bound on the maximal density of packings of congruent nonspherical particles and apply it to Platonic solids, Archimedean solids, superballs, and ellipsoids. Provided that what we term the "asphericity" (ratio of the circumradius to inradius) is sufficiently small, the upper bounds are relatively tight and thus close to the corresponding densities of the optimal lattice packings of the centrally symmetric Platonic and Archimedean solids. Our simulation results, rigorous upper bounds, and other theoretical arguments lead us to the conjecture that the densest packings of Platonic and Archimedean solids with central symmetry are given by their corresponding densest lattice packings. This can be regarded to be the analog of Kepler's sphere conjecture for these solids.The truncated tetrahedron is the only nonchiral Archimedean solid that is not centrally symmetric [corrected], the densest known packing of which is a non-lattice packing with density at least as high as 23/24=0.958 333... . We discuss the validity of our conjecture to packings of superballs, prisms, and antiprisms as well as to high-dimensional analogs of the Platonic solids. In addition, we conjecture that the optimal packing of any convex, congruent polyhedron without central symmetry generally is not a lattice packing. Finally, we discuss the possible applications and generalizations of the ASC scheme in predicting the crystal structures of polyhedral nanoparticles and the study of random packings of hard polyhedra.

Dense packings of the Platonic and Archimedean solids

Dense particle packings have served as useful models of the structures of liquid, glassy and crystalline states of matter, granular media, heterogeneous materials and biological systems. Probing the symmetries and other mathematical properties of the densest packings is a problem of interest in discrete geometry and number theory. Previous work has focused mainly on spherical particles-very little is known about dense polyhedral packings. Here we formulate the generation of dense packings of polyhedra as an optimization problem, using an adaptive fundamental cell subject to periodic boundary conditions (we term this the 'adaptive shrinking cell' scheme). Using a variety of multi-particle initial configurations, we find the densest known packings of the four non-tiling Platonic solids (the tetrahedron, octahedron, dodecahedron and icosahedron) in three-dimensional Euclidean space. The densities are 0.782..., 0.947..., 0.904... and 0.836..., respectively. Unlike the densest tetrahedral packing, which must not be a Bravais lattice packing, the densest packings of the other non-tiling Platonic solids that we obtain are their previously known optimal (Bravais) lattice packings. Combining our simulation results with derived rigorous upper bounds and theoretical arguments leads us to the conjecture that the densest packings of the Platonic and Archimedean solids with central symmetry are given by their corresponding densest lattice packings. This is the analogue of Kepler's sphere conjecture for these solids.