High pressure melt line of nickel using a generalized embedded atomic method potential

Abstract

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.