The OO-CCD orbitals are very similar to the (non-optimized) Brueckner orbitals, and it has been shown previously that the corresponding electron density and associated molecular properties are generally more accurate than those obtained from the Hartree-Fock determinant. Implementation Previously, the db-MBE was implemented in our scripting framework PyAdf 60 in combination with DFT calculations using the Amsterdam Modeling Suite (Ams). For details, we refer to the description in Ref. 42. In this case, we made use of Ams's own tools for exporting electron densities and Coulomb potentials on an integration grid, which had previously been implemented to allow for the export of embedding potentials. While PyAdf is interfaced to numerous quantum-chemical program packages, this interface was specific to Ams. Therefore, we decided to implement a more general, high-level interface that abstracts the handling of the electron densities and Coulomb potentials on a grid from the calculation of these densities. This general interface then allows for an evaluation of the density-based correction [Eq. ( )] that is agnostic to the source of the electron densities and Coulomb potentials.