Following the previous detailed discussion about how to couple domain-decomposition solvation methods to QM host codes, three main modifications were needed inside Psi4: (1) the handling of the additional input parameters for ddX, (2) the computation of solvation energy and Fock matrix contributions as well as (3) the addition of integration routines for electrostatic quantities, which are based on a Becke-type quadrature scheme, which Psi4 uses for DFT computations. This quadrature is used for ddX, since it is able to properly deal with the overlapping atom-centred balls Ω j of the ddX cavity . ( ) and ( ) are strongly dependent on the internals of Psi4, thus of minor interest to implement interfaces with other host codes; we will thus focus our discussion on point (2).