Earth’s inner core is predominantly composed of solid iron (Fe) and displays intriguing properties such as strong shear softening and an ultrahigh Poisson’s ratio. Insofar, physical mechanisms to explain these features coherently remain highly debated. Here, we have studied longitudinal and shear wave velocities of hcp-Fe (hexagonal close-packed iron) at relevant pressure–temperature conditions of the inner core using in situ shock experiments and machine learning molecular dynamics (MLMD) simulations. Our results demonstrate that the shear wave velocity of hcp-Fe along the Hugoniot in the premelting condition, defined as
We conducted shock wave experiments on iron carbide Fe3C up to a Hugoniot pressure of 245 GPa. The correlation between the particle velocity (
- Award ID(s):
- 1619868
- NSF-PAR ID:
- 10446330
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 46
- Issue:
- 20
- ISSN:
- 0094-8276
- Page Range / eLocation ID:
- p. 11018-11024
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
T/Tm (Tm : melting temperature of iron) above 0.96, is significantly reduced by ~30%, while Poisson’s ratio jumps to approximately 0.44. MLMD simulations at 230 to 330 GPa indicate that collective motion with fast diffusive atomic migration occurs in premelting hcp-Fe primarily along [100] or [010] crystallographic direction, contributing to its elastic softening and enhanced Poisson’s ratio. Our study reveals that hcp-Fe atoms can diffusively migrate to neighboring positions, forming open-loop and close-loop clusters in the inner core conditions. Hcp-Fe with collective motion at the inner core conditions is thus not an ideal solid previously believed. The premelting hcp-Fe with collective motion behaves like an extremely soft solid with an ultralow shear modulus and an ultrahigh Poisson’s ratio that are consistent with seismic observations of the region. Our findings indicate that premelting hcp-Fe with fast diffusive motion represents the underlying physical mechanism to help explain the unique seismic and geodynamic features of the inner core. -
Abstract We present ab initio (LDA + U
s c ) studies of high‐temperature and high‐pressure elastic properties of pure as well as iron‐bearing (ferrous, Fe2+, and ferric, Fe3+) and aluminum‐bearing MgSiO3postperovskite, the likely dominant phase in the deep lower mantle of the Earth. Thermal effects are addressed within the quasiharmonic approximation by combining vibrational density of states and static elastic coefficients. Aggregate elastic moduli and sound velocities for the Mg end members are successfully compared with scarce experimental data available. Effects of iron (Fe) and aluminum (Al) substitutions on elastic properties and their pressure and temperature dependence have been thoroughly investigated. At the observed perovskite to postperovskite transition (P = 125 GPa andT = 2,500 K), compressional and shear velocities increase by 0–1% and 1.5–3.75%, respectively. This observation is consistent with some seismic studies of the D′ ′ discontinuity beneath the Caribbean, which suggests that our robust estimates of elastic properties of the postperovskite phase will be very helpful to understand lateral velocity variations in the deep lower mantle region and to constrain its composition and thermal structure. -
Abstract Knowledge of the sound velocity of core materials is essential to explain the observed anomalously low shear wave velocity (
V S) and high Poisson’s ratio (σ ) in the solid inner core. To date, neitherV Snorσ of Fe and Fe-Si alloy have been measured under core conditions. Here, we presentV Sandσ 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 lowV Sand highσ of the inner core. A possible solution would add carbon (C) into the solid inner core that could further decreaseV Sand 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. -
Abstract Iron nitrides are possible constituents of the cores of Earth and other terrestrial planets. Pressure‐induced magnetic changes in iron nitrides and effects on compressibility remain poorly understood. Here we report synchrotron X‐ray emission spectroscopy (XES) and X‐ray diffraction (XRD) results for ε‐Fe7N3and γ′‐Fe4N up to 60 GPa at 300 K. The XES spectra reveal completion of high‐ to low‐spin transition in ε‐Fe7N3and γ′‐Fe4N at 43 and 34 GPa, respectively. The completion of the spin transition induces stiffening in bulk modulus of ε‐Fe7N3by 22% at ~40 GPa, but has no resolvable effect on the compression behavior of γ′‐Fe4N. Fitting pressure‐volume data to the Birch‐Murnaghan equation of state yields
V 0 = 83.29 ± 0.03 (Å3),K 0 = 232 ± 9 GPa,K 0′ = 4.1 ± 0.5 for nonmagnetic ε‐Fe7N3above the spin transition completion pressure, andV 0 = 54.82 ± 0.02 (Å3),K 0 = 152 ± 2 GPa,K 0′ = 4.0 ± 0.1 for γ′‐Fe4N over the studied pressure range. By reexamining evidence for spin transition and effects on compressibility of other candidate components of terrestrial planet cores, Fe3S, Fe3P, Fe7C3, and Fe3C based on previous XES and XRD measurements, we located the completion of high‐ to low‐spin transition at ~67, 38, 50, and 30 GPa at 300 K, respectively. The completion of spin transitions of Fe3S, Fe3P, and Fe3C induces elastic stiffening, whereas that of Fe7C3induces elastic softening. Changes in compressibility at completion of spin transitions in iron‐light element alloys may influence the properties of Earth's and planetary cores. -
Abstract Fe‐Al‐bearing bridgmanite may be the dominant host for ferric iron in Earth's lower mantle. Here we report the synthesis of (Mg0.5Fe3+0.5)(Al0.5Si0.5)O3bridgmanite (FA50) with the highest Fe3+‐Al3+coupled substitution known to date. X‐ray diffraction measurements showed that at ambient conditions, the FA50 adopted the LiNbO3structure. Upon compression at room temperature to 18 GPa, it transformed back into the bridgmanite structure, which remained stable up to 102 GPa and 2,600 K. Fitting Birch‐Murnaghan equation of state of FA50 bridgmanite yields
V 0 = 172.1(4) Å3,K 0 = 229(4) GPa withK 0′ = 4(fixed). The calculated bulk sound velocity of the FA50 bridgmanite is ~7.7% lower than MgSiO3bridgmanite, mainly because the presence of ferric iron increases the unit‐cell mass by 15.5%. This difference likely represents the upper limit of sound velocity anomaly introduced by Fe3+‐Al3+substitution. X‐ray emission and synchrotron Mössbauer spectroscopy measurements showed that after laser annealing, ~6% of Fe3+cations exchanged with Al3+and underwent the high‐ to low‐spin transition at 59 GPa. The low‐spin proportion of Fe3+increased gradually with pressure and reached 17–31% at 80 GPa. Since the cation exchange and spin transition in this Fe3+‐Al3+‐enriched bridgmanite do not cause resolvable unit‐cell volume reduction, and the increase of low‐spin Fe3+fraction with pressure occurs gradually, the spin transition would not produce a distinct seismic signature in the lower mantle. However, it may influence iron partitioning and isotopic fractionation, thus introducing chemical heterogeneity in the lower mantle.