Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

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 kinetics of chlorosilylene reactions II. ClSiH + C2H4: absolute rate measurements and quantum chemical and RRKM calculations for the prototype pi addition reaction

Time-resolved studies of chlorosilylene, ClSiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, are carried out to obtain rate constants for its bimolecular reaction with ethene, C2H4, in the gas-phase. The reaction is studied over the pressure range 0.13-13.3 kPa (with added SF6) at five temperatures in the range 296-562 K. The second order rate constants, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1))=(-10.55 +/- 0.10) + (3.86 +/- 0.70) kJ mol(-1)/RT ln10. The Arrhenius parameters correspond to a loose transition state and the rate constant at room temperature is 43% of that for SiH2 + C2H4, showing that the deactivating effect of Cl-for-H substitution in the silylene is not large. Quantum chemical calculations of the potential energy surface for this reaction at the G3MP2//B3LYP level show that, as well as 1-chlorosilirane, ethylchlorosilylene is a viable product. The calculations reveal how the added effect of the Cl atom on the divalent state stabilisation of ClSiH influences the course of this reaction. RRKM calculations of the reaction pressure dependence suggest that ethylchlorosilylene should be the main product. The results are compared and contrasted with those of SiH2 and SiCl2 with C2H4.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equations where the error limits are single standard deviations: log(k(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

The effect of thermal conductivity of RIM moulds in kinetics cure

Reaction Injection Moulding (RIM) is a moulding technology used for the production of large size and complex plastic parts. The RIM process is characterized essentially by the injection of a highly reactive chemical system (usually polyurethane) and fast cure, in a mould properly closed and thermally controlled. Several studies show that rapid manufacturing moulds obtained in epoxy resins for Thermoplastic Injection Moulding (TIM) affect the moulding process and the final properties of parts. The cycle time and mechanical properties of final parts are reduced, due to a low thermal conductivity of epoxy materials. In contrast, the low conductivity of materials usually applied for the rapid manufacturing of RIM moulds, increase the mechanical properties of final injected parts and reduce the cycle time. This study shows the effect of the rapid manufacturing moulds material during the RIM process. Several materials have been tested for rapid manufacturing of RIM moulds and the analysis of both, temperature profile of moulded parts during injection and the cure data experimentally obtained in a mixing and reaction cell, allow to determine and model the real effect of the mould material on the RIM process.

Atmospheric chemistry of C4F9OC2H5 (HFE-7200), C4F9OCH3 (HFE-7100), C3F7OCH3 (HFE-7000) and C3F7CH2OH: temperature dependence of the kinetics of their reactions with OH radicals, atmospheric lifetimes and global warming potentials

The atmospheric chemistry of several gases used in industrial applications, C4F9OC2H5 (HFE-7200), C4F9OCH3 (HFE-7100), C3F7OCH3 (HFE-7000) and C3F7CH2OH, has been studied. The discharge flow technique coupled with mass-spectrometric detection has been used to study the kinetics of their reactions with OH radicals as a function of temperature. The infrared spectra of the compounds have also been measured. The following Arrhenius expressions for the reactions were determined (in units of cm3 molecule-1 s-1): k(OH + HFE-7200) = (6.9+2.3-1.7) × 10-11 exp(-(2030 ± 190)/T); k(OH + HFE-7100) = (2.8+3.2-1.5) × 10-11 exp(-(2200 ± 490)/T); k(OH + HFE-7000) = (2.0+1.2-0.7) × 10-11 exp(-(2130 ± 290)/T); and k(OH + C3F7CH2OH) = (1.4+0.3-0.2) × 10-11 exp(-(1460 ± 120)/T). From the infrared spectra, radiative forcing efficiencies were determined and compared with earlier estimates in the literature. These were combined with the kinetic data to estimate 100-year time horizon global warming potentials relative to CO2 of 69, 337, 499 and 36 for HFE-7200, HFE-7100, HFE-7000 and CF3CF2CF2CH2OH, respectively.

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).

A temperature-dependent relative-rate study of the OH initiated oxidation of n-butane: the kinetics of the reactions of the 1-and 2-butoxy radicals

