We then present our results including free energies of solvation and adsorption with both RISM water models and CANDLE. At the same time, we include comparison with the static-bilayer approximation, a widely used method where only the first layer of explicit water is considered. We further discuss these results by plotting the 2D averaged solvent charge at the interface in an attempt to explain the disparities between our models. In a second part we explain how we compute workfunctions for solvated systems, and compare them between the models. We then perform calculations outside the PZC at different surface charge densities. For each of these calculations we plot the RISM density distributions for hydrogen and oxygen water sites and show that the model correctly predicts the H-up to H-down flip around the PZC. We continue to discuss the implications of our findings and the need for a comparison with properly sampled Ab-Initio Molecular Dynamics (AIMD) simulations.