Computational Chemistry and Biology Group (CCBG), DETEMA, Facultad de Química, UdelaR, Montevideo, Uruguay.
Laboratorio Universitario de Química y Contaminación Del Aire (L.U.Q.C.A.), Instituto de Investigaciones en Fisicoquímica de Córdoba (I.N.F.I.Q.C.), CONICET, Dpto. de Fisicoquímica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, Ciudad Universitaria, 5000, Córdoba, Argentina.
Chemosphere. 2024 Apr;354:141693. doi: 10.1016/j.chemosphere.2024.141693. Epub 2024 Mar 18.
This work explores theoretically the gas phase oxidation of allyl methyl sulfide (AMS, HCCHCHSCH) 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.
本工作从理论上探讨了 OH 自由基引发烯丙基甲基硫醚(AMS,HCCHCHSCH)的气相氧化,重点研究了 M06-2X-D3/aug-cc-pVTZ 和 MN15/aug-cc-pVTZ 理论水平(m06Tz 和 mn15Tz)下的 H 原子抽提途径。H 原子抽提过程涉及预反应复合物(PRC)的形成,OH 自由基有两种可能的接近方式,分别表示为“α”和“β”。PRC 表现出更高的反应活性,主要是由于硫原子和羟基氢之间的相互作用。提出了一种支持实验鉴定产物并预测形成一些含硫低挥发性产物的 H 原子抽提机制方案。在 AMS 分子中,双键加成和 H 原子抽提途径的势能面(PES)的比较表明,在 m06Tz 理论水平下,C3 上的 H 原子抽提和 C1 的加成具有几乎相同的能量曲线,而在 mn15Tz 理论水平下,最低能量通道是 C1 的加成。在考虑到在到达每个通道的过渡态之前形成 PRC 并假设反应物和 PRC 之间处于热平衡的情况下,计算了每个反应通道的理论速率常数。速率常数在 SVECV-f12/m06Tz 和 SVECV-f12/mn15Tz 理论水平上以多过渡态/多构象的方式进行了计算。SVECV-f12 方法在两个体系中的预测结果一致,仅与实验速率常数存在微小偏差。尽管由于支持 SVECV-f12 计算的 DFT 方法存在一些具体差异,但两种方法都预测 AMS+OH 气相反应中 H 原子抽提的贡献显著,这解释了实验确定的丙烯醛的高生成产率。