We present an ab-initio approach for grand canonical ensembles in thermal equilibrium with local or nonlocal external potentials based on the one-reduced density matrix. We show that equilibrium properties of a grand canonical ensemble are determined uniquely by the eq-1RDM and establish a variational principle for the grand potential with respect to its one-reduced density matrix. We further prove the existence of a Kohn-Sham system capable of reproducing the one-reduced density matrix of an interacting system at finite temperature. Utilizing this Kohn-Sham system as an unperturbed system, we deduce a many-body approach to iteratively construct approximations to the correlation contribution of the grand potential.

Using the discontinuity of the chemical potential as a function of excess charge, the fundamental gaps for finite systems and the band gaps of extended solids are determined within reduced density matrix functional theory. We also present the necessary and sufficient conditions for the one-body reduced density matrix of a system with fractional charge to be ensemble N-representable. The performance of most modern day reduced density matrix functionals is assessed for the gaps and the correlation energy of finite systems. Our results show that for finite systems the PNOF, BBC3, and power functionals yield very accurate correlation energies while for a correct description of the fundamental gap the removal of self-interaction terms is essential. For extended solids we find that the power functional 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.

In this work, we propose a self-consistent minimization procedure for functionals in reduced density matrix functional theory. We introduce an effective noninteracting system at finite temperature which is capable of reproducing the groundstate one-reduced density matrix of an interacting system at zero temperature. By introducing the concept of a temperature tensor the minimization with respect to the occupation numbers is shown to be greatly improved.