Returns the plot of the formation energy of a particle against its
            effect radius

        Args:
            analyzer (SurfaceEnergyPlotter): Analyzer associated with the
                first polymorph
            max_r (float): The maximum radius of the particle to plot up to.
            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.
            label (str): Label of the plot for legend
            increments (int): Number of plot points
            delu_default (float): Default value for all unset chemical potentials
            plt (pylab): Plot
            from_sphere_area (bool): There are two ways to calculate the bulk
                formation energy. Either by treating the volume and thus surface
                area of the particle as a perfect sphere, or as a Wulff shape.
            r_units (str): Can be nanometers or Angstrom
            e_units (str): Can be keV or eV
            normalize (str): Whether or not to normalize energy by volume

def plot_one_stability_map(self, analyzer, max_r, delu_dict=None, label="",
                               increments=50, delu_default=0, plt=None,
                               from_sphere_area=False, e_units="keV",
                               r_units="nanometers", normalize=False,
                               scale_per_atom=False):
        """
        Returns the plot of the formation energy of a particle against its
            effect radius

        Args:
            analyzer (SurfaceEnergyPlotter): Analyzer associated with the
                first polymorph
            max_r (float): The maximum radius of the particle to plot up to.
            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.
            label (str): Label of the plot for legend
            increments (int): Number of plot points
            delu_default (float): Default value for all unset chemical potentials
            plt (pylab): Plot
            from_sphere_area (bool): There are two ways to calculate the bulk
                formation energy. Either by treating the volume and thus surface
                area of the particle as a perfect sphere, or as a Wulff shape.
            r_units (str): Can be nanometers or Angstrom
            e_units (str): Can be keV or eV
            normalize (str): Whether or not to normalize energy by volume
        """

        plt = plt if plt else pretty_plot(width=8, height=7)

        wulffshape = analyzer.wulff_from_chempot(delu_dict=delu_dict,
                                                 delu_default=delu_default,
                                                 symprec=self.symprec)

        gform_list, r_list = [], []
        for r in np.linspace(1e-6, max_r, increments):
            gform, r = self.wulff_gform_and_r(wulffshape, analyzer.ucell_entry,
                                              r, from_sphere_area=from_sphere_area,
                                              r_units=r_units, e_units=e_units,
                                              normalize=normalize,
                                              scale_per_atom=scale_per_atom)
            gform_list.append(gform)
            r_list.append(r)

        ru = "nm" if r_units == "nanometers" else "\AA"
        plt.xlabel(r"Particle radius ($%s$)" %(ru))
        eu = "$%s/%s^3$" %(e_units, ru)
        plt.ylabel(r"$G_{form}$ (%s)" %(eu))

        plt.plot(r_list, gform_list, label=label)

        return plt