As a final simplification of the model, we will henceforth refer to dµ/dR, as the derivative of µ along the direction of the surface normal, as electrochemical barriers are conventionally computed in quasi-2D surface setups, where we are usually concerned with the dipole moment only in the out-of-plane direction. We obtain the dipole derivative in two ways, similar to Ref 39.