Hi Mathieu, Thanks for replying. I used the fourier code in CPMD and dipoles calculated with cp2K (format was modified) to calculated IR spectrum. And I find that the peaks are too noisy and the positions some main peaks are not accurate. (Simulated time >10ps and >600 data points in IR spectrum). Any advise? Regards, Ross, SUN