where q is the proton charge. The sum is performed over the ionic species of valence z i , and (r) = -Z ef f q 4πa 2 δ(r -a) is the uniform charge density at the nanoparticle surface arising from the deprotoneted sites. Application of the Euler-Lagrange condition (7) leads to the general relation