A cell theory for solid solutions: Application to hard sphere mixtures

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.

High antisite defect concentrations in hard-sphere colloidal Laves phases

Binary mixtures of hard spheres can spontaneously self-assemble into binary crystals. Computer simulations have been especially useful in mapping out the phase behaviour of these mixtures, under the assumption that the stoichiometry of the binary crystal is ideal. Here we show that for a size ratio of q = 0.82 this assumption is not valid near the coexistence region between the fluid and the stable binary crystal, the MgZn2 Laves phase. Instead we find a surprisingly high number of antisite defects: up to 2% of the large spheres are replaced by small spheres in equilibrium. We demonstrate that the defect concentration can be estimated using simple approximations, providing an easy way to identify systems where antisite defects play an important role. Our results shed new light on the self-assembly of colloidal Laves phases, and demonstrate the importance of antisite defects in binary crystals.

Low-Temperature Crystal Structures of the Hard Core Square Shoulder Model

In many cases, the stability of complex structures in colloidal systems is enhanced by a competition between different length scales. Inspired by recent experiments on nanoparticles coated with polymers, we use Monte Carlo simulations to explore the types of crystal structures that can form in a simple hard-core square shoulder model that explicitly incorporates two favored distances between the particles. To this end, we combine Monte Carlo-based crystal structure finding algorithms with free energies obtained using a mean-field cell theory approach, and draw phase diagrams for two different values of the square shoulder width as a function of the density and temperature. Moreover, we map out the zero-temperature phase diagram for a broad range of shoulder widths. Our results show the stability of a rich variety of crystal phases, such as body-centered orthogonal (BCO) lattices not previously considered for the square shoulder model.

Hierarchically self-assembled hexagonal honeycomb and kagome superlattices of binary 1D colloids

Synthesis of binary nanoparticle superlattices has attracted attention for a broad spectrum of potential applications. However, this has remained challenging for one-dimensional nanoparticle systems. In this study, we investigate the packing behavior of one-dimensional nanoparticles of different diameters into a hexagonally packed cylindrical micellar system and demonstrate that binary one-dimensional nanoparticle superlattices of two different symmetries can be obtained by tuning particle diameter and mixing ratios. The hexagonal arrays of one-dimensional nanoparticles are embedded in the honeycomb lattices (for AB₂ type) or kagome lattices (for AB₃ type) of micellar cylinders. The maximization of free volume entropy is considered as the main driving force for the formation of superlattices, which is well supported by our theoretical free energy calculations. Our approach provides a route for fabricating binary one-dimensional nanoparticle superlattices and may be applicable for inorganic one-dimensional nanoparticle systems.

Binary mixtures of 1D particles are rarely observed to cooperatively self-assemble into binary superlattices, as the particle types separate into phases. Here, the authors design a system that avoids phase separation, obtaining binary superlattices with different symmetries by simply tuning the particle diameter and mixture composition.

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.

Phase behaviour of deionized binary mixtures of charged colloidal spheres

We review recent work on the phase behaviour of binary charged sphere mixtures as a function of particle concentration and composition. Both size ratios Γ and charge ratios Λ are varied over a wide range. Unlike the case for hard spheres, the long-ranged Coulomb interaction stabilizes the crystal phase at low particle concentrations and shifts the occurrence of amorphous solids to particle concentrations considerably larger than the freezing concentration. Depending on Γ and Λ, we observe upper azeotrope, spindle, lower azeotrope and eutectic types of phase diagrams, all known well from metal systems. Most solids are of body centred cubic structure. Occasionally stoichiometric compounds are formed at large particle concentrations. For very low Γ, entropic effects dominate and induce a fluid-fluid phase separation. Since for charged spheres the charge ratio Λ is also decisive for the type of phase diagram, future experiments with charge variable silica spheres are suggested.

Crystal nucleation in binary hard sphere mixtures: A Monte Carlo simulation study

We present calculations of the nucleation barrier during crystallization in binary hard sphere mixtures under moderate degrees of supercooling using Monte Carlo simulations in the isothermal-isobaric semigrand ensemble in conjunction with an umbrella sampling technique. We study both additive and negatively nonadditive binary hard sphere systems. The solid-fluid phase diagrams of such systems show a rich variety of behavior, ranging from simple spindle shapes to the appearance of azeotropes and eutectics to the appearance of substitutionally ordered solid phase compounds. We investigate the effect of these types of phase behavior upon the nucleation barrier and the structure of the critical nucleus. We find that the underlying phase diagram has a significant effect on the mechanism of crystal nucleation. Our calculations indicate that fractionation of the species upon crystallization increases the difficulty of crystallization of fluid mixtures and in the absence of fractionation (azeotropic conditions) the nucleation barrier is comparable to pure fluids. We also calculate the barrier to nucleation of a substitutionally ordered compound solid. In such systems, which also show solid-solid phase separation, we find that the phase that nucleates is the one whose equilibrium composition is closer to the composition of the fluid phase.

Desorption behavior of quench-condensed argon-neon mixtures

The desorption behavior of quench-condensed rare gas films has been investigated using high frequency surface acoustic waves. Measurements of pure films of argon and neon and of the binary mixture Ar(1-c)/Ne(c) have been carried out. For small and very large neon concentration c(Ne) a behavior is found which indicates the existence of a substitutionally disordered solid. In contrast, in the wide range of concentration 0.25< c(Ne)< 0.92 two discrete temperatures for neon desorption exist. The data clearly indicate the occurrence of two separate phases, one of pure neon, the other of crystallites with either Ar(2)Ne or Ar(3)Ne structure.

Order-disorder transition in the solid phase of a charged hard sphere model

We investigate the solid phases of the restricted primitive model (RPM). Monte Carlo simulations show the existence of an order-disorder transition from a substitutionally disordered face centered cubic lattice (fcc) to a new ordered fcc structure which is proposed as the ground state of the RPM at the close packing density. Our results suggest that the new phase might turn out in a new triple point in the RPM phase diagram involving three solid phases: CsCl, fcc ordered and fcc disordered structures. The order-disorder transition is also studied using the cell theory. The theory shows good agreement with the simulation results and suggests that the transition is weakly first order.