A potential energy surface for the ground state of acetylene, C2H2

A potential energy function has been derived for the ground state surface of C2H2 as a many-body expansion. The 2- and 3-body terms have been obtained by preliminary investigation of the ground state surfaces of CH2( 3B1) and C2H( 2Σ+). A 4-body term has been derived which reproduces the energy, geometry and harmonic force field of C2H2. The potential has a secondary minimum corresponding to the vinylidene structure and the geometry and energy of this are in close agreement with predictions from ab initio calculations. The saddle point for the HCCH-H2CC rearrangement is predicted to lie 2•530 eV above the acetylene minimum.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

Analytical potentials for triatomic molecules VIII. A two valued surface for the lowest 1A¢ surfaces of H2O

Ab initio calculations of the energy have been made at approximately 150 points on the two lowest singlet A' potential energy surfaces of the water molecule, 1A' and 1A', covering structures having D∞h, C∞v, C2v and Cs symmetries. The object was to obtain an ab initio surface of uniform accuracy over the whole three-dimensional coordinate space. Molecular orbitals were constructed from a double zeta plus Rydberg basis, and correlation was introduced by single and double excitations from multiconfiguration states which gave the correct dissociation behaviour. A two-valued analytical potential function has been constructed to fit these ab initio energy calculations. The adiabatic energies are given in our analytical function as the eigenvalues of a 2 2 matrix, whose diagonal elements define two diabatic surfaces. The off-diagonal element goes to zero for those configurations corresponding to surface intersections, so that our adiabatic surface exhibits the correct Σ/II conical intersections for linear configurations, and singlet/triplet intersections of the O + H2 dissociation fragments. The agreement between our analytical surface and experiment has been improved by using empirical diatomic potential curves in place of those derived from ab initio calculations.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Potential energy surface and MULTIMODE vibrational analysis of C2H3+

A full dimensional, ab initio-based semiglobal potential energy surface for C2H3+ is reported. The ab initio electronic energies for this molecule are calculated using the spin-restricted, coupled cluster method restricted to single and double excitations with triples corrections [RCCSD(T)]. The RCCSD(T) method is used with the correlation-consistent polarized valence triple-zeta basis augmented with diffuse functions (aug-cc-pVTZ). The ab initio potential energy surface is represented by a many-body (cluster) expansion, each term of which uses functions that are fully invariant under permutations of like nuclei. The fitted potential energy surface is validated by comparing normal mode frequencies at the global minimum and secondary minimum with previous and new direct ab initio frequencies. The potential surface is used in vibrational analysis using the "single-reference" and "reaction-path" versions of the code MULTIMODE. (c) 2006 American Institute of Physics.

A reaction surface Hamiltonian study of malonaldehyde

We report calculations using a reaction surface Hamiltonian for which the vibrations of a molecule are represented by 3N-8 normal coordinates, Q, and two large amplitude motions, s(1) and s(2). The exact form of the kinetic energy operator is derived in these coordinates. The potential surface is first represented as a quadratic in Q, the coefficients of which depend upon the values of s(1),s(2) and then extended to include up to Q(6) diagonal anharmonic terms. The vibrational energy levels are evaluated by solving the variational secular equations, using a basis of products of Hermite polynomials and appropriate functions of s(1),s(2). Our selected example is malonaldehyde (N=9) and we choose as surface parameters two OH distances of the migrating H in the internal hydrogen transfer. The reaction surface Hamiltonian is ideally suited to the study of the kind of tunneling dynamics present in malonaldehyde. Our results are in good agreement with previous calculations of the zero point tunneling splitting and in general agreement with observed data. Interpretation of our two-dimensional reaction surface states suggests that the OH stretching fundamental is incorrectly assigned in the infrared spectrum. This mode appears at a much lower frequency in our calculations due to substantial transition state character. (c) 2006 American Institute of Physics.

The vibrations and tunnelling of malonaldehyde on a Moller-Plesset surface

The vibrations and tunnelling motion of malonaldehyde have been studied in their full dimensionality using an internal coordinate path Hamiltonian. In this representation there is one large amplitude internal coordinate s and 3N - 7 (=20) normal coordinates Q which are orthogonal to the large amplitude motion at all points. It is crucial that a high accuracy potential energy surface is used in order to obtain a good representation for the tunneling motion; we use a Moller-Plesset (MP2) surface. Our methodology is variational, that is we diagonalize a sufficiently large matrix in order to obtain the required vibrational levels, so an exact representation for the kinetic energy operator is used. In a harmonic valley representation (s, Q) complete convergence of the normal coordinate motions and the internal coordinate motions has been obtained; for the anharmonic valley in which we use two- and three-body terms in the surface (s, Q(1), Q(2)), we also obtain complete convergence. Our final computed stretching fundamentals are deficient because our potential energy surface is truncated at quartic terms in the normal coordinates, but our lower fundamentals are good.

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

