Homogeneous nucleation in a Poiseuille flow

Variation of the critical nucleus size and the corresponding work of formation with average flow speed at axisymmetric axis.

Effects of ensembles on methane hydrate nucleation kinetics

With an exception for its critical nucleus size, the more the thermostat and barostat are used, the faster the methane hydrate nucleates and the lower its crystallinity is.

Nucleation rate analysis of methane hydrate from molecular dynamics simulations

Clathrate hydrates are solid crystalline structures most commonly formed from solutions that have nucleated to form a mixed solid composed of water and gas. Understanding the mechanism of clathrate hydrate nucleation is essential to grasp the fundamental chemistry of these complex structures and their applications. Molecular dynamics (MD) simulation is an ideal method to study nucleation at the molecular level because the size of the critical nucleus and formation rate occur on the nano scale. Various analysis methods for nucleation have been developed through MD to analyze nucleation. In particular, the mean first-passage time (MFPT) and survival probability (SP) methods have proven to be effective in procuring the nucleation rate and critical nucleus size for monatomic systems. This study assesses the MFPT and SP methods, previously used for monatomic systems, when applied to analyzing clathrate hydrate nucleation. Because clathrate hydrate nucleation is relatively difficult to observe in MD simulations (due to its high free energy barrier), these methods have yet to be applied to clathrate hydrate systems. In this study, we have analyzed the nucleation rate and critical nucleus size of methane hydrate using MFPT and SP methods from data generated by MD simulations at 255 K and 50 MPa. MFPT was modified for clathrate hydrate from the original version by adding the maximum likelihood estimate and growth effect term. The nucleation rates calculated by MFPT and SP methods are within 5%, and the critical nucleus size estimated by the MFPT method was 50% higher, than values obtained through other more rigorous but computationally expensive estimates. These methods can also be extended to the analysis of other clathrate hydrates.

Trapping Effect on the Kinetic Critical Radius in Nucleation and Growth Processes

The critical nucleus size—above which nuclei grow, below dissolve—during diffusion controlled nucleation in binary solid-solid phase transformation process is calculated using kinetic Monte Carlo (KMC). If atomic jumps are slower in an A-rich nucleus than in the embedding B-rich matrix, the nucleus traps the A atoms approaching its surface. It doesn’t have enough time to eject A atoms before new ones arrive, even if it would be favourable thermodynamically. In this case the critical nucleus size can be even by an order of magnitude smaller than expected from equilibrium thermodynamics or without trapping. These results were published in [Z. Erdélyi et al., Acta Mater. 58 (2010) 5639]. In a recent paper M. Leitner [M. Leitner, Acta Mater. 60 (2012) 6709] has questioned our results based on the arguments that his simulations led to different results, but he could not point out the reason for the difference. In this paper we summarize our original results and on the basis of recent KMC and kinetic mean field (KMF) simulations we show that Leitner’s conclusions are not valid and we confirm again our original results.