To prevent sampling of the non-physical states between λ i = 0 and λ i = 1 on this flat potential energy surface while still enabling sufficient transitions between the physical end-states to sample both protonation states with the correct thermodynamic weight, we introduce the biasing potential V bias i (λ i ) suggested by Donnini et al. (Figure ).