Phase Diagram App¶
Welcome to the Phase Diagram App Reference!
Here you will find all relevant information regarding the Phase Diagram App (PDApp), including release notes, the manual, etc.
Release Notes¶
Recent changes¶
Version 0.3  May 2, 2011¶
 Bug fixes for tabbing behavior, and unstable phase listing.
 New scroll layout for Data Tables for easier viewing.
Version 0.2  March 28, 2011¶
 New tabbed formatting for PDApp which supports up to 7 PDs.
Version 0.1  March 10, 2011¶
 Initial release of Phase Diagram App
 Key features:
 24 compositional components phase diagrams supported.
 Highlighting of nodes on mouseover.
 Marking of an arbitrary composition in the phase diagram.
 Mouseover popup of details at particular nodes.
 Complete data table with listing of stable and unstable phases.
Manual¶
Introduction¶
Phase diagrams represent the thermodynamic phase equilibria of multicomponent systems and reveal useful insights into fundamental material aspects regarding the processing and reactions of materials. However, the experimental determination of a phase diagram is an extremely timeconsuming process, requiring careful synthesis and characterization of all phases in a chemical system.
Computational tools can accelerate phase diagram construction significantly. By calculating the energies of all known compounds in a given chemical system (e.g. \(\ce{LiFeO}\)), we can determine the phase diagram for that system at 0 K and 0 atm. Furthermore, for systems comprised of predominantly solid phases open with respect to a gaseous element, approximations can be made as to the finite temperature and pressure phase diagrams.
In this manual, we will outline the general usage and thermodynamics methodology of our computational phase diagrams app. The details of the calculations, including computational methodology and accuracy, used in constructing these phase diagrams are available in the Calculations Manual and will not be further covered in this manual.
Using the Computational Phase Diagram App¶
To construct the phase diagram of a 24 component chemical system (e.g. \(\ce{LiFeO}\)), you only need to select the desired elements from the periodic table and click the "Generate" button.
To generate a grand potential phase diagram, for which one of the components of the phase diagram is considered 'open' to the system, select a projected element from the dropdown list, and input a desired chemical potential either manually or using the temperature slider (for supported elements \(\ce{O2}\), \(\ce{N2}\), \(\ce{H2}\), \(\ce{F2}\), and \(\ce{Cl2}\) only). When available, the temperature slider has been calibrated using entropy data from NIST. ^{1}
Once the phase diagram has been generated, there are multiple ways to interact with the data:
 When mousing over nodes in the phase diagram, a popup appears that identifies the phase. The popup lists both stable and unstable crystal structures computed at that composition. You may click the id link to bring up more details of that compound, such as detailed crystal structure information.
 Checking the 'Show Unstable' box will display the compositions at which unstable compounds have been computed as blue circles.
 You can zoom into a portion of the phase diagram by defining a zoom area: click and hold at the upperleft corner of the desired area, and drag down and right, releasing the mouse at the lowerright corner of the desired area. A popup with the zoomedin portion should appear. To return to the original phase diagram, click anywhere outside the popup.
Interpreting Phase Diagrams¶
This section provides a short guide on interpreting phase diagrams generated by the PDApp. A more comprehensive discussion and example of phase diagram interpretation can be found in References 2 and 3.^{2}^{3}
Basic Phase Diagram Information¶
While the PDApp generates several types of phase diagrams, some information is common to all phase diagrams. As an illustrative example, we refer to Figure 1, a ternary \(\ce{LiFeO}\) phase diagram at 0 K and 0 atm calculated using the PDApp.
Figure 1: Calculated LiFeO Phase Diagram
Several pieces of information can be obtained from Figure 1:
 The red nodes on the phase diagram represent phases that are calculated to be stable under the given conditions.
 Blue circles (not shown) represent phases that are calculated to be unstable under the given conditions.
 An arbitrary composition will thermodynamically decompose to stable phases (red nodes), with the relative proportion of each stable phase determined by the lever rule.
