The important question is how to estimate the order of magnitude of possible artifacts introduced by the damping and understand the implications when applying KS-DFT/AMOEBA with PWs model to realistic complex biomolecular systems. As a starting point, we choose here to look into molecular properties and SCF convergence in solute-solvent systems, where a QM solute is surrounded by many MM solvent molecules, reproducing the bulk solution properties.