This embedding potential includes the full electrostatic nuclear and electronic Coulomb potentials of the monomers in the environment as well as non-classical contributions due to the exchange-correlation and kinetic energy, which are evaluated using approximate functionals. For obtaining the electron densities of the fragments in the environment, the simplest variant is the use of those obtained in calculations for the isolated monomers. Since for ion-water clusters large polarization effects can be expected, we also iteratively updated all embedded monomer densities using freeze-andthaw (ft) cycles. While the calculation of the relaxed monomer densities comes with some computational effort, this only needs to be done at the level of the monomers, i.e., the scaling of the computational effort remains unchanged.