I’m trying to use fix bond/react to create a polymer network. The crosslinker has three reaction sites, all with the same topology, and the monomer has two active sites, both with the same topology. The simulation runs successfully but no reactive events occur despite the atoms that I intended to set with BondingIDs being within the Rmin - Rmax range. Based on dumping the bond_react_MASTER_group, I think the dynamic group is never populated. My guess is that I’m not defining the molecule templates correctly but I haven’t been able to figure out what is wrong. The templates both contain entire molecules and as far as I can tell the topology in the templates matches that of the initial data file. Could someone please advise on what my error might be or where I should look to debug? As a new user I can’t upload attachments, but hopefully I can add them once the topic is created.