The kinetics of the reactions of 1-and 2-butoxy radicals have been studied using a slow-flow photochemical reactor with GC-FID detection of reactants and products. Branching ratios between decomposition, CH3CH(O-.)CH2CH3 CH3CHO + C2H5, reaction (7), and reaction with oxygen, CH3CH(O-.)CH2CH3 + O-2 -> CH3C(O)C2H5 + HO2, reaction (6), for the 2-butoxy radical and between isomerization, CH3CH2CH2CH2O. -> CH2CH2CH2CH2OH, reaction (9), and reaction with oxygen, CH3CH2CH2CH2O. + O-2 -> C3H7CHO + HO2, reaction (8), for the 1-butoxy radical were measured as a function of oxygen concentration at atmospheric pressure over the temperature range 250-318 K. Evidence for the formation of a small fraction of chemically activated alkoxy radicals generated from the photolysis of alkyl nitrite precursors and from the exothermic reaction of 2-butyl peroxy radicals with NO was observed. The temperature dependence of the rate constant ratios for a thermalized system is given by k(7)/k(6) = 5.4 x 1026 exp[(-47.4 +/- 2.8 kJ mol(-1))/RT] molecule cm(-3) and k(9)/k(8) = 1.98 x 10(23) exp[(-22.6 +/- 3.9 kJ mol(-1))/RT] molecule cm(-3). The results agree well with the available experimental literature data at ambient temperature but the temperature dependence of the rate constant ratios is weaker than in current recommendations.

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).

Thermochemical kinetics: does it still give insights?

This tutorial review revisits the subject of the seminal book written by Sidney Benson in 1968. A short summary of the nature of the subject is presented, including its place in the wider world of quantitative chemistry. A number of themes are selected to illustrate its previous and continuing usefulness in evaluating numerical values of important quantities, and probing ideas of reaction mechanism. These include strain enthalpies for biradical combination, chain reactions, why some reactions don't occur and the involvement of carbenes in hydrocarbon rearrangements.

Gas-phase kinetics of chlorosilylene reactions. I. ClSiH + Me3SiH: absolute rate measurements and theoretical calculations for prototype Si−H insertion reactions

Time-resolved studies of chlorosilylene, CISiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with trimethylsilane, Me3SiH, in the gas phase. The reaction was studied at total pressures up to 100 torr (with and without added SF6) over the temperature range 297-407 K. The rate constants were found to be pressure independent and gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-13.97 +/- 0.25) + (12.57 +/- 1.64) kJ mol(-1)/RT In 10. The Arrhenius parameters are consistent with a mechanism involving an intermediate complex, whose rearrangement is the rate-determining step. Quantum chemical calculations of the potential energy surface for this reaction and also the reactions of CISiH with SiH4 and the other methylsilanes support this conclusion. Comparisons of both experiment and theory with the analogous Si-H insertion processes of SiH2 and SiMe2 show that the main factor causing the lower reactivity of ClSiH is the secondary energy barrier. The calculations also show the existence of a novel intramolecular H-atom exchange process in the complex of ClSiH with MeSiH3.

Direct time-resolved study of the kinetics of the reaction of silylene with phenylsilane in the gas phase. Does SiH2 react with the aromatic ring?

Laser flash photolysis studies of silylene, SiH2, generated by the 193 nm laser flash photolysis phenylsilane, PhSiH3, have been carried out to obtain rate constants for its bimolecular reaction with PhSiH3 itself, in the gas phase. The reaction was studied in SF6 (mostly at 10 Torr total pressure) over the temperature range 298-595 K. The rate constants (also found to be pressure independent) gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-9.92 +/- 0.04) + (3.31 +/- 0.27) kJ mol(-1)/RT ln 10 Similar investigations of the reaction of silylene with benzene, C6H6, (295-410 K) gave data suggestive of the fact that SiH2 might be reacting with photochemical products of C6H6 as well as with C6H6 itself. However, in the latter system, apparent rate constants were sufficiently low to indicate that in the reaction of SiH2 with PhSiH3 addition to the aromatic ring was unlikely to be in excess of 3% of the total. Quantum chemical calculations of the energy surface for SiH2 + C6H6 indicate that 7-silanorcaradiene and 7-silacycloheptatriene are possible products but that PhSiH3 formation is unlikely. RRKM calculations suggest that 7-silanorcaradiene should be the initial product but that it cannot be collisionally stabilized under experimental conditions

A kinetics and mechanistic study of the atmospherically relevant reaction between molecular chlorine and dimethyl sulfide (DMS)

A gas-phase kinetics study of the atmospherically important reaction between Cl2 and dimethyl sulfide (DMS)Cl2 + CH3SCH3 → productshas been made using a flow-tube interfaced to a photoelectron spectrometer. The rate constant for this reaction has been measured at 1.6 and 3.0 torr at T = (294 ± 2) K as (3.4 ± 0.7) × 10−14 cm3 molecule−1 s−1. Reaction (1) has been found to proceed via an intermediate, (CH3)2SCl2, to give CH3SCH2Cl and HCl as the products. The mechanism of this reaction and the structure of the intermediate were investigated using electronic structure calculations. A comparison of the mechanisms of the reactions between Cl atoms and DMS, and Cl2 and DMS has been made and the relevance of the results to atmospheric chemistry is discussed.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.