Densest binary sphere packings

Structures of LaH10, EuH9, and UH8 superhydrides rationalized by electron counting and Jahn-Teller distortions in a covalent cluster model

The superconducting hydrides LaH₁₀, EuH₉ and UH₈ are studied using chemically intuitive bonding analysis of periodic and molecular models. We find trends in the crystallographic and electronic structures of the materials by focusing on chemically meaningful building blocks in the predicted H sublattices. Atomic charge calculations, using two complementary techniques, allow us to assign oxidation states to the metals and divide the H sublattice into neutral and anionic components. Cubic [H₈]^q- clusters are an important structural motif, and molecular orbital analysis of this cluster in isolation shows the crystal structures to be consistent with our oxidation state assignments. Crystal orbital Hamilton population analysis confirms the applicability of the cluster model to the periodic electronic structure. A Jahn-Teller distortion predicted by MO analysis rationalises the distortion observed in a prior study of EuH₉. The impact of this distortion on superconductivity is determined, and implications for crystal structure prediction in other metal-hydrogen systems are discussed. Additionally, the performance of electronic structure analysis methods at high pressures are tested and recommendations for future studies are given. These results demonstrate the value of simple bonding models in rationalizing chemical structures under extreme conditions.

Ligand Effects in Assembly of Cubic and Spherical Nanocrystals: Applications to Packing of Perovskite Nanocubes

We establish the formula representing cubic nanocrystals (NCs) as hard cubes taking into account the role of the ligands and describe how these results generalize to any other NC shapes. We derive the conditions under which the hard cube representation breaks down and provide explicit expressions for the effective size. We verify the results from the detailed potential of mean force calculations for two nanocubes in different orientations as well as with spherical nanocrystals. Our results explicitly demonstrate the relevance of certain ligand conformations, i.e., "vortices", and show that edges and corners provide natural sites for their emergence. We also provide both simulations and experimental results with single component cubic perovskite nanocrystals assembled into simple cubic superlattices, which further corroborate theoretical predictions. In this way, we extend the Orbifold Topological Model (OTM) accounting for the role of ligands beyond spherical nanocrystals and discuss its extension to arbitrary nanocrystal shapes. Our results provide detailed predictions for recent superlattices of perovskite nanocubes and spherical nanocrystals. Problems with existing united atom force fields are discussed.

Discovery of two-dimensional binary nanoparticle superlattices using global Monte Carlo optimization

Binary nanoparticle (NP) superlattices exhibit distinct collective plasmonic, magnetic, optical, and electronic properties. Here, we computationally demonstrate how fluid-fluid interfaces could be used to self-assemble binary systems of NPs into 2D superlattices when the NP species exhibit different miscibility with the fluids forming the interface. We develop a basin-hopping Monte Carlo (BHMC) algorithm tailored for interface-trapped structures to rapidly determine the ground-state configuration of NPs, allowing us to explore the repertoire of binary NP architectures formed at the interface. By varying the NP size ratio, interparticle interaction strength, and difference in NP miscibility with the two fluids, we demonstrate the assembly of an array of exquisite 2D periodic architectures, including AB-, AB₂-, and AB₃-type monolayer superlattices as well as AB-, AB₂-, A₃B₅-, and A₄B₆-type bilayer superlattices. Our results suggest that the interfacial assembly approach could be a versatile platform for fabricating 2D colloidal superlattices with tunable structure and properties.

Hard-sphere jamming through the lens of linear optimization

The jamming transition is ubiquitous. It is present in granular matter, foams, colloids, structural glasses, and many other systems. Yet, it defines a critical point whose properties still need to be fully understood. Recently, a major breakthrough came about when the replica formalism was extended to build a mean-field theory that provides an exact description of the jamming transition of spherical particles in the infinite-dimensional limit. While such theory explains the jamming critical behavior of both soft and hard spheres, investigating the transition in finite-dimensional systems poses very difficult and different problems, in particular from the numerical point of view. Soft particles are modeled by continuous potentials; thus, their jamming point can be reached through efficient energy minimization algorithms. In contrast, the latter methods are inapplicable to hard-sphere (HS) systems since the interaction energy among the particles is always zero by construction. To overcome these difficulties, here we recast the jamming of hard spheres as a constrained optimization problem and introduce the CALiPPSO algorithm, capable of readily producing jammed HS packings without including any effective potential. This algorithm brings a HS configuration of arbitrary dimensions to its jamming point by solving a chain of linear optimization problems. We show that there is a strict correspondence between the force balance conditions of jammed packings and the properties of the optimal solutions of CALiPPSO, whence we prove analytically that our packings are always isostatic and in mechanical equilibrium. Furthermore, using extensive numerical simulations, we show that our algorithm is able to probe the complex structure of the free-energy landscape, finding qualitative agreement with mean-field predictions. We also characterize the algorithmic complexity of CALiPPSO and provide an open-source implementation of it.

