Thermodynamics of mixing in diopside-jadeite, CaMgSi2O6-NaAlSi2O6, solid solution from staticlattice energy calculations

Static lattice energy calculations (SLEC), based on empirical interatomic potentials, have beenperformed for a set of 800 different structures in a 2 2 4 supercell of C2/c diopside with compositionsbetween diopside and jadeite, and with different states of order of the exchangeable Na/Ca and Mg/Al cations. Excess static energies of these structures have been cluster expanded in a basis set of 37 pair-interaction parameters. These parameters have been used to constrain Monte Carlo simulations of temperature-dependent properties in the range of 273?2,023 K and to calculate a temperature?composition phase diagram. The simulations predict the order?disorder transition in omphacite at1,150 20C in good agreement with the experimental data of Carpenter (Mineral Petrol 78:433?440, 1981). The stronger ordering of Mg/Al within the M1 site than of Ca/Na in the M2 site is attributed to the shorter M1?M1 nearest-neighbor distance, and, consequently, the stronger ordering force. The comparison of the simulated relationship between the order parameters corresponding to M1 and M2 sites with the X-ray refinement data on natural omphacites (Boffa Ballaran et al. in Am Mineral83:419?433, 1998) suggests that the cation ordering becomes kinetically ineffective at about 600C

Atomistic model of diopside–K-jadeite (CaMgSi2O6–KAlSi2O6) solid solution

Atomistic model was proposed to describe the thermodynamics of mixing in the diopside–K-jadeite solid solution (CaMgSi2O6–KAlSi2O6). The simulations were based on minimization of the latticeenergies of 800 structures within a 2 × 2 × 4 supercell of C2/c diopside with the compositions betweenCaMgSi2O6 and KAlSi2O6 and with variable degrees of order/disorder in the arrangement of Ca/K cations in M2 site and Mg/Al in Ml site. The energy minimization was performed with the help of a force-field model. The results of the calculations were used to define a generalized Ising model, which included 37 pair interaction parameters. Isotherms of the enthalpy of mixing within the range of 273–2023 K were calculated with a Monte Carlo algorithm, while the Gibbs free energies of mixing were obtained by thermodynamic integration of the enthalpies of mixing. The calculated T–X diagram for the system CaMgSi2O6–KAlSi2O6 at temperatures below 1000 K shows several miscibility gaps, which are separated by intervals of stability of intermediate ordered compounds. At temperatures above 1000 K a homogeneous solid solution is formed. The standard thermodynamic properties of K-adeite (KAlSi2O6) evaluated from quantum mechanical calculations were used to determine location of several mineral reactions with the participation of the diopside–K-jadeite solid solution. The results of the simulations suggest that the low content of KalSi2O6 in natural clinopyroxenes is not related to crystal chemical factors preventing isomorphism, but is determined by relatively high standard enthalpy of this end member

Order/disorder phase transition in cordierite and its possible relationship to the development of symplectite reaction textures in granulites

Based on a consistent set of empirical interatomic potentials, static structure energy calculations of various Al/Si configurations in the supercell of Mg-cordierite and Monte Carlo simulations the phase transition between the orthorhombic and hexagonal modifications of cordierite (Crd) is predicted at 1623 K. The temperature dependences of the enthalpy, entropy, and free energy of the Al/Si disorder were calculated using the method of thermodynamic integration. The simulations suggest that the commonly observed crystallization of cordierite in the disordered hexagonal form could be related to a tendency of Al to occupy T1 site, which is driven by local charge balance. The increase in the Al fraction in the T1 site over the ratio of 2/3(T1): 1/3(T2), that characterizes the ordered state, precludes formation of the domains of the orthorhombic phase. This intrinsic tendency to the crystallization of the metastable hexagonal phase could have significantly postponed the formation of the association of orthorhombic cordierite and orthopyroxene over the association of quartz and garnet in metapelites subjected to granulite facies metamorphism. The textures of local metasomatic replacement (the formation of Crd + Opx or Spr + Crd symplectites between the grains of garnet and quartz) indicate the thermodynamic instability of the association of Qtz + Grt at the moment of the metasomatic reaction. This instability could have been caused by the difficulty of equilibrium nucleation of orthorhombic cordierite

Modelling in relation to cation ordering

We review the methodology of using computer models to obtain quantitative information about cation ordering. Empirical interactomic potentials or ab initio electronic structure calculations are used to generate the energies for many configurations containing disordered arrangements of cations, and the parameters in model Hamiltonians can be determined from these energies. Monte Carlo simulations are then used to generate ensemble averages as functions of temperature or chemical composition. Analysis of the Monte Carlo ensembles directly yields the temperature dependence of long-range and short-range order, and thermodynamic quantities such as energy and heat capacity. Use of thermodynamic integration allows for the calculation of entropy and free energy. The methods are illustrated by examples showing long-range order/disorder phase transitions (feldspars), short-range order in solid solutions (pyrope-grossular), and non-convergent ordering (magnesium aluminate spinel); where comparisons with experimental data are possible, the model calculations are seen to give results that are reasonably accurate. The example in which ab initio electronic structure calculations are used show that it is now possible to extract accurate thermodynamic data for ordering processes using models that require no prior experimental data

Computational study of tetrahedral Al-Si ordering in muscovite

The nature of Al-Si ordering across the tetrahedral sites in muscovite, computational techniques. Values of the atomic exchange interaction parameters J1 were obtained. From these parameters, a two-dimensional Al-Si ordering scheme was deduced. The transition temperature Tc for this two-dimensional ordering is 1900 K. These are several possible ordering schemes in three dimensions, based on different stacking sequences of ordered sheets of tetrahedral sites. Monte Carlo simulations of both two-dimensional and three-dimensional ordering were performed, but in the three-dimensional simulation only the two-dimensional ordering is seen, implying that three-dimensional ordering is too slow to be attained during the timescale of the simulation. The effect of the three-dimensional interactions is to raise the two-dimensional ordering temperature to 2140 K. From the three-dimensional Monte Carlo simulation, the frequency of occurrence of 4SiOAl, 3Si1Al, 2Si2Al and 1Si3Al clusters was determined, which match those inferred by 29Si MAS-NMR measurements reasonably well. In fact, the match suggests that the cation ordering seen in experiments corresponds to a configuration with considerable short-range order but no long-range order, similar to a state that is at a temperature just above an ordering phase transition

Computational methods for the study of energies of cation distributions: applications to cation-ordering phase transitions and solid solutions

The structural and thermodynamic properties of minerals are strongly affected by cation site-ordering processes. We describe methods to determine the main interatomic interactions that drive the ordering process, which are based on parameterizing model Hamiltonians using empirical interatomic potentials and/or ab initio quantum mechanics methods. The methods are illustrated by a number of case study examples, including Al/Si ordering in aluminosilicates, Mg/Ca ordering in garnets, simultaneous Al/Si and Mg/Al ordering in pyroxenes, micas and amphiboles, and Mg/Al non-convergent ordering in spinel using only quantum mechanical methods