The potential energy surface (PES) for a given vibrational mode at a particular wave vector can be obtained by calculating the internal energies of a set of structures along the corresponding normal-mode coordinate. We have used the program MODEMAP , which computes the displacement of atoms along a given eigenvector e qλj using