Fcc vs. hcp competition in colloidal hard-sphere nucleation: on their relative stability, interfacial free energy and nucleation rate

Modulating the Short-Range Correlation through the Assembly of Anisotropic Colloidal Clusters

Disordered structures with short-range correlation exhibit unique photonic or thermal properties, often arising from the self-assembly of monodisperse microspheres. Colloidal clusters, aggregates of a defined number of microspheres, have anisotropic shapes and offer a unique approach to tailor the short-range correlation compared to traditional isotropic microspheres. However, this potential remains largely unexplored. To this end, we have prepared triangular and tetrahedral colloidal clusters and coassemble them with polystyrene microspheres into both monolayers and thick films (∼10 μm) over 5 mm × 5 mm areas. Through combining scanning electron microscopy images and laser scattering measurements, we show that increasing cluster concentration from 0% to 40% increases disorder in both thin and thick films. Beyond 40%, the short-range order in thin films plateaus, while it increases with triangular clusters but decreases with tetrahedral clusters in thick films. The distinct roles of triangular and tetrahedral clusters offer new insights into design of correlated disordered materials for novel applications.

Statistical mechanics of crystal nuclei of hard spheres

In the study of crystal nucleation via computer simulations, hard spheres are arguably the most extensively explored model system. Nonetheless, even in this simple model system, the complex thermodynamics of crystal nuclei can sometimes give rise to counterintuitive results, such as the recent observation that the pressure inside a critical nucleus is lower than that of the surrounding fluid, seemingly clashing with the strictly positive Young-Laplace pressure we would expect in liquid droplets. Here, we re-derive many of the founding equations associated with crystal nucleation and use the hard-sphere model to demonstrate how they give rise to this negative pressure difference. We exploit the fact that, in the canonical ensemble, a nucleus can be in a (meta)stable equilibrium with the fluid and measure the surface stress for both flat and curved interfaces. Additionally, we explain the effect of defects on the chemical potential inside the crystal nucleus. Finally, we present a simple, fitted thermodynamic model to capture the properties of the nucleus, including the work required to form critical nuclei.

Entropy-Driven Crystallization of Hard Colloidal Mixtures of Polymers and Monomers

The most trivial example of self-assembly is the entropy-driven crystallization of hard spheres. Past works have established the similarities and differences in the phase behavior of monomers and chains made of hard spheres. Inspired by the difference in the melting points of the pure components, we study, through Monte Carlo simulations, the phase behavior of athermal mixtures composed of fully flexible polymers and individual monomers of uniform size. We analyze how the relative number fraction and the packing density affect crystallization and the established ordered morphologies. As a first result, a more precise determination of the melting point for freely jointed chains of tangent hard spheres is extracted. A synergetic effect is observed in the crystallization leading to synchronous crystallization of the two species. Structural analysis of the resulting ordered morphologies shows perfect mixing and thus no phase separation. Due to the constraints imposed by chain connectivity, the local environment of the individual spheres, as quantified by the Voronoi polyhedron, is systematically more spherical and more symmetric compared to that of spheres belonging to chains. In turn, the local environment of the ordered phase is more symmetric and more spherical compared to that of the initial random packing, demonstrating the entropic origins of the phase transition. In general, increasing the polymer content reduces the degree of crystallinity and increases the melting point to higher volume fractions. According to the present findings, relative concentration is another determining factor in controlling the phase behavior of hard colloidal mixtures based on polymers.

Ordered hierarchical superlattice amplifies coated-CeO2 nanoparticles luminescence

Achieving a controlled preparation of nanoparticle superstructures with spatially periodic arrangement, also called superlattices, is one of the most intriguing and open questions in soft matter science. The interest in such regular superlattices originates from the potentialities in tailoring the physicochemical properties of the individual constituent nanoparticles, eventually leading to emerging behaviors and/or functionalities that are not exhibited by the initial building blocks. Despite progress, it is currently difficult to obtain such ordered structures; the influence of parameters, such as size, softness, interaction potentials, and entropy, are neither fully understood yet and not sufficiently studied for 3D systems. In this work, we describe the synthesis and characterization of spatially ordered hierarchical structures of coated cerium oxide nanoparticles in water suspension prepared by a bottom-up approach. Covering the CeO2 surface with amphiphilic molecules having chains of appropriate length makes it possible to form ordered structures in which the particles occupy well-defined positions. In the present case superlattice arrangement is accompanied by an improvement in photoluminescence (PL) efficiency, as an increase in PL intensity of the superlattice structure of up to 400 % compared with that of randomly dispersed nanoparticles was observed. To the best of our knowledge, this is one of the first works in the literature in which the coexistence of 3D structures in solution, such as face-centered cubic (FCC) and Frank-Kasper (FK) phases, of semiconductor nanoparticles have been related to their optical properties.

