I am trying to construct Pourbaix Diagram for some 2-dimensional arbitrary structures. I have checked most references, but I could not understand what inputs are required for constructing Pourbaix diagrams. My main catalytic reaction is Oxygen Evolution, so I guess I have to just calculate the total energies of my structures bonded to OH and O (as intermediates), along with the total energy of my pristine structure. Are these values enough for constructing the Pourbaix diagram by Pymatgen? Are the potential energies of ions such as O, OH, and those related to my structure included in the Pymatgen database?; how can I check which data are already available in Pymatgen and which data are missing?
In addition, I could not find any handy tutorial on using Pymatgen to construct Pourbaix diagrams.