AN EQUATION OF STATE FOR LIQUID-IRON AND IMPLICATIONS FOR THE EARTHS CORE

An equation of state is presented for liquid iron based on published ultrasonic, thermal expansion, and enthalpy data at 1 bar and on pulse-heating and shock wave compression and sound speed data up to 10 Mbar. The equation of state parameters, centered at 1 bar and 1811 K (the normal melting point of iron), are density, rho0 = 7019 kg/m3, isentropic bulk modulus, K(S0) = 109.7 GPa, and the first-and second-pressure derivatives of K(S), K(S0)' = 4.66 and K(S0)'' = -0.043 GPa-1. A parameterization of the Gruneisen parameter gamma as a function of density rho and specific internal energy E is gamma = gamma0 + gamma'(rho/rho0)n(E - E0) where gamma0 = 1.735, gamma' = -0.130 kg/MJ, n = -1.87, and E0 is the internal energy of the liquid at 1 bar and 1811 K. The model gives the temperature dependence of gamma at constant volume as (partial derivative gamma/partial derivative T)V\1bar,1811K = -8.4 x 10(-5) K-1. The constant volume specific heat of liquid Fe at core conditions is 4.0-4.5 R. The model gives excellent agreement with measured temperatures of Fe under shock compression. Comparison with a preliminary reference Earth model indicates that the light component of the core does not significantly affect the magnitude of the isentropic bulk modulus of liquid Fe but does decrease its pressure derivative by approximately 10%. Pure liquid Fe is 3-6% more dense than the inner core, supporting the presence of several percent of light elements in the inner core.