Phase field theory of interfaces and crystal nucleation in a eutectic system of fcc structure: I. Transitions in the one-phase liquid region

The phase field theory (PFT) has been applied to predict equilibrium interfacial properties and nucleation barrier in the binary eutectic system Ag-Cu using double well and interpolation functions deduced from a Ginzburg-Landau expansion that considers fcc (face centered cubic) crystal symmetries. The temperature and composition dependent free energies of the liquid and solid phases are taken from CALculation of PHAse Diagrams-type calculations. The model parameters of PFT are fixed so as to recover an interface thickness of approximately 1 nm from molecular dynamics simulations and the interfacial free energies from the experimental dihedral angles available for the pure components. A nontrivial temperature and composition dependence for the equilibrium interfacial free energy is observed. Mapping the possible nucleation pathways, we find that the Ag and Cu rich critical fluctuations compete against each other in the neighborhood of the eutectic composition. The Tolman length is positive and shows a maximum as a function of undercooling. The PFT predictions for the critical undercooling are found to be consistent with experimental results. These results support the view that heterogeneous nucleation took place in the undercooling experiments available at present. We also present calculations using the classical droplet model [classical nucleation theory (CNT)] and a phenomenological diffuse interface theory (DIT). While the predictions of the CNT with a purely entropic interfacial free energy underestimate the critical undercooling, the DIT results appear to be in a reasonable agreement with the PFT predictions.     

Correlations between surface and interface energies with respect to crystal nucleation

The energetic chemical and structural properties of interfaces between solid and liquid metals are of great interest in numerous applications. However, solid-liquid interfacial energies of metals are often determined by nucleation experiments, requiring particular attention to the measurements so obtained. The purpose of this paper is to conduct an analysis of the level of liquid undercooling of 3d-, 4d-, and 5d-transition metals, as well as of other common undercooled elements, through a critical survey of thermophysical data and experimental results. As already pointed out for the liquid-vapor surface energy sigma(LV), the solid-liquid interface energy sigma(LS) determined from the maximum amount of liquid undercooling is well connected to the position of the element in the periodic table, leading to the individualization of each behavior, including the size of the critical nucleus. The beta ratio of sigma(LS)/sigma(LV) is demonstrated to be an interesting dimensionless number by which to classify the elements into distinctive groups. No significant unexplained anomaly is identified, except for Co, supporting the belief that Turnbull's classical theory furnishes a robust support to describe the crystal nucleation in pure elements.     

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)].     

Atomistic theory of amyloid fibril nucleation

We consider the nucleation of amyloid fibrils at the molecular level when the process takes place by a direct polymerization of peptides or protein segments into β-sheets. Employing the atomistic nucleation theory (ANT), we derive a general expression for the work to form a nanosized amyloid fibril (protofilament) composed of successively layered β-sheets. The application of this expression to a recently studied peptide system allows us to determine the size of the fibril nucleus, the fibril nucleation work, and the fibril nucleation rate as functions of the supersaturation of the protein solution. Our analysis illustrates the unique feature of ANT that the size of the fibril nucleus is a constant integer in a given supersaturation range. We obtain the ANT nucleation rate and compare it with the rates determined previously in the scope of the classical nucleation theory (CNT) and the corrected classical nucleation theory (CCNT). We find that while the CNT nucleation rate is orders of magnitude greater than the ANT one, the CCNT and ANT nucleation rates are in very good quantitative agreement. The results obtained are applicable to homogeneous nucleation, which occurs when the protein solution is sufficiently pure and/or strongly supersaturated.     

Simple correction to the classical theory of homogeneous nucleation

Formation of the new disperse phase via homogeneous nucleation plays a fundamental role wherever the first-order phase transitions occur. Inconsistent temperature dependence of the nucleation rates and poor agreement of theoretical critical supersaturations with experimental data for a number of substances are fundamental problems of the classical nucleation theory (CNT). Here we show that these problems can be solved with a simple empirical correction to CNT. Despite its simplicity, the corrected CNT (CCNT) accurately predicts temperature dependences and absolute values of the critical supersaturations for both organic and inorganic substances with widely varying properties at different ambient conditions and it works surprisingly well in a wide size range down to few molecules. The difference in predictions of CCNT and other versions of the classical nucleation theory commonly used in analyzing experimental data is discussed. It has been found that CCNT consistently gives better agreement with experimental data than other versions of classical nucleation theory.     

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.     

Nucleation rate isotherms of argon from molecular dynamics simulations

We report six nucleation rate isotherms of vapor-liquid nucleation of Lennard-Jones argon from molecular dynamics simulations. The isotherms span three orders of magnitude in nucleation rates, 10(23)相似文献   

氟化钾团簇的相变、成核和重结晶的分子动力学模拟

