Pourbaix App Manual

From Wiki
Jump to: navigation, search


A Pourbaix diagram, also frequently called a potenial-pH diagram, or E-pH diagram, is a representation of aqueous phase electrochemical equilibria. It is a two-dimensional representation of a three-dimensional free energy-pH-potential diagram. In other words, it depicts the stable phase as a function of pH and potential. Potential is defined with respect to the standard hydrogen electrode.

Experimentally determining Pourbaix Diagrams is painstaking work, as we need not only the free energy of aqueous ions, but also of all solid phases that a system can exist in. Since Materials Project offers a very convenient, yet powerful database of materials, it is quite straightforward to utilize this data to generate Pourbaix diagrams in a high-throughput manner.

This manual outlines the usage of the Pourbaix App to calculate Pourbaix diagrams, and the thermodynamic formalism underlying the app.

Using the Computational Pourbaix App

The Pourbaix app is capable of plotting elemental, and multi-element Pourbaix diagrams. To construct an Elemental Pourbaix diagram, enter the element of choice, and click on the Generate button.

To generate a multi-elemental Pourbaix diagram, choose the desired multiple elements from the periodic table, and click the Generate button. Note that oxygen and hydrogen are included by default, since these elements are always "open" in a Pourbaix diagram.

For multi-elemental diagrams, sliders are provided to specify a composition of the elements. Note that for each composition, a new Pourbaix diagram is generated from scratch.

Interfacing with the data

  • To ensure a clutter-free diagram, the domains on the Pourbaix diagram are not labeled by default. To view labels on the plot, select the "Label Domains" check box.
  • Each domain has a mouse-over point located at the center of the domain. Mousing over these points displays the entries corresponding to the domain.
  • Domains which remain solid in solution are shaded. This helps identify passivation regions, especially in multi-elemental systems, were identifying passivation regions is tricky.
  • Zooming into domains is quite simple! Just drag a selection window using your mouse. To reset zoom, click on the reset zoom button which appears on the upper right corner.
  • Unlike the Phase Diagram App, data tables are shown by default. Mousing over rows in the Stable Compounds table highlights the corresponding entries in the Pourbaix diagram. Some stable entries may not get highlighted. This is because, the domains corresponding to these entries lies outside the standard limits of the Pourbaix diagram.
  • The unstable compounds table lists the unstable entries, and their corresponding energies above hull.
  • Links in the data tables lead to more information about the corresponding entries. Unlinked entries are ions. Mousing over the book icon next to the ions shows the reference for the free energy of formation of the aqueous ion.

Thermodynamic formalism of Pourbaix Diagrams

Figure 1: Schematic to reference experimental aqueous ion energies to DFT data

To calculate a Pourbaix diagram, free energies of the solid phases, and of the aqueous ions is required. Calculating free energies of ions is tricky, and time-consuming. To overcome this problem, a methodology utilizing experimentally measured free energies of aqueous ions and the calculated DFT energies for solid phases available in the Materials Project was developed. [1] Note that the correction scheme described below is applied over and above any compatibilities/corrections which are applied to the species

Referencing energies of aqueous ions

Briefly, for each ion, a reference solid is chosen, and the correction term is calculated for the ion as the energy difference between the experimental and the DFT calculated energies of the reference solid. The basic idea behind this scheme is that, if we have a reference energy for an aqueous ion which reproduces the correct dissolution for one solid, then accurate DFT solid-solid energy differences ensure that all other solids dissolve accurately with respect to that ion. The better the solid is represented by DFT, the more transferable the reference aqueous energy becomes. We therefore prefer to choose simple chemical systems (primarily binaries with an uncomplicated electronic structure) as representative solids.

For an aqueous ion i at standard state conditions (e.g., room temperature, atmospheric pressure, and 10^{-6} M concentration) using a representative solid s, we define the chemical potential as:

\begin{align} \mu_{i(aq)}^0 &= \mu_{i(aq)}^{0,{\rm exp}} + \left[\Delta g_s^{0,{\rm DFT}} - \Delta g_s^{0,{\rm exp}}\right] \\
& = \mu_{i(aq)}^{0,{\rm exp}} + \Delta\mu_s^{0, {\rm DFT - exp}}

Figure 1 shows this schematically.

Correction for Water

In an aqueous environment, many chemical and electrochemical reactions are enabled by the breakdown, formation, or incorporation of water molecules. It is therefore important that the free energy of formation of water is captured accurately. This is known accurately from experiments as -2.46 eV. So, at standard state, the free energy of formation of water is set as follows:

\mu^0_{{\rm H_2O}} &= \Delta g^{0, {\rm exp}}_{\rm H_2O} \\
 & = -2.46 eV

Correction for elemental gaseous, and liquid states, except H_2

For all gaseous elements, the experimentally determined entropic contribution at 298 K is added to the DFT/corrected energy of the element as follows:

