Controls on Trace Element Distribution in Oxides and Silicates

ABSTRACT Understanding and quantifying the partitioning of elements at low concentrations is important for petrology, as minor and trace elements are used as tracers of many geological processes. The lattice strain model has consequently been found to be very useful, as it links mineral/melt partition coefficients to the elastic properties of the mineral and to the difference in ionic radius between the trace element and the cation it replaces. However, this model has limitations, particularly in terms of describing crystal strain, due to the form of the equation and from the choice of a hard-sphere model. After a review of the thermodynamics of trace element incorporation into crystal sites and equilibrium between mineral phases, we present classical atomistic modelling using transferable empirical potentials. Following incorporation of one (or more) mismatching element(s), crystal strain appears strongly related to the environment of the site of exchange, with anisotropic, vacancy-rich minerals deforming more than densely-packed minerals, due to anisotropic deformation, rotation and tilting of the surrounding polyhedra. As shown by the computed cation–oxygen length in strained crystal sites, bond strain is smaller than the difference in cation radii, and varies between structures, with densely-packed minerals being less strained. Consequently, computed relaxation energies are smaller in less compressible minerals. This has implications for modelling the partitioning of trace elements, and highlights the limits of using continuum mechanics below crystal cell scale. Neither oxides nor silicates have isotropic elastic properties, and at nanoscale they do not deform like continuous material. Predictions of partitioning between mineral and melt (or fluid) remain hampered by several factors, amongst which are (i) knowledge of the thermodynamics properties of the dissolved species (in melts and fluids); (ii) the precision of estimated defect and strain energies; and (iii) the presence of solid solution in natural systems, even when limited to a few percent. The latter may have orders-of-magnitude effects on the calculated partition coefficients due to interactions between strain fields around minor cations, leading to a chemical mixing regime at low concentration different from solid solutions where phase components are found in high proportions. Partitioning between mineral phases in an assemblage is described with similar equations as for mineral/liquid equil