In the previous implementation of constant pH MD in GROMACS' fork, the smooth interpolation of the force field potential energy function between the protonated and deprotonated states was achieved by linearly interpolating the force field potentials of these states. Thus, for a single titratable site, the λ-dependent potential is given by