New ab initio coupled potential energy surfaces for the Br(2P3/2, 2P1/2) + H2 reaction

Abstract

The three lowest (1A('), 2A('), and 1A('')) adiabatic potential energy surfaces (PESs) for the Br((2)P) + H(2) reactive system have been computed based on the multi-reference configuration interaction (MRCI) method including the Davidson's correction with a large basis set. These three adiabatic PESs have been transformed to a diabatic representation, leading to four coupling potentials. In addition, the spin-orbit matrix elements were also obtained using the Breit-Pauli Hamiltonian and the unperturbed MRCI wavefunctions in the Br + H(2) channel and the transition state region. Consequently, six coupling potentials were obtained and their characteristics were extensively discussed. Nonadiabatic quantum dynamics calculations for this system have been realized with these realistic diabatic potentials instead of previous semi-empirical diabatic potentials. Based on two-state model nonadiabatic calculations for the Br((2)P(3∕2), (2)P(1∕2)) + H(2) reaction, the Br((2)P(1∕2)) + H(2) reaction was found to show less reactivity than the Br((2)P(3∕2)) + H(2) reaction at collision energies beyond the threshold of the Br((2)P(3∕2)) + H(2) reaction. Our results are consistent with the previous studies on the XH(2) (X = F, Cl) system, which indicate that the adiabatically forbidden channel is dominant at low energies in the open-shell halogen atom plus H(2) reactions.