Restraining rotatable bonds can lead to faster convergence in some cases, however, it can also lead to slower convergence of the free energy estimate in other cases. If e.g. the dihedral angle rotates freely in the solvent, releasing dihedral restraints in the solvent can require the addition of λ windows to obtain sufficient overlap of the work distributions. Therefore, our general SepTop protocol does not restrain rotatable bonds.