We present a method for the calculation of photoemission spectra in terms of reduced density matrices. We start from the spectral representation of the one-body Green's function G, whose imaginary part is related to photoemission spectra, and we introduce a frequency-dependent effective energy that accounts for all the poles of G. Simple approximations to this effective energy give accurate spectra in model systems in the weak as well as strong correlation regime. In real systems reduced density matrices can be obtained from reduced density-matrix functional theory. Here we use this approach to calculate the photoemission spectrum of bulk NiO: our method yields a qualitatively correct picture both in the antiferromagnetic and paramagnetic phases, contrary to mean-field methods, in which the paramagnet is a metal.

In this work we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard molecule at 1/4 and 1/2 filling as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis are made concerning the approximations of the observables. In particular we focus on the atomic limit, where the two sites of the molecule are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard molecule at 1/2 filling with or without a spin-symmetry-broken ground state, allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW the signature of strong correlation are present in the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover we show how the spectroscopic properties change from one spin structure to the other. Our findings can be generalized to other situations, which allows us to make connections to real materials and experiment.