Gas-phase rate coefficients for reactions of NO3, OH, O-3 and O(P-3) with unsaturated alcohols and ethers: correlations and structure-activity relations (SARs)

Experimental difficulties sometimes force modellers to use predicted rate coefficients for reactions of oxygenated volatile organic compounds (oVOCs). We examine here methods for making the predictions for reactions of atmospheric initiators of oxidation, NO3, OH, O-3 and O(P-3), with unsaturated alcohols and ethers. Logarithmic correlations are found between measured rate coefficients and calculated orbital energies, and these correlations may be used directly to estimate rate coefficients for compounds where measurements have not been performed. To provide a shortcut that obviates the need to calculate orbital energies, structure-activity relations (SARs) are developed. Our SARs are tested for predictive power against compounds for which experimental rate coefficients exist, and their accuracy is discussed. Estimated atmospheric lifetimes for oVOCs are presented. The SARs for alkenols successfully predict key rate coefficients, and thus can be used to enhance the scope of atmospheric models incorporating detailed chemistry. SARs for the ethers have more limited applicability, but can still be useful in improving tropospheric models. (C) 2008 Elsevier Ltd. All rights reserved.

Gas-phase rate coefficients for the reactions of NO3, OH and O-3 with alpha,beta-unsaturated esters and ketones: structure-activity relations (SARs)

Gas-phase rate coefficients for the atmospherically important reactions of NO3, OH and O-3 are predicted for 55 alpha,beta-unsaturated esters and ketones. The rate coefficients were calculated using a correlation described previously [Pfrang, C., King, M.D., C. E. Canosa-Mas, C.E., Wayne, R.P., 2006. Atmospheric Environment 40, 1170-1179]. These rate coefficients were used to extend structure-activity relations for predicting the rate coefficients for the reactions of NO3, OH or O-3 with alkenes to include alpha,beta-unsaturated esters and ketones. Conjugation of an alkene with an alpha,beta-keto or alpha,beta-ester group will reduce the value of a rate coefficient by a factor of similar to 110, similar to 2.5 and similar to 12 for reaction with NO3, OH or O-3, respectively. The actual identity of the alkyl group, R, in -C(O)R or -C(O)OR has only a small influence. An assessment of the reliability of the SAR is given that demonstrates that it is useful for reactions involving NO3 and OH, but less valuable for those of O-3 or peroxy nitrate esters. (c) 2006 Elsevier Ltd. All rights reserved.

Gas-phase rate coefficients for the reactions of O(3P), S(3P), Se(3P), and Te(3P) with alkenes: application of perturbation frontier molecular orbital theory, correlations, and structure-activity relations (SARs)

The kinetics of the reactions of the atoms O(P-3), S(P-3), Se(P-3), and Te((3)p) with a series of alkenes are examined for correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes. These correlations may be employed to predict rate coefficients from the calculated HOMO energy of any other alkene of interest. The rate coefficients obtained from the correlations were used to formulate structure-activity relations (SARs) for reactions of O((3)p), S(P-3), Se (P-3), and Te((3)p) with alkenes. A comparison of the values predicted by both the correlations and the SARs with experimental data where they exist allowed us to assess the reliability of our method. We demonstrate the applicability of perturbation frontier molecular orbital theory to gas-phase reactions of these atoms with alkenes. The correlations are apparently not applicable to reactions of C(P-3), Si(P-3), N(S-4), and Al(P-2) atoms with alkenes, a conclusion that could be explained in terms of a different mechanism for reaction of these atoms.

Gas-phase rate coefficients for the reactions of nitrate radicals with (Z)-pent-2-ene, (E)-pent-2-ene, (Z)-hex-2-ene, (E)-hex-2-ene, (Z)-hex-3-ene, (E)-hex-3-ene and (E)-3-methylpent-2-ene at room temperature

Rate coefficients for reactions of nitrate radicals (NO3) with (Z)-pent-2-ene, (E)-pent-2-ene, (Z)-hex-2-ene, (E)-hex-2-ene, (Z)-hex-3-ene, (E)-hex-3-ene and (E)-3-methylpent-2-ene were determined to be (6.55 +/- 0.78) x 10(-13) cm(3) molecule(-1) s(-1), (3.78 +/- 0.45) x 10(-13) cm(3) molecule(-1) s(-1), (5.30 +/- 0.73) x 10(-13) cm(3) molecule(-1) s(-1), (3.83 +/- 0.47) x 10(-13) cm(3) molecule(-1) s(-1), (4.37 +/- 0.49) x 10(-13) cm(3) molecule(-1) s(-1), (3.61 +/- 0.40) x 10(-13) cm(3) molecule(-1) s(-1) and (8.9 +/- 1.5) x 10(-12) cm(3) molecule(-1) s(-1), respectively. We performed kinetic experiments at room temperature and atmospheric pressure using a relative-rate technique with GC-FID analysis. The experimental results demonstrate a surprisingly large cis-trans (Z-E) effect, particularly in the case of the pent-2-enes, where the ratio of rate coefficients is ca. 1.7. Rate coefficients are discussed in terms of electronic and steric influences, and our results give some insight into the effects of chain length and position of the double bond on the reaction of NO3 with unsaturated hydrocarbons. Atmospheric lifetimes were calculated with respect to important oxidants in the troposphere for the alkenes studied, and NO3-initiated oxidation is found to be the dominant degradation route for (Z)-pent-2-ene, (Z)-hex-3-ene and (E)-3-methylpent-2-ene.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