我们利用Born-Mayer-Huggins相互作用势函数对(KF)_N(N=108，256，500和864)团簇进行了分子动力学(MD)模拟。为了避免周期性边界条件对相变、成核和重结晶的干扰作用，对体系采用了自由边界。基于MD模拟结果，对团簇的熔化温度、熔化焓、扩散系数、成核速率、固液界面自由能、临界核大小等进行了计算和讨论。在对(KF)₈₆₄双晶团簇的热退火模拟中，观察到了固态的重结晶和晶粒的生长。经典的成核理论成功地解释了(KF)₈₆₄双晶团簇的重结晶MD模拟结果。     

Test of classical nucleation theory on deeply supercooled high-pressure simulated silica

We test classical nucleation theory (CNT) in the case of simulations of deeply supercooled, high density liquid silica, as modeled by the van Beest-Kramer-van Santen potential. We find that at density rho=4.38 gcm(3), spontaneous nucleation of crystalline stishovite occurs in conventional molecular dynamics simulations at temperature T=3000 K, and we evaluate the nucleation rate J directly at this T via "brute force" sampling of nucleation events in numerous independent runs. We then use parallel, constrained Monte Carlo simulations to evaluate DeltaG(n), the free energy to form a crystalline embryo containing n silicon atoms, at T=3000, 3100, 3200, and 3300 K. By comparing the form of DeltaG(n) to CNT, we test the ability of CNT to reproduce the observed behavior as we approach the regime where spontaneous nucleation occurs on simulation time scales. We find that the prediction of CNT for the n dependence of DeltaG(n) fits reasonably well to the data at all T studied. Deltamu, the chemical potential difference between bulk liquid and stishovite, is evaluated as a fit parameter in our analysis of the form of DeltaG(n). Compared to directly determined values of Deltamu extracted from previous work, the fitted values agree only at T=3300 K; at lower T the fitted values increasingly overestimate Deltamu as T decreases. We find that n(*), the size of the critical nucleus, is approximately ten silicon atoms at T=3300 K. At 3000 K, n(*) decreases to approximately 3, and at such small sizes methodological challenges arise in the evaluation of DeltaG(n) when using standard techniques; indeed even the thermodynamic stability of the supercooled liquid comes into question under these conditions. We therefore present a modified approach that permits an estimation of DeltaG(n) at 3000 K. Finally, we directly evaluate at T=3000 K the kinetic prefactors in the CNT expression for J, and find physically reasonable values; e.g., the diffusion length that Si atoms must travel in order to move from the liquid to the crystal embryo is approximately 0.2 nm. We are thereby able to compare the results for J at 3000 K obtained both directly and based on CNT, and find that they agree within an order of magnitude. In sum, our work quantifies how certain predictions of CNT (e.g., for Deltamu) break down in this deeply supercooled limit, while others [the n dependence of DeltaG(n)] are not as adversely affected.     

Critical nucleus size effects on nanoparticle formation in microemulsions: a comparison study between experimental and simulation results

Monte Carlo simulations were performed to study the influence of critical nucleus size on nanoparticle formation in microemulsions. It was found that critical nucleus size strongly affected nucleation and growth rates, as well as final nanoparticle sizes. An increase of critical nucleus leads to a slower nucleation process. In contrast, it gives rise to acceleration of the growth process. Final nanoparticle sizes increase as the critical nucleus value increases. It is predicted that this dependence will be less pronounced when a high reactant concentration is used. We have compared the simulation results with experimental data taken from different authors. Good agreement between the two kinds of results supports the conclusions of this paper.     

Direct measurement of critical nucleus size in confined volumes

In crystallization, the critical nucleus size is of pivotal importance. Above this size, it is favorable for the new crystalline phase to form; below this size, the clusters will tend to dissolve rather than grow. To date, there has been no direct method for measuring the critical nucleus size. Instead, the size is typically calculated from the variation of crystallization rates with temperature. This involves using bulk values of the interfacial tension and enthalpy of fusion, which are inappropriate for small critical nucleus sizes. Here, we present a direct method for measuring the size of the critical nucleus, based on observing crystallization temperatures of materials within microemulsions. Using this approach, the number of molecules in the critical nucleus can be found simply by measuring the droplet size. Data on the freezing of water in water-in-oil microemulsions with and without the nucleating agent, heptacosanol, are presented to support our hypothesis. The results show that the critical nucleus contains 90-350 ice molecules for water pool radii of approximately 1.2-1.8 nm for the heptacosanol-doped microemulsions in which heterogeneous nucleation is initiated at the droplet interface. For the microemulsions without heptacosanol, the critical nucleus contains 70-210 ice molecules for water pool radii of approximately 1.2-1.8 nm. The smaller values arise because homogeneous nucleation occurs and therefore the crystallization temperatures are lower. We can also determine how bulk properties are perturbed at the nanoscale, and we find that the ratio of the ice-water interfacial tension to the enthalpy of fusion decreases significantly for water pool radii that are <2 nm.     

