The elastic properties and stability of fcc-Fe and fcc-FeNi alloys at inner-core conditions

The agreement between shear wave velocities for the Earth's inner core observed from seismology with those derived from mineral physics is considerably worse than for any other region of the Earth. Furthermore, there is still debate as to the phase of iron present in the inner core, particularly when alloying with nickel and light elements is taken into account. To investigate the extent to which the mismatch between seismology and mineral physics is a function of either crystal structure and/or the amount of nickel present, we have used ab initio molecular dynamics simulations to calculate the elastic constants and seismic velocities (V-p and V-s) of face centred cubic (fcc) iron at Earth's inner core pressures (360 GPa) and at temperatures up to similar to 7000 K. We find that V-p for fcc iron (fcc-Fe) is very similar to that for hexagonal close packed (hcp) iron at all temperatures. In contrast, V-s for fcc-Fe is significantly higher than in hcp-Fe, with the difference increasing with increasing temperature; the difference between V-s for the core (from seismology) and V-s for fcc-Fe exceeds 40 per cent. These results are consistent with previous work at lower temperatures. We have also investigated the effect of 6.5 and 13 atm% Ni in fcc-Fe. We find that Ni only slightly reduces V-p and V-s (e.g. by 2 per cent in V-s for 13 atm% Ni at 5500 K), and cannot account for the difference between the velocities observed in the core and those of pure fcc-Fe. We also tried to examine pre-melting behaviour in fcc-Fe, as reported in hcp-Fe by extending the study to very high temperatures (at which superheating may occur). However, we find that fcc-Fe spontaneously transforms to other hcp-like structures before melting; two hcp-like structures were found, both of hexagonal symmetry, which may most easily be regarded as being derived from an hcp crystal with stacking faults. That the structure did not transform to a true hcp phase is likely as a consequence of the limited size of the simulation box (108 atoms). At 360 GPa, in pure fcc-Fe, we find that the transition from fcc to the hcp-like structure occurs at 7000 K, whereas in the Ni bearing system, the transition occurs at higher temperature (7250 K). This reinforces previous work showing that fcc-Fe might transform to hcp-Fe just before melting, and that Ni tends to stabilize the fcc structure with respect to hcp.