The Densest Packing of AB Binary Hard-Sphere Homogeneous Compounds across all Size Ratios

Tuning assembly structures of hard shapes in confinement via interface curvature

Assembly in confinement is a problem of great interest in colloidal structure design, plasmonics, photonics, and industrial packaging. Along with the range of design choices provided by particle shape and attraction or repulsion, confined systems add an additional layer of complexity through the interactions between particles and the container holding them. The range of possible behaviors produced by these systems remains largely unexplored, yet has profound consequences on the resultant assembled structure. Here, we address this problem by exploring how the assembly of hard tetrahedral particles is affected by a spherical container. We simulate particle assemblies in containers holding 4 to 10 000 particles and analyze the range of resultant structures. We find that the presence of a curved wall causes organization into distinct concentric shells in containers holding up to thousands of particles. In addition, we see that wall curvature affects structural motifs in systems as large as 10 000 particles, promoting local environments that maximally conform to the wall and providing a seed for the propagation of these motifs into the interior of the container. Through this work, we show how confining interfaces can be used to promote the assembly of structures markedly distinct from those seen in the more commonly studied bulk systems.

Inherent structure landscape of hard spheres confined to narrow cylindrical channels

The inherent structure landscape for a system of hard spheres confined to a hard cylindrical channel, such that spheres can only contact their first and second neighbors, is studied using an analytical model that extends previous results [Phys. Rev. Lett. 115, 025702 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.025702] to provide a comprehensive picture of jammed packings over a range of packing densities. In the model, a packing is described as an arrangement of k helical sections, separated by defects, that have alternating helical twist directions and where all spheres satisfy local jamming constraints. The structure of each helical section is determined by a single helical twist angle, and a jammed packing is obtained by minimizing the length of the channel per particle with respect to the k helical section angles. An analysis of a small system of N=20 spheres shows that the basins on the inherent structure landscape associated with these helical arrangements split into a number of distinct jammed states separated by low barriers giving rise to a degree of hierarchical organization. The model accurately predicts the geometric properties of packings generated using the Lubachevsky and Stillinger compression scheme (N=10^{4}) and provides insight into the nature of the probability distribution of helical section lengths.

Densest ternary sphere packings

We present our exhaustive exploration of the densest ternary sphere packings (DTSPs) for 45 radius ratios and 237 kinds of compositions, which is a packing problem of three kinds of hard spheres with different radii, under periodic boundary conditions by a random structure searching method. To efficiently explore DTSPs we further develop the searching method based on the piling-up and iterative balance methods [Koshoji et al., Phys. Rev. E 103, 023307 (2021)2470-004510.1103/PhysRevE.103.023307]. The unbiased exploration identifies diverse 38 putative DTSPs appearing on phase diagrams in which 37 DTSPs of them are discovered in the study. The structural trend of DTSPs changes depending especially on the radius of small spheres. In case that the radius of small spheres is relatively small, structures of many DTSPs can be understood as derivatives of densest binary sphere packings (DBSPs), while characteristic structures specific to the ternary system emerge as the radius of small spheres becomes larger. In addition to DTSPs, we reveal a lot of semi-DTSPs (SDTSPs) which are obtained by excluding DBSPs in the calculation of phase diagram, and investigate the correspondence of DTSPs and SDTSPs with real crystals based on the space group, showing a considerable correspondence of SDTSPs having high symmetries with real crystals including Cu_{2}GaSr and ThCr_{2}Si_{2} structures. Our study suggests that the diverse structures of DBSPs, DTSPs, and SDTSPs can be effectively used as structural prototypes for searching complex crystal structures.

Symmetry-derived structure directing agents for two-dimensional crystals of arbitrary colloids

We derive properties of self-assembling rings which can template the organization of an arbitrary colloid into any periodic symmetry in two Euclidean dimensions. By viewing this as a tiling problem, we illustrate how the shape and chemical patterning of these rings are derivable, and are explicitly reflected by the symmetry group's orbifold symbol. We performed molecular dynamics simulations to observe their self-assembly and found 5 different characteristics which could be easily rationalized on the basis of this symbol. These include systems which undergo chiral phase separation, are addressably complex, exhibit self-limiting growth into clusters, form ordered "rods" in only one-dimension akin to a smectic phase, and those from symmetry groups which are pluripotent and allow one to select rings which exhibit different behaviors. We discuss how the curvature of the ring's edges plays an integral role in achieving correct self-assembly, and illustrate how to obtain these shapes. This provides a method for patterning colloidal systems at interfaces without explicitly programming this information onto the colloid itself.

Diverse densest binary sphere packings and phase diagram

