Hidden carbon in Earth's inner core revealed by shear softening in dense Fe7C3

Iron alloys of volatile elements in the deep Earth's interior

Investigations into the compositional model of the Earth, particularly the atypical concentrations of volatile elements within the silicate portion of the early Earth, have attracted significant interest due to their pivotal role in elucidating the planet's evolution and dynamics. To understand the behavior of such volatile elements, an established 'volatility trend' has been used to explain the observed depletion of certain volatile elements. However, elements such as Se and Br remain notably over-depleted in the silicate Earth. Here we show the results from first-principles simulations that explore the potential for these elements to integrate into hcp-Fe through the formation of substitutional alloys, long presumed to be predominant constituents of the Earth's core. Based on our findings, the thermodynamic stability of these alloys suggests that these volatile elements might indeed be partially sequestered within the Earth's core. We suggest potential reservoirs for volatile elements within the deep Earth, augmenting our understanding of the deep Earth's composition.

Superionic iron alloys and their seismic velocities in Earth's inner core

Earth's inner core (IC) is less dense than pure iron, indicating the existence of light elements within it¹. Silicon, sulfur, carbon, oxygen and hydrogen have been suggested to be the candidates^2,3, and the properties of iron-light-element alloys have been studied to constrain the IC composition^4-19. Light elements have a substantial influence on the seismic velocities^4-13, the melting temperatures^14-17 and the thermal conductivities^18,19 of iron alloys. However, the state of the light elements in the IC is rarely considered. Here, using ab initio molecular dynamics simulations, we find that hydrogen, oxygen and carbon in hexagonal close-packed iron transform to a superionic state under the IC conditions, showing high diffusion coefficients like a liquid. This suggests that the IC can be in a superionic state rather than a normal solid state. The liquid-like light elements lead to a substantial reduction in the seismic velocities, which approach the seismological observations of the IC^20,21. The substantial decrease in shear-wave velocity provides an explanation for the soft IC²¹. In addition, the light-element convection has a potential influence on the IC seismological structure and magnetic field.

Inner core composition paradox revealed by sound velocities of Fe and Fe-Si alloy

Knowledge of the sound velocity of core materials is essential to explain the observed anomalously low shear wave velocity (VS) and high Poisson's ratio (σ) in the solid inner core. To date, neither VS nor σ of Fe and Fe-Si alloy have been measured under core conditions. Here, we present VS and σ derived from direct measurements of the compressional wave velocity, bulk sound velocity, and density of Fe and Fe-8.6 wt%Si up to ~230 GPa and ~5400 K. The new data show that neither the effect of temperature nor incorporation of Si would be sufficient to explain the observed low VS and high σ of the inner core. A possible solution would add carbon (C) into the solid inner core that could further decrease VS and increase σ. However, the physical property-based Fe-Si-C core models seemingly conflict with the partitioning behavior of Si and C between liquid and solid Fe.

A pebble accretion model for the formation of the terrestrial planets in the Solar System

Pebbles of millimeter sizes are abundant in protoplanetary discs around young stars. Chondrules inside primitive meteorites-formed by melting of dust aggregate pebbles or in impacts between planetesimals-have similar sizes. The role of pebble accretion for terrestrial planet formation is nevertheless unclear. Here, we present a model where inward-drifting pebbles feed the growth of terrestrial planets. The masses and orbits of Venus, Earth, Theia (which later collided with Earth to form the Moon), and Mars are all consistent with pebble accretion onto protoplanets that formed around Mars' orbit and migrated to their final positions while growing. The isotopic compositions of Earth and Mars are matched qualitatively by accretion of two generations of pebbles, carrying distinct isotopic signatures. Last, we show that the water and carbon budget of Earth can be delivered by pebbles from the early generation before the gas envelope became hot enough to vaporize volatiles.

The carbon content of Earth and its core

Earth's core is likely the largest reservoir of carbon (C) in the planet, but its C abundance has been poorly constrained because measurements of carbon's preference for core versus mantle materials at the pressures and temperatures of core formation are lacking. Using metal-silicate partitioning experiments in a laser-heated diamond anvil cell, we show that carbon becomes significantly less siderophile as pressures and temperatures increase to those expected in a deep magma ocean during formation of Earth's core. Based on a multistage model of core formation, the core likely contains a maximum of 0.09(4) to 0.20(10) wt% C, making carbon a negligible contributor to the core's composition and density. However, this accounts for ∼80 to 90% of Earth's overall carbon inventory, which totals 370(150) to 740(370) ppm. The bulk Earth's carbon/sulfur ratio is best explained by the delivery of most of Earth's volatiles from carbonaceous chondrite-like precursors.

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.

