A statistical analysis of the first stages of freezing and melting of Lennard-Jones particles: Number and size distributions of transient nuclei

TL;DR

This paper analyzes the early stages of freezing and melting in Lennard-Jones particles through molecular simulations, focusing on transient nuclei's number and size distributions and their evolution across phases.

Contribution

It provides a detailed characterization of proto-nuclei distributions in the initial nucleation stages, including their existence in stable phases and their continuous evolution across phases.

Findings

Proto-nuclei exist even in stable phases.

Number and size distributions evolve continuously across phases.

Size distributions relate to free energy barrier calculations.

Abstract

The freezing/melting transition is at the heart of many natural and industrial processes. In the classical picture, the transition proceeds via the nucleation of the new phase, which has to overcome a barrier associated to the free energy cost of the growing nucleus. The total nucleation rate is also influenced by a kinetic factor which somehow depends on the number of attempts to create a nucleus, that translates into a significant density of proto-nuclei in the system. These transient tiny nuclei are not accessible to experiments, but they can be observed in molecular simulations, and their number and size distributions can be acquired and analysed. The number distributions are carefully characterized as a function of the system size, showing the expected behavior, with limited spurious effects due to the finite simulation box. It is also shown that the proto-nuclei do exist even in the stable phase, in agreement with the fact that the (unfavorable) volume contribution to their free energy is negligible in the first stages of nucleation. Moreover, the number and size distributions evolve continuously between the stable and the metastable phases, in particular when crossing the coexistence temperature. The size distributions associated to \textit{any} nucleus and to the \textit{largest} one have also been calculated, and their relationship recently established for bubbles in a liquid [J. Puibasset, J. Chem. Phys. 157, 191102 (2022)] has been shown to apply here. This is an important relation for free energy barrier calculations with biased molecular simulations.