Stresses inside critical nuclei

The usual derivation of classical nucleation theory is inappropriate for crystal nucleation. In particular, it leads to a seriously flawed estimate of the pressure inside a critical nucleus. This has consequences for the prediction of possible metastable phases during the nucleation process. In this paper, we reanalyze the theory for crystal nucleation based on the thermodynamics of small crystals suspended in a liquid, due to Mullins (J. Chem. Phys. 1984, 81, 1436). As an illustration of the difference between the classical picture and the present approach, we consider a numerical study of crystal nucleation in binary mixtures of hard spherical colloids with a size ratio of 1:10. The stable crystal phase of this system can be either dense or expanded. We find that, in the vicinity of the solid-solid critical point where the crystallites are highly compressible, small crystal nuclei are less dense than large nuclei. This phenomenon cannot be accounted for by either classical nucleation theory or by the Gibbsian droplet model.     

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.     

Thermodynamics of crystalline nanonucleus formation from liquid phase

The energy of crystal nucleation from liquid phase was considered, with the following two stages taken into account: (1) the formation of metastable supercooled melt (solution), containing pre-nuclei with intermediate amorphous (quasicrystalline) structure, and (2) the transformation of amorphous clusters into solid crystalline nuclei having different structures. With growth of a nucleus the nucleation energy profile manifests 2–3 maxima corresponding to these stages, and the kinetics of the non-stationary nucleation has five characteristic variations.     

Crystal Nucleation Rates for Particles Experiencing Short-Range Attractions: Applications to Proteins

A kinetic model for the nucleation of a crystalline phase consisting of particles experiencing short-range attractions is developed. Of particular significance is the proximity of the metastable fluid/fluid phase boundary. The model incorporates self-consistent thermodynamics, changes in gradient diffusivity, and density fluctuations in the vicinity of the critical point. Density fluctuations associated with the spinodal of this metastable phase transition greatly enhance nucleation rates, suggesting that experimental conditions may be found where rapid nucleation and slow crystal growth can be achieved by moving the metastable critical point relative to the solubility boundary. Copyright 2000 Academic Press.     

Free-energy landscape of nucleation with an intermediate metastable phase studied using capillarity approximation

Capillarity approximation is used to study the free-energy landscape of nucleation when an intermediate metastable phase exists. The critical nucleus that corresponds to the saddle point of the free-energy landscape as well as the whole free-energy landscape can be studied using this capillarity approximation, and various scenarios of nucleation and growth can be elucidated. In this study, we consider a model in which a stable solid phase nucleates within a metastable vapor phase when an intermediate metastable liquid phase exists. We predict that a composite critical nucleus that consists of a solid core and a liquid wetting layer as well as pure liquid and pure solid critical nuclei can exist depending not only on the supersaturation of the liquid phase relative to that of the vapor phase but also on the wetting behavior of the liquid surrounding the solid. The existence of liquid critical nucleus indicates that the phase transformation from metastable vapor to stable solid occurs via the intermediate metastable liquid phase, which is quite similar to the scenario of nucleation observed in proteins and colloidal systems. By studying the minimum-free-energy path on the free-energy landscape, we can study the evolution of the composition of solid and liquid within nuclei which is not limited to the critical nucleus.     

Can amorphous nuclei grow crystalline clathrates? The size and crystallinity of critical clathrate nuclei

