Hit and miss of classical nucleation theory as revealed by a molecular simulation study of crystal nucleation in supercooled sulfur hexafluoride

Classical nucleation theory pictures the homogeneous nucleation of a crystal as the formation of a spherical crystalline embryo, possessing the properties of the macroscopic crystal, inside a parent supercooled liquid. In this work we study crystal nucleation in moderately supercooled sulfur hexafluoride by umbrella sampling simulations. The nucleation free energy evolves from 5.2kBT at T=170 K to 39.1kBT at T=195 K. The corresponding critical nucleus size ranges from 40 molecules at T=170 K to 266 molecules at T=195 K. Both nucleation free energy and critical nucleus size are shown to evolve with temperature according to the equations derived from the classical nucleation theory. Inspecting the obtained nuclei we show, however, that they present quite anisotropic shapes in opposition to the spherical assumption of the theory. Moreover, even though the critical nuclei possess the structure of the stable bcc plastic phase, the only mechanically stable crystal phase for SF6 in the temperature range investigated, they are shown to be less ordered than the corresponding macroscopic crystal. Their crystalline order is nevertheless shown to increase regularly with their size. This is confirmed by a study of a nucleus growth from a critical size to a size of the order of 10(4) molecules. Similarly to the fact that it does not affect the temperature dependence of the nucleation free energy and of the critical nucleus size, the ordering of the nucleus with size does not affect the growth rate of the nucleus.     

How does a transient amorphous precursor template crystallization

Crystallization through metastable phases, such as polymorphism, plays an important role in chemical manufacture, biomineralization, and protein crystallization. However, the kinetics creating the final stable crystalline phase from metastable phases has so far remained unclear. In this study, crystallization via an amorphous precursor, the so-called multistep crystallization (MSC), is studied quantitatively in a colloidal model system. In MSC, amorphous dense droplets are first nucleated from the mother phase. Subsequently, a few unstable subcrystalline nuclei can be created simultaneously by fluctuation from the tiny dense droplets, which is different from previous theoretical predictions. It is necessary for these crystalline nuclei to reach a critical size N*(crys) to become stable. However, in contrast to subcrystalline nuclei, a stable mature crystalline nucleus is not created by fluctuation but by coalescence of subcrystalline nuclei, which is unexpected. To accommodate a mature crystalline nucleus larger than the critical size N*(crys), the dense droplets have to first acquire a critical size N*. This implies that only a fraction of amorphous dense droplets can serve as a precursor of crystal nucleation. As an outcome, the overall nucleation rate of the crystalline phase is, to a large extent, determined by the nucleation rate of crystals in the dense droplets, which is much lower than the previous theoretical expectation. Furthermore, it is surprising to see that MSC will promote the production of defect-free crystals. The knowledge acquired in this study will also significantly advance our understandings in polymorphism related processes.     

Spinodal for the solution-to-crystal phase transformation

The formation of crystalline nuclei from solution has been shown for many systems to occur in two steps: the formation of quasidroplets of a disordered intermediate, followed by the nucleation of ordered crystalline embryos within these droplets. The rate of each step depends on a respective free-energy barrier and on the growth rate of its near-critical clusters. We address experimentally the relative significance of the free-energy barriers and the kinetic factors for the nucleation of crystals from solution using a model protein system. We show that crystal nucleation is 8-10 orders of magnitude slower than the nucleation of dense liquid droplets, i.e., the second step is rate determining. We show that at supersaturations of three or four k(B)T units, crystal nuclei of five, four, or three molecules transform into single-molecule nuclei, i.e., the significant nucleation barrier vanishes below the thermal energy of the molecules. We show that the main factor, which determines the rate of crystal nucleation, is the slow growth of the near-critical ordered clusters within the quasidroplets of the disordered intermediate. Analogous to the spinodal in supersaturated fluids, we define a solution-to-crystal spinodal from the transition to single-molecule crystalline nuclei. We show that heterogeneous nucleation centers accelerate nucleation not only because of the wettinglike effects that lower the nucleation barrier, as envisioned by classical theory, but by helping the kinetics of growth of the ordered crystalline embryos.     

Liquid Metastable Precursors of Ibuprofen as Aqueous Nucleation Intermediates

