pairs not updated during bond rupture of CG polymer

[caapontes]

Upon bond rupture of a CG polymer model, using the homolysis plugin, the pairs, related to the rupture bond are kept. As a consequence the reducing number of existing bonds (black) do not correlate with the number of pairs (red) (

This is contrary to the binding case, using my specbond plugin, for which pairs and bonds perfectly correlate (

)

The unbinding data can be found here: /hits/fast/mbm/aponteco/BD-shear-flow-polymer/BD/scheme2/BD-crosslinks-kimmdy-unbind/0.001/3/