We revisit the densest binary sphere packings (DBSPs) under periodic boundary conditions and present an updated phase diagram, including newly found 12 putative densest structures over the x-α plane, where x is the relative concentration and α is the radius ratio of the small and large spheres. To efficiently explore the DBSPs, we develop an unbiased random search approach based on both the piling-up method to generate initial structures in an unbiased way and the iterative balance method to optimize the volume of a unit cell while keeping the overlap of hard spheres minimized. With those two methods, we have discovered 12 putative DBSPs and thereby the phase diagram is updated, while our results are consistent with those of a previous study [Hopkins et al., Phys. Rev. E 85, 021130 (2012)]PLEEE81539-375510.1103/PhysRevE.85.021130 with a small correction for the case of 12 or fewer spheres in the unit cell. Five of the discovered 12 DBSPs are identified in the small radius range of 0.42≤α≤0.50, where several structures are competitive to each other with respect to packing fraction. Through the exhaustive search, diverse dense packings are discovered and, accordingly, we find that packing structures achieve high packing fractions by introducing distortion and/or combining a few local dense structural units. Furthermore, we investigate the correspondence of the DBSPs with crystals based on the space group. The result shows that many structural units in real crystals, e.g., LaH_{10} and SrGe_{2-δ} being high-pressure phases, can be understood as DBSPs. The correspondence implies that the densest sphere packings can be used effectively as structural prototypes for searching complex crystal structures, especially for high-pressure phases.

Phase Diagram and Structure Map of Binary Nanoparticle Superlattices from a Lennard-Jones Model

A first-principles prediction of the binary nanoparticle phase diagram assembled by solvent evaporation has eluded theoretical approaches. In this paper, we show that a binary system interacting through the Lennard-Jones (LJ) potential contains all experimental phases in which nanoparticles are effectively described as quasi hard spheres. We report a phase diagram consisting of 53 equilibrium phases, whose stability is quite insensitive to the microscopic details of the potentials, thus giving rise to some type of universality. Furthermore, we show that binary lattices may be understood as consisting of certain particle clusters, i.e., motifs, that provide a generalization of the four conventional Frank-Kasper polyhedral units. Our results show that metastable phases share the very same motifs as equilibrium phases. We discuss the connection with packing models, phase diagrams with repulsive potentials, and the prediction of likely experimental superlattices.

Observation of 9-Fold Coordinated Amorphous TiO2 at High Pressure

Knowledge of the structure in amorphous dioxides is important in many fields of science and engineering. Here we report new experimental results of high-pressure polyamorphism in amorphous TiO₂ (a-TiO₂). Our data show that the Ti coordination number (CN) increases from 7.2 ± 0.3 at ∼16 GPa to 8.8 ± 0.3 at ∼70 GPa and finally reaches a plateau at 8.9 ± 0.3 at ≲86 GPa. The evolution of the structural changes under pressure is rationalized by the ratio (γ) of the ionic radius of Ti to that of O. It appears that the CN ≈ 9 plateau correlates with the two 9-fold coordinated polymorphs (cotunnite, Fe₂P) with different γ values. This CN-γ relationship is compared with those of a-SiO₂ and a-GeO₂, displaying remarkably consistent behavior between CN and γ. The unified CN-γ relationship may be generally used to predict the compression behavior of amorphous AO₂ compounds under extreme conditions.

Influence of Softness on the Stability of Binary Colloidal Crystals

Mixtures of two types of nanoparticles can self-assemble into a wide variety of binary colloidal crystals (also called binary nanoparticle superlattices), which are interesting for their structural diversity and potential applications. Although so-called packing models-which usually treat the particles as "hard" with only excluded volume interactions-seem to explain many reported dense crystalline phases, these models often fail to predict the right structure. Here, we examine the role of soft repulsive interparticle interactions on binary colloidal crystals comprising two sizes of spherical particles; such "softness" can arise due to ligand shells or screened electrostatics. We determine the ground state phase diagram of binary systems of particles interacting with an additive inverse power law potential using a basin hopping algorithm to calculate the enthalpy of an extremely large pool of candidate structures. We find that a surprisingly small amount of softness can destabilize dense packings in favor of less densely packed structures, which provides further evidence that considerations beyond packing are necessary for describing many of the observed phases of binary colloidal crystals. Importantly, we find that several of the phases stabilized by softness, which are characterized by relatively few interparticle contacts and a tendency for local icosahedral order, are more likely to be observed experimentally than those predicted by packing models. We also report a previously unknown dense AB₄ phase and conduct free energy calculations to examine how the stability of several crystals will vary with temperature. Our results further our understanding of why particular binary colloidal crystals form and will be useful as a reference for experimentalists working with softly repulsive colloids.

