To account for anharmonic effects in theoretical vibrational spectroscopy, one needs to go beyond the harmonic approximation to the potential energy surface given by Eq. (3) or Eq. ( ). Commonly, this is done by using a truncated n-mode expansion, which we perform in localized-mode coordinates. When including up to two-mode contributions, this corresponds to 87, 88