where the indices µ and ν refer to the atomic orbitals (AOs), D AO µν is an element of the density matrix in AO basis, Z m is the nuclear charge of the m'th nucleus out of M nuclei, and the integrals ν AO µν (R i ) are defined by