Trying to address the details and the respective strengths and weaknesses of the different classical polarization models and with that the steps required in their parametrization according to points (i) and (ii) exceeds the scope of this review and we point the reader to the original literature, instead. We will focus on general remarks to points (iii) and (iv), which are to a large extent agnostic to the specifics of the MM models used. Duchemin et al. proposed to not calculate the real-space form of the reaction field, v I-I reac (r, r ′ , ω), as in Eq. ( ) but instead exploit the resolution-of-identity (RI) method. Recall that in Eq. ( ), auxiliary basis functions ξ µ were used to calculate two-center Coulomb repulsion integrals