Using symmetry to elucidate the importance of stoichiometry in colloidal crystal assembly

We demonstrate a method based on symmetry to predict the structure of self-assembling, multicomponent colloidal mixtures. This method allows us to feasibly enumerate candidate structures from all symmetry groups and is many orders of magnitude more computationally efficient than combinatorial enumeration of these candidates. In turn, this permits us to compute ground-state phase diagrams for multicomponent systems. While tuning the interparticle potentials to produce potentially complex interactions represents the conventional route to designing exotic lattices, we use this scheme to demonstrate that simple potentials can also give rise to such structures which are thermodynamically stable at moderate to low temperatures. Furthermore, for a model two-dimensional colloidal system, we illustrate that lattices forming a complete set of 2-, 3-, 4-, and 6-fold rotational symmetries can be rationally designed from certain systems by tuning the mixture composition alone, demonstrating that stoichiometric control can be a tool as powerful as directly tuning the interparticle potentials themselves.

Statistics of small length scale density fluctuations in supercooled viscous liquids

Many successful theories of liquids near the melting temperature assume that small length scale density fluctuations follow Gaussian statistics. This paper presents a numerical investigation of density fluctuations in the supercooled viscous regime using an enhanced sampling method. Five model systems are investigated: the single component Lennard-Jones liquid, the Kob-Andersen binary mixture, the Wahnström binary mixture, the Lewis-Wahnström model of ortho-terphenyl, and the TIP4P/Ice model of water. The results show that the Gaussian approximation persists to a good degree into the supercooled viscous regime; however, it is less accurate at low temperatures. The analysis suggests that non-Gaussian fluctuations are related to crystalline configurations. Implications for theories of the glass transition are discussed.

Recent Advances in Multicomponent Particle Assembly

Particle assembly and co-assembly have been research frontiers in chemistry and material science in the past few decades. To achieve a large variety of intricate structures and functional materials, remarkable progress has been made in particle assembly principles and strategies. Essentially, particle assembly is driven by intrinsic interparticle interactions or the external control. In this article, we focus on binary or ternary particle co-assembly and review the principles and feasible strategies. These advances have led to new disciplines of microfabrication technology and material engineering. Although significant achievement on particle-based structures has been made, it is still challenging to fully develop general and facile strategies to precisely control the one-dimensional (1D) co-assembly. This article reviews the recent development on multicomponent particle co-assembly, which significantly increases structural complexity and functional diversity. In particular, we highlight the advances in the particle co-assembly of well-ordered 1D binary superstructures by liquid soft confinement. Finally, prospective outlook for future trends in this field is proposed.

Charging changes contact composition in binary sphere packings

Equal volume mixtures of small and large polytetrafluorethylene spheres are shaken in an atmosphere of controlled humidity which allows one to also control their tribocharging. We find that the contact numbers are charge dependent: As the charge density of the beads increases, the number of same-type contacts decreases and the number of opposite-type contacts increases. This change is not caused by a global segregation of the sample. Hence, tribocharging can be a way to tune the local composition of a granular material.

The stabilization of tubular crystals in mixtures of spherical particles

Novel crystal structures in binary atomic mixtures arise when the attractive well is wide enough to allow double occupancy by small particles. The resulting crystals consist of ordered packings of self assembled linear structures comprised of a cylindrical tube of large particles enclosing a close packed core of small particles that corresponds to a stacking of overlapping icosahedra. We show that the stability of these structures depends on two essential features of the spherically symmetric pairwise interactions: (i) a radius ratio between 0.414 and 0.588, and (ii) a width w of the attractive well in the interaction between unlike particles that satisfies w > σ_SS where σ_SS is the diameter of the small particle.

Self-Assembly of Colloidal Nanocrystals: From Intricate Structures to Functional Materials

Chemical methods developed over the past two decades enable preparation of colloidal nanocrystals with uniform size and shape. These Brownian objects readily order into superlattices. Recently, the range of accessible inorganic cores and tunable surface chemistries dramatically increased, expanding the set of nanocrystal arrangements experimentally attainable. In this review, we discuss efforts to create next-generation materials via bottom-up organization of nanocrystals with preprogrammed functionality and self-assembly instructions. This process is often driven by both interparticle interactions and the influence of the assembly environment. The introduction provides the reader with a practical overview of nanocrystal synthesis, self-assembly, and superlattice characterization. We then summarize the theory of nanocrystal interactions and examine fundamental principles governing nanocrystal self-assembly from hard and soft particle perspectives borrowed from the comparatively established fields of micrometer colloids and block copolymer assembly. We outline the extensive catalog of superlattices prepared to date using hydrocarbon-capped nanocrystals with spherical, polyhedral, rod, plate, and branched inorganic core shapes, as well as those obtained by mixing combinations thereof. We also provide an overview of structural defects in nanocrystal superlattices. We then explore the unique possibilities offered by leveraging nontraditional surface chemistries and assembly environments to control superlattice structure and produce nonbulk assemblies. We end with a discussion of the unique optical, magnetic, electronic, and catalytic properties of ordered nanocrystal superlattices, and the coming advances required to make use of this new class of solids.

