In conclusion, the continuum treatment of the GMPNP equations provides a computationally tractable method to model the reaction environment of CO 2 ER close to highly charged surfaces, although with certain limitations. We believe that our approach can provide a qualitatively sound starting point for further theoretical as well as experimental investigations while bridging some of the gap between quantum mechanical surface reactivity studies and mass transport of species from the bulk solution in electrocatalytic systems.