Fluid-solid transitions in photonic crystals of soft, thermoresponsive microgels

Microgels are often discussed as well-suited model system for soft colloids. In contrast to rigid spheres, the microgel volume and, coupled to this, the volume fraction in dispersion can be manipulated by external stimuli. This behavior is particularly interesting at high packings where phase transitions can be induced by external triggers such as temperature in the case of thermoresponsive microgels. A challenge, however, is the determination of the real volume occupied by these deformable, soft objects and consequently, to determine the boundaries of the phase transitions. Here we propose core-shell microgels with a rigid silica core and a crosslinked, thermoresponsive poly-N-isopropylacrylamide (PNIPAM) shell with a carefully chosen shell-to-core size ratio as ideal model colloids to study fluid-solid transitions that are inducible by millikelvin changes in temperature. Specifically, we identify the temperature ranges where crystallization and melting occur using absorbance spectroscopy in a range of concentrations. Slow annealing from the fluid to the crystalline state leads to photonic crystals with Bragg peaks in the visible wavelength range and very narrow linewidths. Small-angle X-ray scattering is then used to confirm the structure of the fluid phase as well as the long-range order, crystal structure and microgel volume fraction in the solid phase. Thanks to the scattering contrasts and volume ratio of the cores with respect to the shells, the scattering data do allow for form factor analysis revealing osmotic deswelling at volume fractions approaching and also exceeding the hard sphere packing limit.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.

RSeeds: Rigid Seeding Method for Studying Heterogeneous Crystal Nucleation

Heterogeneous nucleation is the dominant form of liquid-to-solid transition in nature. Although molecular simulations are most uniquely suited to studying nucleation, the waiting time to observe even a single nucleation event can easily exceed the current computational capabilities. Therefore, there exists an imminent need for methods that enable computationally fast and feasible studies of heterogeneous nucleation. Seeding is a technique that has proven to be successful at dramatically expanding the range of computationally accessible nucleation rates in simulation studies of homogeneous crystal nucleation. In this article, we introduce a new seeding method for heterogeneous nucleation called Rigid Seeding (RSeeds). Crystalline seeds are treated as pseudorigid bodies and simulated on a surface with metastable liquid above its melting temperature. This allows the seeds to adapt to the surface and identify favorable seed-surface configurations, which is necessary for reliable predictions of crystal polymorphs that form and the corresponding heterogeneous nucleation rates. We demonstrate and validate RSeeds for heterogeneous ice nucleation on a flexible self-assembled monolayer surface, a mineral surface based on kaolinite, and two model surfaces. RSeeds predicts the correct ice polymorph, exposed crystal plane, and rotation on the surface. RSeeds is semiquantitative and can be used to estimate the critical nucleus size and nucleation rate when combined with classical nucleation theory. We demonstrate that RSeeds can be used to evaluate nucleation rates spanning many orders of magnitude.

Crystal Polymorph Selection Mechanism of Hard Spheres Hidden in the Fluid

Nucleation plays a critical role in the birth of crystals and is associated with a vast array of phenomena, such as protein crystallization and ice formation in clouds. Despite numerous experimental and theoretical studies, many aspects of the nucleation process, such as the polymorph selection mechanism in the early stages, are far from being understood. Here, we show that the hitherto unexplained excess of particles in a face-centered-cubic (fcc)-like environment, as compared to those in a hexagonal-close-packed (hcp)-like environment, in a crystal nucleus of hard spheres can be explained by the higher order structure in the fluid phase. We show using both simulations and experiments that in the metastable fluid phase, pentagonal bipyramids, clusters with fivefold symmetry known to be inhibitors of crystal nucleation, transform into a different cluster, Siamese dodecahedra. These clusters are closely similar to an fcc subunit, which explains the higher propensity to grow fcc than hcp in hard spheres. We show that our crystallization and polymorph selection mechanism is generic for crystal nucleation from a dense, strongly correlated fluid phase.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log₁₀ (J/(m³ s)) = -7(5).

Direct Calculation of the Interfacial Free Energy between NaCl Crystal and Its Aqueous Solution at the Solubility Limit

