We have presented an implementation of the PE model that employs the MIC for the induced electrostatics. Thereby we obtain consistent and accurate polarization throughout the molecular system, as opposed to truncated molecular models, such as the droplet model, that necessarily will result in unphysical polarization at the outer boundary. In this work, we studied the effects of the unphysical polarization for some often used cutoff models: a complete cutoff/droplet model where everything outside the cutoff distance is removed, a partial cutoff model that retains permanent charges/multipoles outside the cutoff distance, and a model that includes everything from a given simulation box but without using the MIC. The latter is essentially a complete cutoff model but without the need to define a cutoff distance (as it is already given by the size of the simulation box). We used a set of homo-geneous solute-solvent systems (PNPs in aqueous solution) and heterogeneous biomolecular systems (NR embedded in a BLG protein or a POPC lipid bilayer) for which we analyzed the induced dipoles in the environment and their associated electrostatic potential at the quantum subsystem, as well as the lowest-lying excitation energy, oscillator strength, and TPA cross-section of the quantum subsystem.