An ab initio study of the relative stabilities and equations of state of Fe3S polymorphs

Abstract

AbstractAn investigation of the relative stabilities and equations of state of possible Fe3S polymorphs was conducted using first-principles pseudopotential calculations. These calculations were based on density functional theory and performed using ultrasoft Vanderbilt pseudopotentials within the generalized gradient approximation. In accord with experiment, we found that the tetragonal Fe3P-type polymorph is the only stable phase along the 0 K isotherm as a function of pressure. Fe3S exhibits permanent magnetism at ambient conditions (Fei et al., 2000), but magnetism is suppressed by pressure and temperature, and therefore non-magnetic data are appropriate ones to use for modelling planetary interiors. For this reason, and because the Fe3P-type polymorph of Fe3S contains 32 atoms per unit cell it was impractical to incorporate magnetic properties into the simulations of this phase, we studied the behaviour of the non-magnetic phase. We obtained values of 250 GPa for the bulk modulus, K0, and 4.61 for its first derivative withrespect to pressure, K0′, by fitting a 3rd order Birch-Murnaghan equation of state to the calculated internal energy as a function of volume for the non-magnetic Fe3P-type Fe3S. This suggests that a pressure far greater than that expected in the Martian interior would be needed to achieve a density comparable to that of the Martian core. We therefore conclude that it is unlikely that the core of Mars contains significant amounts of solid Fe3S.