Here, we represent the N λ i -coordinates as an N -dimensional vector λ. The 2 N possible protonation states of the system are represented by the N -dimensional vector l with elements l i equal to 0 or 1 that indicate whether a site i is protonated (λ i = 0) or deprotonated (λ i = 1).