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

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

Toward a mineral physics reference model for the Moon's core

International audienceIron is the main constituent of terrestrial planetary cores, taking on a hexagonal closed packed structure under the conditions of Earth’s inner core, and a face-centered cubic (fcc) structure at the more moderate pressures of smaller bodies, such as the Moon, Mercury, or Mars. Here we present sound velocity and density measurements of fcc iron at pressures and temperatures characteristic of small planetary interiors. The results indicate that the seismic velocities currently proposed for the Moon’s inner core are well below those of fcc iron or plausible iron alloys. Our dataset provides strong constraints to seismic models of the lunar core and cores of small telluric planets, and allows us to build a direct compositional and velocity model of the Moon’s core

Phase relations in the system Fe-Ni-Si to 200 GPa and 3900 K and implications for Earth’s core

Phase relations in Fe–5 wt%Ni–4 wt%Si alloy was examined in an internally resistive heated diamond anvil cell under high pressure (P) and temperature (T) conditions to about 200 GPa and 3900 K by in-situ synchrotron X-ray diffraction. The hexagonal close-packed (hcp) structure was observed to the highest P–T condition, supporting the idea that the stable iron alloy structure in Earth's inner core is hcp. The P–T locations of the phase transition between the face-centred cubic (fcc) and hcp structures were also constrained to 106 GPa. The transition occurs at 15 GPa and 1000 K similar to for pure Fe. The Clausius–Clapeyron slope is however, 0.0480 GPa/K which is larger than reported slopes for Fe (0.0394 GPa/K), Fe–9.7 wt%Ni (0.0426 GPa/K), and Fe–4 wt%Si (0.0394 GPa/K), stabilising the fcc structure towards high pressure. Thus the simultaneous addition of Ni and Si to Fe increases the dP/dT slope of the fcc–hcp transition. This is associated with a small volume change upon transition in Fe–Ni–Si. The triple point, where the fcc, hcp, and liquid phases coexist in Fe–5 wt%Ni–4 wt%Si is placed at 145 GPa and 3750 K. The resulting melting temperature of the hcp phase at the inner core-outer core boundary lies at 550 K lower than in pure Fe

Phase transition boundary between fcc and hcp structures in Fe-Si alloy and its implications for terrestrial planetary cores

International audienc