Review and construction of interatomic potentials for molecular dynamics studies of hydrogen embrittlement in Fe─C based steels

Reducing hydrogen embrittlement in the low‐cost Fe─C based steels have the potential to significantly impact the development of hydrogen energy technologies. Molecular dynamics studies of hydrogen interactions with Fe─C steels provide fundamental information about the behavior of hydrogen at microstructural length scales, although such studies have not been performed due to the lack of an Fe─C─H ternary interatomic potential. In this work, the literature on interatomic potentials related to the Fe─C─H systems are reviewed with the aim of constructing an Fe─C─H potential from the published binary potentials. We found that Fe─C, Fe─H, and C─H bond order potentials exist and can be combined to construct an Fe─C─H ternary potential. Therefore, we constructed two such Fe─C─H potentials and demonstrate that these ternary potentials can reasonably capture hydrogen effects on deformation characteristics and deformation mechanisms for a variety of microstructural variations of the Fe─C steels, including martensite that results from γ to α phase transformation, and pearlite that results from the eutectic formation of the Fe 3 C cementite compound.