A study of the thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) with atomic chlorine (Cl) in the absence and presence of water, using electronic structure methods.

The thermodynamics and mechanisms of the atmospherically relevant reaction dimethyl sulphide (DMS) + atomic chlorine (Cl) were investigated in the absence and presence of a single water molecule, using electronic structure methods. Stationary points on each reaction surface were located using density functional theory (DFT) with the M06-2X functional with aug-cc-pVDZ (aVDZ) and aug-cc-pVTZ (aVTZ) basis sets. Then fixed point calculations were carried out using the UM06-2X/aVTZ optimised stationary point geometries, with aug-cc-pVnZ basis sets (n = T and Q), using the coupled cluster method [CCSD(T)], as well as the domain-based local pair natural orbitals coupled cluster [DLPNO-UCCSD(T)] approach. Four reaction channels are possible, formation of (A) CH3SCH2 + HCl, (B) CH3S + CH3Cl, (C) CH3SCl + CH3, and (C') CH3S(Cl)CH3. The results show that, in the absence of water, channels A and C' are the dominant channels. In the presence of water, the calculations show that the reaction mechanisms for A and C formation change significantly. Channel A occurs via submerged TSs and is expected to be rapid. Channel B occurs via TSs which present significant energy barriers indicating that this channel is not significant in the presence of water relative to CH3SCH2 + HCl and DMS·Cl adduct formation, as is the case in the absence of water. Channel C was not considered as it is endothermic in the absence of water. In the presence of water, pathways which proceed via (a) DMS·H2O + Cl, (b) Cl·H2O + DMS and (c) DMS·Cl + H2O were considered. It was found that under tropospheric conditions, reactions via pathway (b) are of minor importance relative to those that proceed via pathways (a) and (c). This study has shown that water changes the mechanisms of the DMS + Cl reactions significantly but the presence of water is not expected to affect the overall reaction rate coefficient under atmospheric conditions as the DMS + Cl reaction has a rate coefficient at room temperature close to the collisional limit.

Comment on "Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products" by N. T. Tsona, S. Tang and L. Du, Phys. Chem. Chem. Phys., 2019, 21, 20296.

The reaction, BrO + HO2 → HOBr + O2, is exothermic and can produce O2 in both its ground state (X[combining tilde]3∑g-) and its first excited state (ã1Δg). As a result, this reaction can proceed on both a singlet and a triplet potential energy surface. Recently, Tsona, Tang and Du published a paper entitled "Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products (Phys. Chem. Chem. Phys. 2019, 21, 20296-203072). The results of this work showed significant differences from those published earlier on this reaction by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569). Further calculations performed in this present work, combined with higher level calculations published by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569), demonstrate that the work of Tsona et al. is flawed because the integration grid size used in their lowest singlet and triplet calculations is too small, and a closed-shell wavefunction, rather than an open-shell wavefunction, has been used for the singlet surface. The major conclusion in the work of Tsona et al. that the lowest singlet and triplet channels are barrierless is shown to be incorrect. Also, the computed rate coefficients of Tsona et al. showed a positive temperature dependence, which is inconsistent with the experimentally observed negative temperature dependence, whereas the singlet rate coefficients computed by Chow et al. (Phys. Chem. Chem. Phys. 2016, 18, 30554-30569) showed a negative temperature dependence consistent with experiment.

The Atmospherically Important Reaction of Hydroxyl Radicals with Methyl Nitrate: A Theoretical Study Involving the Calculation of Reaction Mechanisms, Enthalpies, Activation Energies, and Rate Coefficients.

A theoretical study, involving the calculation of reaction enthalpies, activation energies, mechanisms, and rate coefficients, was made of the reaction of hydroxyl radicals with methyl nitrate, an important process for methyl nitrate removal in the earth's atmosphere. Four reaction channels were considered: formation of H2O + CH2ONO2, CH3OOH + NO2, CH3OH + NO3, and CH3O + HNO3. For all channels, geometry optimization and frequency calculations were performed at the M06-2X/6-31+G** level, while relative energies were improved at the UCCSD(T*)-F12/CBS level. The major channel is found to be the H abstraction channel, to give the products H2O + CH2ONO2. The reaction enthalpy (ΔH298 KRX) of this channel is computed as -17.90 kcal mol-1. Although the other reaction channels are also exothermic, their reaction barriers are high (>24 kcal mol-1), and therefore these reactions do not contribute to the overall rate coefficient in the temperature range considered (200-400 K). Pathways via three transition states were identified for the H abstraction channel. Rate coefficients were calculated for these pathways at various levels of variational transition state theory including tunneling. The results obtained are used to distinguish between two sets of experimental rate coefficients, measured in the temperature range of 200-400 K, one of which is approximately an order of magnitude greater than the other. This comparison, as well as the temperature dependence of the computed rate coefficients, shows that the lower experimental values are favored. The implications of the results to atmospheric chemistry are discussed.

Direct Measurements of Unimolecular and Bimolecular Reaction Kinetics of the Criegee Intermediate (CH3)2COO.

The Criegee intermediate acetone oxide, (CH3)2COO, is formed by laser photolysis of 2,2-diiodopropane in the presence of O2 and characterized by synchrotron photoionization mass spectrometry and by cavity ring-down ultraviolet absorption spectroscopy. The rate coefficient of the reaction of the Criegee intermediate with SO2 was measured using photoionization mass spectrometry and pseudo-first-order methods to be (7.3 ± 0.5) × 10-11 cm3 s-1 at 298 K and 4 Torr and (1.5 ± 0.5) × 10-10 cm3 s-1 at 298 K and 10 Torr (He buffer). These values are similar to directly measured rate coefficients of anti-CH3CHOO with SO2, and in good agreement with recent UV absorption measurements. The measurement of this reaction at 293 K and slightly higher pressures (between 10 and 100 Torr) in N2 from cavity ring-down decay of the ultraviolet absorption of (CH3)2COO yielded even larger rate coefficients, in the range (1.84 ± 0.12) × 10-10 to (2.29 ± 0.08) × 10-10 cm3 s-1. Photoionization mass spectrometry measurements with deuterated acetone oxide at 4 Torr show an inverse deuterium kinetic isotope effect, kH/kD = (0.53 ± 0.06), for reactions with SO2, which may be consistent with recent suggestions that the formation of an association complex affects the rate coefficient. The reaction of (CD3)2COO with NO2 has a rate coefficient at 298 K and 4 Torr of (2.1 ± 0.5) × 10-12 cm3 s-1 (measured with photoionization mass spectrometry), again similar to rate for the reaction of anti-CH3CHOO with NO2. Cavity ring-down measurements of the acetone oxide removal without added reagents display a combination of first- and second-order decay kinetics, which can be deconvolved to derive values for both the self-reaction of (CH3)2COO and its unimolecular thermal decay. The inferred unimolecular decay rate coefficient at 293 K, (305 ± 70) s-1, is similar to determinations from ozonolysis. The present measurements confirm the large rate coefficient for reaction of (CH3)2COO with SO2 and the small rate coefficient for its reaction with water. Product measurements of the reactions of (CH3)2COO with NO2 and with SO2 suggest that these reactions may facilitate isomerization to 2-hydroperoxypropene, possibly by subsequent reactions of association products.

Simulation of the photodetachment spectrum of HHfO- using coupled-cluster calculations.

The photodetachment spectrum of HHfO- was simulated using restricted-spin coupled-cluster single-double plus perturbative triple {RCCSD(T)} calculations performed on the ground electronic states of HHfO and HHfO-, employing basis sets of up to quintuple-zeta quality. The computed RCCSD(T) electron affinity of 1.67 ± 0.02 eV at the complete basis set limit, including Hf 5s25p6 core correlation and zero-point energy corrections, agrees well with the experimental value of 1.70 ± 0.05 eV from a recent photodetachment study [X. Li et al., J. Chem. Phys. 136, 154306 (2012)]. For the simulation, Franck-Condon factors were computed which included allowances for anharmonicity and Duschinsky rotation. Comparisons between simulated and experimental spectra confirm the assignments of the molecular carrier and electronic states involved but suggest that the experimental vibrational structure has suffered from poor signal-to-noise ratio. An alternative assignment of the vibrational structure to that suggested in the experimental work is presented.

A theoretical study of the atmospherically important radical-radical reaction BrO + HO2; the product channel O2(a1Δg) + HOBr is formed with the highest rate.

A theoretical study has been made of the BrO + HO2 reaction, a radical-radical reaction which contributes to ozone depletion in the atmosphere via production of HOBr. Reaction enthalpies, activation energies and mechanisms have been determined for five reaction channels. Also rate coefficients have been calculated, in the atmospherically important temperature range 200-400 K, for the two channels with the lowest activation energies, both of which produce HOBr: (R1a) HOBr(X1A') + O2(X3Σ) and (R1b) HOBr(X1A') + O2(a1Δg). The other channels considered are: (R2) BrO + HO2 → HBr + O3, (R3) BrO + HO2 → OBrO + OH and (R4) BrO + HO2 → BrOO + OH. For all channels, geometry optimization and frequency calculations were carried out at the M06-2X/AVDZ level, while relative energies of the stationary points on the reaction surface were improved at a higher level (BD(TQ)/CBS or CCSD(T)/CBS). The computed standard reaction enthalpies (ΔH) for channels (R1a), (R1b), (R2), (R3) and (R4) are -47.5, -25.0, -4.3, 14.9 and 5.9 kcal mol-1, and the corresponding computed activation energies (ΔE) are 2.53, -3.07, 11.83, 35.0 and 37.81 kcal mol-1. These values differ significantly from those obtained in earlier work by Kaltsoyannis and Rowley (Phys. Chem. Chem. Phys., 2002, 4, 419-427), particularly for channel (R1b), and reasons for this are discussed. In particular, the importance of obtaining an open-shell singlet wavefunction, rather than a closed-shell singlet wavefunction, for the transition state of this channel is emphasized. Rate coefficient calculations from computed potential energy surfaces were made for BrO + HO2 for the first time. Although channel (R1a) is the most exothermic, channel (R1b) has the lowest barrier height, which is negative (at -3.07 kcal mol-1). Most rate coefficient calculations were therefore made for (R1b). A two transition state model has been used, involving an outer and an inner transition state. The inner transition state was found to be the major bottleneck of the reaction with the outer transition state having essentially no effect on the overall rate coefficient (k) in the temperature range considered. Studying the entropy, enthalpy and free energy of activation of this channel as a function of temperature shows that the main contributor to the magnitude of ln k at a selected temperature is the entropy term (ΔS#/kB) whereas the temperature dependence of ln k is determined mainly by the enthalpy term (-ΔH#/kBT). This compares with reactions with positive barrier heights where the enthalpy term makes a bigger contribution to ln k. Comparison of the computed rate coefficients with available experimental values shows that the computed values have a negative temperature dependence, as observed experimentally, but are too low by approximately an order of magnitude at any selected temperature in the range 200-400 K.

Simulation of the single-vibronic-level emission spectra of HAsO and DAsO.

The single-vibronic-level (SVL) emission spectra of HAsO and DAsO have been simulated by electronic structure/Franck-Condon factor calculations to confirm the spectral molecular carrier and to investigate the electronic states involved. Various multi-reference (MR) methods, namely, NEVPT2 (n-electron valence state second order perturbation theory), RSPT2-F12 (explicitly correlated Rayleigh-Schrodinger second order perturbation theory), and MRCI-F12 (explicitly correlated multi-reference configuration interaction) were employed to compute the geometries and relative electronic energies for the X̃(1)A(') and Ã(1)A(″) states of HAsO. These are the highest level calculations on these states yet reported. The MRCI-F12 method gives computed T0 (adiabatic transition energy including zero-point energy correction) values, which agree well with the available experimental T0 value much better than previously computed values and values computed with other MR methods in this work. In addition, the potential energy surfaces of the X̃(1)A(') and Ã(1)A(″) states of HAsO were computed using the MRCI-F12 method. Franck-Condon factors between the two states, which include anharmonicity and Duschinsky rotation, were then computed and used to simulate the recently reported SVL emission spectra of HAsO and DAsO [R. Grimminger and D. J. Clouthier, J. Chem. Phys. 135, 184308 (2011)]. Our simulated SVL emission spectra confirm the assignments of the molecular carrier, the electronic states involved, and the vibrational structures observed in the SVL emission spectra but suggest a loss of intensity in the reported experimental spectra at the low emission energy region almost certainly due to a loss of responsivity near the cutoff region (∼800 nm) of the detector used. Computed and experimentally derived re (equilibrium) and/or r0 {the (0,0,0) vibrational level} geometries of the two states of HAsO are discussed.

A Study of H2O2 with Threshold Photoelectron Spectroscopy (TPES) and Electronic Structure Calculations: Redetermination of the First Adiabatic Ionization Energy (AIE).

In this work, hydrogen peroxide has been studied with threshold photoelectron (TPE) spectroscopy and photoelectron (PE) spectroscopy. The TPE spectrum has been recorded in the 10.0-21.0 eV ionization energy region, and the PE spectrum has been recorded at 21.22 eV photon energy. Five bands have been observed which have been assigned on the basis of UCCSD(T)-F12/VQZ-F12 and IP-EOM CCSD calculations. Vibrational structure has only been resolved in the TPE spectrum of the first band, associated with the X̃(2)Bg H2O2(+) ← X̃(1)A H2O2 ionization, on its low energy side. This structure is assigned with the help of harmonic Franck-Condon factor calculations that use the UCCSD(T)-F12a/VQZ-F12 computed adiabatic ionization energy (AIE), and UCCSD(T)-F12a/VQZ-F12 computed equilibrium geometric parameters and harmonic vibrational frequencies for the H2O2 X̃(1)A state and the H2O2(+) X̃(2)Bg state. These calculations show that the main vibrational structure on the leading edge of the first TPE band is in the O-O stretching mode (ω3) and the HOOH deformation mode (ω4), and comparison of the simulated spectrum to the experimental spectrum gives the first AIE of H2O2 as (10.685 ± 0.005) eV and ω4 = (850 ± 30) and ω3 = (1340 ± 30) cm(-1) in the X̃(2)Bg state of H2O2(+). Contributions from ionization of vibrationally excited levels in the torsion mode have been identified in the TPE spectrum of the first band and the need for a vibrationally resolved TPE spectrum from vibrationally cooled molecules, as well as higher level Franck-Condon factors than performed in this work, is emphasized.

Can Cyclen Bind Alkali Metal Azides? A DFT Study as a Precursor to Synthesis.

Can cyclen (1,4,7,10-tetraazacyclododecane) bind alkali metal azides? This question is addressed by studying the geometric and electronic structures of the alkali metal azide-cyclen [M(cyclen)N3] complexes using density functional theory (DFT). The effects of adding a second cyclen ring to form the sandwich alkali metal azide-cyclen [M(cyclen)2N3] complexes are also investigated. N3(-) is found to bind to a M(+) (cyclen) template to give both end-on and side-on structures. In the end-on structures, the terminal nitrogen atom of the azide group (N1) bonds to the metal as well as to a hydrogen atom of the cyclen ring through a hydrogen bond in an end-on configuration to the cyclen ring. In the side-on structures, the N3 unit is bonded (in a side-on configuration to the cyclen ring) to the metal through the terminal nitrogen atom of the azide group (N1), and through the other terminal nitrogen atom (N3) of the azide group by a hydrogen bond to a hydrogen atom of the cyclen ring. For all the alkali metals, the N3-side-on structure is lowest in energy. Addition of a second cyclen unit to [M(cyclen)N3] to form the sandwich compounds [M(cyclen)2N3] causes the bond strength between the metal and the N3 unit to decrease. It is hoped that this computational study will be a precursor to the synthesis and experimental study of these new macrocyclic compounds; structural parameters and infrared spectra were computed, which will assist future experimental work.

Reaction between CH3O2 and BrO radicals: a new source of upper troposphere lower stratosphere hydroxyl radicals.

Over the last two decades it has emerged that measured hydroxyl radical levels in the upper troposphere are often underestimated by models, leading to the assertion that there are missing sources. Here we report laboratory studies of the kinetics and products of the reaction between CH3O2 and BrO radicals that shows that this could be an important new source of hydroxyl radicals:BrO + CH3O2 → products (1). The temperature dependent value in Arrhenius form of k(T) is k1 = (2.42­0.72+1.02) × 10­14 exp[(1617 ± 94)/T] cm3 molecule­1 s­1. In addition, CH2OO and HOBr are believed to be the major products. Global model results suggest that the decomposition of H2COO to form OH could lead to an enhancement in OH of up to 20% in mid-latitudes in the upper troposphere and in the lower stratosphere enhancements in OH of 2­9% are inferred from model integrations. In addition, reaction 1 aids conversion of BrO to HOBr and slows polar ozone loss in the lower stratosphere.

A theoretical study of the mechanism of the atmospherically relevant reaction of chlorine atoms with methyl nitrate, and calculation of the reaction rate coefficients at temperatures relevant to the troposphere.

The reaction between atomic chlorine (Cl) and methyl nitrate (CH3ONO2) is significant in the atmosphere, as Cl is a key oxidant, especially in the marine boundary layer, and alkyl nitrates are important nitrogen-containing organic compounds, which are temporary reservoirs of the reactive nitrogen oxides NO, NO2 and NO3 (NOx). Four reaction channels HCl + CH2ONO2, CH3OCl + NO2, CH3Cl + NO3 and CH3O + ClNO2 were considered. The major channel is found to be the H abstraction channel, to give the products HCl + CH2ONO2. For all channels, geometry optimization and frequency calculations were carried out at the M06-2X/6-31+G** level, while relative electronic energies were improved to the UCCSD(T*)-F12/CBS level. The reaction barrier (ΔE(‡)0K) and reaction enthalpy (ΔH(RX)298K) of the H abstraction channel were computed to be 0.61 and -2.30 kcal mol(-1), respectively, at the UCCSD(T*)-F12/CBS//M06-2X/6-31+G** level. Reaction barriers (ΔE(‡)0K) for the other channels are more positive and these pathways do not contribute to the overall reaction rate coefficient in the temperature range considered (200-400 K). Rate coefficients were calculated for the H-abstraction channel at various levels of variational transition state theory (VTST) including tunnelling. Recommended ICVT/SCT rate coefficients in the temperature range 200-400 K are presented for the first time for this reaction. The values obtained in the 200-300 K region are particularly important as they will be valuable for atmospheric modelling calculations involving reactions with methyl nitrate. The implications of the results to atmospheric chemistry are discussed. Also, the enthalpies of formation, ΔHf,298K, of CH3ONO2 and CH2ONO2 were computed to be -29.7 and 19.3 kcal mol(-1), respectively, at the UCCSD(T*)-F12/CBS level.

Simulation of the single-vibronic-level emission spectrum of HPS.

We have computed the potential energy surfaces of the X¹A' and ùA" states of HPS using the explicitly correlated multi-reference configuration interaction (MRCI-F12) method, and Franck-Condon factors between the two states, which include anharmonicity and Duschinsky rotation, with the aim of testing the assignment of the recently reported single-vibronic-level (SVL) emission spectrum of HPS [R. Grimminger, D. J. Clouthier, R. Tarroni, Z. Wang, and T. J. Sears, J. Chem. Phys. 139, 174306 (2013)]. These are the highest level calculations on these states yet reported. It is concluded that our spectral simulation supports the assignments of the molecular carrier, the electronic states involved and the vibrational structure of the experimental laser induced fluorescence, and SVL emission spectra proposed by Grimminger et al. [J. Chem. Phys. 139, 174306 (2013)]. However, there remain questions unanswered regarding the relative electronic energies of the two states and the geometry of the excited state of HPS.

Rate coefficients of the Cl + CH3C(O)OCH3 → HCl + CH3C(O)OCH2 reaction at different temperatures calculated by transition-state theory with ab initio and density functional theory reaction paths.

The complex relationship of computed rate coefficients (k's) with different ab initio/DFT and TST levels was studied. The MEPs, gradients, and Hessians of the title reaction were computed using the MP2 and DFT methods. Electronic energies were improved to the UCCSD(T)-F12x/CBS level, and k's were calculated at the TST, CVT, and ICVT levels with various tunnelling corrections. Although computed microcanonical and tunnelling effects are small, computed k(TST) values are larger than computed k(TST/ZCT) and k(TST/SCT) values by 3 orders of magnitude at low temperatures, because computed κ((TST/CAG)) values are as small as 6 × 10(-4). In some cases, the maximum of the ΔG/s curves at a certain T is far away from the MEP maximum. This raises the question of the range of s to be considered in a VTST calculation and, of a possible scenario, where no maximum on the ΔG curve can be located and hence a breakdown of VTST occurs. For dual-level direct dynamics calculations, different entropic contributions from different lower levels can lead to computed k's, which differ by more than 1 order of magnitude. Matching computed and experimental k values leads to an empirical barrier of 1.34 kcal mol(-1) for the title reaction.

Simulated photodetachment spectra of AlH2(-).

We have carried out high-level ab initio calculations on AlH2 and its anion, as well as Franck-Condon factor calculations, which include anharmonicity and Duschinsky rotation, to simulate the photodetachment spectrum of AlH2(-), with the aim of assigning the very recently reported photodetachment spectrum of AlH2(-) [X. Zhang, H. Wang, E. Collins, A. Lim, G. Ganteför, B. Kiran, H. Schnöckel, B. Eichhorn, and K. Bowen, J. Chem. Phys. 138, 124303 (2013)]. However, our simulated spectra do not support the assignment of the reported experimental spectrum to AlH2(-).

Direct measurements of conformer-dependent reactivity of the Criegee intermediate CH3CHOO.

Although carbonyl oxides, "Criegee intermediates," have long been implicated in tropospheric oxidation, there have been few direct measurements of their kinetics, and only for the simplest compound in the class, CH2OO. Here, we report production and reaction kinetics of the next larger Criegee intermediate, CH3CHOO. Moreover, we independently probed the two distinct CH3CHOO conformers, syn- and anti-, both of which react readily with SO2 and with NO2. We demonstrate that anti-CH3CHOO is substantially more reactive toward water and SO2 than is syn-CH3CHOO. Reaction with water may dominate tropospheric removal of Criegee intermediates and determine their atmospheric concentration. An upper limit is obtained for the reaction of syn-CH3CHOO with water, and the rate constant for reaction of anti-CH3CHOO with water is measured as 1.0 × 10(-14) ± 0.4 × 10(-14) centimeter(3) second(-1).

Rate coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 reaction at different temperatures calculated by transition state theory with ab initio and DFT reaction paths.

The minimum energy path (MEP) of the reaction, CF(3)CHFCF(3) + H → transition state (TS) → CF(3)CFCF(3) + H(2), has been computed at different ab initio levels and with density functional theory (DFT) using different functionals. The computed B3LYP/6-31++G**, BH&HLYP/cc-pVDZ, BMK/6-31++G**, M05/6-31+G**, M05-2X/6-31+G**, UMP2/6-31++G**, PUMP2/6-31++G**//UMP2/6-31++G**, RCCSD(T)/aug-cc-pVDZ//UMP2/6-31++G**, RCCSD(T)/aug-cc-pVTZ(spd,sp)//UMP2//6-31++G**, RCCSD(T)/CBS//M05/6-31+G**, and RCCSD(T)/CBS//UMP2/6-31++G** MEPs, and associated gradients and Hessians, were used in reaction rate coefficient calculations based on the transition state theory (TST). Reaction rate coefficients were computed between 300 and 1500 K at various levels of TST, which include conventional TST, canonical variational TST (CVT) and improved CVT (ICVT), and with different tunneling corrections, namely, Wigner, zero-curvature, and small-curvature (SCT). The computed rate coefficients obtained at different ab initio, DFT and TST levels are compared with experimental values available in the 1000-1200 K temperature range. Based on the rate coefficients computed at the ICVT/SCT level, the highest TST level used in this study, the BH&HLYP functional performs best among all the functionals used, while the RCCSD(T)/CBS//MP2/6-31++G** level is the best among all the ab initio levels used. Comparing computed reaction rate coefficients obtained at different levels of theory shows that, the computed barrier height has the strongest effect on the computed reaction rate coefficients as expected. Variational effects on the computed rate coefficients are found to be negligibly small. Although tunneling effects are relatively small at high temperatures (~1500 K), SCT corrections are significant at low temperatures (~300 K), and both barrier heights and the magnitudes of the imaginary frequencies affect SCT corrections.

Regional and global impacts of Criegee intermediates on atmospheric sulphuric acid concentrations and first steps of aerosol formation.

Carbonyl oxides ("Criegee intermediates"), formed in the ozonolysis of alkenes, are key species in tropospheric oxidation of organic molecules and their decomposition provides a non-photolytic source of OH in the atmosphere (Johnson and Marston, Chem. Soc. Rev., 2008, 37, 699, Harrison et al, Sci, Total Environ., 2006, 360, 5, Gäb et al., Nature, 1985, 316, 535, ref. 1-3). Recently it was shown that small Criegee intermediates, C.I.'s, react far more rapidly with SO2 than typically represented in tropospheric models, (Welz, Science, 2012, 335, 204, ref. 4) which suggested that carbonyl oxides could have a substantial influence on the atmospheric oxidation of SO2. Oxidation of 502 is the main atmospheric source of sulphuric acid (H2SO4), which is a critical contributor to aerosol formation, although questions remain about the fundamental nucleation mechanism (Sipilä et al., Science, 2010, 327, 1243, Metzger et al., Proc. Natl. Acad. Sci. U. S. A., 2010 107, 6646, Kirkby et al., Nature, 2011, 476, 429, ref. 5-7). Non-absorbing atmospheric aerosols, by scattering incoming solar radiation and acting as cloud condensation nuclei, have a cooling effect on climate (Intergovernmental Panel on Climate Change (IPCC), Climate Change 2007: The Physical Science Basis, Cambridge University Press, 2007, ref. 8). Here we explore the effect of the Criegees on atmospheric chemistry, and demonstrate that ozonolysis of alkenes via the reaction of Criegee intermediates potentially has a large impact on atmospheric sulphuric acid concentrations and consequently the first steps in aerosol production. Reactions of Criegee intermediates with SO2 will compete with and in places dominate over the reaction of OH with SO2 (the only other known gas-phase source of H2SO4) in many areas of the Earth's surface. In the case that the products of Criegee intermediate reactions predominantly result in H2SO4 formation, modelled particle nucleation rates can be substantially increased by the improved experimentally obtained estimates of the rate coefficients of Criegee intermediate reactions. Using both regional and global scale modelling, we show that this enhancement is likely to be highly variable spatially with local hot-spots in e.g. urban outflows. This conclusion is however contingent on a number of remaining uncertainties in Criegee intermediate chemistry.

Spectroscopy of the simplest Criegee intermediate CH2OO: simulation of the first bands in its electronic and photoelectron spectra.

CH(2)OO, the simplest Criegee intermediate, and ozone are isoelectronic. They both play very important roles in atmospheric chemistry. Whilst extensive experimental studies have been made on ozone, there were no direct gas-phase studies on CH(2)OO until very recently when its photoionization spectrum was recorded and kinetics studies were made of some reactions of CH(2)OO with a number of molecules of atmospheric importance, using photoionization mass spectrometry to monitor CH(2)OO. In order to encourage more direct studies on CH(2)OO and other Criegee intermediates, the electronic and photoelectron spectra of CH(2)OO have been simulated using high level electronic structure calculations and Franck-Condon factor calculations, and the results are presented here. Adiabatic and vertical excitation energies of CH(2)OO were calculated with TDDFT, EOM-CCSD, and CASSCF methods. Also, DFT, QCISD and CASSCF calculations were performed on neutral and low-lying ionic states, with single energy calculations being carried out at higher levels to obtain more reliable ionization energies. The results show that the most intense band in the electronic spectrum of CH(2) OO corresponds to the B(1)A' ← X(1)A' absorption. It is a broad band in the region 250-450 nm showing extensive structure in vibrational modes involving O-O stretching and C-O-O bending. Evidence is presented to show that the electronic absorption spectrum of CH(2)OO has probably been recorded in earlier work, albeit at low resolution. We suggest that CH(2)OO was prepared in this earlier work from the reaction of CH(2)I with O(2) and that the assignment of the observed spectrum solely to CH(2)IOO is incorrect. The low ionization energy region of the photoelectron spectrum of CH(2)OO consists of two overlapping vibrationally structured bands corresponding to one-electron ionizations from the highest two occupied molecular orbitals of the neutral molecule. In each case, the adiabatic component is the most intense and the adiabatic ionization energies of these bands are expected to be very close, at 9.971 and 9.974 eV at the highest level of theory used.

Development and testing of a compact basis set for use in effective core potential calculations on rhodium complexes.

We present a set of effective core potential (ECP) basis sets for rhodium atoms which are of reasonable size for use in electronic structure calculations. In these ECP basis sets, the Los Alamos ECP is used to simulate the effect of the core electrons while an optimized set of Gaussian functions, which includes polarization and diffuse functions, is used to describe the valence electrons. These basis sets were optimized to reproduce the ionization energy and electron affinity of atomic rhodium. They were also tested by computing the electronic ground state geometry and harmonic frequencies of [Rh(CO)(2) µ-Cl](2) , Rh(CO)(2) ClPy, and RhCO (neutral and its positive, and negative ions) as well as the enthalpy of the reaction of [Rh(CO)(2) µ-Cl](2) with pyridine (Py) to give Rh(CO)(2) ClPy, at different levels of theory. Good agreement with experimental values was obtained. Although the number of basis functions used in our ECP basis sets is smaller than those of other ECP basis sets of comparable quality, we show that the newly developed ECP basis sets provide the flexibility and precision required to reproduce a wide range of chemical and physical properties of rhodium compounds. Therefore, we recommend the use of these compact yet accurate ECP basis sets for electronic structure calculations on molecules involving rhodium atoms.