title: Hole localization in [AlO4]0 defects in silica materials.
abstract: 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.
authors.affiliation: Davy Faraday Research Laboratory, Royal Institution of Great Britain, 21 Albemarle Street, London W1S 4BS, UK. judy@ri.ac.uk
authors.forename: Judy
authors.initials: J
authors.lastname: To
authors.forename: Alexey A
authors.initials: AA
authors.lastname: Sokol
authors.forename: Samuel A
authors.initials: SA
authors.lastname: French
authors.forename: Nikolas
authors.initials: N
authors.lastname: Kaltsoyannis
authors.forename: C Richard A
authors.initials: CR
authors.lastname: Catlow