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.

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.

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.

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.

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.

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.

Carbon chemistry in Galactic bulge planetary nebulae

Galactic bulge planetary nebulae show evidence of mixed chemistry with emission from both silicate dust and polycyclic aromatic hydrocarbons (PAHs). This mixed chemistry is unlikely to be related to carbon dredge-up, as third dredge-up is not expected to occur in the low-mass bulge stars. We show that the phenomenon is widespread and is seen in 30 nebulae out of 40 of our sample, selected on the basis of their infrared flux. Hubble Space Telescope (HST) images and Ultraviolet and Visual Echelle Spectrograph (UVES) spectra show that the mixed chemistry is not related to the presence of emission-line stars, as it is in the Galactic disc population. We also rule out interaction with the interstellar medium (ISM) as origin of the PAHs. Instead, a strong correlation is found with morphology and the presence of a dense torus. A chemical model is presented which shows that hydrocarbon chains can form within oxygen-rich gas through gas-phase chemical reactions. The model predicts two layers, one at A_V~ 1.5, where small hydrocarbons form from reactions with C+, and one at A_V~ 4, where larger chains (and by implication, PAHs) form from reactions with neutral, atomic carbon. These reactions take place in a mini-photon-dominated region (PDR). We conclude that the mixed-chemistry phenomenon occurring in the Galactic bulge planetary nebulae is best explained through hydrocarbon chemistry in an ultraviolet (UV)-irradiated, dense torus.

Gas and grain chemistry in a protoplanetary disk

The chemistry in a protoplanetary accretion disk is modelled between a radius of 100 and 0.1 AU of the central object. We find that interaction of the gas with the dust grains is very important, both by removing a large fraction of the material from the gas in the outer regions and through the chemical reactions which can occur on the dust grain surfaces. In addition, collision with grains neutralises gaseous ions effectively and keeps the ionization fraction low. This results in a chemistry which is dominated by neutral-neutral reactions, even if ionization is provided by cosmic rays or by the decay of radioactive isotopes. We model the effects of two desorption processes with very different efficiencies and find that while these produce similar results over much of the disk for many species, some molecules are extremely sensitive to the nature of the desorption and may one day be used as an observational test for the desorption process.

Galactic Bulge PNe:Carbon molecules in oxygen-rich environments

Galactic bulge planetary nebulae show evidence of mixed chemistry with emission from both silicate dust and PAHs. This mixed chemistry is unlikely to be related to carbon dredge up, as third dredge-up is not expected to occur in the low mass bulge stars. We show that the phenomenon is widespread, and is seen in 30 nebulae out of our sample of 40. A strong correlation is found between strength of the PAH bands and morphology, in particular, the presence of a dense torus. A chemical model is presented which shows that hydrocarbon chains can form within oxygen-rich gas through gas-phase chemical reactions. We conclude that the mixed chemistry phenomenon occurring in the galactic bulge planetary nebulae is best explained through hydrocarbon chemistry in an UV-irradiated, dense torus. © 2012 International Astronomical Union.