skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth's solid inner core: THERMAL EQUATION OF STATE OF HCP-IRON

Abstract

We conducted high-pressure experiments on hexagonal close packed iron (hcp-Fe) in MgO, NaCl, and Ne pressure-transmitting media and found general agreement among the experimental data at 300 K that yield the best fitted values of the bulk modulus K 0 = 172.7(±1.4) GPa and its pressure derivative K 0'= 4.79(±0.05) for hcp-Fe, using the third-order Birch-Murnaghan equation of state. Using the derived thermal pressures for hcp-Fe up to 100 GPa and 1800 K and previous shockwave Hugoniot data, we developed a thermal equation of state of hcp-Fe. The thermal equation of state of hcp-Fe is further used to calculate the densities of iron along adiabatic geotherms to define the density deficit of the inner core, which serves as the basis for developing quantitative composition models of the Earth's inner core. We determine the density deficit at the inner core boundary to be 3.6%, assuming an inner core boundary temperature of 6000 K.

Authors:
 [1];  [1];  [2];  [1];  [3]
  1. Geophysical Laboratory, Carnegie Institution of Washington, Washington District of Columbia USA
  2. Geophysical Laboratory, Carnegie Institution of Washington, Washington District of Columbia USA; Now at Frontier Research Institute for Interdisciplinary Sciences, Tohoku University, Sendai Japan
  3. School of Sciences, Wuhan University of Technology, Wuhan China
Publication Date:
Research Org.:
Argonne National Lab. (ANL), Argonne, IL (United States). Advanced Photon Source (APS)
Sponsoring Org.:
NSFOTHER
OSTI Identifier:
1298276
Resource Type:
Journal Article
Resource Relation:
Journal Name: Geophysical Research Letters; Journal Volume: 43; Journal Issue: 13
Country of Publication:
United States
Language:
ENGLISH
Subject:
58 GEOSCIENCES

Citation Formats

