The first contribution is directly computed from the RISM correlation functions. We find a clear domination of this term in all solvation energies reported in Table , with the residual often accounting only for less than 5 % of the total solvation energy. Interestingly, the excess term should account for both stabilizing and destabilizing effect involved in solvation such as the energy cost to replace water molecules. Lee and Schmidt 100 recently pointed out that implicit solvent calculations should be corrected to account for the latter. This cost is accounted for in our bilayer calculations by explicitly removing a water molecule. For the RISM models, we argue that this cost is directly included in the excess term. This is supported by the fact that RISM needs to compute the total number of sites, and therefore the total number of solvent molecules for a given simulation cell volume. Since the excess term is directly computed as a summation over sites, Additionally, RISM models predict a well-structured, symmetric solvation profile on the clean slab which can be seen through the solvent charge distribution plots in Figure .