The workflow begins with the selection of the two input structures: one representing the cubic modification of Prussian Blue and the other representing the rhombohedral modification. Following the identification of all sodium sites within these structures, the pymatgen library is used to generate permutations of possible sodium orderings for each concentration.