We employ state-of-the-art ab initio simulations within the dynamical mean-field theory to study three likely phases of iron (hexogonal close-packed, hcp, face centered cubic, fcc, and body centered cubic, bcc) at the Earths core conditions. We demon
strate that the correction to the electronic free energy due to correlations can be significant for the relative stability of the phases. The strongest effect is observed in bcc Fe, which shows a non-Fermi liquid behaviour, and where a Curie-Weiss behaviour of the uniform susceptbility hints at a local magnetic moment still existing at 5800 K and 300 GPa. We predict that all three structures have sufficiently high magnetic susceptibility to stabilize the geodynamo.
An accurate and easily extendable method to deal with lattice dynamics of solids is offered. It is based on first-principles molecular dynamics simulations and provides a consistent way to extract the best possible harmonic - or higher order - potent
ial energy surface at finite temperatures. It is designed to work even for strongly anharmonic systems where the traditional quasiharmonic approximation fails. The accuracy and convergence of the method are controlled in a straightforward way. Excellent agreement of the calculated phonon dispersion relations at finite temperature with experimental results for bcc Li and bcc Zr is demonstrated.
