Thermodynamics and Equations of State of Iron to 350 GPa and 6000 K

Abstract

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.3 GPa and 820 K, bcc-fcc-liquid at 5.2 GPa and 1998 K, and fcc-hcp-liquid at 106.5 GPa and 3787 K. At conditions near the fcc-hcp-liquid triple point, the Clapeyron slope of the fcc-liquid curve is dT/dP = 12.8 K/GPa while the slope of the hcp-liquid curve is higher (dT/dP = 13.7 K/GPa). Therefore, the hcp-liquid curve overlaps the metastable fcc-liquid curve at pressures of about 160 GPa. 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.9 GPa at adiabatic temperature conditions started from 5882 K (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 = 5882 K and P = 328.9 GPa) to be 4.4%.