In this article we describe a range of simulations (lattice dynamics and molecular dynamics) of the inelastic incoherent neutron scattering spectra of ices (normal ice, ice lI and ice N ). These simulations use a variety of different intermolecular potentials from simple classic pair-wise (rigid and non-rigid molecule) potentials to sophisticated polarisable potentials. It was found that MCY makes stretching and bending interactions too weak while others do them well. We demonstrate that in order to reproduce the measured neutron spectrum, greater anisotropy (or orientational variation) is required than these potentials presently provide. tion and Raman scattering. These are very powerful techniques, which have been highly refined, and their usage has resulted in extensive and valuable data for ice (Li et al., 1992;Li et al., 1991). For ice, however, the normal selection rules governing the interaction of radiation with matter are broken due to the local structural disorder (or proton disordering in crystalline structures of ice), and analysis of the spectral intensities is in general difficult. On the other hand, although IR and Raman spectra are very sensitive to the intramolecular modes involving the O-H stretching and bending, they are less sensitive to the intermolecular modes involving the vibrations of water molecules against each other. Therefore, under normal circumstances, optical spectroscopy provides only limited data in the translational region which is vitally important for us to obtain direct information about the hydrogen bond interaction. In contrast, Inelastic Neutron Scattering (INS) spectroscopy is a more powerful probe, its results are directly proportional to Phonon Density of States (PDOS), INS provides an excellent opportunity for the construction and testing potential functions. Because optical selection rules are not involved, INS measures all modes (IR/Raman measure the modes at the Brillouin Zone q -0) and is particularly suitable for studying disordered systems. It provides direct information on the hydrogen bond interactions in ice.