Thomas W. Keal

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.

DL-FIND: An Open-Source Geometry Optimizer for Atomistic Simulations

Geometry optimization, including searching for transition states, accounts for most of the CPU time spent in quantum chemistry, computational surface science, and solid-state physics, and also plays an important role in simulations employing classical force fields. We have implemented a geometry optimizer, called DL-FIND, to be included in atomistic simulation codes. It can optimize structures in Cartesian coordinates, redundant internal coordinates, hybrid-delocalized internal coordinates, and also functions of more variables independent of atomic structures. The implementation of the optimization algorithms is independent of the coordinate transformation used. Steepest descent, conjugate gradient, quasi-Newton, and L-BFGS algorithms as well as damped molecular dynamics are available as minimization methods. The partitioned rational function optimization algorithm, a modified version of the dimer method and the nudged elastic band approach provide capabilities for transition-state search. Penalty function, gradient projection, and Lagrange-Newton methods are implemented for conical intersection optimizations. Various stochastic search methods, including a genetic algorithm, are available for global or local minimization and can be run as parallel algorithms. The code is released under the open-source GNU LGPL license. Some selected applications of DL-FIND are surveyed.

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.

Activation of Carbon Dioxide over Zinc Oxide by Localised Electrons

The unique mechanism of activation of carbon dioxide over zinc oxide is unravelled using advanced quantum mechanical methods. The key process is the CO(2) chemisorption catalysed by a highly localized electron carrier trapped at a vacant oxygen interstitial surface site. At the top of the reaction barrier CO(2) pulls the electron from the vacancy and thus becomes active.

Characterization of hydrogen dissociation over aluminium-doped zinc oxide using an efficient massively parallel framework for QM/MM calculations

A task-farm parallelization framework has been implemented in the ChemShell computational chemistry environment to provide a facility for parallelizing common chemical calculations, including finite-difference Hessian evaluation, the nudged elastic band method for reaction path optimization, and population-based methods for global optimization. The optimization methods are provided by a parallel interface to the DL-FIND optimization library. As ChemShell can already exploit parallel external programs for energy and gradient evaluations, the new methods result in a two-level approach to parallelization that gives significantly improved performance for massively parallel calculations. For typical systems, speed-up factors of five to eight times have been observed compared with non-task-farmed calculations. The task-farming version of ChemShell has been used to study the heterolytic dissociation of a hydrogen molecule over a polar oxygen-terminated surface of aluminium-doped zinc oxide using an embedded cluster hybrid QM/MM approach. We calculate a 42 kcal mol−1 heat of reaction and a 30 kcal mol−1 activation energy, which is equivalent to a high backward reaction barrier of 72 kcal mol−1 per H2 molecule, in close agreement with temperature programmed desorption experiments. The dissociation path includes a stable intermediate comprising a hydride ion in an oxygen vacancy and physisorbed atomic hydrogen.

Active-site protein dynamics and solvent accessibility in native Achromobacter cycloclastes copper nitrite reductase

Multiple structures obtained from one crystal of copper nitrite reductase at elevated cryogenic temperature, together with molecular-dynamics simulations, reveal catalyically important protein and solvent dynamics at the active site.

Prediction of multiband luminescence due to the gallium vacancy–oxygen defect complex in GaN

Oxygen is the most common unintentional impurity found in GaN. We study the interaction between substitutional oxygen (ON) and the gallium vacancy (VGa) to form a point defect complex in GaN. The formation energy of the gallium vacancy is largely reduced in n-type GaN by complexing with oxygen, while thermodynamic and optical transition levels remain within the bandgap. We study the spectroscopy of this complex using a hybrid quantum-mechanical molecular-mechanical embedded-cluster approach. We reveal how a single defect center can be responsible for multiband luminescence, including possible contributions to the ubiquitous yellow luminescence signatures observed in n-type GaN, owing to the coexistence of diffuse (extended) and compact (localized) holes.

MSOX crystallography and simulations to capture redox enzyme catalysis

Michael Alexander Hough1, Demet Kekilli1, Sam Horrell1, Kakali Sen1, Chin Yong2, Thomas W.K. Keal2, Svetlana V Antonyuk3, Robert R. Eady3, S Samar Hasnain3, Richard W. Strange1 1School Of Biological Sciences, Colchester, United Kingdom, 2Scientific Computing Department, STFC Daresbury Laboratory, Warrington, United Kingdom, 3Institute of Integrative Biology, University of Liverpool, Liverpool, United Kingdom E-mail: mahough@essex.ac.uk

Modelling the chemistry of Mn-doped MgO for bulk and (100) surfaces

We have investigated the energetic properties of Mn-doped MgO bulk and (100) surfaces using a QM/MM embedding computational method, calculating the formation energy for doped systems, as well as for surface defects, and the subsequent effect on chemical reactivity. Low-concentration Mn doping is endothermic for isovalent species in the bulk but exothermic for higher oxidation states under p-type conditions, and compensated by electrons going to the Fermi level rather than cation vacancies. The highest occupied dopant Mn 3d states are positioned in the MgO band gap, about 4.2 eV below the vacuum level. Surface Mn-doping is more favourable than subsurface doping, and marginally exothermic on a (100) surface at high O2 pressures. For both types of isovalent Mn-doped (100) surfaces, the formation energy for catalytically important oxygen defects is less than for pristine MgO, with F0 and F2+-centres favoured in n- and p-type conditions, respectively. In addition, F+-centres are stabilised by favourable exchange coupling between the Mn 3d states and the vacancy-localised electrons, as verified through calculation of the vertical ionisation potential. The adsorption of CO2 on to the pristine and defective (100) surface is used as a probe of chemical reactivity, with isovalent subsurface Mn dopants mildly affecting reactivity, whereas isovalent surface-positioned Mn strongly alters the chemical interactions between the substrate and adsorbate. The differing chemical reactivity, when compared to pristine MgO, justifies further detailed investigations for more varied oxidation states and dopant species.