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: Density of Fe‐Ni‐C Liquids at High Pressures and Implications for Liquid Cores of Earth and the Moon
Abstract The presence of light elements in the metallic cores of the Earth, the Moon, and other rocky planetary bodies has been widely proposed. Carbon is among the top candidates in light of its high cosmic abundance, siderophile nature, and ubiquity in iron meteorites. It is, however, still controversial whether carbon‐rich core compositional models can account for the seismic velocity observations within the Earth and lunar cores. Here, we report the density and elasticity of Fe90Ni10‐3 wt.% C and Fe90Ni10‐5 wt.% C liquid alloys using synchrotron‐based X‐ray absorption experiments and first‐principles molecular dynamics simulations. Our results show that alloying of 3 wt.% and 5 wt.% C lowers the density of Fe90Ni10liquid by ∼2.9–3.1% at 2 GPa, and ∼3.4–3.6% at 9 GPa. More intriguingly, our experiments and simulations both demonstrate that the bulk moduli of the Fe‐Ni‐C liquids are similar to or slightly higher than those of Fe‐Ni liquids. Thus, the calculated compressional velocities (vp) of Fe‐Ni‐C liquids are higher than that of pure Fe‐Ni alloy, promoting carbon as a possible candidate to explain the elevatedvpin the Earth's outer core. However, the values and slopes of both density andvpof the studied two Fe‐Ni‐C liquids do not match the outer core seismic models, suggesting that carbon may not be the sole principal light element in Earth's outer core. The highvpof Fe‐Ni‐C liquids does not match the presumptivevpof the lunar outer core well, indicating that carbon is less likely to be its dominant light element.  more » « less
Award ID(s):
2017294 1565708 1555388 1829273
PAR ID:
10385980
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
3
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract The viscosity of iron alloy liquids is the key for the core dynamo and core‐mantle differentiation of terrestrial bodies. Here we measured the viscosity of Fe‐Ni‐C liquids up to 7 GPa using the floating sphere viscometry method and up to 330 GPa using first‐principles calculations. We found a viscosity increase at ∼3–5 GPa, coincident with a structural transition in the liquids. After the transition, the viscosity reaches ∼14–27 mPa·s, a factor of 2–4 higher than that of Fe and Fe‐S liquids. Our computational results from 5 to 330 GPa also indicate a high viscosity of the Fe‐Ni‐C liquids. For a carbon‐rich core in large terrestrial body, the level of turbulence in the outer core would be lessened approaching the inner core boundary. It is also anticipated that Fe‐Ni‐C liquids would percolate in Earth's deep silicate mantle at a much slower speed than Fe and Fe‐S liquids. 
    more » « less
  3. 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
  4. 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
  5. 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