The issue of melting of pure iron and iron alloyed with lighter elements at high pressure is critical to the physics of the Earth. The iron melting curve in the relevant pressure range between 3 and 4 Mbar is reasonably well established from the theoretical point of view. However, so far no one attempted a direct atomistic simulation of iron alloyed with light elements. We investigate here the impact of alloying the body-centered cubic (bcc) Fe with Si. We simulate melting of the bcc Fe and ${\text{Fe}}_{0.9375}{\text{Si}}_{0.0625}$ alloy by ab initio molecular dynamics. The addition of light elements to the hexagonal-close-packed (hcp) iron is known to depress its melting temperature $\left(T_m\right)$. We obtain, in marked contrast, that alloying of bcc Fe with Si does not lead to $T_m$ depression; on the contrary, the $T_m$ slightly increases. This suggests that if Si is a typical impurity in the Earth’s inner core, then the stable phase in the core is bcc rather than hcp.

• Received 28 May 2009

DOI:https://doi.org/10.1103/PhysRevB.79.220102