I seriously doubt that the choice of integration fix is the cause of your
problems. It almost never is. More likely are errors in other parts of
your input leading to incorrect dynamics. There is far too little
information in you description for any kind of advice.
There is no other dynamics running in my simulation. Only they are bonded
and angled. So, i am attaching my input script as well as read_data file.
Can you please take a look.
Input Script :
#Input script to equilibrate a random network of cross-linked polymer
boundary p p p
atom_modify sort 0 1.