Here, P and P are TDMs in MO and AO basis, respectively. In our case, the P matrix has a size of 12×12 (since there are 8 occupied and 4 virtual orbitals in the active space), and the P matrix 264×264 (the applied semiempirical method uses 264 AOs to describe the electronic structure of the tetramer).