Nonadiabatic Renner-Teller quantum dynamics of OH(X2?) + H+ reactive collisions

Abstract

Following previous studies on the O(P-3) + H2(+)(X-2 Sigma(+)(g)) collisions, we present the nonadiabatic quantum dynamics of the reactions OH(X-2 Pi) + H'(+) -> OH'(X2P) + H+, exchange (e), -> OH+(X-3 Sigma(-)) + H'(S-2), quenching (q), and -> OH' + (X-3 Sigma(-)) + H(2S), exchange-quenching (eq). The reactants and products correlate via the ground X(2)A and first excited A A(2)A' electronic states of OH2+, which are the degenerate components of linear (2)Pi species. Therefore, they are strongly perturbed by nonadiabatic Renner-Teller (RT) effects, opening the (q) and (eq) channels that are closed in the Born-Oppenheimer approximation. Using accurate potential energy surfaces (PESs) and RT matrix elements, initial-state-resolved reaction probabilities, real-time dynamics, cross sections, and rate constants of the product channels are obtained through the time-dependent real wavepacket (WP) method and full coupled-channel calculations. Owing to the nonadiabatic couplings, the WP jumps from the excited A A(2)A' surface to the A(2)A' ground PES, avoiding any barrier, opening the quenching channels, and giving many collision complexes into the deep minima of both PESs, as it is clearly shown by the oscillations of the reaction probabilities and by the time-dependent WP dynamics. All the results show that the nonadiabatic-RT channels (q) and (eq) are highly reactive, much more than the adiabatic one (e), pointing out large RT effects. The reactivity of the quenching channels is similar, accounting for 97% of the overall reactivity. In fact, the maximum values of the (q) and (eq) cross sections sigma(q) and sigma(eq) are equal to 31.6 A(2), whereas the maximum se value equals 1.34 A(2), and the maximum values of the rate constants k(q), k(eq), and k(e) are 2.07 Chi 10 Chi 10, 2.45 Chi 10 Chi 10, and 0.23 Chi 10 Chi 10 cm(3) s(-1). Some calculations show that the centrifugalsudden and the truncated coupled-channel approximations cannot be employed for the (q) channel. After a sharp increase at the threshold, sq and seq decrease at larger collision energies while se and the rate constants depend slightly on the collision energy and temperature, respectively. These findings are consistent with the barrierless nature of both PESs and the exoergicity of the quenching products, with the small role played by the centrifugal and RT barriers in the reactant channel, and with the large RT couplings in the OH2+ intermediates. Finally, we contrast the present results with those for the opposite reactions O + H-2 + and for the nonadiabatic-induced quenchings NH + H' and OH + H'.