The effective potential energy surfaces of the nonadiabatic collision BPja2 +H2S +g1,n,j ↔BPj'a 2+H2S +g1,n' ,j'

Effective potential energy surfaces (PESs) are calculated for the nonadiabatic collision BPja2 +H2S +g1,n,j ↔BPj'a 2+H2S +g1,n' ,j' . This calculation employed 1 2A', 2 2A', and 1 2A" adiabatic PESs numerically calculated at the state-averaged multiconfigurational self-consistent field (SA-MCSCF)/configuration interaction (CI) level for several values of the H2 bond length, H2 orientation angle, and boron distance. The associated nonadiabatic coupling terms (NACTs) were calculated from the SA-MCSCF/CI wave functions using analytic gradient techniques. A line integral through the NACTs was then used to determine the adiabatic-to-diabatic mixing angle required to transform from the 1 2A' and 2 2A' adiabatic basis to a corresponding diabatic basis. When all nonadiabatic coupling terms between all electronic states are considered, the line integral is path independent. However, only NACTs between the 1 2A' and 2 2A' states were considered in these calculations, and the line integral was therefore path dependent. The path dependence of the line integral was used to characterize the error introduced by employing a truncated set of adiabatic states. A method for reducing the effect of this error through the use of symmetry derived boundary conditions was developed. The resulting diabatic PESs were combined with the total B + H2 rotational kinetic energy and boron spin-orbit coupling to yield diabatic effective PESs. The diabatic effective PESs were diagonalized to yield adiabatic effective PESs. Diabatic effective PESs data was extracted for the equilibrium H2 bond length and used to calculate inelastic scattering matrix elements using the time dependent channel packet method. These matrix elements were compared previous results [D. E. Weeks et al., J. Chem. Phys. 125, 164301 (2006)] to observe the sensitivity of this calculation to the input electronic structure data.