O interstitial diffusion in : Direct approach via master diffusion equations
Phys. Rev. B 112, 024112 – Published 28 July, 2025
DOI: https://doi.org/10.1103/tx7k-l4bq
Abstract
Monoclinic , a promising wide-band-gap semiconducting material, exhibits complex, anisotropic diffusional characteristics and mass transport behavior as a results of its low-symmetry crystal structure. From first-principles calculations combined with master diffusion equations, we determine three-dimensional diffusion tensors for neutral () and charged oxygen interstitials (). Systematic exploration of the configurational space identifies stable configurations in these two dominant charge states and their corresponding formation energies. By connecting every pair of low-energy configurations considering both interstitial or interstitialcy hops, we construct three-dimensional diffusion networks and evaluate hopping barriers of all transition pathways in networks. Combining the collection of (i) defect configurations and their formation energies and (ii) the hopping barriers that link them, we construct and solve the master diffusion equations for and separately through the Onsager approach, resulting in the respective three-dimensional diffusion tensors and . Both and present the fastest diffusion along the axis, demonstrating significant anisotropy. The predicted self-diffusivities along [100] and [] fall within the range of previously reported values from isotopically labeled oxygen tracer experiments, though some quantitative differences remain.