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
This content will become publicly available on October 29, 2025
Phase Transitions and Thermal Equation of State of Fe‐9wt.%Si Applied to the Moon and Mercury
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
- PAR ID:
- 10582455
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Planets
- Volume:
- 129
- Issue:
- 11
- ISSN:
- 2169-9097
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Mechanisms of hexagonal close-packed (HCP) to body-centered cubic (BCC) phase transformation in Mg single crystals are observed using a combination of polychromatic beam Laue diffraction and monochromatic beam powder diffraction techniques under quasi-hydrostatic pressures of up to 58 ± 2 GPa at ambient temperature. Although experiments were performed with both He and Ne pressure media, crystals inevitably undergo plastic deformation upon loading to 40–44 GPa. The plasticity is accommodated by dislocation glide causing local misorientations of up to 1°–2°. The selected crystals are tracked by mapping Laue diffraction spots up to the onset of the HCP to BCC transformation, which is determined to be at a pressure of 56.6 ± 2 GPa. Intensity of the Laue reflections from HCP crystals rapidly decrease but no reflections from crystalline BCC phase are observed with a further increase of pressure. Nevertheless, the powder diffraction shows the formation of 110 BCC peak at 56.6 GPa. The peak intensity increases at 59.7 GPa. Upon the full transformation, a powder-like BCC aggregate is formed revealing the destructive nature of the HCP to BCC transformation in single crystals of Mg.more » « less
-
null (Ed.)Natural kamacite samples (Fe92.5Ni7.5) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the c/a ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal c/a ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth’s core elasticity.more » « less
-
Abstract The motion of liquid iron (Fe) alloy materials in the outer core drives the dynamo, which generates Mercury's magnetic field. The assessment of core models requires laboratory measurements of the melting temperature of Fe alloys at high pressure. Here, we experimentally determined the melting curve of Fe9wt%Si and Fe17wt%Si up to 17 GPa using in situ and ex situ measurements of intermetallic fast diffusion that serves as the melting criterion in a large‐volume press. Our determined melting slopes are comparable with previous studies up to about 17 GPa. However, when extrapolated, our melting slopes significantly deviate from previous studies at higher pressures. For Mercury's core with a model composition of Fe9wt%Si, the melting temperature‐depth profile determined in our study is lower by ∼150–250 K when compared with theoretical calculations. Using the new melting curve of Fe9wt%Si and the electrical resistivity values from a previous study of Fe8.5wt%Si, we estimate that the electronic thermal conductivity of liquid Fe9wt%Si is 30 Wm−1K−1at the Mercury'sCMBpressure of 5 GPa and 37 Wm−1K−1at an assumedICBof 21 GPa, corresponding to heat flux values of 23 mWm−2and 32 mWm−2, respectively. These values provide new constraints on the core models.more » « less