A new activity model for biotite and its application

Area of Science:

• Geochemistry
• Mineral Physics
• Computational Materials Science

Background:

• Biotite is a key mineral in metamorphic petrology, but accurately modeling its thermodynamic properties remains challenging.
• Existing biotite activity models often lack precision due to simplified assumptions about mixing behavior.
• The KFMASHTO system requires new thermodynamic models to account for titanium and ferric end-members in biotite.

Purpose of the Study:

• To develop a novel, physically-based activity model for biotite in the KFMASHTO system.
• To incorporate titanium-biotite, ferric-biotite, and pyrophyllite end-members into the biotite model.
• To utilize density functional theory (DFT) to derive microscopic interaction parameters for enhanced thermodynamic calculations.

Main Methods:

• Formulation of a new biotite activity model within the KFMASHTO system.
• Application of density functional theory (DFT) using Castep software for phonon calculations and entropy determination.
• Parameterization of macroscopic mixing properties (macro-W's) using single-defect DFT and experimental data.
• Derivation of microscopic interaction energies (micro-w's) for Mg-Al and Si-Al mixing in biotite.

Main Results:

• Standard entropies and heat capacity functions for titanium-biotite and ferric-biotite end-members were calculated using DFT.
• Enthalpies of formation for titanium-biotite and natural biotite were constrained using experimental phase-equilibrium data.
• DFT-derived microscopic interaction energies (micro-w's) were obtained for key mixing parameters in biotite.
• The new model, incorporating DFT-based micro-w's, showed improved agreement with measured mineral compositions in test cases.

Conclusions:

• The developed biotite activity model represents a next-generation approach, integrating physically-based parameters from DFT.
• The model demonstrates superior performance compared to existing models in predicting biotite mineral chemistry across various metamorphic conditions.
• This physically-grounded model enhances the accuracy of petrological calculations and thermodynamic modeling of metamorphic processes.