Self-Assembly of Colloidal Nanocrystals into 3D Binary Mesocrystals

ConspectusBiominerals are unique materials found in many living organisms that often display outstanding functionalities attributed to their mesocrystalline structure. Mesocrystals are nanocrystal superstructures with mutual crystallographic alignment of the building units. One could thus imagine these optimized evolutionary systems as archetypes to fabricate advanced materials. The main advantage of such systems relies on their ability to combine the features of the nanocrystals with those of single crystalline microscopic structures, yielding assemblies with directional, enhanced, and potentially emergent properties. Moreover, fueled by the promises of multifunctional materials with unprecedented and tunable properties, the rational design of mesocrystals assembled from two distinct colloidal nanocrystal ensembles has become a recent focus of research. However, the combination of dissimilar nanocrystals into ordered binary superstructures is still a major scientific challenge due to the nature of the coassembly process.We focus this Account on the growth of tridimensional (3D) binary mesocrystals and the understanding of the self-assembly of two colloidal nanocrystal ensembles with the ultimate goal to serve as a basis for more rational mesocrystal syntheses in the future. The formation of mesocrystals demands nanocrystals with defined surface faceting, the primary factor influencing their oriented self-assembly. Notably, such a process cannot be successfully afforded without functionalized nanocrystals with high and, in many cases, tunable colloidal stability. Besides, the nature and solvation degree of the surface ligand shell influences the effective shape of the nanocrystals and the kinetics of self-assembly. If the assembly is triggered by reducing the colloidal stability with nonsolvents, 3D single-component mesocrystals are often grown. Here, the different magnitude of the van der Waals attraction forces between nanocrystals with differing compositions, dimensions, and morphologies generally favors the segregation and growth of single component mesocrystals. This phenomenon was illustrated during the successful preparation of 3D binary mesocrystals composed of iron oxide and platinum nanocubes. Although the building blocks possessed comparable sizes and were stabilized by similar ligands, the amount of the second component could only be arbitrarily tuned up to some extent, even when the assembly conditions were rationally optimized to achieve 3D binary mesocrystals. Only a small amount of it was effectively incorporated into the matrix of the initial mesocrystal. The 3D binary mesocrystal growth process demands a delicate control over the size, shape, and surface chemistry of the nanocrystals, the solvent nature, and the self-assembly process. Hence, the improvement of our ability to control the synthesis of 3D binary mesocrystalline materials is critical to exploit their potential toward technological applications in catalysis, energy storage, or structural materials.

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.

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.

Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids

Computer simulations of the fluid-to-solid phase transition in the hard sphere system were instrumental for our understanding of crystallization processes. But while colloid experiments and theory have been predicting the stability of several binary hard sphere crystals for many years, simulations were not successful to confirm this phenomenon. Here, we report the growth of binary hard sphere crystals isostructural to Laves phases, AlB_{2}, and NaZn_{13} in simulation directly from the fluid. We analyze particle kinetics during Laves phase growth using event-driven molecular dynamics simulations with and without swap moves that speed up diffusion. The crystallization process transitions from nucleation and growth to spinodal decomposition already deep within the fluid-solid coexistence regime. Finally, we present packing fraction-size ratio state diagrams in the vicinity of the stability regions of three binary crystals.

Crystal structures in binary hard-sphere colloid-droplet mixtures with patchy cross interactions

