However, even if the change in the charge dominates the free energy of changing the protonation state, hydrogen bond re-arrangements can contribute as well. Because the effects due to such structural rearrangements are neglected in LR theories, we hypothesize that higher-order fits may be necessary for obtaining sufficiently accurate correction potentials.