where k runs on the occupied MOs in the ground state configuration and c kµ are the expansion coefficients of the k-th MO on the AOs (see ESI, Eqs. S1 and S2). The HF Hamiltonian is solved self-consistently. Typically, 40-50 HF iterations are needed to reach convergence on the density matrix elements.