A binary mixture of droplets and patchy colloids, where patches are arranged in tetrahedral symmetry, is studied with Metropolis Monte Carlo simulations. The colloidal patches attract droplets, while both the colloid-colloid and the droplet-droplet interactions are hard sphere like. We find stable crystal structures with atomic analogs ZnS, CaF_{2}, and fcc or hcp (face centered cubic or hexagonal close packed) of the droplets coexisting with a dispersed fluid of the colloids. The simulated crystal structures agree well with those predicted by close-packing calculations for an intermediate range of droplet-colloid size ratios. A discrepancy between the simulations and theoretical predictions occurs at low and high size ratios. The results of the simulations for mixtures with anisotropic colloid-droplet interactions reveal a richer phase diagram, with ZnS-gas and ZnS-fluid coexistence, as compared to the isotropic case. For the example of a square planar patch arrangement, we find a particular crystal structure, consisting of two interpenetrating fcc or hcp lattices with right bond angles. Such a structure has no known atomic analog. Our study of generic models of anisotropic colloid-droplet mixtures could provide a promising way towards the fabrication of novel and complex colloidal structures.

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.

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.

Population annealing simulations of a binary hard-sphere mixture

Population annealing is a sequential Monte Carlo scheme well suited to simulating equilibrium states of systems with rough free energy landscapes. Here we use population annealing to study a binary mixture of hard spheres. Population annealing is a parallel version of simulated annealing with an extra resampling step that ensures that a population of replicas of the system represents the equilibrium ensemble at every packing fraction in an annealing schedule. The algorithm and its equilibration properties are described, and results are presented for a glass-forming fluid composed of a 50/50 mixture of hard spheres with diameter ratio of 1.4:1. For this system, we obtain precise results for the equation of state in the glassy regime up to packing fractions φ≈0.60 and study deviations from the Boublik-Mansoori-Carnahan-Starling-Leland equation of state. For higher packing fractions, the algorithm falls out of equilibrium and a free volume fit predicts jamming at packing fraction φ≈0.667. We conclude that population annealing is an effective tool for studying equilibrium glassy fluids and the jamming transition.

Highly Elastic and Self-Healing Composite Colloidal Gels

Composite colloidal gels are formed by the pH-induced electrostatic assembly of silica and gelatin nanoparticles. These injectable and moldable colloidal gels are able to withstand substantial compressive and tensile loads, and exhibit a remarkable self-healing efficiency. This study provides new, critical insight into the structural and mechanical properties of composite colloidal gels and opens up new avenues for practical application of colloidal gels.

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.

Quasicrystalline nanocrystal superlattice with partial matching rules

Expanding the library of self-assembled superstructures provides insight into the behaviour of atomic crystals and supports the development of materials with mesoscale order. Here we build on recent findings of soft matter quasicrystals and report a quasicrystalline binary nanocrystal superlattice that exhibits correlations in the form of partial matching rules reducing tiling disorder. We determine a three-dimensional structure model through electron tomography and direct imaging of surface topography. The 12-fold rotational symmetry of the quasicrystal is broken in sublayers, forming a random tiling of rectangles, large triangles and small triangles with 6-fold symmetry. We analyse the geometry of the experimental tiling and discuss factors relevant for the stabilization of the quasicrystal. Our joint experimental-computational study demonstrates the power of nanocrystal superlattice engineering and further narrows the gap between the richness of crystal structures found with atoms and in soft matter assemblies.

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.

Confined disordered strictly jammed binary sphere packings

Disordered jammed packings under confinement have received considerably less attention than their bulk counterparts and yet arise in a variety of practical situations. In this work, we study binary sphere packings that are confined between two parallel hard planes and generalize the Torquato-Jiao (TJ) sequential linear programming algorithm [Phys. Rev. E 82, 061302 (2010)] to obtain putative maximally random jammed (MRJ) packings that are exactly isostatic with high fidelity over a large range of plane separation distances H, small to large sphere radius ratio α, and small sphere relative concentration x. We find that packing characteristics can be substantially different from their bulk analogs, which is due to what we term "confinement frustration." Rattlers in confined packings are generally more prevalent than those in their bulk counterparts. We observe that packing fraction, rattler fraction, and degree of disorder of MRJ packings generally increase with H, though exceptions exist. Discontinuities in the packing characteristics as H varies in the vicinity of certain values of H are due to associated discontinuous transitions between different jammed states. When the plane separation distance is on the order of two large-sphere diameters or less, the packings exhibit salient two-dimensional features; when the plane separation distance exceeds about 30 large-sphere diameters, the packings approach three-dimensional bulk packings. As the size contrast increases (as α decreases), the rattler fraction dramatically increases due to what we call "size-disparity" frustration. We find that at intermediate α and when x is about 0.5 (50-50 mixture), the disorder of packings is maximized, as measured by an order metric ψ that is based on the number density fluctuations in the direction perpendicular to the hard walls. We also apply the local volume-fraction variance σ(τ)(2)(R) to characterize confined packings and find that these packings possess essentially the same level of hyperuniformity as their bulk counterparts. Our findings are generally relevant to confined packings that arise in biology (e.g., structural color in birds and insects) and may have implications for the creation of high-density powders and improved battery designs.

