Thermodynamic modeling of melts in the system Na 2O-NaAlO 2-SiO 2-F 2O −1

Abstract

Fluorine is a common volatile element in magmatic-hydrothermal systems, but its solution mechanisms and thermodynamic description in highly polymerized silicate melts are poorly known. We have developed a thermodynamic model for fluorosilicate liquids that links experimentally determined phase equilibria and spectroscopic information on melt structure. The model is applicable to crystallization of fluoride minerals, fluoride-silicate immiscibility in natural felsic melts, and metallurgical processes. Configurational properties of fluorosilicate melts are described by mixing on three site levels (sublattices): (1) alkali fluoride, polyhedral aluminofluoride and silicofluoride species and nonbridging terminations of the aluminosilicate network, (2) alkali-aluminate and silicate tetrahedra within the network and (3) bridging oxygen, nonbridging oxygen and terminal fluorine atoms on tetrahedral apices of the network. Abundances of individual chemical species are described by a homogeneous equilibrium representing melt depolymerization: F − (free) + O 0 (bridging) = F 0 (terminal) + O − (nonbridging) which corresponds to a replacement of an oxygen bridging two tetrahedra by a pair of terminations, one with F and the other with an O and a charge-balancing Na. In cryolite-bearing systems two additional interaction mechanisms occur: (1) the self-dissociation of octahedral aluminofluoride complexes: [AlF 6] = [AlF 4] + 2 [F], and (2) the short-range order between (O,F)-corners and (Si,NaAl)-centers of tetrahedra: Si-O-Si + 2 [NaAl]-F = [NaAl]-O-[NaAl] + 2 Si-F. Portrayal of these equilibria in ternary Thompson reaction space allows for the decrease in the number of interaction mechanisms by linearly combining melt depolymerization with tetrahedral short-range order. In this formulation, the nonideal thermodynamic properties are represented by reaction energies of homogeneous equilibria, thus defining directly individual chemical species concentrations and configurational properties. Thermodynamic expressions for the activity-composition relationships are simplified if all entities are expressed using symbolic molecular notation (e.g., SiO 2, SiF 4, [NaAl]O 2, [NaAl]F 4, NaF etc.) with corresponding nonfractional site multiplicities (1, 2 or 4). The model has been applied to three subsystems of the Na 2O-NaAlO 2-SiO 2-F 2O −1 compositional space. Activity-composition relationships in the villiaumite-sodium silicate binaries require clustering of silicate tetrahedra and only negligible interaction between fluoride species and silicate polymer. Phase equilibria in the cryolite-albite system with a large depression of albite liquidus are interpreted via complete substitution of O 0 by O − and F 0 in the silicate framework. With increasing fluorine content, initial Al-F and Si-O short-range order evolves into the partial O-F disorder. The present model provides a useful relationship between experimental equilibria, macroscopic thermodynamics and melt speciation, thus it facilitates comparisons with, and interpretations of, spectroscopic and molecular simulation data.