Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0910.2635 · Physical Review A

Physical meaning of the natural orbitals: Analysis of exactly solvable models

N. Helbig, I. V. Tokatly, Ángel Rubio

Published Feb 9, 2010• 19 citations
arXivHTMLPDF

Abstract

We investigate the suitability of natural orbitals as a basis for describing many-body excitations. We analyze to which extent the natural orbitals describe both bound as well as ionized excited states and show that depending on the specifics of the excited state the ground-state natural orbitals may yield a good approximation. We show that the success of reduced density-matrix functional theory in describing molecular dissociation lies in the flexibility provided by fractional occupation numbers while the role of the natural orbitals is minor.

DOI
10.1103/physreva.81.022504

Related papers

Orbitals from local RDMFT: Are they Kohn-Sham or natural orbitals?

Iris Theophilou, Nektarios N. Lathiotakis, Nikitas I. Gidopoulos +2 more

Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.

arXiv:1505.00627•The Journal of Chemical Physics•Aug 5, 2015•12 citations

Orbitals from local RDMFT: Are they Kohn-Sham or natural orbitals?

Theophilou, Iris, Lathiotakis, Nektarios N., Gidopoulos, Nikitas I. +2 more

Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.

•Jan 1, 2015•5 citations

Natural orbitals, overlap functions, and mean-field orbitals in an exactly solvableA-body system

Dimitri Van Neck, A.E.L. Dieperink, Michel Waroquier

We consider a simple but nontrivial many-body system interacting through short-range forces, and confirm a property of such systems that was recently proposed on general grounds, namely, that single-particle overlap functions, spectroscopic factors, and separation energies of the bound (A-1)-particle eigenstates can be derived from the one-body density matrix of the A-particle system in its ground state. The basis of natural orbitals for the system is constructed and its properties are discussed. We also investigate the high-momentum content of the bound-state overlap functions and momentum distribution. It is found that the mean field provides a good approximation for the bound-state overlap functions even in the region of large momenta, where the total momentum distribution is already enhanced by several orders of magnitude over the mean-field result. \textcopyright{} 1996 The American Physical Society.

Physical Review C•May 1, 1996•14 citations

A favorably-scaling natural-orbital functional theory based on higher-order occupation probabilities

Ralph Gebauer, Morrel H. Cohen, Roberto Car

We introduce a novel energy functional for ground-state electronic-structure calculations. Its fundamental variables are the natural spin-orbitals of the implied singlet many-body wave function and their joint occupation probabilities. The functional derives from a sequence of controlled approximations to the two-particle density matrix. Algebraic scaling of computational cost with electron number is obtainable in general, and Hartree-Fock scaling in the seniority-zero version of the theory. Results obtained with the latter version for saturated small molecular systems are compared with those of highly-accurate quantum-chemical computations. The numerical results are variational, capturing most of the correlation energy from equilibrium to dissociation. Their accuracy is considerably greater than that obtainable with current density-functional theory approximations and with current functionals of the one-particle density matrix only.

arXiv:1309.3929•Sep 16, 2013•0 citations

A well-scaling natural orbital theory

Ralph Gebauer, Morrel H. Cohen, Roberto Car

We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix.

arXiv:1611.08692•Proceedings of the National Academy of Sciences•Nov 1, 2016•17 citations

The adiabatic approximation in time-dependent density matrix functional theory: Response properties from dynamics of phase-including natural orbitals

Klaas J. H. Giesbertz, O. V. Gritsenko, Evert Jan Baerends

The adiabatic approximation is problematic in time-dependent density matrix functional theory. With pure density matrix functionals (invariant under phase change of the natural orbitals) it leads to lack of response in the occupation numbers, hence wrong frequency dependent responses, in particular α(ω→0)≠α(0) (the static polarizability). We propose to relinquish the requirement that the functional must be a pure one-body reduced density matrix (1RDM) functional, and to introduce additional variables which can be interpreted as phases of the one-particle states of the independent particle reference system formed with the natural orbitals, thus obtaining so-called phase-including natural orbital (PINO) functionals. We also stress the importance of the correct choice of the complex conjugation in the two-electron integrals in the commonly used functionals (they should not be of exchange type). We demonstrate with the Löwdin-Shull energy expression for two-electron systems, which is an example of a PINO functional, that for two-electron systems exact responses (polarizabilities, excitation energies) are obtained, while writing this energy expression in the usual way as a 1RDM functional yields erroneous responses.

The Journal of Chemical Physics•Nov 4, 2010•32 citations
We investigate the suitability of natural orbitals as a basis for describing many-body excitations. We analyze to which extent the natural orbitals describe both bound as well as ionized excited states and show that depending on the specifics of the excited state the ground-state natural orbitals may yield a good approximation. We show that the success of reduced density-matrix functional theory in describing molecular dissociation lies in the flexibility provided by fractional occupation numbers while the role of the natural orbitals is minor.
Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.
Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.
We consider a simple but nontrivial many-body system interacting through short-range forces, and confirm a property of such systems that was recently proposed on general grounds, namely, that single-particle overlap functions, spectroscopic factors, and separation energies of the bound (A-1)-particle eigenstates can be derived from the one-body density matrix of the A-particle system in its ground state. The basis of natural orbitals for the system is constructed and its properties are discussed. We also investigate the high-momentum content of the bound-state overlap functions and momentum distribution. It is found that the mean field provides a good approximation for the bound-state overlap functions even in the region of large momenta, where the total momentum distribution is already enhanced by several orders of magnitude over the mean-field result. \textcopyright{} 1996 The American Physical Society.
We introduce a novel energy functional for ground-state electronic-structure calculations. Its fundamental variables are the natural spin-orbitals of the implied singlet many-body wave function and their joint occupation probabilities. The functional derives from a sequence of controlled approximations to the two-particle density matrix. Algebraic scaling of computational cost with electron number is obtainable in general, and Hartree-Fock scaling in the seniority-zero version of the theory. Results obtained with the latter version for saturated small molecular systems are compared with those of highly-accurate quantum-chemical computations. The numerical results are variational, capturing most of the correlation energy from equilibrium to dissociation. Their accuracy is considerably greater than that obtainable with current density-functional theory approximations and with current functionals of the one-particle density matrix only.
We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix.
The adiabatic approximation is problematic in time-dependent density matrix functional theory. With pure density matrix functionals (invariant under phase change of the natural orbitals) it leads to lack of response in the occupation numbers, hence wrong frequency dependent responses, in particular α(ω→0)≠α(0) (the static polarizability). We propose to relinquish the requirement that the functional must be a pure one-body reduced density matrix (1RDM) functional, and to introduce additional variables which can be interpreted as phases of the one-particle states of the independent particle reference system formed with the natural orbitals, thus obtaining so-called phase-including natural orbital (PINO) functionals. We also stress the importance of the correct choice of the complex conjugation in the two-electron integrals in the commonly used functionals (they should not be of exchange type). We demonstrate with the Löwdin-Shull energy expression for two-electron systems, which is an example of a PINO functional, that for two-electron systems exact responses (polarizabilities, excitation energies) are obtained, while writing this energy expression in the usual way as a 1RDM functional yields erroneous responses.