Electrochemical nucleation: comparison test of classical and atomistic nucleation models

Classical and atomistic nucleation models have been tested in several aqueous systems dealing with electrocrystallisation. A lot of reported experimental nucleation data have been used and in a wide range of overpotentials. The critical nucleus size has been calculated in those cases not reported in the original work, and the results obtained with the classical and atomistic models have been tabulated, compared and discussed. Small values for the critical nucleus size occur in most of the systems.     

Homogeneous nucleation rates of higher n-alcohols measured in a laminar flow diffusion chamber

Nucleation rate isotherms of n-butanol, n-pentanol, n-hexanol, n-heptanol, and n-octanol were measured in a laminar flow diffusion chamber using helium as carrier gas. The measurements were made at 250-310 K, corresponding to reduced temperatures of 0.43-0.50, and at atmospheric pressure. Experimental nucleation rate range was from 10(3) to 10(7) cm(-3) s(-1). The expression and accuracy of thermodynamic parameters, in particular equilibrium vapor pressure, were found to have a significant effect on calculated nucleation rates. The results were compared to the classical nucleation theory (CNT), the self-consistency corrected classical theory (SCC) and the Hale's scaled model of the CNT. The average ratio between the experimental and theoretical nucleation rates for all alcohols used was 1.5x10(3) when the CNT was used, and 0.2x10(-1) when the SCC was used and 0.7x10(-1) when the Hale's scaled theory was used. The average values represent all the alcohols used at the same reduced temperatures. The average ratio was about the same throughout the temperature range, although J(exp)/J(the) calculated with the Hale's scaled theory increased slightly with increasing temperature. The saturation ratio dependency was predicted closest to experiment with the classical nucleation theory. The nucleation rates were compared to those found in the literature. The measurements were in reasonable agreement with each other. The molecular content of critical alcohol clusters was between 35 and 80 molecules. At a fixed reduced temperature, the number of molecules in a critical cluster decreased as a function of alcohol carbon chain length. The number of molecules in critical clusters was compared to those predicted by the Kelvin equation. The theory predicted the critical cluster sizes well.     

Kinetic theory of heterogeneous nucleation; effect of nonuniform density in the nuclei

The heterogeneous nucleation of a liquid from a vapor in contact with a planar solid surface or a solid surface with cavities is examined on the basis of the kinetic theory of nucleation developed by Nowakowski and Ruckenstein [J. Phys. Chem. 96 (1992) 2313] which is extended to nonuniform fluid density distribution (FDD) in the nucleus. The latter is determined under the assumption that at each moment the FDD in the nucleus is provided by the density functional theory (DFT) for a nanodrop. As a result of this assumption, the theory does not require to consider that the contact angle which the nucleus makes with the solid surface and the density of the nucleus are independent parameters since they are provided by the DFT. For all considered cases, the nucleation rate is higher in the cavities than on a planar surface and increases with increasing strength of the fluid-solid interactions and decreasing cavity radius. The difference is small at high supersaturations (small critical nuclei), but becomes larger at low supersaturations when the critical nucleus has a size comparable with the size of the cavity. The nonuniformity of the FDD in the nucleus decreases the nucleation rate when compared to the uniform FDD.     

Formation mechanism of critical nucleus during nucleation process of liquid metal sodium

To deeply understand the formation mechanism of a critical nucleus during the nucleation process of liquid metal sodium, a system consisting of 10 000 Na atoms has been simulated by using molecular dynamics method. The evolutions of nuclei are traced directly, adopting the cluster-type index method. It is found that the energies of clusters and their geometrical constraints interplay to form the favorable microstructures during the nucleation process. The nucleus can be formed through many different pathways, and the critical size of the nucleus would be different for each pathway. It is also found that the critical nucleus is nonspherical and may include some metastable structures. Furthermore, the size of the cluster and its internal structure both play a crucial role in determining whether it is a critical nucleus, and this is in agreement with the simulations by computing the free energy of the Lennard-Jones system [D. Moroni, P. R. ten Wolde, and P. G. Bolhuis, Phys. Rev. Lett. 94, 235703 (2005)].     

Molecular dynamics simulations of nucleation from vapor to solid composed of Lennard-Jones molecules

