Simple bond-order-type interatomic potential for an intermixed Fe-Cr-C system of metallic and covalent bondings in heat-resistant ferritic steels

Abstract

It is known that M23C6(M = Cr/Fe) behavior in heat-resistant ferritic steels affects the strength of the material at high temperature. The ability to garner direct information regarding the atomic motion using classical molecular dynamics simulations is useful for investigating the M23C6 behavior in heat-resistant ferritic steels. For such classical molecular dynamics calculations, a suitable interatomic potential is needed. To satisfy this requirement, an empirical bond-order-type interatomic potential for Fe-Cr-C systems was developed because the three main elements to simulate the M23C6 behavior in heat-resistant ferritic steels are Fe, Cr, and C. The angular-dependent term, which applies only in non-metallic systems, was determined based on the similarity between a Finnis-Sinclair-type embedded-atom-method interatomic potential and a Tersoff-type bond-order potential. The potential parameters were determined such that the material properties of Fe-Cr-C systems were reproduced. These properties include the energy and lattice constants of 89 crystal structures; the elastic constants of four realistic precipitates; the bulk moduli of B1, B2, and B3 crystals; the surface energies of B1 and B2 crystals; and the defect-formation energies and atomic configurations of 66 Fe-Cr-C complexes. Most of these material properties were found to be reproduced by our proposed empirical bond-order potentials. The formation energies and lattice constants of randomly mixed Fe-Cr alloys calculated using the interatomic potentials were comparable to those obtained through experiments and first-principles calculations. Furthermore, the energies and structures of interfaces between Cr carbide and α-Fe as predicted through first-principles calculations were well reproduced using these interatomic potentials.