Note that both functions ξ µ and ξ ν are auxiliary functions in the QM region (subsystem I). These integrals can be computed by considering the basis function ξ µ as a charge distribution that polarizes the embedding region according to the adopted MM model. Then, the action of this response field on the basis function ξ ν is calculated.