Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5338021

Formats

Article sections

- Abstract
- Thermodynamic model for solid and liquid iron
- EoSs for solid and liquid Fe to 350GPa
- Discussion and geophysical implications
- Summary
- Additional Information
- Supplementary Material
- References

Authors

Related links

Sci Rep. 2017; 7: 41863.

Published online 2017 March 6. doi: 10.1038/srep41863

PMCID: PMC5338021

Received 2016 October 10; Accepted 2016 December 29.

Copyright © 2017, The Author(s)

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

The equations of state for solid (with bcc, fcc, and hcp structures) and liquid phases of Fe were defined via simultaneous optimization of the heat capacity, bulk moduli, thermal expansion, and volume at room and higher temperatures. The calculated triple points at the phase diagram have the following parameters: bcc–fcc–hcp is located at 7.3GPa and 820K, bcc–fcc–liquid at 5.2GPa and 1998K, and fcc–hcp–liquid at 106.5GPa and 3787K. At conditions near the fcc–hcp–liquid triple point, the Clapeyron slope of the fcc–liquid curve is d*T*/d*P*=12.8K/GPa while the slope of the hcp–liquid curve is higher (d*T*/d*P*=13.7K/GPa). Therefore, the hcp–liquid curve overlaps the metastable fcc–liquid curve at pressures of about 160GPa. At high-pressure conditions, the metastable bcc–hcp curve is located inside the fcc-Fe or liquid stability field. The density, adiabatic bulk modulus and P-wave velocity of liquid Fe calculated up to 328.9GPa at adiabatic temperature conditions started from 5882K (outer/inner core boundary) were compared to the PREM seismological model. We determined the density deficit of hcp-Fe at the inner core boundary (*T*=5882K and *P*=328.9GPa) to be 4.4%.

Iron is a major component of the Earth’s core, therefore knowledge of its *P–V–T* relations and thermodynamic properties is extremely important^{1}^{,2}^{,3}^{,4}^{,5}^{,6}. The phase diagram of Fe is relatively complex due to the existence of several polymorphic modifications^{7}^{,8}. At the standard conditions (*T*=298.15K and *P*=1bar), iron is a ferromagnet and has a body-centred cubic (bcc) structure (α-Fe or bcc-Fe). The Curie temperature (*T*_{C}) of 1043K marks the transition to the paramagnetic state with the same structure. At this transition, the heat capacity of Fe has a characteristic λ-shape form with a maximum at *T*_{C}^{9}^{,10}. At 1185–1667K the crystal structure of iron changes to a face-centred cubic (fcc) cell (γ-Fe or fcc-Fe), however, above 1667K and up to the melting temperature of 1811K, iron again has the bcc structure (δ-Fe).

At 10.5GPa and 753K^{7} (or by more precise and recent measurements at 8.2GPa and 678K^{11}), there is a triple point between bcc-Fe, fcc-Fe, and the high-pressure phase hcp-Fe, which has the hexagonal close-packed (hcp) structure (hcp-Fe). It was argued that hcp-Fe is likely a stable phase in the inner core of the Earth^{12}^{,13}^{,14}. However, bcc-Fe is also suggested as a reliable candidate in the inner core^{15}^{,16}. There are two triple points along the melting line of Fe. The bcc-Fe and fcc-Fe phases are stable with liquid at the first triple point at *Р*=5.2GPa and *T*=1991K^{7}, whereas hcp-Fe and fcc-Fe were found in equilibrium with liquid at the second triple point, whose location is variable in different works (e.g. *Р*=88GPa and *Т*=2800K^{17} or *Р*=98.5GPa and *Т*=3712K^{14}).

The most recent thermodynamic calculations of the Fe phase diagram including solid phases and liquid were performed in a few studies^{18}^{,19}^{,20}^{,21} using the CALPHAD thermodynamic formalism^{10}^{,22}. In these calculations, the volume contribution to the Gibbs energy was calculated using the Anderson–Grüneisen parameter^{23}. In addition, a wide-range multiphase equation of state (EoS) to 10 TPa and 10^{5}K was calculated^{24}.

During the last several years, a significant amount of new *P–V–T* data for fcc-Fe and hcp-Fe, especially at very high temperatures, have appeared^{6}^{,13}^{,25}^{,26}^{,27}^{,28} and, in addition, the melting curve of Fe was shifted to higher temperatures according to measurements in ref. ^{14}. In these works, various pressure scales were used based on the EoSs of Au, MgO, NaCl, KCl, and hcp-Fe. Thus, the problem of consistent EoSs for solid phases (bcc-Fe, fcc-Fe, and hcp-Fe) and liquid iron remains extremely important. These EoSs should be consistent with the phase diagram of iron, including the melting curve and with *P–V–T*, thermochemical, and sound velocity data at 0.1MPa.

In this work, we propose updated EoSs for Fe phases based on a formalism from our recent publications^{29}^{,30}^{,31}^{,32}^{,33}^{,34}^{,35}^{,36}. In addition, the magnetic contribution to the Helmholtz free energy was calculated according to ref. ^{10} and ^{37}, whereas the EoS of Fe liquid was calculated using the standard approach^{38}^{,39}^{,40}.

The Helmholtz free energy of solid phases can be expressed in its classical form^{41} as:

where *U*_{0} is the reference energy, *E*_{0}(*V*) is the potential (cold) part of the free energy at the reference isotherm *T*_{0}=298.15K, which depends only on *V, F*_{th}(*V, T*) is the thermal part of the free energy, which depends on *V* and *T, F*_{e}(*V, T*) is the free electrons’ contribution to the Helmholtz free energy, which also depends on *V* and *Т*, and *F*_{mag}(*T*) is the magnetic contribution, depending on *T* only.

The pressure at 298K isotherm is calculated from the Vinet–Rydberg equation^{42}:

