Metal/water interfaces play an essential role in a wide variety of (electro-)chemical processes. Computational modeling of such interfaces is essential for developing a detailed, atomistic understanding of the structure, dynamics, and reactivity that occur at the interface during solutionphase heterogeneous catalysis and/or electrocatalysis. As such, there have been many attempts to develop methodologies to model these interfaces. These range from those based upon an implicit treatment of the surrounding water, to fully atomistic ab initio molecular dynamics (AIMD) or hybrid quantum mechanics/molecule mechanics (QM/MM) methods that utilize an explicit treatment of the surrounding electrolyte at QM and MM levels, respectively.