Mixing effects in the crystallization of supercooled quantum binary liquids

By means of Raman spectroscopy of liquid microjets, we have investigated the crystallization process of supercooled quantum liquid mixtures composed of parahydrogen (pH2) or orthodeuterium (oD2) diluted with small amounts of neon. We show that the introduction of the Ne impurities affects the crystallization kinetics in terms of a significant reduction of the measured pH2 and oD2 crystal growth rates, similarly to what found in our previous work on supercooled pH2-oD2 liquid mixtures [Kühnel et al., Phys. Rev. B 89, 180201(R) (2014)]. Our experimental results, in combination with path-integral simulations of the supercooled liquid mixtures, suggest in particular a correlation between the measured growth rates and the ratio of the effective particle sizes originating from quantum delocalization effects. We further show that the crystalline structure of the mixtures is also affected to a large extent by the presence of the Ne impurities, which likely initiate the freezing process through the formation of Ne-rich crystallites.

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.

Generation of random microstructures and prediction of sound velocity and absorption for open foams with spherical pores

This paper proposes and discusses an approach for the design and quality inspection of the morphology dedicated for sound absorbing foams, using a relatively simple technique for a random generation of periodic microstructures representative for open-cell foams with spherical pores. The design is controlled by a few parameters, namely, the total open porosity and the average pore size, as well as the standard deviation of pore size. These design parameters are set up exactly and independently, however, the setting of the standard deviation of pore sizes requires some number of pores in the representative volume element (RVE); this number is a procedure parameter. Another pore structure parameter which may be indirectly affected is the average size of windows linking the pores, however, it is in fact weakly controlled by the maximal pore-penetration factor, and moreover, it depends on the porosity and pore size. The proposed methodology for testing microstructure-designs of sound absorbing porous media applies the multi-scale modeling where some important transport parameters-responsible for sound propagation in a porous medium-are calculated from microstructure using the generated RVE, in order to estimate the sound velocity and absorption of such a designed material.

Attractive particle interaction forces and packing density of fine glass powders

We study the packing of fine glass powders of mean particle diameter in the range (4-52) μm both experimentally and by numerical DEM simulations. We obtain quantitative agreement between the experimental and numerical results, if both types of attractive forces of particle interaction, adhesion and non-bonded van der Waals forces are taken into account. Our results suggest that considering only viscoelastic and adhesive forces in DEM simulations may lead to incorrect numerical predictions of the behavior of fine powders. Based on the results from simulations and experiments, we propose a mathematical expression to estimate the packing fraction of fine polydisperse powders as a function of the average particle size.

Simple effective rule to estimate the jamming packing fraction of polydisperse hard spheres

A recent proposal in which the equation of state of a polydisperse hard-sphere mixture is mapped onto that of the one-component fluid is extrapolated beyond the freezing point to estimate the jamming packing fraction ϕJ of the polydisperse system as a simple function of M1M3/M22, where Mk is the kth moment of the size distribution. An analysis of experimental and simulation data of ϕJ for a large number of different mixtures shows a remarkable general agreement with the theoretical estimate. To give extra support to the procedure, simulation data for seventeen mixtures in the high-density region are used to infer the equation of state of the pure hard-sphere system in the metastable region. An excellent collapse of the inferred curves up to the glass transition and a significant narrowing of the different out-of-equilibrium glass branches all the way to jamming are observed. Thus, the present approach provides an extremely simple criterion to unify in a common framework and to give coherence to data coming from very different polydisperse hard-sphere mixtures.

Dense periodic packings of tori

