nMOLDYN has the capability to use the velocity auto correlation function to obtain the phonon density of state for an NVT simulation. There are problems with normalization of the spectra. Does anybody know a better package for this purpose?
You get the velocities of each atom from MD, then use Matlab to read the velocities and calculate the VACF, make the fourrier transform and plot it as a function of frequency.