Homogeneous nucleation and growth in simple fluids. II. Scaling behavior, instabilities, and the (n,v) order parameter

Calculation of critical nucleation rates by the persistent embryo method: application to quasi hard sphere models

We study the crystal nucleation of the Weeks-Chandler-Andersen (WCA) model, using the recently introduced persistent embryo method (PEM). The method provides detailed characterization of pre-critical, critical and post-critical nuclei, as well as nucleation rates that compare favorably with those obtained using other methods (umbrella sampling, forward flux sampling or seeding). We further map our results to a hard sphere model allowing comparison with other existing predictions. Implications for experiments are also discussed.

Understanding the sensitivity of nucleation free energies: The role of supersaturation and temperature

Simulations are used to investigate the vapor-to-liquid nucleation of water for several different force fields at various sets of physical conditions. The nucleation free energy barrier is found to be extremely sensitive to the force field at the same absolute conditions. However, when the results are compared at the same supersaturation and reduced temperature or the same metastability parameter and reduced temperature, then the differences in the nucleation free energies of the different models are dramatically reduced. This finding suggests that comparisons of experimental data and computational predictions are most meaningful at the same relative conditions and emphasizes the importance of knowing the phase diagram of a given computational model, but such information is usually not available for models where the interaction energy is determined directly from electronic structure calculations.

Toward a molecular theory of homogeneous bubble nucleation: II. Calculation of the number density of critical nuclei and the rate of nucleation

In the present paper, we develop a method to calculate the rate of homogeneous bubble nucleation within a superheated L-J liquid based on the (n,v) equilibrium embryo free energy surface introduced in the first paper (DOI: 10.1021/jp404149n). We express the nucleation rate as the product of the concentration of critical nuclei within the metastable liquid phase and the relevant forward rate coefficient. We calculate the forward rate coefficient of the critical nuclei from their average lifetime as determined from MD simulations of a large number of embryo trajectories initiated from the transitional region of the metastable liquid configuration space. Therefore, the proposed rate coefficient does not rely on any predefined reaction coordinate. In our model, the critical nuclei belong to the region of the configuration space where the committor probability is about one-half, guaranteeing the dynamical relevance of the proposed embryos. One novel characteristic of our approach is that we define a limit for the configuration space of the equilibrium metastable phase and do not include the configurations that have zero committor probability in the nucleation free energy surface. Furthermore, in order to take into account the transitional degrees of freedom of the critical nuclei, we develop a simulation-based approach for rigorously mapping the free energy of the (n,v) equilibrium embryos to the concentration of the critical nuclei within the bulk metastable liquid phase.

Homogeneous nucleation and growth in simple fluids. I. Fundamental issues and free energy surfaces of bubble and droplet formation

The free energy of forming a droplet and a bubble with a given particle number n and volume v within the pure-component Lennard-Jones supercooled vapor and superheated liquid, respectively, are further explored using density-functional theory. Similar to what was found previously [M. J. Uline and D. S. Corti, Phys. Rev. Lett. 99, 076102 (2007); M. J. Uline and D. S. Corti, J. Chem. Phys. 129, 234507 (2008)], the limits of stability again appear within both free energy surfaces evaluated at two other metastability conditions, one closer to the binodal and one closer to the spinodal. Furthermore, an ad hoc bond connectivity criterion is also applied in an attempt, however approximately, to eliminate certain configurational redundancies that arise from the chosen droplet and bubble definitions. What results are free energy surfaces describing the formation of equilibrium embryos that should be an improved representation of the fluctuations that are relevant to those nonequilibrium embryos seen in an actual nucleation event. Finally, we discuss in some detail the use of the (n,v) reaction coordinate within the framework of an equilibrium-based theory and its relation to other descriptions of nucleation.