Dense packings of nonoverlapping bodies in three-dimensional Euclidean space R(3) are useful models of the structure of a variety of many-particle systems that arise in the physical and biological sciences. Here we investigate the packing behavior of congruent ring tori in R(3), which are multiply connected nonconvex bodies of genus 1, as well as horn and spindle tori. Specifically, we analytically construct a family of dense periodic packings of unlinked tori guided by the organizing principles originally devised for simply connected solid bodies [Torquato and Jiao, Phys. Rev. E 86, 011102 (2012)]. We find that the horn tori as well as certain spindle and ring tori can achieve a packing density not only higher than that of spheres (i.e., π/sqrt[18] = 0.7404...) but also higher than the densest known ellipsoid packings (i.e., 0.7707...). In addition, we study dense packings of clusters of pair-linked ring tori (i.e., Hopf links), which can possess much higher densities than corresponding packings consisting of unlinked tori.

Shape alloys of nanorods and nanospheres from self-assembly

Mixtures of anisotropic nanocrystals promise a great diversity of superlattices and phase behaviors beyond those of single-component systems. However, obtaining a colloidal shape alloy in which two different shapes are thermodynamically coassembled into a crystalline superlattice has remained a challenge. Here we present a joint experimental-computational investigation of two geometrically ubiquitous nanocrystalline building blocks-nanorods and nanospheres-that overcome their natural entropic tendency toward macroscopic phase separation and coassemble into three intriguing phases over centimeter scales, including an AB2-type binary superlattice. Monte Carlo simulations reveal that, although this shape alloy is entropically stable at high packing fraction, demixing is favored at experimental densities. Simulations with short-ranged attractive interactions demonstrate that the alloy is stabilized by interactions induced by ligand stabilizers and/or depletion effects. An asymmetry in the relative interaction strength between rods and spheres improves the robustness of the self-assembly process.

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.

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.

Binary superlattices from colloidal nanocrystals and giant polyoxometalate clusters

We report a new kind of long-range ordered binary superlattices comprising atomically defined inorganic clusters and colloidally synthesized nanocrystals. In a model system, we combined surfactant-encapsulated, nearly spherical giant polyoxometalate clusters containing 2.9 nm polyoxomolybdate or 2.5 nm polyoxovanadomolybdate cores with monodisperse colloidal semiconductor nanocrystals (PbS, CdSe, PbS/CdS; 4-11 nm). The results are rationalized on the basis of dense packing principles of sterically stabilized particles with predominantly hard-spherelike interparticle interactions. By varying the size-ratios and relative concentrations of constituents, we obtained known thermodynamically stable binary packings of hard-spheres such as NaCl, AlB2, and NaZn13 lattices and also CaCu5-type lattice and aperiodic quasicrystals with 12-fold symmetry. These results suggest that other kinds of cluster materials such as fullerenes and magic-sized metallic and semiconductor clusters can also be integrated into supramolecular assemblies with nanocrystals. Furthermore, synergistic effects are expected from the combination of redox and catalytic properties of polyoxometalates with excitonic and plasmonic properties of inorganic nanocrystals.

Electron tomography resolves a novel crystal structure in a binary nanocrystal superlattice

The self-assembly of different nanocrystals into a binary superlattice is of interest for both colloidal science and nanomaterials science. New properties may emerge from the interaction between the nanocrystal building blocks that are ordered in close contact in three dimensions. Identification of the superlattice structure including its defects is of key interest in understanding the electrical and optical properties of these systems. Transmission electron microscopy (TEM) has been very instrumental to reach this goal but fails for complex crystal structures and buried defects. Here, we use electron tomography to resolve the three-dimensional crystal structure of a binary superlattice that could not be resolved by TEM only. The structure with a [PbSe]6[CdSe]19 stoichiometry has no analogue in the atomic world. Moreover we will show how tomography can overcome the clouding effects of planar defects on structure identification by TEM.

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.

[N]pT ensemble and finite-size-scaling study of the critical isostructural transition in the generalized exponential model of index 4

First-order transitions of system where both lattice site occupancy and lattice spacing fluctuate, such as cluster crystals, cannot be efficiently studied by traditional simulation methods, which necessarily fix one of these two degrees of freedom. The difficulty, however, can be surmounted by the generalized [N]pT ensemble [J. Chem. Phys. 136, 214106 (2012)]. Here we show that histogram reweighting and the [N]pT ensemble can be used to study an isostructural transition between cluster crystals of different occupancy in the generalized exponential model of index 4 (GEM-4). Extending this scheme to finite-size scaling studies also allows us to accurately determine the critical point parameters and to verify that it belongs to the Ising universality class.