Machine learning metadynamics simulation of reconstructive phase transition

Simulating reconstructive phase transition requires an accurate description of potential energy surface (PES). Density-functional-theory (DFT) based molecular dynamics can achieve the desired accuracy but it is computationally unfeasible for large systems and/or long simulation times. Here we introduce an approach that combines the metadynamics simulation and machine learning representation of PES at the accuracy close to the DFT calculations, but with the computational cost several orders of magnitude less, and scaling with system size approximately linear. The high accuracy of the method is demonstrated in the simulation of pressure-induced $B4\text{\ensuremath{-}}B1$ phase transition in gallium nitride (GaN). The large-scale simulation using a 4096-atom simulation box reveals the phase transition with excellent detail, revealing different simulated transition paths under particular stress conditions. With well-trained machine learning potentials, this method can be easily applied to all types of systems for accurate scalable simulations of solid-solid reconstructive phase transition.