where *X*=(*V*/*V*_{0})^{1/3}, and *η*=3*K*_{0}′/2 − 3/2. Differentiating eq. (2) with respect to volume, we obtain the bulk modulus at reference isotherm *T*_{0}=298.15K and its pressure derivative:

The potential energy at the 298K isotherm is:

The thermal part of the Helmholtz free energy can be expressed via the Einstein model, which coincides with the Debye model at the high-temperature limit^{43}:

where Θ is the characteristic temperature depending on volume, which is related to the Debye temperature (Θ_{D}) via Θ=0.75 Θ_{D}; *n* is the number of atoms in the chemical formula of the compound, and *R* is the gas constant.

Differentiating eq. (6) with respect to temperature at constant volume one can obtain the entropy, the thermal part of the free energy and the heat capacity at constant volume:

Differentiating eq. (6) with respect to volume at constant temperature, one can obtain the thermal pressure and isothermal bulk modulus:

In eqs (10, 11) *γ* is the Grüneisen parameter, , and . Differentiation of the thermal pressure with respect to temperature at constant volume gives the pressure slope at constant volume:

The volume dependence of *γ* and *q* was accepted in the Altshuler form^{44}:

The volume dependence of the Einstein temperature can be expressed as:

In eqs (12,13,14,15) *γ*_{0} is the Grüneisen parameter at reference conditions, *γ*_{∞} is the Grüneisen parameter at infinite compression (*x*=0), and β is a fitted parameter.

The contribution of free electrons to the Helmholtz free energy in the simplest form can be expressed as^{41}:

where the parameter *e*_{0} denotes the electronic contribution to the Helmholtz free energy, *g* is an electronic analogue of the Grüneisen parameter, and *x*=*V*/*V*_{0}. The contribution of this part to entropy, internal energy, heat capacity, pressure, isothermal bulk modulus, and pressure slope at constant volume can be estimated as:

It should be emphasized that if the electronic contribution to the Helmholtz free energy is not equal to zero, the thermal Grüneisen parameter:

will be different from the Grüneisen parameter calculated from eq. (13).

The formalism for the magnetic contribution to the Helmholtz free energy was adapted from ref. ^{10} and ^{37} with modifications from ref. ^{45} to obtain the correct limit of entropy at 0K. The magnetic contribution to the Helmholtz free energy can be expressed as^{45}:

where *B*_{0} is an average magnetic moment per atom, τ=*T*/*T**, *T** is the critical temperature (*T*_{C} for ferromagnetic materials and Neel temperature *T*_{N} for paramagnetic materials). The magnetic moment^{10}^{,37} of bcc-Fe is *B*_{0}=2.22. The function *f*(τ) can be written as:

The value of the parameter *p* changes from *p*=0.4 (for bcc-Fe) to *p*=0.28 for other Fe polymorphs. However, it was shown that at 18GPa bcc-Fe transforms into hcp-Fe and the value of the magnetic moment becomes zero^{46}. Based on another study^{47}, the magnet moment of both fcc-Fe and hcp-Fe approaches zero with increasing temperature and pressure. In our study, the EoS of fcc-Fe and hcp-Fe were constrained with the magnetic moment fixed to zero and this assumption is common for thermodynamic calculations of the iron phase diagram at high pressures^{6}^{,18}^{,19}^{,20}^{,21}^{,48}^{,49}.

The equation for the Helmholtz free energy for liquid metal can be rewritten^{38}^{,39}^{,40}. The entropic parameter *a*_{S}, which characterizes residual entropy of the liquid at 0K is introduced. Accordingly, in equation (1) for liquid Fe, the magnetic part is excluded, whereas the entropy parameter is added:

The reference temperature for liquid Fe at 0.1MPa is *T*_{0}=1811K. Along the reference isotherm, the pressure was calculated using eq. (2). In the first approximation, the entropic parameter is independent of volume. Calibration of fitted parameters for liquid iron is considered below.

The full solution of the equations allowed us to find all necessary parameters for the thermodynamic description of the EoS for a solid phase. In addition, from the obtained parameters one can calculate the coefficient of thermal expansion *α*=(*P*/*T*)_{V}/*K*_{T}, heat capacity at constant pressure *C*_{P}=*C*_{V}+*α*^{2}*TVK*_{T}, and adiabatic bulk modulus *K*_{S}=*K*_{T}+*VT*(*αK*_{T})^{2}/*C*_{V}, which can also be compared with direct experimental measurements. The enthalpy and the Gibbs free energy can be found from the following relations: *H*=*E*+*PV, G*=*F*+*PV*.

The fitted parameters of the EoSs for bcc-Fe, fcc-Fe, hcp-Fe, and liquid iron (Table 1) were obtained by simultaneous optimization of the experimental measurements of the heat capacity, volume and thermal expansion, adiabatic bulk modulus at room pressure, and *P–V–T* data (all references are included in the Supplementary Information). Pressures in the *P–V–T* dataset were corrected based on the self-consistent ruby pressure scale or EoSs of Au, MgO, and W^{31}^{,35}^{,36}. Supplementary Figures S1–S4 show the temperature dependence of the isobaric heat capacity, the volume coefficient of thermal expansion, molar volume, adiabatic and isothermal bulk modulus for bcc-Fe, fcc-Fe, and hcp-Fe at 0.1MPa calculated from our EoSs in comparison with direct experimental measurements and reference data. Supplementary Figures S5–S8 show the differences between calculated pressures at different temperatures with direct experimental measurements of pressure. The pressure scales are shown in the legends to the figures.

Supplementary Figure S5 shows that our EoS for bcc-Fe is reliably consistent with experimental measurements^{49}^{,50}^{,51}, which were obtained in quasihydrostatic conditions in He or Ne pressure media. The pressure in these works was calculated using ruby^{30}^{,35} and Au^{52} pressure scales and in the pressure range of 0–15GPa these scales give consistent results. The measurements in Ne and Ar pressure media show an alternative deviation from our data^{53}. In another paper^{54}, the measurement was performed in non-hydrostatic conditions and deviates from our data at higher pressures. In ref. ^{55} and ^{56}, pressure was calculated using the NaCl EoS^{57}^{,58}. As was shown recently^{59}, this scale underestimates pressure by 0.5GPa (at 10–15GPa) in comparison with another NaCl EoS^{60}. If one recalculates the data from ref. ^{55}^{, }^{56} and ^{60} using the NaCl scale^{60}, a better consistency with our EoS for bcc-Fe is obtained, which confirms the relevance of the experimental data mentioned above^{55}^{,56}. The pressure obtained based on the compressional wave sound velocity and density measurements of bcc-Fe^{61} is also shown in Supplementary Figure S5 for comparison. The calculated thermodynamic properties (*P, T, x*=*V*/*V*_{0}, α, *S, C*_{V}, *C*_{P}, *K*_{T}, *K*_{S}, *γ*_{th}, *K*′, and Gibbs energy, *G*) for α-Fe (bcc-Fe) are listed in Supplementary Table S1 (see also bcc-Fe.xls in the Supplementary Information).

