Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Gibbs free energy of homogeneous nucleation: from atomistic nuclei to the planar limit
66
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper we discuss how the information contained in atomistic simulations of homogeneous nucleation should be used when fitting the parameters in macroscopic nucleation models. We show how the number of solid and liquid atoms in such simulations can be determined unambiguously by using a Gibbs dividing surface and how the free energy as a function of the number of solid atoms in the nucleus can thus be extracted. We then show that the parameters of a model based on classical nucleation theory can be fit using the information contained in these free-energy profiles but that the parameters in such models are highly correlated. This correlation is unfortunate as it ensures that small errors in the computed free energy surface can give rise to large errors in the extrapolated properties of the fitted model. To resolve this problem we thus propose a method for fitting macroscopic nucleation models that uses simulations of planar interfaces and simulations of three-dimensional nuclei in tandem. We show that when the parameters of the macroscopic model are fitted in this way the numerical errors for the final fitted model are smaller and that the extrapolated predictions for large nuclei are thus more reliable.
rate research
Read More
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 to extract the surface free energy of the nucleus. Neither the knowledge of the (in general non-spherical) 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 $Delta F^*=(4pi/3)^{1/3}bar{gamma}(V^*)^{2/3}+const.$ with the nucleus volume $V^*$, and the effective surface tension $bar{gamma}$ that accounts implicitly for the nonspherical shape can be precisely estimated.
Macroscopic models of nucleation provide powerful tools for understanding activated phase transition processes. These models do not provide atomistic insights and can thus sometime lack material-specific descriptions. Here we provide a comprehensive framework for constructing a continuum picture from an atomistic simulation of homogeneous nucleation. We use this framework to determine the shape of the equilibrium solid nucleus that forms inside bulk liquid for a Lennard-Jones potential. From this shape, we then extract the anisotropy of the solid-liquid interfacial free energy, by performing a reverse Wulff construction in the space of spherical harmonic expansions. We find that the shape of the nucleus is nearly spherical and that its anisotropy can be perfectly described using classical models.
The homogeneous and heterogeneous nucleation of a Lennard-Jones liquid is investigated using the umbrella sampling method. The free energy cost of forming a nucleating droplet is determined as a function of the quench depth, and the saddle point nature of the droplets is verified using an intervention technique. The structure and symmetry of the nucleating droplets is found for a range of temperatures. We find that for deep quenches the nucleating droplets become more anisotropic and diffuse with no well defined core or surface. The environment of the nucleating droplets form randomly stacked hexagonal planes. This behavior is consistent with a spinodal nucleation interpretation. We also find that the free energy barrier for heterogeneous nucleation is a minimum when the lattice spacing of the impurity equals the lattice spacing of the equilibrium crystalline phase. If the lattice spacing of the impurity is different, the crystal grows into the bulk instead of wetting the impurity.
Molecular dynamics simulation is used to study the time-scales involved in the homogeneous melting of a superheated crystal. The interaction model used is an embedded-atom model for Fe developed in previous work, and the melting process is simulated in the microcanonical $(N, V, E)$ ensemble. We study periodically repeated systems containing from 96 to 7776 atoms, and the initial system is always the perfect crystal without free surfaces or other defects. For each chosen total energy $E$ and number of atoms $N$, we perform several hundred statistically independent simulations, with each simulation lasting for between 500 ps and 10 ns, in order to gather statistics for the waiting time $tau_{rm w}$ before melting occurs. We find that the probability distribution of $tau_{rm w}$ is roughly exponential, and that the mean value $<tau_{rm w} >$ depends strongly on the excess of the initial steady temperature of the crystal above the superheating limit identified by other researchers. The mean $<tau_{rm w}>$ also depends strongly on system size in a way that we have quantified. For very small systems of $sim 100$ atoms, we observe a persistent alternation between the solid and liquid states, and we explain why this happens. Our results allow us to draw conclusions about the reliability of the recently proposed Z method for determining the melting properties of simulated materials, and to suggest ways of correcting for the errors of the method.
Estimating the homogeneous ice nucleation rate from undercooled liquid water is at the same time crucial for understanding many important physical phenomena and technological applications, and challenging for both experiments and theory. From a theoretical point of view, difficulties arise due to the long time scales required, as well as the numerous nucleation pathways involved to form ice nuclei with different stacking disorders. We computed the homogeneous ice nucleation rate at a physically relevant undercooling for a single-site water model, taking into account the diffuse nature of ice-water interfaces, stacking disorders in ice nuclei, and the addition rate of particles to the critical nucleus.We disentangled and investigated the relative importance of all the terms, including interfacial free energy, entropic contributions and the kinetic prefactor, that contribute to the overall nucleation rate.There has been a long-standing discrepancy for the predicted homogeneous ice nucleation rates, and our estimate is faster by 9 orders of magnitude compared with previous literature values. Breaking down the problem into segments and considering each term carefully can help us understand where the discrepancy may come from and how to systematically improve the existing computational methods.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
