This dataset includes the full computational workflows of a density functional theory based thermodynamics model for the overpotential of the oxygen evolution reaction (OER) on an oxidized palladium surface. The model assumes an oxygen bridge vacancy as an active site on the 110 surface of the tetragonal PdO2 (rutile type structure). The critical OER potential has been computed with variation of the Fe modifyer type, either dopant or adsorbate, and modifyer position. Furthermore, an alternative bifunctional pathway of OER has been considered by adding an H atom to an auxiliary O-bridge site from which the proton−electron pair for second OER reaction step is released rather than from the hydroxylated active site. The computed OER overpotential on the Fe-free surface via this bifunctional route is 0.42 V. A substitution of Pd with Fe directly at this active site further reduces the calculated OER overpotential, over the same route, to 0.35 V. This 70 mV decrease in overpotential is in good agreement with the experimentally measured decrease of 60 mV in the onset potential.