In order to find the linear combination coefficients a and b from Eq. (1), a generalized eigenvalue problem in the basis of constructed quasi-diabatic states Φ 1 and Φ 2 needs to be solved. This requires a construction of 2 × 2 overlap S and Hamilton H matrices. The elements of S are equal to determinants of MO overlap matrices with elements