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.

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

References in corpus (18)

• Predicting Band Gaps with Hybrid Density Functionals
• A direct view at excess electrons in TiO${}_2$ rutile and anatase
• Self-Energy and Excitonic Effects in the Electronic and Optical Properties of TiO2 Crystalline Phases
• On the oxidation state of titanium in titanium dioxide
• First-principles DFT+GW study of oxygen vacancies in rutile TiO2
• New generation of effective core potentials from correlated calculations: 3d transition metal series
• Wave functions for quantum Monte Carlo calculations in solids: Orbitals from density functional theory with hybrid exchange-correlation functionals
• Point-Defect Optical Transitions and Thermal Ionization Energies from Quantum Monte Carlo Methods: Application to F-center Defect in MgO
• 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
• A comparison between quantum chemistry and quantum Monte Carlo techniques for the adsorption of water on the (001) LiH surface
• 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
• Straightforward and accurate automatic auxiliary basis set generation for molecular calculations with atomic orbital basis sets
• Diffusion quantum Monte Carlo calculations of SrFeO${}_3$ and LaFeO${}_3$
• First-principles embedded-cluster calculations of the neutral and charged oxygen vacancy at the rutile TiO${}_2$(110) surface
• The color center singlet state of oxygen vacancies in TiO${}_2$