Search papers, labs, and topics across Lattice.
This paper addresses the challenge of accurately modeling bond breaking in reduced density matrix functional theory (RDMFT) by refining the natural orbital functional through the inclusion of 4-index cumulant contributions. The authors correct the cumulant contribution from the PNOF5 functional to ensure that it aligns with the known physical properties of dissociated fragments, such as local spin and delocalization index. Benchmark results demonstrate that the updated functional can reproduce complete active space self-consistent field (CASSCF) energies in the dissociation regime for various molecular systems, significantly enhancing performance in strongly correlated scenarios.
Correcting cumulant contributions in natural orbital functionals leads to exact energy predictions for bond dissociation, bridging a critical gap in RDMFT accuracy.
Accurate modeling of bond breaking remains a central challenge for reduced density matrix functional theory (RDMFT). Although some modern functionals can yield reasonably accurate dissociation energies, they often fail to reproduce key properties of the dissociated fragments, such as a vanishing fragment population covariance (also known as the delocalization index) and the correct total spin angular momentum of each fragment (local spin). In this work, we revisit the construction of natural orbital functionals by correcting the cumulant contribution produced by the PNOF5 functional. Our method enforces known contributions of the cumulant to local spin fragments and the delocalization index at the dissociation limit. We obtain the closest cumulant consistent with these physically motivated constraints and subsequently purify the corresponding one- and two-electron reduced density matrices by imposing the standard $P, Q, \text{and } G$ $N$-representability conditions. The resulting functional yields improved behavior in strongly correlated regimes. Benchmarking on the dissociation of the singlet states of \ce{N2}, \ce{NO+}, \ce{O2}, \ce{S2}, and \ce{CO} shows that in the dissociation regime the energies computed from the updated cumulant exactly reproduce the complete active space self-consistent field (CASSCF) energies. We further analyze the limitations of the approach and identify scenarios in which the current approach performs poorly. This work provides a pathway for systematically improving natural orbital functionals to achieve reliable bond-breaking calculations within RDMFT.