Fei, Yingwei, Murphy, Caitlin, Shibazaki, Yuki, Shahar, Anat, and Huang, Haijun. Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth's solid inner core: THERMAL EQUATION OF STATE OF HCP-IRON. United States: N. p., 2016. Web. doi:10.1002/2016GL069456.
Fei, Yingwei, Murphy, Caitlin, Shibazaki, Yuki, Shahar, Anat, & Huang, Haijun. Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth's solid inner core: THERMAL EQUATION OF STATE OF HCP-IRON. United States. doi:10.1002/2016GL069456.
Fei, Yingwei, Murphy, Caitlin, Shibazaki, Yuki, Shahar, Anat, and Huang, Haijun. 2016. "Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth's solid inner core: THERMAL EQUATION OF STATE OF HCP-IRON". United States. doi:10.1002/2016GL069456.
@article{osti_1298276,
title = {Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth's solid inner core: THERMAL EQUATION OF STATE OF HCP-IRON},
author = {Fei, Yingwei and Murphy, Caitlin and Shibazaki, Yuki and Shahar, Anat and Huang, Haijun},
abstractNote = {We conducted high-pressure experiments on hexagonal close packed iron (hcp-Fe) in MgO, NaCl, and Ne pressure-transmitting media and found general agreement among the experimental data at 300 K that yield the best fitted values of the bulk modulus K0 = 172.7(±1.4) GPa and its pressure derivative K0'= 4.79(±0.05) for hcp-Fe, using the third-order Birch-Murnaghan equation of state. Using the derived thermal pressures for hcp-Fe up to 100 GPa and 1800 K and previous shockwave Hugoniot data, we developed a thermal equation of state of hcp-Fe. The thermal equation of state of hcp-Fe is further used to calculate the densities of iron along adiabatic geotherms to define the density deficit of the inner core, which serves as the basis for developing quantitative composition models of the Earth's inner core. We determine the density deficit at the inner core boundary to be 3.6%, assuming an inner core boundary temperature of 6000 K.},
doi = {10.1002/2016GL069456},
journal = {Geophysical Research Letters},
number = 13,
volume = 43,
place = {United States},
year = 2016,
month = 7
}
  • No abstract prepared.
  • The outer core of the Earth contains several weight percent of one or more unknown light elements, which may include silicon. Therefore it is critical to understand the high pressure–temperature properties and behavior of an iron–silicon alloy with a geophysically relevant composition (16 wt% silicon). We experimentally determined the melting curve, subsolidus phase diagram, and equations of state of all phases of Fe–16 wt%Si to 140 GPa, finding a conversion from the D0 3 crystal structure to a B2+hcp mixture at high pressures. The melting curve implies that 3520 K is a minimum temperature for the Earth's outer core, ifmore » it consists solely of Fe–Si alloy, and that the eutectic composition in the Fe–Si system is less than 16 wt% silicon at core–mantle boundary conditions. Comparing our new equation of state to that of iron and the density of the core, we find that for an Fe–Ni–Si outer core, 11.3±1.5 wt% silicon would be required to match the core's observed density at the core–mantle boundary. We have also performed first-principles calculations of the equations of state of Fe 3Si with the D0 3 structure, hcp iron, and FeSi with the B2 structure using density-functional theory.« less
  • Understanding the thermal-chemical state of the Earth's core requires knowledge of the thermal expansion of iron-rich alloys at megabar pressures and high temperatures. Our survey of literature revealed a significant lack of such data. We have determined the unit-cell parameters of the iron-sulfur compound Fe{sub 3}S by using synchrotron x-ray diffraction techniques and externally heated diamond-anvil cells at pressures up to 42.5 GPa and temperatures up to 900 K. The zero-pressure thermal expansivity of Fe{sub 3}S is determined in the form {alpha} = a{sub 1} + a{sub 2} T, where a{sub 1} = 3.0 {+-} 1.3 x 10{sup -5} K{supmore » -1} and {alpha}{sub 2} = 2.8 {+-} 1.5 x 10{sup -8} K{sup -2}. The temperature dependence of isothermal bulk modulus ({partial_derivative}K{sub T,0}/{partial_derivative}T){sub P} is estimated at -3.75 {+-} 1.80 x 10{sup -2} GPa K{sup -1}. Our data at 42.5 GPa and 900 K suggest that {approx}2.1 at. % (1.2 wt. %) sulfur produces 1% density deficit in iron. We have also carried out energy-dispersive x-ray diffraction measurements on pure iron and Fe{sub 0.864}Si{sub 0.136} alloy samples that were placed symmetrically in the same multianvil cell assemblies, using the SPring-8 synchrotron facility in Japan. Based on direct comparison of unit cell volumes under presumably identical pressures and temperatures, our data suggest that at most 3.2 at. % (1.6 wt. %) silicon is needed to produce 1% density deficit with respect to pure iron.« less
  • The theory of the earth's thermal history is investigated, taking account of the nonuniform energy density distribution from both the initial energy of accretion and the subsequent energy released in gravitational reorganization. The origin of a solid inner core is discussed.
  • The electronic structure of the clusters [Fe{sup m}O{sub n}]{sup m{minus}2n}, which model local structure defects arising in doping the tetragonal phase of yttrium-barium ceramics with iron atoms, has been calculated using the X{sub {alpha}}-discrete variation (DV) method. According to the results of calculations of the hyperfine magnetic field and the quadrupole splitting of {sup 57}Fe Moessbauer spectra in YBa{sub 2}(Cu{sub 1-x}Fe{sub x}){sub 3}O{sub 6.0}, the gamma resonance (GR) doublet with quadrupole splitting {vert_bar}{Delta}{vert_bar} {approximately} 2 mm/s is most likely attributable to the Fe atoms localized in the Cu(1) position and having the valence of 2+ and a dumbbell configuration ofmore » the nearest environment cluster.« less