Recent studies reveal that amorphous intermediates are involved in the formation of clathrate hydrates under conditions of high driving force, raising two questions: first, how could amorphous nuclei grow into crystalline clathrates and, second, whether amorphous nuclei are intermediates in the formation of clathrate crystals for temperatures close to equilibrium. In this work, we address these two questions through large-scale molecular simulations. We investigate the stability and growth of amorphous and crystalline clathrate nuclei and assess the thermodynamics and kinetic factors that affect the crystallization pathway of clathrates. Our calculations show that the dissociation temperature of amorphous clathrates is just 10% lower than for the crystals, facilitating the formation of metastable amorphous intermediates. We find that, at any temperatures, the critical crystalline nuclei are smaller than critical amorphous nuclei. The temperature dependence of the critical nucleus size is well described by the Gibbs-Thomson relation, from which we extract a liquid-crystal surface tension in excellent agreement with experiments. Our analysis suggests that at high driving force the amorphous nuclei may be kinetically favored over crystalline nuclei because of lower free energy barriers of formation. We investigated the role of the initial structure and size of the nucleus on the subsequent growth of the clathrates, and found that both amorphous and sI crystalline nuclei yield crystalline clathrates. Interestingly, growth of the metastable sII crystal polymorph is always favored over the most stable sI crystal, revealing kinetic control of the growth and indicating that a further step of ripening from sII to sI is needed to reach the stable crystal phase. The latter results are in agreement with the observed metastable formation of sII CO(2) and CH(4) clathrate hydrates and their slow conversion to sI under experimental conditions.     

Nucleation of self-associating fluids: free versus activated association

We use the density functional theory of statistical mechanics in a square gradient approximation to analyze the structure, size, and work of formation of critical nuclei in self-associating fluids where association reduces the strength of the interactions between bonded particles. This effect is expected in systems of strongly dipolar particles that associate into chains. In this work we analyze the nucleation behavior of two types of self-associating fluids: a system comprised of particles that can freely associate, and a system in which the association process involves a thermally activated initiation step. For the first case, we explore the properties of critical nuclei in fluids that exhibit a metastable critical point between a vapor phase and a highly associated liquid phase. In fluids where the association dynamics involves an initiation step, we investigate the nucleation behavior in the vicinity of the polymerization transition. In both cases critical nuclei undergo a structural transition that shares many of the features of the coil-globule transition reported in Monte Carlo simulations of strongly dipolar Stockmayer fluids. Our results suggest that the sharp structural transition observed in these simulations is evidence of the existence of a second-order or nearly second-order association transition in these model fluids.     

A metastable prerequisite for the growth of lumazine synthase crystals

Dense liquid phases, metastable with respect to a solid phase, form in solutions of proteins and small-molecule materials. They have been shown to serve as a prerequisite for the nucleation of crystals and other ordered solid phases. Here, using crystals of the protein lumazine synthase from Bacillus subtilis, which grow by the generation and spreading of layers, we demonstrate that within a range of supersaturations the only mechanism of generation of growth layers involves the association of submicrometer-size droplets of the dense liquid to the crystal surface. The dense liquid is metastable not only with respect to the crystals, but also with respect to the low-concentration solution: dynamic light scattering reveals that the droplets' lifetime is limited to several seconds, after which they decay into the low-concentration solution. The short lifetime does not allow growth to detectable dimensions so that liquid-liquid phase separation is not observed within a range of conditions broader than the one used for crystallization. If during their lifetime the droplets encounter a crystal surface, they lower their free energy not by decay, but by transformation into crystalline matter, ensuring perfect registry with the substrate. These observations illustrate two novel features of phase transformations in solutions: the existence of doubly metastable, short-lifetime dense phases and their crucial role for the growth of an ordered solid phase.     

Crystal Memory near Discontinuous Triacylglycerol Phase Transitions: Models,Metastable Regimes,and Critical Points

It is proposed that “crystal memory”, observed in a discontinuous solid-liquid phase transition of saturated triacylglycerol (TAG) molecules, is due to the coexistence of solid TAG crystalline phases and a liquid TAG phase, in a superheated metastable regime. Such a coexistence has been detected. Solid crystals can act as heterogeneous nuclei onto which molecules can condense as the temperature is lowered. We outlined a mathematical model, with a single phase transition, that shows how the time-temperature observations can be explained, makes predictions, and relates them to recent experimental data. A modified Vogel-Fulcher-Tammann (VFT) equation is used to predict time-temperature relations for the observation of “crystal memory” and to show boundaries beyond which “crystal memory” is not observed. A plot of the lifetime of a metastable state versus temperature, using the modified VFT equation, agrees with recent time-temperature data. The model can be falsified through its predictions: the model possesses a critical point and we outline a procedure describing how it could be observed by changing the hydrocarbon chain length. We make predictions about how thermodynamic functions will change as the critical point is reached and as the system enters a crossover regime. The model predicts that the phenomenon of “crystal memory” will not be observed unless the system is cooled from a superheated metastable regime associated with a discontinuous phase transition.     

