What is the problem with high viscosity?


I have a lammps input file with viscosity. When I have viscosity = 1000, I have no problem but when I increase the viscosity to 10000 I get to the following error:

ERROR on proc 0: Bond atoms 1 2 missing on proc 0 at step 10 (…/ntopo_bond_all.cpp:63)
I do not know why this happens!
It would be great if somebody could give me a hint.


Can be caused by a variety of things. Bad initial geometry, bad initial velocities, bad force field, bad simulation parameters, bad choice of boundary conditions, bad neighborlist settings, and probably some more issues, that I cannot recall at the moment.