A comparison of two classical trajectory surface hopping methods for Na(2P)+H2,D2 quenching

Abstract

Trajectory surface hopping calculations were performed on optimized diatomics-in-molecules surfaces to study Na(2P) collisions with H2 and D2 molecules (v=0, j=1) at four different translational energies (0.039, 0.062, 0.101, and 0.140 eV). Two methods were used to predict surface hopping: (1) transformation of the multidimensional surface intersection to a local one-dimensional curve crossing and calculation of the Landau–Zener transition probability, and (2) integration of the coefficients of the adiabatic electronic states to determine transition probability. For all initial conditions used in this work, we found that method (2) gave significantly larger quenching cross sections. Also in this paper we present results that show nonadiabatic coupling terms calculated by the diatomics-in-molecules method are in good agreement with ab initio values.