Modeling the mechanisms that control in-stream dissolved organic carbon dynamics in upland and forested catchments

[1] We present a new, process-based model of soil and stream water dissolved organic carbon (DOC): the Integrated Catchments Model for Carbon (INCA-C). INCA-C is the first model of DOC cycling to explicitly include effects of different land cover types, hydrological flow paths, in-soil carbon biogeochemistry, and surface water processes on in-stream DOC concentrations. It can be calibrated using only routinely available monitoring data. INCA-C simulates daily DOC concentrations over a period of years to decades. Sources, sinks, and transformation of solid and dissolved organic carbon in peat and forest soils, wetlands, and streams as well as organic carbon mineralization in stream waters are modeled. INCA-C is designed to be applied to natural and seminatural forested and peat-dominated catchments in boreal and temperate regions. Simulations at two forested catchments showed that seasonal and interannual patterns of DOC concentration could be modeled using climate-related parameters alone. A sensitivity analysis showed that model predictions were dependent on the mass of organic carbon in the soil and that in-soil process rates were dependent on soil moisture status. Sensitive rate coefficients in the model included those for organic carbon sorption and desorption and DOC mineralization in the soil. The model was also sensitive to the amount of litter fall. Our results show the importance of climate variability in controlling surface water DOC concentrations and suggest the need for further research on the mechanisms controlling production and consumption of DOC in soils.

Laboratory kinetic studies of the reactions of Cl atoms with species of biogenic origin: δ3-Carene, isoprene, methacrolein and methyl vinyl ketone

The rate coefficients for the reaction between atomic chlorine and a number of naturally occurring species have been measured at ambient temperature and atmospheric pressure using the relative rate technique. The values obtained were (4.0 ± 0.8) × 10-10, (2.1 ± 0.5) × 10-10, (3.2 ± 0.5) × 10-10, and (4.9 ± 0.5) × 10-10 cm3 molecule-1 s-1, for reactions with isoprene, methyl vinyl ketone, methacrolein and δ3-carene, respectively. The value obtained for isoprene compares favourably with previously reported values. No values have been reported to date for the rate constants of the other reactions.

A comparison of the reactivity of germylene and dimethylgermylene with some methylgermanes. Direct kinetic and quantum chemical studies

