Theoretical investigation on isomerization and decomposition reactions of pentanol radicals–part II: linear pentanol isomers

Abstract

High-level ab initio calculations are conducted for studying the kinetics of three linear pentanol radicals generated through H-atom abstraction reactions. The species involved are optimized using the M06-2X/6-311++G(d,p) level of theory, while a relaxed scan at the M06-2X/6-31g level of theory with 10° increments is used for the hindrance potential for low-frequency torsional modes. Single-point energies for all stationary points are obtained through the QCISD(T) and MP2 methods in combination with cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets, which can be extrapolated to the complete basis set (CBS) limit. The rate constants and branching ratios for isomerization and decomposition reactions are computed over a temperature range of 250–2000 K and a pressure range of 0.01–100 atm. Isomerization reactions are dominant at low temperatures, while decomposition reactions are more dominant at high temperatures. The branching ratio of the isomerization reaction exhibits a slight decrease with increasing pressure, while the trend for decomposition reactions depends on the type of the breaking bond. Based on the calculations for five branched pentanol radicals in part I, kinetics of linear and branched pentanol radicals are compared in this work and the results reveal that, for the same kind of β-scission reaction at similar positions of linear and branched pentanol radicals, the rate constants of branched ones are faster than those of linear ones at low temperatures. The hydroxyl group adjacent to the breaking bond can increase the β-scission reaction rate constants, while the effect can be ignored when the hydroxyl group is not adjacent to the breaking bond. Moreover, compared to when the hydroxyl group is located in the middle of the carbon chain, its positioning at the chain's end yields a more noticeable impact on the products and rate constants of C–O bond and O–H bond β-scission reactions. Besides, when incorporating calculated rate constants into the CRECK model, the updated mechanism shows a better performance for ignition delay times of 1-pentanol in the NTC range but exhibits lower reactivity at higher temperatures. The simulation of speciation profiles also shows better agreement with the experimental data obtained using a flow reactor.