I try to simulate FTIR spectra for polyethylene from molecular dynamics (namely LAMMPS) by Fourier transform of dipole moment.I have gained FTIR for nonreactive force fields,but reactive potentials are the target of my study.My previous efforts gave incorrect results.Does anybody has any idea?