Unexpected high yield of acrolein underlies the importance of the hydrogen-abstraction mechanism in photooxidation of allyl methyl sulfide (AMS)

This work explores theoretically the gas phase oxidation of allyl methyl sulfide (AMS, H2CCHCH2SCH3) initiated by •OH radicals, focusing on the H-abstraction pathway at the M06-2X-D3/aug-cc-pVTZ and MN15/aug-cc-pVTZ levels of theory (m06Tz and mn15Tz). The formation of a prereactive complex (PRC) is involved in H-abstraction processes with two potential directions of approach for the OH radical, denoted as "α" and "β". The PRCs, demonstrate increased reactivity, primarily due to the interaction between the sulfur atoms and the hydroxyl hydrogen. A scheme for the H-abstraction mechanism that supports the experimentally identified products and predicts the formation of some S-containing low volatility products is proposed. The comparison of the potential energy surface (PES) between the double bond addition and H-abstraction paths in the AMS molecule shows that at the m06Tz level of theory, the H-abstraction on C3 and the addition to C1 have nearly the same profile of energy, while at the mn15Tz level, the minimum energy channel is the addition to C1. The theoretical rate coefficient for each reaction channel was calculated, considering the formation of a PRC prior to reaching the transition state of each channel and assuming thermal equilibrium between reactants and the PRC. The rate constants were calculated in a multi-TS/multi-conformer way at the SVECV-f12/m06Tz and SVECV-f12/mn15Tz levels of theory. The SVECV-f12 method is consistent in its predictions in both systems and exhibits only minor deviations from the experimental rate constants. Despite some specific differences due to the DFT method supporting the SVECV-f12 calculations, both methodologies predict a significant H-abstraction contribution in the AMS + OH gas phase reaction, which explains the high formation yield for acrolein determined experimentally.

Dark transformation from 17β-estradiol to estrone initiated by hydroxyl radical in dissolved organic matter

The occurrence and fate of 17β-estradiol (E₂) in natural water have gained extensive attention owing to its high ecotoxic risk to wildlife. Dissolved organic matter (DOM) is a ubiquitous water constituent and contributes significantly to E₂ removal, although the reaction mechanism is rarely clarified. The present study aims to investigate E₂ transformation in water containing fresh or aged DOM surrogates at environmentally relevant concentrations in the dark. Experiments along with radical probes of benzene and furfuryl alcohol reveal that reactive radicals, particularly hydroxyl radical (·OH), formed non-photochemically at higher concentrations in aged DOM than in fresh DOM. The contribution of ·OH in E₂ removal is indicated by the decreases in the removal of radical probes in the presence of E₂; moreover, E₂ removal is inhibited in the presence of radical scavengers. The dose-dependent inhibitive effect of substrate concentrations, including E₂ and coexistent propylparaben, shows that the radical concentration is a limiting factor for E₂ removal, which could be enhanced by increasing DOM concentration, dissolved oxygen, and light supply. As the main byproduct, estrone (E₁) is persistent in the current DOM water in the dark, but it can be easily photodegraded when exposed to light. Theoretical analysis reveals that the initial step is ·OH-initiated H- abstraction on the hydroxyl group in the cyclopentane ring of E₂. The formed singlet excited state of E₂ undergoes further intramolecular rearrangement and oxidative dehydrogenation to generate E₁ and the hydroperoxy radical (·HO₂). Considering the universal occurrence of E₂ in DOM-rich aquatic matrices, the present findings have special implications for the biogeochemical cycle and risk assessment of this pollutant in natural aquatic environments, particularly those beyond the photic zone.

Computational investigations on kinetics of reaction between t-butanol and OH radical and ozone formation potential

The kinetics of the gas-phase atmospheric reaction of t-butanol with OH radicals is computationally studied using the CCSD(T)/aug-cc-pVTZ//M06-2X/6-311++G(d,p) level of calculation. The rate coefficients are evaluated for a wide temperature range of 250-1200 K and the calculated rate coefficient value of 0.83×10^-12cm³molecule^-1s^-1 at 298K is in close agreement with experimental results. The H-abstraction from the -CH₃ group is predicted to be the main reaction channel. A weak negative temperature dependence of rate coefficient is observed in 250-300 K. The study also highlighted the possibility of re-generation of OH radicals at higher temperature. The ozone formation potential of t-butanol in the troposphere has also been estimated and discussed.

