Toward a Consistent Prediction of Defect Chemistry in CeO₂

Abstract

Polarizable shell-model potentials are widely used for atomic-scale modeling of charged defects in solids using the Mott–Littleton approach and hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) embedded-cluster techniques. However, at the pure MM level of theory, the calculated defect energetics may not satisfy the requirement of quantitative predictions and are limited to only certain charged states. Here, we proposed a novel interatomic potential development scheme that unifies the predictions of all relevant charged defects in CeO2 based on the Mott–Littleton approach and QM/MM electronic-structure calculations. The predicted formation energies of oxygen vacancies accompanied by different excess electron localization patterns at the MM level of theory reach the accuracy of density functional theory (DFT) calculations using hybrid functionals. The new potential also accurately reproduces a wide range of physical properties of CeO2, showing excellent agreement with experimental and other computational studies. These findings provide opportunities for accurate large-scale modeling of the partial reduction and nonstoichiometry in CeO2, as well as a prototype for developing robust interatomic potentials for other defective crystals.

Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments

“…Our recently proposed SKZCAM protocol is particularly suited for tackling this challenge. It is based upon the electrostatic embedding approach ,− (top panel of Figure b) and provides the design rubrics to generate a series of quantum clusters of systematically increasing size (middle panel of Figure b). We have shown previously and here (bottom panel of Figure b) that these clusters converge smoothly and rapidly to the bulk limit.…”

Level-Shifted Embedded Cluster Method for Modeling the Chemistry of Metal Oxides

“…The key idea of the LSEC method is the alignment of a target energy level with a reference result. This differs from the previously proposed approach of tailoring boundary ECPs, in which the parameters of the boundary ECPs, which usually take the form of a linear combination of three Gaussian functions, are adjusted to minimize the gradients on the relaxed ions and the spread of deep core levels in the energy spectrum. − In addition, note that there are no serious restrictions in terms of the target system, the reference model, the QM method, and the boundary ECPs. The target system can be metal oxides or other ionically bonded systems.…”

“…In (QM/MM) embedded cluster calculations, effective core potentials (ECPs) can be placed adjacent to the boundary of the cluster (i.e., QM region in QM/MM calculations), more specifically, at the positions of the cations surrounding the cluster, to better reproduce the electrostatic feature at the central region of the cluster and to prevent charge spreading . In most cases, a large-core ECP of the corresponding metal in the metal oxide is employed for boundary ECPs. − Alternatively, specially tailored ECPs can be used to improve the results. − …”

Multiscale QM/MM modelling of catalytic systems with ChemShell

“…217 The ionic QM/MM embedded cluster approach can be used to model both bulk defects and surface reactivity. In recent ChemShell studies it has been used to characterise native point defects in GaN, 218 to create optimised models of rutile TiO 2 surfaces, 219 to study oxygen vacancies in TiO 2 using DFT and highlevel wavefunction methods, 220 to improve interatomic potential descriptions of CeO 2 , 221 to investigate bulk and surface vacancies in MnO with potential catalytic applications to the transformation of CO 2 , 222 and to study the defect properties of Cu in ZnO, an important industrial catalyst for the synthesis of methanol. 223 Sainna et al used ChemShell to study the reaction of glycerol over a catalytic MgO surface to form methanol.…”