Extension of Frozen Natural Orbital Approximation to Open-Shell References: Theory, Implementation, and Application to Single-Molecule Magnets

25 November 2019, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


Natural orbitals are often used in quantum chemistry to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual-virtual block of the state density matrix instead of the canonical Hartree-Fock molecular orbitals results in smaller errors when the same fraction of virtual orbitals is frozen. This strategy, termed frozen natural orbitals (FNO) approach, has been successfully used to reduce the cost of state-specific coupled-cluster (CC) calculations, such as ground-state CC, as well as some multi-state methods, i.e., EOM-IP-CC (equation-of-motion CC method for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip), which has been developed to describe certain multi-configurational wave-functions within the single-reference framework. In contrast to EOM-IP-CCSD, which describes open-shell target states by using a closed-shell reference, EOM-SF-CCSD relies on high-spin open-shell references (triplets, quartets, etc). Consequently, straightforward application of FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties, which can be attributed to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. The OSFNO algo-rithm first identifies corresponding orbitals by means of singular value decomposition (SVD) of the overlap matrix of the α and β molecular orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference thus facilitating a safe truncation of the virtual space. The performance of the OSFNO approximation in combination with different choices of reference orbitals is benchmarked for a set of diradicals and triradicals. For a set of di-copper single-molecule magnets, a conservative truncation criterion corresponding to a two-fold reduction of the virtual space in a triple-zeta basis leads to errors of 5–18 cm-1 in the singlet–triplet gaps and errors of ∼2-3 cm-1 in the spin–orbit coupling constants.


Frozen Natural Orbital
Open Shell
coupled cluster theory
Reduced Cost
Single Molecule Magnets

Supplementary materials



Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.