• Citation: Y. Sun, M.I. Mendelev, F. Zhang, X. Liu, B. Da, C.-Z. Wang, R.M. Wentzcovitch, and K.-M. Ho (2024), "Unveiling the effect of Ni on the formation and structure of Earth’s inner core", Proceedings of the National Academy of Sciences, 121(4), e2316477121. DOI: 10.1073/pnas.2316477121.
    Abstract: Ni is the second most abundant element in the Earth’s core. Yet, its effects on the inner core’s structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe’s melting point at inner core pressures. The melting temperature of Ni is shown to be 700 to 800 K higher than that of Fe at 323 to 360 GPa. hcp, bcc, and liquid phase relations differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe’s crystallization at core pressures. These results suggest that Ni may substantially impact the structure and formation process of the solid inner core.

    Notes: The potential was employed in the TI calculations in the above reference. It can be used as an initial approximation for MD simulations under the Earth’s inner core conditions.

  • LAMMPS pair_style eam/fs (2024--Sun-Y--Fe-Ni--LAMMPS--ipr1)
    See Computed Properties
    Notes: This file was provided by Mikhail Mendelev on February 16, 2024.
    File(s):
  • Citation: Y. Mishin, M.J. Mehl, and D.A. Papaconstantopoulos (2005), "Phase stability in the Fe-Ni system: Investigation by first-principles calculations and atomistic simulations", Acta Materialia, 53(15), 4029-4041. DOI: 10.1016/j.actamat.2005.05.001.
    Abstract: First-principles calculations of the energy of various crystal structures of Fe, Ni and ordered Fe–Ni compounds with different stoichiometries have been performed by the linearized augmented plane wave (LAPW) method in the generalized gradient approximation. The most stable compounds are L12–Ni3Fe, L10–FeNi, C11f–Ni2Fe and C11f–Fe2Ni. The L12-Ni3Fe compound has the largest negative formation energy, which is consistent with the experimental Fe–Ni phase diagram. The L10–FeNi compound has also been observed experimentally in meteorite samples as a metastable phase. It is suggested here that the C11f compounds could also form in Fe–Ni alloys at low temperatures. A new semi-empirical interatomic potential has been developed for the Fe–Ni system by fitting to experimental data and the results of the LAPW calculations. Recognizing the significance of the covalent component of bonding in this system, the potential is based on the embedded-atom method (EAM) but additionally includes a bond-angle dependence. In comparison with the existing modified EAM method, our potential form is simpler, extends interactions to several (3–5) coordination shells and replaces the screening procedure by a smooth cutoff of the potential functions. The potential reproduces a variety of properties of Fe and Ni with a reasonable accuracy. It also reproduces all stability trends across the Fe–Ni system established by the LAPW calculations. The potential can be useful in atomistic simulations of the phases of the Fe–Ni system.

  • ADP tabulated functions (2005--Mishin-Y--Fe-Ni--table--ipr1)
    Notes: These files were provided by Yuri Mishin (George Mason University) and posted on 22 Dec. 2009. Prof. Mishin requested the following note be included: "The equation appearing in the Appendix on page 4040 contains a typing error: the sign before 1/3 in the last line must be negative." He provided the corrected equation for the angular-dependent force contributions in ADP_Forces.jpg or ADP_Forces.pdf.
    File(s):
    Fe F(ρ): F_Fe.plt
    Ni F(ρ): F_Ni.plt
    Fe ρ(r): fFe.plt
    Ni ρ(r): fNi.plt
    Fe φ(r): pFe.plt
    Ni φ(r): pNi.plt
    Fe-Ni φ(r): pFeNi.plt
    Fe u(r): dFe.plt
    Ni u(r): dNi.plt
    Fe-Ni u(r): dFeNi.plt
    Fe w(r): qFe.plt
    Ni w(r): qNi.plt
    Fe-Ni w(r): qFeNi.plt
    ADP_Forces.jpg
    ADP_Forces.pdf

Date Created: October 5, 2010 | Last updated: March 20, 2025