[gmx-users] Constraint algorithm error

[nicolas Dinter]

Dear Gromacs users,

I am simulting a 64 diphenyl carbonate molecule box, at 1bar and 298K. After
building a big box (15*15*15) I am reducing the box size by NPT. So far it
runs well but now the calculation crash after several steps.  In the md.log
file I found the error: "constraint error in algorithm step 84910...".
I wonder if it's not a Gromacs problem, because in my mdp file I specified:
constraints= none, and I tried to use the shake algorithm with a 0.0001
shake_toll and I also tried the lincs algorithm with lincs_order = 4.
Any idea?
Thanks.

Nicolas Dinter