5- and 6-membered rings: A natural orbital functional study

The Global Natural Orbital Functional (GNOF) provides a straightforward approach to capture most electron correlation effects without needing perturbative corrections or limited active spaces selection. In this work, we evaluate both the original GNOF and its modified variant GNOFm on a set of twelve 5- and 6-membered molecular rings, systems characterized primarily by dynamic correlation. This reference set is vital as it comprises essential substructures of more complex molecules. We report complete-basis-set limit correlation energies for GNOF, GNOFm, and the benchmark CCSD(T) method. Across the Dunning basis sets, both functionals deliver a balanced and accurate description of the molecular set, with GNOFm showing small but systematic improvements while preserving the overall robustness of the original formulation. These results confirm the reliability of the GNOF family and its ability to capture dynamic correlation effects.

💡 Research Summary

This paper presents a comprehensive benchmark of the Global Natural Orbital Functional (GNOF) and its modified variant GNOFm on a set of twelve five‑ and six‑membered aromatic and hetero‑aromatic rings. These molecules—cyclopentadiene, furan, imidazole, pyrrole, thiophene, benzene, pyridine, pyrazine, pyridazine, s‑tetrazine, s‑triazine, among others—are dominated by dynamic electron correlation, making them an ideal testbed for assessing the ability of a density‑matrix‑based method to capture dynamic correlation without resorting to perturbative corrections or active‑space selection.

The authors first review the theoretical foundations of natural orbital functional theory (NOFT), emphasizing the electron‑pairing framework that partitions the orbital space into paired (Ω II) and unpaired (Ω I) subspaces. In this context, GNOF is expressed as a sum of intra‑pair, inter‑pair Hartree‑Fock, static inter‑pair, and dynamic inter‑pair contributions (Eqs. 7‑14). The modification leading to GNOFm re‑introduces interactions between strongly occupied orbitals in antiparallel‑spin blocks, a feature originally proposed in PNOF7, and results in a compact static inter‑pair term (Eq. 16). This adjustment is intended to improve the description of systems where both static and dynamic correlation are present.

All calculations were performed with the open‑source DoNOF code, employing all‑electron Dunning correlation‑consistent basis sets ranging from double‑ζ (cc‑pVDZ) to quintuple‑ζ (cc‑pV5Z) as well as weighted core‑valence variants (cc‑pwCVXZ). For each basis, complete‑basis‑set (CBS) extrapolations were carried out to obtain basis‑set‑limit correlation energies. Benchmark reference data consist of coupled‑cluster singles‑doubles with perturbative triples