23 d-block metal surfaces were modeled with GC-DFT using JDFTx. The CANDLE implicit solvent model was used to account for the effects of the aqueous electrolyte at the electrochemical interface. All bulk crystal structures were taken from the Materials Project. All bulk lattice constants were first optimized in JDFTx, from which surface facets were obtained using Pymatgen. Following the methods of Montoya, 8 the (110) surface facets were cut from body-centered cubic bulk structures and the (111) and (211) surfaces were cut from face-centered cubic structures. The top and bottom faces of each slab are separated by 10 Å of space to provide sufficient distance for charge screening by the solvation model. The flat (110) and (111) surfaces have 333 atoms per unit cell. The stepped (211) surfaces have 334 atoms per unit cell. The bottom layer of surface atoms was frozen to the bulk geometry. Surfaces (with and without adsorbates) are converged to 0.1 meV (0.0001 eV) by relaxing the ionic positions. All of the converged structures in this work can be found at , which also includes the standardized calculation parameter files. All surfaces and their corresponding reaction energy diagrams are provided in the SI.