In an attempt to tackle the various challenges previously cited, an important amount of work has been dedicated to the development of implicit solvation models. While these models have proved helpful for small, isolated, solutes, their accuracy is often assumed transferable to the electrochemical interface without extensive justifications. Additional terms such as dispersion and repulsion corrections or entropic considerations are often heavily parametrized for specific molecular datasets. For this reason, their usefulness for different environments such as electrochemical interfaces is unknown. Additional issues arise when considering electrochemical systems such as the lack of quadrupole interactions which can considerably affect energy levels. Conventional Poisson-Boltzmann models depict electrolytes as continuous charge-density with exponential decay, therefore ignoring angular and orientation contributions. This leads to a simplified model of the double layer with important shortcomings when considering effects arising far from the Potential of Zero Charge (PZC) or when the size of electrolytes becomes significant.