\mu_i^{\rm ref} = E_i^{0, {\rm DFT}} + \Delta E_i^{\rm correction} - Ts_i^{\rm exp}

This is implemented for the following elements: O_2, F_2, Cl_2, Br_2, Hg

Correction for H_2

In an aqueous environment, O_2 and H_2 in their gaseous states are in equilibrium with water through the reaction

{\rm H}_{2(g)} + \frac{1}{2}{\rm O}_{2(g)} \leftrightarrow {\rm H}_2{\rm O}_{(l)}

Hence, the hydrogen energy is corrected such that the experimental free energy of formation of H_2O is reproduced.

\mu_{\rm H}^{\rm ref} & = \frac{1}{2}\left[g_{\rm H_2O}^0 - \mu_O^0 - \Delta g_{\rm H_2O}^{0, {\rm exp}} \right] \\
& \approx \frac{1}{2}\left[ E_{\rm H_2O}^{0, {\rm DFT}} -Ts_{\rm H_2O}^0 - \mu_O^0 - \mu_{\rm H_2O}^{0, {\rm exp}} \right]

Demonstration of the Pourbaix App

Figure 2: Pourbaix diagram for elemental Fe.
Figure 3: Orange lines are the H2 and O2 lines in the Pourbaix diagram.
Figure 4: Multi-elemental Fe-Cr Pourbaix diagram at a composition of 18% Cr, 82% Fe.

This section briefly demonstrates both elemental, and multi-elemental Pourbaix diagrams.

Elemental Pourbaix Diagram

Figure 2 shows the elemental Pourbaix diagram for Fe. Note that in the current version of the Pourbaix App, the concentration of all ions is fixed at 10^{-6}M, and does not offer a way to change this default setting.

The two orange lines are the hydrogen reduction line, and the line denoting water oxidation to O2. These are clearly labeled in Figure 3. These lines provide a good reference to predict whether hydrogen/oxygen gas is evolved at a particular set of conditions. For example, water is unstable below the H2 line, and so, hydrogen gas evolves at the cathode at conditions below this line. Similarly, above the O2 line, oxygen gas is evolved at the anode.

Multi-elemental Pourbaix diagram

The stability of multiple-elements in aqueous environments is predicted using multi-elemental Pourbaix diagrams like the one shown in Figure 4. The composition slider bar can be seen above the plot. Here, the small white bar separating the two colors can be clicked on and dragged to change the ratio of Fe to Cr. This may, or may not have any effect on the Pourbaix diagram. More information about how multi-elemental Pourbaix diagrams vary as a function of composition can be found elsewhere. [2] [3]

For an n-element diagram, there are n coexisting phases. Note that these can be any mixture of solid, and aqueous phases. So, for the two-element diagram shown in Figure 4, there are two coexisting phases in each domain. As mentioned above, domains shaded blue indicate purely solid domains. For multi-elemental diagrams, shaded domains indicate those in which purely solid phases coexist.

Literature references for ions

The free energies of ions in the aqueous phase have been taken from standard references/recent publications. The acronyms which show up on the tool-tips associated with the aqueous ions, and their corresponding references are as follows.

  • NBS Tables: NBS Thermodynamic tables.[4]
  • M. Pourbaix (1974): Atlas of Electrochemical Equilibria in Aqueous Solutions. [5]
  • Barin Knacke Kubaschewski: Thermochemical Properties of Inorganic Substances [6]
  • Barner and Scheuerman (1978): Handbook of thermochemical data for compounds and aqueous species [7]
  • Beverskog and Puigdomenech (1997): Beverskog and Puigdomenech, Corr. Sci. (1997) [8]


  1. K. A. Persson, B. Waldwick, P. Lazic, and G. Ceder, Phys. Rev. B, 85, 235438 (2012)
  2. S. Jayaraman, S. P. Ong, K. A. Persson, and G. Ceder, in preparation (2013)
  3. Pourbaix Diagrams for Multielement Systems, Thompson, W. T., Kaye, M. H., Bale, C. W. and Pelton, A. D. (2011), in Uhlig's Corrosion Handbook, Third Edition (ed R. W. Revie), John Wiley & Sons, Inc., Hoboken, NJ, USA.
  4. NBS Technical Note 270-1 to 270-8. D. D. Wagman et. al, U. S. Department of Commerce (1973)
  5. Atlas of Electrochemical Equilibria in Aqueous Solutions, M. Pourbaix,, NACE (1974)
  6. Thermochemical Properties of Inorganic Substances, I. Barin, O. Knacke, and O. Kubaschewski, Springer-Verlag, Berlin (1977)
  7. Handbook of thermochemical data for compounds and aqueous species, H. E. Barner, and R. V. Scheuerman, Wiley, New York, 1978
  8. Revised Pourbaix diagrams for Ni at 25-300^oC, B. Beverskog and I. Puigdomenech, Corr. Sci., 39, 969-980 (1997)


  1. Sai Jayaraman