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. 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
  3. Using a combination of two-dimensional infrared (2D IR) and variable temperature Fourier transform infrared (FTIR) spectroscopies the rapid structural isomerization of a five-coordinate ruthenium complex is investigated. In methylene chloride, three exchanging isomers were observed: (1) square pyramidal equatorial, ( 1 ); (2) trigonal bipyramidal, ( 0 ); and (3) square pyramidal apical, ( 2 ). Exchange between 1 and 0 was found to be an endergonic process (Δ H = 0.84 (0.08) kcal mol −1 , Δ S = 0.6 (0.4) eu) with an isomerization time constant of 4.3 (1.5) picoseconds (ps, 10 −12 s). Exchange between 0 and 2 however was found to be exergonic (Δ H = −2.18 (0.06) kcal mol −1 , Δ S = −5.3 (0.3) eu) and rate limiting with an isomerization time constant of 6.3 (1.6) ps. The trigonal bipyramidal complex was found to be an intermediate, with an activation barrier of 2.2 (0.2) kcal mol −1 and 2.4 (0.2) kcal mol −1 relative to the equatorial and apical square pyramidal isomers respectively. This study provides direct validation of the mechanism of Berry pseudorotation – the pairwise exchange of ligands in a five-coordinate complex – a process that was first described over fifty years ago. This study also clearly demonstrates that the rate of pseudorotation approaches the frequency of molecular vibrations. 
    more » « less
  4. Abstract

    Reaction of {LiC6H2−2,4,6‐Cyp3⋅Et2O}2(Cyp=cyclopentyl) (1) of the new dispersion energy donor (DED) ligand, 2,4,6‐triscyclopentylphenyl with SnCl2afforded a mixture of the distannene {Sn(C6H2−2,4,6‐Cyp3)2}2(2), and the cyclotristannane {Sn(C6H2−2,4,6‐Cyp3)2}3(3).2is favored in solution at higher temperature (345 K or above) whereas3is preferred near 298 K. Van't Hoff analysis revealed the3to2conversion has a ΔH=33.36 kcal mol−1and ΔS=0.102 kcal mol−1 K−1, which gives a ΔG300 K=+2.86 kcal mol−1, showing that the conversion of3to2is an endergonic process. Computational studies show that DED stabilization in3is −28.5 kcal mol−1per {Sn(C6H2−2,4,6‐Cyp3)2unit, which exceeds the DED energy in2of −16.3 kcal mol−1per unit. The data clearly show that dispersion interactions are the main arbiter of the3to2equilibrium. Both2and3possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results).

     
    more » « less
  5. Abstract

    Reaction of {LiC6H2−2,4,6‐Cyp3⋅Et2O}2(Cyp=cyclopentyl) (1) of the new dispersion energy donor (DED) ligand, 2,4,6‐triscyclopentylphenyl with SnCl2afforded a mixture of the distannene {Sn(C6H2−2,4,6‐Cyp3)2}2(2), and the cyclotristannane {Sn(C6H2−2,4,6‐Cyp3)2}3(3).2is favored in solution at higher temperature (345 K or above) whereas3is preferred near 298 K. Van't Hoff analysis revealed the3to2conversion has a ΔH=33.36 kcal mol−1and ΔS=0.102 kcal mol−1 K−1, which gives a ΔG300 K=+2.86 kcal mol−1, showing that the conversion of3to2is an endergonic process. Computational studies show that DED stabilization in3is −28.5 kcal mol−1per {Sn(C6H2−2,4,6‐Cyp3)2unit, which exceeds the DED energy in2of −16.3 kcal mol−1per unit. The data clearly show that dispersion interactions are the main arbiter of the3to2equilibrium. Both2and3possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results).

     
    more » « less