Reactions of Cl Atoms with Dimethyl Sulfide: A Theoretical Calculation and an Experimental Study with Cavity Ring-Down Spectroscopy

Atmospheric Gas-Phase Formation of Methanesulfonic Acid

Despite its impact on the climate, the mechanism of methanesulfonic acid (MSA) formation in the oxidation of dimethyl sulfide (DMS) remains unclear. The DMS + OH reaction is known to form methanesulfinic acid (MSIA), methane sulfenic acid (MSEA), the methylthio radical (CH3S), and hydroperoxymethyl thioformate (HPMTF). Among them, HPMTF reacts further to form SO2 and OCS, while the other three form the CH3SO2 radical. Based on theoretical calculations, we find that the CH3SO2 radical can add O2 to form CH3S(O)2OO, which can react further to form MSA. The branching ratio is highly temperature sensitive, and the MSA yield increases with decreasing temperature. In warmer regions, SO2 is the dominant product of DMS oxidation, while in colder regions, large amounts of MSA can form. Global modeling indicates that the proposed temperature-sensitive MSA formation mechanism leads to a substantial increase in the simulated global atmospheric MSA formation and burden.

Measurement of the Intramolecular Hydrogen-Shift Rate Coefficient for the CH3SCH2OO Radical between 314 and 433 K

