Hello.
I need to simulate the phonon dispersion curves of bccFe. Can lammps do this?

I don’t know how this calculation is generally done. LAMMPS does
not do it directly  my guess is that this
is better as a postprocessing task. Maybe someone
else can comment.
Steve
Sialon, as far as I know lammps does not do this directly, but it can be done by postprocessing a dump file. See Lukes & Zhong, J. Heat Transfer 129(6) 2007 for the necessary theory.
Jeremy
Jeremy and Sialon,
As far as I know, what have done in Lukes & Zhong’s paper is plotting phonon density of state not dispersion curve via fast Fourier transformation of velocity autocorrelation (common method). For CNT, to obtain the dispersion curve, one can basically get the phonon dispersion curve by manupilating the velocity of atoms. Then one may need some knowledge on lattice dynamics to get the dispersion curve by connecting these discrete points to be lines in dispersion curve.
I’m not sure how to do with bulk materials like bccFe. What people usually do with ab initio is first get the force on each atoms in the unit cell, and diagonalize the lattice dynamic matrix to get the eigenvalues. This requires some sort of postprocessing. I’m sure the paper publishing the potential should have the info of dispersion curves therein.
Best,
Ajing