The thermodynamic properties of 154 mineral end-members, 13 silicate liquid end-members and 22 aqueous fluid species are presented in a revised and updated data set. The use of a temperature-dependent thermal expansion and bulk modulus, and the use of high-pressure equations of state for solids and fluids, allows calculation of mineral–fluid equilibria to 100 kbar pressure or higher. A pressure-dependent Landau model for order–disorder permits extension of disordering transitions to high pressures, and, in particular, allows the alpha–beta quartz transition to be handled more satisfactorily. Several melt end-members have been included to enable calculation of simple phase equilibria and as a first stage in developing melt mixing models in NCKFMASH. The simple aqueous species density model has been extended to enable speciation calculations and mineral solubility determination involving minerals and aqueous species at high temperatures and pressures. The data set has also been improved by incorporation of many new phase equilibrium constraints, calorimetric studies and new measurements of molar volume, thermal expansion and compressibility. This has led to a significant improvement in the level of agreement with the available experimental phase equilibria, and to greater flexibility in calculation of complex mineral equilibria. It is also shown that there is very good agreement between the data set and the most recent available calorimetric data.

An internally consistent thermodynamic data set for phases of petrological interest

The thermodynamic properties of 154 mineral end-members, 13 silicate liquid end-members and 22 aqueous fluid species are presented in a revised and updated data set. The use of a temperature-dependent thermal expansion and bulk modulus, and the use of high-pressure equations of state for solids and fluids, allows calculation of mineral-fluid equilibria to 100 kbar pressure or higher. A pressure-dependent Landau model for order-disorder permits extension of disordering transitions to high pressures, and, in particular, allows the alpha-beta quartz transition to be handled more satisfactorily. Several melt end- members have been included to enable calculation of simple phase equilibria and as a first stage in developing melt mixing models in NCKFMASH. The simple aqueous species density model has been extended to enable speciation calculations and mineral solubility determination involving minerals and aqueous species at high temperatures and pressures. The data set has also been improved by incorporation of many new phase equilibrium constraints, calorimetric studies and new measurements of molar volume, thermal expansion and compressibility. This has led to a significant improvement in the level of agreement with the available experimental phase equilibria, and to greater flexibility in calculation of complex mineral equilibria. It is also shown that there is very good agreement between the data set and the most recent available calorimetric data. kinetics which apply to determining directly the greatest majority of such equilibria in the laboratory, for forming solid solutions, and inclusion of aqueous and silicate melt species), and provides uncertainties especially at lower temperatures, as well as the diYculty of establishing reversals of reactions involving solid allowing the likely uncertainties on the results of thermodynamic calculations to be estimated. This is a solutions. The levels of precision and accuracy required of thermodynamic data in order to be able to forward- critical issue in that calculations using data sets should always involve uncertainty propagation to help evalu- model synthetic and natural mineral assemblages mean that the continuing upgrading and expansion of the ate the results. Because the experimental phase equilib- ria involve overlapping subsets of compositional space, data set by incorporation of new phase equilibrium constraints, calorimetry and new measurements of the derived thermodynamic data are highly correlated, and it is only the inclusion of the correlations which molar volume, thermal expansion and compressibility are more than justified. enables the reliable calculation of uncertainties on mineral reactions to be performed. Earlier work on mineral thermodynamic data sets for rock-forming minerals includes compilations of The thermodynamic data extraction involves using weighted least squares on the diVerent types of data

The thermodynamic properties of 254 end-members, including 210 mineral end-members, 18 silicate liquid end-members and 26 aqueous fluid species are presented in a revised and updated internally consistent thermodynamic data set. The PVT properties of the data set phases are now based on a modified Tait equation of state (EOS) for the solids and the Pitzer & Sterner (1995) equation for gaseous components. Thermal expansion and compressibility are linked within the modified Tait EOS (TEOS) by a thermal pressure formulation using an Einstein temperature to model the temperature dependence of both the thermal expansion and bulk modulus in a consistent way. The new EOS has led to improved fitting of the phase equilibrium experiments. Many new end-members have been added, including several deep mantle phases and, for the first time, sulphur-bearing minerals. Silicate liquid end-members are in good agreement with both phase equilibrium experiments and measured heat of melting. The new dataset considerably enhances the capabilities for thermodynamic calculation on rocks, melts and aqueous fluids under crustal to deep mantle conditions. Implementations are already available in thermocalc to take advantage of the new data set and its methodologies, as illustrated by example calculations on sapphirine-bearing equilibria, sulphur-bearing equilibria and calculations to 300 kbar and 2000 °C to extend to lower mantle conditions.

An improved and extended internally consistent thermodynamic dataset for phases of petrological interest, involving a new equation of state for solids

This paper is the fourth in a series relating the lattice vibrational properties to the thermodynamic properties of minerals. The temperature dependence of the harmonic lattice heat capacity is calculated from a model which uses only elastic, crystallographic, and spectroscopic data for the following minerals: calcite, zircon, forsterite, grossular, pyrope, almandine, spessartine, andradite, kyanite, andalusite, sillimanite, clinoenstatite, orthoenstatite, jadeite, diopside, tremolite, talc, and muscovite. The heat capacities of these minerals reflect structural and compositional differences. The ‘excess’ entropy of pyrope—compared with that of grossular—is shown to arise from low-frequency optic modes of vibration. The entropy differences between kyanite, andalusite, and sillimanite are well reproduced by the model, although the absolute values calculated are systematically about 3% high. Model values of the heat capacity and entropy are compared with experimental values at 298.15, 700, and 1000°K for the 32 minerals included in papers 1–4 of this series. The average deviation of the entropies at 298°K from well-determined calorimetric values is ±1.5%. A method is given for obtaining greater accuracy in the model thermodynamic functions by fitting one parameter to experimental data when partial calorimetric data (such as the heat capacity at a single temperature in the range 50–100°K) are available; such a method should permit accurate extrapolation of calorimetric data beyond the range of experiment.

J. K. Winter

Papers

A high-temperature structural study of high albite, monalbite, and the analbite-->monalbite phase transition

A high-temperature study of the thermal expansion and the anisotropy of the sodium atom in low albite