Workfunctions play a crucial role in electrochemistry due to their direct connection with electrode potentials. Computing workfunctions in vacuum is relatively straightforward, as it is calculated by taking the difference between the Fermi energy of the electrode and the electrostatic potential in vacuum close to the surface. However, determining solvated workfunctions requires additional considerations. In this section, we first explain the methodology used to compute the PZC for each model and then compare the obtained workfunctions with those reported in the literature. The absolute value of the Standard Hydrogen Electrode (SHE) can be directly computed by considering the equilibrium of reaction 7.