An official website of the United States government
Abstract Nitrogen has been proposed to be stored within planetary cores, but its effects on the structure and density of molten Fe–alloys have not been explored experimentally. Using energy‐dispersive X‐ray diffraction, we determined the structure of Fe–N(–C) liquids at core conditions (1–7 GPa and 1700–1900°C) within a Paris‐Edinburgh press. Variation of N up to 7 wt.% and C up to 1.5 wt.% results in near‐linear changes in Fe–Fe atom distances and structure factor with increasing light element content. We did not observe a significant pressure‐driven structural transition in Fe–N(–C) liquids. We model the expansion of the Fe–Fe bonds using a modified Birch‐Murnaghan equation of state. With this model, we demonstrate that N or C contamination could lead to an overestimation of the Fe–Fe distances of pure Fe. We observe that the incorporation of 1 wt.% N or C into Fe results in a change in Fe–Fe distances that is twice as significant as the effect of 1 GPa. By approximating the change in volume, we infer that N and C incorporated in liquid iron could contribute to the density deficit observed in the cores of terrestrial bodies.
more »
« less
Effect of Carbon on the Volume of Solid Iron at High Pressure: Implications for Carbon Substitution in Iron Structures and Carbon Content in the Earth’s Inner Core
Understanding the effect of carbon on the density of hcp (hexagonal-close-packed) Fe-C alloys is essential for modeling the carbon content in the Earth’s inner core. Previous studies have focused on the equations of state of iron carbides that may not be applicable to the solid inner core that may incorporate carbon as dissolved carbon in metallic iron. Carbon substitution in hcp-Fe and its effect on the density have never been experimentally studied. We investigated the compression behavior of Fe-C alloys with 0.31 and 1.37 wt % carbon, along with pure iron as a reference, by in-situ X-ray diffraction measurements up to 135 GPa for pure Fe, and 87 GPa for Fe-0.31C and 109 GPa for Fe-1.37C. The results show that the incorporation of carbon in hcp-Fe leads to the expansion of the lattice, contrary to the known effect in body-centered cubic (bcc)-Fe, suggesting a change in the substitution mechanism or local environment. The data on axial compressibility suggest that increasing carbon content could enhance seismic anisotropy in the Earth’s inner core. The new thermoelastic parameters allow us to develop a thermoelastic model to estimate the carbon content in the inner core when carbon is incorporated as dissolved carbon hcp-Fe. The required carbon contents to explain the density deficit of Earth’s inner core are 1.30 and 0.43 wt % at inner core boundary temperatures of 5000 K and 7000 K, respectively.
more »
« less
- Award ID(s):
- 1619868
- PAR ID:
- 10224472
- Date Published:
- Journal Name:
- Minerals
- Volume:
- 9
- Issue:
- 12
- ISSN:
- 2075-163X
- Page Range / eLocation ID:
- 720
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Carbon and nitrogen are considered as candidate light elements present in planetary cores. However, there is limited understanding regarding the structure and physical properties of Fe‐C‐N alloys under extreme conditions. Here diamond anvil cell experiments were conducted, revealing the stability of hexagonal‐structured Fe7(N0.75C0.25)3up to 120 GPa and 2100 K, without undergoing any structural transformation or dissociation. Notably, the thermal expansion coefficient and Grüneisen parameter of the alloy exhibit a collapse at 55–70 GPa. First‐principles calculations suggest that such anomaly is associated with the spin transition of iron within Fe7(N0.75C0.25)3. Our modeling indicates that the presence of ∼1.0 wt% carbon and nitrogen in liquid iron contributes to 9–12% of the density deficit of the Earth's outer core. The thermoelastic anomaly of the Fe‐C‐N alloy across the spin transition is likely to affect the density and seismic velocity profiles of (C,N)‐rich planetary cores, thereby influencing the dynamics of such cores.more » « less
-
Abstract Hexagonal close‐packed (hcp) structured Fe‐Ni alloy is believed to be the dominant phase in the Earth's inner core. This phase is expected to contain 4%–5% light elements, such as Si and H. While the effects of individual light element candidates on the equation of state (EoS) of the hcp Fe metal have been studied, their combined effects remain largely unexplored. In this study, we report the equations of state for two hcp‐structured Fe‐Si‐H alloys, namely Fe0.83Si0.17H0.07and Fe0.83Si0.17H0.46, using synchrotron X‐ray diffraction measurements up to 125 GPa at 300 K. These alloys were synthesized by cold compression of Fe‐9wt%Si in either pure H2or Ar‐H2mixture medium in diamond‐anvil cells. The volume increase caused by a H atom in hcp Fe‐Si‐H alloys is approximately eight times greater than that by a Si atom. We used the improved data set to develop a composition‐dependent EoS that covers a wide range of compositions. Our calculated density and bulk sound velocity of hcp Fe‐Si‐H alloys suggest a large trade‐off between Si and H contents in fitting the seismic properties of the inner core. Combining our new EoS with geophysical and geochemical constraints, we propose 1.6–3 wt% Si and 0.15–0.6 wt% H in the Earth's inner core.more » « less
-
Abstract Accurate knowledge of the phase transitions and thermoelastic properties of candidate iron alloys, such as Fe‐Si alloys, is essential for understanding the nature and dynamics of planetary cores. The phase diagrams of some Fe‐Si alloys between 1 atm and 16 GPa have been back‐extrapolated from higher pressures, but the resulting phase diagram of Fe83.6Si16.4(9 wt.% Si) is inconsistent with temperature‐induced changes in its electrical resistivity between 6 and 8 GPa. This study reports in situ synchrotron X‐ray diffraction (XRD) measurements on pre‐melted and powder Fe83.6Si16.4samples from ambient conditions to 60 GPa and 900 K using an externally heated diamond‐anvil cell. Upon compression at 300 K, thebccphase persisted up to ∼38 GPa. Thehcpphase appeared near 8 GPa in the pre‐melted sample, and near 17 GPa in the powder sample. The appearance of thehcpphase in the pre‐melted sample reconciles the reported changes in electrical resistivity of a similar sample, thus resolving the low‐pressure region of the phase diagram. The resulting high‐temperature Birch‐Murnaghan equation of state (EoS) and thermal EoS based on the Mie‐Gruneisen‐Debye model of thebccandhcpstructures are consistent with, and complement the literature data at higher pressures. The calculated densities based on the thermal EoS of Fe‐9wt.%Si indicate that bothbccandhcpphases agree with the reported core density estimates for the Moon and Mercury.more » « less
-
Abstract Evaluating carbon’s candidacy as a light element in the Earth’s core is critical to constrain the budget and planet-scale distribution of this life-essential element. Here we use first principles molecular dynamics simulations to estimate the density and compressional wave velocity of liquid iron-carbon alloys with ~4-9 wt.% carbon at 0-360 gigapascals and 4000-7000 kelvin. We find that for an iron-carbon binary system, ~1-4 wt.% carbon can explain seismological compressional wave velocities. However, this is incompatible with the ~5-7 wt.% carbon that we find is required to explain the core’s density deficit. When we consider a ternary system including iron, carbon and another light element combined with additional constraints from iron meteorites and the density discontinuity at the inner-core boundary, we find that a carbon content of the outer core of 0.3-2.0 wt.%, is able to satisfy both properties. This could make the outer core the largest reservoir of terrestrial carbon.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3390/min9120720