New high-pressure phases of Fe7N3 and Fe7C3 stable at Earth's core conditions: evidences for carbon–nitrogen isomorphism in Fe-compounds

We carried out ab initio calculations on the crystal structure prediction and determination of P–T diagrams within the quasi-harmonic approximation for Fe₇N₃ and Fe₇C₃. Two new isostructural phases Fe₇N₃-C2/m and Fe₇C₃-C2/m which are dynamically and thermodynamically stable under the Earth's core conditions were predicted. The Fe₇C₃-C2/m phase stabilizes preferentially to the known h-Fe₇C₃ at 253–344 GPa in the temperature range of 0–5000 K, and the Fe₇N₃-C2/m stabilizes preferentially relative to the β-Fe₇N₃ – at ∼305 GPa over the entire temperature range. This indicate that carbon and nitrogen can mutually coexist and replace each other in the Earth's and other planetary cores similarly to low pressure phases of the same compounds.

We carried out ab initio calculations on the crystal structure prediction and determination of P–T diagrams within the quasi-harmonic approximation for Fe₇N₃ and Fe₇C₃.

Abnormal Elasticity of Single-Crystal Magnesiosiderite across the Spin Transition in Earth's Lower Mantle

Brillouin light scattering and impulsive stimulated light scattering have been used to determine the full elastic constants of magnesiosiderite [(Mg_{0.35}Fe_{0.65})CO_{3}] up to 70 GPa at room temperature in a diamond-anvil cell. Drastic softening in C_{11}, C_{33}, C_{12}, and C_{13} elastic moduli associated with the compressive stress component and stiffening in C_{44} and C_{14} moduli associated with the shear stress component are observed to occur within the spin transition between ∼42.4 and ∼46.5  GPa. Negative values of C_{12} and C_{13} are also observed within the spin transition region. The Born criteria constants for the crystal remain positive within the spin transition, indicating that the mixed-spin state remains mechanically stable. Significant auxeticity can be related to the electronic spin transition-induced elastic anomalies based on the analysis of Poisson's ratio. These elastic anomalies are explained using a thermoelastic model for the rhombohedral system. Finally, we conclude that mixed-spin state ferromagnesite, which is potentially a major deep-carbon carrier, is expected to exhibit abnormal elasticity, including a negative Poisson's ratio of -0.6 and drastically reduced V_{P} by 10%, in Earth's midlower mantle.

Pressure-dependent isotopic composition of iron alloys

Our current understanding of Earth's core formation is limited by the fact that this profound event is far removed from us physically and temporally. The composition of the iron metal in the core was a result of the conditions of its formation, which has important implications for our planet's geochemical evolution and physical history. We present experimental and theoretical evidence for the effect of pressure on iron isotopic composition, which we found to vary according to the alloy tested (FeO, FeH(x), or Fe3C versus pure Fe). These results suggest that hydrogen or carbon is not the major light-element component in the core. The pressure dependence of iron isotopic composition provides an independent constraint on Earth's core composition.

Carbon-depleted outer core revealed by sound velocity measurements of liquid iron-carbon alloy

The relative abundance of light elements in the Earth's core has long been controversial. Recently, the presence of carbon in the core has been emphasized, because the density and sound velocities of the inner core may be consistent with solid Fe₇C₃. Here we report the longitudinal wave velocity of liquid Fe₈₄C₁₆ up to 70 GPa based on inelastic X-ray scattering measurements. We find the velocity to be substantially slower than that of solid iron and Fe₃C and to be faster than that of liquid iron. The thermodynamic equation of state for liquid Fe₈₄C₁₆ is also obtained from the velocity data combined with previous density measurements at 1 bar. The longitudinal velocity of the outer core, about 4% faster than that of liquid iron, is consistent with the presence of 4–5 at.% carbon. However, that amount of carbon is too small to account for the outer core density deficit, suggesting that carbon cannot be a predominant light element in the core.

The composition of the Earth's core, particularly the light elements present, is not well constrained. Here, the authors report sound velocities of liquid iron-carbon alloy as measured at very high pressures using inelastic X-ray scattering and suggest that carbon cannot be predominant in the outer core.