All-electron study of the electronic properties of quartz with Al substitutional impurity

Modeling doped and defective oxides in catalysis with density functional theory methods: room for improvements

Due to the well-known problem of the self-interaction, standard density functional theory (DFT) methods tend to produce delocalized holes and electrons in defective oxide materials even when there is ample experimental evidence of a strong localization. For late transition metal compounds or rare earth oxides, this results in the incorrect description of the electronic structure of the system (e.g., magnetic insulators are predicted to be metallic). Practical ways to correct this deficiency are based on the use of hybrid functionals or of the DFT+U approach. In this way, most of the limitations related to the self-interaction are removed, and the electronic structure is properly described. What is less clear is to what extent hybrid functionals, DFT+U approaches, or standard DFT functionals can properly describe the strength of the chemical bonds at the surface of an oxide. This is a crucial question if one is interested in the catalytic properties of oxide surfaces. Oxidation reactions often involve oxygen detachment from the surface and incorporation into an organic substrate. Oxides are doped with heteroatoms to create defects and facilitate oxygen removal from the surface, with formation of oxygen vacancies. Do standard DFT calculations provide a good binding energy of the missing oxygen despite the failure in giving the right electronic structure? Can hybrid functionals or the DFT+U approach provide a simple yet reliable way to get accurate reaction enthalpies and energy barriers? In this essay, we discuss these problems by analyzing some case histories and the relatively scarce data existing in the literature. The conclusion is that while modern electronic structure methods accurately reproduce and predict a wide range of electronic, optical, and magnetic properties of oxides, the description of the strength of chemical bonds still needs considerable improvements.

Hole localization in [AlO4]0 defects in silica materials

First-principles calculations based on cluster models have been performed to investigate the ground state and the optically excited states of the [AlO(4)](0) hole in alpha-quartz and in the siliceous zeolite ZSM-5. The structure and spectroscopic properties of this defect have been studied using the recently developed Becke88-Becke95 one-parameter model for kinetics (BB1K) functional of Zhao et al., [J. Phys. Chem. A 108, 2715 (2004)]. Our results show that the BB1K method is significantly more reliable and more accurate than the standard density-functional theory (DFT) functionals at reproducing the localized spin density on one oxygen atom and the hyperfine coupling constants associated with the hole. Furthermore, we find that the BB1K results are in close agreement with experiments, and with the self-interaction-free unrestricted Hartree-Fock (UHF) and unrestricted second-order Møller-Plesset perturbation theory (UMP2) calculations. For the first time, we present results of the ground-state paramagnetic properties of the Al defect in ZSM-5. Similar to the theoretical work for defective alpha-quartz, we find that the BB1K, UHF, UHFLee-Yang-Parr, and UMP2 calculations show a localized hole on one oxygen neighboring the Al, while even the best to date thermochemically derived hybrid generalized gradient approximation density-functional, B97-2, predicts a different model where the hole is distributed over two oxygen. We have further considered the optical transitions of the [AlO(4)](0) center in alpha-quartz and ZSM-5. In both systems, our BB1K time-dependent density-functional theory (TDDFT) and configuration interaction singles (CIS) calculations predict that the most likely transition involves electron transfer from the hole-bearing oxygen to other neighboring oxygen ions. This reinforces the experimental conclusions obtained for defective alpha-quartz. Notably, the two lowest, most dominant excitation energies calculated by BB1K-TDDFT (1.99 and 3.03 eV) show excellent agreement with experiment (1.96 and 2.85 eV [B. K. Meyer, J.M. Spaeth, and J.A. Weil, J. Phys. C: Solid State Phys. 17, L31 (1987)]) clearly outperforming the CIS method and other DFT calculations available in the literature.

Hole trapping at Al impurities in silica: a challenge for density functional theories

The atomic geometry and electronic structure around a neutral substitutional Al impurity in silica is investigated using either the unrestricted Hartree-Fock (UHF) approximation, or Beckes three-parameter hybrid functional (B3LYP). It is found that the B3LYP functional fails to describe the structural distortions around the Al impurity, while the UHF results are consistent with experimental information. We argue that the failure of the B3LYP functional is caused by the incomplete self-interaction cancellation usually present in density functional theories.