code: def remove_bond(self, particle_pair):
          
          from mbuild.port import Port
          if self.root.bond_graph is None or not self.root.bond_graph.has_edge(
                  *particle_pair):
              warn("Bond between {} and {} doesn't exist!".format(*particle_pair))
              return
          self.root.bond_graph.remove_edge(*particle_pair)
          bond_vector = particle_pair[0].pos - particle_pair[1].pos
          if np.allclose(bond_vector, np.zeros(3)):
              warn("Particles {} and {} overlap! Ports will not be added."
                   "".format(*particle_pair))
              return
          distance = np.linalg.norm(bond_vector)
          particle_pair[0].parent.add(Port(anchor=particle_pair[0],
                                           orientation=-bond_vector,
                                           separation=distance / 2), 'port[$]')
          particle_pair[1].parent.add(Port(anchor=particle_pair[1],
                                           orientation=bond_vector,
                                           separation=distance / 2), 'port[$]')
docs: Deletes a bond between a pair of Particles
  
          Parameters
          ----------
          particle_pair : indexable object, length=2, dtype=mb.Compound
              The pair of Particles to remove the bond between