At a surface charge density σ = 0, the calculated potential U M|S (SHE) corresponds to the PZC of the system vs. SHE. For the RISM (TIP5P) model at 1 M, we obtain a value of 0.39 V vs. SHE, which is in good agreement with previous implicit solvation work: 0.54 V vs. RHE 27 and 0.57 V vs. SHE. 117 Moreover, our result is consistent with experimental values of 0.35 V vs. RHE 118 and 0.37 V vs. SHE. In contrast, the SPC model yields an underestimated value of 0.10 V vs. SHE, most likely due to the poor solvation energies for H 3 O + . The CANDLE model gives a higher value of 0.58 V, in good agreement with previous studies using implicit solvation models.