Why do you use position restrains during your minimization ? The problem seems to be for a water molecule (because of the atom, number given by the LINCS warning). Try to minimize using SPC flexible: define = -DFLEX_SPC XAvier