code: def check_dissociated(self, cutoff=1.2):
          
          dissociated = False
          if not len(self.B) > self.nslab + 1:  # only one adsorbate
              return dissociated
  
          adsatoms = [atom for atom in self.B[self.nslab:]]
          ads0, ads1 = set(atom.symbol for atom in adsatoms)
          bond_dist = get_ads_dist(self.B, ads0, ads1)
  
          Cradii = [cradii[atom.number]
                    for atom in [ase.Atom(ads0), ase.Atom(ads1)]]
          bond_dist0 = sum(Cradii)
  
          if bond_dist > cutoff * bond_dist0:
              print('DISSOCIATED: {} Ang > 1.2 * {} Ang'
                    .format(bond_dist, bond_dist0))
              dissociated = True
  
              return dissociated
docs: Check if adsorbate dissociates