The most fundamental struggle of modern electronic structure theory is the accurate description of electron correlation, which is commonly divided into its dynamic and multireference or strong parts. The former is readily captured with various modern electronic structure methods at polynomial computational cost, while the latter exhibits exponential scaling in traditional wave function based theories. For the precise and tractable computation of molecular electronic properties and processes it is essential to develop methods that capture both dynamic and strong correlation effects at affordable computational cost. Reduced density matrix (RDM) based theories, which exploit the fact that fermions interact pairwise to express the electronic en...