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.


This content will become publicly available on December 1, 2025

Title: High Pressure Melting Curve of Fe‐Si: Implication for the Thermal Properties in Mercury's Core
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
Award ID(s):
2022492
PAR ID:
10597565
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Journal of Geophysical Research: Planets
Volume:
129
Issue:
12
ISSN:
2169-9097
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The high-pressure melting curve of FeO controls key aspects of Earth’s deep interior and the evolution of rocky planets more broadly. However, existing melting studies on wüstite were conducted across a limited pressure range and exhibit substantial disagreement. Here we use an in-situ dual-technique approach that combines a suite of >1000 x-ray diffraction and synchrotron Mössbauer measurements to report the melting curve for Fe1-xO wüstite to pressures of Earth’s lowermost mantle. We further observe features in the data suggesting an order-disorder transition in the iron defect structure several hundred kelvin below melting. This solid-solid transition, suggested by decades of ambient pressure research, is detected across the full pressure range of the study (30 to 140 GPa). At 136 GPa, our results constrain a relatively high melting temperature of 4140 ± 110 K, which falls above recent temperature estimates for Earth’s present-day core-mantle boundary and supports the viability of solid FeO-rich structures at the roots of mantle plumes. The coincidence of the defect order-disorder transition with pressure-temperature conditions of Earth’s mantle base raises broad questions about its possible influence on key physical properties of the region, including rheology and conductivity. 
    more » « less
  2. 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
  3. 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 yieldsV0 = 172.1(4) Å3,K0 = 229(4) GPa withK0′ = 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. 
    more » « less
  4. Abstract The heat extracted from the core by the overlying mantle across the core‐mantle boundary controls the thermal evolution of the core. This in turn leads to the solidification of the inner core in association with the exsolution of light alloying elements into the liquid outer core. Although the temperature (T) at the inner core boundary (ICB) would be adjusted to account for the effects of the light elements, the melting T of Fe places an upper bound at the ICB and it is a vital point in the thermal profile of the core. Here, we determine the melting T of Fe in the multi‐anvil press by characterizing the interface of Fe‐W interaction. Our data place a tighter constraint on the melting curve of Fe between 8 and 21 GPa, that is directly applicable to small planetary bodies and serves as an anchor for melting curve of Fe at higher pressure. 
    more » « less
  5. We report experimental constraints on the melting curve of potassium chloride (KCl) between 3.2 and 9 GPa from in situ ionic conduction measurements using a multi-anvil apparatus. On the basis of concurrent measurements of KCl and sodium chloride (NaCl) at 1 bar using the differential thermal analysis (DTA) method and Pt sphere marker, we show that the peak rate of increase in ionic current with temperature upon heating coincides with latent heat ledge and fall of Pt sphere, thus establishing the criterion for melting detection from ionic conduction measurements. Applying this criterion to high pressures, we found that the melting point of KCl rose steeply with increasing pressure to exceed 2443 ± 100 K at 9 GPa. Fitting the results of this study together with existing data at pressures below 4 GPa and above 20 GPa, we obtained the Simon’s melting equation for KCl in the simple cubic B2 structure between 1.8 and 50 GPa: T m = 1323 ( P − 1.87 2.2 ( 1 ) + 1 ) 1 2.7 ( 1 ) , where T is in K and P is in GPa. Starting at 1 bar, the melting point of KCl increases at an average rate of ~150 K/GPa to cross that of Pt near 9 GPa. The highly refractory nature of KCl makes it a sensitive pressure calibrant for the large-volume pressure at moderate pressures and a potential sample container for experiments at moderate pressures and very high temperatures. 
    more » « less