Optimal packings of superdisks and the role of symmetry

Computational design of anisotropic stealthy hyperuniform composites with engineered directional scattering properties

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here we generalize the Fourier-space-based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)1359-645410.1016/j.actamat.2017.09.053] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ[over ̃]_{_{V}}(k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform microstructures with χ[over ̃]_{_{V}}(k)=0 for k∈Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped, and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions and are nonhyperuniform along others. We find that while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the global symmetry of the resulting media, they can still possess structures at a local level that are nearly isotropic. Both the isotropic and anisotropic hyperuniform microstructures associated with the elliptical-disk, square, and rectangular Ω possess phase-inversion symmetry over certain range of volume fractions and a percolation threshold ϕ_{c}≈0.5. On the other hand, the directionally hyperuniform microstructures associated with the butterfly-shaped and lemniscate-shaped Ω do not possess phase-inversion symmetry and percolate along certain directions at much lower volume fractions. We also apply our general procedure to construct stealthy nonhyperuniform systems. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures via the shape, size, and symmetries of Ω, which is crucial to engineering directional optical, transport, and mechanical properties of two-phase composite media.

Design rules for 2D field mediated assembly of different shaped colloids into diverse microstructures

Assembling different shaped particles into ordered microstructures is an open challenge in creating multifunctional particle-based materials and devices. Here, we report the two-dimensional (2D) AC electric field mediated assembly of different shaped colloidal particles into amorphous, liquid crystalline, and crystalline microstructures. Particle shapes investigated include disks, ellipses, squares, and rectangles, which show how systematic variations in anisotropy and corner curvature determine the number and type of resulting microstructures. AC electric fields induce dipolar interactions to control particle positional and orientational order. Microstructural states are determined via particle tracking to compute order parameters, which agree with computer simulations and show how particle packing and dipolar interactions together produce each structure. Results demonstrate how choice of particle shape and field conditions enable kinetically viable routes to assemble nematic, tetratic, and smectic liquid crystal structures as well as crystals with stretched 4- and 6-fold symmetry. Results show it is possible to assemble all corresponding hard particle phases, but also show how dipolar interactions influence and produce additional microstructures. Our findings provide design rules for the assembly of diverse microstructures of different shaped particles in AC electric fields, which could enable scalable and reconfigurable particle-based materials, displays, and printing technologies.

Finding the grain size distribution that produces the densest arrangement in frictional sphere packings: Revisiting and rediscovering the century-old Fuller and Thompson distribution

By means of discrete-element methods, we investigate the joint effects of the grain size distribution (GSD) and contact friction on the structure of three-dimensional samples composed of spherical grains. Specifically, we compress these systems isotropically until jamming and then analyze their structure in terms of density, connectivity, coefficients of uniformity and curvature, and parameters of grading entropy. Our study focuses on power-law GSDs and particularly on the Fuller and Thompson distribution, proposed over a century ago. First, we show that, among the set of GSDs investigated, this particular distribution produces the densest and best-connected systems, falsifying a conjecture recently posed in the literature. Second, we find that the jamming packing fraction can be accurately predicted as a function of simple descriptors of the GSD, but among these descriptors the granular entropy concept proves to be the most useful. This allows for an alternative interpretation of both jamming and grading entropy concepts. Finally, we compare the Fuller and Thompson distribution with two well-known GSDs: that of the Apollonian sphere packing and that towards which granular systems evolve after intensive grain fracturing. Surprisingly, we find that these three GSDs are practically coincident in the limit of large size spans, despite having been introduced or discovered in different scientific contexts (i.e., engineering, mathematics, and earth sciences, respectively).

Self-assembly of colloidal cube superstructures with critical Casimir attractions

The structure of self-assembled materials is determined by the shape and interactions of the building blocks. Here, we investigate the self-assembly of colloidal 'superballs', i.e. cubes with rounded corners, by temperature-tunable critical Casimir forces to obtain insight into the coupling of a cubic shape and short range attractions. The critical Casimir force is a completely reversible and controllable attraction that arises in a near-critical solvent mixture. Using confocal microscopy and particle tracking, we follow the self-assembly dynamics and structural transition in a quasi-2D system. At low attraction, we observe the formation of small clusters with square symmetry. When the attraction is increased, a transition to a rhombic Λ₁-lattice is observed. We explain our findings by the change in contact area at faces and corners of the building blocks combined with the increase in attraction strength and range of the critical Casimir force. Our results show that the coupling between the rounded cubic shape and short-range attraction plays a crucial role for the superstructures that form and provide new insights for the active assembly control of micro and nanocubes.

Hard superellipse phases: particle shape anisotropy & curvature

We report computer simulations of two-dimensional convex hard superellipse particle phases vs. particle shape parameters including aspect ratio, corner curvature, and sidewall curvature. Shapes investigated include disks, ellipses, squares, rectangles, and rhombuses, as well as shapes with non-uniform curvature including rounded squares, rounded rectangles, and rounded rhombuses. Using measures of orientational order, order parameters, and a novel stretched bond orientational order parameter, we systematically identify particle shape properties that determine liquid crystal and crystalline phases including their coarse boundaries and symmetry. We observe phases including isotropic, nematic, tetratic, plastic crystals, square crystals, and hexagonal crystals (including stretched variants). Our results catalog known benchmark shapes, but include new shapes that also interpolate between known shapes. Our results indicate design rules for particle shapes that determine two-dimensional liquid, liquid crystalline, and crystalline microstructures that can be realized via particle assembly.

Three-step melting of hard superdisks in two dimensions

We explore the link between the melting scenarios of two-dimensional systems of hard disks and squares through replica-exchange Monte Carlo simulations of hard superdisks. The well-known melting scenarios are observed in the disk and square limits, while we observe an unusual three-step scenario for dual shapes. We find that two mesophases mediate the melting: a hexatic phase and another fluid phase with a D_{2} local symmetry, we call it rhombatic, where both bond and particle orientational orders are quasi-long-range. Our results show that not only can the melting process of liquid-crystal forming molecules be complicated, where elongated shapes stabilize several mesophases, but also that of anisotropic quasispherical molecules.

Dense packings of hard circular arcs

This work investigates dense packings of congruent hard infinitesimally thin circular arcs in the two-dimensional Euclidean space. It focuses on those denotable as major whose subtended angle θ∈(π,2π]. Differently than those denotable as minor whose subtended angle θ∈[0,π], it is impossible for two hard infinitesimally thin circular arcs with θ∈(π,2π] to arbitrarily closely approach once they are arranged in a configuration, e.g., on top of one another, replicable ad infinitum without introducing any overlap. This makes these hard concave particles, in spite of being infinitesimally thin, most densely pack with a finite number density. This raises the question as to what are these densest packings and what is the number density that they achieve. Supported by Monte Carlo numerical simulations, this work shows that one can analytically construct compact closed circular groups of hard major circular arcs in which a specific, θ-dependent, number of them (counter) clockwise intertwine. These compact closed circular groups then arrange on a triangular lattice. These analytically constructed densest-known packings are compared to corresponding results of Monte Carlo numerical simulations to assess whether they can spontaneously turn up.

Demixing and tetratic ordering in some binary mixtures of hard superellipses

We examine the fluid phase behavior of binary mixtures of hard superellipses using the scaled particle theory. The superellipse is a general two-dimensional convex object that can be tuned between the elliptical and rectangular shapes continuously at a given aspect ratio. We find that the shape of the particle affects strongly the stability of isotropic, nematic, and tetratic phases in the mixture even if the side lengths of both species are fixed. While the isotropic-isotropic demixing transition can be ruled out using the scaled particle theory, the first order isotropic-nematic and the nematic-nematic demixing transition can be stabilized with strong fractionation between the components. It is observed that the demixing tendency is strongest in small rectangle-large ellipse mixtures. Interestingly, it is possible to stabilize the tetratic order at lower densities in the mixture of hard squares and rectangles where the long rectangles form a nematic phase, while the squares stay in the tetratic order.

Intrinsic Effect of Nanoparticles on the Mechanical Rupture of Doubled-Shell Colloidal Capsule via In Situ TEM Mechanical Testing and STEM Interfacial Analysis

The discovery of Pickering emulsion templated assembly enables the design of a hybrid colloidal capsule with engineered properties. However, the underlying mechanisms by which nanoparticles affect the mechanical properties of the shell are poorly understood. Herein, in situ mechanical compression on the transmission electron microscope and aberration-corrected scanning transmission microscope are unprecedentedly implemented to study the intrinsic effect of nanoparticles on the mechanical properties of the calcium carbonate (CaCO₃ )-decorated silica (SiO₂ ) colloidal capsule. The stiff and brittle nature of the colloidal capsule is due to the interfacial chemical bonding between the CaCO₃ nanoparticles and SiO₂ inner shell. Such bonding strengthens the mechanical strength of the SiO₂ shell (166 ± 14 nm) from the colloidal capsule compared to the thicker single SiO₂ shell (310 ± 70 nm) from the silica hollow sphere. At elevated temperature, this interfacial bonding accelerates the formation of the single calcium silicate shell, causing shell morphology transformation and yielding significantly enhanced mechanical strength by 30.9% and ductility by 94.7%. The superior thermal durability of the heat-treated colloidal capsule holds great potential for the fabrication of the functional additives that can be applied in the wide range of applications at elevated temperatures.

Scattering from colloidal cubic silica shells: Part I, particle form factors and optical contrast variation

HYPOTHESIS

Colloidal cubic silica shells, prepared from cuprous oxide cubes, with a typical size of 100 nm are promising model particles for scattering studies on dilute, as well as concentrated fluids, of non-spherical colloids.

EXPERIMENTS

Small angle X-ray scattering, and static light scattering are employed to determine form factors of cubic silica shells and silica covered cuprous oxide cubes. Contrast variation experiments are performed to assess the refractive index and optical homogeneity of the cubic silica shells, which is important for the extension of the scattering study to concentrated dispersions of cubic shells in Part II (Dekker, submitted for publication).

RESULTS

The experimental form factors, which compare well to theoretical form factors, manifest cubic silica shells that are dispersed as single stable colloids with a shape intermediate between a sphere and a perfect cube. Contrast variation demonstrates that the silica shells are optically homogeneous, with a refractive index that is independent of the shell thickness. The results presented here open up the possibility to extract structure factors from light scattering measurements on concentrated cube dispersions in Part II.

Random sequential adsorption of particles with tetrahedral symmetry

We study random sequential adsorption (RSA) of a class of solids that can be obtained from a cube by specific cutting of its vertices, in order to find out how the transition from tetrahedral to octahedral symmetry affects the densities of the resulting jammed packings. We find that in general solids of octahedral symmetry form less dense packing; however, the lowest density was obtained for the packing built of tetrahedra. The densest packing is formed by a solid close to a tetrahedron but with vertices and edges slightly cut. Its density is θ_{max}=0.41278±0.00059 and is higher than the mean packing fraction of spheres or cuboids but is lower than that for the densest RSA packings built of ellipsoids or spherocylinders. The density autocorrelation function of the studied packings is typical for random media and vanishes very quickly with distance.

Convectively Assembled Monolayers of Colloidal Cubes: Evidence of Optimal Packings

We employ a system of cubic colloids with rounded corners to study the close-packed monolayers that form via convective assembly. We show that by controlled solvent evaporation large densely packed monolayers of colloidal cubes are obtained. Using scanning electron microscopy and particle-tracking algorithms, we investigate the local order in detail and show that the obtained monolayers possess their predicted close-packed optimal packings, the Λ₀-lattice and the Λ₁-lattice, as well as the simple square-lattice and disordered packings. We further show that shape details of the cube corners are important for the final packing symmetry, where the frequency of the Λ₁-lattice increases with decreasing roundness of the corners, whereas the frequency of the Λ₀-lattice is unaffected. The formation of both optimal packings is found to be a consequence of the out-of-equilibrium formation process, which leads to small shifts in rows of cubes, thereby transforming the Λ₁-lattice into the Λ₀-lattice.

Continuum percolation-based tortuosity and thermal conductivity of soft superball systems: shape dependence from octahedra via spheres to cubes

Understanding the effect of particle shape on the percolation threshold, tortuosity and thermal conductivity of soft (geometrical overlapping) particle systems is very crucial for the design and optimization of such materials, including colloids, polymers, and porous and fracture media. In this work, we first combine the excluded-volume theory with the Monte Carlo simulations to determine the percolation threshold for a family of soft superballs, the shape of which interpolates between octahedra and cubes via spheres. Then, we propose two continuum percolation-based models to respectively obtain the tortuosity and effective thermal conductivity of soft superball systems considering their percolation behavior, where monodisperse overlapping superballs are uniformly distributed in a homogeneous solid matrix. Specifically, both models cover the whole feasible range of superball volume fractions, including near the percolation threshold. Comparison with extensive experimental, numerical and theoretical results confirms that the present models are capable of precisely predicting the percolation threshold, tortuosity and thermal conductivity of such systems. Furthermore, we apply the proposed models to probe the influence of particle shape on these important parameters. Our results show that the decreasing percolation threshold and tortuosity as soft particles become more anisotropic is consistent with increasing conductivity. It suggests that the anisotropic-shaped inclusion phase is more conducting than the spherical inclusion phase. The present theoretical strategies and conclusions may provide sound guidance for the synthesis of colloidal and polymer superballs.

Transition of Dielectrophoresis-Assembled 2D Crystals to Interlocking Structures under a Magnetic Field

Aspherical cubic hematite colloids with cylindrical arms protruding from each face, referred to as "hexapods", were assembled via negative dielectrophoresis and then manipulated using an applied magnetic field. Upon application of an ac electric field, the hexapods aligned in close-packed linear chains parallel to the field direction. The chains then aggregated to the center of the device, with adjacent chains separated by distances approximately equal to twice the arm length. The resulting open packing structure exhibited cmm plane group symmetry due to the obstruction of arms, with a high density of incorporated defects. Subsequent application of a magnetic field to the dielectrophoresis (DEP)-assembled structure was found to anneal the colloidal crystal by reorienting the hexapods to align their intrinsic magnetic dipoles with the magnetic field direction. During reorganization, the colloidal packing density was found to decrease by more than 10% at both the center and edges of the crystal, accompanied by a significant loss of ordering, prior to redensification of the 2D lattice with fewer defects. Reorganization at the edge was 1.5 times faster than at the center, consistent with the need for cooperative colloidal motion to remove defects at the centers of the crystals.

Observation of solid-solid transitions in 3D crystals of colloidal superballs

Self-organization in anisotropic colloidal suspensions leads to a fascinating range of crystal and liquid crystal phases induced by shape alone. Simulations predict the phase behaviour of a plethora of shapes while experimental realization often lags behind. Here, we present the experimental phase behaviour of superball particles with a shape in between that of a sphere and a cube. In particular, we observe the formation of a plastic crystal phase with translational order and orientational disorder, and the subsequent transformation into rhombohedral crystals. Moreover, we uncover that the phase behaviour is richer than predicted, as we find two distinct rhombohedral crystals with different stacking variants, namely hollow-site and bridge-site stacking. In addition, for slightly softer interactions we observe a solid-solid transition between the two. Our investigation brings us one step closer to ultimately controlling the experimental self-assembly of superballs into functional materials, such as photonic crystals.

Depletion-driven crystallization of cubic colloids sedimented on a surface

Cubic colloids, sedimented on a surface and immersed in a solution of depletant molecules, were modeled with a family of shapes which smoothly varies from squares to circles. Using Wang-Landau simulations with expanded ensembles, we observe the formation of rhombic lattices, square lattices, hexagonal lattices, and a fluid phase. This systematic investigation includes locating transitions between all combinations of the three lattice structures upon changing the shape and transitions between the fluid and crystal upon changing the depletant concentration. The rhombic lattice deforms smoothly between square-like and hexagonal-like angles, depending on both the shape and the depletant concentration. Our results on the effect of the depletant concentration, depletant size, and colloid shape to influence the stability of the fluid and the lattice structures may help guide experimental studies with recently synthesized cubic colloids.

The Synthesis of Organic Molecules of Intrinsic Microporosity Designed to Frustrate Efficient Molecular Packing

Efficient reactions between fluorine-functionalised biphenyl and terphenyl derivatives with catechol-functionalised terminal groups provide a route to large, discrete organic molecules of intrinsic microporosity (OMIMs) that provide porous solids solely by their inefficient packing. By altering the size and substituent bulk of the terminal groups, a number of soluble compounds with apparent BET surface areas in excess of 600 m(2)  g(-1) are produced. The efficiency of OMIM structural units for generating microporosity is in the order: propellane>triptycene>hexaphenylbenzene>spirobifluorene>naphthyl=phenyl. The introduction of bulky hydrocarbon substituents significantly enhances microporosity by further reducing packing efficiency. These results are consistent with findings from previously reported packing simulation studies. The introduction of methyl groups at the bridgehead position of triptycene units reduces intrinsic microporosity. This is presumably due to their internal position within the OMIM structure so that they occupy space, but unlike peripheral substituents they do not contribute to the generation of free volume by inefficient packing.

Shape-sensitive crystallization in colloidal superball fluids

Guiding the self-assembly of materials by controlling the shape of the individual particle constituents is a powerful approach to material design. We show that colloidal silica superballs crystallize into canted phases in the presence of depletants. Some of these phases are consistent with the so-called "Λ1" lattice that was recently predicted as the densest packing of superdisks. As the size of the depletant is reduced, however, we observe a transition to a square phase. The differences in these entropically stabilized phases result from an interplay between the size of the depletants and the fine structure of the superball shape. We find qualitative agreement of our experimental results both with a phase diagram computed on the basis of the volume accessible to the depletants and with simulations. By using a mixture of depletants, one of which is thermosensitive, we induce solid-to-solid phase transitions between square and canted structures. The use of depletant size to leverage fine features of the shape of particles in driving their self-assembly demonstrates a general and powerful mechanism for engineering novel materials.

Random sequential adsorption of starlike particles

Random packing of surfaceless starlike particles built of 3 to 50 line segments was studied using random sequential adsorption algorithm. Numerical simulations allow us to determine saturated packing densities as well as the first two virial expansion coefficients for such objects. Measured kinetics of the packing growth supports the power law known to be valid for particles with a finite surface; however, the dependence of the exponent in this law on the number of star arms is unexpected. The density autocorrelation function shows fast superexponential decay as for disks, but the typical distance between closest stars is much smaller than between disks of the similar size, especially for a small number of arms.

Understanding shape entropy through local dense packing

Many natural systems are structured by the ordering of repeated, distinct shapes. Understanding how this happens is difficult because shape affects structure in two ways. One is how the shape of a cell or nanoparticle, for example, affects its surface, chemical, or other intrinsic properties. The other is an emergent, entropic effect that arises from the geometry of the shape itself, which we term “shape entropy,” and is not well understood. In this paper, we determine how shape entropy affects structure. We quantify the mechanism and determine when shape entropy competes with intrinsic shape effects. Our results show that in a wide class of systems, shape affects bulk structure because crowded particles optimize their local packing.

Understanding shape entropy through local dense packing

Entropy drives the phase behavior of colloids ranging from dense suspensions of hard spheres or rods to dilute suspensions of hard spheres and depletants. Entropic ordering of anisotropic shapes into complex crystals, liquid crystals, and even quasicrystals was demonstrated recently in computer simulations and experiments. The ordering of shapes appears to arise from the emergence of directional entropic forces (DEFs) that align neighboring particles, but these forces have been neither rigorously defined nor quantified in generic systems. Here, we show quantitatively that shape drives the phase behavior of systems of anisotropic particles upon crowding through DEFs. We define DEFs in generic systems and compute them for several hard particle systems. We show they are on the order of a few times the thermal energy ([Formula: see text]) at the onset of ordering, placing DEFs on par with traditional depletion, van der Waals, and other intrinsic interactions. In experimental systems with these other interactions, we provide direct quantitative evidence that entropic effects of shape also contribute to self-assembly. We use DEFs to draw a distinction between self-assembly and packing behavior. We show that the mechanism that generates directional entropic forces is the maximization of entropy by optimizing local particle packing. We show that this mechanism occurs in a wide class of systems and we treat, in a unified way, the entropy-driven phase behavior of arbitrary shapes, incorporating the well-known works of Kirkwood, Onsager, and Asakura and Oosawa.

Random packing of regular polygons and star polygons on a flat two-dimensional surface

Random packing of unoriented regular polygons and star polygons on a two-dimensional flat continuous surface is studied numerically using random sequential adsorption algorithm. Obtained results are analyzed to determine the saturated random packing ratio as well as its density autocorrelation function. Additionally, the kinetics of packing growth and available surface function are measured. In general, stars give lower packing ratios than polygons, but when the number of vertexes is large enough, both shapes approach disks and, therefore, properties of their packing reproduce already known results for disks.

Particle shape anisotropy in pickering emulsions: cubes and peanuts

We have investigated the effect of particle shape in Pickering emulsions by employing, for the first time, cubic and peanut-shaped particles. The interfacial packing and orientation of anisotropic microparticles are revealed at the single-particle level by direct microscopy observations. The uniform anisotropic hematite microparticles adsorb irreversibly at the oil-water interface in monolayers and form solid-stabilized o/w emulsions via the process of limited coalescence. Emulsions were stable against further coalescence for at least 1 year. We found that cubes assembled at the interface in monolayers with a packing intermediate between hexagonal and cubic and average packing densities of up to 90%. Local domains displayed densities even higher than theoretically achievable for spheres. Cubes exclusively orient parallel with one of their flat sides at the oil-water interface, whereas peanuts preferentially attach parallel with their long side. Those peanut-shaped microparticles assemble in locally ordered, interfacial particle stacks that may interlock. Indications for long-range capillary interactions were not found, and we hypothesize that this is related to the observed stable orientations of cubes and peanuts that marginalize deformations of the interface.

Kinetics of random sequential adsorption of nearly spherically symmetric particles

Kinetics of random sequential adsorption (RSA) of disks on flat, two-dimensional surfaces is governed by a power law with exponent -1/2. The study has shown that for RSA of nearly spherically symmetric particles this exponent is -1/3, whereas other characteristics typically measured in RSA simulations approach values known for disks with the increase of symmetry of the particles.

Kinetics of deposition of oriented superdisks

We use numerical Monte Carlo simulation to study the kinetics of the deposition of oriented superdisks, bounded by the Lame curves of the form |x|(2p)+|y|(2p)=1 on a regular planar substrate. Recently, it was shown that the maximum packing density as well as jamming density ρ(J) exhibit a discontinuous derivative at p=0.5 when the shape changes from convex to concave form. By careful examination of the late-stage approach to the jamming limit, we find that the leading term in the temporal development is also nonanalytic at p=0.5 and offer heuristic excluded-area arguments for this behavior.

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.

Maximally dense packings of two-dimensional convex and concave noncircular particles

Dense packings of hard particles have important applications in many fields, including condensed matter physics, discrete geometry, and cell biology. In this paper, we employ a stochastic search implementation of the Torquato-Jiao adaptive-shrinking-cell (ASC) optimization scheme [Nature (London) 460, 876 (2009)] to find maximally dense particle packings in d-dimensional Euclidean space R(d). While the original implementation was designed to study spheres and convex polyhedra in d≥3, our implementation focuses on d=2 and extends the algorithm to include both concave polygons and certain complex convex or concave nonpolygonal particle shapes. We verify the robustness of this packing protocol by successfully reproducing the known putative optimal packings of congruent copies of regular pentagons and octagons, then employ it to suggest dense packing arrangements of congruent copies of certain families of concave crosses, convex and concave curved triangles (incorporating shapes resembling the Mercedes-Benz logo), and "moonlike" shapes. Analytical constructions are determined subsequently to obtain the densest known packings of these particle shapes. For the examples considered, we find that the densest packings of both convex and concave particles with central symmetry are achieved by their corresponding optimal Bravais lattice packings; for particles lacking central symmetry, the densest packings obtained are nonlattice periodic packings, which are consistent with recently-proposed general organizing principles for hard particles. Moreover, we find that the densest known packings of certain curved triangles are periodic with a four-particle basis, and we find that the densest known periodic packings of certain moonlike shapes possess no inherent symmetries. Our work adds to the growing evidence that particle shape can be used as a tuning parameter to achieve a diversity of packing structures.

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.

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.

Deposition of general ellipsoidal particles

We present a systematic overview of granular deposits composed of ellipsoidal particles with different particle shapes and size polydispersities. We study the density and anisotropy of such deposits as functions of small to moderate size polydispersity and two shape parameters that fully describe the shape of a general ellipsoid. Our results show that, while shape influences significantly the macroscopic properties of the deposits, polydispersity in the studied range plays apparently a secondary role. The density attains a maximum for a particular family of nonsymmetrical ellipsoids, larger than the density observed for prolate or oblate ellipsoids. As for anisotropy measures, the contact forces are increasingly preferred along the vertical direction as the shape of the particles deviates from a sphere. The deposits are constructed by means of a molecular dynamics method, where the contact forces are efficiently and accurately computed. The main results are discussed in the light of applications for porous media models and sedimentation processes.

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.

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.

Morphologically controlled synthesis of colloidal upconversion nanophosphors and their shape-directed self-assembly

We report a one-pot chemical approach for the synthesis of highly monodisperse colloidal nanophosphors displaying bright upconversion luminescence under 980 nm excitation. This general method optimizes the synthesis with initial heating rates up to 100 °C/minute generating a rich family of nanoscale building blocks with distinct morphologies (spheres, rods, hexagonal prisms, and plates) and upconversion emission tunable through the choice of rare earth dopants. Furthermore, we employ an interfacial assembly strategy to organize these nanocrystals (NCs) into superlattices over multiple length scales facilitating the NC characterization and enabling systematic studies of shape-directed assembly. The global and local ordering of these superstructures is programmed by the precise engineering of individual NC's size and shape. This dramatically improved nanophosphor synthesis together with insights from shape-directed assembly will advance the investigation of an array of emerging biological and energy-related nanophosphor applications.

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.

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.

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.

Random close packing of disks and spheres in confined geometries

Studies of random close packing of spheres have advanced our knowledge about the structure of systems such as liquids, glasses, emulsions, granular media, and amorphous solids. In confined geometries, the structural properties of random-packed systems will change. To understand these changes, we study random close packing in finite-sized confined systems, in both two and three dimensions. Each packing consists of a 50-50 binary mixture with particle size ratio of 1.4. The presence of confining walls significantly lowers the overall maximum area fraction (or volume fraction in three dimensions). A simple model is presented, which quantifies the reduction in packing due to wall-induced structure. This wall-induced structure decays rapidly away from the wall, with characteristic length scales comparable to the small particle diameter.

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.

Random sequential adsorption of oriented superdisks

In this work we extend recent study of the properties of the dense packing of "superdisks," by Y. Jiao [Phys. Rev. Lett. 100, 245504 (2008)] to the jammed state formed by these objects in random sequential adsorption. The superdisks are two-dimensional shapes bound by the curves of the form |x|2p+|y|2p=1, with p>0. We use Monte Carlo simulations and theoretical arguments to establish that p=1/2 is a special point at which the jamming density, rhoJ(p), has a discontinuous derivative as a function of p . The existence of this point can be also argued for by a phenomenological excluded-area argument.

Optimal packings of superballs

Dense hard-particle packings are intimately related to the structure of low-temperature phases of matter and are useful models of heterogeneous materials and granular media. Most studies of the densest packings in three dimensions have considered spherical shapes, and it is only more recently that nonspherical shapes (e.g., ellipsoids) have been investigated. Superballs (whose shapes are defined by |x1|2p+|x2|2p+|x3|2p<or=1) provide a versatile family of convex particles (p>or=0.5) with both cubic-like and octahedral-like shapes as well as concave particles (0<p<0.5) with octahedral-like shapes. In this paper, we provide analytical constructions for the densest known superball packings for all convex and concave cases. The candidate maximally dense packings are certain families of Bravais lattice packings (in which each particle has 12 contacting neighbors) possessing the global symmetries that are consistent with certain symmetries of a superball. We also provide strong evidence that our packings for convex superballs (p>or=0.5) are most likely the optimal ones. The maximal packing density as a function of p is nonanalytic at the sphere point (p=1) and increases dramatically as p moves away from unity. Two more nontrivial nonanalytic behaviors occur at pc*=1.150 9... and po*=ln 3/ln 4=0.792 4... for "cubic" and "octahedral" superballs, respectively, where different Bravais lattice packings possess the same densities. The packing characteristics determined by the broken rotational symmetry of superballs are similar to but richer than their two-dimensional "superdisk" counterparts [Y. Jiao, Phys. Rev. Lett. 100, 245504 (2008)] and are distinctly different from that of ellipsoid packings. Our candidate optimal superball packings provide a starting point to quantify the equilibrium phase behavior of superball systems, which should deepen our understanding of the statistical thermodynamics of nonspherical-particle systems.