Item

Abstract

Observing nonclassical nucleation pathways remains challenging in simulations of complex materials with technological interests. This is because it requires very accurate force fields that can capture the whole complexity of their underlying interatomic interactions and an advanced structural analysis able to discriminate between competing crystalline phases. Here, we first report the construction and particularly thorough validation of a machine learning force field for zinc oxide interactions using the Physical LassoLars Interaction Potentials approach which allows us to be predictive even for high-temperature dynamical systems such as ZnO melt. Then, we carried out several types of crystallization simulations and followed the formation of ZnO crystals with atomistic precision. Our results, which were analyzed using a data-driven approach based on bond order parameters, demonstrate the presence of both prenucleation clusters and two-step nucleation scenarios, thus retrieving seminal predictions of nonclassical nucleation pathways made on much simpler models. Dedicated calculations of high temperature ZnO free energy within a newly developed automated nonequilibrium thermodynamic integration method revealed the existence of a thermodynamic bias for the predicted nonclassical nucleation scenarios.