For solutes 1 (Figure ) and 2 (Figure ), hydrogen bound water molecules were reintroduced into the QM region. For the first snapshot for each solute, this means that there are, respectively, 4 and 8 water molecules included in the QM region. For the negatively charged solute 2 the increased number of water molecules is due to more favorable interactions with the surrounding water molecules. Again, the remaining water molecules are treated with either TIP3P, WAT or M2P2 embedding potentials. Compared to the results obtained for the small QM region (Figure ) there are definite improvements for solute 1 also when using TIP3P as an embedding potential. Specifically, progress is observed around the central peak at 5.9 eV where all three embedding potentials are now reproducing this part of the spectrum as opposed to only the polarizable models in the previous section. This again confirms that the peak position and intensity are dominated by mutual polarization of the solutes and the solvents. Considering that even the peak around 6.7 eV is observed to improve for all three embedding potentials, it is suggested that the effect has character of energy transfer instead of being influenced by polarization effects. It is noted that the position of the peak around 5.0 eV is still poorly reproduced by the TIP3P potential. For solute 2, similar conclusions can be made where an overall improvement across the entire spectrum is observed. This is valid for both the peak position (3.3 eV) in the low-energy regime and also for the higher energy excitations. However, at the peak around 5.3 eV a shift of 0.1 eV is observed for the polarizable models. This discrepancy is most likely due to over-polarization effects which could probably be ameliorated by including some form of damping when solving for the induced moments. When computing the deviations of the spectra to the reference spectrum for solute 1 one obtains 0.38, 0.27 and 0.24 for W1, W2 and W3, respectively. This is, for W1 and W2 almost an improvement by a factor of two. However, because the spectrum for model W3 with the small QM region was very good (δ dev = 0.16), an increase in deviation is observed for the large QM region to 0.24. Similar results are obtained for solute 2 where the improvements for all three water models are almost a factor of three for TIP3P (W1) down to a more modest improvement for the M2P2 potential (W3) of a factor 1.4.