aiqbench analysis
This research advances Reduced Density Matrix Functional Theory (RDMFT) by incorporating 4-index cumulant contributions to accurately model bond breaking and strongly correlated electronic states. By enforcing physical constraints on local spin and delocalization indices, the new functional achieves exact agreement with CASSCF energies for dissociated molecules like N2 and O2. This methodological improvement addresses a critical limitation in current quantum chemistry approaches where fragment properties are often inaccurately reproduced during chemical reactions.
Towards more accurate natural orbital functional approximations: including 4-index cumulant contributions
arXiv:2607.21399v1 Announce Type: cross Abstract: 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.
Abstract: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.
Chemical Physics (physics.chem-ph); Quantum Physics (quant-ph)
arXiv:2607.21399 [physics.chem-ph]
(or arXiv:2607.21399v1 [physics.chem-ph] for this version)
https://doi.org/10.48550/arXiv.2607.21399 Focus to learn more arXiv-issued DOI via DataCite (pending registration)
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.
Full text under the publisher's open licence · 317 words
Related coverage
Running Variational Quantum Eigensolver with Qiskit Aer on AMD Instinct
AMD · May 29, 2026
News Archive 2026
Quantum Computing Report · Jul 25, 2026

News
Quantum Computing Report · Jul 25, 2026
J++ Lens: Jacobian Filtering Enables More Faithful Workspace Lenses
LessWrong · Oct 8, 2026