| for domcumentation of Perple_X thermodynamic data file structure and content refer to: | www.perplex.ethz.ch/perplex_thermodynamic_data_file.html | comments are indicated by the | character. | check for warnings at the end of the header section. | this file can only be read by Perple_X 6.7.8+ | for Harrison & Sverjensky's comments on the DEW data refer to the individual entries in DEW13ver.dat Harrison & Sverjensky HKF '13, Holland & Powell '02, Connolly & Kerrick '02 |<= data base title begin_standard_variables |<= name (<9 characters), reference value, tolerance P(bar) 1.00 0.1E-3 T(K) 298.15 0.1E-4 Y(CO2) 0.00 0.1E-6 mu(C1) 0.00 0.1E-2 mu(C2) 0.00 0.1E-2 end_standard_variables tolerance -.1E-2 |<= DTOL for unconstrained minimization, energy units HSC_conversion |<= tag enabling HSC to SUP apparent energy conversion, requires elemental entropies in the component list below begin_components | < 6 chars, molar weight (g), elemental entropy (R) Na2O 61.9790 205.175 MgO 40.3040 135.255 Al2O3 101.9610 364.425 SiO2 60.0840 223.96 K2O 94.1960 231.935 CaO 56.0770 144.205 TiO2 79.8660 235.87 MnO 70.9370 134.795 FeO 71.8440 129.855 NiO 74.6930 132.375 ZrO2 123.2200 244.33 Cl2 70.9060 223.08 O2 31.9990 205.15 H2O 18.0150 233.255 CO2 44.0100 210.89 CuO 79.5450 135.725 Cr2O3 151.9900 358.811 S2 64.130 64.1 F2 37.9968 202.79 end_components begin_makes | the make definitions segment of the thermodynamic | data file is preceeded by the "begin_makes" keyword, | and terminated by the "end_makes" keyword. Within | this segment comments must be preceeded by a "|" | character. | the make definitions define new thermodynamic entities | (e.g., fO2 buffers, endmembers etc) as a linear combination of | the "real" entities in the thermodynamic file. the definition | includes a "DQF" correction that is used to increment the | G of the linearly combined G's as a function of pressure and | temperature, i.e., the total G computed for a made entity is | the G of the right-hand-side of the make reaction + Gdqf. | the format assumes data on one line of less than 240 characters, | the expected format is | name = num_1 * name_1 + num_2 * name_2 ....+ num_int * name_int | dnum_1 dnum_2 dnum_3 | where i num_j is a number or fraction (i.e., two numbers separated | by a '/') and name_j are the names of the j existing entities. | and the dqf correction to the entity 'name' is | Gdqf(J/mol) = dnum_1 + T[K]*dnum_2 + P[bar]*dnum_3 | made entities are excluded from calculations in vertex if they | are included in the excluded phase list in the computational option | file (e.g., "in"); however, real entities to the right of the equals | sign in a make definition are not effected by the excluded phase list. | for example, if qL (liquid silica) is excluded, the make defintion | q8L = 4 qL | will function, but qL will not be considered as a possible phase in | the vertex calculation. | kyco silica buffer | kyco = 1 ky - 1 cor | 0 0 0 | sphene + cc + ru + q CO2 buffer | tcrq = 1 cc + 1 ru + 1 q - 1 sph | 0. 0. 0. | Ti phengite endmember for Auzanneau et al. (CMP 159:1-24, 2010) mica model | Gdqf corrected from -482876. to -13469.8 J/mol, D. Dolejs, Mar 6, 2010. | tip = 1 fcel + 1 geik - 1/2 fs | -13469.8 -14.694 .84 | mthm = 6 hem - 4 mt | mthm O2 buffer | 0 0 0 | qfm = 2 mt + 3 q -3 fa | qfm O2 buffer | 0 0 0 |emog = -1 fo - 2 gph + 1 mag + 1/2 en | EMOG O2 buffer |0 0 0 |ffm = -6 fa + 3 fs + 2 mt | ffm O2 buffer |0 0 0 |fff = 2 fa - 2 iron - 1 fs | iron fa fs O2 buffer |0 0 0 tbi = 1 phl - 1 fo + 1 ru + 1 q | Ti-biotite endmember for TiBio(HP) model 45e3 0 0 | tbi1 = 1 phl - 1 br + 1 ru | Ti-biotite endmember for TiBio(WPH) model | 78e3 0 0 | L Tajcamanova, 5/11/2006. | fbi = 1 east - 1/2 cor + 1/2 hem | Fe3+-biotite endmember for TiBio(WPH) model, definition corrected for | | presumed typo in the WPH (JMG, 2007) paper by L Tajcamanova, 5/11/2007. | 13.45e3 0 0 | fbi = 1 east - 1/2 cor + 1 hem, 8/10/2007, L Tajcamanova | Fe endmembers for ultramafic mineral solution | models, based on the talc Fe-Mg exchange. fbr = 1 br - 1/3 ta + 1/3 fta 0 0 0 fchum = 1 chum - 3 ta + 3 fta 0 0 0 fphA = 1 phA - 7/3 ta + 7/3 fta 0 0 0 fper = 1 per - 1/3 ta + 1/3 fta 0 0 0 fatg = 1 atg - 16 ta + 16 fta 0 0 0 atgts = 4 clin + 9/17 atg - 24/17 br | Tshermaks antigorite for Atg(PN) model (Padrón-Navarta et al., 2013, Lithos) -2e3 46.1 0 | ettore's makes: | odo_ef = 1/2 mag 1/2 cc | -1000. 0 0 | oank = 1/4 sid 1/2 cc 1/4 mag | -750 0 0 | odo = 1/2 mag + 1/2 cc | ordered dolmite species for H&P dolomite model | -13.5d3 0 0 | endmembers for the melt(HP) model, the dqf corrections | made in the solution model could be applied here, but | should not be done in both the thermodynamic data file | and the solution model. The small dqf corrections here | are made to avoid degeneracy between sil8L and silL etc. | These may not be necessary, and are definitely superfluous | if silL, foL and faL are excluded in calculations with | vertex. sil8L = 8/5 silL | sil8L is dqf corrected in the melt(HP) solution model -20 0 0 fo8L = 2 foL | fa8L is dqf corrected in the melt(HP) solution model -20 0 0 fa8L = 2 faL | fa8L is dqf corrected in the melt(HP) solution model -20 0 0 q8L = 4 qL | q8L has no dqf correction -20 0 0 | endmembers for the Diener et al. (JMG, 2007) amphibole models mrb = 1 gl - 2 jd +2 acm | magnesio riebekite, stoichiometry corrected 3/30/08 JADC 8d3 0. 0. ts_dqf = 1 ts | tschermakite 10000. 0. 0. parg_dqf = 1 parg | pargasite, Corrected from 10 kJ, 10/14/11, JADC 15000. 0. 0. gl_dqf = 1 gl | glaucophane, Corrected from 15 kJ, 10/14/11, MJC 3000. 0. 0. mpa = 1 parg -1 tr +1 anth | ortho-pargasite, stoichiometry corrected 3/30/08 JADC 27d3 0. 0. ged_dqf = 1 ged | gedrite 22000. 0. 0. ogl_dqf = 1 gl | ortho-glaucophane 15000. 0. 0. fanth_dq = 1 fanth | ferro-anthophyllite 7000. 0. 0. omrb_dqf = 1 gl -2 jd +2 acm | ortho-magnesio riebekite (mrb dqf + 25000 J/mol), stoichiometry corrected 3/30/08 JADC 33000. 0. 0. | Make definitions for the revised (Diener et al, JMG 2011 25:631-656) | ortho-amphibole model. Oct 31, 2011. mpa2 = 1 parg - 1 tr +1 anth 25d3 0. 0. ged_dqf2 = 1 ged 20000. 0. 0. | make definition for coma of the unpublished | thermocalc scapolite model coma = 3 ab + 1 cc 25000. 0. 0. | titano-clinohumite to fit weiss experiments. JADC, 2/08. | tichum = 4 fo + 1/2 br + 1/2 ru | -33581. 32. 0. tbit = 1 phl - 1 br + 1 ru |Ti-Fe3+ Biotite after Tajcmanova et al. (2009) 84e3 -11.5 0 |Model entered by Lucie Tajcmanova, September, 2008. fbit = 1 east - 1/2 cor + 1/2 hem |Ti-Fe3+ Biotite after Tajcmanova et al. (2009) 6e3 0 0 | Make definitions for the Smye et al (JMG, 2011, 28:753-768) | models for Mica, Carpholite and Chloritoid, to use these | make definitions delete the leading "|" comment characters: ocar = 1 fcar + 1/2 hem - 1/2 cor 45d3 0 0 mncar = 1 mcar + 1 mang - 1 per | corrected from 1 mcar + 1 mang - 1/2 cor, G Hoschek, 10/24/11. 30d3 0 0 octd = 1 fctd + 1/4 hem - 1/4 cor 125d2 0 0 fmu = 1 mu + 1/2 hem - 1/2 cor 30d3 0 0 | corrected from -30 kJ/mol, G Hoschek, 10/24/11. ma_dqf = 1 ma 3d3 0 0 | Make definition for the White et al (JMG, 2005) garnet, Gt(WPPH), JADC, 10/24/11. kho = 1 py - 1 gr + 1 andr 40d3 0 0 | Make definitions for the White et al (JMG, 2000) ilmenite, Ilm(WHP), JADC, 10/29/11. | Corrected coeffiecients for TC un-landau, formerly, the TC DQF was used without | modification. NOTE: the ilm and ilm_nol endmembers must be excluded from calculations | with the Ilm(WPH) model and the ilm_nol, oilm and dilm must be excluded from calculations | without the Ilm(WPH) model. JADC, 6/19/13. dilm = 1 ilm_nol | 1993 -2.1 0 => the TC DQF coeffecients, the values below are | oilm + {DH,R*ln(4),0} 15789.27763 -12.19977769 .1836386612d-1 | made definitions for Sapphirine, O/D non-ideal, Taylor-Jones & Powell | JMG 2010, 28:615-633. JADC, 6/26/2013. | ospr = 1 spr4 - 1 jd + 1 acm | 5d3 0 0 | Make definitions for Kelsey & Powell (JMG, 2011) Zr melting relations. JADC, 2/6/12 zr8L = 2 zrc | Zr-melt endmember 69e3 0 0 zrg = 1 gr + 3 zrc - 6 q | Zr garnet endmember 96e3 0.165e3 2.7 zrru = 1 bdy | Zr rutile endmember 52.84e3 -0.01829e3 -0.074 mots = 1 mgts - 1/2 cor + 1/2 hem | ferric opx endmember MgFe3+AlSiO6 for Opx(HP) 22d3 0 0 esn = 1 cats - 1/2 cor + 1/2 hem | essenite, ferric cpx endmember CaFe3+AlSiO6 for Cpx(HP) 22d3 0 0 end_makes th pdata produced at 16.35 on Wed 13 Feb,2002 (with sigma fit = 1.089) The Thermocalc data base includes data for ionic aqueous species, this data does not appear to be maintained and is present in hpaqver.dat JADC Mar 9, 2002. The Holland & Powell data base reproduced here has been augmented by shear moduli with sources as specified in Connolly & Kerrick EPSL '02. JADC Mar 9, 2004. see comments in hp96ver.dat for additional information on changes from the earlier hp94ver.dat and hp90ver.dat files. comments can be written in the data file between entries, provided the first character of comment card is left blank. WARNINGS: 1) The Holland & Powell thermodynamic data herein has been augmented by data for the Ghiorso et al. (2002, G3) pMELTS model. These data are not necessarily consistent and results obtained using the mixed data sources should be viewed with caution. The Ghirso et al data consists of the following melt (liquid) endmembers: Name Composition qGL Si4O8 modified 5/16/2012, see entry coGL Al4O6 modified 5/16/2012, see entry faGL Fe2SiO4 modified 5/16/2012, see entry foGL Mg2SiO4 modified 5/16/2012, see entry woGL Ca2Si2O6 modified 5/16/2012, see entry nasGL NaSi1/2O3/2 kalGL KAlSiO4 modified 5/16/2012, see entry h2oGL H2O tiGl TiO2 Notation (not alphabetical!) Entity Symbol Formula akermanite ak Ca2MgSi2O7 almandine alm Fe3Al2Si3O12 andalusite and Al2SiO5 andradite andr Ca3Fe2Si3O12 clinohumite chum Mg9Si4O16(OH)2 clinozoisite cz Ca2Al3Si3O12(OH) cordierite crd Mg2Al4Si5O18 epidote(ordered) ep Ca2FeAl2Si3O12(OH) fayalite fa Fe2SiO4 Fe-chloritoid fctd FeAl2SiO5(OH)2 Fe-cordierite fcrd Fe2Al4Si5O18 Fe-epidote fep Ca2Fe2AlSi3O12(OH) Fe-osumilite fosm KFe2Al5Si10O30 Fe-staurolite fst Fe4Al18Si7.5O48H4 forsterite fo Mg2SiO4 gehlenite geh Ca2Al2SiO7 grossular gr Ca3Al2Si3O12 hydrouscordierite hcrd Mg2Al4Si5O18H2O hydroxy-topaz tpz Al2SiO4(OH)2 kyanite ky Al2SiO5 larnite-bredigite larn Ca2SiO4 lawsonite law CaAl2Si2O7(OH)2H2O merwinite merw Ca3MgSi2O8 Mg-chloritoid mctd MgAl2SiO5(OH)2 Mg-staurolite mst Mg4Al18Si7.5O48H4 Mn-chloritoid mnctd MnAl2SiO5(OH)2 Mn-cordierite mncrd Mn2Al4Si5O18 Mn-staurolite mnst Mn4Al18Si7:5O48H4 monticellite mont CaMgSiO4 osumilite(1) osm1 KMg2Al5Si10O30 osumilite(2) osm2 KMg3Al3Si11O30 phase A phA Mg7Si2O8(OH)6 pumpellyite pump Ca4MgAl5Si6O21(OH)7 pyrope py Mg3Al2Si3O12 rankinite rnk Ca3Si2O7 sillimanite sill Al2SiO5 spessartine spss Mn3Al2Si3O12 sphene sph CaTiSiO5 spurrite spu Ca5Si2O8(CO3) tephroite teph Mn2SiO4 tilleyite ty Ca5Si2O7(CO3)2 vesuvianite vsv Ca19Mg2Al11Si18O69(OH)9 zircon zrc ZrSiO4 zoisite zo Ca2Al3Si3O12(OH) acmite acm NaFeSi2O6 Ca-tschermaks pyroxene cats CaAl2SiO6 Diopside di CaMgSi2O6 enstatite en Mg2Si2O6 ferrosilite fs Fe2Si2O6 hedenbergite hed CaFeSi2O6 jadeite jd NaAlSi2O6 mg-tschermak mgts MgAl2SiO6 pseudowollastonite pswo CaSiO3 pyroxmangite pxmn MnSiO3 rhodonite rhod MnSiO3 wollastonite wo CaSiO3 anthophyllite anth Mg7Si8O22(OH)2 cummingtonite cumm Mg7Si8O22(OH)2 Fe-anthophyllite fanth Fe7Si8O22(OH)2 Fe-glaucophane fgl Na2Fe3Al2Si8O22(OH)2 ferroactinolite ftr Ca2Fe5Si8O22(OH)2 gedrite(Na-free) ged Mg5Al4Si6O22(OH)2 glaucophane gl Na2Mg3Al2Si8O22(OH)2 grunerite grun Fe7Si8O22(OH)2 pargasite parg NaCa2Mg4Al3Si6O22(OH)2 riebeckite rieb Na2Fe5Si8O22(OH)2 tremolite tr Ca2Mg5Si8O22(OH)2 tschermakite ts Ca2Mg3Al4Si6O22(OH)2 deerite deer Fe18Si12O40(OH)10 fe-carpholite fcar FeAl2Si2O6(OH)4 fe-sapphirine(793) fspr Fe3.5Al9Si1.5O20 mg-carpholite mcar MgAl2Si2O6(OH)4 sapphirine(442) spr4 Mg4Al8Si2O20 sapphirine(793) spr7 Mg3.5Al9Si1.5O20 annite ann KFe3AlSi3O10(OH)2 celadonite cel KMgAlSi4O10(OH)2 eastonite east KMg2Al3Si2O10(OH)2 Fe-celadonite fcel KFeAlSi4O10(OH)2 margarite ma CaAl4Si2O10(OH)2 Mn-biotite mnbi KMn3AlSi3O10(OH)2 muscovite mu KAl3Si3O10(OH)2 Na-phlogopite naph NaMg3AlSi3O10(OH)2 paragonite pa NaAl3Si3O10(OH)2 phlogopite phl KMg3AlSi3O10(OH)2 Al-free chlorite afchl Mg6Si4O10(OH)8 amesite(14Ang) ames Mg4Al4Si2O10(OH)8 clinochlore(ordered) clin Mg5Al2Si3O10(OH)8 daphnite daph Fe5Al2Si3O10(OH)8 Fe-sudoite fsud Fe2Al4Si3O10(OH)8 Mn-chlorite mnchl Mn5Al2Si3O10(OH)8 Sudoite sud Mg2Al4Si3O10(OH)8 antigorite atg Mg48Si34O85(OH)62 chrysotile chr Mg3Si2O5(OH)2 Fe-talc fta Fe3Si4O10(OH)2 Kaolinite kao Al2Si2O5(OH)4 prehnite pre Ca2Al2Si3O10(OH)2 pyrophyllite prl Al2Si4O10(OH)2 talc ta Mg3Si4O10(OH)2 tschermak-talc tats Mg2Al2Si3O10(OH)2 albite ab NaAlSi3O8 analcite anl NaAlSi2O6H2O anorthite an CaAl2Si2O8 coesite coe SiO2 cristobalite crst SiO2 heulandite heu CaAl2Si7O186H2O highalbite abh NaAlSi3O8 kalsilite kals KAlSiO4 laumontite lmt CaAl2Si4O124H2O leucite lc KAlSi2O6 meionite me Ca4Al6Si6O24(CO3) microcline mic KAlSi3O8 nepheline ne NaAlSiO4 quartz q SiO2 sanidine san KAlSi3O8 stilbite stlb CaAl2Si7O187H2O stishovite stv SiO2 tridymite trd SiO2 wairakite wrk CaAl2Si4O12H2O3 baddeleyite bdy ZrO2 corundum cor Al2O3 geikielite geik MgTiO3 hematite hem Fe2O3 hercynite herc FeAl2O4 ilmenite ilm FeTiO3 lime lime CaO magnesioferrite mft MgFe2O4 magnetite mt Fe3O4 manganosite mang MnO nickel oxide NiO periclase per MgO pyrophanite pnt MnTiO3 rutile ru TiO2 spinel sp MgAl2O4 ulvospinel usp Fe2TiO4 brucite br Mg(OH)2 diaspore dsp AlO(OH) goethite gth FeO(OH) ankerite ank CaFe(CO3)2 aragonite arag CaCO3 calcite cc CaCO3 dolomite dol CaMg(CO3)2 magnesite mag MgCO3 rhodochrosite rhc MnCO3 siderite sid FeCO3 diamond diam C graphite gph C iron iron Fe nickel Ni Ni carbon dioxide CO2 carbon monoxide CO hydrogen H2 methane CH4 oxygen O2 water fluid H2O albite liquid abL NaAlSi3O8 anorthite liquid anL CaAl2Si2O8 diopside liquid diL CaMgSi2O6 enstatite liquid enL Mg2Si2O6 fayalite liquid faL Fe2SiO4 Fe-liquid (in KFMASH) fliq K3Fe0:5Al4Si19:5O47 Forsterite liquid foL Mg2SiO4 H2O liquid h2oL H2O H2O liquid (in KFMASH) hliq H2O K-feldspar liquid kspL KAlSi3O8 Mg liquid (in KFMASH) mliq K3Mg0:5Al4Si19:5O47 Silica liquid qL SiO2 Sillimanite liquid silL Al2SiO5 H+(aq) H+ Cl(aq) Cl- OH(aq) OH- Na+(aq) Na+ K+(aq) K+ Ca2+(aq) Ca++ Mg2+(aq) Mg++ Fe2+(aq) Fe++ Al3+(aq) Al+++ CO3--(aq) CO3 Al(OH)3(aq) AlOH3 Al(OH)4----(aq) AlOH4- KOH(aq) KOH HCl(aq) HCL KCl(aq) KCL NaCl(aq) NaCl CaCl(aq) CaCl2 CaCl+(aq) CaCl+ MgCl2(aq) MgCl2 MgCl+(aq) MgCl FeCl(aq) FeCl2 Aqueous silica aqSi SiO2 the following "end" card marks the end of the header section of the data file end fo EoS = 2 | H= -2172240. MgO(2)SiO2(1) G0 = -2053138 S0 = 95.1 V0 = 4.366 c1 = 233.3 c2 = .1494E-2 c3 = -603800 c5 = -1869.7 b1 = .613E-4 b5 = -.613E-3 b6 = 1250000. b7 = -187.5 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 end fa EoS = 2 | H= -1478160. SiO2(1)FeO(2) G0 = -1378945 S0 = 151 V0 = 4.631 c1 = 201.1 c2 = .1733E-1 c3 = -1960600 c5 = -900.9 b1 = .505E-4 b5 = -.505E-3 b6 = 1330000. b7 = -199.5 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end teph EoS = 2 | H= -1732300. SiO2(1)MnO(2) G0 = -1631725 S0 = 155.9 V0 = 4.899 c1 = 219.6 c3 = -1292700 c5 = -1308.3 b1 = .505E-4 b5 = -.505E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end lrn EoS = 2 | H= -2307160. SiO2(1)CaO(2) G0 = -2192411 S0 = 127.6 V0 = 5.16 c1 = 247.5 c2 = -.3206E-2 c5 = -2051.9 b1 = .505E-4 b5 = -.505E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 transition = 1 type = 4 t1 = 1710 t2 = 10.03 t3 = .5E-1 end mont EoS = 2 | H= -2252950. MgO(1)SiO2(1)CaO(1) G0 = -2135055 S0 = 108.1 V0 = 5.148 c1 = 250.7 c2 = -.10433E-1 c3 = -797200 c5 = -1996.1 b1 = .563E-4 b5 = -.563E-3 b6 = 1120000. b7 = -168 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end chum EoS = 2 | H= -9607380. MgO(9)SiO2(4)H2O(1) G0 = -9040340 S0 = 445 V0 = 19.785 c1 = 1070 c2 = -.16533E-1 c3 = -7899600 c5 = -7373.9 b1 = .61E-4 b5 = -.61E-3 b6 = 1290000. b7 = -193.5 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end py EoS = 2 | H= -6284750. MgO(3)Al2O3(1)SiO2(3) G0 = -5934105 S0 = 266.3 V0 = 11.318 c1 = 633.5 c3 = -5196100 c5 = -4315.2 b1 = .436E-4 b5 = -.436E-3 b6 = 1737000. b7 = -260.55 b8 = 4 m0 = 926900 m1 = 1.47 m2 = -87 end alm EoS = 2 | H= -5263630. Al2O3(1)SiO2(3)FeO(3) G0 = -4939788 S0 = 340 V0 = 11.511 c1 = 677.3 c3 = -3772700 c5 = -5044 b1 = .403E-4 b5 = -.403E-3 b6 = 1690000. b7 = -253.5 b8 = 4 m0 = 955400 m1 = 1.4 m2 = -109 end spss EoS = 2 | H= -5646410. Al2O3(1)SiO2(3)MnO(3) G0 = -5326387 S0 = 367 V0 = 11.792 c1 = 584.6 c2 = -.1593E-2 c3 = -7516700 c5 = -2750.1 b1 = .462E-4 b5 = -.462E-3 b6 = 1790000. b7 = -268.5 b8 = 4 m0 = 955400 m1 = 1.4 m2 = -109 end gr EoS = 2 | H= -6644130. Al2O3(1)SiO2(3)CaO(3) G0 = -6282110 S0 = 255 V0 = 12.535 c1 = 626 c3 = -5779200 c5 = -4002.9 b1 = .393E-4 b5 = -.393E-3 b6 = 1680000. b7 = -252 b8 = 4 m0 = 1070000 m1 = 1.6 m2 = -120 end andr EoS = 2 | H= -5768250. SiO2(3)CaO(3)FeO(2)O2(.5) G0 = -5425652 S0 = 318 V0 = 13.204 c1 = 638.6 c3 = -4955100 c5 = -3989.2 b1 = .393E-4 b5 = -.393E-3 b6 = 1590000. b7 = -238.5 b8 = 4 m0 = 1070000 m1 = 1.6 m2 = -120 end osm1 EoS = 2 | H= -.1496858E8 MgO(2)Al2O3(2.5)SiO2(10)K2O(.5) G0 = -14122760 S0 = 701 V0 = 37.893 c1 = 1625.8 c2 = -.35548E-1 c3 = -8063500 c5 = -13490.9 b1 = .76E-5 b5 = -.76E-4 b6 = 810000 b7 = -121.5 b8 = 4 m0 = 1070000 m1 = 1.6 m2 = -120 end osm2 EoS = 2 | H= -.1481075E8 MgO(3)Al2O3(1.5)SiO2(11)K2O(.5) G0 = -13973340 S0 = 724 V0 = 38.44 c1 = 1610.6 c2 = -.34457E-1 c3 = -8262099 c5 = -13128.8 b1 = .76E-5 b5 = -.76E-4 b6 = 810000 b7 = -121.5 b8 = 4 m0 = 1070000 m1 = 1.6 m2 = -120 end fosm EoS = 2 | H= -.1424892E8 Al2O3(2.5)SiO2(10)K2O(.5)FeO(2) G0 = -13424510 S0 = 762 V0 = 38.32 c1 = 1656 c2 = -.34163E-1 c3 = -6497700 c5 = -14114.3 b1 = .8E-5 b5 = -.8E-4 b6 = 800000 b7 = -120 b8 = 4 m0 = 1070000 m1 = 1.6 m2 = -120 end vsv EoS = 2 | H= -.4235251E8 MgO(2)Al2O3(5.5)SiO2(18)CaO(19)H2O(4.5) G0 = -39905380 S0 = 1890 V0 = 85.2 c1 = 4488 c2 = -.57952E-1 c3 = -22269300 c5 = -33478 b1 = .5E-4 b5 = -.5E-3 b6 = 1670000. b7 = -250.5 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end and EoS = 2 | H= -2588750. Al2O3(1)SiO2(1) G0 = -2440924 S0 = 92.7 V0 = 5.153 c1 = 277.3 c2 = -.6588E-2 c3 = -1914100 c5 = -2265.6 b1 = .411E-4 b5 = -.411E-3 b6 = 1334000. b7 = -200.1 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end ky EoS = 2 | H= -2593110. Al2O3(1)SiO2(1) G0 = -2442541 S0 = 83.5 V0 = 4.414 c1 = 279.4 c2 = -.7124E-2 c3 = -2055600 c5 = -2289.4 b1 = .404E-4 b5 = -.404E-3 b6 = 1590000. b7 = -238.5 b8 = 4 m0 = 966000 m1 = 2 m2 = -15 end sill EoS = 2 | H= -2585860. Al2O3(1)SiO2(1) G0 = -2438869 S0 = 95.5 V0 = 4.986 c1 = 280.2 c2 = -.69E-2 c3 = -1375700 c5 = -2399.4 b1 = .221E-4 b5 = -.221E-3 b6 = 1320000. b7 = -198 b8 = 4 m0 = 929000 m1 = 2 m2 = -15 transition = 1 type = 4 t1 = 2200 t2 = 4 t3 = .35E-1 end tpz EoS = 2 | H= -2904980. Al2O3(1)SiO2(1)H2O(1) G0 = -2689921 S0 = 100.5 V0 = 5.339 c1 = 387.7 c2 = -.712E-2 c3 = -857200 c5 = -3744.2 b1 = .404E-4 b5 = -.404E-3 b6 = 1315000. b7 = -197.25 b8 = 4 m0 = 966000 m1 = 2 m2 = -15 end mst EoS = 2 | H= -.2510146E8 MgO(4)Al2O3(9)SiO2(7.5)H2O(2) G0 = -23593330 S0 = 910 V0 = 44.26 c1 = 2820.5 c2 = -.59366E-1 c3 = -13774000 c5 = -24126 b1 = .12E-4 b5 = -.12E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end fst EoS = 2 | H= -.2375423E8 Al2O3(9)SiO2(7.5)FeO(4)H2O(2) G0 = -22282350 S0 = 1010 V0 = 44.88 c1 = 2880 c2 = -.56595E-1 c3 = -10642000 c5 = -25373 b1 = .12E-4 b5 = -.12E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end mnst EoS = 2 | H= -.2420382E8 Al2O3(9)SiO2(7.5)MnO(4)H2O(2) G0 = -22730480 S0 = 1024 V0 = 45.46 c1 = 2873.3 c2 = -.89064E-1 c3 = -12688000 c5 = -24749 b1 = .12E-4 b5 = -.12E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end mctd EoS = 2 | H= -3551430. MgO(1)Al2O3(1)SiO2(1)H2O(1) G0 = -3307815 S0 = 140 V0 = 6.875 c1 = 464.4 c2 = -.12654E-1 c3 = -1147200 c5 = -4341 b1 = .542E-4 b5 = -.542E-3 b6 = 1465000. b7 = -219.75 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end fctd EoS = 2 | H= -3215460. Al2O3(1)SiO2(1)FeO(1)H2O(1) G0 = -2977927 S0 = 155 V0 = 6.98 c1 = 484.6 c2 = -.13808E-1 c3 = -198900 c5 = -4762.2 b1 = .542E-4 b5 = -.542E-3 b6 = 1465000. b7 = -219.75 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end mnctd EoS = 2 | H= -3329280. Al2O3(1)SiO2(1)MnO(1)H2O(1) G0 = -3093616 S0 = 166 V0 = 7.175 c1 = 464.4 c2 = -.12654E-1 c3 = -1147200 c5 = -4341 b1 = .542E-4 b5 = -.542E-3 b6 = 1465000. b7 = -219.75 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end merw EoS = 2 | H= -4546260. MgO(1)SiO2(2)CaO(3) G0 = -4318804 S0 = 253.1 V0 = 9.847 c1 = 417.5 c2 = .8117E-2 c3 = -2923000 c5 = -2320.3 b1 = .615E-4 b5 = -.615E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end spu EoS = 2 | H= -5849370. SiO2(2)CaO(5)CO2(1) G0 = -5536279 S0 = 330 V0 = 14.697 c1 = 614.1 c2 = -.3508E-2 c3 = -2493100 c5 = -4168 b1 = .65E-4 b5 = -.65E-3 b6 = 950000 b7 = -142.5 b8 = 4 m0 = 575000 m1 = 2 m2 = -15 end zo EoS = 2 | H= -6898520. Al2O3(1.5)SiO2(3)CaO(2)H2O(.5) G0 = -6502908 S0 = 297 V0 = 13.575 c1 = 595.7 c2 = .62297E-1 c3 = -5921300 c5 = -3394.7 b1 = .67E-4 b5 = -.67E-3 b6 = 1120000. b7 = -168 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end cz EoS = 2 | H= -6898050. Al2O3(1.5)SiO2(3)CaO(2)H2O(.5) G0 = -6503630 S0 = 301 V0 = 13.63 c1 = 567 c2 = .18063E-1 c3 = -7034000 c5 = -2603 b1 = .46E-4 b5 = -.46E-3 b6 = 1120000. b7 = -168 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end fep EoS = 2 | H= -6002200. Al2O3(.5)SiO2(3)CaO(2)FeO(2)O2(.5)H2O(.5) G0 = -5625115 S0 = 357 V0 = 14.19 c1 = 520.1 c2 = .31499E-1 c3 = -15426000 c5 = 218.8 b1 = .505E-4 b5 = -.505E-3 b6 = 1294000. b7 = -194.1 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end ep EoS = 2 | H= -6463180. Al2O3(1)SiO2(3)CaO(2)FeO(1)O2(.25)H2O(.5) G0 = -6077129 S0 = 328 V0 = 13.91 c1 = 544.6 c2 = .24781E-1 c3 = -11230000 c5 = -1192.1 b1 = .505E-4 b5 = -.505E-3 b6 = 1233000. b7 = -184.95 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end law EoS = 2 | H= -4869070. Al2O3(1)SiO2(2)CaO(1)H2O(2) G0 = -4513272 S0 = 230 V0 = 10.132 c1 = 687.8 c2 = .1566E-2 c3 = 375900 c5 = -7179.2 b1 = .582E-4 b5 = -.582E-3 b6 = 1014000. b7 = -152.1 b8 = 4 m0 = 520000 m1 = .9613 m2 = -91 end pump EoS = 2 | H= -.143894E8 MgO(1)Al2O3(2.5)SiO2(6)CaO(4)H2O(3.5) G0 = -13448720 S0 = 629 V0 = 29.55 c1 = 1720.8 c2 = -.24928E-1 c3 = -5998700 c5 = -14620.3 b1 = .5E-4 b5 = -.5E-3 b6 = 1615000. b7 = -242.25 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end geh EoS = 2 | H= -3986850. Al2O3(1)SiO2(1)CaO(2) G0 = -3785608 S0 = 202 V0 = 9.024 c1 = 405.7 c2 = -.7099E-2 c3 = -1188300 c5 = -3174.4 b1 = .417E-4 b5 = -.417E-3 b6 = 1080000. b7 = -162 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 transition = 1 type = 4 t1 = 700 t2 = 11 t3 = .97E-1 end ak EoS = 2 | H= -3866300. MgO(1)SiO2(2)CaO(2) G0 = -3669742 S0 = 212.5 V0 = 9.254 c1 = 385.4 c2 = .3209E-2 c3 = -247500 c5 = -2889.9 b1 = .508E-4 b5 = -.508E-3 b6 = 1420000. b7 = -213 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end rnk EoS = 2 | H= -3944440. SiO2(2)CaO(3) G0 = -3744468 S0 = 210 V0 = 9.651 c1 = 372.3 c2 = -.2893E-2 c3 = -2462400 c5 = -2181.3 b1 = .65E-4 b5 = -.65E-3 b6 = 950000 b7 = -142.5 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end ty EoS = 2 | H= -6368190. SiO2(2)CaO(5)CO2(2) G0 = -6010097 S0 = 390 V0 = 17.039 c1 = 741.7 c2 = -.5345E-2 c3 = -1434600 c5 = -5878.5 b1 = .65E-4 b5 = -.65E-3 b6 = 950000 b7 = -142.5 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end crd EoS = 2 | H= -9163560. MgO(2)Al2O3(2)SiO2(5) G0 = -8653094 S0 = 407.5 V0 = 23.322 c1 = 821.3 c2 = .43339E-1 c3 = -8211200 c5 = -5000 b1 = .76E-5 b5 = -.76E-4 b6 = 810000 b7 = -121.5 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 transition = 1 type = 4 t1 = 1800 t2 = 20 t3 = .2 end hcrd EoS = 2 | H= -9447170. MgO(2)Al2O3(2)SiO2(5)H2O(1) G0 = -8890938 S0 = 487.3 V0 = 23.322 c1 = 869.7 c2 = .51995E-1 c3 = -7723700 c5 = -5251.2 b1 = .76E-5 b5 = -.76E-4 b6 = 810000 b7 = -121.5 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 transition = 1 type = 4 t1 = 1800 t2 = 20 t3 = .2 end fcrd EoS = 2 | H= -8436360. Al2O3(2)SiO2(5)FeO(2) G0 = -7949239 S0 = 475 V0 = 23.71 c1 = 851.5 c2 = .44724E-1 c3 = -6645000 c5 = -5623.4 b1 = .76E-5 b5 = -.76E-4 b6 = 810000 b7 = -121.5 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 transition = 1 type = 4 t1 = 1800 t2 = 20 t3 = .2 end mncrd EoS = 2 | H= -8681360. Al2O3(2)SiO2(5)MnO(2) G0 = -8191419 S0 = 475 V0 = 24.027 c1 = 847.7 c2 = .2849E-1 c3 = -7668200 c5 = -5311.4 b1 = .76E-5 b5 = -.76E-4 b6 = 810000 b7 = -121.5 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 transition = 1 type = 4 t1 = 1800 t2 = 20 t3 = .2 end | phase as modified for lars (commented out) | phA 1 0 0 0 H= -7130190. | 0. 7. 0. 2. 0. 0. 0. 0. 0. 0. 0. 0. 0. 3. 0. | -6580054. 390. 15.3 964. -.11521E-1 -4517800. | 0. -7724.7 0. 0. 0.826E-4 0. | 0. 0. -.826E-3 1514848. -217.5 4. | original HP phase A data: phA EoS = 2 |H= -7130190. MgO(7)SiO2(2)H2O(3) G0 = -6609954 S0 = 350 V0 = 15.442 c1 = 964 c2 = -.11521E-1 c3 = -4517800 c5 = -7724.7 b1 = .826E-4 b5 = -.826E-3 b6 = 1450000. b7 = -217.5 b8 = 4 m0 = 226000 m1 = 3.53 m2 = -104.8 end sph EoS = 2 | H= -2595490. SiO2(1)CaO(1)TiO2(1) G0 = -2454504 S0 = 131.2 V0 = 5.565 c1 = 233.7 c2 = .4043E-2 c3 = -2306500 c5 = -1207.6 b1 = .42E-4 b5 = -.42E-3 b6 = 1100000. b7 = -165 b8 = 4 m0 = 1121400 m1 = .78 m2 = -21 end zrc EoS = 2 | H= -2031820. SiO2(1)ZrO2(1) G0 = -1917278 S0 = 84.03 V0 = 3.926 c1 = 237 c2 = -.1788E-1 c3 = -149600 c5 = -2267.8 b1 = .222E-4 b5 = -.222E-3 b6 = 1160000. b7 = -174 b8 = 4 m0 = 1121400 m1 = .78 m2 = -21 end en EoS = 2 | H= -3090500. MgO(2)SiO2(2) G0 = -2915760 S0 = 132.5 V0 = 6.262 c1 = 356.2 c2 = -.299E-2 c3 = -596900 c5 = -3185.3 b1 = .505E-4 b5 = -.505E-3 b6 = 1070000. b7 = -160.5 b8 = 4 m0 = 757000 m1 = 1.5 m2 = -115 end | 0.757e6 2.0 -90. 0. 0. 0. 0. 0. 0. 0. => C/K values fs EoS = 2 | H= -2388650. SiO2(2)FeO(2) G0 = -2234453 S0 = 190.6 V0 = 6.592 c1 = 398.7 c2 = -.6579E-2 c3 = 1290100 c5 = -4058 b1 = .632E-4 b5 = -.632E-3 b6 = 1010000. b7 = -151.5 b8 = 4 m0 = 514000 m1 = 1.5 m2 = -115 end | 0.514e6 2.0 -90. 0. 0. 0. 0. 0. 0. 0. => C/K values mgts EoS = 2 | H= -3189350. MgO(1)Al2O3(1)SiO2(1) G0 = -3012610 S0 = 131 V0 = 5.9 c1 = 371.4 c2 = -.4082E-2 c3 = -398400 c5 = -3547.1 b1 = .508E-4 b5 = -.508E-3 b6 = 1143999. b7 = -171.6 b8 = 4 m0 = 757000 m1 = 1.5 m2 = -115 end di EoS = 2 | H= -3202720. MgO(1)SiO2(2)CaO(1) G0 = -3028353 S0 = 142.7 V0 = 6.619 c1 = 314.5 c2 = .41E-4 c3 = -2745900 c5 = -2020.1 b1 = .57E-4 b5 = -.57E-3 b6 = 1223000. b7 = -183.45 b8 = 4 m0 = 670000 m1 = 1.7 m2 = -100 end hed EoS = 2 | H= -2844360. SiO2(2)CaO(1)FeO(1) G0 = -2680995 S0 = 174.2 V0 = 6.795 c1 = 340.2 c2 = .812E-3 c3 = -1047800 c5 = -2646.7 b1 = .57E-4 b5 = -.57E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 610000 m1 = 1.7 m2 = -100 end jd EoS = 2 | H= -3027850. Na2O(.5)Al2O3(.5)SiO2(2) G0 = -2849148 S0 = 133.5 V0 = 6.04 c1 = 301.1 c2 = .10143E-1 c3 = -2239300 c5 = -2055.1 b1 = .466E-4 b5 = -.466E-3 b6 = 1284000. b7 = -192.6 b8 = 4 m0 = 840000 m1 = 1.7 m2 = -130 end acm EoS = 2 | H= -2586610. Na2O(.5)SiO2(2)FeO(1)O2(.25) G0 = -2419288 S0 = 170.6 V0 = 6.459 c1 = 307.1 c2 = .16758E-1 c3 = -1685500 c5 = -2125.8 b1 = .466E-4 b5 = -.466E-3 b6 = 1060000. b7 = -159 b8 = 4 m0 = 580000 m1 = 2 m2 = -10 end cats EoS = 2 | H= -3307000. Al2O3(1)SiO2(1)CaO(1) G0 = -3129678 S0 = 138 V0 = 6.356 c1 = 347.6 c2 = -.6974E-2 c3 = -1781600 c5 = -2757.5 b1 = .443E-4 b5 = -.443E-3 b6 = 1140000. b7 = -171 b8 = 4 m0 = 740000 m1 = 2 m2 = -20 end rhod EoS = 2 | H= -1321760. SiO2(1)MnO(1) G0 = -1244801 S0 = 100.5 V0 = 3.494 c1 = 138.4 c2 = .4088E-2 c3 = -1936000 c5 = -538.9 b1 = .508E-4 b5 = -.508E-3 b6 = 1250000. b7 = -187.5 b8 = 4 m0 = 740000 m1 = 2 m2 = -20 end pxmn EoS = 2 | H= -1322540. SiO2(1)MnO(1) G0 = -1245224 S0 = 99.3 V0 = 3.472 c1 = 138.4 c2 = .4088E-2 c3 = -1936000 c5 = -538.9 b1 = .508E-4 b5 = -.508E-3 b6 = 1250000. b7 = -187.5 b8 = 4 m0 = 740000 m1 = 2 m2 = -20 end wo EoS = 2 | H= -1634110. SiO2(1)CaO(1) G0 = -1548917 S0 = 82.5 V0 = 3.993 c1 = 159.3 c3 = -967300 c5 = -1075.4 b1 = .46E-4 b5 = -.46E-3 b6 = 795000 b7 = -119.25 b8 = 4 m0 = 620000 m1 = 1.7 m2 = -100 end pswo EoS = 2 | H= -1627740. SiO2(1)CaO(1) G0 = -1544246 S0 = 88.2 V0 = 4.008 c1 = 157.8 c3 = -967300 c5 = -1075.4 b1 = .539E-4 b5 = -.539E-3 b6 = 1050000. b7 = -157.5 b8 = 4 m0 = 620000 m1 = 1.7 m2 = -100 end tr EoS = 2 | H= -.1231023E8 MgO(5)SiO2(8)CaO(2)H2O(1) G0 = -11582670 S0 = 550 V0 = 27.27 c1 = 1260.2 c2 = .383E-2 c3 = -11455000 c5 = -8237.601 b1 = .534E-4 b5 = -.534E-3 b6 = 762000.1 b7 = -114.3 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end ftr EoS = 2 | H= -.1050391E8 SiO2(8)CaO(2)FeO(5)H2O(1) G0 = -9833597 S0 = 715 V0 = 28.28 c1 = 1290 c2 = .29991E-1 c3 = -8447500 c5 = -8947 b1 = .534E-4 b5 = -.534E-3 b6 = 760000 b7 = -114 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end ts EoS = 2 | H= -.1254101E8 MgO(3)Al2O3(2)SiO2(6)CaO(2)H2O(1) G0 = -11808850 S0 = 545 V0 = 26.8 c1 = 1244.8 c2 = .24348E-1 c3 = -11965000 c5 = -8112.1 b1 = .534E-4 b5 = -.534E-3 b6 = 760000 b7 = -114 b8 = 4 m0 = 524500 m1 = 1.012085 m2 = -81.5 end parg EoS = 2 | H= -.127206E8 Na2O(.5)MgO(4)Al2O3(1.5)SiO2(6)CaO(2)H2O(1) G0 = -11988550 S0 = 601 V0 = 27.19 c1 = 1280.2 c2 = .22997E-1 c3 = -12359500 c5 = -8065.8 b1 = .534E-4 b5 = -.534E-3 b6 = 912000 b7 = -136.8 b8 = 4 m0 = 432000 m1 = .8266 m2 = -67 end gl EoS = 2 | H= -.1196964E8 Na2O(1)MgO(3)Al2O3(1)SiO2(8)H2O(1) G0 = -11234430 S0 = 535 V0 = 26.05 c1 = 1717.5 c2 = -.12107 c3 = 7075000 c5 = -19272 b1 = .53E-4 b5 = -.53E-3 b6 = 882999.9 b7 = -132.45 b8 = 4 m0 = 696000 m1 = 1.20833 m2 = -94 end fgl EoS = 2 | H= -.1089164E8 Na2O(1)Al2O3(1)SiO2(8)FeO(3)H2O(1) G0 = -10187790 S0 = 624 V0 = 26.59 c1 = 1762.9 c2 = -.118992 c3 = 9423700 c5 = -20207.1 b1 = .53E-4 b5 = -.53E-3 b6 = 890000 b7 = -133.5 b8 = 4 m0 = 696000 m1 = 1.20833 m2 = -94 end rieb EoS = 2 | H= -.1003511E8 Na2O(1)SiO2(8)FeO(5)O2(.5)H2O(1) G0 = -9351875 S0 = 691 V0 = 27.49 c1 = 1746.9 c2 = -.113572 c3 = 9370300 c5 = -19468.7 b1 = .53E-4 b5 = -.53E-3 b6 = 890000 b7 = -133.5 b8 = 4 m0 = 696000 m1 = 1.20833 m2 = -94 end anth EoS = 2 | H= -.1206926E8 MgO(7)SiO2(8)H2O(1) G0 = -11342860 S0 = 536 V0 = 26.54 c1 = 1277.3 c2 = .25825E-1 c3 = -9704600 c5 = -9074.7 b1 = .5E-4 b5 = -.5E-3 b6 = 700000.1 b7 = -105 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end fanth EoS = 2 | H= -9627250. SiO2(8)FeO(7)H2O(1) G0 = -8968476 S0 = 725 V0 = 27.87 c1 = 1383.1 c2 = .30669E-1 c3 = -4224700 c5 = -11257.6 b1 = .5E-4 b5 = -.5E-3 b6 = 700000.1 b7 = -105 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end cumm EoS = 2 | H= -.1205757E8 MgO(7)SiO2(8)H2O(1) G0 = -11332370 S0 = 540 V0 = 26.33 c1 = 1277.3 c2 = .25825E-1 c3 = -9704600 c5 = -9074.7 b1 = .5E-4 b5 = -.5E-3 b6 = 700000.1 b7 = -105 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end grun EoS = 2 | H= -9613260. SiO2(8)FeO(7)H2O(1) G0 = -8955976 S0 = 730 V0 = 27.84 c1 = 1383.1 c2 = .30669E-1 c3 = -4224700 c5 = -11257.6 b1 = .5E-4 b5 = -.5E-3 b6 = 648000 b7 = -97.2 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end ged EoS = 2 | H= -.1231988E8 MgO(5)Al2O3(2)SiO2(6)H2O(1) G0 = -11584120 S0 = 515 V0 = 25.8 c1 = 1307.7 c2 = .23642E-1 c3 = -9307400 c5 = -9799 b1 = .48E-4 b5 = -.48E-3 b6 = 770000 b7 = -115.5 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end spr4 EoS = 2 | H= -.1101977E8 MgO(4)Al2O3(4)SiO2(2) G0 = -10422240 S0 = 443 V0 = 19.905 c1 = 1160.3 c2 = -.24324E-1 c3 = -7706600 c5 = -8974.2 b1 = .49E-4 b5 = -.49E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 end spr7 EoS = 2 | H= -.1106577E8 MgO(3.5)Al2O3(4.5)SiO2(1.5) G0 = -10468950 S0 = 448 V0 = 19.775 c1 = 1167.9 c2 = -.2487E-1 c3 = -7607300 c5 = -9155.3 b1 = .49E-4 b5 = -.49E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 810000 m1 = 1.82 m2 = -140 end fspr EoS = 2 | H= -9835830. Al2O3(4.5)SiO2(1.5)FeO(3.5) G0 = -9275354 S0 = 551 V0 = 20.391 c1 = 1257.8 c2 = -.22171E-1 c3 = -1664000 c5 = -11348.4 b1 = .49E-4 b5 = -.49E-3 b6 = 1200000. b7 = -180 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end mcar EoS = 2 | H= -4781270. MgO(1)Al2O3(1)SiO2(2)H2O(2) G0 = -4422178 S0 = 210 V0 = 10.59 c1 = 667.8 c2 = -.12559E-1 c3 = -1167100 c5 = -6440 b1 = .5E-4 b5 = -.5E-3 b6 = 525000 b7 = -78.75 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end fcar EoS = 2 | H= -4413440. Al2O3(1)SiO2(2)FeO(1)H2O(2) G0 = -4069375 S0 = 255 V0 = 10.69 c1 = 674.8 c2 = -.10092E-1 c3 = -715800 c5 = -6554.5 b1 = .5E-4 b5 = -.5E-3 b6 = 525000 b7 = -78.75 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end deer EoS = 2 | H= -.1834767E8 SiO2(12)FeO(18)O2(1.5)H2O(5) G0 = -16901570 S0 = 1650 V0 = 55.74 c1 = 3164.4 c2 = -.27883E-1 c3 = -5039100 c5 = -26721 b1 = .5E-4 b5 = -.5E-3 b6 = 630000 b7 = -94.5 b8 = 4 m0 = 510000 m1 = .62 m2 = -108 end mu EoS = 2 | H= -5984180. Al2O3(1.5)SiO2(3)K2O(.5)H2O(1) G0 = -5603722 S0 = 292 V0 = 14.083 c1 = 756.4 c2 = -.1984E-1 c3 = -2170000 c5 = -6979.2 b1 = .596E-4 b5 = -.596E-3 b6 = 490000 b7 = -73.5 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end cel EoS = 2 | H= -5842450. MgO(1)Al2O3(.5)SiO2(4)K2O(.5)H2O(1) G0 = -5462950 S0 = 290 V0 = 13.957 c1 = 741.2 c2 = -.18748E-1 c3 = -2368800 c5 = -6616.9 b1 = .596E-4 b5 = -.596E-3 b6 = 700000.1 b7 = -105 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end fcel EoS = 2 | H= -5477590. Al2O3(.5)SiO2(4)K2O(.5)FeO(1)H2O(1) G0 = -5111327 S0 = 329 V0 = 14.07 c1 = 756.3 c2 = -.19147E-1 c3 = -1586100 c5 = -6928.7 b1 = .596E-4 b5 = -.596E-3 b6 = 700000.1 b7 = -105 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end pa EoS = 2 | H= -5946330. Na2O(.5)Al2O3(1.5)SiO2(3)H2O(1) G0 = -5565091 S0 = 276 V0 = 13.211 c1 = 803 c2 = -.3158E-1 c3 = 217000 c5 = -8151 b1 = .774E-4 b5 = -.774E-3 b6 = 550000 b7 = -82.5 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end ma EoS = 2 | H= -6241150. Al2O3(2)SiO2(2)CaO(1)H2O(1) G0 = -5857267 S0 = 267 V0 = 12.964 c1 = 744.4 c2 = -.168E-1 c3 = -2074400 c5 = -6783.2 b1 = .487E-4 b5 = -.487E-3 b6 = 1300000. b7 = -195 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end phl EoS = 2 | H= -6219210. MgO(3)Al2O3(.5)SiO2(3)K2O(.5)H2O(1) G0 = -5837160 S0 = 328 V0 = 14.964 c1 = 770.3 c2 = -.36939E-1 c3 = -2328900 c5 = -6531.6 b1 = .579E-4 b5 = -.579E-3 b6 = 513000 b7 = -76.95 b8 = 4 m0 = 274000 m1 = 2 m2 = -50 end ann EoS = 2 | H= -5151860. Al2O3(.5)SiO2(3)K2O(.5)FeO(3)H2O(1) G0 = -4801474 S0 = 418 V0 = 15.432 c1 = 815.7 c2 = -.34861E-1 c3 = 19800 c5 = -7466.7 b1 = .579E-4 b5 = -.579E-3 b6 = 513000 b7 = -76.95 b8 = 4 m0 = 274000 m1 = 2 m2 = -50 end mnbi EoS = 2 | H= -5462860. Al2O3(.5)SiO2(3)K2O(.5)MnO(3)H2O(1) G0 = -5112715 S0 = 433 V0 = 15.264 c1 = 809.9 c2 = -.59213E-1 c3 = -1514400 c5 = -6998.7 b1 = .579E-4 b5 = -.579E-3 b6 = 513000 b7 = -76.95 b8 = 4 m0 = 274000 m1 = 2 m2 = -50 end east EoS = 2 | H= -6338210. -5951625. MgO(2)Al2O3(1.5)SiO2(2)K2O(.5)H2O(1) G0 = -5951625 S0 = 318 V0 = 14.738 c1 = 785.5 c2 = -.38031E-1 c3 = -2130300 c5 = -6893.7 b1 = .579E-4 b5 = -.579E-3 b6 = 513000 b7 = -76.95 b8 = 4 m0 = 274000 m1 = 2 m2 = -50 end naph EoS = 2 | H= -6172870. Na2O(.5)MgO(3)Al2O3(.5)SiO2(3)H2O(1) G0 = -5791828 S0 = 318 V0 = 14.45 c1 = 773.5 c2 = -.40229E-1 c3 = -2597900 c5 = -6512.6 b1 = .579E-4 b5 = -.579E-3 b6 = 513000 b7 = -76.95 b8 = 4 m0 = 274000 m1 = 2 m2 = -50 end clin EoS = 2 | H= -8912480. MgO(5)Al2O3(1)SiO2(3)H2O(4) G0 = -8251890 S0 = 430.5 V0 = 21.09 c1 = 1161.8 c2 = .10133E-1 c3 = -7657300 c5 = -9690.899 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end ames EoS = 2 | H= -9034470. MgO(4)Al2O3(2)SiO2(2)H2O(4) G0 = -8366214 S0 = 410 V0 = 20.52 c1 = 1177 c2 = .9041E-2 c3 = -7458700 c5 = -10053 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end afchl EoS = 2 | H= -8727500. MgO(6)SiO2(4)H2O(4) G0 = -8067718 S0 = 428 V0 = 21.66 c1 = 1146.6 c2 = .11225E-1 c3 = -7855900 c5 = -9328.8 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end daph EoS = 2 | H= -7135350. Al2O3(1)SiO2(3)FeO(5)H2O(4) G0 = -6522911 S0 = 565 V0 = 21.34 c1 = 1237.4 c2 = .13594E-1 c3 = -3743000 c5 = -11250 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end mnchl EoS = 2 | H= -7666560. Al2O3(1)SiO2(3)MnO(5)H2O(4) G0 = -7056014 S0 = 595 V0 = 22.59 c1 = 1227.8 c2 = -.2699E-1 c3 = -6299800 c5 = -10469.4 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end sud EoS = 2 | H= -8626230. MgO(2)Al2O3(2)SiO2(3)H2O(4) G0 = -7970064 S0 = 404 V0 = 20.3 c1 = 1436.1 c2 = -.48749E-1 c3 = -2748500 c5 = -13764 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end fsud EoS = 2 | H= -7904130. Al2O3(2)SiO2(3)FeO(2)H2O(4) G0 = -7268477 S0 = 462 V0 = 20.4 c1 = 1466.3 c2 = -.47365E-1 c3 = -1182800 c5 = -14388 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end prl EoS = 2 | H= -5640620. Al2O3(1)SiO2(4)H2O(1) G0 = -5266608 S0 = 239.4 V0 = 12.81 c1 = 784.5 c2 = -.42948E-1 c3 = 1251000 c5 = -8495.9 b1 = .75E-4 b5 = -.75E-3 b6 = 525000 b7 = -78.75 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end ta EoS = 2 | H= -5897150. MgO(3)SiO2(4)H2O(1) G0 = -5516955 S0 = 260 V0 = 13.625 c1 = 622.2 c3 = -6385500 c5 = -3916.3 b1 = .37E-4 b5 = -.37E-3 b6 = 480000 b7 = -72 b8 = 4 m0 = 226000 m1 = 3.53 m2 = -104.8 end fta EoS = 2 | H= -4803160. SiO2(4)FeO(3)H2O(1) G0 = -4455225 S0 = 352 V0 = 14.225 c1 = 579.7 c2 = .39494E-1 c3 = -6459300 c5 = -3088.1 b1 = .37E-4 b5 = -.37E-3 b6 = 480000 b7 = -72 b8 = 4 m0 = 226000 m1 = 3.53 m2 = -104.8 end tats EoS = 2 | H= -5987760. MgO(2)Al2O3(1)SiO2(3)H2O(1) G0 = -5605713 S0 = 259 V0 = 13.51 c1 = 549.5 c2 = .36324E-1 c3 = -8606600 c5 = -2515.3 b1 = .37E-4 b5 = -.37E-3 b6 = 480000 b7 = -72 b8 = 4 m0 = 226000 m1 = 3.53 m2 = -104.8 end kao EoS = 2 | H= -4122140. Al2O3(1)SiO2(2)H2O(2) G0 = -3801504 S0 = 203.7 V0 = 9.934 c1 = 436.7 c2 = -.34295E-1 c3 = -4055900 c5 = -2699.1 b1 = .51E-4 b5 = -.51E-3 b6 = 645000 b7 = -96.75 b8 = 4 m0 = 226000 m1 = 3.53 m2 = -104.8 end pre EoS = 2 | H= -6203140. Al2O3(1)SiO2(3)CaO(2)H2O(1) G0 = -5825834 S0 = 292.8 V0 = 14.026 c1 = 724.9 c2 = -.13865E-1 c3 = -2059000 c5 = -6323.9 b1 = .51E-4 b5 = -.51E-3 b6 = 835000 b7 = -125.25 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end chr EoS = 2 | H= -4359130. MgO(3)SiO2(2)H2O(2) G0 = -4031415 S0 = 221.3 V0 = 10.746 c1 = 624.7 c2 = -.2077E-1 c3 = -1721800 c5 = -5619.4 b1 = .47E-4 b5 = -.47E-3 b6 = 525000 b7 = -78.75 b8 = 4 m0 = 146364 m1 = 5.37 m2 = -88.8 end atg EoS = 2 | H= -.7142492E8 MgO(48)SiO2(34)H2O(31) G0 = -66132440 S0 = 3591 V0 = 175.48 c1 = 9621 c2 = -.91183E-1 c3 = -35941600 c5 = -83034.2 b1 = .47E-4 b5 = -.47E-3 b6 = 525000 b7 = -78.75 b8 = 4 m0 = 146364 m1 = 5.37 m2 = -88.8 end ab EoS = 2 | H= -3934530. Na2O(.5)Al2O3(.5)SiO2(3) G0 = -3711877 S0 = 210.1 V0 = 10.006 c1 = 452 c2 = -.13364E-1 c3 = -1275900 c5 = -3953.6 b1 = .456E-4 b5 = -.456E-3 b6 = 593000 b7 = -88.95 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 transition = 1 type = 4 t1 = 950 t2 = 16 t3 = .124 end abh EoS = 2 | H= -3924770. Na2O(.5)Al2O3(.5)SiO2(3) G0 = -3706083 S0 = 223.4 V0 = 10.109 c1 = 452 c2 = -.13364E-1 c3 = -1275900 c5 = -3953.6 b1 = .456E-4 b5 = -.456E-3 b6 = 593000 b7 = -88.95 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 end mic EoS = 2 | H= -3975100. Al2O3(.5)SiO2(3)K2O(.5) G0 = -3750217 S0 = 216 V0 = 10.892 c1 = 448.8 c2 = -.10075E-1 c3 = -1007300 c5 = -3973.1 b1 = .335E-4 b5 = -.335E-3 b6 = 574000 b7 = -86.1 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 end san EoS = 2 | H= -3964950. Al2O3(.5)SiO2(3)K2O(.5) G0 = -3744241 S0 = 230 V0 = 10.9 c1 = 448.8 c2 = -.10075E-1 c3 = -1007300 c5 = -3973.1 b1 = .335E-4 b5 = -.335E-3 b6 = 574000 b7 = -86.1 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 end an EoS = 2 | H= -4233420. Al2O3(1)SiO2(2)CaO(1) G0 = -4007795 S0 = 200 V0 = 10.079 c1 = 371.6 c2 = .12615E-1 c3 = -4110200 c5 = -2038.4 b1 = .238E-4 b5 = -.238E-3 b6 = 919000 b7 = -137.85 b8 = 4 m0 = 421000 m1 = 3.48 m2 = -43 transition = 1 type = 4 t1 = 2300 t2 = 11 t3 = .5E-1 end q EoS = 2 | H= -910840. SiO2(1) G0 = -856424.6 S0 = 41.5 V0 = 2.2688 c1 = 110.7 c2 = -.5189E-2 c5 = -1128.3 b1 = .65E-5 b5 = -.65E-4 b6 = 750000 b7 = -112.5 b8 = 4 m0 = 431250 m1 = .46 m2 = -14 transition = 1 type = 4 t1 = 847 t2 = 4.95 t3 = .1188 end trd EoS = 2 | H= -906690. SiO2(1) G0 = -853646.1 S0 = 46.1 V0 = 2.7 c1 = 97.9 c2 = -.335E-2 c3 = -636200 c5 = -774 b1 = .5E-5 b5 = -.5E-4 b6 = 750000 b7 = -112.5 b8 = 4 m0 = 431250 m1 = .46 m2 = -14 end crst EoS = 2 | H= -906000. SiO2(1) G0 = -853075.4 S0 = 46.5 V0 = 2.61 c1 = 97.9 c2 = -.335E-2 c3 = -636200 c5 = -774 b1 = .81E-5 b5 = -.81E-4 b6 = 600000 b7 = -90 b8 = 4 m0 = 431250 m1 = .46 m2 = -14 end coe EoS = 2 | H= -905490. SiO2(1) G0 = -850865.9 S0 = 40.8 V0 = 2.064 c1 = 96.5 c2 = -.577E-3 c3 = -444800 c5 = -798.2 b1 = .18E-4 b5 = -.18E-3 b6 = 999999.5 b7 = -150 b8 = 4 m0 = 620000 m1 = 1.5 m2 = -150 end | 0.59d6 1.40097 -24. 0. 0. 0. 0. 0. 0. 0. helffrich www1.gly.bris.ac.uk/~ge stv EoS = 2 | H= -875670. SiO2(1) G0 = -816186.1 S0 = 24.5 V0 = 1.4014 c1 = 68.1 c2 = .601E-2 c3 = -1978200 c5 = -82.1 b1 = .25E-4 b5 = -.25E-3 b6 = 3160000. b7 = -474 b8 = 4 m0 = 2170000 m1 = 1.8 m2 = -180 end | 2.26E+06 1.8 -24.0 0. 0. 0. 0. 0. 0. 0. CK '02 ne EoS = 2 | H= -2095070. Na2O(.5)Al2O3(.5)SiO2(1) G0 = -1980443 S0 = 124.4 V0 = 5.419 c1 = 272.7 c2 = -.12398E-1 c5 = -2763.1 b1 = .81E-4 b5 = -.81E-3 b6 = 600000 b7 = -90 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 transition = 1 type = 4 t1 = 467 t2 = 10 t3 = .8E-1 end kals EoS = 2 | H= -2121900. Al2O3(.5)SiO2(1)K2O(.5) G0 = -2006146 S0 = 134 V0 = 6.04 c1 = 242 c2 = -.4482E-2 c3 = -895800 c5 = -1935.8 b1 = .576E-4 b5 = -.576E-3 b6 = 589999.9 b7 = -88.5 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 end lc EoS = 2 | H= -3029170. Al2O3(.5)SiO2(2)K2O(.5) G0 = -2866305 S0 = 200 V0 = 8.828 c1 = 369.8 c2 = -.16332E-1 c3 = 684700 c5 = -3683.1 b1 = .367E-4 b5 = -.367E-3 b6 = 630000 b7 = -94.5 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 transition = 1 type = 4 t1 = 938 t2 = 18 t3 = .482 end me EoS = 2 | H= -.1384337E8 Al2O3(3)SiO2(6)CaO(4)CO2(1) G0 = -13105920 S0 = 752 V0 = 33.985 c1 = 1359 c2 = .36442E-1 c3 = -8594700 c5 = -9598.2 b1 = .316E-4 b5 = -.316E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 end wrk EoS = 2 | H= -6666380. Al2O3(1)SiO2(4)CaO(1)H2O(2) G0 = -6220237 S0 = 375 V0 = 19.04 c1 = 838.3 c2 = -.2146E-1 c3 = -2272000 c5 = -7292.3 b1 = .238E-4 b5 = -.238E-3 b6 = 999999.5 b7 = -150 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end lmt EoS = 2 | H= -7268410. Al2O3(1)SiO2(4)CaO(1)H2O(4) G0 = -6707599 S0 = 457 V0 = 20.37 c1 = 1013.4 c2 = -.21413E-1 c3 = -2235800 c5 = -8806.7 b1 = .238E-4 b5 = -.238E-3 b6 = 999999.5 b7 = -150 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end heu EoS = 2 | H= -.105991E8 Al2O3(1)SiO2(7)CaO(1)H2O(6) G0 = -9762014 S0 = 669 V0 = 31.8 c1 = 1504.8 c2 = -.33224E-1 c3 = -2959300 c5 = -13297.2 b1 = .238E-4 b5 = -.238E-3 b6 = 999999.5 b7 = -150 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end stlb EoS = 2 | H= -.1089895E8 Al2O3(1)SiO2(7)CaO(1)H2O(7) G0 = -10004530 S0 = 710 V0 = 32.87 c1 = 1588.4 c2 = -.32043E-1 c3 = -3071600 c5 = -13966.9 b1 = .238E-4 b5 = -.238E-3 b6 = 999999.5 b7 = -150 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end anl EoS = 2 | H= -3309960. Na2O(.5)Al2O3(.5)SiO2(2)H2O(1) G0 = -3091067 S0 = 232 V0 = 9.74 c1 = 643.5 c2 = -.16067E-1 c3 = 9302300 c5 = -9179.6 b1 = .5E-4 b5 = -.5E-3 b6 = 400000 b7 = -60 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end lime EoS = 2 | H= -634940. CaO(1) G0 = -603297.3 S0 = 38.1 V0 = 1.6764 c1 = 52.4 c2 = .3673E-2 c3 = -750700 c5 = -51 b1 = .665E-4 b5 = -.6649999E-3 b6 = 1160000. b7 = -174 b8 = 4 m0 = 612000 m1 = .9714 m2 = -106 end ru EoS = 2 | H= -944190. TiO2(1) G0 = -888963.7 S0 = 50.6 V0 = 1.882 c1 = 90.4 c2 = .29E-2 c5 = -623.8 b1 = .443E-4 b5 = -.443E-3 b6 = 2224999. b7 = -333.75 b8 = 4 m0 = 1121400 m1 = .78 m2 = -21 end per EoS = 2 | H= -601590. MgO(1) G0 = -569276.5 S0 = 26.9 V0 = 1.125 c1 = 60.5 c2 = .362E-3 c3 = -535800 c5 = -299.2 b1 = .62E-4 b5 = -.62E-3 b6 = 1650000. b7 = -247.5 b8 = 4 m0 = 1121400 m1 = .78 m2 = -21 end mang EoS = 2 | H= -385150. MnO(1) G0 = -362815.6 S0 = 59.7 V0 = 1.3221 c1 = 59.8 c2 = .36E-2 c3 = -31400 c5 = -282.6 b1 = .63E-4 b5 = -.63E-3 b6 = 1640000. b7 = -246 b8 = 4 m0 = 1121400 m1 = .78 m2 = -21 end cor EoS = 2 | H= -1675210. Al2O3(1) G0 = -1581710 S0 = 50.9 V0 = 2.558 c1 = 139.5 c2 = .589E-2 c3 = -2460600 c5 = -589.2 b1 = .419E-4 b5 = -.419E-3 b6 = 2520000. b7 = -378 b8 = 4 m0 = 1620000 m1 = 1.8 m2 = -190 end hem EoS = 2 | H= -825720. FeO(2)O2(.5) G0 = -743740.7 S0 = 87.39999 V0 = 3.0274 c1 = 163.9 c3 = -2257200 c5 = -657.6 b1 = .599E-4 b5 = -.599E-3 b6 = 1996000. b7 = -299.4 b8 = 4 m0 = 915000 m1 = 1.2 m2 = -130 transition = 1 type = 4 t1 = 955 t2 = 15.6 end NiO EoS = 2 | H= -239440. NiO(1) G0 = -211272.7 S0 = 38 V0 = 1.097 c1 = 47.7 c2 = .7824E-2 c3 = -392500 b1 = .62E-4 b5 = -.62E-3 b6 = 1650000. b7 = -247.5 b8 = 4 m0 = 915000 m1 = 1.2 m2 = -130 transition = 1 type = 4 t1 = 520 t2 = 5.7 end pnt EoS = 2 | H= -1359180. TiO2(1)MnO(1) G0 = -1280009 S0 = 104.9 V0 = 3.288 c1 = 141.9 c2 = .3373E-2 c3 = -1940700 c5 = -407.6 b1 = .495E-4 b5 = -.495E-3 b6 = 1770000. b7 = -265.5 b8 = 4 m0 = 915000 m1 = 1.2 m2 = -130 end geik EoS = 2 | H= -1567150. MgO(1)TiO2(1) G0 = -1478745 S0 = 74.60001 V0 = 3.086 c1 = 151 c3 = -1890400 c5 = -652.2 b1 = .495E-4 b5 = -.495E-3 b6 = 1770000. b7 = -265.5 b8 = 4 m0 = 1320000 m1 = 1.7 m2 = -170 end ilm EoS = 2 | H= -1231270. TiO2(1)FeO(1) G0 = -1154702 S0 = 108.9 V0 = 3.169 c1 = 138.9 c2 = .5081E-2 c3 = -1288800 c5 = -463.7 b1 = .495E-4 b5 = -.495E-3 b6 = 1770000. b7 = -265.5 b8 = 4 m0 = 910000 m1 = 1.7 m2 = -170 transition = 1 type = 4 t1 = 1900 t2 = 11 t3 = .2E-1 end ilm_nol EoS = 2 | created from ilm for Ilm(WPH) by suppressing the Landau model TiO2(1)FeO(1) G0 = -1154702 S0 = 108.9 V0 = 3.169 c1 = 138.9 c2 = .5081E-2 c3 = -1288800 c5 = -463.7 b1 = .495E-4 b5 = -.495E-3 b6 = 1770000. b7 = -265.5 b8 = 4 m0 = 910000 m1 = 1.7 m2 = -170 end bdy EoS = 2 | H= -1101310. ZrO2(1) G0 = -1043530 S0 = 50.4 V0 = 2.115 c1 = 90.7 c3 = -813300 c5 = -438.8 b1 = .376E-4 b5 = -.376E-3 b6 = 2224999. b7 = -333.75 b8 = 4 m0 = 910000 m1 = 1.7 m2 = -170 end sp EoS = 2 | H= -2300740. MgO(1)Al2O3(1) G0 = -2176030 S0 = 81.5 V0 = 3.978 c1 = 242.7 c2 = -.6037E-2 c3 = -2315100 c5 = -1678.1 b1 = .431E-4 b5 = -.431E-3 b6 = 1945000. b7 = -291.75 b8 = 4 m0 = 1080000 m1 = .5 m2 = -90 end | 0.64d6 0.2 0. 0. 0. 0. 0. 0. 0. 0. => chang & barsch '73, G'? no dG/dT herc EoS = 2 | H= -1959200. Al2O3(1)FeO(1) G0 = -1843852 S0 = 107.5 V0 = 4.075 c1 = 283.3 c2 = -.5376E-2 c3 = 609800 c5 = -2713.6 b1 = .395E-4 b5 = -.395E-3 b6 = 2120000. b7 = -318 b8 = 4 m0 = 840000 m1 = .4 m2 = -70 end | 0.64d6 0.2 0. 0. 0. 0. 0. 0. 0. 0. => sp value, chang & barsch '73, no dG/dT mt EoS = 2 | H= -1115490. FeO(3)O2(.5) G0 = -1012288 S0 = 146.1 V0 = 4.452 c1 = 262.5 c2 = -.7204E-2 c3 = -1926200 c5 = -1655.7 b1 = .696E-4 b5 = -.696E-3 b6 = 1850000. b7 = -277.5 b8 = 4 m0 = 915000 m1 = 1.2 m2 = -130 transition = 1 type = 4 t1 = 848 t2 = 35 end mft EoS = 2 | H= -1440650. MgO(1)FeO(2)O2(.5) G0 = -1329995 S0 = 126.5 V0 = 4.457 c1 = 217.9 c2 = .355E-3 c3 = -3108000 c5 = -745.9 b1 = .696E-4 b5 = -.696E-3 b6 = 1850000. b7 = -277.5 b8 = 4 m0 = 915000 m1 = 1.2 m2 = -130 transition = 1 type = 4 t1 = 665 t2 = 12.9 end usp EoS = 2 | H= -1497470. TiO2(1)FeO(2) G0 = -1401886 S0 = 175 V0 = 4.682 c1 = -102.6 c2 = .14252 c3 = -9144500 c5 = 5270.7 b1 = .69E-4 b5 = -.69E-3 b6 = 1850000. b7 = -277.5 b8 = 4 m0 = 915000 m1 = 1.2 m2 = -130 end br EoS = 2 | H= -924820. MgO(1)H2O(1) G0 = -834158.6 S0 = 64.5 V0 = 2.463 c1 = 158.4 c2 = -.4076E-2 c3 = -1052300 c5 = -1171.3 b1 = .13E-3 b5 = -.13E-2 b6 = 485000 b7 = -72.75 b8 = 4 m0 = 352000 m1 = 3.4 m2 = -105 end dsp EoS = 2 | H= -999450. Al2O3(.5)H2O(.5) G0 = -920768.2 S0 = 35 V0 = 1.776 c1 = 145.1 c2 = .8709E-2 c3 = 584400 c5 = -1741.1 b1 = .797E-4 b5 = -.797E-3 b6 = 2300000. b7 = -345 b8 = 4 m0 = 352000 m1 = 3.4 m2 = -105 end gth EoS = 2 | H= -561660. FeO(1)O2(.25)H2O(.5) G0 = -490870.2 S0 = 60.4 V0 = 2.082 c1 = 139.3 c2 = .147E-3 c3 = -212700 c5 = -1077.8 b1 = .797E-4 b5 = -.797E-3 b6 = 2300000. b7 = -345 b8 = 4 m0 = 352000 m1 = 3.4 m2 = -105 end cc EoS = 2 | H= -1207470. CaO(1)CO2(1) G0 = -1129155 S0 = 92.5 V0 = 3.689 c1 = 140.9 c2 = .5029E-2 c3 = -950700 c5 = -858.4 b1 = .44E-4 b5 = -.44E-3 b6 = 760000 b7 = -114 b8 = 4 m0 = 327000 m1 = .8829 m2 = -143.88 transition = 1 type = 4 t1 = 1240 t2 = 10 t3 = .4E-1 end CACO3L EoS = 2 CaO(1)CO2(1) G0 = -1047060 S0 = 130.505 V0 = 3.9858 c1 = 131.156 b1 = .256015E-3 b5 = -.256015E-2 b6 = 74728.65 b7 = -11.2093 b8 = 7 m0 = 327000 m1 = .8829 m2 = -143.88 end MGCO3L EoS = 2 MgO(1)CO2(1) G0 = -947706 S0 = 131.128 V0 = 3.11738 c1 = 99.9876 b1 = .256015E-3 b5 = -.256015E-2 b6 = 74728.65 b7 = -11.2093 b8 = 7 m0 = 680000 m1 = .8829 m2 = -143.88 end FECO3L EoS = 2 FeO(1)CO2(1) G0 = -615465.807755467 S0 = 141.295014032831 V0 = 3.51057462761850 c1 = 126.023 b1 = .256015E-3 b5 = -.256015E-2 b6 = 74728.65 b7 = -11.2093 b8 = 7 m0 = 680000 m1 = .8829 m2 = -143.88 end arag EoS = 2 | H= -1207580. CaO(1)CO2(1) G0 = -1128370 S0 = 89.5 V0 = 3.415 c1 = 192.3 c2 = -.3052E-2 c3 = 1149700 c5 = -2118.3 b1 = .115E-3 b5 = -.115E-2 b6 = 650000 b7 = -97.5 b8 = 4 m0 = 385000 m1 = .8829 m2 = -143.88 transition = 1 type = 4 t1 = 1240 t2 = 9 t3 = .4E-1 end mag EoS = 2 | H= -1111320. MgO(1)CO2(1) G0 = -1027504 S0 = 65.1 V0 = 2.803 c1 = 186.4 c2 = -.3772E-2 c5 = -1886.2 b1 = .648E-4 b5 = -.648E-3 b6 = 1460000. b7 = -219 b8 = 4 m0 = 680000 m1 = .8829 m2 = -143.88 end sid EoS = 2 | H= -761450. FeO(1)CO2(1) G0 = -688158.8 S0 = 95 V0 = 2.938 c1 = 168.4 c5 = -1483.6 b1 = .11E-3 b5 = -.11E-2 b6 = 1200000. b7 = -180 b8 = 4 m0 = 680000 m1 = .8829 m2 = -143.88 end rhc EoS = 2 | H= -891090. MnO(1)CO2(1) G0 = -817283 S0 = 98 V0 = 3.107 c1 = 169.5 c5 = -1534.3 b1 = .65E-4 b5 = -.65E-3 b6 = 800000 b7 = -120 b8 = 4 m0 = 680000 m1 = .8829 m2 = -143.88 end dol EoS = 2 | H= -2324460. MgO(1)CaO(1)CO2(2) G0 = -2161852 S0 = 156 V0 = 6.434 c1 = 358.9 c2 = -.4905E-2 c5 = -3456.2 b1 = .635E-4 b5 = -.6350001E-3 b6 = 900000 b7 = -135 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 transition = 1 type = 4 t1 = 1373 t2 = 13 t3 = .15E-1 end ank EoS = 2 | H= -1971300. CaO(1)FeO(1)CO2(2) G0 = -1819545 S0 = 187 V0 = 6.606 c1 = 341 c2 = -.1161E-2 c5 = -3054.8 b1 = .635E-4 b5 = -.6350001E-3 b6 = 900000 b7 = -135 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 transition = 1 type = 4 t1 = 1273 t2 = 9 t3 = .1E-1 end syv EoS = 2 | H= -436500. K2O(.5)Cl2(.5)O2(-.25) G0 = -408605.5 S0 = 82.6 V0 = 3.752 c1 = 46.2 c2 = .1797E-1 b1 = .247E-3 b5 = -.247E-2 b6 = 170000 b7 = -25.5 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 end hlt EoS = 2 | H= -411300. Na2O(.5)Cl2(.5)O2(-.25) G0 = -384264.2 S0 = 72.1 V0 = 2.702 c1 = 45.2 c2 = .1797E-1 b1 = .269E-3 b5 = -.269E-2 b6 = 240000 b7 = -36 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 end iron EoS = 2 | H= 0. FeO(1)O2(-.5) G0 = -11.9257 S0 = 27.32 V0 = .7092 c1 = 46.2 c2 = .5158E-2 c3 = 723100 c5 = -556.2 b1 = .746E-4 b5 = -.746E-3 b6 = 1680000. b7 = -252 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 transition = 1 type = 4 t1 = 1042 t2 = 8.3 end Ni EoS = 2 | H= 0. NiO(1)O2(-.5) G0 = 1.121429 S0 = 29.87 V0 = .6588 c1 = 49.8 c3 = 585900 c5 = -533.9 b1 = .886E-4 b5 = -.886E-3 b6 = 1870000. b7 = -280.5 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 transition = 1 type = 4 t1 = 631 t2 = 3 end gph EoS = 2 | H= 0. O2(-1)CO2(1) G0 = -32.79654 S0 = 5.85 V0 = .5298 c1 = 51 c2 = -.4428E-2 c3 = 488600 c5 = -805.5 b1 = .484E-4 b5 = -.484E-3 b6 = 390000 b7 = -58.5 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 end diam EoS = 2 | H= 2100. O2(-1)CO2(1) G0 = 3125.636 S0 = 2.3 V0 = .3417 c1 = 24.3 c2 = .6272E-2 c3 = -377400 c5 = -273.4 b1 = .165E-4 b5 = -.165E-3 b6 = 5800000. b7 = -870 b8 = 4 m0 = 457000 m1 = .8829 m2 = -143.88 end H2O EoS = 101 | H= -241810. H2O(1) G0 = -228542.3 S0 = 188.8 c1 = 40.1 c2 = .8656001E-2 c3 = 487500 c5 = -251.2 m0 = .1E-15 end CO2 EoS = 102 | H= -393510. CO2(1) G0 = -394332.9 S0 = 213.7 c1 = 87.8 c2 = -.2644E-2 c3 = 706400 c5 = -998.9 m0 = .1E-15 end CO EoS = 103 | H= -110530. O2(-.5)CO2(1) G0 = -137163.7 S0 = 197.67 c1 = 45.7 c2 = -.97E-4 c3 = 662700 c5 = -414.7 m0 = .1E-15 end CH4 EoS = 104 | H= -74810. O2(-2)H2O(2)CO2(1) G0 = -50695.63 S0 = 186.26 c1 = 150.1 c2 = .2062E-2 c3 = 3427700 c5 = -2650.4 m0 = .1E-15 end O2 EoS = 107 | H= 0. O2(1) S0 = 205.2 c1 = 48.3 c2 = -.691E-3 c3 = 499200 c5 = -420.7 m0 = .1E-15 end H2 EoS = 105 | H= 0. O2(-.5)H2O(1) G0 = -.4549408E-2 S0 = 130.7 c1 = 23.3 c2 = .4627E-2 c5 = 76.3 m0 = .1E-15 end syvL EoS = 2 | H= -416520. K2O(.5)Cl2(.5)O2(-.25) G0 = -392412 S0 = 95.3 V0 = 3.857 c1 = 66.9 b1 = .5E-3 b5 = -.5E-2 b6 = 59000 b7 = -8.85 b8 = 4 m0 = .1E-15 end hltL EoS = 2 | H= -393450. Na2O(.5)Cl2(.5)O2(-.25) G0 = -368680.2 S0 = 79.7 V0 = 2.965 c1 = 72 c2 = -.3223E-2 b1 = .5E-3 b5 = -.5E-2 b6 = 66000.01 b7 = -9.9 b8 = 4 m0 = .1E-15 end foL EoS = 2 | H= -2225200. MgO(2)SiO2(1) G0 = -2061346 S0 = -55 V0 = 4.243 c1 = 267.9 b1 = .145E-3 b5 = -.145E-2 b6 = 730000 b7 = -109.5 b8 = 4 m0 = .1E-15 end faL EoS = 2 | H= -1458790. SiO2(1)FeO(2) G0 = -1345114 S0 = 102.5 V0 = 4.695 c1 = 239.7 b1 = .169E-3 b5 = -.169E-2 b6 = 410000 b7 = -61.5 b8 = 4 m0 = .1E-15 end silL EoS = 2 | H= -2583280. Al2O3(1)SiO2(1) G0 = -2419443 S0 = 39 V0 = 6.419 c1 = 237.6 b1 = .1E-4 b5 = -.1E-3 b6 = 300000.1 b7 = -45 b8 = 4 m0 = .1E-15 end anL EoS = 2 | H= -4257680. Al2O3(1)SiO2(2)CaO(1) G0 = -3987929 S0 = 52 V0 = 10.206 c1 = 417.5 b1 = .49E-4 b5 = -.49E-3 b6 = 200000 b7 = -30 b8 = 4 m0 = .1E-15 end h2oL EoS = 2 | H= -296060. H2O(1) G0 = -240067.4 S0 = 45.5 V0 = 1.414 c1 = 80 b1 = .1079E-2 b5 = -.1079E-1 b6 = 40000 b7 = -6 b8 = 4 m0 = .1E-15 end enL EoS = 2 | H= -3091500. MgO(2)SiO2(2) G0 = -2876659 S0 = -2 V0 = 6.9 c1 = 354.9 b1 = .129E-3 b5 = -.129E-2 b6 = 260000 b7 = -39 b8 = 4 m0 = .1E-15 end kspL EoS = 2 | H= -3992580. Al2O3(.5)SiO2(3)K2O(.5) G0 = -3741907 S0 = 129.5 V0 = 11.468 c1 = 367.3 b1 = .6E-4 b5 = -.6E-3 b6 = 260000 b7 = -39 b8 = 4 m0 = .1E-15 end qL EoS = 2 | H= -920860. SiO2(1) G0 = -858990.9 S0 = 16.5 V0 = 2.64 c1 = 82.5 b1 = -.5E-5 b5 = .5E-4 b6 = 470000 b7 = -70.5 b8 = 4 m0 = .1E-15 end diL EoS = 2 | H= -3207640. MgO(1)SiO2(2)CaO(1) G0 = -2997823 S0 = 23.8 V0 = 7.34 c1 = 345.3 b1 = .129E-3 b5 = -.129E-2 b6 = 300000.1 b7 = -45 b8 = 4 m0 = .1E-15 end abL EoS = 2 | H= -3934330. Na2O(.5)Al2O3(.5)SiO2(3) G0 = -3692268 S0 = 145 V0 = 10.71 c1 = 358.5 b1 = .45E-4 b5 = -.45E-3 b6 = 390000 b7 = -58.5 b8 = 4 m0 = .1E-15 end kal3o EoS = 1 Al2O3(1.5)K2O(.5) G0 = 100 V0 = .1E-3 m0 = .1E-15 end kalo2 EoS = 1 Al2O3(.5)K2O(.5) G0 = 100 V0 = .1E-3 m0 = .1E-15 end k2o EoS = 1 K2O(1) G0 = 100 V0 = .1E-3 m0 = .1E-15 end Na2O EoS = 1 Na2O(1) G0 = 100 V0 = .1E-3 m0 = .1E-15 end qGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -3.555755e6 | modified 5/16/2012 from -3383116. to match 1 bar liquidus at 1999 K SiO2(4) G0 = -3363173.2 S0 = 107.93 V0 = 10.76 c1 = 330.4 b2 = -.756E-4 b3 = .52E-7 b8 = -5 m0 = .1E-15 end coGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -3.307937e6 | modified 5/16/2012 from -3034014. to match 1 bar liquidus at 2320 K Al2O3(2) G0 = -3135674.4 S0 = -44.57 V0 = 7.422 c1 = 340.6 b1 = .524E-3 b2 = -.452E-4 b3 = .54E-7 b8 = -5 m0 = .1E-15 end faGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -1.465621e6 | modified 5/16/2012 from -1336321. to match fa 1 bar liquidus at 1490 K SiO2(1)FeO(2) G0 = -1339255. S0 = 106.41 V0 = 5.42 c1 = 240.2 b1 = .584E-3 b2 = -.279E-4 b3 = -.23E-7 b8 = -5 m0 = .1E-15 end foGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -2.166508e6 | modified 5/16/2012 from -2059270. to match fo 1 bar liquidus at 2163 K MgO(2)SiO2(1) G0 = -2039445. S0 = -47.34 V0 = 4.98 c1 = 271 b1 = .524E-3 b2 = -.135E-4 b3 = -.13E-7 b8 = -5 m0 = .1E-15 end woGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -3.276717e6 | modified 5/16/2012 from -3108527. to match 1 bar liquidus at 1817 K SiO2(2)CaO(2) G0 = -3117049.6 S0 = 20.31 V0 = 8.694 c1 = 344.8 b1 = .584E-3 b2 = -.31E-4 b3 = -.32E-7 b8 = -5 m0 = .1E-15 end nasGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -8.00742e5 Na2O(.5)SiO2(.5) G0 = -739767.9 S0 = 50.01 V0 = 2.784 c1 = 90.1 b1 = .371E-3 b2 = -.215E-4 b3 = -.265E-7 b8 = -5 m0 = .1E-15 end kalGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) -2.127829e6 | modified 5/16/2012 from -2014274. to match fo 1 bar liquidus at 2023 K Al2O3(.5)SiO2(1)K2O(.5) G0 = -2019205.3 S0 = 49.59 V0 = 6.8375 c1 = 217 b1 = .7265E-3 b2 = -.64E-4 b3 = -.46E-7 b8 = -5 m0 = .1E-15 end h2oGL EoS = 3 | PMELT liquid endmember (Ghiorso et al, G3, 2002) 6.64e2 H2O(1) G0 = 67460.3 S0 = 9.22 V0 = 2.775 b1 = .1086E-2 b2 = -.6E-4 b8 = -3.5 m0 = .1E-15 end Si EoS = 2 | Silicon metal, not from HP data base, JADC 1/2010 SiO2(1)O2(-1) S0 = 18.82 V0 = 1.21 c1 = 36.12 c2 = -.398E-3 c3 = .8987E-2 c5 = -285.3 b1 = .9E-5 b5 = -.9E-4 b6 = 933470.1 b7 = -146 b8 = 4 m0 = .1E-15 end SiC EoS = 2 | Silicon Carbide, not from HP data base, JADC 1/2010 SiO2(1)O2(-2)CO2(1) G0 = -69139 S0 = 16.485 V0 = 2.1652 c1 = 64.21 c2 = .3308E-2 c3 = .1991E-1 c5 = -631.7 b1 = .831E-5 b5 = -.831E-4 b6 = 2139822. b7 = -336 b8 = 4 m0 = .1E-15 end fpum EoS = 2 | ferro-pumpellyite from Massonne & Willner, EJM, 2008 Al2O3(2.5)SiO2(6)CaO(4)FeO(1)H2O(3.5) G0 = -13099500 S0 = 654 V0 = 29.72 c1 = 1568.5 c3 = -26900000 c5 = -10207.2 c7 = .3560562E10 b1 = .5E-4 b5 = -.5E-3 b6 = 1615000. b7 = -242.25 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end ffpu EoS = 2 | ferro-pumpellyite from Massonne & Willner, EJM, 2008 Al2O3(2)SiO2(6)CaO(4)FeO(2)O2(.25)H2O(3.5) G0 = -12672500 S0 = 687 V0 = 30.1 c1 = 1574.92 c3 = -24553960 c5 = -10268.96 c7 = .3154003E10 b1 = .5E-4 b5 = -.5E-3 b6 = 1615000. b7 = -242.25 b8 = 4 m0 = 454000 m1 = .97133 m2 = -68 end acti EoS = 2 | actinolite from Massonne & Willner, EJM, 2008 MgO(3)SiO2(8)CaO(2)FeO(2)H2O(1) G0 = -10903000 S0 = 614 V0 = 27.67 c1 = 1258.028 c3 = -24442220 c5 = -7489.36 c7 = .3162433E10 b1 = .534E-4 b5 = -.534E-3 b6 = 761253.7 b7 = -114 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end mrie EoS = 2 | magnesio-riebeckite from Massonne & Willner, EJM, 2008 Na2O(1)MgO(3)SiO2(8)FeO(2)O2(.5)H2O(1) G0 = -10370000 S0 = 602 V0 = 27.13 c1 = 1246.525 c3 = -29546270 c5 = -6823.08 c7 = .3998917E10 b1 = .534E-4 b5 = -.534E-3 b6 = 761253.7 b7 = -114 b8 = 4 m0 = 617000 m1 = 1.19757 m2 = -96 end stlp EoS = 2 | stilpnomelane from Massonne & Willner, EJM, 2008 Al2O3(2.5)SiO2(67)K2O(2.5)FeO(48)H2O(60) G0 = -91150000 S0 = 7364 V0 = 303.7 c1 = 15936.51 c3 = -181141800 c5 = -124844 c7 = .2860103E11 b1 = .596E-4 b5 = -.596E-3 b6 = 490000 b7 = -73.5 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end mstl EoS = 2 | magnesio-stilpnomelane from Massonne & Willner, EJM, 2008 MgO(48)Al2O3(2.5)SiO2(67)K2O(2.5)H2O(60) G0 = -107600000 S0 = 5911 V0 = 296.5 c1 = 15032.19 c3 = -248561600 c5 = -104513 c7 = .3374243E11 b1 = .596E-4 b5 = -.596E-3 b6 = 490000 b7 = -73.5 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end mnsp EoS = 2 | manganese-stilpnomelane from Massonne & Willner, EJM, 2008 Al2O3(2.5)SiO2(67)K2O(2.5)MnO(48)H2O(60) G0 = -97200000 S0 = 7642 V0 = 313.3 c1 = 12238.78 c3 = -181141800 c5 = -96778 c7 = .2833413E11 b1 = .596E-4 b5 = -.596E-3 b6 = 490000 b7 = -73.5 b8 = 4 m0 = 353000 m1 = 1.21 m2 = -54 end kcym EoS = 2 | potassium whatever from Massonne & Willner, EJM, 2008 Al2O3(.5)SiO2(3)K2O(.5)H2O(1) G0 = -3966800 S0 = 278.6 V0 = 11.459 c1 = 481.225 c3 = -9931330 c5 = -2980.73 c7 = 1416000000 b1 = .335E-4 b5 = -.335E-3 b6 = 574000 b7 = -86.1 b8 = 4 m0 = 301000 m1 = 4.26 m2 = -70 end |Fe-BCC EoS = 600 | Komabayashi & Fei JGR 2010 Special EoS |FeO(1)O2(-.5) |0 = -6.5 V0 = .709 |c1 = 1400 c2 = 124.06 c3 = -25.5143 c4 = -.439752E-2 c5 = -.58927E-7 c6 = 77359 |b1 = .262E-4 b2 = .18E-7 b6 = 1640000 b8 = -5.29 |m0 = 454000 m1 = .97133 m2 = -68 |end |Fe-FCC EoS = 601 | Komabayashi & Fei JGR 2010 Special EoS |FeO(1)O2(-.5) |S0 = -6.5 V0 = .6835 |c1 = 16300.92 c2 = 381.4716 c3 = -52.2754 c4 = .177578E-3 c6 = -395355.4 |b1 = .64E-4 b6 = 1653000 b8 = -5.5 |m0 = 454000 m1 = .97133 m2 = -68 |end |Fe-HCP EoS = 602 | Komabayashi & Fei JGR 2010 Special EoS (dummy alpha & delta) |FeO(1)O2(-.5) |S0 = -4.97 V0 = .6765 |c1 = 14405.92 c2 = 384.8716 c3 = -52.2754 c4 = .177578E-3 c6 = -395355.4 |b1 = .64E-4 b6 = 16500000 b8 = -4.97 |m0 = 454000 m1 = .97133 m2 = -68 |end |Fe-Liq EoS = 603 | Komabayashi & Fei JGR 2010 Special EoS |FeO(1)O2(-.5) |S0 = -5.22 V0 = .6921 |c1 = -9007.34 c2 = 290.2987 c3 = -46 |b1 = .92E-4 b6 = 12400000 b8 = -5.77 |m0 = 454000 m1 = .97133 m2 = -68 |end skiag EoS = 2 | Malaspina et al., 2009. Skiagite garnet SiO2(3)FeO(5)O2(.5) G0 = -4072.02e3 S0 = 369.871 V0 = 12.135 c1 = 870.4 c2 = -.14476 c3 = -400450 c4 = .43207E-4 c5 = -8652.6 b1 = .2336E-4 b2 = .70286E-8 b4 = -.295 b6 = 1574000 b8 = 6.7 end tiGL EoS = 3 | PMELT liquid endmember, G/S adjusted to match 1840 K 1 bar fusion for H&P rutile TiO2(1) G0 = -867802 S0 = 13.0582 V0 = 2.316 c1 = 109.2 b1 = .7246E-3 b2 = -.231E-4 b8 = -5 m0 = .1E-15 end sud_dqf EoS = 2 | H= -8647122. MgO(2)Al2O3(2)SiO2(3)H2O(4) G0 = -7995205 S0 = 417.39 V0 = 20.3 c1 = 1104.9 c2 = .11217E-1 c3 = -6811700 c5 = -10110.5 b1 = .398E-4 b5 = -.398E-3 b6 = 870000 b7 = -130.5 b8 = 4 m0 = 405000 m1 = 1.20881 m2 = -54 end | EoS 610-633 yields HSC convention apparent Gibbs energies | it is corrected here to GSC by adding G0 = 298*sum(S_elements) Fe-BCC EoS = 610 FeO(1)O2(-.5) G0 = 8134. V0 = .705 b1 = .27919E-4 b6 = 1700000 b8 = 6.2 end Si-BCC EoS = 611 SiO2(1)O2(-1) G0 = 5608. V0 = .93 b1 = .608E-5 b6 = 710000 b8 = 4 end Fe-FCC EoS = 612 FeO(1)O2(-.5) G0 = 8134. V0 = .6826 b1 = .5797E-4 b5 = -.5797E-3 b6 = 1400000 b8 = 8 end Si-FCC EoS = 613 SiO2(1)O2(-1) G0 = 5608. V0 = .896 b1 = .608E-5 b6 = 707000 b8 = 4 end Fe_LIQ EoS = 614 FeO(1)O2(-.5) G0 = 8134. V0 = .706 b1 = .512E-4 b5 = -.512E-3 b6 = 1400000 b8 = 5.8 end Si_LIQ EoS = 615 SiO2(1)O2(-1) G0 = 5608. V0 = .951 b1 = .12067E-3 b5 = -.12067E-2 b6 = 660000 b8 = 4 end |fe2si EoS = 616 |SiO2(.3333333)FeO(.6666667)O2(-.6666667) |G0 = 7292. |end |fe5si3 EoS = 617 |SiO2(.375)FeO(.625)O2(-.6875) |G0 = 7187. |end fesi EoS = 618 SiO2(.5)FeO(.5)O2(-.75) G0 = 6871 V0 = .6745 b1 = .485E-4 b6 = 1870000 b8 = 5.2 end |fesi2 EoS = 619 |SiO2(.6666667)FeO(.3333333)O2(-.8333333) |G0 = 6450 |end |fe3si7 EoS = 620 |SiO2(.7)FeO(.3)O2(-.85) |G0 = 6366. |end Si-diam EoS = 621 SiO2(1)O2(-1) G0 = 5608. V0 = 1.2 b1 = .608069E-5 b5 = -.608069E-4 b6 = 1000000 b8 = 4 end FeC-BCC EoS = 622 FeO(1)CO2(3)O2(-3.5) G0 = 13268. V0 = .72 b6 = 1700000 b8 = 6.2 end SiC-BCC EoS = 623 SiO2(1)CO2(3)O2(-4) G0 = 10742. V0 = .93 b6 = 710000 b8 = 4 end FeC-FCC EoS = 624 FeO(1)CO2(1)O2(-1.5) G0 = 13268. V0 = .726 b6 = 1400000 b8 = 8 end SiC-FCC EoS = 625 SiO2(1)CO2(1)O2(-2) G0 = 10742. V0 = .896 b6 = 707000 b8 = 4 end C_LIQ EoS = 626 CO2(1)O2(-1) G0 = 1711. V0 = .778 b1 = .544E-3 b5 = -.544E-2 b6 = 100000 b8 = 4 end GPH EoS = 627 CO2(1)O2(-1) G0 = 1711. V0 = .5298 b1 = .484E-4 b5 = -.484E-3 b6 = 390000 b8 = 4 end SiC-lac EoS = 628 SiO2(1)CO2(1)O2(-2) G0 = 10742. V0 = 1.247 b1 = .104105E-4 b5 = -.104105E-3 b6 = 2200000 b8 = 4 end Fe-CBCC EoS = 629 FeO(1)O2(-.5) G0 = 8134. V0 = .6677 b1 = .5678E-4 b6 = 1700000 b8 = 5.5 end Si-CBCC EoS = 630 SiO2(1)O2(-1) G0 = 5608. V0 = .9 b1 = .608E-5 b6 = 700000 b8 = 4 end |FeC-CBCC EoS = 631 |FeO(1)CO2(1)O2(-1.5) |G0 = 13268. |end |SiC-CBCC EoS = 632 |SiO2(1)CO2(1)O2(-2) |G0 = 10742. |end CEM EoS = 633 FeO(3)CO2(1)O2(-2.5) G0 = 7091. V0 = 2.323 b1 = .97E-4 b5 = -.97E-3 b6 = 1740000 b8 = 4.8 end Fe-HCP EoS = 634 FeO(1)O2(-.5) G0 = 8134. V0 = .6677 b1 = .5678E-4 b6 = 1700000 b8 = 5.5 end Si-HCP EoS = 635 SiO2(1)O2(-1) G0 = 5608. V0 = .9 b1 = .608E-5 b6 = 700000 b8 = 4 end FeC-HCP EoS = 636 FeO(1)CO2(.5)O2(-1) G0 = 13268. V0 = .6677 b1 = .5678E-4 b5 = -.5678E-3 b6 = 1700000 b8 = 5.5 end SiC-HCP EoS = 637 SiO2(1)CO2(.5)O2(-1.5) G0 = 10742. V0 = .9 b1 = .608E-5 b5 = -.608E-4 b6 = 700000 b8 = 4 end |FeSiC EoS = 638 |SiO2(.1818)FeO(.7273)CO2(.0909)O2(-.63635) |G0 = 7091. |end Cdiam EoS = 639 CO2(1)O2(-1) G0 = 1711. V0 = .342 b1 = .54E-5 b5 = -.54E-4 b6 = 4420000 b8 = 4 end FeC4H6O4 EoS = 16 | Gf = -844331. Fe(C2H6C2O2O2)(0) Sverjensky et al. (1997) with revised volume increased in order that a1 of the complex is the sum of the a1 valu FeO(1)H2O(3)CO2(4)O2(-4) G0 = -844331.2 S0 = 48.24152 w = -15899.2 a1 = 4.58056 a2 = 7929.308 a3 = -7.119494 a4 = -149050.8 c1 = 476.0103 c2 = 1442932. end FeC2H3O2 EoS = 16 | Gf = -468190. Fe(CH3COO)(+) Sverjensky et al. (1997) with revised volume increased in order that a1 of the complex is the sum of the a1 values of the ions FeO(1)H2O(1.5)CO2(2)O2(-2.25) G0 = -468189.6 S0 = -6.276 w = 240831 q = 1 a1 = 1.644898 a2 = 759.8562 a3 = 21.08652 a4 = -119411.4 c1 = 247.097 c2 = 565082.7 end NaCO3- EoS = 16 | Gf = -797311. NaCO3(-) "Fit to Stefansson et al. (2013; 80 to 200 C) & Garrels et al., (1961; at 25 C) to get G, S and Cp; with a1 pred. from the sum of the ions and Na2O(.5)CO2(1)O2(.25) G0 = -797311.4 S0 = 33.472 w = 541577 q = -1 a1 = 3.197706 a2 = 1252.48 a3 = 9.74077 a4 = -121453.2 c1 = 81.47503 c2 = -135490.5 end CO3-2 EoS = 16 | Gf = -527983. CO3(-2) Facq et al. (2014) CO2(1)O2(.5) G0 = -527983.1 S0 = -54.1828 w = 1924640. q = -2 a1 = 2.42672 a2 = 2092 a3 = -8.368 a4 = -451872 c1 = 75.312 c2 = -836800 end CO2,aq EoS = 16 | Gf = -385974. CO2(0) Facq et al. (2014) CO2(1) G0 = -385974 S0 = 117.5704 w = -83680 a1 = 3.259336 a2 = 1313.776 a3 = 9.24664 a4 = -121754.4 c1 = 154.808 c2 = 271960 end BENZENE EoS = 16 | Gf = 133888. C6H6(0) "Plyasunov & Shock (2001) for c1, c2 & omega; a1 predicted with new correlation for complex species from Sverjensky et al. (2014)" H2O(3)CO2(6)O2(-7.5) G0 = 133888 S0 = 149.3688 w = -104600 a1 = 7.330326 a2 = 5187.156 a3 = -23.25091 a4 = -137716.4 c1 = 279.9096 c2 = 410032 end ACETATE EoS = 16 | Gf = -369322. CH3COO- Shock & Helgeson (1990); Shock (1995) H2O(1.5)CO2(2)O2(-1.75) G0 = -369321.7 S0 = 86.1904 w = 551534.9 q = -1 a1 = 3.243646 a2 = 3639.913 a3 = 31.72518 a4 = -131314.8 c1 = 110.0392 c2 = -161502.4 end ACETIC-A EoS = 16 | Gf = -396476. CH3COOH Plyasunov & Shock (2001) H2O(1.5)CO2(2)O2(-1.75) G0 = -396475.8 S0 = 178.6568 w = -31798.4 a1 = 4.8116 a2 = 2301.2 a3 = 6.98728 a4 = -120080.8 c1 = 187.8616 c2 = -110039.2 end Al+3 EoS = 16 | Gf = -483708. Al(+3) Shock et al. (1997) Al2O3(.5)O2(-.75) G0 = -483708.1 S0 = -325.0968 w = 1151855. q = 3 a1 = -1.414276 a2 = -7115.771 a3 = 60.7454 a4 = -86851.47 c1 = 44.7688 c2 = -337230.4 end AlO2- EoS = 16 | Gf = -831332. AlO2(-) Sverjensky et al. (2014) Al2O3(.5)O2(.25) G0 = -831331.5 S0 = -30.20848 w = 728769.1 q = -1 a1 = 1.798283 a2 = 1044.996 a3 = -7.639147 a4 = -120595.4 c1 = 79.9144 c2 = -259408 end Ca(HCO3) EoS = 16 | Gf = -1145705. Ca(HCO3)(+) Facq et al. (2014) CaO(1)H2O(.5)CO2(1)O2(-.25) G0 = -1145705. S0 = 66.944 w = 292880 q = 1 a1 = 3.824176 a2 = 1849.328 a3 = 4.72792 a4 = -123846.4 c1 = 275.7256 c2 = 648520 end Ca(OH)+ EoS = 16 | Gf = -716719. Ca(OH)(+) Shock et al. (1997) with new a1 prediction CaO(1)H2O(.5)O2(-.25) G0 = -716719.2 S0 = 28.0328 w = 188112.6 q = 1 a1 = -.3369375 a2 = -4077.224 a3 = 40.05427 a4 = -99416.02 c1 = 34.9774 c2 = -161054.7 end Ca+2 EoS = 16 | Gf = -552790. Ca(+2) Facq et al. (2014) CaO(1)O2(-.5) G0 = -552790.1 S0 = -56.484 w = 517393.4 q = 2 a1 = -.102508 a2 = -3034.237 a3 = 20.92 a4 = -103763.2 c1 = 37.656 c2 = -105436.8 end CaCl+ EoS = 16 | Gf = -682410. CaCl(+) Sverjensky et al. (1997) with new V from the sum of a1 values of ions CaO(1)Cl2(.5)O2(-.5) G0 = -682410.4 S0 = 18.828 w = 195602 q = 1 a1 = 1.198298 a2 = -328.4022 a3 = 25.31948 a4 = -114913.6 c1 = 89.47693 c2 = 17238.08 end CaCl2,aq EoS = 16 | Gf = -811696. CaCl(0) Sverjensky et al. (1997) with new V from the sum of a1 values of ions CaO(1)Cl2(.5)O2(-.5) G0 = -811696 S0 = 25.104 w = -15899.2 a1 = 2.685459 a2 = 3302.264 a3 = 11.05957 a4 = -129921.6 c1 = 97.33574 c2 = 126754.3 end CaCO3,aq EoS = 16 | Gf = -1099346. CaCO3(0) Facq et al. (2014) CaO(1)CO2(1) G0 = -1099346 w = 836800 a1 = 2.786544 a2 = 857.72 a3 = 13.05408 a4 = -119662.4 c1 = 137.2352 c2 = -7531.2 end CaSO4,aq EoS = 16 | Gf = -1309299. CaSO4(0) Sverjensky et al. (1997) revised using exptl V and pred. a1 from new correlation CaO(1)S2(.5)O2(1.5) G0 = -1309299. S0 = 20.92 w = -15899.2 a1 = 1.162524 a2 = -417.4377 a3 = 25.70148 a4 = -114545.4 c1 = -36.96564 c2 = -340033.7 end Cl- EoS = 16 | Gf = -131290. Cl(-) Shock & Helgeson (1988) Cl2(.5) G0 = -131289.7 S0 = 56.73504 w = 609190.4 q = -1 a1 = 1.686989 a2 = 2008.738 a3 = 23.27559 a4 = -119118.5 c1 = -18.4096 c2 = -239073.8 end CO,aq EoS = 16 | Gf = -120005. CO(0) "Shock & McKibben (1993) for c1, c2 and omega with revised a1 from correlation for complex species " CO2(1)O2(-.5) G0 = -120005.5 S0 = 102.6335 w = -155435.6 a1 = 2.609686 a2 = 3116.996 a3 = 11.79219 a4 = -129155.9 c1 = 168.4156 c2 = 418492 end Cr+2 EoS = 16 | Gf = -164850. Cr(+2) Shock et al. (1997) revised with new predicted a1 for cations Cr2O3(.5)O2(-.75) G0 = -164849.6 S0 = -101.2528 w = 597768.1 q = 2 a1 = -.3362262 a2 = -4075.216 a3 = 40.03586 a4 = -99420.21 c1 = 62.97004 c2 = -189179.6 end Cr+3 EoS = 16 | Gf = -206271. Cr(+3) Shock et al. (1997) revised with new predicted a1 for cations Cr2O3(.5)O2(-.75) G0 = -206271.2 S0 = -322.168 w = 1146542. q = 3 a1 = -2.972439 a2 = -1395.155 a3 = 55.50285 a4 = -110507.8 c1 = 75.27602 c2 = -322134.5 end Cr2O7-2 EoS = 16 | Gf = -1301224. Cr2O7(-2) Shock et al. (1997) revised with new predicted a1 for cations Cr2O3(1)O2(2) G0 = -1301224 S0 = 261.9184 w = 947843.4 q = -2 a1 = 5.238368 a2 = 9534.876 a3 = -13.41809 a4 = -155690.8 c1 = 63.10309 c2 = -300833.8 end CrO4-2 EoS = 16 | Gf = -731363. CrO4(-2) Shock et al. (1997) revised with new predicted a1 for cations Cr2O3(.5)O2(1.25) G0 = -731363.2 S0 = 57.7392 w = 1256204. q = -2 a1 = 2.296639 a2 = 2352.37 a3 = 14.80425 a4 = -125992.8 c1 = -11.96959 c2 = -660490.4 end Cu+ EoS = 16 | Gf = 49999. Cu(+) Shock et al. (1997) revised with new predicted a1 for cations CuO(1)O2(-.5) G0 = 49998.8 S0 = 40.5848 w = 140205.8 q = 1 a1 = .3278164 a2 = -2455.255 a3 = 33.7084 a4 = -106123 c1 = 72.31249 c2 = -10204.78 end Cu+2 EoS = 16 | Gf = 65584. Cu(+2) Shock et al. (1997) revised with new predicted a1 for cations CuO(1)O2(-.5) G0 = 65584.2 S0 = -97.0688 w = 617935 q = 2 a1 = -.4611186 a2 = -4381.736 a3 = 41.28018 a4 = -98160.82 c1 = 84.9352 c2 = -183677.6 end ETHANE,A EoS = 16 | Gf = -16259. C2H6 Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(3)CO2(2)O2(-3.5) G0 = -16259.02 S0 = 112.173 w = -169870.4 a1 = 3.612466 a2 = 5565.18 a3 = 2.177772 a4 = -139277 c1 = 226.6703 c2 = 625595.9 end ETHANOL, EoS = 16 | Gf = -181293. C2H5OH Plyasunov & Shock (2001) H2O(2.5)CO2(2)O2(-2.75) G0 = -181292.7 S0 = 150.2056 w = -85228.08 a1 = 3.863213 a2 = 4166.469 a3 = 50.81259 a4 = -133494.7 c1 = 251.1132 c2 = 6305.288 end ETHYLENE EoS = 16 | Gf = 81379. C2H4 "Plyasunov & Shock (2001) for c1, c2 & omega; a1 predicted with new correlation for complex species from Sverjensky et al. (2014)" H2O(2)CO2(2)O2(-3) G0 = 81378.8 S0 = 119.6624 w = -129704 a1 = 4.311779 a2 = 2313.166 a3 = .8468416 a4 = -125838 c1 = 184.096 c2 = 539736 end Fe+2 EoS = 16 | Gf = -91504. Fe(+2) Shock et al. (1997) revised with new predicted a1 for cations FeO(1)O2(-.5) G0 = -91504.08 S0 = -105.8552 w = 609776.2 q = 2 a1 = -.3264775 a2 = -4052.915 a3 = 39.98774 a4 = -99520.62 c1 = 62.60603 c2 = -194296.6 end Fe+3 EoS = 16 | Gf = -17238. Fe(+3) Shock et al. (1997) revised with new predicted a1 for cations FeO(1)O2(-.5) G0 = -17238.08 S0 = -277.3992 w = 1079974. q = 3 a1 = -1.014871 a2 = -5730.448 a3 = 46.50139 a4 = -92579.37 c1 = 79.68805 c2 = -285486.9 end FeCl+ EoS = 16 | Gf = -221878. FeCl(+) Sverjensky et al. (1997) with revised volume increased in order that a1 of the complex is the sum of the a1 values of the ions FeO(1)Cl2(.5)O2(-.5) G0 = -221877.5 S0 = -42.09104 w = 293005.5 q = 1 a1 = .868431 a2 = -1135.077 a3 = 28.51145 a4 = -111574.7 c1 = 103.2348 c2 = 48346.12 end |FeCl+ FeCl(+) -53030 -61260 -10.060 5.8 20.6 0.7003 1 2.0756 -2.7129 6.8144 -2.6667 24.6737 1.1555 Sverjensky et al. (1997) with revised volume increased in order that a1 of the complex is the sum of the a1 values of the ions |FeCl++ FeCl(++) -37573 -50875 -28.700 -22.8 3.1 1.4967 2 -0.7191 -9.5343 9.4904 -2.3848 21.7718 -2.4031 Sverjensky et al. (1997) with revised volume increased in order that a1 of the complex is the sum of the a1 values of the ions FeCl2,aq EoS = 16 | Gf = -307440. FeCl2(0) Sverjensky et al. (1997) with revised volume increased in order that a1 of the complex is the sum of the a1 values of the ions FeO(1)Cl2(1)O2(-.5) G0 = -307440.3 S0 = 179.912 w = -15899.2 a1 = 2.391156 a2 = 2582.532 a3 = 13.91515 a4 = -126946.7 c1 = 93.33835 c2 = 112863.4 end FORMATE, EoS = 16 | Gf = -350879. HCOO- Shock & Helgeson (1990) revised with new predicted a1 for ions H2O(.5)CO2(1)O2(-.25) G0 = -350878.6 S0 = 90.7928 w = 544045.5 q = -1 a1 = 2.204299 a2 = 3214.818 a3 = 14.52099 a4 = -129561.7 c1 = 71.128 c2 = -518816 end FORMIC-A EoS = 16 | Gf = -372301. HCOOH Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(.5)CO2(1)O2(-.25) G0 = -372300.7 S0 = 162.7576 w = -144013.3 a1 = 3.385902 a2 = 1431.639 a3 = 8.238714 a4 = -122193.7 c1 = 92.853 c2 = -130524.1 end GLUTARAT EoS = 16 | Gf = -714878. C5H7O4- Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(3.5)CO2(5)O2(-6.25) G0 = -714878.2 S0 = 221.3336 w = 346393.4 q = -1 a1 = 5.77392 a2 = 10271.72 a3 = -4.10032 a4 = -158573.6 c1 = 148.8667 c2 = -161502.4 end GLUTARIC EoS = 16 | Gf = -739648. C5H8O4 Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(4)CO2(5)O2(-5) G0 = -739647.5 S0 = 306.2688 w = -41840 a1 = 6.40152 a2 = 11673.36 a3 = -7.28016 a4 = -164431.2 c1 = 252.7136 c2 = 57320.8 end GLYCOLAT EoS = 16 | Gf = -506975. C2H3O3- Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(1.5)CO2(2)O2(-2.25) G0 = -506975.3 S0 = 109.6208 w = 516054.6 q = -1 a1 = 3.194442 a2 = 4626.249 a3 = 4.114546 a4 = -135394.2 c1 = 108.9777 c2 = -168493.9 end GLYCOLIC EoS = 16 | Gf = -528858. C2H4O3 Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(2)CO2(2)O2(-1.5) G0 = -528857.6 S0 = 180.3304 w = -126231.3 a1 = 3.658699 a2 = 5083.393 a3 = 16.82888 a4 = -137285.4 c1 = 182.7023 c2 = -19476.52 end H+ EoS = 16 | Gf = 0. H(+) "Perfectly known, by convention!" H2O(.5)O2(-.25) q = 1 HOH = 1 end H2,aq EoS = 16 | Gf = 17723. H2(0) Shock et al. (1989) H2O(1)O2(-.5) G0 = 17723.42 S0 = 57.7392 w = -87445.6 a1 = 2.151706 a2 = 1998.195 a3 = 16.20421 a4 = -124532.6 c1 = 115.5834 c2 = 213091.1 end H2S,aq EoS = 16 | Gf = -27920. H2S(0) Shock et al. (1989) with revised prediction of a1 using new correlation H2O(1)S2(.5)O2(-.5) G0 = -27919.83 S0 = 125.52 w = -41840 a1 = 2.723658 a2 = 2833.572 a3 = 24.95589 a4 = -127988.6 c1 = 135.1432 c2 = 197903.2 end HAlO2,aq EoS = 16 | Gf = -868180. HAlO2(0) Sverjensky et al. (2014) Al2O3(.5)H2O(.5) G0 = -868180 S0 = -27.196 w = 209200 a1 = 1.792844 a2 = -85.10256 a3 = 20.9564 a4 = -115921.9 c1 = 161.548 c2 = 277867.8 end HCl,aq EoS = 16 | Gf = -127235. HCl(0) Regression of Sverjensky et al. (1991) & Tagirov et al. (1997) with new full correlations H2O(.5)Cl2(.5)O2(-.25) G0 = -127235.4 S0 = 13.3888 w = -83680 a1 = .5249665 a2 = -1973.886 a3 = 31.81639 a4 = -108114.6 c1 = 69.92887 c2 = 120193.8 end HCO3- EoS = 16 | Gf = -586940. HCO3(-) Facq et al. (2014) H2O(.5)CO2(1)O2(.25) G0 = -586939.9 S0 = 98.44952 w = 532748.7 q = -1 a1 = 3.20076 a2 = 384.928 a3 = 2.5104 a4 = -117988.8 c1 = 46.024 c2 = -158992 end HCrO4- EoS = 16 | Gf = -764835. HCrO4(-) Shock et al. (1997) revised with new predicted a1 for complex species Cr2O3(.5)H2O(.5)O2(1) G0 = -764835.2 S0 = 184.096 w = 402584.5 q = -1 a1 = 3.393726 a2 = 5030.8 a3 = 4.284834 a4 = -137072 c1 = 112.9144 c2 = 46898.46 end HEXANE,A EoS = 16 | Gf = 18493. C6H14 "Plyasunov & Shock (2001) for c1, c2 & omega; a1 predicted with new correlation for complex species from Sverjensky et al. (2014)" H2O(.5)CO2(6)O2(-6.25) G0 = 18493.28 S0 = 222.5888 w = -196648 a1 = 9.805162 a2 = 7543.417 a3 = -43.00817 a4 = -147456.7 c1 = 451.872 c2 = 962320 end HO2- EoS = 16 | Gf = -67362. HO2(-) Shock et al. (1997) revised with new predicted a1 for complex species H2O(.5)O2(.75) G0 = -67362.4 S0 = 23.8488 w = 646386.2 q = -1 a1 = 1.283275 a2 = -122.3402 a3 = 24.53874 a4 = -115767.1 c1 = 11.29973 c2 = -384354.8 end HS- EoS = 16 | Gf = 11966. HS(-) Shock et al. (1997) revised with new predicted a1 for complex species H2O(.5)S2(.5)O2(-.25) G0 = 11966.24 S0 = 68.1992 w = 602914.4 q = -1 a1 = 2.096979 a2 = 2083.59 a3 = 14.54568 a4 = -124888.2 c1 = 14.30928 c2 = -262336.8 end HSiO3- EoS = 16 | Gf = -1015879. HSiO3(-) "Sverjensky et al. (1997), revised to predict a1 using new correlation for complex species in Sverjensky et al. (2014)" SiO2(1)H2O(.5)O2(.25) G0 = -1015879. S0 = 20.92 w = 648980.2 q = -1 a1 = 1.519838 a2 = 455.6794 a3 = 22.25511 a4 = -118152 c1 = 1.803304 c2 = -408216.1 end HSO3- EoS = 16 | Gf = -527728. HSO3(-) Sverjensky et al. (1997) revised using published estimated V and pred. a1 from new correlation H2O(.5)S2(.5)O2(1.25) G0 = -527727.9 S0 = 139.7456 w = 469988.7 q = -1 a1 = 2.803866 a2 = 3590.541 a3 = 9.945786 a4 = -131118.2 c1 = 65.66746 c2 = -138900.4 end HSO4- EoS = 16 | Gf = -755756. HSO4(-) Shock et al. (1997) revised with new predicted a1 for complex species H2O(.5)S2(.5)O2(1.75) G0 = -755755.9 S0 = 125.52 w = 491536.3 q = -1 a1 = 2.91993 a2 = 3873.966 a3 = 8.831587 a4 = -132289.7 c1 = 84.08208 c2 = -81797.2 end HSO5- EoS = 16 | Gf = -637516. HSO5(-) Shock et al. (1997) revised with new predicted a1 for complex species H2O(.5)S2(.5)O2(2.25) G0 = -637516.1 S0 = 212.1288 w = 360284.2 q = -1 a1 = 3.763006 a2 = 5932.494 a3 = .7409864 a4 = -140800 c1 = 130.5935 c2 = 121896.7 end K+ EoS = 16 | Gf = -282462. K(+) Shock & Helgeson (1988) K2O(.5)O2(-.25) G0 = -282461.8 S0 = 101.0436 w = 80625.68 q = 1 a1 = 1.489086 a2 = -616.3032 a3 = 22.74004 a4 = -113470.1 c1 = 30.9616 c2 = -74935.44 end KCl,aq EoS = 16 | Gf = -413379. KCl(0) Fitted to Ho & Palmer (1997) data with new a1 pred. from the volume K2O(.5)Cl2(.5)O2(-.25) G0 = -413379.2 S0 = 184.096 w = 167360 a1 = 1.819203 a2 = 1186.331 a3 = 19.39535 a4 = -121177 c1 = -10.4554 c2 = -246282.8 end KOH,aq EoS = 16 | Gf = -430115. KOH(0) Fitted to Ho & Palmer (1997) data with a1 pred. from the sum of the ions and used to predict the volume K2O(.5)H2O(.5) G0 = -430115.2 S0 = 151.8792 w = 200832 a1 = 2.471572 a2 = 561.1162 a3 = 15.5377 a4 = -118591.3 c1 = 8.116542 c2 = -242027.7 end KSO4- EoS = 16 | Gf = -1031649. KSO4(-) "Revised G, S and Cp to fit Noyes et al. (1906) & Sharygin et al. (2006) with a1 of the complex equal to the sum of the a1 values of the i K2O(.5)S2(.5)O2(1.75) G0 = -1031649. S0 = 142.256 w = 460072.6 q = -1 a1 = 2.542868 a2 = 2952.649 a3 = 12.46455 a4 = -128473.9 c1 = 41.21324 c2 = -220714.4 end LACTATE, EoS = 16 | Gf = -512666. C3H5O3- Shock & Helgeson (1990) revised with new predicted a1 for ions H2O(2.5)CO2(3)O2(-3.75) G0 = -512665.5 S0 = 133.4696 w = 479863 q = -1 a1 = 4.121156 a2 = 5365.729 a3 = 33.88622 a4 = -138452.7 c1 = 188.0783 c2 = -158251.4 end LACTIC-A EoS = 16 | Gf = -534715. C3H6O3 Shock & Helgeson (1990) revised with new predicted a1 for complex species H2O(3)CO2(3)O2(-3) G0 = -534715.2 S0 = 208.3632 w = -107612.5 a1 = 4.674365 a2 = 7214.555 a3 = 15.93393 a4 = -146092.7 c1 = 255.3073 c2 = 68157.36 end METHANE, EoS = 16 | Gf = -34451. CH4 "Plyasunov & Shock (2001) for c1, c2 & omega; a1 predicted with new correlation for complex species from Sverjensky et al. (2014)" H2O(2)CO2(1)O2(-2) G0 = -34451.06 S0 = 87.82216 w = -167360 a1 = 3.481088 a2 = 1522.976 a3 = 7.48936 a4 = -122591.2 c1 = 171.1256 c2 = 269868 end METHANOL EoS = 16 | Gf = -175937. CH3OH Plyasunov & Shock (2001) H2O(1.5)CO2(1)O2(-1.25) G0 = -175937.2 S0 = 132.2144 w = -62760 a1 = 3.51456 a2 = 1715.44 a3 = -1.00416 a4 = -59831.2 c1 = 143.9296 c2 = 40166.4 end Mg(HCO3) EoS = 16 | Gf = -1046837. Mg(HCO3)(+) Sverjensky et al. (1997) revised V increased so that a1 of the complex is the sum of the a1 values of the ions MgO(1)H2O(.5)CO2(1)O2(-.25) G0 = -1046837. S0 = -12.552 w = 250412.4 q = 1 a1 = 1.368712 a2 = 86.65064 a3 = 23.71115 a4 = -116629 c1 = 197.8337 c2 = 390794 end MgHSiO3+ EoS = 16 | Gf = -1477057. Mg(HSiO3)(+) Sverjensky et al. (1997) with new a1 from the sum of the ions and V predicted from a1 value MgO(1)SiO2(1)H2O(.5)O2(-.25) G0 = -1477057. S0 = -99.49552 w = 383965.7 q = 1 a1 = .2631318 a2 = -2611.988 a3 = 34.29499 a4 = -105474.5 c1 = 153.9218 c2 = 195401.2 end Mg+2 EoS = 16 | Gf = -453985. Mg(+2) Shock & Helgeson (1988) MgO(1)O2(-.5) G0 = -453984.9 S0 = -138.072 w = 643164.5 q = 2 a1 = -.3437993 a2 = -3597.822 a3 = 35.10376 a4 = -99997.6 c1 = 87.0272 c2 = -246521.3 end MgCl+ EoS = 16 | Gf = -584505. MgCl(+) Sverjensky et al. (1997) revised V increased so that a1 of the complex is the sum of the a1 values of the ions MgO(1)Cl2(.5)O2(-.5) G0 = -584504.8 S0 = -79.496 w = 353506.2 q = 1 a1 = 1.343064 a2 = -2058.528 a3 = 26.23368 a4 = -107947.2 c1 = 119.6624 c2 = 86106.72 end MgCO3,aq EoS = 16 | Gf = -998972. MgCO3(0) "Sverjensky et al. (1997) revised volume increased in order that a1 of the complex is the sum of the a1 values of the ions; and MgO(1)CO2(1) G0 = -998971.8 S0 = -100.416 w = -15899.2 a1 = -.2281535 a2 = -3810.118 a3 = 38.96852 a4 = -100516.4 c1 = -43.92991 c2 = -364238.1 end MgOH+ EoS = 16 | Gf = -624357. MgOH(+) Shock et al. (1997) with revised a1 equal to the sum of the ions used to predict new values of V using new correlation in Sverjensky et al MgO(1)H2O(.5)O2(-.25) G0 = -624357.4 S0 = -79.9144 w = 353506.2 q = 1 a1 = .9667132 a2 = -893.9116 a3 = 27.54202 a4 = -112574.7 c1 = 133.8913 c2 = 135536.5 end MgSO4,aq EoS = 16 | Gf = -1211055. MgSO4 "S, Cp, V, omega from fit to Akilan et al. (2006a) and Frantz et al. (1994)" MgO(1)S2(.5)O2(1.5) G0 = -1211055. S0 = -4.184 w = 20920 a1 = .810399 a2 = -1278.086 a3 = 29.10139 a4 = -110984.8 c1 = -28.61438 c2 = -312080.4 end Mn+2 EoS = 16 | Gf = -228028. Mn(+2) Shock et al. (1997) revised with new predicted a1 for cations MnO(1)O2(-.5) G0 = -228028 S0 = -73.6384 w = 586011 q = 2 a1 = -.4250944E-1 a2 = -3359.543 a3 = 37.2627 a4 = -102386.7 c1 = 69.7364 c2 = -161912.4 end MnCl+ EoS = 16 | Gf = -361037. MnCl(+) Sverjensky et al. (1997) with revised a1 equal to the sum of the ions used to predict new values of V using new correlation in Sverjensky et a MnO(1)Cl2(.5)O2(-.5) G0 = -361037.4 S0 = 50.208 w = 154222.2 q = 1 a1 = 1.230305 a2 = -250.4124 a3 = 25.01363 a4 = -115235.7 c1 = 100.4026 c2 = 82947.8 end MnO4- EoS = 16 | Gf = -447270. MnO4(-) Shock et al. (1997) revised with new predicted a1 for ions MnO(1)O2(1.5) G0 = -447269.6 S0 = 191.2088 w = 391957.1 q = -1 a1 = 3.275612 a2 = 4742.397 a3 = 5.41828 a4 = -135879.6 c1 = 57.49946 c2 = -142310.4 end MnO4-2 EoS = 16 | Gf = -500825. MnO4(-2) Shock et al. (1997) revised with new predicted a1 for ions MnO(1)O2(1.5) G0 = -500824.8 S0 = 58.576 w = 1256204. q = -2 a1 = 2.731775 a2 = 3414.521 a3 = 10.6374 a4 = -130390.2 c1 = -3.877313 c2 = -632369.8 end MnSO4,aq EoS = 16 | Gf = -985918. MnSO4(0) Sverjensky et al. (1997) with revised a1 equal to the sum of the ions used to predict new values of V using new correlation in Sver MnO(1)S2(.5)O2(1.5) G0 = -985917.8 S0 = 20.92 w = -15899.2 a1 = .9724453 a2 = -880.2299 a3 = 27.49014 a4 = -112629.1 c1 = -28.33405 c2 = -310034.4 end Na+ EoS = 16 | Gf = -261881. Na(+) Shock & Helgeson (1988) Na2O(.5)O2(-.25) G0 = -261880.7 S0 = 58.40864 w = 138323 q = 1 a1 = .7694376 a2 = -956.044 a3 = 13.6231 a4 = -114055.8 c1 = 76.06512 c2 = -124725 end NaCl,aq EoS = 16 | Gf = -388735. NaCl(0) Sverjensky et al. (1997) with new predicted a1 value Na2O(.5)Cl2(.5)O2(-.25) G0 = -388735.4 S0 = 117.152 w = -15899.2 a1 = 2.107188 a2 = 1981.752 a3 = 14.29003 a4 = -124465.6 c1 = 45.1872 c2 = -54392 end NaHCO3,a EoS = 16 | Gf = -849720. NaHCO3(0) "Fit to Stefansson et al. (2013; 80 to 200 C) & Nakayama (1971; at 25 C) to get G, S and Cp; with a1 pred. from the sum of the io Na2O(.5)H2O(.5)CO2(1) G0 = -849720.2 S0 = 154.808 w = -15899.2 a1 = 2.092 a2 = 4184 a3 = 8.368 a4 = -83680 c1 = 83.68 c2 = 209200 end NaHSiO3, EoS = 16 | Gf = -1285074. NaHSiO3(0) Sverjensky et al. (1997) with revised a1 equal to the sum of the ions used to predict new values of V using new correlation in Sve Na2O(.5)SiO2(1)H2O(.5) G0 = -1285074. S0 = 79.496 w = -15899.2 a1 = 1.747071 a2 = 1009.097 a3 = 20.11207 a4 = -120440.6 c1 = 30.07334 c2 = -107018.4 end NaOH,aq EoS = 16 | Gf = -422166. NaOH(0) Fitted to Ho & Palmer (1997) data with a1 pred. from the sum of the ions and used to predict the volume Na2O(.5)H2O(.5) G0 = -422165.6 S0 = 27.196 w = 83680 a1 = 1.651383 a2 = -219.7855 a3 = 22.08566 a4 = -115365.4 c1 = 136.498 c2 = 230990.3 end Ni+2 EoS = 16 | Gf = -45606. Ni(+2) Shock & Helgeson (1988) revised with new predicted a1 for ions NiO(1)O2(-.5) G0 = -45605.6 S0 = -128.8672 w = 630403.3 q = 2 a1 = -.7088533 a2 = -4986.533 a3 = 43.65753 a4 = -95658.79 c1 = 55.18905 c2 = -226684.9 end O2,aq EoS = 16 | Gf = 16544. O2(0) Shock et al. (1989) O2(1) G0 = 16543.54 S0 = 108.9514 w = -164975.1 a1 = 2.422076 a2 = 2658.346 a3 = 13.60972 a4 = -127264.7 c1 = 147.917 c2 = 350309.6 end OH- EoS = 16 | Gf = -157297. OH(-) "August 16th, 2013 fit to Bandura & Lvov (2005)" H2O(.5)O2(.25) G0 = -157297.5 S0 = -20.92 w = 719648 q = -1 HOH = 2 a1 = .87864 a2 = -2092 a3 = 4.184 a4 = -112968 c1 = 50.71008 c2 = -314636.8 end PROPANE, EoS = 16 | Gf = -8213. C3H8 Shock & Helgeson (1990) with new a1-a4 based on revised correlation to predict a1 in Sverjensky et al. (2013) H2O(4)CO2(3)O2(-5) G0 = -8213.192 S0 = 139.6201 w = -211417.5 a1 = 4.50303 a2 = 7738.224 a3 = -6.331647 a4 = -148260 c1 = 277.5222 c2 = 815654.1 end PROPANOA EoS = 16 | Gf = -363046. C2H5COO- Revised a1 from new delVn correlation for -1 ions; kept original regression c1 and c2 plus omega from Shock & Helgeson (1990) H2O(2.5)CO2(3)O2(-3.25) G0 = -363045.7 S0 = 110.876 w = 513627.8 q = -1 a1 = 3.505606 a2 = 8090.392 a3 = .192464 a4 = -149720.3 c1 = 218.8232 c2 = -175728 end PROPANOI EoS = 16 | Gf = -390995. C3H6O2 Plyasunov & Shock (2001) H2O(3)CO2(3)O2(-3.5) G0 = -390994.8 S0 = 206.6896 w = -37656 a1 = 6.10864 a2 = 3305.36 a3 = 28.8696 a4 = -179912 c1 = 264.8472 c2 = -58576 end PROPANOL EoS = 16 | Gf = -175351. C3H7OH Plyasunov & Shock (2001) H2O(3.5)CO2(3)O2(-4.25) G0 = -175351.4 S0 = 169.8704 w = -83680 a1 = 6.10864 a2 = 3807.44 a3 = 51.0448 a4 = -230120 c1 = 310.8712 c2 = 138072 end S2-2 EoS = 16 | Gf = 79496. S2(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1) G0 = 79496 S0 = 28.4512 w = 1300513. q = -2 a1 = 2.334546 a2 = 2444.544 a3 = 14.44986 a4 = -126377.7 c1 = -14.01473 c2 = -681803.7 end S2O3-2 EoS = 16 | Gf = -522582. S2O3(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1)O2(1.5) G0 = -522581.6 S0 = 66.944 w = 1242397. q = -2 a1 = 2.7901 a2 = 5227.95 a3 = -32.33437 a4 = -137883.7 c1 = -.2414168 c2 = -615324.1 end S2O4-2 EoS = 16 | Gf = -600404. S2O4(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1)O2(2) G0 = -600404 S0 = 92.048 w = 1203820. q = -2 a1 = 3.1721 a2 = 4489.683 a3 = 6.411562 a4 = -134833.6 c1 = 6.990209 c2 = -577823 end S2O5-2 EoS = 16 | Gf = -790776. S2O5(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1)O2(2.5) G0 = -790776 S0 = 104.6 w = 1185871. q = -2 a1 = 3.337786 a2 = 4894.276 a3 = 4.821642 a4 = -136507.2 c1 = 11.22316 c2 = -557371.6 end S2O6-2 EoS = 16 | Gf = -966504. S2O6(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1)O2(3) G0 = -966504 S0 = 125.52 w = 1154240. q = -2 a1 = 3.607654 a2 = 5553.172 a3 = 2.231746 a4 = -139231 c1 = 18.11714 c2 = -523276.1 end S2O8-2 EoS = 16 | Gf = -1115036. S2O8(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1)O2(4) G0 = -1115036 S0 = 244.3456 w = 974077 q = -2 a1 = 5.590744 a2 = 10395.32 a3 = -16.80002 a4 = -159247.2 c1 = 54.23803 c2 = -340037.9 end S3-2 EoS = 16 | Gf = 73638. S3(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1.5) G0 = 73638.4 S0 = 66.1072 w = 1244698. q = -2 a1 = 2.830936 a2 = 3656.649 a3 = 9.68596 a4 = -131390.2 c1 = -1.504148 c2 = -620437 end S3O6-2 EoS = 16 | Gf = -958136. S3O6(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(1.5)O2(3) G0 = -958136 S0 = 138.072 w = 1135161. q = -2 a1 = 3.778696 a2 = 5970.819 a3 = .589944 a4 = -140954.8 c1 = 22.24591 c2 = -502824.8 end S4-2 EoS = 16 | Gf = 69036. S4(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(2) G0 = 69036 S0 = 103.3448 w = 1187838. q = -2 a1 = 3.321301 a2 = 4853.942 a3 = 4.979797 a4 = -136339.8 c1 = 10.91229 c2 = -559074.4 end S4O6-2 EoS = 16 | Gf = -1040561. S4O6(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(2)O2(3) G0 = -1040561. S0 = 257.316 w = 954161.2 q = -2 a1 = 5.349244 a2 = 9805.706 a3 = -14.4825 a4 = -156812.1 c1 = 61.47677 c2 = -308503.1 end S5-2 EoS = 16 | Gf = 65689. S5(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(2.5) G0 = 65688.8 S0 = 140.5824 w = 1131814. q = -2 a1 = 3.811917 a2 = 6051.947 a3 = .2715416 a4 = -141293.7 c1 = 23.16304 c2 = -498561.3 end S5O6-2 EoS = 16 | Gf = -958136. S5O6(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(2.5)O2(3) G0 = -958136 S0 = 167.36 w = 1091020. q = -2 a1 = 4.164544 a2 = 6912.972 a3 = -3.112896 a4 = -144850.1 c1 = 31.90969 c2 = -455093.7 end Si2O4,aq EoS = 16 | Gf = -1676742. Si2O4(0) "Fit to Raman speciation and quartz solubility data (Sverjensky et al., 2014)." SiO2(2) G0 = -1676742. S0 = 75.312 w = 41840 a1 = 4.214292 a2 = 1126.793 a3 = 2.809556 a4 = -120930.2 c1 = 149.8077 c2 = 290649.9 end SiO2,aq EoS = 16 | Gf = -834946. SiO2(0) "Fit to Raman speciation and quartz solubility data (Sverjensky et al., 2014)." SiO2(1) G0 = -834946.5 S0 = 22.1752 w = 150624 a1 = 2.050076 a2 = 589.5256 a3 = 18.43721 a4 = -118708.4 c1 = 107.6003 c2 = 109114.5 end SO2,aq EoS = 16 | Gf = -301164. SO2(0) Shock et al. (1989) with revised a1 predicted as a complex from delVn S2(.5)O2(1) G0 = -301164.3 S0 = 161.9208 w = -102968.2 a1 = 2.907964 a2 = 3844.678 a3 = 8.946647 a4 = -132168.4 c1 = 130.5831 c2 = 270194.4 end SO3-2 EoS = 16 | Gf = -486599. SO3(-2) Shock & Helgeson (1988) revised with new predicted a1 for ions S2(.5)O2(1.5) G0 = -486599.2 S0 = -29.288 w = 1389506. q = -2 a1 = 1.528708 a2 = 133.5114 a3 = 30.9001 a4 = -116825.6 c1 = -32.78917 c2 = -775554.6 end SO4-2 EoS = 16 | Gf = -744459. SO4(-2) Shock & Helgeson (1988) S2(.5)O2(2) G0 = -744459.1 S0 = 18.828 w = 1316412. q = -2 a1 = 3.473306 a2 = -830.3566 a3 = -25.99184 a4 = -112842.5 c1 = 6.86176 c2 = -753036.3 end TOLUENE, EoS = 16 | Gf = 126608. C6H5CH3 Plyasunov & Shock (2001) with new a1 to a4 predicted with revised a1 consistent with Sverjensky et al. (2013) H2O(4)CO2(7)O2(-9) G0 = 126607.8 S0 = 184.096 w = -121336 a1 = 8.500549 a2 = 6301.271 a3 = -32.59294 a4 = -142322.9 c1 = 341.8328 c2 = 598312 end