To investigate this further, we consider the SAPT0 + MBD energy components. These are plotted for each of the ions in Fig. , again separating bulk and interfacial environments and ensemble-averaging over either data set. In considering the energy decomposition in Eq. ( ), we have opted to group first-order electrostatics and exchange together,