Hydroxyl radical oxidation of cyclic methylsiloxanes D4 ∼ D6 in aqueous phase

Cyclic methylsiloxanes (CMS) were listed as candidates of substances of very high concerns in 2018 by the REACH. These compounds can enter environmental waters, and potentially cause harmful effects to aquatic organisms and human beings. Until now, reaction mechanisms of these pollutants with hydroxyl radicals (HO) in aqueous phase were unknown. In this study, reaction mechanisms of three typical CMS (D4 ∼ D6) with HO in aqueous phase were investigated by employing both UV/H₂O₂ experiments and density functional theoretical calculations. Bimolecular reaction rate constants (k_HO·) of D4 ∼ D6 with HO were determined as k_HO·(D4) = 8k_HO·(D5) = 12k_HO·(D6) = 6.6 × 10⁸ L mol^-1 s^-1. Half-lives of HO oxiding D4 ∼ D6 ranged from 12 to 140 days at [HO] = 10^-15 mol L^-1 in sunlit surface water, and were comparable to (D4, D5) or much shorter (D6) than hydrolytic half-lives. The reactivity to HO decreased with the increasing size of siloxane ring in aqueous phase, in an order totally opposite to that in gaseous phase. Calculation results indicated that HO oxidation of the three CMS proceeded spontaneously through an exothermic H-abstraction process at the first step. Water molecules participated into H-abstraction of CMS and caused energy barrier of D5 higher than that of D4. Thus, H-bonds formed by water molecules were responsible for the reverse reactivity of CMS in aqueous phase. This work provided basic evidences suggesting environmental persistence of CMS in aqueous phase completely different from that in gaseous phase.

Atmospheric oxidation of HFE-7300 [n-C2F5CF(OCH3)CF(CF3)2] initiated by •OH/Cl oxidants and subsequent degradation of its product radical: a DFT approach

To understand the atmospheric chemistry of hydrofluoroethers, we have studied the oxidation of a highly fluorinated compound n-C₂F₅CF(OCH₃)CF(CF₃)₂ (HFE-7300) by OH/Cl oxidants. Here, we have employed M06-2X functional along with a 6-31 + G(d,p) basis set to obtain the optimized structures, various forms of energies, and different modes of frequencies for all species. We have characterized energies of all species on the potential energy surface, and it indicates that H-abstraction from n-C₂F₅CF(OCH₃)CF(CF₃)₂ by Cl atom is kinetically more dominant than the H-abstraction reaction initiated by OH radical. In contrast, the calculated energy change (Δ_rH°₂₉₈ and Δ_rG°₂₉₈) results govern that OH-initiated H-abstraction reaction is highly exothermic and spontaneous compared to the Cl-initiated H-abstraction reaction. Rate constants are estimated using transition state theory as well as canonical variation transition state theory at the temperature range 200-1000 K and 1 atm pressure. The calculated rate constants of the H-abstraction channels are found to be in good agreement with the reported experimental rate constant at 298 K. Moreover, we have estimated the atmospheric lifetimes of HFE-7300 for the reaction with OH radical and Cl atom and are found to be 1.75 and 153.93 years, respectively. Additionally, the global warming potentials for HFE-7300 molecule are also estimated for 20-, 100-, and 500-year time horizons. Further, subsequent aerial oxidation of product radical (n-C₂F₅CF(OCH₂)CF(CF₃)₂) in the presence of NO radical is performed, and it produced alkoxy radical via formation of peroxy radical. This alkoxy radical undergoes unimolecular decompositions via two different ways and formed n-C₂F₅CF(OCHO)CF(CF₃)₂ and n-C₂F₅CF(OH) CF(CF₃)₂ products.

Reaction of 1-propanol with Ozone in Aqueous Media

