code: def wulff_from_chempot(self, delu_dict=None, delu_default=0, symprec=1e-5,
                             no_clean=False, no_doped=False):
          
  
          latt = SpacegroupAnalyzer(self.ucell_entry.structure). \
              get_conventional_standard_structure().lattice
  
          miller_list = self.all_slab_entries.keys()
          e_surf_list = []
          for hkl in miller_list:
              # For all configurations, calculate surface energy as a
              # function of u. Use the lowest surface energy (corresponds
              # to the most stable slab termination at that particular u)
              gamma = self.get_stable_entry_at_u(hkl, delu_dict=delu_dict,
                                                 delu_default=delu_default,
                                                 no_clean=no_clean,
                                                 no_doped=no_doped)[1]
              e_surf_list.append(gamma)
  
          return WulffShape(latt, miller_list, e_surf_list, symprec=symprec)
docs: Method to get the Wulff shape at a specific chemical potential.
  
          Args:
              delu_dict (Dict): Dictionary of the chemical potentials to be set as
                  constant. Note the key should be a sympy Symbol object of the
                  format: Symbol("delu_el") where el is the name of the element.
              delu_default (float): Default value for all unset chemical potentials
              symprec (float): See WulffShape.
              no_doped (bool): Consider stability of clean slabs only.
              no_clean (bool): Consider stability of doped slabs only.
  
          Returns:
              (WulffShape): The WulffShape at u_ref and u_ads.