deltaphi phiM-phi "Electrode-OHP potential difference" rho_surf epsilon0_const*eps_S*deltaphi/xS "Surface charge density"