Hi Teo, I checked the numerical differentiation step by doing a force field calculation for the diamond structure, where any QM related noise is absent. The step size was 0.01 and I got zeros of the order of 0.1 cm-1. Therefore I concluded it is not the step size which is causing the absence of nice zeros. An