We performed molecular dynamics (MD) simulations of nucleation from vapor at temperatures below the triple point for systems consisting of 10(4)-10(5) Lennard-Jones (L-J) type molecules in order to test nucleation theories at relatively low temperatures. Simulations are performed for a wide range of initial supersaturation ratio (S(0) ? 10-10(8)) and temperature (kT = 0.2-0.6ε), where ε and k are the depth of the L-J potential and the Boltzmann constant, respectively. Clusters are nucleated as supercooled liquid droplets because of their small size. Crystallization of the supercooled liquid nuclei is observed after their growth slows. The classical nucleation theory (CNT) significantly underestimates the nucleation rates (or the number density of critical clusters) in the low-T region. The semi-phenomenological (SP) model, which corrects the CNT prediction of the formation energy of clusters using the second virial coefficient of a vapor, reproduces the nucleation rate and the cluster size distributions with good accuracy in the low-T region, as well as in the higher-T cases considered in our previous study. The sticking probability of vapor molecules onto the clusters is also obtained in the present MD simulations. Using the obtained values of sticking probability in the SP model, we can further refine the accuracy of the SP model.     

Homogeneous nucleation of n-nonane and n-propanol mixtures: a comparison of classical nucleation theory and experiments

The homogeneous nucleation rates for n-nonane-n-propanol vapor mixtures have been calculated as a function of vapor-phase activities at 230 K using the classical nucleation theory (CNT) with both rigorous and approximate kinetic prefactors and compared to previously reported experimental data. The predicted nucleation rates resemble qualitatively the experimental results for low n-nonane gas phase activity. On the high nonane activity side the theoretical nucleation rates are about three orders of magnitude lower than the experimental data when using the CNT with the approximate kinetics. The accurate kinetics improves the situation by reducing the difference between theory and experiments to two orders of magnitude. Besides the nucleation rate comparison and the experimental and predicted onset activities, the critical cluster composition is presented. The total number of molecules is approximated by CNT with reasonable accuracy. Overall, the classical nucleation theory with rigorous kinetic prefactor seems to perform better. The thermodynamic parameters needed to calculate the nucleation rates are revised extensively. Up-to-date estimates of liquid phase activities using universal functional activity coefficient Dortmund method are presented together with the experimental values of surface tensions obtained in the present study.     

Sensitivity of nucleation phenomena on range of interaction potential

Theoretical and computational investigations of nucleation have been plagued by the sensitivity of the phase diagram to the range of the interaction potential. As the surface tension depends strongly on the range of interaction potential and as the classical nucleation theory (CNT) predicts the free energy barrier to be directly proportional to the cube of the surface tension, one expects a strong sensitivity of nucleation barrier to the range of the potential; however, CNT leaves many aspects unexplored. We find for gas-liquid nucleation in Lennard-Jones system that on increasing the range of interaction the kinetic spinodal (KS) (where the mechanism of nucleation changes from activated to barrierless) shifts deeper into the metastable region. Therefore the system remains metastable for larger value of supersaturation and this allows one to explore the high metastable region without encountering the KS. On increasing the range of interaction, both the critical cluster size and pre-critical minima in the free energy surface of kth largest cluster, at respective kinetic spinodals, shift towards smaller cluster size. In order to separate surface tension contribution to the increase in the barrier from other non-trivial factors, we introduce a new scaling form for surface tension and use it to capture both the temperature and the interaction range dependence of surface tension. Surprisingly, we find only a weak non-trivial contribution from other factors to the free energy barrier of nucleation.     

Molecular dynamics investigation of the transient regime in the freezing of salt clusters

Molecular dynamics (MD) investigations of the freezing of supercooled liquids can identify nuclei far smaller than can be detected in laboratory experiments, to date, and consequently can provide information about nucleation so far inaccessible to experiment. In a recent MD study of the freezing of clusters of SeF6, a new method of recording nucleation events was introduced. It involved the observation of times of first appearance of nuclei of the size of n. An advantage of the new approach is that it provides information about the size of the critical nucleus. For nuclei smaller than the critical size, it also avoids the overshoots of nucleation rates that precluded the application of the Miloshev-Wu method in the subcritical region. Kinetic information in the transient regime can be characterized by three parameters, the time lag, the reduced moment, and the steady state nucleation rate. To get some idea of how general the new approach is, a very different system was investigated, that of clusters of NaCl. Two different fitting functions were used to analyze the results. The first one adopted the log-normal probability distribution of Wu. The second function was a modification of Shneidman's analytical solution appropriate for large nuclei. The second function gave a rather good account of MD data for all nuclear sizes and temperatures and gave more stable results in the subcritical region. Several inferences of the sizes of critical nuclei were made. Applying the criterion for n* based on the Zeldovich solution of the Becker-D?ring equations, we estimated the critical nucleus sizes to be 14, 18, and 24 ions for quench temperatures of 640, 690, and 740 K, respectively. Even though the interionic interactions initiating nucleation in salt are very different from the van der Waals interactions in clusters of SeF6, the characteristic aspects of the transient regimes of the two systems were quite similar.