The challenges of hydrogen and metals.

The Royal Society Scientific Discussion Meeting 'The challenges of hydrogen and metals' was held in Carlton House Terrace, London, UK, on 16-18 January 2017. This is the introductory article to the discussion meeting issue which includes contributed papers and seven discussion papers. Here, we introduce the motivation to hold the Meeting and give a brief overview of the contents. We conclude with acknowledgements.This article is part of the themed issue 'The challenges of hydrogen and metals'.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules.

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water.

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice--leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point.

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water.

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster.

A tight binding model for water.

We demonstrate for the first time a tight binding model for water incorporating polarizable oxygen atoms. A novel aspect is that we adopt a "ground up" approach in that properties of the monomer and dimer only are fitted. Subsequently we make predictions of the structure and properties of hexamer clusters, ice-XI and liquid water. A particular feature, missing in current tight binding and semiempirical hamiltonians, is that we reproduce the almost two-fold increase in molecular dipole moment as clusters are built up toward the limit of bulk liquid. We concentrate on properties of liquid water, particularly dielectric constant and self diffusion coefficient, which are very well rendered in comparison with experiment. Finally we comment on the question of the contrasting densities of water and ice which is central to an understanding of the subtleties of the hydrogen bond.

Bandstructure approach to near edge structure.

We review the current state of the art in EELS fingerprinting by computer simulation, focusing on the bandstructure approach to the problem. Currently calculations are made using a one electron theory, but we describe in principle the way to go beyond this to include final state effects. We include these effects within the one electron framework using the Slater transition state formula and assess the errors involved. Two examples are then given which illustrate the use of the one electron approximation within density functional theory. Our approach is to combine predicted atomic structure with predicted electronic structure to assist in fingerprinting of complex crystal structures.

Energization of the transport systems for arabinose and comparison with galactose transport in Escherichia coli.

1. Strains of Escherichia coli were obtained containing either the AraE or the AraF transport system for arabinose. AraE+,AraF- strains effected energized accumulation and displayed an arabinose-evoked alkaline pH change indicative of arabinose-H+ symport. In contrast, AraE-,AraF+ strains accumulated arabinose but did not display H+ symport. 2. The ability of different sugars and their derivatives to elicit sugar-H+ symport in AraE+ strains was examined. Only L-arabinose and D-fucose were good substrates, and arabinose was the only inducer. 3. Membrane vesicles prepared from an AraE+,AraF+ strain accumulated the sugar, energized most efficiently by the respiratory substrates ascorbate + phenazine methosulphate. Addition of arabinose or fucose to an anaerobic suspension of membrane vesicles caused an alkaline pH change indicative or sugar-H+ symport on the membrane-bound transport system. 4. Kinetic studies and the effects of arsenate and uncoupling agents in intact cells and membrane vesicles gave further evidence that AraE is a low-affinity membrane-bound sugar-H+ symport system and that AraF is a binding-protein-dependent high-affinity system that does not require a transmembrane protonmotive force for energization. 5. The interpretation of these results is that arabinose transport into E. coli is energized by an electrochemical gradient of protons (AraE system) or by phosphate bond energy (AraF system). 6. In batch cultures the rates of growth and carbon cell yields on arabinose were lower in AraE-,AraF+ strains than in AraE+,AraF- or AraE+,AraF+ strains. The AraF system was more susceptible to catabolite repression than was the AraE system. 7. The properties of the two transport systems for arabinose are compared with those of the genetically and biochemically distinct transport systems for galactose, GalP and MglP. It appears that AraE is analogous to GalP, and AraF to MglP.