Bond atoms missing (delete bonds)

Dear all,

I receive the following error when trying to delete bonds and atoms in the
polymer system:

Dear all,

I receive the following error when trying to delete bonds and atoms in
the
polymer system:

----------
read_restart ...

run...

delete_bonds center multi any
delete_atoms region void

timestep 0.0
run 1
-----------
output:
-----------

Deleting bonds ...
51832 total bonds, 48405 turned on, 3427 turned off
91200 total angles, 84875 turned on, 6325 turned off
6384 total dihedrals, 5916 turned on, 468 turned off
36480 total impropers, 33880 turned on, 2600 turned off
Deleted 3078 atoms, new total = 45866
Setting up run ...
ERROR on proc 7: Bond atoms 20855 20854 missing on proc 7 at step
500200 (neigh_bond.cpp:87)

-----------

Does anyone have any idea what is wrong?

please (re-)read the documentation! unless you use the
"remove" option to the delete_bonds command, bonds
are not removed but only disabled. hence, when you
delete the corresponding atoms, you *have* to get the
exact error message that you see.

axel.

Dear Axel,

Thank you for your reply.
I agree, that the input script command I wrote lacks "remove",
nevertheless, before writing to the mailing list I actually tried
it with "remove" as well, the problem is the same:

Dear Axel,

Thank you for your reply.
I agree, that the input script command I wrote lacks "remove",
nevertheless, before writing to the mailing list I actually tried
it with "remove" as well, the problem is the same:

do you have proof, that you didn't delete
any atoms that are still bonded? please
note that delete_bonds will only delete
"wholly contained" entities.

if you believe LAMMPS is at fault, please
produce the smallest possible input that
reproduces the issue and post it to the
list for further debugging.

thanks,
    axel.

Dear Axel,

You comment made me realize my mistake,
the group "center" was defined using region "void",
however in between I run calculations, and thus the
two do not (necessary) contain the same atoms any longer.

The following works for me:

delete_bonds center multi any remove special
delete_atoms group center

Thank you for your help.

Kind regards,
Denis.