code: def defect_concentrations(self, chemical_potentials, temperature=300, fermi_level=0.):
          
          concentrations = []
          for dfct in self.all_stable_entries:
              concentrations.append({
                  'conc':
                  dfct.defect_concentration(
                      chemical_potentials=chemical_potentials, temperature=temperature, fermi_level=fermi_level),
                  'name':
                  dfct.name,
                  'charge':
                  dfct.charge
              })
  
          return concentrations
docs: Give list of all concentrations at specified efermi in the DefectPhaseDiagram
          args:
              chemical_potentials = {Element: number} is dictionary of chemical potentials to provide formation energies for
              temperature = temperature to produce concentrations from
              fermi_level: (float) is fermi level relative to valence band maximum
                  Default efermi = 0 = VBM energy
          returns:
              list of dictionaries of defect concentrations