Iron isotope fractionation between metal and silicate during core-mantle differentiation in rocky bodies

Abstract

Fe isotope variations in rocky bodies reveal fundamental information about planetary evolution. However, experimental results have come to contradictory conclusions on the equilibrium Fe isotope fractionation between metal and silicate during core-mantle differentiation. Many different processes, including evaporation, core formation, partial melting and disproportion of mantle silicate, have been consequently proposed to explain the observed Fe isotope variations in rocky solar system bodies. Here we perform ab initio molecular dynamics simulations and find that the anharmonicity in iron strongly decreases the force constant of Fe at low pressures (<∼50 GPa), which even reverses the equilibrium Fe isotope fractionation between metal and silicate. We conclude that pyrolitic melt is always enriched in heavy Fe isotopes relative to liquid Fe-alloys, no matter what pressure. Therefore core-mantle differentiation will play a significant role in explaining the heavy Fe isotope compositions of the mantles of some rocky bodies (e.g., Earth, the ureilite parent body, and possibly the asteroid Vesta). As all previously proposed processes for Fe isotope fractionation can only enrich the mantle-derived rocks in heavy Fe isotopes, the near/sub-chondritic Fe isotope signatures of Mars and the aubrite parent body thus imply that iron sulfide enriched in light Fe isotopes may significantly contribute to the iron components of those meteoritic samples.