Abstract
The knowledge of diffusion mechanisms in materials is crucial for predicting their high-temperature performance and stability, yet accurately capturing the underlying physics like thermal effects remains challenging. In particular, the origin of the experimentally observed non-Arrhenius diffusion behavior has remained elusive, largely due to the lack of effective computational tools. Here we propose an efficient ab initio framework to compute the Gibbs energy of the transition state in vacancy-mediated diffusion including the relevant thermal excitations at the density-functional-theory level. With the aid of a bespoke machine-learning interatomic potential, the temperature-dependent vacancy formation and migration Gibbs energies of the prototype system body-centered cubic (BCC) tungsten are shown to be strongly affected by anharmonicity. This finding explains the physical origin of the experimentally observed non-Arrhenius behavior of tungsten self-diffusion. A remarkable agreement between the calculated and experimental temperature-dependent self-diffusivity and, in particular, its curvature is revealed. The proposed computational framework is robust and broadly applicable, as evidenced by first tests for a hexagonal close-packed (HCP) multicomponent high-entropy alloy. The successful applications underscore the attainability of an accurate ab initio diffusion database.