The main aim of this work is to substantiate the mechanism of 1-propanol oxidation by ozone in aqueous solution when the substrate is present in large excess. Further goals are assessment of the products, their formation yields as well as the kinetic parameters of the considered reaction. The reaction of ozone with 1-propanol in aqueous solution occurs via hydride transfer, H-abstraction and insertion. Of these three mechanisms, the largest share is for hydride transfer. This implies the extraction of an hydride ion from the activated C-H group by O₃ according to reaction: (C₂H₅)(H)(HO)C-H + O₃ → [(C₂H₅)(H)(HO)C⁺ + HO₃^-]_cage → (C₂H₅)(H)(HO)C⁺ + HO₃^-. The experimentally determined products and their overall formation yields with respect to ozone are: propionaldehyde-(60 ± 3)%, propionic acid-(27.4 ± 1.0)%, acetaldehyde-(4.9 ± 0.3)%, acetic acid-(0.3 ± 0.1)%, formaldehyde-(1.0 ± 0.1)%, formic acid-(4.6 ± 0.3)%, hydrogen peroxide-(11.1 ± 0.3)% and hydroxyl radical-(9.8 ± 0.3)%. The reaction of ozone with 1-propanol in aqueous media follows a second order kinetics with a reaction rate constant of (0.64 ± 0.02) M^-1·s^-1 at pH = 7 and 23 °C. The dependence of the second order rate constant on temperature is described by the equation: l n   k I I = ( 27.17 ± 0.38 ) - ( 8180 ± 120 ) × T - 1 , which gives the activation energy, Ea = (68 ± 1) kJ mol^-1 and pre-exponential factor, A = (6.3 ± 2.4) × 10¹¹ M^-1 s^-1. The nature of products, their yields and the kinetic data can be used in water treatment. The fact that the hydride transfer is the main pathway in the 1-propanol/ozone system can probably be transferred on other systems in which the substrate is characterized by C-H active sites only.

Application of one-dimensional semiclassical transition state theory to the CH3OH + H ⇌ CH2OH/CH3O + H2 reactions

The rate constants of the two branches of H-abstractions from CH₃OH by the H-atom and the corresponding reactions in the reverse direction are calculated using the one-dimensional semiclassical transition state theory (1D SCTST). In this method, only the reaction mode vibration of the transition state (TS) is treated anharmonically, while the remaining internal degrees of freedom are treated as they would have been in a standard TS theory calculation. A total of eight ab initio single-point energy calculations are performed in addition to the computational cost of a standard TS theory calculation. This allows a second-order Richardson extrapolation method to be employed to improve the numerical estimation of the third- and fourth-order derivatives, which in turn are used in the calculation of the anharmonic constant. Hindered-rotor (HR) vibrations are identified in the equilibrium states of CH₃OH and CH₂OH, and the TSs of the reactions. The partition function of the HRs are calculated using both a simple harmonic oscillator model and a more sophisticated one-dimensional torsional eigenvalue summation (1D TES) method. The 1D TES method can be easily adapted in 1D SCTST computation. The resulting 1D SCTST with 1D TES rate constants show good agreement to previous theoretical and experimental works. The effects of the HR on rate constants for different reactions are also investigated.This article is part of the theme issue 'Modern theoretical chemistry'.

Solvent Degradation in Nonaqueous Li-O2 Batteries: Oxidative Stability versus H-Abstraction

Developing rechargeable Li-O2 batteries hinges on identifying stable solvents resistant to decomposition. Here, we focus on solvent stability against adsorption-induced H-abstraction during discharge. Using a detailed thermodynamic analysis, we show that a solvent's propensity to resist H-abstraction is determined by its acid dissociation constant, pKa, in its own environment. Upon surveying hundreds of solvents for their pKa values in different media, we find linear correlations between the pKa values across various classes of solvents in any two given media. Utilizing these correlations, we choose DMSO as the common standard to compare the relative stability trends. We construct a stability plot based on the solvent's HOMO level and its pKa in DMSO, which reveals that most solvents obey a correlation where solvents with lower HOMO levels tend to have lower pKa values in DMSO. However, this is at odds with the stability requirement that demands deep HOMO levels and high pKa values. Thus, stable solvents need to be outliers to this observed correlation.