{"title":"A new activity model for biotite and its application","authors":"Edgar Dachs, Artur Benisek","doi":"10.1007/s00410-024-02173-6","DOIUrl":null,"url":null,"abstract":"<div><p>A new activity model for biotite is formulated in the system K<sub>2</sub>O-FeO-MgO-Al<sub>2</sub>O<sub>3</sub>-SiO<sub>2</sub>-H<sub>2</sub>O-TiO<sub>2</sub>-O<sub>2</sub> (KFMASHTO), which extends that for the KFMASH system by introducing a titanium-biotite and a ferric-biotite end-member (<i>tbio</i>: K(TiMg<sub>2</sub>)[(O)<sub>2</sub>(AlSi<sub>3</sub>)O<sub>10</sub>] and <i>fbio</i>: K(Fe<sup>3+</sup>Mg<sub>2</sub>)[(OH)<sub>2</sub>(Al<sub>2</sub>Si<sub>2</sub>)O<sub>10</sub>]), as well as a pyrophyllite end-member (<i>pyp</i>: Al<sub>2</sub>[(OH)<sub>2</sub>Si<sub>4</sub>O<sub>10</sub>]) that accounts for the presence of octahedral excess-Al in natural biotites. Phonon calculations applying density functional theory (DFT) using the software <i>Castep</i> yielded the standard entropies of <i>tbio</i> and <i>fbio</i> as <i>S</i><sup>o</sup><sub>tbio</sub> = 328.06 J/(mol·K) and <i>S</i><sup>o</sup><sub>fbio</sub> = 301.69 J/(mol·K), and their heat capacity functions. From experimental phase-equilibrium data, the enthalpy of formation value of <i>tbio</i> was constrained as <span>\\(\\Delta {H}_{f,tbio}^{o}\\)</span> = −6124.68 ± 3.33 kJ/mol. Natural data were used to derive <span>\\(\\Delta {H}_{f,fbio}^{o}\\)</span>= −5935.3 ± 6.6 kJ/mol. The single-defect DFT method was applied to parameterize important macroscopic mixing properties (macro-<i>W</i>’s) involving <i>tbio</i> and <i>pyp</i> end-members in the model (<i>fbio</i> was treated ideal). <i>Castep</i>-derived microscopic interaction energies (micro-<i>w</i>’s) are presented herein for KFMASH-biotite. The octahedral same-site (M1) Mg–Al mixing micro-<i>w</i> (<i>w</i><sub>MgAl(M1)</sub>), the same-site tetrahedral Si-Al mixing parameter (<i>w</i><sub>SiAl(T1)</sub>) and the related cross-site term are: <i>w</i><sub>MgAl(M1)</sub> = 82.5 kJ/mol, <i>w</i><sub>SiAl(T1)</sub> = 95.6 kJ/mol (two T1-sites) and <span>\\({w}_{MgAlAlSi(M1T1)}=\\)</span> 175.1 kJ/mol. The linear combination of these micro-<i>w</i>’s gives a macroscopic <i>W</i><sub>phleas</sub> = 18.8 kJ/mol, that is not transferable to other mineral groups. Micro <i>w</i>’s for Mg-Fe mixing in biotite (<i>w</i><sub>MgFe(M1)</sub>, <i>w</i><sub>MgFe(M2),</sub> <span>\\({w}_{MgMgFeFe(M1M2)}\\)</span>), are all close to ideality. The biotite activity model of this study is thus a first example of next-generation activity models that use DFT- and thus physically based micro-<i>w</i>’s and reassembled macro-<i>W</i>’s for petrological calculations. Test calculations on 5 samples from low- to high-grade metamorphic environments covering metapelite to greywacke bulk-compositions using <i>Perple_X</i> suite of programs illustrate the performance of the new biotite activity model. Computed mineral-chemistries are in all cases in better agreement with measured compositions than resulting from published activity models of biotite.</p></div>","PeriodicalId":526,"journal":{"name":"Contributions to Mineralogy and Petrology","volume":"179 10","pages":""},"PeriodicalIF":3.5000,"publicationDate":"2024-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11452188/pdf/","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Contributions to Mineralogy and Petrology","FirstCategoryId":"89","ListUrlMain":"https://link.springer.com/article/10.1007/s00410-024-02173-6","RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
引用次数: 0
Abstract
A new activity model for biotite is formulated in the system K2O-FeO-MgO-Al2O3-SiO2-H2O-TiO2-O2 (KFMASHTO), which extends that for the KFMASH system by introducing a titanium-biotite and a ferric-biotite end-member (tbio: K(TiMg2)[(O)2(AlSi3)O10] and fbio: K(Fe3+Mg2)[(OH)2(Al2Si2)O10]), as well as a pyrophyllite end-member (pyp: Al2[(OH)2Si4O10]) that accounts for the presence of octahedral excess-Al in natural biotites. Phonon calculations applying density functional theory (DFT) using the software Castep yielded the standard entropies of tbio and fbio as Sotbio = 328.06 J/(mol·K) and Sofbio = 301.69 J/(mol·K), and their heat capacity functions. From experimental phase-equilibrium data, the enthalpy of formation value of tbio was constrained as \(\Delta {H}_{f,tbio}^{o}\) = −6124.68 ± 3.33 kJ/mol. Natural data were used to derive \(\Delta {H}_{f,fbio}^{o}\)= −5935.3 ± 6.6 kJ/mol. The single-defect DFT method was applied to parameterize important macroscopic mixing properties (macro-W’s) involving tbio and pyp end-members in the model (fbio was treated ideal). Castep-derived microscopic interaction energies (micro-w’s) are presented herein for KFMASH-biotite. The octahedral same-site (M1) Mg–Al mixing micro-w (wMgAl(M1)), the same-site tetrahedral Si-Al mixing parameter (wSiAl(T1)) and the related cross-site term are: wMgAl(M1) = 82.5 kJ/mol, wSiAl(T1) = 95.6 kJ/mol (two T1-sites) and \({w}_{MgAlAlSi(M1T1)}=\) 175.1 kJ/mol. The linear combination of these micro-w’s gives a macroscopic Wphleas = 18.8 kJ/mol, that is not transferable to other mineral groups. Micro w’s for Mg-Fe mixing in biotite (wMgFe(M1), wMgFe(M2),\({w}_{MgMgFeFe(M1M2)}\)), are all close to ideality. The biotite activity model of this study is thus a first example of next-generation activity models that use DFT- and thus physically based micro-w’s and reassembled macro-W’s for petrological calculations. Test calculations on 5 samples from low- to high-grade metamorphic environments covering metapelite to greywacke bulk-compositions using Perple_X suite of programs illustrate the performance of the new biotite activity model. Computed mineral-chemistries are in all cases in better agreement with measured compositions than resulting from published activity models of biotite.
在 K2O-FeO-MgO-Al2O3-SiO2-H2O-TiO2-O2 体系(KFMASHTO)中提出了一种新的生物黄铁矿活性模型,该模型通过引入钛-生物黄铁矿和铁-生物黄铁矿末端分子(tbio:K(TiMg2)[(O)2(AlSi3)O10]和 fbio:K(Fe3+Mg2)[(OH)2(Al2Si2)O10]),以及辉绿岩端粒(pyp:Al2[(OH)2Si4O10]),说明天然生物岩中存在八面体过量铝。利用密度泛函理论(DFT)软件 Castep 进行的声子计算得出了 tbio 和 fbio 的标准熵 S o tbio = 328.06 J/(mol-K) 和 S o fbio = 301.69 J/(mol-K),以及它们的热容函数。根据相平衡实验数据,tbio 的形成焓值被确定为 Δ H f , t b i o o = -6124.68 ± 3.33 kJ/mol。利用自然数据得出 Δ H f , f b i o o = -5935.3 ± 6.6 kJ/mol。应用单缺陷 DFT 方法对模型中涉及 tbio 和 pyp 末端分子(fbio 被视为理想)的重要宏观混合特性(宏观-W)进行了参数化。本文介绍了 KFMASH 生物沸石的 Castep 衍生微观相互作用能(微观-W's)。八面体同位(M1)镁-铝混合微w(w MgAl(M1))、同位四面体硅-铝混合参数(w SiAl(T1))和相关的跨位项分别为:w MgAl(M1) = 82.5 kJ/mol,w SiAl(T1) = 95.6 kJ/mol(两个 T1 位)和 w M g A l A l S i ( M 1 T 1 ) = 175.1 kJ/mol。这些微观 w 的线性组合得出了宏观 W phleas = 18.8 kJ/mol,但这并不能用于其他矿物组。生物岩中镁铁混合的微观 W(w MgFe(M1)、w MgFe(M2)、w M g M g F e F e ( M 1 M 2 ) )都接近理想状态。因此,本研究中的黑云母活性模型是下一代活性模型的第一个实例,这些模型使用基于 DFT 的微观-w 和重新组装的宏观-w 进行岩石学计算,因此也是基于物理的微观-w 和重新组装的宏观-w 进行岩石学计算。使用 Perple_X 套装程序对从低级到高级变质环境的 5 个样本进行了测试计算,涵盖了从元青石到灰岩的体成分,说明了新的生物活度模型的性能。在所有情况下,计算出的矿物化学成分与实测成分的一致性都优于已发表的生物橄榄岩活性模型:在线版本包含补充材料,可查阅 10.1007/s00410-024-02173-6。
期刊介绍:
Contributions to Mineralogy and Petrology is an international journal that accepts high quality research papers in the fields of igneous and metamorphic petrology, geochemistry and mineralogy.
Topics of interest include: major element, trace element and isotope geochemistry, geochronology, experimental petrology, igneous and metamorphic petrology, mineralogy, major and trace element mineral chemistry and thermodynamic modeling of petrologic and geochemical processes.