Performance of one-body reduced density matrix functionals for the homogeneous electron gas

The subject of this study is the exchange-correlation-energy functional of reduced density matrix functional theory. Approximations of this functional are tested by applying them to the homogeneous electron gas. We find that two approximations recently proposed by Gritsenko, Pernal, and Baerends, J. Chem. Phys., {\bf 122}, 204102 (2005), yield considerably better correlation energies and momentum distributions than previously known functionals. We introduce modifications to these functionals which, by construction, reproduce the exact correlation energy of the homogeneous electron gas

Benchmark calculations for reduced density-matrix functional theory

Reduced density-matrix functional theory (RDMFT) is a promising alternative approach to the problem of electron correlation. Like standard density functional theory, it contains an unknown exchange-correlation functional, for which several approximations have been proposed in the last years. In this article, we benchmark some of these functionals in an extended set of molecules with respect to total and atomization energies. Our results show that the most recent RDMFT functionals give very satisfactory results compared to more involved quantum chemistry and density functional approaches.Comment: 17 pages, 1 figur

Discontinuity of the chemical potential in reduced-density-matrix-functional theory

We present a novel method for calculating the fundamental gap. To this end, reduced-density-matrix-functional theory is generalized to fractional particle number. For each fixed particle number, $M$, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to a function, $E_{\mathrm{tot}}^M$, whose derivative with respect to the particle number has a discontinuity identical to the gap. In contrast to density functional theory, the energy minimum is generally not a stationary point of the total-energy functional. Numerical results, presented for alkali atoms, the LiH molecule, the periodic one-dimensional LiH chain, and solid Ne, are in excellent agreement with CI calculations and/or experimental data.Comment: 9 pages, 3 figures, version as publishe

Local reduced-density-matrix-functional theory: Incorporating static correlation effects in Kohn-Sham equations

We propose a novel scheme to bring reduced density matrix functional theory (RDMFT) into the realm of density functional theory (DFT) that preserves the accurate density functional description at equilibrium, while incorporating accurately static and left-right correlation effects in molecules and keeping the good computational performance of DFT-based schemes. The key ingredient is to relax the requirement that the local potential is the functional derivative of the energy with respect to the density. Instead, we propose to restrict the search for the approximate natural orbitals within a domain where these orbitals are eigenfunctions of a single-particle hamiltonian with a local effective potential. In this way, fractional natural occupation numbers are accommodated into Kohn-Sham equations allowing for the description of molecular dissociation without breaking spin symmetry. Additionally, our scheme provides a natural way to connect an energy eigenvalue spectrum to the approximate natural orbitals and this spectrum is found to represent accurately the ionization potentials of atoms and small molecules

Open shells in reduced-density-matrix-functional theory

Reduced-density-matrix-functional theory is applied to open-shell systems. We introduce a spin-restricted formulation by appropriately expressing approximate correlation-energy functionals in terms of spin-dependent occupation numbers and spin-independent natural orbitals. We demonstrate that the additional constraint of total-spin conservation is indispensable for the proper treatment of open-shell systems. The formalism is applied to the first-row open-shell atoms. The obtained ground-state energies are in very good agreement with the exact values as well as other state of the art quantum chemistry calculationsComment: 4 pages, 2 figures, corrected typo

Reduced Density Matrix Functional for Many-Electron Systems

Reduced density matrix functional theory for the case of solids is presented and a new exchange correlation functional based on a fractional power of the density matrix is introduced. We show that compared to other functionals, this produces more accurate results for both finite systems. Moreover, it captures the correct band gap behavior for conventional semiconductors as well as strongly correlated Mott insulators, where a gap is obtained in absence of any magnetic ordering.Comment: 4 figs and 1 tabl

Conditions for describing triplet states in reduced density matrix functional theory

We consider necessary conditions for the one-body-reduced density matrix (1RDM) to correspond to a triplet wave-function of a two electron system. The conditions concern the occupation numbers and are different for the high spin projections, $S_z=\pm 1$, and the $S_z=0$ projection. Hence, they can be used to test if an approximate 1RDM functional yields the same energies for both projections. We employ these conditions in reduced density matrix functional theory calculations for the triplet excitations of two-electron systems. In addition, we propose that these conditions can be used in the calculation of triplet states of systems with more than two electrons by restricting the active space. We assess this procedure in calculations for a few atomic and molecular systems. We show that the quality of the optimal 1RDMs improves by applying the conditions in all the cases we studied

Generalized Pauli constraints in reduced density matrix functional theory

Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman's ensemble $N$-representability conditions. Recently, the topic of pure-state $N$-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble $N$-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone