skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Thermal conductivity of Fe-Si alloys and thermal stratification in Earth’s core
Light elements in Earth’s core play a key role in driving convection and influencing geodynamics, both of which are crucial to the geodynamo. However, the thermal transport properties of iron alloys at high-pressure and -temperature conditions remain uncertain. Here we investigate the transport properties of solid hexagonal close-packed and liquid Fe-Si alloys with 4.3 and 9.0 wt % Si at high pressure and temperature using laser-heated diamond anvil cell experiments and first-principles molecular dynamics and dynamical mean field theory calculations. In contrast to the case of Fe, Si impurity scattering gradually dominates the total scattering in Fe-Si alloys with increasing Si concentration, leading to temperature independence of the resistivity and less electron–electron contribution to the conductivity in Fe-9Si. Our results show a thermal conductivity of ∼100 to 110 W⋅m −1 ⋅K −1 for liquid Fe-9Si near the topmost outer core. If Earth’s core consists of a large amount of silicon (e.g., > 4.3 wt %) with such a high thermal conductivity, a subadiabatic heat flow across the core–mantle boundary is likely, leaving a 400- to 500-km-deep thermally stratified layer below the core–mantle boundary, and challenges proposed thermal convection in Fe-Si liquid outer core.  more » « less
Award ID(s):
1901801
PAR ID:
10315202
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
119
Issue:
1
ISSN:
0027-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. unknown (Ed.)
    The Earth’s core–mantle boundary presents a dramatic change in materials, from silicate to metal. While little is known about chemical interactions between them, a thin layer with a lower velocity has been proposed at the topmost outer core (Eʹ layer) that is difficult to explain with a change in concentration of a single light element. Here we perform high-temperature and -pressure laser-heated diamond-anvil cell experiments and report the formation of SiO2 and FeHx from a reaction between water from hydrous minerals and Fe–Si alloys at the pressure–temperature conditions relevant to the Earth’s core–mantle boundary. We suggest that, if water has been delivered to the core–mantle boundary by subduction, this reaction could enable exchange of hydrogen and silicon between the mantle and the core. The resulting H-rich, Si-deficient layer formed at the topmost core would have a lower density, stabilizing chemical stratification at the top of the core, and a lower velocity. We suggest that such chemical exchange between the core and mantle over gigayears of deep transport of water may have contributed to the formation of the putative Eʹ layer. 
    more » « less
  2. Seismic observations suggest that the uppermost region of Earth’s liquid outer core is buoyant, with slower velocities than the bulk outer core. One possible mechanism for the formation of a stably stratified layer is immiscibility in molten iron alloy systems, which has yet to be demonstrated at core pressures. We find immisci- bility between liquid Fe-Si and Fe-Si-O persisting to at least 140 GPa through a combination of laser-heated diamond-anvil cell experiments and first-principles molecular dynamics simulations. High-pressure immiscibility in the Fe-Si-O system may explain a stratified layer atop the outer core, complicate differentiation and evolution of the deep Earth, and affect the structure and intensity of Earth’s magnetic field. Our results support silicon and oxy- gen as coexisting light elements in the core and suggest that SiO2 does not crystallize out of molten Fe-Si-O at the core-mantle boundary. 
    more » « less
  3. Abstract The motion of liquid iron (Fe) alloy materials in the outer core drives the dynamo, which generates Mercury's magnetic field. The assessment of core models requires laboratory measurements of the melting temperature of Fe alloys at high pressure. Here, we experimentally determined the melting curve of Fe9wt%Si and Fe17wt%Si up to 17 GPa using in situ and ex situ measurements of intermetallic fast diffusion that serves as the melting criterion in a large‐volume press. Our determined melting slopes are comparable with previous studies up to about 17 GPa. However, when extrapolated, our melting slopes significantly deviate from previous studies at higher pressures. For Mercury's core with a model composition of Fe9wt%Si, the melting temperature‐depth profile determined in our study is lower by ∼150–250 K when compared with theoretical calculations. Using the new melting curve of Fe9wt%Si and the electrical resistivity values from a previous study of Fe8.5wt%Si, we estimate that the electronic thermal conductivity of liquid Fe9wt%Si is 30 Wm−1K−1at the Mercury'sCMBpressure of 5 GPa and 37 Wm−1K−1at an assumedICBof 21 GPa, corresponding to heat flux values of 23 mWm−2and 32 mWm−2, respectively. These values provide new constraints on the core models. 
    more » « less
  4. Abstract Electrical resistivity experiments were conducted on three alloys in the iron-rich side of the Fe-Ni(-S) system (Fe-5 wt% Ni, Fe-10 wt% Ni, Fe-10 wt% Ni-5 wt% S) at 4.5 and 8 GPa and up to 1900 K using the multi-anvil apparatus and the 4-electrode technique. For all samples, increasing temperature increases resistivity. At a specified temperature, Fe-Ni(-S) alloys are more resistive than Fe by a factor of about 3. Fe-Ni alloys containing 5 and 10 wt% Ni present comparable electrical resistivity values. The resistivity of Fe-Ni(-S) alloys is comparable to the one of Fe = 5 wt% S at 4.5 GPa and is about three times higher than the resistivity of Fe = 5 wt% S at 8 GPa, due to a different pressure dependence of electrical resistivity between Fe-Ni and Fe-S alloys. Based on these electrical results and experimentally determined thermal conductivity values from the literature, lower and upper bounds of thermal conductivity were calculated. For all Ni-bearing alloys, thermal conductivity estimates range between ~12 and 20 W/(m⋅K) over the considered pressure and temperature ranges. Adiabatic heat fluxes were computed for both Ganymede's core and the Lunar core, and heat flux values suggest a significant dependence to both core composition and the adiabatic temperature. Comparison with previous thermochemical models of the cores of Ganymede and the Moon suggests that some studies may have overestimated the thermal conductivity and hence, the heat flux along the adiabat in these planetary cores. 
    more » « less
  5. Abstract We present a joint experimental‐modeling investigation of core cooling in small terrestrial bodies. Significant amounts of light elements (S, O, Mg, Si) may compose the metallic cores of terrestrial planets and moons. However, the effect of multiple light elements on transport properties, in particular, electrical resistivity and thermal conductivity, is not well constrained. Electrical experiments were conducted at 10 GPa and up to 1850 K on high‐purity powder mixtures in the Fe‐S‐O(±Mg, ±Si) systems using the multianvil apparatus and the four‐electrode technique. The sample compositions contained 5 wt.% S, up to 3 wt.% O, up to 2 wt.% Mg, and up to 1 wt.% Si. We observe that above the eutectic temperature, electrical resistivity is significantly sensitive to the nature and amount of light elements. For each composition, thermal conductivity‐temperature equations were estimated using the experimental electrical results and a modified Wiedemann‐Franz law. These equations were implemented in a thermochemical core cooling model to study the evolution of the dynamo. Modeling results suggest that bulk chemistry significantly affects the entropy available to power dynamo action during core cooling. In the case of Mars, the presence of oxygen would delay the dynamo cessation by up to 1 Gyr compared to an O‐free, Fe‐S core. Models with 3 wt% O can be reconciled with the inferred cessation time of the Martian dynamo if the core‐mantle boundary heat flow falls from >2 TW to ~0.1 TW in the first 0.5 Gyr following core formation. 
    more » « less