Alexey A. Sokol

Band alignment of rutile and anatase TiO2

The most widely used oxide for photocatalytic applications owing to its low cost and high activity is TiO₂. The discovery of the photolysis of water on the surface of TiO₂ in 1972 launched four decades of intensive research into the underlying chemical and physical processes involved. Despite much collected evidence, a thoroughly convincing explanation of why mixed-phase samples of anatase and rutile outperform the individual polymorphs has remained elusive. One long-standing controversy is the energetic alignment of the band edges of the rutile and anatase polymorphs of TiO₂ (ref. ). We demonstrate, through a combination of state-of-the-art materials simulation techniques and X-ray photoemission experiments, that a type-II, staggered, band alignment of ~ 0.4 eV exists between anatase and rutile with anatase possessing the higher electron affinity, or work function. Our results help to explain the robust separation of photoexcited charge carriers between the two phases and highlight a route to improved photocatalysts.

QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis

Abstract We describe the work of the European project QUASI (Quantum Simulation in Industry, project EP25047) which has sought to develop a flexible QM/MM scheme and to apply it to a range of industrial problems. A number of QM/MM approaches were implemented within the computational chemistry scripting system, ChemShell, which provides the framework for deploying a variety of independent program packages. This software was applied in several large-scale QM/MM studies which addressed the catalytic decomposition of N 2 O by Cu-containing zeolites, the methanol synthesis reaction catalysed by Cu clusters supported on ZnO surfaces, and the modelling of enzyme structure and reactivity.

ChemShell—a modular software package for QM/MM simulations

ChemShell is a modular computational chemistry package with a particular focus on hybrid quantum mechanical/molecular mechanical (QM/MM) simulations. A core set of chemical data handling modules and scripted interfaces to a large number of quantum chemistry and molecular modeling packages underpin a flexible QM/MM scheme. ChemShell has been used in the study of small molecules, molecular crystals, biological macromolecules such as enzymes, framework materials including zeolites, ionic solids, and surfaces. We outline the range of QM/MM coupling schemes and supporting functions for system setup, geometry optimization, and transition‐state location (including those from the open‐source DL‐FIND optimization library). We discuss recently implemented features allowing a more efficient treatment of long range electrostatic interactions, X‐ray based quantum refinement of crystal structures, free energy methods, and excited‐state calculations. ChemShell has been ported to a range of parallel computers and we describe a number of options including parallel execution based on the message‐passing capabilities of the interfaced packages and task‐farming for applications in which a number of individual QM, MM, or QM/MM calculations can performed simultaneously. We exemplify each of the features by brief reference to published applications.

Computational approaches to the determination of active site structures and reaction mechanisms in heterogeneous catalysts

We apply quantum chemical methods to the study of active site structures and reaction mechanisms in mesoporous silica and metal oxide catalysts. Our approach is based on the use of both molecular cluster and embedded cluster (QM/MM) techniques, where the active site and molecular complex are described using density functional theory (DFT) and the embedding matrix simulated by shell model potentials. We consider three case studies: alkene epoxidation over the microporous TS-1 catalyst; methanol synthesis on ZnO and Cu/ZnO and C–H bond activation over Li-doped MgO.

Advances in computational studies of energy materials

We review recent developments and applications of computational modelling techniques in the field of materials for energy technologies including hydrogen production and storage, energy storage and conversion, and light absorption and emission. In addition, we present new work on an Sn2TiO4 photocatalyst containing an Sn(II) lone pair, new interatomic potential models for SrTiO3 and GaN, an exploration of defects in the kesterite/stannite-structured solar cell absorber Cu2ZnSnS4, and report details of the incorporation of hydrogen into Ag2O and Cu2O. Special attention is paid to the modelling of nanostructured systems, including ceria (CeO2, mixed CexOy and Ce2O3) and group 13 sesquioxides. We consider applications based on both interatomic potential and electronic structure methodologies; and we illustrate the increasingly quantitative and predictive nature of modelling in this field.

Point defects in ZnO

We have investigated intrinsic point defects in ZnO and extended this study to Li, Cu and Al impurity centres. Atomic and electronic structures as well as defect energies have been obtained for the main oxidation states of all defects using our embedded cluster hybrid quantum mechanical/molecular mechanical approach to the treatment of localised states in ionic solids. With these calculations we were able to explain the nature of a number of experimentally observed phenomena. We show that in zinc excess materials the energetics of zinc interstitial are very similar to those for oxygen vacancy formation. Our results also suggest assignments for a number of bands observed in photoluminescence and other spectroscopic studies of the material.

Zinc oxide: A case study in contemporary computational solid state chemistry

Computational techniques have been applied to study a broad range of chemical and physical properties of zinc oxide. Both interatomic‐potential and density functional theory methods are used to investigate structural, thermodynamic, surface, and defect properties. We survey the structures and energies of nano‐particulate zinc oxide.

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.

The effect of local environment on photoluminescence: A time-dependent density functional theory study of silanone groups on the surface of silica nanostructures

The optical absorption spectrum and lowest photoluminescence (PL) signal for silanone terminated silica nanostructures are studied using time-dependent density functional theory calculations on a range of realistic low energy silica nanocluster models. We show that the broad experimental absorption spectrum for silanone centers [V. A. Radtsig and I. M. Senchenya Russ. Chem. Bull. 45, 1849 (1996)] is most likely the result of a synergetic combination of inhomogeneous broadening, thermal broadening and the small energy differences between different excitations. We further demonstrate that upon relaxation of the excited state the excited electron and hole localize on only one silanone center, and that there is a clear and distinct link between the local environment of a silanone center and its absorption and PL spectra. Finally, we provide strong evidence that the silanone center does not have a double bond between the constituent silicon and oxygen atoms but rather can be probably more aptly described as the =Si(+)-O(-) charge-transfer species.

Assignment of the complex vibrational spectra of the hydrogenated ZnO polar surfaces using QM/MM embedding

Hydrogenated zinc oxide gives rise to complex vibrational spectra with many prominent features that remain unexplained. Our calculations have unambiguously shown that the presence of vacant oxygen and zinc interstitial surface sites is the only way to rationalize the observed spectra, notably the 1710 cm−1 zinc hydride stretching mode. The large number of such sites, which expose low-coordinated surface ions, are inherent at ionically reconstructed polar surfaces. The thermal stability of the sorbed hydrogen and the infrared activity of the resulting species are correlated with site coordination and coverage.