Modeling equilibrium for an arbitrary rock bulk composition requires thermodynamic data and solution models for all the possible minerals, and software to compare the data and solve (minimize Gibbs energy) for the equilibrium mineral assemblage, including mineral compositions, as a function of temperature and pressure. Many geoscientists have contributed to the efforts to make this possible by collecting data, building solution models, and writing software. You can find citations to some of this work on the web pages, listed here, describing available software.
An example of what is possible with whole rock equilibrium calculations is shown in Figure 6.14, a T-P diagram of chemical reactions
Figure 6.14. Metamorphosed Shale Reactions. Petrogenetic grid based on the thermodynamic model dataset SPaC (version 4/2000) developed originally by Spear and Cheney (1989). Click on the image to see a larger version.
(K2O-FeO-MgO-Al2O3-SiO2-H2O). And continuous reactions are not shown. To see the continuous reactions that occur at equilibrium between the discontinuous reactions shown click here, to open Figure 15, mouse over the T-P diagram, and watch equilibrium assemblages and mineral compositions change on a Thompson AFM diagram.
Although equilibrium information for all bulk compositions is shown in Figure 6.14 with the associated AFM diagrams in Figure 6.15,
Figure 6.16. KFMASH Grid w/AFM Topologies and MAD Examples. Mineral assemblage diagram (MAD) based on the thermodynamic model dataset SPaC (version 4/2000) developed originally by Spear and Cheney (1989). Click on the image to see a larger, interactive version.
Using model calculations like those shown in Figures 6.15 and 6.16, equilibrium conditions for a metamorphic rock can be determined. But what do those temperatures and pressures mean in terms of the total history of the rock?