The accuracy of the vibrational self-consistent field (VSCF) method for the computation of anharmonic vibrational frequencies in the infrared (IR) spectrum of formamide and
thioformamide is investigated. The importance of triple potentials
in the commonly used hierarchical expansion of the potential energy surface (PES) is studied in detail. The PES is expanded in terms of Cartesian as well as internal coordinate normal mode displacements. It is found that triples play an important role when using rectilinear coordinates. A VSCF computation based on rectilinear displacements exhibits serious shortcomings which are only remedied by a large vibrational configuration interaction (VCI) treatment including triple potentials. These limitations are partially removed when using curvilinear coordinates. The merits and disadvantages of either type of displacements for the generation of the PES are discussed.