POURBAIX_DIAGRAM Module

This command or module is not available with GES6.

GES6 is a re-write and new implementation of the Gibbs Energy System module of Thermo‑Calc and it is the default engine. The previous version, GES5, is still available and will continue to work in the background of the program where necessary. The main purpose of GES6 is to support faster development of new features. However, not all GES5 functionality is implemented in GES6. In such cases (except for the POURBAIX module), GES6 falls back silently and uses functionality from GES5 in the background.

Getting started with the Pourbaix module

After you enter the POURBAIX_DIAGRAM module, you are asked if you want to have step-by-step instuctions to help walk you through the module.

Need any instruction on the Pourbaix module? /N/: Y

If you answer Y a series of notes are displayed. The following is a more comprehensive summary of these same notes.

What is a Pourbaix diagram?

The Pourbaix diagram is a calculated equilibrium phase diagram mapped and plotted with the independently-varied electropotential (Eh, as defined with regard to the standard hydrogen electrode as its reference) and acidity (pH), that represents all the equilibrated phase boundaries among aqueous solution, gaseous mixture, and various primary and secondary solids (modeled as either complex solution or simple stoichiometric phases) in a certain multicomponent heterogeneous interaction system, under a defined T-P condition and a specific initial bulk composition (which is, by default, always set as 1 kg of water dissolving a specified amount of metals/alloys and other acids/alkalines/salts).

What database should I use?

At least one database is required that contains an aqueous solution phase (with thermodynamic data for water and various aqueous solutes), that shall be selected from TCAQ (PAQ) or AQS (PAQS) [which use the SIT Model or the Complete Revised HKF Model, respectively] within the Thermo‑Calc databases.

Choosing a custom aqueous solution database

You can also use a custom database but it must be in the Thermo‑Calc TDB format. Due to the restrictions of aqueous solution models used within Thermo‑Calc, such a database must be designed in the same format as in the default TCAQ (PAQ2) for aqueous solution phases. Among other things, keep the following requirements in mind:

The ELECTRON is defined as an special element (ZE) and as the only constituent in its reference phase REFERENCE_ELECTRODE (for determining the electro-potential that is defined as Eh with the standard hydrogen electrode as the reference), but it is not defined as an aqueous species;

The vacancy (VA) is unnecessary for AQUEOUS solution phase and it should be avoided in the definition of phase-constituents in the AQUEOUS phase;

The AQUEOUS solution phase should always be defined as a constitutional solution phase, implying that all the aqueous solution species must be included in a single site, rather than in two or multiple sublattices.

Beside the AQUEOUS solution phase, there shall exist a GAS mixture phase containing at least H2O1, O2 and H2; and for multicomponent systems, normally there shall also contain some solid (stoichiometric or solution) phases. Of course, if desired, you could also choose to calculate and generate a Pourbaix diagram without considering the GAS mixture phase entirely; however, such a plot is not really a complete Pourbaix diagram, due to the fact that thermodynamically-stable phase boundaries between the AQUEOUS solution phase and GAS mixture phase will then not be calculated at all!

All the required thermodynamic data for calculations of Pourbaix diagrams or other diagrams must be retrieved either from one (Single) database which consists of an AQUEOUS solution phase, a GASeous mixture phase, a REF_ELECTRODE phase, and some SOLID phases (being solutions and/or stoichiometric compounds; for primary metals/alloys and for secondary products formed from heterogeneous chemical/electrochemical interactions, or from several (Multiple) databases that respectively contain various solutions/compounds (as listed above).

Such databases suitable for calculations of aqueous-bearing heterogeneous interaction systems can be those default-prompted ones [i.e., in the Single-Database case, the PAQ or PAQS; and in the Multiple-Database case, the PAQ2 or as primarily-switched database, plus the SSUBx as firstly-appended one and the SSOLx as secondly-appended one if it is necessary; even more databases can be appended). Of course, you could also choose to append required data from other appropriate databases (such as TCFE, TCSLD, TCNI/TTNi, TCAL/TTAl, TCMG/TTMg, TTZr, TCMP, SLAG, etc.) for GASeous mixture phase and for varous solid solution and stoichiometric compound phases. Furthermore, an experienced user can also utilize his/her own USER-specified databases in various cases.

The advanced POURBAIX module is designed so that it requires you to just answer some necessary questions, rather than to go through the other Console Mode modules (i.e., TDB, GES5, POLY3 and POST) although you still enter commands at the prompts. Beside the default plotted Pourbaix diagram, it also allows the user to easily and quickly plot many different properties of the system, stable phases and aqueous species, varied along the calculated phase boundaries for the same defined heterogeneous interaction system. Moreover, it permits the user to directly change some plotting settings and manipulate all kinds of plotted diagrams.

The current advanced POURBAIX-Module has been extended so that it is additionally able to directly perform a normal STEPPING calculation (varied with a specified independent variable) and to easily generate various types of property diagrams, for the same heterogeneous interaction system that has been defined in a previous POURBAIX or TDB-GES5-POLY3-POST calculation. 

Copy and Rename the POLY3 file

At the end of a Pourbaix module calculation, a new POLY3 file is automatically saved.

A reminder that Linux and Mac are case sensitive; when applicable ensure you enter file extensions with capital letters.

If you are using a POLY3 file more than once (as is done in the extended examples), it is recommended that you copy and rename the previously-generated POLY3 file. This is so that the required POLY3 file structure and Pourbaix results are not lost. Do this outside of Thermo‑Calc i.e. do not save from within Thermo‑Calc.

The reason you do not save from within Thermo‑Calc is because of how the file is saved. When the command sequence (GO POLY3>SAVE <file name>.POLY3>Y) is used, the POLY3 workspace of the POURBAIX calculation results in the POLY3 module monitor and all the initial equilibrium points and the mapping/stepping calculation results are lost as well as parts of the POLY3 file structure (e.g., some definitions of the previously defined POURBAIX-type calculation system).

Enforce a Pause

After the initial set of information is displayed, you are asked:

Enforce a PAUSE after plotting when running a MACRO? /N/

  • If you choose Y a pause is always enforced after each diagram is plotted;
  • If you choose N, all the command-lines (answers) for the session continue without stopping or until the SET_INTERACTIVE_MODE line in the macro file is reached.

Whenever running a Thermo‑Calc macro (TCM) file, you may want to pause after a specific diagram is plotted on screen. This gives you the flexibility to:

  • Print it (of EMF/PS format) to a printer
  • Convert it (of EMF format) to (PDF) graphical files;
  • Save it as an EMF graphical file;
  • Dump it as a PNG or BMP graphical file;
  • Set background colour for the current diagram
  • You can also set the default font/size and change plot layers for sequential plots.

Pourbaix Main Options

If you answer Y to the Pause prompt above, these options are displayed:

  1. Start a completely new Pourbaix diagram calculation.
  2. Open an old file and plot other property diagrams.
  3. Open an old file and make another Pourbaix calculation.
  4. Open an old file and make another stepping calculation.

POURBAIX_DIAGRAM Commands for Options 1 to 4