Routinely used cutoff strategies for PE calculations introduce unphysical polarization at the outer boundaries of the molecular system. This treatment of the environment polarization may propagate to the quantum subsystem and thus affect its properties. In this work, we have calculated induced dipoles and the associated electrostatic potentials, as well as excitation energies, oscillator strengths, and TPA cross-sections, to scrutinize the effects of different treatments of the induced electrostatics within the PE model for homogeneous solvent environments and heterogeneous biomolecular environments. We assess the effects by applying models A (complete cutoff/droplet), B (partial cutoff), C (full box), and D (full box with MIC), which are presented in Figure . First, we examine the error in the induced dipoles obtained from model C compared to those obtained using model D. Then we examine the environmental electrostatic potentials, i.e., embedding potentials, at the quantum subsystem calculated using models A-D with varying cutoff distances for models A and B. Then we investigate the effects on excitation energies, oscillator strengths, and TPA cross-sections.