Check if adsorbate dissociates

def check_dissociated(self, cutoff=1.2):
        """Check if adsorbate dissociates"""
        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