Hard sphere packings within cylinders

Arrangements of identical hard spheres confined to a cylinder with hard walls have been used to model experimental systems, such as fullerenes in nanotubes and colloidal wire assembly. Finding the densest configurations, called close packings, of hard spheres of diameter σ in a cylinder of diameter D is a purely geometric problem that grows increasingly complex as D/σ increases, and little is thus known about the regime for D > 2.873σ. In this work, we extend the identification of close packings up to D = 4.00σ by adapting Torquato-Jiao's adaptive-shrinking-cell formulation and sequential-linear-programming (SLP) technique. We identify 17 new structures, almost all of them chiral. Beyond D ≈ 2.85σ, most of the structures consist of an outer shell and an inner core that compete for being close packed. In some cases, the shell adopts its own maximum density configuration, and the stacking of core spheres within it is quasiperiodic. In other cases, an interplay between the two components is observed, which may result in simple periodic structures. In yet other cases, the very distinction between the core and shell vanishes, resulting in more exotic packing geometries, including some that are three-dimensional extensions of structures obtained from packing hard disks in a circle.

A Geometric-Structure Theory for Maximally Random Jammed Packings

Maximally random jammed (MRJ) particle packings can be viewed as prototypical glasses in that they are maximally disordered while simultaneously being mechanically rigid. The prediction of the MRJ packing density ϕMRJ, among other packing properties of frictionless particles, still poses many theoretical challenges, even for congruent spheres or disks. Using the geometric-structure approach, we derive for the first time a highly accurate formula for MRJ densities for a very wide class of two-dimensional frictionless packings, namely, binary convex superdisks, with shapes that continuously interpolate between circles and squares. By incorporating specific attributes of MRJ states and a novel organizing principle, our formula yields predictions of ϕMRJ that are in excellent agreement with corresponding computer-simulation estimates in almost the entire α-x plane with semi-axis ratio α and small-particle relative number concentration x. Importantly, in the monodisperse circle limit, the predicted ϕMRJ = 0.834 agrees very well with the very recently numerically discovered MRJ density of 0.827, which distinguishes it from high-density "random-close packing" polycrystalline states and hence provides a stringent test on the theory. Similarly, for non-circular monodisperse superdisks, we predict MRJ states with densities that are appreciably smaller than is conventionally thought to be achievable by standard packing protocols.

Binary nanoparticle superlattices of soft-particle systems

The solid-phase diagram of binary systems consisting of particles of diameter σA = σ and σB = γσ (γ ≤ 1) interacting with an inverse p = 12 power law is investigated as a paradigm of a soft potential. In addition to the diameter ratio γ that characterizes hard-sphere models, the phase diagram is a function of an additional parameter that controls the relative interaction strength between the different particle types. Phase diagrams are determined from extremes of thermodynamic functions by considering 15 candidate lattices. In general, it is shown that the phase diagram of a soft repulsive potential leads to the morphological diversity observed in experiments with binary nanoparticles, thus providing a general framework to understand their phase diagrams. Particular emphasis is given to the two most successful crystallization strategies so far: evaporation of solvent from nanoparticles with grafted hydrocarbon ligands and DNA programmable self-assembly.

Favoured local structures in liquids and solids: a 3D lattice model

We investigate the connection between the geometry of Favoured Local Structures (FLS) in liquids and the associated liquid and solid properties. We introduce a lattice spin model - the FLS model on a face-centered cubic lattice - where this geometry can be arbitrarily chosen among a discrete set of 115 possible FLS. We find crystalline groundstates for all choices of a single FLS. Sampling all possible FLS's, we identify the following trends: (i) low symmetry FLS's produce larger crystal unit cells but not necessarily higher energy groundstates, (ii) chiral FLS's exhibit peculiarly poor packing properties, (iii) accumulation of FLS's in supercooled liquids is linked to large crystal unit cells, and (iv) low symmetry FLS's tend to find metastable structures on cooling.

On the phase behavior of binary mixtures of nanoparticles

