code: def has_antibnd_states_below_efermi(self, spin=None, limit=0.01):
          
          warnings.warn("This method has not been tested on many examples. Check the parameter limit, pls!")
  
          populations = self.cohp
          number_energies_below_efermi = len([x for x in self.energies if x <= self.efermi])
  
          if populations is None:
              return None
          elif spin is None:
              dict_to_return = {}
              for spin, cohpvalues in populations.items():
                  if (max(cohpvalues[0:number_energies_below_efermi])) > limit:
                      dict_to_return[spin] = True
                  else:
                      dict_to_return[spin] = False
          else:
              dict_to_return = {}
              if isinstance(spin, int):
                  spin = Spin(spin)
              elif isinstance(spin, str):
                  s = {"up": 1, "down": -1}[spin.lower()]
                  spin = Spin(s)
              if (max(populations[spin][0:number_energies_below_efermi])) > limit:
                  dict_to_return[spin] = True
              else:
                  dict_to_return[spin] = False
  
          return dict_to_return
docs: Returns dict indicating if there are antibonding states below the Fermi level depending on the spin
              spin: Spin
              limit: -COHP smaller -limit will be considered.