Atomistic modeling of thermodynamic equilibrium and polymorphism of iron

We develop two new modified embedded-atom method (MEAM) potentials for elemental iron, intended to reproduce the experimental phase stability with respect to both temperature and pressure. These simple interatomic potentials are fitted to a wide variety of material properties of bcc iron in close agreement with experiments. Numerous defect properties of bcc iron and bulk properties of the two close-packed structures calculated with these models are in reasonable agreement with the available first-principles calculations and experiments. Performance at finite temperatures of these models has also been examined using Monte Carlo simulations. We attempt to reproduce the experimental iron polymorphism at finite temperature by means of free energy computations, similar to the procedure previously pursued by Müller et al (2007 J. Phys.: Condens. Matter 19 326220), and re-examine the adequacy of the conclusion drawn in the study by addressing two critical aspects missing in their analysis: (i) the stability of the hcp structure relative to the bcc and fcc structures and (ii) the compatibility between the temperature and pressure dependences of the phase stability. Using two MEAM potentials, we are able to represent all of the observed structural phase transitions in iron. We discuss that the correct reproductions of the phase stability among three crystal structures of iron with respect to both temperature and pressure are incompatible with each other due to the lack of magnetic effects in this class of empirical interatomic potential models. The MEAM potentials developed in this study correctly predict, in the bcc structure, the self-interstitial in the (110) orientation to be the most stable configuration, and the screw dislocation to have a non-degenerate core structure, in contrast to many embedded-atom method potentials for bcc iron in the literature.