Regarding polymorphic compounds (compounds that change crystal structure with respect to temperature or pressure), only limited information can be obtained because phase diagrams generated by the PDApp do not include an explicit temperature axis. For most solid materials, our experience indicates that temperature and pressure contribute fairly small changes to differences in total energy. Thus, we expect that crystal structures for which the PDApp reports similar formation energies (i.e., within about 50 meV/atom) might be candidates for polymorphism.
The PDApp models only the thermodynamics of a system. Thus, compounds that are metastable (kinetically stabilized) may not show up as stable on the generated phase diagram. We expect that compounds for which the decomposition energy into stable phases is low are potential candidates for metastability; however, compounds with extremely high decomposition energies into stable phases (i.e., over about 200 meV/atom) are less likely to be metastable. The thermodynamic scale of inorganic crystalline metastability is discussed in Reference 4. ^{4}
Compositional Phase Diagrams¶
In a compositional phase diagram, the system is closed with respect to the environment. We display compositional phase diagrams involving different numbers of components in slightly different ways.
Binary Compositional Phase Diagrams¶
Figure 2 is an example of a calculated binary \(\ce{FeP}\) phase diagram at 0 K and 0 atm using the PDApp. Binary phase diagrams show the complete convex hull for the system, where the yaxis is the formation energy per atom and the xaxis is the composition (e.g. in the \(\ce{FeP}\) phase diagram, the xaxis is the fraction of \(\ce{P}\)).
Figure 2: Calculated FeP Phase Diagram
The black lines show the convex hull construction, which connects stable phases. Unstable phases will always appear above the convex hull line; one measure of the thermodynamic stability of an arbitrary compound is its distance from the convex hull line, which predicts the decomposition energy of that phase into the most stable phases.
Ternary Compositional Phase Diagrams¶
Figure 1 is an example of a calculated ternary \(\ce{LiFeO}\) phase diagram at 0 K and 0 atm using the PDApp. The presentation of ternary phase diagrams differs from that of binary phase diagrams; the energy axis is removed so that the entire compositional space can be represented. (The energy axis, not shown, can be thought of as coming 'out of the page').
The black lines in the ternary phase diagrams are projections of the convex hull construction into compositional space. The lines form Gibbs triangles, which can be used to find stable phases at an arbitrary composition. At any point in the phase diagram other than the stable nodes, the equilibrium phases are given by vertices of the triangle bounding that composition. For example, the equilibrium phases for a composition with \(\ce{Li:Fe:O}\) ratio of 1:1:1 (i.e., in the center of the phase diagram), is predicted to be \(\ce{Li5FeO4}\), \(\ce{LiFeO2}\) and \(\ce{Fe}\).
Quaternary Compositional Phase Diagrams¶
Quaternary phase diagrams are presented in the same general way as ternary phase diagrams, but an additional axis is needed to represent the fourth composition. Therefore, quaternary phase diagrams are shown in three dimensions rather than as a 2D plot.
The lines in quaternary phase diagrams define polyhedra with 4 vertices each rather than triangles. At any point in the phase diagram other than the stable nodes, the equilibrium phases are given by the vertices of the polyhedron bounding that composition.
Grand Potential Phase Diagrams¶
Grand potential phase diagrams are phase diagrams representing phase equilibria that are open to a particular element. In these phase diagrams, the chemical potential of an element is an external input variable that the user must specify. In some cases (e.g. a gaseous element such as oxygen), the chemical potential can be related to temperature and partial pressure via thermodynamic relations and known experimental data.^{6}
A grand potential phase diagram involving n components is presented in the same way as a compositional phase diagram involving \(n1\) components.
Figure 3: Calculated LiFeO grand potential phase diagram at oxygen chemical potential consistent with approximately 700 K and 0 atm.
Figure 3 shows the calculated grand potential phase diagram for the \(\ce{LiFeO}\) system at an oxygen chemical potential consistent with 700 K and 0 atm. The phase diagram provides information on the equilibrium phases under these conditions for any \(\ce{Li:Fe}\) composition (note that the O fraction is controlled by interactions with the environment). For example, we may observe that a compositon with a relatively high Fe ratio (near composition = 1) will comprise \(\ce{Fe2O3}\) and \(\ce{LiFeO2}\). In addition, by examining the changes in the equilibrium phases as the chemical potential changes, we can obtain insights into how the phase equilibria are modified with changes in the chemical potential. A detailed discussion of this application can be found in References 2 and 3. ^{2}^{3}
Accuracy of Calculated Phase Diagrams¶
Figure 4: Experimental LiFeO Phase Diagram
Figure 4 shows the experimentally determined phase diagram for the \(\ce{LiFeO}\) system at 673K.^{5} Even though the experimental phase diagram is at a much higher temperature, we can see that our calculated phase diagram reproduces the features very well. The phases \(\ce{Li2O}\), \(\ce{Fe3O4}\), \(\ce{Fe2O3}\), \(\ce{LiFeO2}\), \(\ce{Li5FeO4}\) are stable in both the experimental diagram and our calculated phase diagram. The only additional phase in our calculated diagram is \(\ce{FeO}\), which is well known to be difficult to obtain in stoichiometric proportion under normal conditions. Our calculations correctly reproduce the \(\ce{FeO}\) formation enthalpy to within experimental accuracy, and our calculated phase diagram is fullly consistent with known experimental thermochemical data along the \(\ce{FeO}\) line at 1 atm and 298K.
In general, we can expect that compositional phase diagrams comprising of predominantly solid phases to be reproduced fairly well by our calculations. However, it should be noted that there are inherent limitations in accuracy in the DFT calculated energies. Furthermore, our calculated phase diagrams are at 0 K and 0 atm, and differences with nonzero temperature phase diagrams are to be expected.
For grand potential phase diagrams, further approximations are made as to the entropic contributions.^{2}^{3} They are therefore expected to be less accurate, but nonetheless provide useful insights on general trends.
Thermodynamics Methodology¶
In the remainder of this section, we will use an \(\ce{LiFeO}\) system to illustrate the methodology without loss of generality.
Compositional Phase Diagrams¶
To construct a phase diagram, one would need to compare the relative thermodynamic stability of phases belonging to the system using an appropriate free energy model. For an isothermal, isobaric, closed system, the relevant thermodynamic potential is the Gibbs free energy, \(G\), which can be expressed as a Legendre transform of the enthalpy, \(H\), and internal energy, \(E\), as follows:
where \(T\) is the temperature of the system, \(S\) is the entropy of the system, \(P\) is the pressure of the system, \(V\) is the volume of the system, and \(N_i\) is the number of atoms of species \(i\) in the system.
For systems comprising primarily of condensed phases, the \(PV\) term can be neglected and at 0K, the expression for \(G\) simplifies to just \(E\). Normalizing \(E\) with respect to the total number of particles in the system, we obtain \(\bar{E}(0,P,x_{Li},x_{Fe},x_O)\), where \(x_i = \frac{N_i}{N_{Li} + N_{Fe} + N_{O}}\). By taking the convex hull^{6} of \(E\) for all phases belonging to the Mcomponent system and projecting the stable nodes into the \((M  1)\)  dimension composition space, one can obtain the 0 K phase diagram for the closed system at constant pressure. The convex hull of a set of points is the smallest convex set containing the points. For instance, to construct a 0 K, closed \(\ce{LiFeO}\) system phase diagram, the convex hull is taken on the set of points in \((\bar{E},x_{Li},x_{Fe})\) space with \(x_O\) being related to the other composition variables by \(x_O = 1  x_{Li}  x_{Fe}\).
Grand Potential Phase Diagrams¶
In many scientific applications, the phase equilibria of interest is not that of a closed system, but rather one which is open to a particular element. For instance, synthesis might be carried out in flowing Argon gas at a particular temperature, which would mean the system is open to gaseous elements such as oxygen. Other experiments may be carried out in air, which essentially serves as an infinite reservoir of atmospheric elements such as oxygen, nitrogen and others.
In environments that are open to a particular element (say oxygen), the relevant control variable is the chemical potential of that element, \(\mu_{O}\). The relevant thermodynamic potential to study phase equilibria in an open system is the grand potential, which is defined as the following for an \(\ce{LiFeO}\) system that is open with respect to oxygen:
where the \(PV\) term is again ignored.
Normalizing with respect to \(\ce{LiFe}\) composition and dropping the explicit expression of the functional dependence of \(E\), \(S\), and \(N_{O}\) on the righthand side henceforth for brevity, we can obtain
where \(x_i = \frac{N_i}{N_{Li}+N_{Fe}}\) is the fraction of component \(i\) in \(\ce{LiFe}\) composition space.
To formally introduce temperature into ab initio phase stability calculations, one would usually need to take into account all the relevant excitations (e.g. vibrational, configurational, and electronic) that contribute to entropy. However, for systems where phase equilibria changes takes place primarily between solid phases with absorption or loss of as gas, a few simplifying assumptions can be made that allow us to obtain a useful approximate phase diagram with less effort. In such reactions, the reaction entropy is dominated by the entropy of the gas, and the effect of temperature is mostly captured by changes in the chemical potential of the gas. Since the \(TS\) term is the entropy contribution of the condensed system, it can be neglected compared to the entropy effect of the gas, and the expression for \(\bar{\phi}\) simplifies to:
Using the above assumption, the effect of temperature and partial pressure can be fully captured in a single chemical potential variable, with a more negative value corresponding to higher temperatures or lower partial pressures. The chemical potential is then treated as an external variable to obtain a grand potential phase diagram under a particular condition by taking the convex hull of \(\bar{\phi}\) for all phases and projecting the stable nodes into the \(\ce{LiFe}\) composition space.
Citation¶
To cite the Computational Phase Diagram App, please reference the following works:
 S. P. Ong, L. Wang, B. Kang, G. Ceder., The LiFePO2 Phase Diagram from First Principles Calculations, Chemistry of Materials, vol. 20, Mar. 2008, pp. 17981807.
 S.P. Ong, A. Jain, G. Hautier, B. Kang, and G. Ceder, Thermal stabilities of delithiated olivine MPO4 (M=Fe, Mn) cathodes investigated using first principles calculations, Electrochemistry Communications, vol. 12, 2010, pp. 427430.
Authors¶
 Shyue Ping Ong
 Anubhav Jain
 Matthew McDermott
Release Roadmap¶
Other Resources¶
References¶

https://doi.org/10.1021/je000236i ↩

https://doi.org/10.1126/sciadv.1600225 ↩

V. Raghavan, FeLiO Phase Diagram, ASM Alloy Phase Diagrams Center, P. Villars, editorinchief; H. Okamoto and K. Cenzual, section editors; http://www1.asminternational.org/AsmEnterprise/APD, ASM International, Materials Park, OH, 2006. ↩