In Eq. 29, D jk is the excited state transition density matrix between the excited states j and k. If the ϕ n (r) are expressed in an atomic orbital basis {χ α (r)} according to ϕ n (r) = α d n α χ α (r), Eq. ( ) can be rewritten as