The calculated thermodynamic properties for γ-Fe (fcc-Fe) are listed in Supplementary Table S2 (see also fcc-Fe.xls in Supplementary Information). Supplementary Figure S6 shows that our EoS is reliably consistent with experimental *P–V–T* data^{26}^{,28}^{,62}^{,63} calibrated using Au, MgO, and NaCl pressure scales. At high pressures, the present EoS of fcc-Fe averages out the experimental measurements^{14}^{,17}, where pressures were calculated using the EoS of hcp-Fe^{49}. One more experimental dataset^{64} overestimates pressures by up to 6GPa compared with our data for fcc-Fe. The calculated molar volume of fcc-Fe at 0.1MPa is consistent with measurements in ref. ^{65} and ^{66} (Supplementary Figure S6). The calculated volumes of fcc-Fe^{26}^{,28} are plotted between bcc and fcc phases in Supplementary Figure S2. Their isothermal bulk modulus is consistent with our calculations (Supplementary Figure S4). The Grüneisen parameter for fcc-Fe (eq. 13) is almost independent of volume (Table 1); however, the thermal Grüneisen parameter (eq. 18) decreases with increase of temperature and pressure up to 1.5–1.6 at 100GPa and 3000–4000K (Supplementary Table S2).

The calculated thermodynamic properties for ε-Fe (hcp-Fe) are listed in Table S3 (see also hcp-Fe.xls in Supplementary data). The 298K isotherm for ε-Fe was calculated using the compressibility curve from ref. ^{49} corrected based on pressure scales^{31}^{,35}. The other parameters were calculated by optimization of the *P–V–T* data (Supplementary Figures S7, S8). Most measurements are scattered for less than ±2GPa in the 80GPa pressure range (Supplementary Figure S7). However, the more scattered data^{62} deviate by −2 to +4GPa at 298K isotherm and significantly overestimate pressures at high temperatures. It can be emphasized that in the Supplementary Figure S8 the measurements for hcp-Fe^{14} deviate in the same way, but with larger scattering at 140GPa than data for fcc-Fe in Supplementary Figure S6. Two sets of data^{13}^{,67} at 250–300GPa are significantly different (Supplementary Figures S8). This may indicate that either the pressure scales used were incorrect or significant stress was accumulated in the samples.

The calibration of the Gibbs energy of hcp-Fe was performed using the α–γ–ε triple point^{11}, the α–ε transition at room temperature^{68}^{,69}^{,70}, and considering the slope of the γ–ε transition^{14} (Fig. 1). One more important marker for calibration of the Gibbs energy for hcp-Fe and liquid Fe was the γ–ε–liquid triple point^{14}. Recently, the *P–V–T* relationships for hcp-Fe in different pressure-transmitting media at 300K up to 205GPa and at 1800K up to 100GPa have been investigated^{6}. To calculate the pressure, the Ne, NaCl-B2, and Pt pressure scales^{52}, and the MgO pressure scale^{71} were used. Supplementary Figures S7, S8 show that values obtained in ref. ^{6} are in reasonable agreement with our EoS of hcp-Fe. The authors^{6} used a third-order Birch–Murnaghan EoS to fit the parameters on the room temperature isotherm with a fixed initial density *ρ*_{0}=8.2695gcm^{−3} and obtained *K*_{0}=172.7GPa and *K*_{0}′=4.79. The pressure calculated based on these parameters is 3GPa higher than that obtained by our EoS at 50–150GPa. At higher pressures, the room temperature isotherm^{6} is consistent with our EoS of hcp-Fe (Supplementary Figure S8). At the condition of the inner core boundary (~330GPa and ~6000K), our EoS is in good agreement with the EoS obtained in ref. ^{6}. These authors determined the density deficit at the inner core boundary^{72} to be 3.6%, whereas based on our data it would be 4.4% at *T*=5882K and *P*=328.9GPa.

The position of the γ–ε–liquid triple point in the phase diagram of iron has been determined many times and revealed highly controversial results (Fig. 2). The position of the triple point was estimated at 75GPa and 2500K^{73}. This point was also placed at 100GPa and 2700K^{1}. Later, the point was shifted back to 75GPa and 2700K^{74}. Another scientific group^{75} placed the γ–ε–liquid triple point at 60GPa and 2800K using pressures calculated from the EoS of hcp-Fe^{12} and EoS of MgO^{76}. The melting line of Fe determined in ref. ^{77} is located at a higher temperature than most determined triple points (Fig. 2). The pressure in this work was estimated before laser heating using the old ruby scale^{78}. Melting of Fe was also determined at 105GPa and 3510K^{79} with pressures estimated by the EoS of hcp-Fe^{12} before laser heating. Estimations in ref. ^{17} revealed the triple point at 88GPa and 2800K. In addition, the melting temperature of Fe at the core–mantle boundary (*P*=135GPa) at 3500±100K with careful estimation of thermal pressure and anharmonicity was determined^{80}. A novel method was presented for detecting the solid–liquid phase boundary of compressed iron at high temperatures using synchrotron Mössbauer spectroscopy^{81}. The melting points are shown in Fig. 2 and the pressure was determined using the ruby scale^{82} with corrections for thermal pressure (note that at Fig. 2 their data are shown without these corrections). Finally, the melting line of Fe measured using a synchrotron-based fast X-ray diffraction method was also proposed^{14}. The authors fixed the triple point at *P*=98.5GPa and *T*=3712K. The pressures were estimated from the EoS of hcp-Fe^{49} and KCl^{83}. Thus, direct experimental measurements of the melting line of Fe and the γ–ε–liquid triple point indicate about 40GPa uncertainty in pressure and about 1000K uncertainty in temperature. The triple point calculated from our data is located at *P*=106.5GPa and *T*=3787K.

