Assigning signs to the electronic nonadiabatic coupling terms: the H2,O system as a case study

This paper is devoted to a specific difficulty related to the electronic nonadiabatic coupling terms (NACT), namely, how to determine correctly their signs. It is well known that correct NACTs, including their signs, are crucial for any numerical treatment of the nuclear Schrodinger equation [see, i.e., A. Kuppermaan and R. Abrol, Adv. Chem. Phys. 124, 283 (2003)]. In most cases the derivation of the correct sign of the nonadiabatic coupling matrix (NACM) is done employing various continuity procedures. However, there are cases where these procedures do not suffice and for these cases we suggest to apply an additional procedure based on a mathematical lemma which asserts that the exponentiated line integral which yields the D matrix is invariant with respect to the initial point of the integration [M. Baer, J. Phys. Chem. A 104, 3181 (2000)]. In the numerical study we apply this lemma to determine the signs of the 3x3 NACM elements for the three excited states of the {H(2),O} system (some of these NACTs are presented here for the first time). It turns out that the ab initio treatment yields results from which one can form eight different 3x3 NACMs. However the application of this lemma (which does not require any significant additional numerical effort) reduces this number to two. The final selection is done by an enhanced numerical study which requires more accurate calculations.