Born-Oppenheimer and Renner-Teller Quantum Dynamics of CH(X(2)Π) + D((2)S) Reactions on Three CHD Potential Surfaces.

Abstract

The quantum dynamics of three CH(X(2)Π) + D((2)S) reactions is studied by means of the coupled-channel time-dependent real-wavepacket (WP) and flux methods at collision energy Ecol ≤ 0.6 eV and on three potential energy surfaces (PESs): the Born-Oppenheimer (BO) ground PES X̃(3)A″ and the excited ones ã(1)A' and b̃(1)A″, coupled by nonadiabatic (NA) Renner-Teller (RT) effects. This three-state model is suitable for obtaining initial-state-resolved observables, is based on a complete analysis of the correlation diagram of the lowest electronic states of the CHD intermediate and of their NA interactions, and neglects the smaller coupling effects due to the asymptotic electronic angular momenta that become important in state-to-state dynamics. WPs are propagated on each PES at total angular momentum values J ≤ 70, with CH in the two lowest vibrational states v0 and in the ground rotational state j0 = 1. Reaction probabilities are obtained for three possible final products (f): (dP) CH decay and C((3)P) + HD(X(1)Σ(+)) formation that occurs on the uncoupled ground PES, (dD) CH decay and C((1)D) + HD(X(1)Σ(+)) formation that depends on the RT-coupled singlet species, and (ex) exchange to CD(X(2)Π) + H((2)S) available adiabatically from the X̃(3)A″ PES and nonadiabatically from ã(1)A' and b̃(1)A″. Observable cross sections σf,v0j0 and rate constants kf,v0j0 in the temperature range T = 100-2000 K are obtained for (dP), (dD), and (ex) channels. Comparing BO with RT probabilities, we show that NA effects are important at high J values for the (ex) channel at v0 = 1. Real time mechanisms on the three PESs show that RT couplings are opened after some time and clearly point out the formation of the product channels. Both cross sections and rate constants present the same sequence, for example σex,11 > σdP,01 ∼ σex,01 > σdP,11 ≫ σdD,11 ≫ σdD,01, and the CH vibrational excitation enhances the total removal CH+D reactivity by a factor of ∼1.7, mainly due to the increase of the (ex) channel contribution from ∼47% at v0 = 0 to ∼76% at v0 = 1. This fact implies a considerable vibrational enhancement of combustion processes at high temperature. In agreement with the probability results, the ã(1)A'/b̃(1)A″ RT coupling increases both σex,11 and kex,11 up to ∼30%. Moreover, including the three PESs in the dynamics simulation of CH+D increase by far the (ex)/(dP) branching ratio with respect to the CH + H' reaction. Thus, at room temperature, kdP,01 changes from 10.8 × 10(-11) to 3.4 × 10(-11) cm(3) s(-1) substituting H atom by D.