The EoS of liquid Fe was estimated using the following approach (eq. 22). Standard conditions for liquid Fe were chosen at *T*_{0}=1811K and *P*=0.1MPa. For these conditions, the following parameters were recommended^{84}: *V*_{0}=7.957cm^{3}mol^{−1}, *K*_{S0}=109.5GPa, *α*(*V*)=92E–6K^{−1}, *γ*_{0}=1.735, and *C*_{P}=46.632Jmol^{−1}K^{−1}. Fitting parameters for our EoS for liquid iron were calculated using the following assumptions (as a starting point for calculations): (i) the Gibbs energy of α-Fe and liq-Fe must be the same at 1811K and 0.1MPa; (ii) the entropy of liquid Fe increases by 7.626Jmol^{−1}K^{−1} upon melting^{9} in comparison with the entropy of α-Fe; (iii) the melting line of Fe was taken from ref. ^{14} as more accurate. The fitted parameters of EoS for liquid Fe are listed in Table 1. The calculated thermodynamic properties for liquid Fe are listed in Supplementary Table S4 (see also liquid-Fe.xls in Supplementary Information). The calculated melting lines of iron from the bcc, fcc, and hcp phases are shown in Figs 1, ,22 and and33.

The thermodynamic properties of liquid iron at a pressure of 0.1MPa calculated from our EoS (Supplementary Table S4) are in good agreement with the measured values of density^{84}, the sound velocity (*v*_{P}) and adiabatic bulk modulus^{84}^{,85}^{,86}. The calculated entropy at pressure 0.1MPa is very close to the reference data^{9}.

We applied the new EoS data for all Fe phases to calculate the phase diagram up to 350GPa (Figs 1, ,22 and and3).3). Our melting line of Fe is in close agreement with that from ref. ^{14} because our EoS of liquid iron is based mainly on these data. The calculated triple points have the following parameters: α–γ–liquid is located at 5.2GPa and 1998K, α–γ–ε at 7.3GPa and 820K, and γ–ε–liquid at 106.5GPa and 3787K.

Most of the previous estimations plot the melting temperature of iron below the melting line obtained from our EoSs (Fig. 3). For example estimations in refs ^{1} and ^{18} indicate melting temperatures of 5000K and 5600±200K at 330GPa, respectively. The melting line of Fe in ref. ^{19} was calculated using the revised thermodynamic properties of Fe. Their triple point γ–ε–liquid is located at 81GPa and 3200K and at 330GPa the melting line is located at 5400K. Another estimation indicates a melting temperature of 6100K at 330GPa^{21}. The newly measured melting line of Fe extrapolated to 330GPa reveals a temperature of about 6230±500 K^{14}. Most of the *ab initio* estimations provide higher melting temperatures of Fe at 330 GPa compared with our data: 6700±600K^{87}, 6370±100K^{88}, 6900±400K^{89}, 7100–7200K^{90}, 6325K^{91}, and 6345K^{92}.

Our thermodynamic calculations indicate a melting temperature of hcp-Fe at 5882K at 328.9GPa, and it is in a good agreement with the estimated melting temperature *T*=5800±500K obtained from the shock wave data^{93} and extrapolation from the experiments^{14}, which is *T*=6230±500K. Figures 2 and and33 show the fcc–liquid melting curve and its extrapolation to high pressures. The melting curve was found to be very close to the hcp–liquid curve. At the conditions near the fcc–hcp–liquid triple point the Clapeyron slope of the fcc–liquid curve is d*T*/d*P*=12.8K/GPa while the slope of the hcp–liquid curve is higher (d*T*/d*P*=13.7K/GPa). Therefore, the hcp–liquid curve overlaps the metastable fcc–liquid curve at a pressure of ≈160GPa. At higher pressures, the fcc phase is in the stability field of the liquid. These transitions can also be confirmed by the temperature dependence of the Gibbs energy of the fcc, hcp, and liquid iron at 125GPa (Fig. 4). At 125GPa the fcc-Fe melts at 4019K, however, this point is in the stability field of the hcp-Fe. The hcp-Fe melts at 4033K and at higher temperatures the liquid is the only stable phase.

Meanwhile, this assumption is made based on the equilibrium thermodynamics, the shock data^{93} on the iron shows two phase transitions on the Hugoniot curve at pressures from 77 to 400GPa. A discontinuity in sound velocities of iron at 200GPa may mark the transition of hcp-iron to fcc iron. The authors described the second discontinuity at 243GPa as melting of the fcc-Fe. The study^{94} based on new measurements of sound velocity and reassessment of previously obtained shock data established the melting temperature of the iron between 5100±500K at *P*=225±3GPa and 6100±500K and *P*=260±3GPa on the Hugoniot curve. According to the small temperature difference between the melting lines of hcp–liquid and fcc–liquid, there might be a possible fcc–hcp transition in a very narrow temperature range. The following transition can be caused by the kinetics of the transition process or the effect of fcc-Fe and hcp-Fe magnetic moment.

The magnetic moment of the hcp-Fe is *В*_{0}=0.28 at ambient conditions^{10}. At high pressure and temperature, iron was found to be paramagnetic^{47}. However, under the Earth’s core conditions (*P* ≈ 360GPa, *T* ≈ 6000K) Fe acquires a substantial local magnetic moment^{47}, up to 1.3μ_{B}. The authors used a microscopic phenomenological model for longitudinal spin fluctuations (LSFs) based on density functional theory calculations. In our model (eq. 19) it means that *B*_{0}=0.3 and this value seems to be too high for the EoS of hcp-Fe. If *B*_{0}=0.14, then melting of the hcp-Fe is observed at 7150K and 328.9GPa and is in agreement with the computer calculations^{90}. Whereas the melting curve goes 500K higher than the experimentally determined one^{14} at 100–200GPa. The lower value of the hcp-Fe magnetic moment seems to be more realistic. If the magnetic moment for hcp-Fe is *В*_{0}=0.04 and *T*_{N}=67K (see eqs 19,20,21), then the hcp–liquid curve shifts to higher temperatures (red squares on Fig. 3) and crosses the melting curve of fcc-Fe at a temperature of ≈5300K and pressure of ≈230GPa. This result is consistent with the data from the shock wave experiments^{93}^{,94}. The metastable line bcc–hcp-iron was also calculated at high pressures. The extrapolation of this line from the low-pressure data (Fig. 1) is presented in Fig. 3. The bcc–hcp line is located in the stability field of either the fcc-Fe or liquid at high pressures. Thus, the bcc-Fe could not be stable at high pressures based on our thermodynamic model.

