RRKM and master equation kinetic analysis of parallel addition reactions of isomeric radical intermediates in hydrocarbon flames

We have calculated the temperature-dependent rate coefficients of the addition reactions of butadien2- yl (C4H5) and acroylyl (C3H3O) radicals with ethene (C2H4), carbon monoxide (CO), formaldehyde (H2CO), hydrogen cyanide (HCN), and ketene (H2CCO), in order to explore the balance between kinetic and thermodynamic control in these combustion-related reactions. For the C4H5 radical, the 1,3-diene form of the addition products is more stable than the 1,2-diene, but the 1,2-diene form of the radical intermediate is stabilized by an allylic delocalization, which may influence the relative activation energies. For the reactions combiningC(3)H(3)OwithC(2)H(4), CO, and HCN, the opposite is true: the 1,2-enone form of the addition products is more stable than the 1,3-enone, whereas the 1,3-enone is the slightly more stable radical species. Optimized geometries and vibrational modes were computed with the QCISD/aug-cc-pVDZ level and basis, followed by single-point CCSD(T)-F12a/cc-pVDZF12 energy calculations. Our findings indicate that the kinetics in all cases favor reaction along the 1,3 pathway for both the C4H5 and C3H3O systems. The Rice-Ramsperger-Kassel-Marcus (RRKM) microcanonical rate coefficients and subsequent solution of the chemical master equation were used to predict the time-evolution of our system under conditions from 500 K to 2000 K and from 10 fi 5 bar to 10 bars. Despite the 1,3 reaction pathway being more favorable for theC(4)H(5) system, our results predict branching ratios of the 1,2 to 1,3 product as high as 0.48 at 1 bar. Similar results hold for the acroylyl system under these combustion conditions, suggesting that under kinetic control the branching of these reactions may be much more significant than the thermodynamics would suggest. This effect may be partly attributed to the low energy difference between 1,2 and 1,3 forms of the radical intermediate. No substantial pressure-dependence is found for the overall forward reaction rates until pressures decrease below 0.1 bar. Published by AIP Publishing.