The first solvation shell of a water molecule comprises 4-5 water molecules (Fig. shows a typical structure from the equilibrium simulations). Thus, our minimal QM system for computing bulk IEs comprises 5 water molecules. As we show below, the convergence of the core IEs with respect to the QM size is slow and much larger QM systems are needed for converged results.