The assembly of mixtures of nanoparticles with different properties into a binary nanoparticle superlattice (BNSL) provides a route to fabricate novel classes of materials with properties emerging from the choice of the building blocks. The common theoretical approach based on the hard-spheres model predicts crystallization of only a few metastable binary superstructures (NaCl, AlB₂ or the AB₁₃). Recently [Shevchenko, E. V.; Talapin, D. V.; O'Brien, S.; Murray, C. B. Nature 2006; 439, 55.)], it has been demonstrated that with the use of a combination of semiconducting, metallic, and magnetic nanoparticles, a variety of novel BNSL structures were formed, where at least 10 were low density structures that have not been previously reported. While some of the structures can be explained by the addition of electrostatic interactions, it is clear that at the nanometer scale one needs to consider other influences, such as van der Waals forces, steric effects, etc. Motivated by those experiments, we study, using Monte Carlo simulations, the phase behavior of binary mixtures of nanoparticles interacting via a combination of hard-core electrostatics and van der Waals forces. We include a tuning parameter that can be used to balance between electrostatic and dispersion interactions and study the phase behavior as a function of the different charges and size ratios of the nanoparticles. The results indicate that at the nanoscale, both electrostatic and dispersion interactions are necessary to explain the experimental observed BNSL structures.

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.

Geometrical frustration and static correlations in a simple glass former

We study the geometrical frustration scenario of glass formation for simple hard-sphere models. We find that the dual picture in terms of defects brings little insight and no theoretical simplification for the understanding of the slowing down of relaxation, because of the strong frustration characterizing these systems. The possibility of a growing static length is furthermore found to be physically irrelevant in the regime that is accessible to computer simulations.

Phase diagram and structural diversity of the densest binary sphere packings

The densest binary sphere packings have historically been very difficult to determine. The only rigorously known packings in the α-x plane of sphere radius ratio α and relative concentration x are at the Kepler limit α=1, where packings are monodisperse. Utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [S. Torquato and Y. Jiao, Phys. Rev. E 82, 061302 (2010)], we present the most comprehensive determination to date of the phase diagram in (α,x) for the densest binary sphere packings. Unexpectedly, we find many distinct new densest packings.

Structural searches using isopointal sets as generators: densest packings for binary hard sphere mixtures

Algorithms to search for crystal structures that optimize some extensive property (energy, volume, etc) typically make use of random particle reorganizations in the context of one or more numerical techniques such as simulated annealing, genetic algorithms or biased random walks, applied to the coordinates of every particle in the unit cell, together with the cell angles and lengths. In this paper we describe the restriction of such searches to predefined isopointal sets, breaking the problem into countable sub-problems which exploit crystal symmetries to reduce the dimensionality of the search space. Applying this method to the search for maximally packed mixtures of hard spheres of two sizes, we demonstrate that the densest packed structures can be identified by searches within a couple of isopointal sets. For the A(2)B system, the densest known packings over the entire tested range 0.2 < r(A)/r(B) < 2.5, including some improvements on previous optima, can all be identified by searches within a single isopointal set. In the case of the AB composition, searches of two isopointal sets generate the densest packed structures over the radius ratio range 0.2 < r(A)/r(B) < 5.0.

Efficient method for predicting crystal structures at finite temperature: variable box shape simulations

We present an efficient and robust method based on Monte Carlo simulations for predicting crystal structures at finite temperature. We apply this method, which is surprisingly easy to implement, to a variety of systems, demonstrating its effectiveness for hard, attractive, and anisotropic interactions, binary mixtures, semi-long-range soft interactions, and truly long-range interactions where the truly long-range interactions are treated using Ewald sums. In the case of binary hard-sphere mixtures, star polymers, and binary Lennard-Jones mixtures, the crystal structures predicted by this algorithm are consistent with literature, providing confidence in the method. Finally, we predict new crystal structures for hard asymmetric dumbbell particles, bowl-like particles and hard oblate cylinders and present the phase diagram for the oblate cylinders based on full free energy calculations.

Prediction of binary hard-sphere crystal structures

We present a method based on a combination of a genetic algorithm and Monte Carlo simulations to predict close-packed crystal structures in hard-core systems. We employ this method to predict the binary crystal structures in a mixture of large and small hard spheres with various stoichiometries and diameter ratios between 0.4 and 0.84. In addition to known binary hard-sphere crystal structures similar to NaCl and AlB2, we predict additional crystal structures with the symmetry of CrB, gammaCuTi, alphaIrV, HgBr2, AuTe2, Ag2Se, and various structures for which an atomic analog was not found. In order to determine the crystal structures at infinite pressures, we calculate the maximum packing density as a function of size ratio for the crystal structures predicted by our GA using a simulated annealing approach.