Performance of density functionals for calculating barrier heights of chemical reactions relevant to astrophysics

The performance of 39 different LDA, GGA, meta-GGA, and hybrid density functionals has been evaluated, for calculating forward and reverse barrier heights of 10 gas-phase reactions involving hydrogen. The reactions are all relevant to astrochemistry. Special focus is put on the applicability of DFT for calculating the rates of corresponding surface hydrogenation reactions that are relevant to the chemistry of ice-coated interstellar grains. General trends in the performance of the density functionals for reactions involving H atoms, H-2, and OH are discussed. The OH+CO reaction is shown to be a very problematic case for DFT. The best overall performance is found for the hybrid density functionals, such as MPW1K, B97-1, B97-2, and B1B95. For several reactions, the HCTH GGA functionals and the VS98 and OLAP3 meta-GGA functionals also give results that are almost as good as those of the hybrid functionals.

Metal Oxide Nanoparticle Mediated Enhanced Raman Scattering and Its Use in Direct Monitoring of Interfacial Chemical Reactions

Metal oxide nanoparticles (MONPs) have widespread usage across many disciplines, but monitoring molecular processes at their surfaces in situ has not been possible. Here we demonstrate that MONPs give highly enhanced (X10(4)) Raman scattering signals from molecules at the interface permitting direct monitoring of their reactions, when placed on top of flat metallic surfaces. Experiments with different metal oxide materials and molecules indicate that the enhancement is generic and operates at the single nanoparticle level. Simulations confirm that the amplification is principally electromagnetic and is a result of optical modulation of the underlying plasmonic metallic surface by MONPs, which act as scattering antennae and couple light into the confined region sandwiched by the underlying surface. Because of additional functionalities of metal oxides as magnetic, photoelectrochemical and catalytic materials, enhanced Raman scattering mediated by MONPs opens up significant opportunities in fundamental science, allowing direct tracking and understanding of application-specific transformations at such interfaces. We show a first example by monitoring the MONP-assisted photocatalytic decomposition reaction of an organic dye by individual nanoparticles.

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. (C) 2014 AIP Publishing LLC.

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.

CO oxidation and NO reduction on metal surfaces: density functional theory investigations

This article reviews the accumulated theoretical results, in particular density functional theory calculations, on two catalytic processes, CO oxidation and NO reduction on metal surfaces. Owing to their importance in automotive emission control, these two reactions have generated a lot of interest in the last 20 years. Here the pathways and energetics of the involved elementary reactions under different catalytic conditions are described in detail and the understanding of the reactions is generalized. It is concluded that density functional theory calculations can be applied to catalysis to elucidate mechanisms of complex surface reactions and to understand the electronic structure of chemical processes in general. The achieved molecular knowledge of chemical reactions is certainly beneficial to new catalyst design.

The Validation of a Two-Dimensional Transient Catalyst Model for Direct Injection Two-Stroke Applications

This paper describes the detailed validation of a computer model designed to simulate the transient light-off in a two-stroke oxidation catalyst. A plug flow reactor is employed to provide measurements of temperature and gas concentration at various radial and axial locations inside the catalyst. These measurements are recorded at discrete intervals during a transient light-off in which the inlet temperature is increased from ambient to 300oC at rates of up to 6oC/sec. The catalyst formulation used in the flow reactor, and its associated test procedures, are then simulated by the computer and a comparison made between experimental readings and model predictions. The design of the computer model to which this validation exercise relates is described in detail in a separate technical paper. The first section of the paper investigates the warm-up characteristics of the substrate and examines the validity of the heat transfer predictions between the wall and the gas in the absence of chemical reactions. The predictions from a typical single-component CO transient light-off test are discussed in the second section and are compared with experimental data. In particular the effect of the temperature ramp on the light-off curve and reaction zone development is examined. An analysis of the C3H6 conversion is given in the third section while the final section examines the accuracy of the light-off curves which are produced when both CO and C3H6 are present in the feed gas. The analysis shows that the heat and mass transfer calculations provided reliable predictions of the warm-up behaviour and post light-off gas concentration profiles. The self-inhibition and cross-inhibition terms in the global rate expressions were also found to be reasonably reliable although the surface reaction rates required calibration with experimental data.

Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

