Model A is a complete cutoff model, also known as the droplet model. In this approach, a smaller-often spherical-part of the simulation box is cut out based on a cutoff distance from the quantum subsystem. The environment outside the cutoff distance is thus completely neglected. Model B is a partial cutoff model, where the same cutoff distance is used, but only for the induced electrostatics, i.e., polarizabilities outside the cutoff distance are removed while permanent multipoles (or force-field charges) are included for all sites in the entire simulation box. In model C, the full simulation box is used, i.e., the permanent multipoles and polarizabilities for all sites are included. Note that even though this model does not depend on a cutoff distance measured from the quantum part (as in models A and B), there is a complete cutoff at the edges of the simulation box, i.e., there are no periodic boundaries, and it is, in that respect, similar to model A. Finally, in model D (i.e., the PE-MIC approach), the permanent multipoles and polarizabilities for all sites in the simulation box are included as in model C, but the MIC is used for the induced electrostatics as described in the Theory section.