We investigate the equation of state and elastic properties of hcp iron at high pressures and high temperatures using first principles linear response linear-muffin-tin-orbital method in the generalized-gradient approximation. We calculate the Helmholtz free energy as a function of volume, temperature, and volume-conserving strains, including the electronic excitation contributions from band structures and lattice vibrational contributions from quasi-harmonic lattice dynamics. We perform detailed investigations
IntroductionIron is one of the most abundant elements in the Earth, and is fundamental to our world. The study of iron at high pressures and high temperatures is of great geophysical interest, since both the Earth's liquid outer core and solid inner core are composed mostly of this element. Although the crystal structure of iron at the extremely high temperature (4000 to 8000 K) and high pressure (330 to 360 GPa) conditions found in the inner core is still under intensive debate, 1-10 the hexagonal-close-packed phase (ε-Fe) is commonly believed to have a wide stability field extending from deep mantle to core conditions, and serves as a starting point for understanding the nature of the inner core. 11 Significant experimental and theoretical efforts have been recently devoted to investigate various properties of hcp iron at high pressures and high temperatures. New high-pressure diamond-anvil-cell techniques have been developed or significantly improved, which makes it possible to reach higher pressures and provide more valuable information on material properties in these extreme states. First-principles based theoretical techniques have been improved in reliability and accuracy, and have been widely used to predicate the high pressure-temperature behavior and provide fundamental understandings to the experiment.Despite intensive investigations, numerous fundamental problems remain unresolved, and many of the current results are mutually inconsistent. 11 The melting line at very high pressures has been one of the most difficult and controversial problems.
12-19Other major problems include possible subsolidus phase transitions 2,4,5,11,20 and the magnetic structure of the dense hexagonal iron. In section II we detail the theoretical methods to perform the first-principles calculations and obtain the thermal properties and elastic moduli. We present the results and related discussions about the thermal equation of state in section III, and about the thermoelasticity in section IV. We conclude with a brief summary in Section V.
II. Theoretical methodsThe Helmholtz free energy F for many metals has three major contributionswith V as the volume, T as the temperature, and δ as the strain. E static is the zerotemperature energy of a static lattice, F el is the thermal free energy arising from electronic excitations, and F ph is the lattice vibrational energy contribution. We obtain both E static and F el from first-principles calculations directly, assuming that the eigenvalues for given lattice and nuclear pos...