General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

Abstract

The O vacancy (Ov) formation energy, $E_{\text{Ov}}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_{\text{Ov}}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_{\text{Ov}}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_{\text{Ov}}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

Accepted for publication in Journal of Chemical Physics on 2022-03-10

References in corpus (9)

• Predicting Band Gaps with Hybrid Density Functionals
• First-principles DFT+GW study of oxygen vacancies in rutile TiO2
• Wave functions for quantum Monte Carlo calculations in solids: Orbitals from density functional theory with hybrid exchange-correlation functionals
• Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo
• Local embedding of Coupled Cluster theory into the Random Phase Approximation using plane-waves
• Embedded-Cluster Calculations in a Numeric Atomic Orbital Density-Functional Theory Framework
• Self-consistent hybrid functional calculations: Implications for structural, electronic, and optical properties of oxide semiconductors
• A periodic equation-of-motion coupled-cluster implementation applied to $F$-centers in alkaline earth oxides
• The color center singlet state of oxygen vacancies in TiO${}_2$