Christian de Capitani Version 11.03.2020

Example: Fig 2c from Green et al. (2016)

Green, E.C.R., White, R.W., Diener, J.F.A., Powell, R., Holland, T.J.B. & Palin, R.M. (2016).
Activity-composition relations for the calculation of partial melting equilibria in metabasic rocks. Journal of Metamorphic Geology, 34, 845-869.


Fig. 2c of Green et al. (2016)

Equilibrium phase diagram calculated with Domino (Version 11.03.2020)

The bulk is SM89, taken from Table 1 of Green et al. (2016), with excess H2O. So the comparison is only for the subsolidus region. 15 solutions were included, and all Fe3+ phases were considered.
The database is tcds62cid from the 11.03.2020 distribution, (iron and microcline excluded), If you want to test different solution models or databases, the THERIN-file is: (for copy/paste):

! -----Version: 03.11.2020
    500     18000
1   SI(52.47)TI(1.05)AL(18.20)FE(8.15)MG(12.71)CA(12.21)NA(5.22)K(0.46)O(1.47) H(100)SI(10)O(?)   *    SM89, Green, 2016

Minor differences are probably due to slightly different solution models used. (I have not checked the models, but the results are not different enough for MORB Geochemistry to be rewritten)

Version 11.03.2020 is faster than previous versions, because it stops the calculation, when minimum is only approximately reached. Further tests will show if this is reasonable.
There are only few obvious failures (point clouds. stuttering lines) These are mainly close to a narrow 3-amphibole field. (see clipping)
For labeling the code '-1' was used: labeling solutions with the dominant endmember. This will produce several (many) "apparent reactions" representing compositional variations of a solution, which is a minor inconvenience, but is much more robust when exsolutions are involved. e.g. reactions 20 and 5 (bottom center and left) are simply the change in epidote composition from cz-rich to ep-rich and reaction 128 (top right) is the CPXo composition varying from om-rich to di-rich and 147 the amphibole from parg to gl etc.

For more details, this example, including scripts and all output files can be downloaded here

The calculation with Domino took about 24 hours (MacBookPro14,3, 3.1 GHz Quad-Core Intel Core i7) and about 26 hours (MacBookPro13,1, 2 GHz Dual-Core Intel Core i5)
While you wait for the program to finish, here are some suggestions of what you might do in the meantime.