diff --git a/pymatgen/analysis/phase_diagram.py b/pymatgen/analysis/phase_diagram.py
index 68902be86e8..d3ebfc909da 100644
--- a/pymatgen/analysis/phase_diagram.py
+++ b/pymatgen/analysis/phase_diagram.py
@@ -1608,7 +1608,7 @@ def _get_2d_plot(self, label_stable=True, label_unstable=True,
             plt.ylim((miny - ybuffer, ybuffer))
             center = (0.5, miny / 2)
             plt.xlabel("Fraction", fontsize=28, fontweight='bold')
-            plt.ylabel("Formation energy (eV/fu)", fontsize=28,
+            plt.ylabel("Formation energy (eV/atom)", fontsize=28,
                        fontweight='bold')
 
         for coords in sorted(labels.keys(), key=lambda x: -x[1]):