Although results obtained in the gas phase are interesting, we care more about Aβ-Zn(II) systems in the physiological environment. Looking at our reference systems (see Figure ), we see that the metal centers of these proteins are fully exposed to the solvent, therefore, we have performed calculations in the solvent phase at ε = 78. Results obtained in the solvent phase are shown in Figure and. Once again, we start comparing the relative stabilities of the different clusters (see Figure ). Considering solvent effects, there is an attenuation of the difference in stabilities among the different clusters, but with a shift in the relative stabilities of clusters. Thus, unlike in the gas phase, the q+1 rule is no longer fulfilled, and instead, the most stable clusters have fewer acetate ligands, e.g., either [Zn(Im) 2 (Ac) 2 ] or [Zn(Im) 3 (Ac) 1 ] + .