Úkol: Vykreslete solvus mezi muskovitem a paragonitem při tlaku 2000 bar za použití SS modelu Coggona a Hollanda (2002). "Uvažujte, že na míšení se nepodílejí další komponenty." Postup: Princip solvní termobarometrie pro pár muskovit-paragonit je obdobný (respektive jednodušší) solvní živcové termometrii. Znovu se využívá rovnosti aktivit (respektive rovnosti chemických potenciálů) komponent v obou slídách při jejich rovnováze. Ze dvou rovnic pro rovnost chem. potenciálů lze vyjádřit teplotu - "jednu pro aktivitu muskovitu, druhou pro aktivitu paragonitu. Za rovnováhy se tyto teploty rovnají. Výpočet teplot z dvojice slídje navržen tak, že vložíme rovnovážná " složení obou koexistujících slíd a vypočtené teploty by měly být shodné. Pro vykreslení solvu zvolíme odlišný přístup za použití stejného výpočtu. "Využijeme doplňku řešitel, který nám umožní iterovat hodnotu určité buňky tak, abychom dosáhli cílové hodnoty v buňce jiné. Konkrétně budeme měnit složení" "muskovitu, zatímco rovnovážné složení paragonitu necháme iterovat. Jako cílovou buňku zvolíme F12, která obsahuje vzorec pro rozdíl dvou teplot" " (viz výše), zatímco měněná buňka bude I34. Program necháme proiterovat různá složení muskovitu." amu T (°C) T (K) R 8.3144598 199 472 P (bar) 2000 0 apa T (°C) T (K) 199 472 ACTIVITIES MUSCOVITE MUSCOVITE proportions site occupancy prl ma phl fcel cel ideal activity W(mu-pa) 12430.62631 activity x y z na ca p Na Ca K Fe Mg Si Al ppa pprl pma pphl pmu pfcel pcel "K,A" "Na,A" "Ca,A" "vac,M1" "Mg,M1" "Al,M2A" "Mg,M2A" "Fe,M2A" "Al,M2B" "Mg,M2B" "Al,T1" "Si,T1" sum(A) sum(M1) sum(M2A) sum(M2B) sum(T1) *RTlnγpa RTlnγprl RTlnγma *RTlnγmu RTlnγfcel RTlnγcel mu pa RTlnγmu RTlnγpa γmu γpa mu pa 0 0.5 0 0.35 0 0 0.35 0.65 3 3 0.35 0 0 0 0.65 0 0 0.65 0.35 0 1 0 1 0 0 1 0 0.5 0.5 1 1 1 1 1 4624.129885 18887.66833 787.4093625 13555.79819 13555.79819 0.65 0.35 904.1189305 5309.516962 1.259112941 3.86935494 0.818423411 1.354274229 pa prl ma phl mu fcel cel 19349.34153 assymetr Ф 0.24025974 0 0 0 0.75974026 0 0 WHITE MICA formalism α 0.37 0.5 0.37 0.63 0.63 0.63 0.63 "(T),P-dependent" p*α 0.1295 0 0 0 0.4095 0 0 W=WH+PWV pa ma cel fcel prl WS mu pa mu 10826 30000 400 400 20400 mu-pa -3.4 -0.247292798 -1.45224844 pa 14500 52000 52000 20400 ma 30000 30000 30400 cel 0 25400 fcel 25400 ACTIVITIES PARAGONITE PARAGONITE proportions site occupancy prl ma phl fcel cel ideal activity W(mu-pa) 12430.62624 activity x y z na ca p Na Ca K Fe Mg Si Al ppa pprl pma pphl pmu pfcel pcel "K,A" "Na,A" "Ca,A" "vac,M1" "Mg,M1" "Al,M2A" "Mg,M2A" "Fe,M2A" "Al,M2B" "Mg,M2B" "Al,T1" "Si,T1" sum(A) sum(M1) sum(M2A) sum(M2B) sum(T1) *RTlnγpa RTlnγprl RTlnγma *RTlnγmu RTlnγfcel RTlnγcel mu pa RTlnγmu RTlnγpa γmu γpa mu pa 0 0.5 0 0.349794191 0 0 0.349794191 0.650205809 3 3 0.349794191 0 0 0 0.650205809 0 0 0.650205809 0.349794191 0 1 0 1 0 0 1 0 0.5 0.5 1 1 1 1 1 4626.140034 18889.62709 786.3274642 13546.21631 13546.21631 0.650205809 0.349794191 902.8766678 5311.825021 1.258714407 3.871631784 0.818423419 1.354274308 pa prl ma phl mu fcel cel 19348.45071 assymetr Ф 0.240094625 0 0 0 0.759905375 0 0 WHITE MICA formalism α 0.37 0.5 0.37 0.63 0.63 0.63 0.63 "(T),P-dependent" Řešení: p*α 0.129423851 0 0 0 0.40962966 0 0 W=WH+PWV pa ma cel fcel prl WS mu pa mu 10826 30000 400 400 20400 mu-pa -3.4 -0.246953018 -1.452879745 XNa pa 14500 52000 52000 20400 mus par T ma 30000 30000 30400 0.050 0.995 62 cel 0 25400 0.100 0.987 146 fcel 25400 0.150 0.975 213 0.200 0.962 270 0.250 0.946 321 0.300 0.929 367 0.350 0.909 408 0.400 0.886 446 0.450 0.860 479 0.500 0.832 508 0.550 0.799 531 0.600 0.762 549 0.650 0.720 560 0.700 0.700 561