An ab initio RRKM-based master equation study for kinetics of OH-initiated oxidation of 2-methyltetrahydrofuran and its implications in kinetic modeling
Tam V-T Mai, Thanh Q Bui, Nguyen Thi Ai Nhung, Phan Tu Quy, Krishna Prasad Shrestha, Fabian Mauss, Binod Raj Giri, Lam K Huynh
Abstract
Cyclic ethers (CEs) can be promising future biofuel candidates. Most CEs possess physico-chemical and combustion indicators comparable to conventional fuels, making them suitable for internal combustion engines. This work computationally investigates the kinetic behaviors of hydrogen abstraction from 2-methyl tetrahydrofuran (2MTHF), one of the promising CEs, by hydroxyl radicals under combustion and atmospheric relevant conditions. The various reaction pathways were explored using the CCSD(T)/cc-pVTZ//M06-2X/aug-cc-pVTZ level of theory. The Rice–Ramsperger–Kassel–Marcus-based master equation (RRKM-ME) rate model, including treatments for hindered internal rotation and tunneling, was employed to describe time-dependent species profiles and pressure and temperature-dependent rate coefficients. Our kinetic model revealed that the H-abstraction proceeds via an addition-elimination mechanism forming reaction complexes at both the entrance and exit channels. Eight different reaction channels yielding five radical products were located. The reaction exhibited complex kinetics yielding a U-shaped Arrhenius behavior. An unusual occurrence of negative temperature dependence was observed at low temperatures, owing to the negative barrier height for the hydrogen abstraction reaction from the C-H bond at the vicinity of the O-atom. A shift in the reaction mechanism was observed with the dominance of the abstraction at Cα-H of 2MTHF ring (causing negative-T dependence) and at CH3 (positive-T dependence) at low and high temperatures, respectively. Interestingly, the pressure effect was observed at low temperatures, revealing the kinetic significance of the pre-reaction complex. Under atmospheric pressure, our theoretical rate coefficients showed excellent agreement with the available literature data. Our model nicely captured the negative temperature-dependent behaviors at low temperatures. Our predicted global rate coefficients can be expressed as k (T, 760 Torr) = 3.55 × 101 × T−4.72 × exp [−340.0 K/T] + 8.21 × 10−23 × T3.49 × exp [918.8 K/T] (cm3/molecule/s). Our work provides a detailed kinetic picture of the OH-initiated oxidation kinetics of 2MTHF. Hence, this information is useful for building a kinetic me chanism for methylated cyclic ethers.