Recently, thermodynamic functions for fcc, hcp, and liquid iron were tabulated^{20}. It is important to compare the thermodynamics of these phases calculated by different methods. The thermodynamic formalism in ref. ^{18} and ^{20} is based on the Gibbs free energy calculations depending on temperature and pressure. First, two partial derivatives of the Gibbs energy allow calculation of the entropy and volume . The Helmholtz free energy depends on temperature and volume; thus, its derivatives yield entropy and pressure . Figures S9–S11 show comparisons of the thermodynamic functions calculated by these two methods for fcc, hcp, and liquid iron. These figures show that the entropy and molar volume of fcc- and hcp-Fe are consistent, especially at the conditions of their stability field. At the same time, our EoSs cannot be extrapolated to very high temperature at 0.1MPa. At moderate and high pressures, such extrapolation provides reliable results with a close agreement with reference data^{20}. At 0.1MPa the calculated entropy^{20} of solid phases of Fe decreases much faster than in our EoSs. The volumes for liquid iron from our EoS and that from ref. ^{20} are very consistent; however, entropy is consistent only at temperatures close to the melting line (Supplementary Figure S11).

Supplementary Figure S12 shows isochores of liquid Fe calculated from our EoS and plotted in comparison with *P–V–T* data calculated by *ab initio* molecular dynamics simulations^{95}^{,96}. The comparison with earlier *ab initio* calculations can be found in the cited papers.

Comparisons of other thermodynamic functions, which are second derivatives of the Gibbs energy or the Helmholtz free energy, are not necessary. The values of these functions can be found in Tables S1–S4. If entropy and molar volume (pressure) are comparable in different EoSs, then the Gibbs energy will be similar and deviations in calculated lines in equilibrium phase diagram would be negligible. The differences would be defined by reference points chosen by different authors.