Phase separations in liquid crystal-colloid mixtures

We present a mean-field theory to describe phase separations in mixtures of a nematic liquid crystal and a colloidal particle. The theory takes into account an orientational ordering of liquid crystals and a crystalline ordering of colloidal particles. We calculate phase diagrams on the temperature-concentration plane, depending on interactions between a liquid crystal and a colloidal surface and a coupling between nematic and crystalline ordering. We find various phase separation processes, such as a nematic-crystal phase separation and nematic-isotropic-crystal triple point. Inside binodal curves, we find new unstable and metastable regions which are important in phase ordering dynamics. We also find a stable nematic-crystalline (NC) phase, where colloidal particles dispersed in a nematic phase can form a crystalline structure. The coexistence between two NC phases with different concentrations can be appear though the coupling between nematic and crystalline ordering.     

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.     

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.     

Towards understanding the nucleation mechanism for multi-component systems: an atomistic simulation of the ternary nucleation of water/n-nonane/1-butanol

Inspired by the previous finding of some unusual vapour/liquid nucleation results on the ternary water/n-nonane/1-butanol system, atomistic simulations were carried out for a detailed investigation of this mixture. These simulations reproduced the experimentally-reported non-ideal nucleation behaviour for this system, including both onset activities and the average compositions of the critical nuclei. Close examination of the nucleation free energy data and the structure of the critical nuclei reveals two types of phase separation. One occurs internally inside the cluster via formation of a multi-layered structure. The other takes place externally, leading to the coexistence of multiple nucleation channels, characterized by critical clusters of different compositions. Such mechanistic and structural heterogeneity is the microscopic origin of the complex nucleation behaviour observed for this ternary mixture.     

Interface induced crystal structures of dioctyl-terthiophene thin films

Temperature dependent structural and morphological investigations on semiconducting dioctyl-terthiophene (DOTT) thin films prepared on silica surfaces reveals the coexistence of surface induce order and distinct crystalline/liquid crystalline bulk polymorphs. X-ray diffraction and scanning force microscopy measurements indicate that at room temperature two polymorphs are present: the surface induced phase grows directly on the silica interface and the bulk phase on top. At elevated temperatures the long-range order gradually decreases, and the crystal G (340 K), smectic F (348 K), and smectic C (360 K) phases are observed. Indexation of diffraction peaks reveals that an up-right standing conformation of DOTT molecules is present within all phases. A temperature stable interfacial layer close to the silica-DOTT interface acts as template for the formation of the different phases. Rapid cooling of the DOTT sample from the smectic C phase to room temperature results in freezing into a metastable crystalline state with an intermediated unit cell between the room temperature crystalline phase and the smectic C phase. The understanding of such interfacial induced phases in thin semiconducting liquid crystal films allows tuning of crystallographic and therefore physical properties within organic thin films.     

Phase transitions and conformational changes in an antiferroelectric liquid crystal 4-(1-methylheptyloxycarbonyl)phenyl 4'-octyloxybiphenyl-4-carboxylate (MHPOBC)

The thermodynamic properties and conformational structures in the crystalline and liquid crystalline phases of 4-(1-methylheptyloxycarbonyl)phenyl 4'-octyloxybiphenyl-4-carboxylate (MHPOBC) have been investigated by Raman scattering, X-ray diffraction and differential scanning calorimetry (DSC). MHPOBC can assume two different crystalline states at room temperature depending on crystallization conditions, and heating transforms the metastable crystal to the stable one before melting. In the liquid crystalline phases, Raman scattering spectra have revealed not only the flip-flop twisting vibration of the biphenyl group but also the internal rotation around C-C bonds between the carbonyl groups and the corresponding benzene rings.     

Nucleation on cylindrical plates: sharp transitions and double barriers

