The structure of Fe–Si alloy in Earth's inner core
Introduction
Crystal structure in the Earth's inner core is of great importance to understand its elastic and transport properties, and dynamics (de Koker et al., 2012, Pozzo et al., 2012, Gomi et al., 2013). While a variety of different crystal structures have been proposed for iron at high pressure and temperature (P–T) (Saxena et al., 1995, Belonoshko et al., 2003, Mikhaylushkin et al., 2007), recent static experiments (Tateno et al., 2010) revealed that the hcp structure is a stable form of iron up to 377 GPa and 5700 K. Such experimental result is also supported by the most recent ab initio calculations (Stixrude, 2012). Nevertheless, the terrestrial core is not pure iron. Silicon is likely to be a major light component in the core (∼6 wt.%, Hirose et al., 2013), which is strongly favored by cosmochemical and geochemical arguments based on 1) the depletion in silicon in the Earth's mantle relative to solar abundance and primitive meteorites (e.g., Allègre et al., 1995) and 2) the difference in 30Si/28Si ratio between meteorites and terrestrial rocks (Georg et al., 2007, Fitoussi et al., 2009, Shahar et al., 2009, Ziegler et al., 2010). Such high concentration of silicon in the core has been supported by metal–silicate partitioning experiments (Takafuji et al., 2005, Sakai et al., 2006, Ozawa et al., 2009 Rubie et al., 2011), although some experimental and theoretical studies argued for a subordinate role of silicon (<4.5 wt.%) (Antonangeli et al., 2010, Siebert et al., 2013, Badro et al., 2014). Since only small fractionation of silicon is known to occur between liquid and solid iron, the inner core may possess similar amount of silicon to that in the bulk core (Kuwayama and Hirose, 2004; Alfè et al., 2002; Fischer et al., 2013). Indeed, previous density and sound velocity measurements of Fe and Fe–Si alloy at high pressure suggested that iron alloyed with silicon can explain seismic observations in the inner core (Badro et al., 2007, Mao et al., 2012).
Phase relations in Fe–Si alloys, however, have been examined by earlier experiments only to 257 GPa and 2400 K (Asanuma et al., 2008), while the inner core is under >330 GPa and ∼5000 K. The phase diagram of Fe–9 wt.% Si was extensively studied under pressure in a DAC (Kuwayama et al., 2009, Lin et al., 2009, Fischer et al., 2013). These previous experiments consistently demonstrated that hcp-Fe containing 9 wt.% silicon decomposes into a mixture of hcp and B2 phases with increasing temperature. It suggests a possibility of two-phase mixture in the inner core (Brosh et al., 2009, Fischer et al., 2013), although the boundary between hcp and hcp + B2 has never been constrained above 240 GPa. On the other hand, ab initio calculations argued that the incorporation of small amount of silicon in iron stabilizes body-centered cubic (bcc) or fcc phase with respect to the hcp phase (Vočadlo et al., 2003; Côté et al., 2008a, Côté et al., 2010), but it has never been verified by experiments.
In this study, we examine the phase diagram and the equation of state (EoS) of Fe–Si alloys at multimegabar pressures by a combination of laser-heated diamond–anvil cell (DAC) techniques and synchrotron X-ray diffraction (XRD) measurements. The phase relation in Fe–9 wt.% Si (Fe0.84Si0.16 in atomic ratio) was determined up to 407 GPa and 5960 K, the highest static P–T ever reported. The pressure–volume data for hcp Fe–9 wt.% Si was also obtained to 305 GPa at 300 K. In addition, the stabilities of hcp and fcc Fe–6.5 wt.% Si (Fe0.88Si0.12 in atomic ratio) were examined up to melting temperature at relatively low pressures (<70 GPa). Based on these new data, we discuss the stable crystal structure of Fe–Si alloy in the Earth's inner core.
Section snippets
High-P/room-T experiments for Fe–9 wt.% Si
High-pressure/room-temperature experiments were performed in a DAC to obtain the EoS of hcp Fe–9 wt.% Si (runs #1–4) (Table S1). We used a powder sample (Goodfellow) as starting material. We examined unheated portion of a recovered sample under scanning transmission electron microscope and confirmed the starting composition and the chemical homogeneity at micrometric scale. It was loaded into a hole in a pre-indented rhenium gasket together with MgO (pressure standard) and pressure medium (Ar
Equation of state of hcp Fe–9 wt.% Si
Four separate sets of experiments were performed to determine the EoS of hcp-structured Fe–9 wt.% Si at 300 K between 13 and 305 GPa (runs #1–4). The observed unit-cell parameters and volumes of hcp Fe–9 wt.% Si and MgO are summarized in Table S1. A representative XRD pattern collected at 305 GPa at room temperature is shown in Fig. 1a, in which the diffraction peaks from hcp Fe–9 wt.% Si and NaCl-type (B1)-MgO were found. The diffraction peaks of SiO2 were not observed through the present
Crystal structure of Fe–Si alloy in the inner core
The density of the inner core deduced from seismological observations is smaller by 5 to 6% than that of pure iron (Dziewonski and Anderson, 1981, Dewaele et al., 2006) when ICB temperature is 4900–5700 K (Alfè et al., 2007; Anzellini et al., 2013, Nomura et al., 2014) (Fig. 10). The EoS of Fe–9 wt.% Si obtained in this study indicates that such inner core density deficit is reconciled with an hcp iron-alloy containing 6 to 8 wt.% Si at the ICB. Additionally, the high Mg/Si ratio (∼1.3) in the
Conclusions
The phase relations and the EoS of Fe–9 wt.% Si were examined to a multimegabar pressure range. Experiments were performed up to 407 GPa and 5960 K in a laser-heated DAC, the highest static P–T condition ever reported. The results demonstrate that iron with an hcp structure contains up to 9 wt.% Si at 4800 K and 330 GPa of ICB pressure. The Earth's inner core may include about 7 wt.% Si to account for its density deficit from pure iron according to the present density measurements, which is
Acknowledgments
We thank N. Hirao, H. Ozawa, and Y. Kudo for their technical assistance. T. Komabayashi is thanked for fruitful discussions. The comments from two anonymous reviewers were helpful to improve the manuscript. The synchrotron X-ray diffraction measurements were conducted at BL10XU of SPring-8 (proposal Nos. 2011B0087, 2011B1405, 2012A0087, 2012B1212, and 2013B0087).
References (56)
- et al.
Composition and temperature of the Earth's core constrained by combining ab initio calculations and seismic data
Earth Planet. Sci. Lett.
(2002) - et al.
The chemical composition of the Earth
Earth Planet. Sci. Lett.
(1995) - et al.
Composition of the Earth's inner core from high-pressure sound velocity measurements in Fe–Ni–Si alloys
Earth Planet. Sci. Lett.
(2010) - et al.
Effect of light elements on the sound velocities in solid iron: implications for the composition of Earth's core
Earth Planet. Sci. Lett.
(2007) - et al.
Thermodynamic analysis of high-pressure phase equilibria in Fe–Si alloys, implications for the inner-core
Phys. Earth Planet. Inter.
(2009) - et al.
The effect of silicon impurities on the phase diagram of iron and possible implications for the Earth's core structure
J. Phys. Chem. Solids
(2008) - et al.
Ab initio lattice dynamics calculations on the combined effect of temperature and silicon on the stability of different iron phases in the Earth's inner core
Phys. Earth Planet. Inter.
(2010) Structure and dynamics of Earth's inner core
Earth Planet. Sci. Lett.
(2012)- et al.
Preliminary reference Earth model
Phys. Earth Planet. Inter.
(1981) - et al.
Phase relations in the Fe–FeSi system at high pressures and temperatures
Earth Planet. Sci. Lett.
(2013)
Si isotope systematics of meteorites and terrestrial peridotites: implications for Mg/Si fractionation in the solar nebula and for Si in the Earth's core
Earth Planet. Sci. Lett.
The axial ratio of hcp iron at the conditions of the Earth's inner core
Phys. Earth Planet. Inter.
The high conductivity of iron and thermal evolution of the Earth's core
Phys. Earth Planet. Inter.
Melting relationships in the Fe–Fe3S system up to the outer core conditions
Earth Planet. Sci. Lett.
In-situ X-ray diffraction measurements of the γ–ε transition boundary of iron in an internally-heated diamond anvil cell
Earth Planet. Sci. Lett.
Experimentally determined Si isotope fractionation between silicate and Fe metal and implications for Earth's core formation
Earth Planet. Sci. Lett.
Metal–silicate silicon isotope fractionation in enstatite meteorites and constraints on Earth's core formation
Earth Planet. Sci. Lett.
Temperature and composition of the Earth's core
Contemp. Phys.
Melting of iron at Earth's inner core boundary based on fast X-ray diffraction
Science
Phase relations of Fe–Si alloy up to core conditions: implications for the Earth inner core
Geophys. Res. Lett.
A seismically consistent compositional model of the core
Proc. Natl. Acad. Sci. USA
Stability of the body-centered-cubic phase of iron in the Earth's inner core
Nature
Light elements in the core: effects of impurities on the phase diagram of iron
Geophys. Res. Lett.
Electrical resistivity and thermal conductivity of liquid Fe alloys at high P and T, and heat flux in Earth's core
Proc. Natl. Acad. Sci. USA
Quasihydrostatic equation of state of iron above 2 Mbar
Phys. Rev. Lett.
Equation of state in the Fe–FeSi system at high pressure and temperatures
J. Geophys. Res.
Silicon in the Earth's core
Nature
Equation of state of iron–silicon alloys to megabar pressure
Phys. Chem. Miner.
Cited by (80)
High-pressure fluid-phase equilibria: Experimental methods, developments and systems investigated (2013–2016)
2024, Fluid Phase EquilibriaEquations of state for B2 and bcc Fe<inf>1-x</inf>Si<inf>x</inf>
2023, Physics of the Earth and Planetary InteriorsShear softening of earth's inner core as indicated by its high poisson ratio and elastic anisotropy
2022, Fundamental ResearchHigh-pressure experimental constraints of partitioning behavior of Si and S at the Mercury's inner core boundary
2021, Earth and Planetary Science Letters