Metadynamics for Vacancy Dynamics in Crystals

We propose a metadynamics-based (MetaD-based) approach for constructing the free energy surface (FES) of vacancy dynamics in crystals. In this approach, the vacancy FES can be constructed without explicitly defining a unique vacancy coordinate or introducing a set of parameters that strongly govern the FES topology, enabled by parallel bias MetaD with partitioned families (PB MetaDPF). In addition, the proposed approach is made more efficient and effective through a multi-hill strategy that exploits crystallographic symmetry. We demonstrate the validity of the proposed approach through applications to self-diffusion and impurity diffusion via monovacancies and divacancies in metallic and ionic crystals.