The nucleation mechanism of crystals of small organic molecules, postulated based on computer simulations, still lacks experimental evidence. In this study we designed an experimental approach to monitor the early stages of the crystallization of ibuprofen as a model system for small organic molecules. Ibuprofen undergoes liquid–liquid phase separation prior to nucleation. The binodal and spinodal limits of the corresponding liquid–liquid miscibility gap were analyzed and confirmed. An increase in viscosity sustains the kinetic stability of the dense liquid intermediate. Since the distances between ibuprofen molecules within the dense liquid phase are similar to those in the crystal forms, this dense liquid phase is identified as a precursor phase in the nucleation of ibuprofen, in which densification is followed by generation of structural order. This discovery may make it possible to enrich poorly soluble pharmaceuticals beyond classical solubility limitations in aqueous environments.     

Calculation of solid-liquid interfacial free energy: a classical nucleation theory based approach

We present a simple approach to calculate the solid-liquid interfacial free energy. This new method is based on the classical nucleation theory. Using the molecular dynamics simulation, we employ spherical crystal nuclei embedded in the supercooled liquids to create an ideal model of a homogeneous nucleation. The interfacial free energy is extracted by fitting the relation between the critical nucleus size and the reciprocal of the critical undercooling temperature. The orientationally averaged interfacial free energy is found to be 0.302+/-0.002 (in standard LJ unit). The temperature dependence of the interfacial free energy is also obtained in this work. We find that the interfacial free energy increases slightly with increasing temperature. The positive temperature coefficient of the interfacial free energy is in qualitative agreement with Spaepen's analysis [Solid State Phys. 47, FS181 (1994)] and Turnbull's empirical estimation [J. Appl. Phys. 21, 1022 (1950)].     

Crystal nucleation and the solid-liquid interfacial free energy

We present the results of molecular dynamics simulation of crystal nucleation in a supercooled Lennard-Jones liquid. Temperature and baric dependences of the nucleation rate, the Zeldovich factor, nucleus size diffusion coefficient, the radius, and the pressure in a critical crystal nucleus are defined in computer simulation. The data obtained have been used in the framework of classical nucleation theory to calculate the effective surface energy of crystal nuclei γ(e). It is shown that the value of γ(e) at T = const exceeds the value of the interfacial free energy at a flat crystal-liquid interface γ(∞) and γ(e) < γ(∞) at p = const.     

Nucleation of Organic Crystals—A Molecular Perspective

The outcome of synthetic procedures for crystalline organic materials strongly depends on the first steps along the molecular self‐assembly pathway, a process we know as crystal nucleation. New experimental techniques and computational methodologies have spurred significant interest in understanding the detailed molecular mechanisms by which nuclei form and develop into macroscopic crystals. Although classical nucleation theory (CNT) has served well in describing the kinetics of the processes involved, new proposed nucleation mechanisms are additionally concerned with the evolution of structure and the competing nature of crystallization in polymorphic systems. In this Review, we explore the extent to which CNT and nucleation rate measurements can yield molecular‐scale information on this process and summarize current knowledge relating to molecular self‐assembly in nucleating systems.     

A study of nucleation in supersaturated ibuprofen vapor in a flow diffusion chamber

Isothermal nucleation of supersaturated ibuprofen racemate vapor has been experimentally studied in a flow diffusion chamber at 293.3 and 301.2 K. Nucleation rates have been measured in the range of 10⁴?10⁴ cm^?3 s^?1 as functions of supersaturation. According to the first nucleation theorem, the numbers of molecules in critical nuclei have been found and used to determine the nucleation rate and supersaturation values as depending on the sizes of critical nuclei. The comparison of the experimental data with theoretical predictions has shown that the nucleation rates measured as functions of the numbers of molecules in critical nuclei are higher than the rates predicted by the classical theory by six to seven decimal orders of magnitude but, within one order of magnitude, coincide with the rates predicted by a theory previously proposed in a work by one of the authors, in which nucleation clusters were considered to be microscopic objects.     

Secondary nucleation and growth of ZnO