We apply methods of density-functional theory in statistical mechanics to study the properties of droplets and bubbles formed on a single cylindrical plate or between two such disks immersed in a metastable fluid. Our approach allows us to analyze the properties of different types of aggregates and investigate the effect of disk size, disk separation, and solid-fluid interactions on the dynamics of a liquid-vapor phase transition. The finite size of disks induces nucleation phenomena that are not observed in the cases of either a planar wall or a slit pore. Heterogeneous nucleation on a single disk is characterized by the existence of two distinct types of critical nuclei that control the phase-transition dynamics at different supersaturations. Asymmetric droplets or bubbles formed on one side of the disk are the preferred nucleation path at high supersaturations. However, these types of aggregates become unstable close to the binodal, where they abruptly collapse into nuclei that engulf the cylindrical plates. Droplet or bubble nucleation in between two disks may occur through a free-energy barrier with one or two maxima depending on the value of the system parameters and the supersaturation. Metastable droplets or bubbles corresponding to local minima of the free energy are observed forming between two plates only after density fluctuations in the system achieve a critical size. These types of aggregates only exist for cylindrical plates larger than a minimum size given a fixed distance between the disks. The stability of these droplets and bubbles decreases when the plates are separated.     

Enhanced protein crystallization around the metastable critical point

We report on a computer-simulation study of homogeneous crystal nucleation in a model for globular proteins. We find that the presence of a metastable vapour-liquid critical point drastically changes the pathway for the formation of a critical nucleus. But what is more important, the large density fluctuations near the critical point also lowers the free-energy barrier to nucleation and hence increases the nucleation rate. As␣the location of the vapour-liquid critical point can be controlled by changing the solvent conditions, our simulation results suggest a guided approach to protein crystallization. Received: 4 June 1998/Accepted: 3 September 1998 / Published online: 10 December 1998     

Dumbbells and onions in ternary nucleation

Molecular simulations for a ternary nucleation system (water/n-nonane/1-butanol) demonstrate a more complex nucleation mechanism than previously thought, where critical nuclei with different compositions are present even for a given vapour-phase composition; the spatial distribution in these critical nuclei is heterogeneous and dumbbell and onion motifs are found; in the former, water and nonane nano-droplets are connected through a butanol handle, while in the latter a water core is surrounded by a nonane corona with an interfacial butanol shell.     

Phase diagram of a two-dimensional system with anomalous liquid properties

Using Monte Carlo simulation techniques, we calculate the phase diagram for a square-shoulder square-well potential in two dimensions that has been previously shown to exhibit liquid anomalies consistent with a metastable liquid-liquid critical point. We consider the liquid, gas, and five crystal phases, and find that all the melting lines are first order, despite a small range of metastability. One melting line exhibits a temperature maximum, as well as a pressure maximum that implies inverse melting over a small range in pressure.     

Molecular dynamics studies of the kinetics of phase changes in clusters III: structures, properties, and crystal nucleation of iron nanoparticle Fe331

Molecular dynamics (MD) computer simulations have been carried out to study the structures, properties, and crystal nucleation of iron nanoparticles with 331 Fe atoms or with diameter around 2 nm. Structure information for the nanoparticles was analyzed from the MD simulations. Three crystalline phases and one amorphous phase were obtained by cooling the nanoparticles from their molten droplets at different cooling rates or with different lengths of cooling time periods. Molten droplets froze into three different solid phases and a solid-solid transition from a disordered body-centered cubic (BCC) phase to an ordered BCC phase were observed during the slow cooling and the quenching processes. Properties of nanoparticle Fe₃₃₁, such as melting point, freezing temperature, heat capacity, heat of fusion, heat of crystallization, molar volume, thermal expansion coefficient, and diffusion coefficient, have been estimated. Nucleation rates of crystallization to two solid phases for Fe₃₃₁ at temperatures of 750, 800, and 850 K are presented. Both classical nucleation theory and diffuse interface theory are used to interpret our observed nucleation results. The interfacial free energy and the diffuse interface thickness between the liquid phase and two different solid phases are estimated from these nucleation theories.