Nucleus-size pinning for determination of nucleation free-energy barriers and nucleus geometry

Effect of particle anisotropy on the thermodynamics and kinetics of ordering transitions in hard faceted particles

Monte Carlo simulations were used to study the influence of particle aspect ratio on the kinetics and phase behavior of hard gyrobifastigia (GBF). First, the formation of a highly anisotropic nucleus shape in the isotropic-to-crystal transition in regular GBF is explained by the differences in interfacial free energies of various crystal planes and the nucleus geometry predicted by the Wulff construction. GBF-related shapes with various aspect ratios were then studied, mapping their equations of state, determining phase coexistence conditions via interfacial pinning, and computing nucleation free-energy barriers via umbrella sampling using suitable order parameters. Our simulations reveal a reduction of the kinetic barrier for isotropic-crystal transition upon an increase in aspect ratio, and that for highly oblate and prolate aspect ratios, an intermediate nematic phase is stabilized. Our results and observations also support two conjectures for the formation of the crystalline state from the isotropic phase: that low phase free energies at the ordering phase transition correlate with low transition barriers and that the emergence of a mesophase provides a steppingstone that expedites crystallization.

Comparing four hard-sphere approximations for the low-temperature WCA melting line

By combining interface-pinning simulations with numerical integration of the Clausius–Clapeyron equation, we accurately determine the melting-line coexistence pressure and fluid/crystal densities of the Weeks–Chandler–Andersen system, covering four decades of temperature. The data are used for comparing the melting-line predictions of the Boltzmann, Andersen–Weeks–Chandler, Barker–Henderson, and Stillinger hard-sphere approximations. The Andersen–Weeks–Chandler and Barker–Henderson theories give the most accurate predictions, and they both work excellently in the zero-temperature limit for which analytical expressions are derived here.

A route to hierarchical assembly of colloidal diamond

Photonic crystals, appealing for their ability to control light, are constructed from periodic regions of different dielectric constants. Yet, the structural holy grail in photonic materials, diamond, remains challenging to synthesize at the colloidal length scale. Here we explore new ways to assemble diamond using modified gyrobifastigial (mGBF) nanoparticles, a shape that resembles two anti-aligned triangular prisms. We investigate the parameter space that leads to the self-assembly of diamond, and we compare the likelihood of defects in diamond self-assembled via mGBF vs. the nanoparticle shape that is the current focus for assembling diamond, the truncated tetrahedra. We introduce a potential route for realizing mGBF particles by dimerizing triangular prisms using attractive patches, and we report the impact of this superstructure on the photonic properties.

Facile self-assembly of colloidal diamond from tetrahedral patchy particles via ring selection

The self-assembly of colloidal diamond–a classic example of an open crystal with the low coordination number of four and much sought after due to its applications in visible light management–from designer spherical colloidal particles has proved challenging over the years. The formation of the diamond lattice from tetrahedral patchy particles is hampered by the propensity to form competing open periodic structures for narrow patches or dynamically arrested states for wider patches, leaving a narrow window in design space where diamond crystals may be realized. Our two-component system of designer tetrahedral patchy particles supports a significantly wider range for patch sizes for programmed self-assembly, thus facilitating experimental fabrication, and offers fundamental insight into crystallization into open lattices.

Diamond-structured crystals, particularly those with cubic symmetry, have long been attractive targets for the programmed self-assembly of colloidal particles, due to their applications as photonic crystals that can control the flow of visible light. While spherical particles decorated with four patches in a tetrahedral arrangement—tetrahedral patchy particles—should be an ideal building block for this endeavor, their self-assembly into colloidal diamond has proved elusive. The kinetics of self-assembly pose a major challenge, with competition from an amorphous glassy phase, as well as clathrate crystals, leaving a narrow widow of patch widths where tetrahedral patchy particles can self-assemble into diamond crystals. Here we demonstrate that a two-component system of tetrahedral patchy particles, where bonding is allowed only between particles of different types to select even-member rings, undergoes crystallization into diamond crystals over a significantly wider range of patch widths conducive for experimental fabrication. We show that the crystallization in the two-component system is both thermodynamically and kinetically enhanced, as compared to the one-component system. Although our bottom-up route does not lead to the selection of the cubic polytype exclusively, we find that the cubicity of the self-assembled crystals increases with increasing patch width. Our designer system not only promises a scalable bottom-up route for colloidal diamond but also offers fundamental insight into crystallization into open lattices.

On the calculation of free energies over Hamiltonian and order parameters via perturbation and thermodynamic integration

In this work, complementary formulas are presented to compute free-energy differences via perturbation (FEP) methods and thermodynamic integration (TI). These formulas are derived by selecting only the most statistically significant data from the information extractable from the simulated points involved. On the one hand, commonly used FEP techniques based on overlap sampling leverage the full information contained in the overlapping macrostate probability distributions. On the other hand, conventional TI methods only use information on the first moments of those distributions, as embodied by the first derivatives of the free energy. Since the accuracy of simulation data degrades considerably for high-order moments (for FEP) or free-energy derivatives (for TI), it is proposed to consider, consistently for both methods, data up to second-order moments/derivatives. This provides a compromise between the limiting strategies embodied by common FEP and TI and leads to simple, optimized expressions to evaluate free-energy differences. The proposed formulas are validated with an analytically solvable harmonic Hamiltonian (for assessing systematic errors), an atomistic system (for computing the potential of mean force with coordinate-dependent order parameters), and a binary-component coarse-grained model (for tracing a solid-liquid phase diagram in an ensemble sampled through alchemical transformations). It is shown that the proposed FEP and TI formulas are straightforward to implement, perform similarly well, and allow robust estimation of free-energy differences even when the spacing of successive points does not guarantee them to have proper overlapping in phase space.

