Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: G. Bonny, R.C. Pasianot, D. Terentyev, and L. Malerba (2011), "Iron chromium potential to model high-chromium ferritic alloys", Philosophical Magazine, 91(12), 1724-1746. DOI: 10.1080/14786435.2010.545780.
Abstract: We present an Fe–Cr interatomic potential to model high-Cr ferritic alloys. The potential is fitted to thermodynamic and point-defect properties obtained from density functional theory (DFT) calculations and experiments. The developed potential is also benchmarked against other potentials available in literature. It shows particularly good agreement with the DFT obtained mixing enthalpy of the random alloy, the formation energy of intermetallics and experimental excess vibrational entropy and phase diagram. In addition, DFT calculated point-defect properties, both interstitial and substitutional, are well reproduced, as is the screw dislocation core structure. As a first validation of the potential, we study the precipitation hardening of Fe–Cr alloys via static simulations of the interaction between Cr precipitates and screw dislocations. It is concluded that the description of the dislocation core modification near a precipitate might have a significant influence on the interaction mechanisms observed in dynamic simulations.
EAM tabulated functions (2011--Bonny-G--Fe-Cr--table--ipr1)
Notes: These files were sent by Dr. Giovanni Bonny (Nuclear Materials Science Institute of SCK-CEN, Belgium) on 2 November 2017 and posted with his permission. File(s):
See Computed Properties Notes: These files were sent by Dr. Giovanni Bonny (Nuclear Materials Science Institute of SCK-CEN, Belgium) on 2 November 2017 and posted with his permission. Giovanni Bonny also included Caution.pdf file, which explains why a large number of grid points for the s-embedding function are necessary. Giovanni Bonny noted that this warning is in fact valid for all known two-band model (2BM) potentials. Giovanni Bonny thanks Junlei Zhao (Department of Physics, University of Helsinki, Finland) for help in preparation of the LAMMPS files. Update March 15, 2020: This version was identified to not be compatible with LAMMPS versions after 7 Aug 2019 due to more rigorous format checks. File(s): superseded
See Computed Properties Notes: This is a modification to the above version posted by Lucas Hale on March 15, 2020. Missing pair function tables of all zeros were added to the FeCr_s.eam.fs file to make the files compatible with LAMMPS versions after 7 Aug 2019. Update May 26 2021: This version is not compatible for LAMMPS versions starting with 29 Oct 2020 due to Infinify and NaN values no longer allowed. File(s): superseded
See Computed Properties Notes: This is a modification to the above version posted by Lucas Hale on May 26, 2021. To make the files compatible with LAMMPS versions after 29 Oct 2020, the Infinity and NaN values associated with the Cr elemental tables at r=0 are replaced with 1e+8 and 0.0 values, respectively. File(s):