title: Adiabatic quantum dynamics of CH(X2Π) + H(2S) reactions on the CH2(X̃3A″) surface and role of the excited electronic states.
abstract: We present the Born-Oppenheimer (BO) quantum mechanical (QM) dynamics of the CH decay (d) CH(X2Π) + H(2S) → C(3P) + H2(X1Σ(g)(+)) and of the H exchange reaction (e) CH(X2Π) + H′(2S) → CH′(X2Π) + H(2S) on the CH2 X̃3A″ adiabatic potential energy surface (PES) of Harding et al. (J. Phys. Chem. 1993, 97, 5472). A thorough analysis of the correlation diagram of the four lowest CH2 electronic states, as well as Renner-Teller and spin–orbit nonadiabatic test calculations on the X̃3A″, ã1A′, and b̃1A″ coupled PESs, validate the X̃3A″ BO results, confirming that these reactions occur essentially on the uncoupled X̃3A″ ground surface. We consider the CH molecule in the ground vibrational state and in the four lowest rotational states j0. Thus, we obtain initial-state resolved reaction probabilities, cross sections, and rate constants by propagating coupled-channel real wave packets and performing flux analyses. If J is the total angular momentum quantum number and K is its projection along the body-fixed z axis, CH + H gives essentially the C + H2 products via a barrierless K-inhibited insertion, CH2 resonances at low J, and large cross sections near the threshold. These cross sections decrease strongly with collision energy and depend slightly on j0. On the other hand, the small cross sections obtained for the (e) channel are nearly independent of energy. From initial-state resolved rate constants and Boltzmann populations at temperature T, we obtain QM thermal rate constants from 100 to 400 K: at 300 K, k(d) = (9.57 ± 0.96) × 10(-11) and k(e) = (1.41 ± 0.14) × 10(-11) cm(3) s(-1) for (d) and (e) reactions, respectively. The k(d) value is in good agreement with previous quasi-classical trajectory (QCT) results on the same PES, but it is larger than that observed at 297 K by a factor of 7. On the contrary, and in agreement with the small role of CH2 excited electronic states, X̃3A″ QCT and experimental rate constants agree at high temperatures. Thus, the discrepancy obtained at room T between theory and experiment should be due to an experimental error or to some theoretical effects that we have not been considered in this work. At the present state of the art, an experimental error is more likely and suggests a new measurement of the rate constant.
authors: Departament de Química Institut de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona, Spain.
authors: P
authors: P
authors: Gamallo
authors: P
authors: P
authors: Defazio
authors: S
authors: S
authors: Akpinar
authors: C
authors: C
authors: Petrongolo