We also introduced the first comprehensive convex hull for the cubic Prussian Blue system, considering all 24 available 24d sites. From this convex hull, we determined an upper bound of 25 meV for the mixing energy difference between the most stable sodium configuration and the next metastable arrangement at equivalent sodium concentrations. A significant clustering of configurations was observed across all sodium concentrations in terms of mixing energy, suggesting considerable degeneracy in sodium configurations and weak