Recently we discovered that under certain conditions new crystal growth (branch) can be induced on specific crystalline planes of the same material. This is a new phenomenon and is in sharp contrast to typical nucleation and growth in which a crystal will simply grow larger in preferred directions depending on the surface energy of the specific crystalline planes. Based on our observation, we developed a sequential nucleation and growth technique offering the power to assemble complex hierarchical crystals step-by-step. However, the key questions of when and how the secondary nucleation takes place have not been answered. Here we systematically study secondary ZnO crystal growth using organic diamine additives with a range of chain lengths and concentration. We found that ZnO branches form for a narrow diamine concentration range with a critical lower and upper critical nucleation concentration limit, which increases by about a factor of 5 for each additional carbon in the diaminoalkane chain. Our results suggest that the narrow window for secondary growth is dictated by the solubility of the ZnO crystals, where the low critical nucleation concentration is determined by slight etching of the surface to produce new nucleation sites, and the upper critical concentration is determined by the supersaturation concentration. Kinetic measurements show that the induction time and growth rate increase with increasing diamine concentration and follow classical nucleation and growth theory. Observations of branch morphological evolution reveal the mechanisms guiding the tunable crystal size and morphology.     

Heterogeneous multicomponent nucleation theorems for the analysis of nanoclusters

In this paper we present a new form of the nucleation theorems applicable to heterogeneous nucleation. These heterogeneous nucleation theorems allow, for the first time, direct determination of properties of nanoclusters formed on pre-existing particles from measured heterogeneous nucleation probabilities. The theorems can be used to analyze the size (first theorem) and the energetics (second theorem) of heterogeneous clusters independent of any specific nucleation model. We apply the first theorem to the study of small water and n-propanol clusters formed at the surface of 8 nm silver particles. According to the experiments the size of the two-component critical clusters is found to be below 90 molecules, and only less than 20 molecules for pure water, less than 300 molecules for pure n-propanol. These values are drastically smaller than the ones predicted by the classical nucleation theory, which clearly indicates that the nucleating clusters are too small to be quantitatively described using a macroscopic theory.     

Universal Asymptotics of the Thermodynamic Characteristics of the Nucleation on Small Macroscopic Nuclei of Soluble Surfactants

Asymptotic behavior of thermodynamic characteristics of nucleation on small macroscopic nuclei of soluble surfactants at their complete dissolution in a nucleating droplet is studied. It is taken into account that, in the region of small sizes of nuclei and corresponding small sizes of critical nuclei of liquid phase, the chemical potential of condensate and the work of droplet formation are affected by the presence of dense surfactant adsorption monolayer on the droplet surface. It is shown that, as the limiting surface area per surfactant molecule in adsorption monolayer increases, the behavior of thermodynamic characteristics of nucleation in the region of small nucleus sizes is characterized by the transition from asymptotics at the adsorption of almost all substance comprising nucleus in a monolayer to the asymptotics at constant adsorption. The study performed is not limited by the selection of specific adsorption isotherms; therefore, the obtained asymptotic dependences of thermodynamic characteristics on the nucleus size can be considered as universal for the heterogeneous nucleation on the nuclei of soluble surfactants.     

Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: II. Nucleation in the metastable liquid immiscibility region

In the second part of our paper, we address crystal nucleation in the metastable liquid miscibility region of eutectic systems that is always present, though experimentally often inaccessible. While this situation resembles the one seen in single component crystal nucleation in the presence of a metastable vapor-liquid critical point addressed in previous works, it is more complex because of the fact that here two crystal phases of significantly different compositions may nucleate. Accordingly, at a fixed temperature below the critical point, six different types of nuclei may form: two liquid-liquid nuclei: two solid-liquid nuclei; and two types of composite nuclei, in which the crystalline core has a liquid "skirt," whose composition falls in between the compositions of the solid and the initial liquid phases, in addition to nuclei with concentric alternating composition shells of prohibitively high free energy. We discuss crystalline phase selection via exploring/identifying the possible pathways for crystal nucleation.     

Rate of homogeneous crystal nucleation in molten NaCl

