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 Giri, Lam K. Huynh

Research output: Contribution to journalArticlepeer-review

2 Scopus citations


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.
Original languageEnglish (US)
Pages (from-to)3730
Issue number9
StatePublished - Apr 27 2023

Bibliographical note

KAUST Repository Item: Exported on 2023-05-02
Acknowledgements: This work was funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 104.06-2020.45. We want to thank Tuyn Phan and Tri Pham (ICST) for their technical assistance.

ASJC Scopus subject areas

  • General Computer Science


Dive into the research topics of 'An Ab Initio RRKM-Based Master Equation Study for Kinetics of OH-Initiated Oxidation of 2-Methyltetrahydrofuran and Its Implications in Kinetic Modeling'. Together they form a unique fingerprint.

Cite this