Time-resolved studies of germylene, GeH2, and dimethygermylene, GeMe2, generated by the 193 nm laser flash photolysis of appropriate precursor molecules have been carried out to try to obtain rate coefficients for their bimolecular reactions with dimethylgermane, Me2GeH2, in the gas-phase. GeH2 + Me2GeH2 was studied over the pressure range 1-100 Torr with SF6 as bath gas and at five temperatures in the range 296-553 K. Only slight pressure dependences were found (at 386, 447 and 553 K). RRKM modelling was carried out to fit these pressure dependences. The high pressure rate coefficients gave the Arrhenius parameters: log(A/cm(3) molecule(-1)s(-1)) = -10.99 +/- 0.07 and E-a = -(7.35 +/- 0.48) kJ mol(-1). No reaction could be found between GeMe2 + Me2GeH2 at any temperature up to 549 K, and upper limits of ca. 10(-14) cm(3) molecule(-1)s(-1) were set for the rate coefficients. A rate coefficient of (1.33 +/- 0.04) x 10(-11)cm(3) molecule(-1)s(-1) was also obtained for GeH2 + MeGeH3 at 296 K. No reaction was found between GeMe2 and MeGeH3. Rate coefficient comparisons showed, inter alia, that in the substrate germane Me-for-H substitution increased the magnitudes of rate coefficients significantly, while in the germylene Me-for-H substitution decreased the magnitudes of rate coefficients by at least four orders of magnitude. Quantum chemical calculations (G2(MP2,SVP)// B3LYP level) supported these findings and showed that the lack of reactivity of GeMe2 is caused by a positive energy barrier for rearrangement of the initially formed complexes. Full details of the structures of intermediate complexes and the discussion of their stabilities are given in the paper.

The kinetics of the gas-phase decomposition of the 2-methyl-2-butoxyl and 2-methyl-2-pentoxyl radicals

The kinetics of the title reactions have been studied by relative-rate methods as a function of temperature. Relative-rate coefficients for the two decomposition channels of 2-methyl-2-butoxyl have been measured at five different temperatures between 283 and 345 K and the observed temperature dependence is consistent with the results of some previous experimental studies. The kinetics of the two decomposition channels of 2-methyl-2-pentoxyl have also been investigated, as a function of temperature, relative to the estimated rate of isomerisation of this radical. Room-temperature rate coefficient data for the two decomposition channels of both 2-methyl-2-pentoxyl and 2-methyl-2-butxoyl (after combining the relative rate coefficient for this latter with a value for the rate coefficient of the major channel, extrapolated from the data presented by Batt et al., Int. J. Chem. Kinet., 1978, 10, 931) are shown to be consistent with a non-linear kinetic correlation, for alkoxyl radical decomposition rate data, previously presented by this laboratory (Johnson et al., Atmos. Environ., 2004, 38, 1755-1765).

Correlation-type structure activity relationships for the kinetics of gas-phase RO2 self-reactions and reaction with HO2

Structure activity relationships (SARs) are presented for the gas-phase reactions of RO2 with HO2, and the self- and cross-reactions of RO2. For RO2+HO2 the SAR is based upon a correlation between the logarithm of the measured rate coefficient and a calculated ionisation potential for the molecule R-CH=CH2, R being the same group in both the radical and molecular analogue. The correlation observed is strong and only for one RO2 species does the measured rate coefficient deviate by more than a factor of two from the linear least-squares regression line. For the self- and cross-reactions of RO2 radicals, the SAR is based upon a correlation between the logarithm of the measured rate coefficient and the calculated electrostatic potential (ESP) at the equivalent carbon atom in the RH molecule to which oxygen is attached in RO2, again R being the same group in the molecule and the radical. For cases where R is a simple alkyl-group, a strong linear correlation observed. For RO2 radicals which contain lone pair-bearing substituents and for which the calculated ESP<-0.05 self-reaction rate coefficients appear to be insensitive to the value of the ESP. For RO2 of this type with ESP>-0.05 a linear relationship between log k and the ESP is again observed. Using the relationships, 84 out of the 85 rate coefficients used to develop the SARs are predicted to within a factor of three of their measured values. A relationship is also presented that allows the prediction of the Arrhenius parameters for the self-reactions of simple alkyl RO2 radicals. On the basis of the correlations, predictions of room-temperature rate coefficients are made for a number of atmospherically important peroxyl-peroxyl radical reactions. (C) 2003 Elsevier Ltd. All rights reserved.

Reactions of NO3 with the man-made emissions 2-methylpent-2-ene, (Z)-3-methylpent-2-ene, ethyl vinyl ether, and the stress-induced plant emission ethyl vinyl ketone

Rate coefficients for reactions of nitrate radicals (NO3) with the anthropogenic emissions 2-methylpent-2-ene, (Z)-3-methylpent-2-ene.. ethyl vinyl ether, and the stress-induced plant emission ethyl vinyl ketone (pent-1-en-3-one) were determined to be (9.3 +/- 1.1) x 10(-12), (9.3 +/- 3.2) x 10(-12), (1.7 +/- 1.3) x 10(-12) and (9.4 + 2.7) x 10(-17) cm(3) molecule(-1) s(-1). We performed kinetic experiments at room temperature and atmospheric pressure using a relative-rate technique with GC-FID analysis. Experiments with ethyl vinyl ether required a modification of our established procedure that might introduce additional uncertainties, and the errors suggested reflect these difficulties. Rate coefficients are discussed in terms of electronic and steric influences. Atmospheric lifetimes with respect to important oxidants in the troposphere were calculated. NO3-initiated oxidation is found to be the strongly dominating degradation route for 2-methylpent-2-ene, (Z)-3-methylpent-2-ene and ethyl vinyl ether. Atmospheric concentrations of the alkenes and their relative contribution to the total NMHC emissions from trucks can be expected to increase if plans for the introduction of particle filters for diesel engines are implemented on a global scale. Thus more kinetic data are required to better evaluate the impact of these emissions.

Gas-phase reactions of NO3 and N2O5 with (Z)-hex-4-en-1-ol, (Z)-hex-3-en-1-ol ('leaf alcohol'), (E)-hex-3-en-1-ol, (Z)-hex-2-en-1-ol and (E)-hex-2-en-1-ol

The night-time atmospheric chemistry of the biogenic volatile organic compounds (Z)-hex-4-en-1-ol, (Z)-hex-3-en-1-ol ('leaf alcohol'), (E)-hex-3-en-1-ol, (Z)-hex-2-en-1-ol and (E)-hex-2-en-1-ol, has been studied at room temperature. Rate coefficients for reactions of the nitrate radical (NO3) with these stress-induced plant emissions were measured using the discharge-flow technique. We employed off-axis continuous-wave cavity-enhanced absorption spectroscopy (CEAS) for the detection of NO3, which enabled us to work in excess of the hexenol compounds over NO3. The rate coefficients determined were (2.93 +/- 0.58) x 10(-13) cm(3) molecule(-1) s(-1), (2.67 +/- 0.42) x 10(-13) cm(3) molecule(-1) s(-1), (4.43 +/- 0.91) x 10(-13) cm(3) molecule(-1) s(-1), (1.56 +/- 0.24) x 10(-13) cm(3) molecule(-1) s(-1), and (1.30 +/- 0.24) x 10(-13) cm(3) molecule(-1) s(-1) for (Z)-hex-4-en-1-ol, (Z)-hex-3en-1-ol, (E)-hex-3-en-1-ol, (Z)-hex-2-en-1-ol and (E)-hex-2-en-1-ol. The rate coefficient for the reaction of NO3 with (Z)-hex-3-en-1-ol agrees with the single published determination of the rate coefficient using a relative method. The other rate coefficients have not been measured before and are compared to estimated values. Relative-rate studies were also performed, but required modification of the standard technique because N2O5 (used as the source of NO3) itself reacts with the hexenols. We used varying excesses of NO2 to determine simultaneously rate coefficients for reactions of NO3 and N2O5 with (E)-hex-3-en-1-ol of (5.2 +/- 1.8) x 10(-13) cm(3) molecule(-1) s(-1) and (3.1 +/- 2.3) x 10(-18) cm(3) molecule(-1) s(-1). Our new determinations suggest atmospheric lifetimes with respect to NO3-initiated oxidation of roughly 1-4 h for the hexenols, comparable with lifetimes estimated for the atmospheric degradation by OH and shorter lifetimes than for attack by O-3. Recent measurements of [N2O5] suggest that the gas-phase reactions of N2O5 with unsaturated alcohols will not be of importance under usual atmospheric conditions, but they certainly can be in laboratory systems when determining rate coefficients.

Structure–activity relations (SARs) for gas-phase reactions of NO3, OH and O3 with alkenes: an update

We present updated structure-activity relations (SARs) for the prediction of rate coefficients for gas-phase reactions with alkenes of the major atmospheric oxidants NO3, OH and O-3. Such SARs provide one way of incorporating essential information about reactivity into atmospheric models. Rate coefficients obtained from correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes were used to refine the SARs. SARs have an advantage for the user over the direct application of the correlations in that knowledge of the structure of the alkene of interest is sufficient to estimate rate coefficients, and no quantum-mechanical calculations need to be performed. A comparison of the values predicted by the SARs with experimental data where they exist allowed us to assess the reliability of our method.

Correlations for gas-phase reactions of NO3, OH and O3 with alkenes: an update

Methods are developed for predicting rate coefficients for reactions of initiators of tropospheric oxidation with unsaturated compounds that are abundant in the atmosphere; prognostic tools of this kind are essential for atmospheric chemists and modellers. To pursue the aim of exploring such tools, the kinetics of reactions of NO3, OH and O-3 with a series of alkenes are examined for correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes. A comparison of the values predicted by the correlations with experimental data (where the latter exist) allowed us to assess the reliability of our method. We used a series of theoretical methods to calculate the HOMO energies, and found that higher computational effort improves the agreement of the predicted rate coefficients with experimental values, especially for reactions of NO3 with alkenes that possess vinyllic halogen substituents. As a consequence, it is expedient to suggest new correlations to replace those presented by us and others that were based on the lower level of theory. We propose the following correlations for the reactions of NO3, OH and O-3 with alkenes: ln(k(NO3)/cm(3) molecule(-1) s(-1)) = 6.40(E-HOMO/eV) + 31.69, ln(k(OH)/cm(3) molecule(-1) s(-1)) = 1.21 (E-HOMO/eV)-12.34 and ln(k(O3)/cm(3) molecule(-1) s(-1)) = 3.28(E-HOMO/eV)-6.78. These new correlations have been developed using the larger experimental data sets now available, and the impact of the extended data on the quality of the correlations is examined in the paper. Atmospheric lifetimes have been calculated from both experimental and estimated rate coefficients to provide an overview of removal efficiencies for different classes of alkenes with respect to oxidative processes initiated by NO3, OH and O-3. A figure is presented to show the spatial scales over which alkenes may survive transport in competition with attack by NO3, OH and O-3. Removal by NO3 or OH is always more important than removal by O-3, and reactions with NO3 dominate for scales up to a few hundred metres.

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.