Static calculations were performed using the Vienna Ab Initio Simulation Package (VASP). Valence electrons were described using Projector-Augmented Wave (PAW) pseudopotentials. The exchange-correlation terms of the Hamiltonian were computed using the Generalized Gradient Approximation (GGA) functional of Perdew, Burke, and Ernzerhof. The solvation of the electrode surface was considered using the VASPsol package. The first Brillouin zone was sampled using a 4x4x1 Gamma-centered k-point mesh. The basis-set cutoff was set to 400 eV to ensure proper treatment of O-metal bonds.