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: Viscosity of Fe‐Ni‐C Liquids up to Core Pressures and Implications for Dynamics of Planetary Cores
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
Award ID(s):
2017294
PAR ID:
10369679
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
4
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. Abstract The occurrences and cycling of slab‐originated carbon and hydrogen are considered to be controlled by their reactions with metallic iron from mantle disproportionation and slab serpentinization, to form Fe alloys containing carbon and hydrogen. Here we show experimental results on the phase relations and melting of the Fe‐C‐H system using laser‐heated diamond anvil cell and X‐ray diffraction techniques up to 72 GPa. The incorporation of hydrogen was found to lower the eutectic melting temperatures of Fe‐C alloy by ∼50–178 K at 20–70 GPa, facilitating the formation of metallic liquids in the deep mantle and thus enhancing the mobility and deep cycling of subducted carbon and hydrogen. Hydrogen also substitutes with carbon in Fe‐C metal to form hydride and diamond at relatively high‐temperature conditions (e.g., 42.6 GPa, >1885 K and 71.8 GPa, >1798 K). The hydrogen‐carbon‐enriched metallic liquids provide the necessary fluid environment for superdeep diamond growth. 
    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