Low Interfacial Free Energy Describes the Bulk Ordering Transition in Colloidal Cubes

Many hard faceted nanoparticles are known to undergo disorder-to-order phase transitions following a classical nucleation and growth mechanism. In a previous study [J. Phys. Chem. B 2018, 122, 9264-9273], it was shown that hard cubes undergo a nonclassical phase transition with a bulk character instead of originating from consolidated nuclei. Significantly, an unusually high fraction of ordered particles was observed in the metastable basin of the disordered phase, even for very low degrees of supersaturation. This work aims to substantiate the conjecture that these unique properties originate from a comparatively low interfacial free energy between the disordered and ordered phases for hard cubes relative to other hard particle systems. Using the cleaving wall method to directly measure the interfacial free energy for cubes, it is found that its values are indeed small; e.g., at phase coexistence conditions, it is only one-fifth that for hard spheres. A theoretical nucleation model is used to explore the broader implications of low interfacial tension values and how this could result in a bulk ordering mechanism.

Self-Assembly of Patchy Colloidal Rods into Photonic Crystals Robust to Stacking Faults

Diamond-structured colloidal photonic crystals are much sought-after for their applications in visible light management because of their ability to support a complete photonic band gap (PBG). However, their realization via self-assembly pathways is a long-standing challenge. This challenge is rooted in three fundamental problems: the design of building blocks that assemble into diamond-like structures, the sensitivity of the PBG to stacking faults, and ensuring that the PBG opens at an experimentally attainable refractive index. Here we address these problems simultaneously using a multipronged computational approach. We use reverse engineering to establish the design principles for the rod-connected diamond structure (RCD), the so-called "champion" photonic crystal. We devise two distinct self-assembly routes for designer triblock patchy colloidal rods, both proceeding via tetrahedral clusters to yield a mixed phase of cubic and hexagonal polymorphs closely related to RCD. We use Monte Carlo simulations to show how these routes avoid a metastable amorphous phase. Finally, we show that both the polymorphs support spectrally overlapping PBGs. Importantly, randomly stacked hybrids of these polymorphs also display PBGs, thus circumventing the requirement of polymorph selection in a scalable fabrication method.

Parasitic crystallization of colloidal electrolytes: growing a metastable crystal from the nucleus of a stable phase

Colloidal particles have been extensively used to comprehend the main principles governing liquid-crystal nucleation. Multiple mechanisms and frameworks have been proposed, through either experiments or computational approaches, to rationalise the ubiquitous formation of colloidal crystals. In this work, we elucidate the nucleation scenario behind the crystallization of oppositely charged colloids. By performing molecular dynamics simulations of colloidal electrolytes in combination with the Seeding technique, we evaluate the fundamental factors, such as the nucleation rate, free energy barrier, surface tension and kinetic pre-factor, that determine the liquid-to-solid transition of several crystalline polymorphs. Our results show that at a high packing fraction, there is a cross-over between the nucleation of the CsCl structure and that of a substitutionally disordered fcc phase, despite the CuAu crystal being the most stable phase. We demonstrate that the crucial factor in determining which phase nucleates the fastest is the free energy cost of the cluster formation rather than their kinetic ability to grow from the liquid. While at a low packing fraction, the stable phase, CsCl, is the one that nucleates and subsequently grows, we show how at moderate and high packing fractions, a disordered fcc phase subsequently grows regardless of the nature of the nucleating phase, termed parasitic crystallization. Taken together, our results provide a panoramic perspective of the complex nucleation scenario of oppositely charged colloids at moderate temperature and rationalise the different thermodynamic and kinetic aspects behind it.

Interfacial free energy of a liquid-solid interface: Its change with curvature

We analyze the changes in the interfacial free energy between a spherical solid cluster and a fluid due to the change of the radius of the solid. Interfacial free energies from nucleation studies using the seeding technique for four different systems, being hard spheres, Lennard-Jones, and two models of water (mW and TIP4P/ICE), were plotted as a function of the inverse of the radius of the solid cluster. In all cases, the interfacial free energy was a linear function of the inverse of the radius of the solid cluster and this is consistent with Tolman's equation. This linear behavior is shown not only in isotherms but also along isobars. The effect of curvature on the interfacial free energy is more pronounced in water, followed by hard spheres, and smaller for Lennard-Jones particles. We show that it is possible to estimate nucleation rates of Lennard-Jones particles at different pressures by using information from simple NpT simulations and taking into account the variation of the interfacial free energy with the radius of the solid cluster. Neglecting the effects of the radius on the interfacial free energy (capillarity approximation) leads to incorrect values of the nucleation rate. For the Lennard-Jones system, the homogeneous nucleation curve is not parallel to the melting curve as was found for water in previous work. This is due to the increase in the interfacial free energy along the coexistence curve as the pressure increases. This work presents a simple and relatively straightforward way to approximately estimate nucleation rates.