Viscosity of hcp iron at Earth's inner core conditions from density functional theory

The inner core, extending to 1,221 km above the Earth's center at pressures between 329 and 364 GPa, is primarily composed of solid iron. Its rheological properties influence both the Earth's rotation and deformation of the inner core which is a potential source of the observed seismic anisotropy. However, the rheology of the inner core is poorly understood. We propose a mineral physics approach based on the density functional theory to infer the viscosity of hexagonal close packed (hcp) iron at the inner core pressure (P) and temperature (T). As plastic deformation is rate-limited by atomic diffusion under the extreme conditions of the Earth's center, we quantify self-diffusion in iron non-empirically. The results are applied to model steady-state creep of hcp iron. Here, we show that dislocation creep is a key mechanism driving deformation of hcp iron at inner core conditions. The associated viscosity agrees well with the estimates from geophysical observations supporting that the inner core is significantly less viscous than the Earth's mantle. Such low viscosity rules out inner core translation, with melting on one side and solidification on the opposite, but allows for the occurrence of the seismically observed fluctuations in inner core differential rotation.