The elastic properties of hcp-Fe alloys under the conditions of the Earth's inner core

Geophysical and cosmochemical constraints suggest the inner-core is mainly composed of iron with a few percent of light elements. However, despite extensive studies over many years, no single alloying light-element has been found that is able to simultaneously match the observed inner-core density and both seismic velocities. This has motivated a number of suggestions of other mechanism to lower velocities, such as anelasticity or premelting. However, an unexplored possibility is that a combination of two or more light-elements might produce the desired reduction in velocities and densities of the inner core. In order to test this, we use ab initio molecular dynamics calculations to map the elastic property space of hcp-Fe alloyed with S, Si and C at 360 GPa up to the melting temperature. Based on a mixing solid solution model together with direct simulations on the ternaries, we found a number of compositions which are able to match the observed properties of the inner core. This is the first time that the density, VP, Vs and the Poisson's ratio of the inner core have been matched directly with an hcp-Fe alloy.