Variability in surface atmospheric electric field measurements

An electrical current of the order one picoamp per metre squared flows vertically in the Earth's atmosphere, between the ionosphere at approximately 50km altitude and the surface. This current is generated by global thunderstorm activity and is modulated by galactic cosmic rays and atmospheric aerosol. In fair weather conditions, this current cause a vertical atmospheric electric field, commonly measured as a potential gradient. For circumstances other than fair weather conditions, the potential gradient varies, from small steady enhancements in fog to large fluctuations in thunderstorms. The atmospheric potential gradient is continuously monitored at the Reading University Atmospheric Observatory. An account of the variability of the potential gradient on a variety of time scales will be presented.

Adsorption and dissociation of benzene on bimetallic surfaces - the influence of surface geometry and electronic structure

This topical review discusses the influence of the surface geometry (e.g. lattice parameters and termination) and electronic structure of well-defined bimetallic surfaces on the adsorption and dissociation of benzene. The available data can be divided into two categories with combinations of non-transition metals and transition metals on the one side and combinations of two transition metals on the other. The main effect of non-transition metals in surface alloys is site blocking which can suppress chemisorption and dissociation of the molecules completely. When two transition metals are combined, the effects are less dramatic. They mainly affect the strength of the chemisorption bond and the degree of dissociation due to electronic and template effects.

Charge transfer in Cr adsorption and reaction at the rutile TiO2(110) surface

The rutile TiO2(110) surface has been doped with sub-monolayer metallic Cr, which oxidises and donates charge to specific surface Ti ions. X-Ray and ultra violet photoemission spectroscopy and first principles density functional theory with Hubbard U are used to assign the oxidation states of Cr and surface Ti and we find that Cr2+ forms on bridging oxygen ions and a 5-fold coordinated surface Ti atom is reduced to Ti3+ and the Cr ions readily react with oxygen (to Cr3+), which leads to depletion of surface Ti3+ 3d electrons.

Solving surface structures from normal incidence X-ray standing wave data

A program is provided to determine structural parameters of atoms in or adsorbed on surfaces by refinement of atomistic models towards experimentally determined data generated by the normal incidence X-ray standing wave (NIXSW) technique. The method employs a combination of Differential Evolution Genetic Algorithms and Steepest Descent Line Minimisations to provide a fast, reliable and user friendly tool for experimentalists to interpret complex multidimensional NIXSW data sets.

Coupling aerosol surface and bulk chemistry with a kinetic double layer model (K2-SUB): oxidation of oleic acid by ozone

We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an approximate upper limit of similar to 10(-11) cm(2) s(-1) for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance.

Kinetic multi-layer model of aerosol surface and bulk chemistry (KM-SUB): the influence of interfacial transport and bulk diffusion on the oxidation of oleic acid by ozone

We present a novel kinetic multi-layer model that explicitly resolves mass transport and chemical reaction at the surface and in the bulk of aerosol particles (KM-SUB). The model is based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007), and it includes reversible adsorption, surface reactions and surface-bulk exchange as well as bulk diffusion and reaction. Unlike earlier models, KM-SUB does not require simplifying assumptions about steady-state conditions and radial mixing. The temporal evolution and concentration profiles of volatile and non-volatile species at the gas-particle interface and in the particle bulk can be modeled along with surface concentrations and gas uptake coefficients. In this study we explore and exemplify the effects of bulk diffusion on the rate of reactive gas uptake for a simple reference system, the ozonolysis of oleic acid particles, in comparison to experimental data and earlier model studies. We demonstrate how KM-SUB can be used to interpret and analyze experimental data from laboratory studies, and how the results can be extrapolated to atmospheric conditions. In particular, we show how interfacial and bulk transport, i.e., surface accommodation, bulk accommodation and bulk diffusion, influence the kinetics of the chemical reaction. Sensitivity studies suggest that in fine air particulate matter oleic acid and compounds with similar reactivity against ozone (carbon-carbon double bonds) can reach chemical lifetimes of many hours only if they are embedded in a (semi-)solid matrix with very low diffusion coefficients (< 10(-10) cm(2) s(-1)). Depending on the complexity of the investigated system, unlimited numbers of volatile and non-volatile species and chemical reactions can be flexibly added and treated with KM-SUB. We propose and intend to pursue the application of KM-SUB as a basis for the development of a detailed master mechanism of aerosol chemistry as well as for the derivation of simplified but realistic parameterizations for large-scale atmospheric and climate models.