“…For a Nth excited state, we proceed from the stationarity condition for the functional where E UHF is the total energy of a system in the UHF formalism; Q α is the orthoprojector on the subspace of virtual orbitals; and λ s and λ o are the Lagrange multipliers, which, in accordance with the AP-method [12,13], ensure, in the limit λ s , λ o → ∞, the fulfillment of conditions for the spin purity and orthogonality of states, respectively (in practice, values λ s = 100 au and λ o = 1000 au ensure the required accuracy). Finally, variations of orbitals in δL (N) = 0, their independence and arbitrariness yield the equations sought for the MOs of α-and β-clusters: (13) Here, f α and f β are the standard Fock operators in the UHF formalism and P β is the orthoprojector on the subspace of occupied orbitals of the β-cluster. Equations (13) , P β = is the β-spin density matrix from occupied orbitals of the α-cluster, Q α = is the density matrix from virtual orbitals of the α-cluster, and P HOMO, k = .…”