The suggestion was to use "analyze set topo_part_sync" command.
It actually works, but in a strange way. For example it interchanged the "mol" numbers for all the particles and now the virtual particles become the COMs of some artificial molecules, that are sets of unbonded particles from different polymers.
1) if someone has ever dealt with the same problem and found a way to use more than one virtual particles; and