No Arabic abstract
The Earth acts as a gigantic heat engine driven by decay of radiogenic isotopes and slow cooling, which gives rise to plate tectonics, volcanoes, and mountain building. Another key product is the geomagnetic field, generated in the liquid iron core by a dynamo running on heat released by cooling and freezing to grow the solid inner core, and on chemical convection due to light elements expelled from the liquid on freezing. The power supplied to the geodynamo, measured by the heat-flux across the core-mantle boundary (CMB), places constraints on Earths evolution. Estimates of CMB heat-flux depend on properties of iron mixtures under the extreme pressure and temperature conditions in the core, most critically on the thermal and electrical conductivities. These quantities remain poorly known because of inherent difficulties in experimentation and theory. Here we use density functional theory to compute these conductivities in liquid iron mixtures at core conditions from first principles- the first directly computed values that do not rely on estimates based on extrapolations. The mixtures of Fe, O, S, and Si are taken from earlier work and fit the seismologically-determined core density and inner-core boundary density jump. We find both conductivities to be 2-3 times higher than estimates in current use. The changes are so large that core thermal histories and power requirements must be reassessed. New estimates of adiabatic heat-flux give 15-16 TW at the CMB, higher than present estimates of CMB heat-flux based on mantle convection; the top of the core must be thermally stratified and any convection in the upper core driven by chemical convection against the adverse thermal buoyancy or lateral variations in CMB heat flow. Power for the geodynamo is greatly restricted and future models of mantle evolution must incorporate a high CMB heat-flux and explain recent formation of the inner core.
The electronic state and transport properties of hot dense iron are of the utmost importance to geophysics. Combining the density functional and dynamical mean field theories we study the impact of electron correlations on electrical and thermal resistivity of hexagonal close-packed $epsilon$-Fe at Earths core conditions. $epsilon$-Fe is found to behave as a nearly perfect Fermi liquid. The quadratic dependence of the scattering rate in Fermi liquids leads to a modification of the Wiedemann-Franz law with suppression of the thermal conductivity as compared to the electrical one. This significantly increases the electron-electron thermal resistivity which is found to be of comparable magnitude to the electron-phonon one. The implications of this effect on the dynamics of Earths core is discussed.
The transport properties of iron under Earths inner core conditions are essential input for the geophysical modelling but are poorly constrained experimentally. Here we show that the thermal and electrical conductivity of iron at those conditions remains high even if the electron-electron-scattering (EES) is properly taken into account. This result is obtained by ab initio simulations taking into account consistently both thermal disorder and electronic correlations. Thermal disorder suppresses the non-Fermi-liquid behavior of the body-centered cubic iron phase, hence, reducing the EES; the total calculated thermal conductivity of this phase is 220 Wm$^{-1}$K$^{-1}$ with the EES reduction not exceeding 20%. The EES and electron-lattice scattering are intertwined resulting in breaking of the Matthiessens rule with increasing EES. In the hexagonal close-packed iron the EES is also not increased by thermal disorder and remains weak. Our main finding thus holds for the both likely iron phases in the inner core.
The crystal structure of iron in the Earths inner core remains debated. Most recent experiments suggest a hexagonal-close-packed (hcp) phase. In simulations, it has been generally agreed that the hcp-Fe is stable at inner core pressures and relatively low temperatures. At high temperatures, however, several studies suggest a body-centered-cubic (bcc) phase at the inner core condition. We have examined the crystal structure of iron at high pressures over 2 million atmospheres (>200GPa) and at high temperatures over 5000 kelvin in a laser-heated diamond cell using microstructure analysis combined with $textit{in-situ}$ x-ray diffraction. Experimental evidence shows a bcc-Fe appearing at core pressures and high temperatures, with an hcp-bcc transition line in pressure-temperature space from about 95$pm$2GPa and 2986$pm$79K to at least 222$pm$6GPa and 4192$pm$104K. The trend of the stability field implies a stable bcc-Fe at the Earths inner core condition, with implications including a strong candidate for explaining the seismic anisotropy of the Earths inner core.
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 demonstrate 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.
Some Bravais lattices have a particular geometry that can slow down the motion of Bloch electrons by pre-localization due to the band-structure properties. Another known source of electronic localization in solids is the Coulomb repulsion in partially filled d- or f-orbitals, which leads to the formation of local magnetic moments. The combination of these two effects is usually considered of little relevance to strongly correlated materials. Here we show that it represents, instead, the underlying physical mechanism in two of the most important ferromagnets: nickel and iron. In nickel, the van Hove singularity has an unexpected impact on the magnetism. As a result, the electron-electron scattering rate is linear in temperature, in violation of the conventional Landau theory of metals. This is true even at Earths core pressures, at which iron is instead a good Fermi liquid. The importance of nickel in models of geomagnetism may have therefore to be reconsidered.