Any rigorous approach to first-order reduced density (Γ) matrix functional theory faces the phase dilemma, that is, having to deal with a large number of possible combinations of signs in terms of the electron-electron interaction energy. This problem was discovered by reducing a ground-state energy generated from an approximate N-particle wavefunction into a functional of Γ, known as the top-down method. Here, we show that the phase dilemma also appears in the bottom-up method, in which the functional E [Γ] is generated by progressive inclusion of N-representability conditions on the reconstructed two-particle reduced density matrix. It is shown that an adequate choice of signs is essential to accurately describe model systems with strong non-dynamic (static) electron correlation, specifically, the one-dimensional Hubbard model with periodic boundary conditions and hydrogen rings. For the latter, the Piris natural orbital functional 7 (PNOF7), with phases equal to -1 for the inter-pair energy terms containing the exchange-time-inversion integrals, agrees with exact diagonalization results.

No abstract

Strongly correlated materials are now under intense development, and natural orbital functional (NOF) methods seem to be able to capture the physics of these systems. We present a benchmark based on the Hubbard model for a class of commonly used NOF approximations (also known as reduced density matrix functional approximations). Our findings highlight the importance of imposing ensemble N-representability conditions in order to obtain consistent results in systems with either weak or strong electronic correlation, such as the Hubbard system with a varying two-particle interaction parameter. Based on the accuracy of the results obtained using PNOF7, which retrieves a large amount of the total strong nondynamic correlation, the Hubbard model points out that N-representability gives solid foundations for NOF development.