A spectroscopic study of the chemistry and reactivity of SO2 on Pt{111}: Reactions with O-2, CO and C3H6

The chemisorption and reactivity of SO2 on Pt{111} have been studied by HREELS, XPS, NEXAFS and temperature-programmed desorption. At 160 K SO2 adsorbs intact at high coverages, with eta(2) S-O coordination to the surface. On annealing to 270 K, NEXAFS indicates the SO2 molecular plane essentially perpendicular to the surface. Preadsorbed O-a reacts with SO2 to yield adsorbed SO4, identified as the key surface species responsible for SO2-promoted catalytic alkane oxidation. Coadsorbed CO or propene efficiently reduce SO2 overlayers to deposit S-a, and the implications of this for catalytic systems are discussed.

Examining the redox and formate mechanisms for water-gas shift reaction on Au/CeO2 using density functional theory

We perform DFT calculations to investigate the redox and formate mechanisms of water-gas-shift (WGS) reaction on Au/CeO2 catalysts. In the redox mechanism, we analyze all the key elementary steps and find that the OH cleavage is the key step. Three possible pathways of OH cleavage are calculated: (1) OHad '' + *'--> H-ad' + O-ad"; (2) H-ad' + OHad '' --> H-2(g) + O-ad '' + *'; and (3) OHad" + OHad '' --> 2O(ad '') + H-2(g) (*': the free adsorption sites on the oxides; ad': adsorption on the metal; ad": adsorption on the oxide, respectively). In the formate mechanism, we identify all the possible pathways for the formation and decomposition of surface formates in the WGS reaction. It is found that there is a shortcoming in the redox and formate mechanisms which is related to surface oxygen reproduction. Four possible pathways for producing surface oxygen are studied, and all the barriers of the four pathways are more than 1 eV. Our results suggest that the processes to reproduce surface oxygen in the reaction circle are not kinetically easy. (C) 2008 Elsevier B.V. All rights reserved.

Dioxygenase-catalysed mono-, di- and tri-oxygenation of dialkyl sulfides and thioacetals: chemoenzymatic synthesis of enantiopure cis-diol sulfoxides

Toluene dioxygenase (TDO)-catalysed monooxygenation of methylsulfanylmethyl phenyl sulfide 1 and methylsulfanylmethyl 2-pyridyl sulfide 4, using whole cells of Pseudomonas putida UV4, occurred exclusively at the alkyl aryl sulfur centre to yield the alkyl aryl sulfoxides 2 and 5 respectively. These sulfoxides, accompanied by the dialkyl sulfoxides 3 and 6, were also obtained from naphthalene dioxygenase (NDO)-catalysed sulfoxidation of thioacetals 1 and 4 using intact cells of P. putida NCIMB 8859. Enzymatic oxidation of methyl benzyl sulfide 7, 2-phenyl-1,3-dithiane 19, and 2-phenyl-1,3-dithiolane 23, using TDO, gave the corresponding dialkyl sulfoxides 8, 20 and 24 as minor bioproducts. TDO-catalysed dioxygenation of the alkyl benzyl sulfides 7, 15 and 17 and the thioacetals 19 and 23, with P. putida UV4, yielded the corresponding enantiopure cis-dihydrodiols 9, 16, 18, 21 and 25 as major metabolites and cis-dihydrodiol sulfoxides 14, 22 and 26 as minor metabolites, resulting from a tandem trioxygenation of substrates 7, 19 and 23 respectively. Chemical oxidation, of the enantiopure cis-dihydrodiol sulfides 9, 16, 18 and 21 with dimethyldioxirane (DMD), gave separable mixtures of the corresponding pairs of cis-dihydrodiol sulfoxide diastereoisomers 14 and 27, 28 and 29, 30 and 31, 22 and 32. While dialkyl sulfoxide bioproducts 3, 6, 20 and 24 were of variable enantiopurity (27-greater than or equal to 98% ee), alkyl aryl monosulfoxides 2 and 5, cis-dihydrodiols 9, 16, 18, 21 and 25 and cis-dihydrodiol sulfoxide bioproducts 14, 22 and 26 were all single enantiomers (greater than or equal to 98% ee). The absolute configurations of the products, obtained from enzyme-catalysed (TDO and NDO) and chemical (DMD) oxidation methods, were determined by stereochemical correlation, circular dichroism, and X-ray crystallographic methods.

Effect of heat treatment on the antioxidant activity, color, and free phenolic acid profile of malt

Green malt was kilned at 95 degrees C following two regimens: a standard regimen (SKR) and a rapid regimen (RKR). Both resulting malts were treated further in a tray dryer heated to 120 degrees C, as was green malt previously dried to 65 degrees C (TDR). Each regimen was monitored by determining the color, antioxidant activity (by both ABTS(center dot+) and FRAP methods), and polyphenolic profile. SKR and RKR malts exhibited decreased L* and increased b* values above approximately 80 degrees C. TDR malts changed significantly less, and color did not develop until 110 degrees C, implying that different chemical reactions lead to color in those malts. Antioxidant activity increased progressively with each regimen, although with TDR malts this became significant only at 110-120 degrees C. The RKR malt ABTS(center dot+) values were higher than those of the SKR malt. The main phenolics, that is, ferulic, p-coumaric, and vanillic acids, were monitored throughout heating. Ferulic acid levels increased upon heating to 80 degrees C for SKR and to 70 degrees C for RKR, with subsequent decreases. However, the levels for TDR malts did not increase significantly. The increase in free phenolics early in kilning could be due to enzymatic release of bound phenolics and/or easier extractability due to changes in the matrix. The differences between the kilning regimens used suggest that further modification of the regimens could lead to greater release of bound phenolics with consequent beneficial effects on flavor stability in beer and, more generally, on human health.

Plasma chloriding of thin-film silver - A novel process in silver-silver chloride reference electrode fabrication

Silver thin films were modified using a novel plasma modification process for the development of thin-film silver-silver chloride reference electrodes. The surface, physical, and electrochemical properties of these electrodes were investigated by atomic force microscopy, thickness and resistivity measurement techniques, as well as impedance spectroscopy and potentiometry. After plasma treatment, thin-film growth was observed and the electrodes, in general, exhibited low interface impedance and a roughened surface. Evidence of a complex surface reorganization was found. Correlating plasma conditions with film properties suggested that increasing pressure and exposure duration increased species availability, therefore governing the reaction rates, while input power appeared to influence the type of surface chemical reactions. Results also indicated that Ar/Cl-2 mixtures should be employed rather than pure chlorine plasmas. (C) 2002 The Electrochemical Society.

Cyanopolyynes in hot cores: modelling G305.2+0.2

We present results from a time-dependent gas-phase chemical model of a hot core based on the physical conditions of G305.2+0.2. While the cyanopolyyne HC3N has been observed in hot cores, the longer chained species, HC5N, HC7N and HC9N, have not been considered as the typical hot-core species. We present results which show that these species can be formed under hot core conditions. We discuss the important chemical reactions in this process and, in particular, show that their abundances are linked to the parent species acetylene which is evaporated from icy grain mantles. The cyanopolyynes show promise as ‘chemical clocks’ which may aid future observations in determining the age of hot core sources. The abundance of the larger cyanopolyynes increases and decreases over relatively short time-scales, ~10^2.5 yr. We present results from a non-local thermodynamic equilibrium statistical equilibrium excitation model as a series of density, temperature and column density dependent contour plots which show both the line intensities and several line ratios. These aid in the interpretation of spectral-line data, even when there is limited line information available. In particular, non-detections of HC5N and HC7N in Walsh et al. are analysed and discussed.