I am running in single precision, my time-step is 1 fs (0.001 ps, usually works fine). However, I achieved very good energy conservation when I changed integrator from "md" to "md-vv". It sounds quite strange to me, but it works, so... Thanks, however! Francesco