Compression experiments to 126 GPa and 2500 K and thermal equation of state of Fe3S: Implications for sulphur in the Earth’s core

Pressure-volume-temperature (P-V-T) experiments on tetragonal Fe3S were conducted to 126 GPa and 2500 K in laser-heated diamond anvil cells (DAC) with in-situ X-ray diffraction (XRD). Seventy nine high-T data as well as four 300-K data were collected, based on which new thermal equations of state (EoS) for Fe3S were established. The room-T data together with existing data were fitted to the third order Birch-Murnaghan EoS, which yielded, GPa and with fixed at 377.0 Å3. A constant term in the thermal pressure equation, Pth = , fitted the high-T data well to the highest temperature, which implies that the contributions from the anharmonic and electronic terms should be minor in the thermal pressure term. The high-T data were also fitted to the Mie-Grüneisen-Debye model; with and q fixed at 417 K and 1 respectively. Calculations from the EoS show that crystalline Fe3S at 4000-5500 K is denser than the Earth's outer core and less dense than the inner core. Assuming a density reduction due to melting, liquid Fe3S would meet the outer core density profile, which however suggests that no less than 16 wt%S is needed to reconcile the observed outer core density deficit. The S-rich B2 phase, which was suggested to be a potential liquidus phase of an Fe3S-outer core above 250 GPa, namely the main constituent of its solid inner core, would likely be less dense than the Earth's inner core. As such, while the outer core density requires as much sulphur as 16 wt%, the resulting liquidus phase cannot meet the density of the inner core. Any sulphur-rich composition should therefore be rejected for the Earth's core

High-pressure melting experiments of Fe3S and a thermodynamic model of the Fe-S liquids for the Earth’s core

Melting experiments on Fe3S were conducted to 75 GPa and 2800 K in laser-heated and internally resistive-heated diamond anvil cells with in-situ x-ray diffraction and/or post-mortem textural observation. From the constrained melting curve, we assessed the thermal equation of state for Fe3S liquid. Then we constructed a thermodynamic model of melting of the system Fe–Fe3S including the eutectic relation under high pressures based on our new experimental data. The mixing properties of Fe–S liquids under high pressures were evaluated in order to account for existing experimental data on eutectic temperature. The results demonstrate that the mixing of Fe and S liquids are nonideal at any core pressure. The calculated sulphur content in eutectic point decreases with increasing pressure to 120 GPa and is fairly constant of 8 wt% at greater pressures. From the Gibbs free energy, we derived the parameters to calculate the crystallising point of an Fe–S core and its isentrope, and then we calculated the density and the longitudinal seismic wave velocity (Vp) of these liquids along each isentrope. While Fe3S liquid can account for the seismologically constrained density and Vp profiles over the outer core, the density of the precipitating phase is too low for the inner core. On the other hand, a hypothetical Fe–S liquid core with a bulk composition on the Fe-rich side of the eutectic point cannot represent the density and Vp profiles of the Earth's outer core. Therefore, Earth's core cannot be approximated by the system Fe–S and it should include another light element