The seeding method: A test case for classical nucleation theory in small systems
Thomas Philippe, Yijian Wu, Aymane Graini
TL;DR
This paper tests classical nucleation theory (CNT) in small, confined systems by applying the NVT seeding MD approach to Lennard-Jones condensation. By inserting a pre-formed liquid seed into a vapor in a finite box, the authors identify two CNT states and extract stable cluster radii $R^*$ and vapor densities $\rho_v^*$, then compare these with CNT predictions using several thermodynamic models. CNT predictions based on the Johnson–Zollweg–Gubbins (JZG) equation of state show very good agreement with seeded MD across a range of temperatures; the KN equation of state also performs well at low temperature, while density functional theory (DFT) and ideal-gas approximations become unreliable at higher temperatures, though the latter remains practically useful for initializing seeds. Overall, the seeding method provides an efficient, robust framework to benchmark CNT in finite systems and offers guidance on which thermodynamic models are most reliable for predicting nucleation in confined environments.
Abstract
Molecular dynamics simulations are widely used to investigate nucleation in first-order phase transitions. Brute-force simulations, though popular, are limited to conditions of high metastability, where the critical cluster and the nucleation barrier are small. The seeding method has recently emerged as a powerful alternative for exploring lower supersaturation regimes by initiating simulations with a pre-formed nucleus. In confined systems (NVT ensemble), the seeded simulations are particularly effective for determining stable cluster properties and provide a stringent test case for classical nucleation theory (CNT). In this work, we perform NVT seeded simulations of Lennard-Jones condensation in small systems and compare them with CNT predictions based on several thermodynamic models, including equations of state, perturbation theory, and ideal gas approximation. We find that CNT accurately predicts stable cluster radii across a wide range of conditions. Notably, even the very simple ideal gas approximation proves useful for initializing seeded simulations. Furthermore, seeded simulation results correspond to the critical cluster radii of infinite systems: CNT predictions with good equations of state show very good agreement with simulations, while the perturbation theory and the ideal gas approximation perform well at low temperatures but deviate significantly at high temperatures.
