I am simulating the polymerization of EDOT monomers using LAMMPS’s fix bond/create command. Polymerization should occur between one Alpha Carbon of one monomer and another alpha Carbon of another monomer. But the problem is that each monomer has two unbonded alpha Carbons. So, the bond creation was done within each monomer rather than between different monomers simply because the two alpha Carbons of each monomer satisfy the bond distance criterion. Can you please suggest to me how to exempt bond creation from two identical atoms (same atom type) within each monomer from bonding?