Salty water is the most abundant electrolyte aqueous mixture on Earth, however, very little is known about the NaCl-saturated solution interfacial free energy (γ_{s}). Here, we provide the first direct estimation of γ_{s} for several NaCl crystallographic planes by means of the mold integration technique, a highly efficient computational method to evaluate interfacial free energies with anisotropic crystal resolution. Making use of the JC-SPC/E model, one of the most benchmarked force fields for NaCl water solutions, we measure γ_{s} of four different crystal planes, (100), (110), (111), and (112[over ¯]) with the saturated solution at normal conditions. We find high anisotropy between the different crystal orientations with values ranging from 100 to 150  mJ m^{-2}, and the average value of the distinct planes being γ[over ¯]_{s}=137(20)  mJ m^{-2}. This value for the coexistence interfacial free energy is in reasonable agreement with previous extrapolations from nucleation studies. Our Letter represents a milestone in the computational calculation of interfacial free energies between ionic crystals and aqueous solutions.

Polymorph Stability and Free Energy of Crystallization of Freely-Jointed Polymers of Hard Spheres

The free energy of crystallization of monomeric hard spheres as well as their thermodynamically stable polymorph have been known for several decades. In this work, we present semianalytical calculations of the free energy of crystallization of freely-jointed polymers of hard spheres as well as of the free energy difference between the hexagonal closed packed (HCP) and face-centered cubic (FCC) polymorphs. The phase transition (crystallization) is driven by an increase in translational entropy that is larger than the loss of conformational entropy of chains in the crystal with respect to chains in the initial amorphous phase. The conformational entropic advantage of the HCP polymer crystal over the FCC one is found to be ΔschHCP-FCC≈0.331×10-5k per monomer (expressed in terms of Boltzmann's constant k). This slight conformational entropic advantage of the HCP crystal of chains is by far insufficient to compensate for the larger translational entropic advantage of the FCC crystal, which is predicted to be the stable one. The calculated overall thermodynamic advantage of the FCC over the HCP polymorph is supported by a recent Monte Carlo (MC) simulation on a very large system of 54 chains of 1000 hard sphere monomers. Semianalytical calculations using results from this MC simulation yield in addition a value of the total crystallization entropy for linear, fully flexible, athermal polymers of Δs≈0.93k per monomer.

Polymorphism and Perfection in Crystallization of Hard Sphere Polymers

We present results on polymorphism and perfection, as observed in the spontaneous crystallization of freely jointed polymers of hard spheres, obtained in an unprecedentedly long Monte Carlo (MC) simulation on a system of 54 chains of 1000 monomers. Starting from a purely amorphous configuration, after an initial dominance of the hexagonal closed packed (HCP) polymorph and a transitory random hexagonal close packed (rHCP) morphology, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. An analysis of chain conformational characteristics, of the spatial distribution of monomers and of the volume accessible to them shows that the phase transition is caused by an increase in translational entropy that is larger than the loss of conformational entropy of the chains in the crystal, compared to the amorphous state. In spite of the significant local re-arrangements, as reflected in the bending and torsion angle distributions, the average chain size remains unaltered during crystallization. Polymers in the crystal adopt ideal random walk statistics as their great length renders local conformational details, imposed by the geometry of the FCC crystal, irrelevant.

Kinetic Phase Diagram for Nucleation and Growth of Competing Crystal Polymorphs in Charged Colloids

We determine the kinetic phase diagram for nucleation and growth of crystal phases in a suspension of charged colloids. Exploiting the seeding approach in extensive simulations, we calculate nucleation barrier heights for face-centered cubic (fcc) and body-centered cubic (bcc) phases for varying screening lengths and supersaturations. We determine for the entire metastable fluid region the crystal polymorph with the lowest nucleation barrier, and find a regime close to the triple point where metastable bcc can form due to a lower nucleation barrier, even though fcc is the stable phase. For higher supersaturation, we find that the difference in barrier heights decreases and we observe a mix of hexagonal close-packed, fcc, and bcc structures in the growth of crystalline seeds as well as in spontaneously formed crystals. Our kinetic phase diagram rationalizes the different crystallization mechanisms observed in previous work.

Homogeneous ice nucleation rates for mW and TIP4P/ICE models through Lattice Mold calculations

Water freezing is the most common liquid-to-crystal phase transition on Earth, however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240K. We employ two different water models, mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic non-polarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih, Ic and a stacking mixture of ice Ih/Ic; reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner which contributes to benchmarking the freezing behaviour of two popular water models.