Monte Carlo simulations of crystal nuclei coexisting with the fluid phase in thermal equilibrium in finite volumes are presented and analyzed, for fluid densities from dense melts to the vapor. Generalizing the lever rule for two-phase coexistence in the canonical ensemble to finite volume, “measurements” of the nucleus volume together with the pressure and chemical potential of the surrounding fluid allows us to extract the surface free energy of the nucleus. Neither the knowledge of the (in general nonspherical) nucleus shape nor of the angle-dependent interface tension is required for this task. The feasibility of the approach is demonstrated for a variant of the Asakura-Oosawa model for colloid-polymer mixtures, which form face-centered cubic colloidal crystals. For a polymer to colloid size ratio of 0.15, the colloid packing fraction in the fluid phase can be varied from melt values to zero by the variation of an effective attractive potential between the colloids. It is found that the approximation of spherical crystal nuclei often underestimates actual nucleation barriers significantly. Nucleation barriers are found to scale as $\mathrm{\Delta }{F}^{*}={(4\pi /3)}^{1/3}\overline{\gamma }{\left(V^{*}\right)}^{2/3}+\text{const}$ with the nucleus volume $V^{*}$, and the effective surface tension $\overline{\gamma }$ that accounts implicitly for the nonspherical shape can be precisely estimated.

• Received 30 May 2017

DOI:https://doi.org/10.1103/PhysRevE.96.042609