We report a numerical simulation of the rate of crystal nucleation of sodium chloride from its melt at moderate supercooling. In this regime nucleation is too slow to be studied with "brute force" molecular-dynamics simulations. The melting temperature of ("Tosi Fumi") NaCl is approximately 1060 K. We studied crystal nucleation at T = 800 and 825 K. We observe that the critical nucleus formed during the nucleation process has the crystal structure of bulk NaCl. Interestingly, the critical nucleus is clearly faceted, the nuclei have a cubical shape. We have computed the crystal-nucleation rate using two completely different approaches, one based on an estimate of the rate of diffusive crossing of the nucleation barrier, the other based on the forward flux sampling and transition interface sampling methods. We find that the two methods yield the same result within an order of magnitude. However, when we compare the extrapolated simulation data with the only available experimental results for NaCl nucleation, we observe a discrepancy of nearly five orders of magnitude. We discuss the possible causes for this discrepancy.     

How a Facet of a Nanocrystal Is Formed: The Concept of the Symmetry Based Kinematic Theory

Conventional nanocrystal (NC) growth mechanisms have overwhelmingly focused on the final exposed facets to explain shape evolution. However, how the final facets are formed from the initial nuclei or seeds, has not been specifically interrogated. In this concept paper, we would like to concentrate on this specific topic, and introduce the symmetry based kinematic theory (SBKT) to explain the formation and evolution of crystal facets. It is a crystallographic theory based on the classical crystal growth concepts developed to illustrate the shape evolution during the NC growth. The most important principles connecting the basic NC growth processes and morphology evolution are the preferential growth directions and the properties of kinematic waves. On the contrary, the final facets are just indications of how the crystal growth terminates, and their formation and evolution rely on the NC growth processes: surface nucleation and layer advancement. Accordingly, the SBKT could even be applied to situations where non-faceted NCs such as spheres are formed.     

Numerical prediction of absolute crystallization rates in hard-sphere colloids

Special computational techniques are required to compute absolute crystal nucleation rates of colloidal suspensions. Using crystal nucleation of hard-sphere colloids as an example, we describe in some detail the novel computational tools that are needed to perform such calculations. In particular, we focus on the definition of appropriate order parameters that distinguish liquid from crystal, and on techniques to compute the kinetic prefactor that enters in the expression for the nucleation rate. In addition, we discuss the relation between simulation results and theoretical predictions based on classical nucleation theory.     

Thermodynamics and Kinetics of Prenucleation Clusters,Classical and Non‐Classical Nucleation

Recent observations of prenucleation species and multi‐stage crystal nucleation processes challenge the long‐established view on the thermodynamics of crystal formation. Here, we review and generalize extensions to classical nucleation theory. Going beyond the conventional implementation as has been used for more than a century now, nucleation inhibitors, precursor clusters and non‐classical nucleation processes are rationalized as well by analogous concepts based on competing interface and bulk energy terms. This is illustrated by recent examples of species formed prior to/instead of crystal nucleation and multi‐step nucleation processes. Much of the discussed insights were obtained from molecular simulation using advanced sampling techniques, briefly summarized herein for both nucleation‐controlled and diffusion‐controlled aggregate formation.     

Molecular simulation of the homogeneous crystal nucleation of carbon dioxide

We report on a molecular simulation study of the homogeneous nucleation of CO2 in the supercooled liquid at low pressure (P = 5 MPa) and for degrees of supercooling ranging from 32% to 60%. In all cases, regardless of the degree of supercooling, the structure of the crystal nuclei is that of the Pa3 phase, the thermodynamically stable phase. For the more moderate degree of supercooling of 32%, the nucleation is an activated process and requires a method to sample states of high free energy. In this work, we apply a series of bias potentials, which promote the ordering of the centers of mass of the molecules and allow us to gradually grow crystal nuclei. The reliability of the results so obtained is assessed by studying the evolution of the nuclei in the absence of any bias potential, and by determining their probability of growth. We estimate that the size of the critical nucleus, for which the probability of growth is 0.5, is approximately 240 molecules. Throughout the nucleation process, the crystal nuclei clearly exhibit a Pa3 structure, in apparent contradiction with Ostwald's rule of stages. The other polymorphs have a much larger free energy. This makes their formation highly unlikely and accounts for the fact that the nucleation of CO2 proceeds directly in the stable Pa3 structure.