Theoretical study of the oxidation mechanisms of thiophene initiated by hydroxyl radicals

Abstract

The mechanisms for the oxidation of thiophene by OH radicals under inert conditions (Ar) have been studied using density functional theory in conjunction with various exchange-correlation functionals. These results were compared with benchmark CBS-QB3 theoretical results. Kinetic rate constants were estimated by means of variational transition state theory (VTST) and the statistical Rice–Ramsperger–Kassel–Marcus (RRKM) theory. Effective rate constants were calculated via a steady-state analysis based upon a two-step model reaction mechanism. In line with experimental results, the computed branching ratios indicate that the most kinetically efficient process involves OH addition to a carbon atom adjacent to the sulfur atom. Due to the presence of negative activation energies, pressures larger than 104 bar are required to reach the high-pressure limit. Nucleus-independent chemical shift indices and natural bond orbital analysis show that the computed activation energies are dictated by changes in aromaticity and charge-transfer effects due to the delocalization of lone pairs from sulfur to empty π* orbitals.