Now I'm trying to use LIGGGHTS-PFM version to simulate the Bonded Breakage and Particles Replacement Model at the same time. I ran the simulation several times, and the bonded breakage was well-performed, while the PRM had a problem. Every time that the PRM occurred, I got an error;
ERROR on proc 0: Bond atoms 2 1 missing on proc 0 at step 26002 (../neigh_bond.cpp:67)
I suppose that it might caused when atom 2 which is bonded with atom 1 was broken and left atom 1 alone, the bond still remained.
Does anyone have any idea or suggestion about how to fix this error? Thank you in advance for every kind assistance.
The following is my code and visualization of my hypothesis. (Please ignore about the values for each variable which will be fixed later)