The successful wavefunction optimization using the KS-DFT/AMOEBA model with PWs for the model dimer systems points out that the modified interaction tensor with atomic covalent radii is effectively screening the electrostatic potential, thus enabling us to perform calculations even when the MM region is overlapping with the real-space grid. Still, we can not completely avoid unphysical artefacts that can be introduced with the ad hoc damping parameters, since the repulsion is not treated explicitly. However, we suggest our model as a cheap and straightforward way to control the screening of electrostatic potential of QM/MM interactions and thus mitigate electron spill-out for short range electrostatic and polarization interactions.