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 butadien-2-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 combining C3H3O with C2H4, 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-pVDZ-F12 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-5 bar to 10 bars. Despite the 1,3 reaction pathway being more favorable for the C4H5 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.