The intramolecular hydrogen-shift rate coefficient of the CH₃SCH₂O₂ (methylthiomethylperoxy, MSP) radical, a product formed in the oxidation of dimethyl sulfide (DMS), was measured using a pulsed laser photolysis flow tube reactor coupled to a high-resolution time-of-flight chemical ionization mass spectrometer that measured the formation of the DMS degradation end product HOOCH₂SCHO (hydroperoxymethyl thioformate). Measurements performed over the temperature range of 314-433 K yielded a hydrogen-shift rate coefficient of k₁(T) = (2.39 ± 0.7) × 10⁹ exp(-(7278 ± 99)/T) s^-1 Arrhenius expression and a value extrapolated to 298 K of 0.06 s^-1. The potential energy surface and the rate coefficient have also been theoretically investigated using density functional theory at the M06-2X/aug-cc-pVTZ level combined with approximate CCSD(T)/CBS energies yielding k₁(273-433 K) = 2.4 × 10¹¹ × exp(-8782/T) s^-1 and k₁(298 K) = 0.037 s^-1 in fair agreement with the experimental results. Present results are compared with the previously reported values of k₁(293-298 K).

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) CH₃SCH₂ + HCl, (B) CH₃S + CH₃Cl, (C) CH₃SCl + CH₃, and (C') CH₃S(Cl)CH₃. 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 CH₃SCH₂ + 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·H₂O + Cl, (b) Cl·H₂O + DMS and (c) DMS·Cl + H₂O 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.

A Gas-Phase Kinetic Study of the Reaction between Bromine Monoxide and Methylperoxy Radicals at Atmospheric Temperatures

The rate constant of the reaction of BrO with CH(3)O(2) was determined to be k1 = (6.2 +/- 2.5) x 10(-12) cm3 molecule(-1) s(-1) at 298 K and 100-200 Torr of O2 diluent. Quoted uncertainty was two standard deviations. No significant pressure dependence of the rate constants was observed at 100-200 Torr total pressure of N2 or O2 diluents. Temperature dependence of the rate constants was further investigated over the range 233-333 K, and an Arrhenius type expression was obtained for k1 = 4.6 x 10(-13) exp[(798 +/- 76)/T] cm3 molecule(-1) s(-1). The product branching ratios were evaluated and the atmospheric implications were discussed.

Rotational spectrum and inversion motions in the neon-dimethyl sulfide complex

The rotational spectra of the (20)Ne and (22)Ne isotopomers of the Ne-dimethyl sulfide (DMS) rare gas dimer have been measured by Fourier transform microwave spectroscopy. MP2/6-311++G(2d,2p) calculations, and the experimental spectroscopic data, suggest a structure of C(s) symmetry in which the Ne atom lies above the heavy atom plane of the DMS (in the sigma(v) plane which bisects the CSC angle). Experimental rotational constants are consistent with a S...Ne distance of 3.943(6) Angstroms and a (cm...S...Ne) angle of 63.2(6) degrees (where cm is the center of mass of DMS). A motion of the Ne atom from one side of the DMS to the other gives rise to inversion splittings of around 3 MHz in the c-type transitions. An ab initio potential energy surface calculation has allowed examination of several possible tunneling pathways, and suggests a barrier of between 20 and 40 cm(-1) for the inversion motion, depending on the tunneling pathway taken by the Ne. Dipole moment measurements are consistent with both the experimental and ab initio structures.

Kinetics, Mechanism, and Thermochemistry of the Gas Phase Reaction of Atomic Chlorine with Dimethyl Sulfoxide

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the reaction of chlorine atoms with dimethyl sulfoxide (CH3S(O)CH3; DMSO) as a function of temperature (270-571 K) and pressure (5-500 Torr) in nitrogen bath gas. At T = 296 K and P > or = 5 Torr, measured rate coefficients increase with increasing pressure. Combining our data with literature values for low-pressure rate coefficients (0.5-3 Torr He) leads to a rate coefficient for the pressure independent H-transfer channel of k1a = 1.45 x 10(-11) cm3 molecule(-1) s(-1) and the following falloff parameters for the pressure-dependent addition channel in N2 bath gas: k(1b,0) = 2.53 x 10(-28) cm6 molecule(-2) s(-1); k(1b,infinity) = 1.17 x 10(-10) cm3 molecule(-1) s(-1), F(c) = 0.503. At the 95% confidence level, both k1a and k1b(P) have estimated accuracies of +/-30%. At T > 430 K, where adduct decomposition is fast enough that only the H-transfer pathway is important, measured rate coefficients are independent of pressure (30-100 Torr N2) and increase with increasing temperature. The following Arrhenius expression adequately describes the temperature dependence of the rate coefficients measured at over the range 438-571 K: k1a = (4.6 +/- 0.4) x 10(-11) exp[-(472 +/- 40)/T) cm3 molecule(-1) s(-1) (uncertainties are 2sigma, precision only). When our data at T > 430 K are combined with values for k1a at temperatures of 273-335 K that are obtained by correcting reported low-pressure rate coefficients from discharge flow studies to remove the contribution from the pressure-dependent channel, the following modified Arrhenius expression best describes the derived temperature dependence: k1a = 1.34 x 10(-15)T(1.40) exp(+383/T) cm3 molecule(-1) s(-1) (273 K < or = T < or = 571 K). At temperatures around 330 K, reversible addition is observed, thus allowing equilibrium constants for Cl-DMSO formation and dissociation to be determined. A third-law analysis of the equilibrium data using structural information obtained from electronic structure calculations leads to the following thermochemical parameters for the association reaction: delta(r)H(o)298 = -72.8 +/- 2.9 kJ mol(-1), deltaH(o)0 = -71.5 +/- 3.3 kJ mol(-1), and delta(r)S(o)298 = -110.6 +/- 4.0 J K(-1) mol(-1). In conjunction with standard enthalpies of formation of Cl and DMSO taken from the literature, the above values for delta(r)H(o) lead to the following values for the standard enthalpy of formation of Cl-DMSO: delta(f)H(o)298 = -102.7 +/- 4.9 kJ mol(-1) and delta(r)H(o)0 = -84.4 +/- 5.8 kJ mol(-1). Uncertainties in the above thermochemical parameters represent estimated accuracy at the 95% confidence level. In agreement with one published theoretical study, electronic structure calculations using density functional theory and G3B3 theory reproduce the experimental adduct bond strength quite well.

Observation of Adducts in the Reaction of Cl Atoms with XCH2I (X = H, CH3, Cl, Br, I) Using Cavity Ring-Down Spectroscopy

The reactions of Cl atoms with XCH2I (X = H, CH3, Cl, Br, I) have been studied using cavity ring-down spectroscopy in 25-125 Torr total pressure of N2 diluent at 250 K. Formation of the XCH2I-Cl adduct is the dominant channel in all reactions. The visible absorption spectrum of the XCH2I-Cl adduct was recorded at 405-632 nm. Absorption cross-sections at 435 nm are as follows (in units of 10(-18) cm2 molecule(-1)): 12 for CH3I, 21 for CH3CH2I, 3.7 for CH2ICl, 7.1 for CH2IBr, and 3.7 for CH2I2. Rate constants for the reaction of Cl with CH3I were determined from rise profiles of the CH3I-Cl adduct. k(Cl + CH3I) increases from (0.4 +/- 0.1) x 10(-11) at 25 Torr to (2.0 +/- 0.3) x 10(-11) cm3 molecule(-1) s(-1) at 125 Torr of N2 diluent. There is no discernible reaction of the CH3I-Cl adduct with 5-10 Torr of O2. Evidence for the formation of an adduct following the reaction of Cl atoms with CF3I and CH3Br was sought but not found. Absorption attributable to the formation of the XCH2I-Cl adduct following the reaction of Cl atoms with XCH2I (X = H, CH3, Br, I) was measured as a function of temperature over the range 250-320 K.