If during the deprotonation move the new coordinate falls into the interior of a colloidal particle, the ∆E is counted as infinite, and the move is rejected. The change in electrostatic energy during each move is calculated using Ewald summation with tin foil boundary condition. The Coulomb energy of a periodically replicated charged system is: