Effects of electron correlations on transport properties of iron at Earth’s core conditions

Electronic correlations in dense iron: from moderate pressure to Earth's core conditions

We discuss the role of dynamical many-electron effects in the physics of iron and iron-rich solid alloys under applied pressure on the basis of recent ab initio studies employing the dynamical mean-field theory (DMFT). We review in detail two particularly interesting regimes: first, a moderate pressure range up to 60 GPa and, second, the ultra-high pressure of about 360 GPa expected inside the solid inner core of Earth. Electronic correlations in iron under the moderate pressure of several tens GPa are discussed in the first section. DMFT-based methods predict an enhancement of electronic correlations at the pressure-induced body-centered cubic α to hexagonal close-packed [Formula: see text] phase transition. In particular, the electronic effective mass, scattering rate and electron-electron contribution to the electrical resistivity undergo a step-wise increase at the transition point. One also finds a significant many-body correction to the [Formula: see text]-Fe equation of state, thus clarifying the origin of discrepancies between previous DFT studies and experiment. An electronic topological transition is predicted to be induced in [Formula: see text]-Fe by many-electron effects; its experimental signatures are analyzed. The next section focuses on the geophysically relevant pressure-temperature regime of the Earth's inner core (EIC) corresponding to the extreme pressure of 360 GPa combined with temperatures up to 6000 K. The three iron allotropes ([Formula: see text], [Formula: see text] and face-centered-cubic [Formula: see text]) previously proposed as possible stable phases at such conditions are found to exhibit qualitatively different many-electron effects as evidenced by a strongly non-Fermi-liquid metallic state of [Formula: see text]-Fe and an almost perfect Fermi liquid in the case of [Formula: see text]-Fe. A recent active discussion on the electronic state and transport properties of [Formula: see text]-Fe at the EIC conditions is reviewed in details. Estimations for the dynamical many-electron contribution to the relative phase stability are presented. We also discuss the impact of a Ni admixture, which is expected to be present in the core matter. We conclude by outlining some limitation of the present DMFT-based framework relevant for studies of iron-base systems as well as perspective directions for further development.

Thermal Conductivity and Electrical Resistivity of Solid Iron at Earth's Core Conditions from First Principles

We compute the thermal conductivity and electrical resistivity of solid hcp Fe to pressures and temperatures of Earth's core. We find significant contributions from electron-electron scattering, usually neglected at high temperatures in transition metals. Our calculations show a quasilinear relation between the electrical resistivity and temperature for hcp Fe at extreme high pressures. We obtain thermal and electrical conductivities that are consistent with experiments considering reasonable error. The predicted thermal conductivity is reduced from previous estimates that neglect electron-electron scattering. Our estimated thermal conductivity for the outer core is 77±10  W m^{-1} K^{-1} and is consistent with a geodynamo driven by thermal convection.

Size-dependent electrical transport properties in Co nanocluster-assembled granular films

A series of Co nanocluster-assembled films with cluster sizes ranging from 4.5 nm to 14.7 nm were prepared by the plasma-gas-condensation method. The size-dependent electrical transport properties were systematically investigated. Both of the longitudinal resistivity (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{xx}$$\end{document}ρxx) and saturated anomalous Hall resistivity (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{xy}^{A}$$\end{document}ρxyA) continuously increased with the decrease of the cluster sizes (d). The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{xx}$$\end{document}ρxx firstly increased and then decreased with increasing the temperature for all samples, which could be well described by involving the thermally fluctuation-induced tunneling (FIT) process and scattering. The tunneling effect was verified to result in the invalidation of classical anomalous Hall effect (AHE) scaling relation. After deducting the contribution from tunneling effect to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{xx}$$\end{document}ρxx, the AHE scaling relation between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{xy}^{A}$$\end{document}ρxyA and the scattering resistivity (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{S}$$\end{document}ρS) by varying the temperature was reconstructed. The value of scaling exponent γ increased with increasing Co cluster sizes. The size dependence of γ might be qualitatively interpreted by the interface and surface-induced spin flip scattering. We also determined the scaling relation between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{xy}^{A}$$\end{document}ρxyA and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{S}$$\end{document}ρS at 5 K by changing the Co cluster sizes, and a large value of γ = 3.6 was obtained which might be ascribed to the surface and interfacial scattering.

Powering Earth's dynamo with magnesium precipitation from the core

Earth's global magnetic field arises from vigorous convection within the liquid outer core. Palaeomagnetic evidence reveals that the geodynamo has operated for at least 3.4 billion years, which places constraints on Earth's formation and evolution. Available power sources in standard models include compositional convection (driven by the solidifying inner core's expulsion of light elements), thermal convection (from slow cooling), and perhaps heat from the decay of radioactive isotopes. However, recent first-principles calculations and diamond-anvil cell experiments indicate that the thermal conductivity of iron is two or three times larger than typically assumed in these models. This presents a problem: a large increase in the conductive heat flux along the adiabat (due to the higher conductivity of iron) implies that the inner core is young (less than one billion years old), but thermal convection and radiogenic heating alone may not have been able to sustain the geodynamo during earlier epochs. Here we show that the precipitation of magnesium-bearing minerals from the core could have served as an alternative power source. Equilibration at high temperatures in the aftermath of giant impacts allows a small amount of magnesium (one or two weight per cent) to partition into the core while still producing the observed abundances of siderophile elements in the mantle and avoiding an excess of silicon and oxygen in the core. The transport of magnesium as oxide or silicate from the cooling core to underneath the mantle is an order of magnitude more efficient per unit mass as a source of buoyancy than inner-core growth. We therefore conclude that Earth's dynamo would survive throughout geologic time (from at least 3.4 billion years ago to the present) even if core radiogenic heating were minimal and core cooling were slow.

Palaeomagnetic field intensity variations suggest Mesoproterozoic inner-core nucleation

The Earth's inner core grows by the freezing of liquid iron at its surface. The point in history at which this process initiated marks a step-change in the thermal evolution of the planet. Recent computational and experimental studies have presented radically differing estimates of the thermal conductivity of the Earth's core, resulting in estimates of the timing of inner-core nucleation ranging from less than half a billion to nearly two billion years ago. Recent inner-core nucleation (high thermal conductivity) requires high outer-core temperatures in the early Earth that complicate models of thermal evolution. The nucleation of the core leads to a different convective regime and potentially different magnetic field structures that produce an observable signal in the palaeomagnetic record and allow the date of inner-core nucleation to be estimated directly. Previous studies searching for this signature have been hampered by the paucity of palaeomagnetic intensity measurements, by the lack of an effective means of assessing their reliability, and by shorter-timescale geomagnetic variations. Here we examine results from an expanded Precambrian database of palaeomagnetic intensity measurements selected using a new set of reliability criteria. Our analysis provides intensity-based support for the dominant dipolarity of the time-averaged Precambrian field, a crucial requirement for palaeomagnetic reconstructions of continents. We also present firm evidence for the existence of very long-term variations in geomagnetic strength. The most prominent and robust transition in the record is an increase in both average field strength and variability that is observed to occur between a billion and 1.5 billion years ago. This observation is most readily explained by the nucleation of the inner core occurring during this interval; the timing would tend to favour a modest value of core thermal conductivity and supports a simple thermal evolution model for the Earth.

Saturation of electrical resistivity of solid iron at Earth's core conditions

We report on the temperature dependence of the electrical resistivity of solid iron at high pressure, up to and including conditions likely to be found at the centre of the Earth. We have extended some of the calculations of the resistivities of pure solid iron we recently performed at Earth's core conditions (Pozzo et al. in Earth Planet Sci Lett 393:159-164, 2014) to lower temperature. We show that at low temperature the resistivity increases linearly with temperature, and saturates at high temperature. This saturation effect is well known as the Mott-Ioffe-Regel limit in metals, but has been largely ignored to estimate the resistivity of iron at Earth's core conditions. Recent experiments (Gomi et al. in Phys Earth Planet Int 224:88-103, 2013) coupled new high pressure data and saturation to predict the resitivity of iron and iron alloys at Earth's core conditions, and reported values up to three times lower than previous estimates, confirming recent first principles calculations (de Koker et al. in Proc Natl Acad Sci 109:4070-4073, 2012; Pozzo et al. in Nature 485:355-358, 2012, Phys Rev B 87:014110-10, 2013, Earth Planet Sci Lett 393:159-164, 2014; Davies et al. in Nat Geosci 8:678-685, 2015). The present results support the saturation effect idea.

PALEOMAGNETISM. A Hadean to Paleoarchean geodynamo recorded by single zircon crystals

Knowing when the geodynamo started is important for understanding the evolution of the core, the atmosphere, and life on Earth. We report full-vector paleointensity measurements of Archean to Hadean zircons bearing magnetic inclusions from the Jack Hills conglomerate (Western Australia) to reconstruct the early geodynamo history. Data from zircons between 3.3 billion and 4.2 billion years old record magnetic fields varying between 1.0 and 0.12 times recent equatorial field strengths. A Hadean geomagnetic field requires a core-mantle heat flow exceeding the adiabatic value and is suggestive of plate tectonics and/or advective magmatic heat transport. The existence of a terrestrial magnetic field before the Late Heavy Bombardment is supported by terrestrial nitrogen isotopic evidence and implies that early atmospheric evolution on both Earth and Mars was regulated by dynamo behavior.

The electrical conductivity of Al2O3 under shock-compression

Sapphire (Al2O3) crystals are used below 100 GPa as anvils and windows in dynamic-compression experiments because of their transparency and high density. Above 100 GPa shock pressures, sapphire becomes opaque and electrically conducting because of shock-induced defects. Such effects prevent temperature and dc conductivity measurements of materials compressed quasi-isentropically. Opacities and electrical conductivities at ~100 GPa are non-equilibrium, rather than thermodynamic parameters. We have performed electronic structure calculations as a guide in predicting and interpreting shock experiments and possibly to discover a window up to ~200 GPa. Our calculations indicate shocked sapphire does not metallize by band overlap at ~300 GPa, as suggested previously by measured non-equilibrium data. Shock-compressed Al2O3 melts to a metallic liquid at ~500 GPa and 10,000 K and its conductivity increases rapidly to ~2000 Ω(-1)cm(-1) at ~900 GPa. At these high shock temperatures and pressures sapphire is in thermal equilibrium. Calculated conductivity of Al2O3 is similar to those measured for metallic fluid H, N, O, Rb, and Cs. Despite different materials, pressures and temperatures, and compression techniques, both experimental and theoretical, conductivities of all these poor metals reach a common end state typical of strong-scattering disordered materials.