A molecular force field dedicated to molecular dynamics simulation of biomembranes was developed. It was parameterized on model compounds related to phospholipids and was able to reproduce at the same time structures, energies, and vibrational spectra. Cross terms in the potential energy function were introduced by solving the redundancy problem among internal coordinates. This force field was used in the 400-ps molecular dynamics simulation of a hydrated bilayer in the gel and liquid crystal phases. The conformational properties of the polar head groups were in particular agreement with the experimental observations using Raman scattering.