skip to main content


Title: Reaction mechanism, energetics, and kinetics of the water-assisted thioformaldehyde + ˙OH reaction and the fate of its product radical under tropospheric conditions
The reactions of thioformaldehyde (H 2 CS) with OH radicals and assisted by a single water molecule have been investigated using high level ab initio quantum chemistry calculations. The H 2 CS + ˙OH reaction can in principle proceed through: (1) abstraction, and (2) addition pathways. The barrier height for the addition reaction in the absence of a catalyst was found to be −0.8 kcal mol −1 , relative to the separated reactants, which has a ∼1.0 kcal mol −1 lower barrier than the abstraction channel. The H 2 CS + ˙OH reaction assisted by a single water molecule reduces the barrier heights significantly for both the addition and abstraction channels, to −5.5 and −6.7 kcal mol −1 respectively, compared to the un-catalyzed H 2 CS + ˙OH reaction. These values suggest that water lowers the barriers by ∼6.0 kcal mol −1 for both reaction paths. The rate constants for the H 2 CS⋯H 2 O + ˙OH and OH⋯H 2 O + H 2 CS bimolecular reaction channels were calculated using Canonical Variational Transition state theory (CVT) in conjunction with the Small Curvature Tunneling (SCT) method over the atmospherically relevant temperatures between 200 and 400 K. Rate constants for the H 2 CS + ˙OH reaction paths for comparison with the H 2 CS + ˙OH + H 2 O reaction in the same temperature range were also computed. The results suggest that the rate of the H 2 CS + ˙OH + H 2 O reaction is slower than that of the H 2 CS + ˙OH reaction by ∼1–4 orders of magnitude in the temperatures between 200 and 400 K. For example, at 300 K, the rates of the H 2 CS + ˙OH + H 2 O and H 2 CS + ˙OH reactions were found to be 2.2 × 10 −8 s −1 and 6.4 × 10 −6 s −1 , respectively, calculated using [OH] = 1.0 × 10 6 molecules cm −3 , and [H 2 O] = 8.2 × 10 17 molecules cm −3 (300 K, RH 100%) atmospheric conditions. Electronic structure calculations on the H 2 C(OH)S˙ product in the presence of 3 O 2 were also performed. The results show that H 2 CS is removed from the atmosphere primarily by reacting with ˙OH and O 2 to form thioformic acid, HO 2 , formaldehyde, and SO 2 as the main end products.  more » « less
Award ID(s):
1710221
NSF-PAR ID:
10197843
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
18
ISSN:
1463-9076
Page Range / eLocation ID:
10027 to 10042
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The acetylperoxy + HO 2 reaction has multiple impacts on the troposphere, with a triplet pathway leading to peracetic acid + O 2 (reaction (1a)) competing with singlet pathways leading to acetic acid + O 3 (reaction (1b)) and acetoxy + OH + O 2 (reaction (1c)). A recent experimental study has reported branching fractions for these three pathways ( α 1a , α 1b , and α 1c ) from 229 K to 294 K. We constructed a theoretical model for predicting α 1a , α 1b , and α 1c using quantum chemical and Rice–Ramsperger–Kassel–Marcus/master equation (RRKM/ME) simulations. Our main quantum chemical method was Weizmann-1 Brueckner Doubles (W1BD) theory; we combined W1BD and equation-of-motion spin-flip coupled cluster (SF) theory to treat open-shell singlet structures. Using RRKM/ME simulations that included all conformers of acetylperoxy–HO 2 pre-reactive complexes led to a 298 K triplet rate constant, k 1a = 5.11 × 10 −12 cm 3 per molecule per s, and values of α 1a in excellent agreement with experiment. Increasing the energies of all singlet structures by 0.9 kcal mol −1 led to a combined singlet rate constant, k 1b+1c = 1.20 × 10 −11 cm 3 per molecule per s, in good agreement with experiment. However, our predicted variations in α 1b and α 1c with temperature are not nearly as large as those measured, perhaps due to the inadequacy of SF theory in treating the transition structures controlling acetic acid + O 3 formation vs. acetoxy + OH + O 2 formation. 
    more » « less
  2. Abstract

    The oxidation of thiram (THR) pesticide initiated by HOradicals in water was investigated using the density functional theory (DFT) approach at the M06‐2X/6‐311++G(3df,3pd)//M06‐2X/6‐31+G(d,p) level of theory. Three oxidation mechanisms were evaluated: radical addition, hydrogen abstraction, and single electron transfer. The results showed that all these reactions are favorable and spontaneous. Among them, the addition reactions at S‐positions have the most prominent rate constants and branching ratios. In addition, the oxidation of THR is temperature‐dependent; the higher the temperature, the faster the oxidation and the higher the removal efficiency. The overall rate constant isk(T) = 7.816 × 10−6×T5.953× exp[929.0 (J mol−1)/(RT)] (M−1s−1), showing a rate constant of 6.09 × 109 m−1s−1at 298.15 K which agrees with the previously reported experimental result. In general, thiram can be easily decomposed in a wastewater treatment process with a high concentration of HO. Still, it may exist in a natural‐like environment with a very low HOconcentration from a few days to several years.

     
    more » « less
  3. Nitroxyl (HNO) and hydrogen peroxide have both been implicated in a variety of reactions relevant to environmental and physiological processes and may contribute to a unique, unexplored, pathway for the production of nitrous acid (HONO) in soil. To investigate the potential for this reaction, we report an in-depth investigation of the reaction pathway of H 2 O 2 and HNO forming HONO and water. We find the breaking of the peroxide bond and a coupled proton transfer in the first step leads to hydrogen nitryl (HNO 2 ) and an endogenous water, with an extrapolated NEVPT2 (multireference perturbation theory) barrier of 29.3 kcal mol −1 . The first transition state is shown to possess diradical character linking the far peroxide oxygen to the bridging, reacting, peroxide oxygen. The energy of this first step, when calculated using hybrid density functional theory, is shown to depend heavily on the amount of Hartree–Fock exchange in the functional, with higher amounts leading to a higher barrier and more diradical character. Additionally, high amounts of spin contamination cause CCSD(T) to significantly overestimate the TS1 barrier with a value of 36.2 kcal mol −1 when using the stable UHF wavefunction as the reference wavefunction. However, when using the restricted Hartree–Fock reference wavefunction, the TS1 CCSD(T) energy is lowered to yield a barrier of 31.2 kcal mol −1 , in much better agreement with the NEVPT2 result. The second step in the reaction is the isomerization of HNO 2 to trans -HONO through a Grotthuss-like mechanism accepting a proton from and donating a proton to the endogenous water. This new mechanism for the isomerization of HNO 2 is shown to have an NEVPT2 barrier of 23.3 kcal mol −1 , much lower than previous unimolecular estimates not including an explicit water. Finally, inclusion of an additional explicit water is shown to lower the HNO 2 isomerization barrier even further. 
    more » « less
  4. The kinetics of the combustion-relevant reaction of the allyl radical, a-C 3 H 5 , with molecular oxygen has been studied in a flow tube reactor at the vacuum ultraviolet (VUV) beamline of the Swiss Light Source storage ring, using the CRF-PEPICO (Combustion Reactions Followed by Photoelectron Photoion Coincidence Spectroscopy) setup. The ability to measure threshold photoelectron spectra enables a background-free detection of reactive species as well as an isomer-specific analysis of reaction products. Allyl was generated by direct photodissociation of allyl iodide at 266 nm and 213 nm and indirectly by the reaction of propene with Cl atoms, which were generated by photolysis from oxalyl chloride at 266 nm. Experiments were conducted at room temperature at low pressures between 0.8 and 3 mbar using Ar as the buffer gas and with excess O 2 to maintain nearly pseudo-first-order reaction conditions. Whereas allyl was detected by photoionisation using synchrotron radiation, the main reaction product allyl peroxy was not observed due to dissociative ionisation of this weakly bound species. From the concentration–time profiles of the allyl signal, second-order rate constants between 1.35 × 10 11 cm 3 mol −1 s −1 at 0.8 mbar and 1.75 × 10 11 cm 3 mol −1 s −1 at 3 mbar were determined. The rates obtained for the different allyl radical generation schemes agree well with each other, but are about a factor of 2 higher than the ones reported previously using He as a buffer gas. The discrepancy is partly attributed to the higher collision efficiency of Ar causing a varying fall-off behavior. When allyl is produced by the reaction of propene with Cl atom, an unexpected product is observed at m / z = 68, which was identified as 1,3-butadienal in the threshold photoelectron spectrum. It is formed in a secondary reaction of allyl with the OCCl radical, which is generated in the 266 nm photolysis of oxalyl chloride. 
    more » « less
  5. The formation and reactivities of [Cu–O–M] 2+ species (M = Ti–Cu, Zr–Mo and Ru–Ag) in metal-exchanged zeolites, as well as stabilities of these species towards autoreduction by O 2 elimination are investigated with density functional theory. These species were investigated in zeolite mordenite in search of insights into active site formation mechanisms, the relationship between stability and reactivity as well as discovery of heterometallic species useful for isothermal methane-to-methanol conversion (MMC). Several [Cu–O–M] 2+ species (M = Ti–Cr and Zr–Mo) are substantially more stable than [Cu 2 O] 2+ . Other [Cu–O–M] 2+ species, (M = Mn–Ni and Ru–Ag) have similar formation energies to [Cu 2 O] 2+ , to within ±10 kcal mol −1 . Interestingly, only [Cu–O–Ag] 2+ is more active for methane activation than [Cu 2 O] 2+ . [Cu–O–Ag] 2+ is however more susceptible to O 2 elimination. By considering the formation energies, autoreduction, cost and activity towards the methane C–H bond, we can only conclude that [Cu 2 O] 2+ is best suited for MMC. Formation of [Cu 2 O] 2+ is initiated by proton transfer from aquo ligands to the framework and proceeds mostly via dehydration steps. Its μ-oxo bridge is formed via water-assisted condensation of two hydroxo groups. To evaluate the relationship between [Cu 2 O] 2+ and other active sites, we also examined the formation energies of other species. The formation energies follow the trend: isolated [Cu–OH] + < paired [Cu–OH] + < [Cu 2 O] 2+ < [Cu 3 O 3 ] 2+ . Inclusion of Gibbs free-energy corrections indicates activation temperatures of 257, 307 and 327 and 331 °C for isolated [Cu–OH] + , paired [Cu–OH] + , [Cu 2 O] 2+ and [Cu 3 O 3 ] 2+ , respectively. The provocative nature of the lower-than-expected activation temperature for isolated [Cu–OH] + species is discussed. 
    more » « less