Density (*ρ* at given *P* and *T*), adiabatic bulk modulus (*K*_{S}) and the P-wave velocity (*v*_{P}=(*K*_{S}/ρ)^{1/2}) of the liquid iron estimated by our formalism can be compared with the values of PREM^{72}. The following comparison has also been assessed by previous studies^{84}^{,95} on temperature calculations of liquid iron in the Earth’s core at isentropic conditions and various pressures. The initial temperature in our model is 5882K at the inner core boundary (ICB) (328.9GPa) (Fig. 3). At these conditions, the entropy for liquid iron, *S*=114.15Jmol^{−1}K^{−1} and hcp-Fe, *S*=107.2Jmol^{−1}K^{−1}, were determined from eqs 7 and 17. Under isentropic conditions, the entropy at the pressures of the mantle–core boundary (MCB, 135.8GPa) is equal to the entropy at the ICB. Figure 5a shows the calculated isentropic temperature profile started from ICB conditions for liquid iron. The calculated isentropic temperature profile for hcp-iron and liquid iron started from 6000K based on the reference data^{95} are also presented in Fig. 5a. The Clapeyron slope (*T*/*P*)_{S} for hcp-Fe differs significantly for both the liquid Fe and the calculated model^{95}. This inconsistency is probably due to the differences in the EoSs.

Figure 5b shows adiabatic bulk modulus and density for solid and liquid iron calculated from our EoSs in comparison with PREM^{72} and the calculated model^{95}. The deviations from the PREM are presented as (*K*_{S PREM} − *K*_{S cal})/*K*_{S PREM}) × 100. The adiabatic bulk modulus for liquid Fe is about 3.0–4.5% lower than the PREM and the calculated^{95}
*K*_{S} is 3.1–9.0% higher than the PREM (Fig. 5b). *K*_{S} for hcp-Fe is almost identical to the PREM at the Earth’s inner core conditions. The calculated isentrope for liquid iron from 5000K to 8000K is presented in ref. ^{84}. Based on these data, the adiabatic bulk modulus at the ICB conditions is consistent with isentropes of 7000K and at the ICB with isentropes of 5000K (see Fig. 10 in ref. ^{84}). The liquid iron density calculated in our study is 7.6–8.2% higher than PREM and the calculated one^{95} is 8.9–7.7% higher than PREM (Fig. 5c). Figure 5d shows that the calculated P-wave velocity for liquid iron is 5.4–5.7% lower than PREM and consistent with the calculated estimations^{84}. The calculations^{95} indicate that the P-wave velocity is 2.5% lower than PREM at the ICB conditions and very close to PREM at the CMB.

The EoSs of the liquid and solid iron phases were constrained based on the Helmholtz free energy. The obtained EoSs allow calculation of *P–V–T* relations and thermodynamic properties of Fe at the Earth’s core conditions (up to 6000K and 350GPa). The calculated triple points have the following parameters: bcc–fcc–hcp is located at 7.3GPa and 820K, bcc–fcc–liquid at 5.2GPa and 1998K, and fcc–hcp–liquid at 106.5GPa and 3787K. The melting temperature is 5882К at ICB pressure (328.9GPa). An extrapolation to high-pressure conditions revealed the very close position of the melting curves of fcc- and hcp-Fe. If the magnetic moment for hcp-Fe is *В*_{0}=0.04 and *T*_{N}=67K (see eqs 19,20,21), the hcp–liquid curve shifts to higher temperatures and crosses the fcc–liquid curve at the temperature of ≈5300K and pressure of ≈230GPa. This estimation is consistent with the shock wave data^{93}^{,94}. At higher pressure conditions, the metastable bcc–hcp curve is in the fcc-Fe or liquid stability field. The iron density, adiabatic bulk modulus and P-wave velocity calculated up to 328.9GPa at adiabatic temperature conditions started from 5882K (outer/inner core boundary) were compared with the PREM and calculated model^{95}. We determined the density deficit of hcp-Fe at the inner core boundary (*T*=5882K and *P*=328.9GPa) to be 4.4%.

**How to cite this article**: Dorogokupets, P. I. *et al*. Thermodynamics and Equations of State of Iron to 350 GPa and 6000 K. *Sci. Rep.*
**7**, 41863; doi: 10.1038/srep41863 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Click here to view.^{(2.3M, doc)}

Click here to view.^{(362K, xls)}

Click here to view.^{(310K, xls)}

Click here to view.^{(284K, xls)}

Click here to view.^{(212K, xls)}

We thank the reviewers for critical comments and suggestions. This work was financially supported by Russian Science Foundation (No 14-17-00601) and performed under program of Ministry of education and science of Russian Federation (Project No 14.B25.31.0032).

The authors declare no competing financial interests.

**Author Contributions** P.D. initiated the project. P.D., T.S., A.D. and K.L. analyzed the data and wrote the manuscript text. A.D., T.S. and K.L. reviewed the manuscript.

- Boehler R. Temperatures in the Earth’s core from melting point measurements of iron at high static pressures. Nature 363, 534–536 (1993).
- Funtikov A. Phase diagram of iron: Implications for the state of the Earth’s core. Izv. Phys. Solid Earth 36, 958–964 (2000).
- Funtikov A. Phase diagram and melting curve of iron obtained using the data of static and shock-wave measurements. High Temp. 41, 850–864 (2003).
- Li J. & Fei Y. In Treatise on Geochemistry Vol. 2 Experimental constraints on core composition (eds Holland H. D. & Turekian K. K.) 521–546 (Elsevier – Pergamon, Oxford, 2003).
- Hirose K., Labrosse S. & Hernlund J.
Composition and state of the core.
*Annual Rev*. Earth Planet. Sciences 41, 657–691 (2013). - Fei Y., Murphy C., Shibazaki Y., Shahar A. & Huang H. Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth’s solid inner core. Geophys. Res. Lett. 43, 6837–6843 (2016).
- Swartzendruber L. The Fe (Iron) System. J. Phase Equilibria 3, 161–165 (1982).
- Tonkov E. Y. & Ponyatovsky E. G. In Advances in Metallic Alloys V.4 Phase transformations of elements under high pressure (eds Fridlyander J. N. & Eskin D. G.) 377 (CRC Press Boca Raton, Florida, 2005).
- Desai P. D. Thermodynamic properties of iron and silicon. J. Phys. Chem. Ref. Data 15, 967–983 (1986).
- Dinsdale A. T. SGTE data for pure elements. Calphad 15, 317–425 (1991).
- Klotz S., Le Godec Y., Strassle T. & Stuhr U. The α–γ–ε triple point of iron investigated by high pressure–high temperature neutron scattering. App. Phys. Lett. 93, 091904 (2008).
- Mao H., Wu Y., Chen L., Shu J. & Jephcoat A. P. Static compression of iron to 300 GPa and Fe0.8Ni0.2 alloy to 260 GPa: Implications for composition of the core. J. Geophys. Res. 95, 21737–21742 (1990).
- Tateno S., Hirose K., Ohishi Y. & Tatsumi Y. The structure of iron in Earth’s inner core. Science 330, 359–361 (2010). [PubMed]
- Anzellini S., Dewaele A., Mezouar M., Loubeyre P. & Morard G. Melting of iron at Earth’s inner core boundary based on fast X-ray diffraction. Science 340, 464–466 (2013). [PubMed]
- Belonoshko A. B., Ahuja R. & Johansson B. Stability of the body-centred-cubic phase of iron in the Earth’s inner core. Nature 424, 1032–1034 (2003). [PubMed]
- Dubrovinsky L. et al. . Body-centered cubic iron-nickel alloy in Earth’s core. Science 316, 1880–1883 (2007). [PubMed]
- Komabayashi T., Fei Y., Meng Y. & Prakapenka V.
*In-situ*X-ray diffraction measurements of the γ-ε transition boundary of iron in an internally-heated diamond anvil cell. Earth Planet. Science Lett. 282, 252–257 (2009). - Komabayashi T. & Fei Y. W. Internally consistent thermodynamic database for iron to the Earth’s core conditions. J. Geophys. Res. 115, B03202 (2010).
- Fei Y. & Brosh E. Experimental study and thermodynamic calculations of phase relations in the Fe–C system at high pressure. Earth Planet. Science Lett. 408, 155–162 (2014).
- Komabayashi T. Thermodynamics of melting relations in the system Fe-FeO at high pressure: Implications for oxygen in the Earth’s core. J. Geophys. Res. 119, 4164–4177 (2014).
- Saxena S. K. & Eriksson G. Thermodynamics of iron at extreme pressures and temperatures. J. Phys. Chem. Solids 84, 70–74 (2015).
- Lu X.-G., Selleby M. & Sundman B. Assessments of molar volume and thermal expansion for selected bcc, fcc and hcp metallic elements. Calphad 29, 68–89 (2005).
- Anderson O. L. & Isaak D. G. The dependence of the Anderson-Grüneisen parameter δT upon compression at extreme conditions. J. Phys. Chem. Solids 54, 221–227 (1993).
- Medvedev A. B. Wide-range multiphase equation of state for iron. Combustion, Explosion, and Shock Waves 50, 582–598 (2014).
- Sakai T., Ohtani E., Hirao N. & Ohishi Y. Stability field of the hcp-structure for Fe, Fe-Ni, and Fe-Ni-Si alloys up to 3 Mbar. Geophys. Res. Lett. 38 (2011).
- Nishihara Y. et al. . Isothermal compression of face-centered cubic iron. Amer. Miner. 97, 1417–1420 (2012).
- Yamazaki D. et al. . P-V-T equation of state for epsilon-iron up to 80GPa and 1900K using the Kawai-type high pressure apparatus equipped with sintered diamond anvils. Geophys. Res. Lett. 39, L20308 (2012).
- Tsujino N. et al. . Equation of state of γ-Fe: Reference density for planetary cores. Earth Planet. Science Lett. 375, 244–253 (2013).
- Dorogokupets P. I. & Dewaele A. Equations of state of MgO, Au, Pt, NaCl-B1, and NaCl-B2: Internally consistent high-temperature pressure scales. High Press. Res. 27, 431–446 (2007).
- Dorogokupets P. I. & Oganov A. R. Ruby, metals, and MgO as alternative pressure scales: A semiempirical description of shock-wave, ultrasonic, x-ray, and thermochemical data at high temperatures and pressures. Phys. Rev. B 75, 024115 (2007).
- Dorogokupets P. I., Sokolova T. S., Danilov B. S. & Litasov K. D. Near-absolute equations of state of diamond, Ag, Al, Au, Cu, Mo, Nb, Pt, Ta, and W for quasi-hydrostatic conditions. Geodyn. Tectonophys. 3, 129–166 (2012).
- Dorogokupets P. I., Sokolova T. S. & Litasov K. D. Thermodynamic properties of bcc-Fe to melting temperature and pressure to 15 GPa. Geodyn. Tectonophys. 5, 1033–1044 (2014).
- Dorogokupets P. I., Dymshits A. M., Sokolova T. S., Danilov B. S. & Litasov K. D.
The equations of state of forsterite, wadsleyite, ringwoodite, akimotoite, MgSiO
_{3}-perovskite and postperovskite and phase diagram of the Mg_{2}SiO_{4}system at pressures to 130 GPa. Russ. Geol. Geophys. 56, 172–189 (2015). - Dorogokupets P. I., Sokolova T. S., Dymshits A. M. & Litasov K. D.
Thermodynamic properties of rock-forming oxides, α-Al
_{2}O_{3}, Cr_{2}O_{3}, α-Fe_{2}O_{3}, and Fe_{3}O_{4}at high temperatures and pressures. Geodyn. Tectonophys. 7, 459–476 (2016). - Sokolova T. S., Dorogokupets P. I. & Litasov K. D. Self-consistent pressure scales based on the equations of state for ruby, diamond, MgO, B2-NaCl, as well as Au, Pt and other metals to 4 Mbars and 3000 K. Russ. Geol. Geophys. 54, 181–199 (2013).
- Sokolova T. S., Dorogokupets P. I., Dymshits A. M., Danilov B. S. & Litasov K. D. Microsoft excel spreadsheets for calculation of P-V-T relations and thermodynamic properties from equations of state of MgO, diamond and nine metals as pressure markers in high-pressure and high-temperature experiments. Comput. Geoscience 94, 162–169 (2016).
- Hillert M. & Jarl M. A model for alloying in ferromagnetic metals. Calphad 2, 227–238 (1978).
- Molodets A., Molodets M. & Nabatov S. Helmholtz’free energy of molten metals. High Temp. 36, 891–896 (1998).
- Molodets A., Molodets M. & Nabatov S. Free energy of liquid diamond. Combustion, Explosion and Shock Waves 35, 185–190 (1999).
- Molodets A. Thermodynamic potentials and non-monotonic melting curve of sodium at high pressure. High Press. Res. 30, 325–331 (2010).
- Zharkov V. N. & Kalinin V. A. Equations of State of Solids at High Pressures and Temperatures (ed. Tybulewicz A.) 257 (Springer Science+Business Media, New York, 1971).
- Vinet P., Ferrante J., Rose J. H. & Smith J. R. Compressibility of solids. J. Geophys. Res. 92, 9319–9325 (1987).
- Dorogokupets P. I. P-V-T Equations of state of MgO and thermodynamics. Phys. Chem. Miner. 37, 677–684 (2010).
- Al’tshuler L. V., Brusnikin S. E. & Kuz’menkov E. A. Isotherms and Gruneisen functions for 25 metals. J. Appl. Mechan. Technic. Phys. 28, 129–141 (1987).
- Jacobs M. H. & Schmid-Fetzer R. Thermodynamic properties and equation of state of fcc aluminum and bcc iron, derived from a lattice vibrational method. Phys. Chem. Miner. 37, 721–739 (2010).
- Iota V. et al. . Electronic structure and magnetism in compressed 3d transition metals. Appl. Phys. Lett. 90, 042505 (2007).
- Ruban A. V., Belonoshko A. B. & Skorodumova N. V. Impact of magnetism on Fe under Earth’s core conditions. Phys. Rev. B 87, 014405 (2013).
- Brosh E., Makov G. & Shneck R. Z. Application of CALPHAD to high pressures.
*Calphad*31, 173–185 (2007). - Dewaele A. et al. . Quasihydrostatic equation of state of iron above 2 Mbar. Phys. Rev. Lett. 97, 215504 (2006). [PubMed]
- Dewaele A. & Garbarino G. Low temperature equation of state of epsilon-iron. Report No. HC-1679 (ESRF, Grenoble, 2014).
- Liu J., Lin J.-F., Alatas A. & Bi W. Sound velocities of bcc-Fe and Fe0.85Si0.15 alloy at high pressure and temperature. Phys. Earth Planet. Inter. 233, 24–32 (2014).
- Fei Y. et al. . Toward an internally consistent pressure scale. Proc. Natl. Acad. Sci. USA 104, 9182–9186 (2007). [PubMed]
- Jephcoat A. P., Mao H. K. & Bell P. M. Static compression of iron to 78 GPa with rare gas solids as pressure-transmitting media. J. Geophys. Res. 91, 4677–4684 (1986).
- Huang E., Bassett W. A. & Tao P. Pressure-temperature-volume relationship for hexagonal close packed iron determined by synchrotron radiation. J. Geophys. Res. 92, 8129–8135 (1987).
- Mao H. K., Bassett W. A. & Takahashi T. Effect of pressure on crystal structure and lattice parameters of iron up to 300 kbar. J. Appl. Phys. 38, 272–276 (1967).
- Zhang J. & Guyot F. Thermal equation of state of iron and Fe0.91Si0.09. Phys. Chem. Miner. 26, 206–211 (1999).
- Decker D. L. Equation of State of NaCl and Its Use as a Pressure Gauge in High‐Pressure Research. J. Appl. Phys. 36, 157–161 (1965).
- Decker D. L. High-pressure equation of state for NaCl, KCl, and CsCl. J. Appl. Phys. 42, 3239–3244 (1973).
- Strässle T., Klotz S., Kunc K., Pomjakushin V. & White J. S. Equation of state of lead from high-pressure neutron diffraction up to 8.9 GPa and its implication for the NaCl pressure scale. Phys. Rev. B. 014101 (2014).
- Brown M. J. The NaCl pressure standard. J. Appl. Phys. 86, 5801–5808 (1999).
- Antonangeli D. et al. . Toward a mineral physics reference model for the Moon’s core. PNAS 115, 3916–3919 (2015). [PubMed]
- Funamori N., Yagi T. & Uchida T.
High-pressure and high-temperature
*in situ*x-ray diffraction study of iron to above 30 GPa using MA8-type apparatus. Geophys. Res. Lett. 23, 953–956 (1996). - Campbell A. J. et al. . High pressure effects on the iron–iron oxide and nickel–nickel oxide oxygen fugacity buffers. Earth Planet. Science Lett. 286, 556–564 (2009).
- Boehler R., Vonbargen N. & Chopelas A. Melting, thermal-expansion, and phase-transitions of iron at high-pressures. J. Geophys. Res. 95, 21731–21736 (1990).
- Basinski Z. S., Hume-Rothery W. & Sutton A. L. The lattice expansion of Iron. Proceedings of the Royal Society of London A 229, 459–467 (1955).
- Kohlhaas R., Dunner P. & Schmitz P. N. Uber die temperaturabhangigkeit der gitterparameter von Eisen, Kobalt und Nickel im bereich hoher temperaturen. Zeitschrift fur Angewandte Physik 23, 245–249 (1967).
- Sakai T. et al. . Equation of state of pure iron and Fe0.9Ni0.1 alloy up to 3 Mbar. Phys. Earth Planet. Inter. 228, 114–126 (2014).
- Ishimatsu N. et al. . a-e transition pathway of iron under quasihydrostatic pressure conditions. Phys. Rev. B 90, 014422 (2014).
- Dewaele A. et al. . Mechanism of the α-d phase transformation in iron. Phys. Rev. B 91, 174105 (2015).
- Zarkevich N. A. & Johnson D. D. Coexistence pressure for a martensitic transformation from theory and experiment: Revisiting the bcc-hcp transition of iron under pressure. Phys. Rev. B 91, 174194 (2015).
- Tange Y., Nishihara Y. & Tsuchiya T. Unified analyses for P-V-T equation of state of MgO: A solution for pressure-scale problems in high P-T experiments. J. Geophys. Res. 114, B03208 (2009).
- Dziewonski A. M. & Anderson D. L. Preliminary reference Earth model. Phys. Earth Planet. Inter. 25, 297–356 (1981).
- Boehler R. The phase diagram of iron to 430 kbar. Geophys. Res. Lett. 13, 1153–1156 (1986).
- Boehler R., Santamaría-Pérez D., Errandonea D. & Mezouar M. Melting, density, and anisotropy of iron at core conditions: new X-ray measurements to 150 GPa. J. Phys. 121, 022018 (2008).
- Shen G., Mao H.-k., Hemley R. J., Duffy T. S. & Rivers M. L. Melting and crystal structure of iron at high pressures and temperatures. Geophys. Res. Lett. 25, 373–376 (1998).
- Jamieson J. C., Fritz J. N. & Manghnani M. H. In High Pressure Research in Geophysics Pressure measurement at high temperature in X-ray diffraction studies: gold as a primary standard (eds Akimoto S. & Manghnani M. H.) 27–48 (Center for Academic Publications, Tokyo, 1982).
- Williams Q., Jeanloz R., Bass J., Svendsen B. & Ahrens T. J. The melting curve of iron to 250 gigapascals: A constraint on the temperature at Earth’s center. Science 236, 181–182 (1987). [PubMed]
- Mao H. K., Bell P. M., Shaner J. W. & Steinberg D. J. Specific volume measurements of Cu, Mo, Pd, and Ag and calibration of the ruby R1 fluorescence pressure gauge from 0.06 to 1 Mbar. J. Appl. Phys. 49, 3276–3283 (1978).
- Ma Y.
et al. .
*In situ*X-ray diffraction studies of iron to Earth-core conditions. Phys. Earth Planet. Inter. 143–144, 455–467 (2004). - Murphy C. A., Jackson J. M., Sturhahn W. & Chen B. Melting and thermal pressure of hcp-Fe from the phonon density of states. Phys. Earth Planet. Inter. 188, 114–120 (2011).
- Jackson J. M. et al. . Melting of compressed iron by monitoring atomic dynamics. Earth Planet. Science Lett. 362, 143–150 (2013).
- Mao H., Xu J.-A. & Bell P. Calibration of the ruby pressure gauge to 800 kbar under quasi-hydrostatic conditions. J. Geophys. Res. 91, 4673–4676 (1986).
- Dewaele A. et al. . High-pressure–high-temperature equation of state of KCl and KBr. Phys. Rev. B 85, 214105 (2012).
- Anderson W. W. & Ahrens T. J. An equation of state for liquid iron and implications for the Earths core. J. Geophys. Res. 99, 4273–4284 (1994).
- Jing Z. et al. . Sound velocity of Fe–S liquids at high pressure: implications for the Moon’s molten outer core. Earth Planet. Science Lett. 396, 78–87 (2014).
- Kuwabara S. Sound velocity and elastic properties of Fe–Ni and Fe–Ni–C liquids at high pressure. Phys. Chem. Miner. 43, 229–236 (2016).
- Alfe D., Gillan M. J. & Price G. D. The melting curve of iron at the pressures of the Earth’s core from ab initio calculations. Nature 401, 462–464 (1999).
- Alfe D. Temperature of the inner-core boundary of the Earth: Melting of iron at high pressure from first-principles coexistence simulations. Phys. Rev. B 79, 060101 (2009).
- Sola E. & Alfe D. Melting of Iron under Earth’s Core Conditions from Diffusion Monte Carlo Free Energy Calculations. Phys. Rev. Lett. 103, 078501 (2009). [PubMed]
- Belonoshko A. B., Ahuja R. & Johansson B. Quasi - Ab initio molecular dynamic study of Fe melting. Phys. Rev. Lett. 84, 3638–3641 (2000). [PubMed]
- Belashchenko D. K. Estimation of the thermodynamic characteristics of the Earth’s core using the embedded atom model. Geochem. Internat. 52, 456–466 (2014).
- Zhang W.-J., Liu Z.-Y., Liu Z.-L. & Cai L.-C. Melting curves and entropy of melting of iron under Earth’s core conditions. Phys. Earth Planet. Inter. 244, 69–77 (2015).
- Brown J. M. & McQueen R. G. Phase transitions, Grüneisen parameter, and elasticity for shocked iron between 77 GPa and 400 GPa. J. Geophys. Res. 91, 7485–7494 (1986).
- Nguyen J. H. & Holmes N. C. Melting of iron at the physical conditions of the Earth’s core. Nature 427, 339–342 (2004). [PubMed]
- Ichikawa H., Tsuchiya T. & Tange Y. The P-V-T equation of state and thermodynamic properties of liquid iron. J. Geophys. Res. 119 (2014).
- Umemoto K. et al. . Liquid iron-sulfur alloys at outer core conditions by first-principles calculations. Geophys. Res. Lett. 41, 2014GL061233 (2014).
- Zhang D. et al. . Temperature of Earth’s core constrained from melting of Fe and Fe0.9Ni0.1 at high pressures. Earth Planet. Science Lett. 447, 72–83 (2016).
- Aitta A. Iron melting curve with a tricritical point. J. Statist. Mechanic. 12, P12015 (2006).
- Alfè D., Vočadlo L., Price G. D. & Gillan M. J. Melting curve of materials: theory versus experiments. J. Phys. Condens. Matter. 16, S973–S982 (2004).

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |