× Updated! Potentials that share interactions are now listed as related models.
 
Citation: M.I. Mendelev (2019), "to be published".

Notes: This potential was designed to simulate stainless steels. All pure components potentials here are original and are not the same as previously published potentials from Mendelev.

See Computed Properties
Notes: This file was provided by Mikhail Mendelev (Ames Laboratory) on 8 October 2019. Update 19 July 2021: The contact email in the file's header has been changed.
File(s):
See Computed Properties
Notes: Listing found at https://openkim.org. This KIM potential is based on the files from 2019--Mendelev-M-I--Fe-Ni-Cr--LAMMPS--ipr1.
Link(s):
Citation: X.W. Zhou, M.E. Foster, and R.B. Sills (2018), "An Fe-Ni-Cr embedded atom method potential for austenitic and ferritic systems", Journal of Computational Chemistry, 39(29), 2420-2431. DOI: 10.1002/jcc.25573.
Abstract: Fe‐Ni‐Cr stainless‐steels are important structural materials because of their superior strength and corrosion resistance. Atomistic studies of mechanical properties of stainless‐steels, however, have been limited by the lack of high‐fidelity interatomic potentials. Here using density functional theory as a guide, we have developed a new Fe‐Ni‐Cr embedded atom method potential. We demonstrate that our potential enables stable molecular dynamics simulations of stainless‐steel alloys at high temperatures, accurately reproduces the stacking fault energy-known to strongly influence the mode of plastic deformation (e.g., twinning vs. dislocation glide vs. cross‐slip)-of these alloys over a range of compositions, and gives reasonable elastic constants, energies, and volumes for various compositions. The latter are pertinent for determining short‐range order and solute strengthening effects. Our results suggest that our potential is suitable for studying mechanical properties of austenitic and ferritic stainless‐steels which have vast implementation in the scientific and industrial communities. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.

LAMMPS pair_style eam/alloy (2018--Zhou-X-W--Fe-Ni-Cr--LAMMPS--ipr1)
See Computed Properties
Notes: This file was provided by Xiaowang Zhou (Sandia National Laboratories) on 1 January 2019 and posted with his permission. The function tabulations are identical to 2018--Zhou-X-W--Fe-Ni-Cr--LAMMPS--ipr2 below, only the file format is different.
File(s):
See Computed Properties
Notes: This file was provided by Xiaowang Zhou (Sandia National Laboratories) on 1 January 2019 and posted with his permission. The function tabulations are identical to 2018--Zhou-X-W--Fe-Ni-Cr--LAMMPS--ipr1 above, only the file format is different.
File(s):
Citation: L.K. Béland, A. Tamm, S. Mu, G.D. Samolyuk, Y.N. Osetsky, A. Aabloo, M. Klintenberg, A. Caro, and R.E. Stoller (2017), "Accurate classical short-range forces for the study of collision cascades in Fe–Ni–Cr", Computer Physics Communications, 219, 11-19. DOI: 10.1016/j.cpc.2017.05.001.
Abstract: The predictive power of a classical molecular dynamics simulation is largely determined by the physical validity of its underlying empirical potential. In the case of high-energy collision cascades, it was recently shown that correctly modeling interactions at short distances is necessary to accurately predict primary damage production. An ab initio based framework is introduced for modifying an existing embedded-atom method FeNiCr potential to handle these short-range interactions. Density functional theory is used to calculate the energetics of two atoms approaching each other, embedded in the alloy, and to calculate the equation of state of the alloy as it is compressed. The pairwise terms and the embedding terms of the potential are modified in accordance with the ab initio results. Using this reparametrized potential, collision cascades are performed in Ni50Fe50, Ni80Cr20 and Ni33Fe33Cr33. The simulations reveal that alloying Ni and NiCr to Fe reduces primary damage production, in agreement with some previous calculations. Alloying Ni and NiFe to Cr does not reduce primary damage production, in contradiction with previous calculations.

Notes: Prof. Béland notes that "The potential takes the 2011--Bonny-G-Terentyev-D-Pasianot-R-C-et-al--Fe-Ni-Cr potential and re-parameterizes the short-distance interactions based on DFT calculations, as explained in the paper and https://doi.org/10.1021/acs.jctc.5b01194. We recommend using this potential for simulating collision cascades."

LAMMPS pair_style eam/alloy (2017--Beland-L-K--Fe-Ni-Cr--LAMMPS--ipr1)
See Computed Properties
Notes: This file was provided by Laurent Béland on 7 Nov 2019 and posted with his permission.
File(s):
Citation: G. Bonny, A. Bakaev, P. Olsson, C. Domain, E.E. Zhurkin, and M. Posselt (2017), "Interatomic potential to study the formation of NiCr clusters in high Cr ferritic steels", Journal of Nuclear Materials, 484, 42-50. DOI: 10.1016/j.jnucmat.2016.11.017.
Abstract: Under irradiation NiSiPCr clusters are formed in high-Cr ferritic martensitic steels as well as in FeCr model alloys. In the literature little is known about the origin and contribution to the hardening of these clusters. In this work we performed density functional theory (DFT) calculations to study the stability of small substitutional NiCr-vacancy clusters and interstitial configurations in bcc Fe. Based on DFT data and experimental considerations a ternary potential for the ferritic FeNiCr system was developed. The potential was applied to study the thermodynamic stability of NiCr clusters by means of Metropolis Monte Carlo (MMC) simulations. The results of our simulations show that Cr and Ni precipitate as separate fractions and suggest only a limited synergetic effect between Ni and Cr. Therefore our results suggest that the NiCrSiP clusters observed in experiments must be the result of other mechanisms than the synergy of Cr and Ni at thermal equilibrium.

LAMMPS pair_style hybrid/overlay eam/alloy eam/fs (2017--Bonny-G--Fe-Cr-Ni--LAMMPS--ipr1)
See Computed Properties
Notes: These files were provided by Giovanni Bonny on October 31, 2023. The files were modified to be compatible with newer versions of LAMMPS by replacing inf and nan values at r=0 in the d file, and adding values of 0.0 to the s file for the missing and unused rphi tables.
File(s):
Citation: C. Wu, B.-J. Lee, and X. Su (2017), "Modified embedded-atom interatomic potential for Fe-Ni, Cr-Ni and Fe-Cr-Ni systems", Calphad, 57, 98-106. DOI: 10.1016/j.calphad.2017.03.007.
Abstract: A semi-empirical interatomic potential formalism, the second-nearest-neighbor modified embedded-atom method (2NN MEAM), has been applied to obtaining interatomic potentials for the Fe-Ni, Cr-Ni and Fe-Cr-Ni systems using previously developed MEAM potentials of Fe and Ni and a newly revised potential of Cr. The potential parameters were determined by fitting the experimental data on the enthalpy of formation or mixing, lattice parameter and elastic constant. The present potentials generally reproduced the fundamental physical properties of the Fe-Ni and Cr-Ni alloys. The enthalpy of formation or mixing of the disordered phase at finite temperature and the enthalpy of mixing of the liquid phase are reasonable in agreements with experiment data and CALPHAD calculations. The potentials can be combined with already-developed MEAM potentials to describe Fe-Cr-Ni-based multicomponent alloys. Moreover, the average diffusivities in the unary, some binary and ternary alloys were simulated based on present potential. Good agreement is obtained in comparison with experimental data.

See Computed Properties
Notes: These potential files were obtained from http://cmse.postech.ac.kr/home_2nnmeam, accessed Nov 9, 2020.
File(s):
Citation: G. Bonny, N. Castin, and D. Terentyev (2013), "Interatomic potential for studying ageing under irradiation in stainless steels: the FeNiCr model alloy", Modelling and Simulation in Materials Science and Engineering, 21(8), 085004. DOI: 10.1088/0965-0393/21/8/085004.
Abstract: The degradation of austenitic stainless steels in a radiation environment is a known problem for the in-core components of nuclear light water reactors. For a better understanding of the prevailing mechanisms responsible for the materials' degradation, large-scale atomistic simulations are desirable. In this framework and as a follow-up on Bonny et al (2011 Modelling Simul. Mater. Sci. Eng. 19 085008), we developed an embedded atom method type interatomic potential for the ternary FeNiCr system to model the production and evolution of radiation defects. Special attention has been drawn to the Fe10Ni20Cr alloy, whose properties were ensured to be close to those of 316L austenitic stainless steels. The potential is extensively benchmarked against density functional theory calculations and the potential developed in our earlier work. As a first validation, the potential is used in AKMC simulations to simulate thermal annealing experiments in order to determine the self-diffusion coefficients of the components in FeNiCr alloys around the Fe10Ni20Cr composition. The results from these simulations are consistent with experiments, i.e., DCr > DNi > DFe.

Notes: Notes from Giovanni Bonny: "The present potential was developed to model POINT DEFECTS near the Fe-10Ni-20Cr composition.

LAMMPS pair_style eam/alloy (2013--Bonny-G--Fe-Ni-Cr--LAMMPS--ipr1)
See Computed Properties
Notes: This file was provided by Giovanni Bonny (Nuclear Materials Science Institute of SCK-CEN, Belgium) on 13 Jan. 2014.
File(s):
EAM tabulated functions (2013--Bonny-G--Fe-Ni-Cr--table--ipr1)
Notes: These files were provided by Giovanni Bonny on 13 Jan. 2014.
File(s):
See Computed Properties
Notes: Listing found at https://openkim.org. This KIM potential is based on the files from 2013--Bonny-G--Fe-Ni-Cr--LAMMPS--ipr1.
Link(s):
Citation: G. Bonny, D. Terentyev, R.C. Pasianot, S. Poncé, and A. Bakaev (2011), "Interatomic potential to study plasticity in stainless steels: the FeNiCr model alloy", Modelling and Simulation in Materials Science and Engineering, 19(8), 085008. DOI: 10.1088/0965-0393/19/8/085008.
Abstract: Austenitic stainless steels are commonly used materials for in-core components of nuclear light water reactors. In service, such components are exposed to harsh conditions: intense neutron irradiation, mechanical and thermal stresses, and aggressive corrosion environment which all contribute to the components' degradation. For a better understanding of the prevailing mechanisms responsible for the materials degradation, large-scale atomistic simulations are desirable. In this framework we developed an embedded atom method type interatomic potential for the ternary FeNiCr system to model movement of dislocations and their interaction with radiation defects. Special attention has been drawn to the Fe-10Ni-20Cr alloy, whose properties were ensured to be close to those of 316L austenitic stainless steel. In particular, the stacking fault energy and elastic constants are well reproduced. The fcc phase for the Fe–10Ni-20Cr random alloy was proven to be stable in the temperature range 0–900 K and under shear strain up to 5%. For the same alloy the stable glide of screw dislocations and stability of Frank loops was confirmed.

Notes: Notes from Giovanni Bonny: "The present potential was developed to model dislocations around the Fe-10Ni-20Cr composition."

LAMMPS pair_style eam/alloy (2011--Bonny-G--Fe-Ni-Cr--LAMMPS--ipr1)
See Computed Properties
Notes: This file was provided by Giovanni Bonny (Nuclear Materials Science Institute of SCK-CEN, Belgium) on 2 Sept. 2013. 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


EAM tabulated functions (2011--Bonny-G--Fe-Ni-Cr--table--ipr1)
Notes: These files were provided by Giovanni Bonny on 2 Sept. 2013.
File(s):
Fe F(ρ): F_Fe.spt
Ni F(ρ): F_Ni.spt
Cr F(ρ): F_Cr.spt
Fe ρ(r): rhoFe.spt
Ni ρ(r): rhoNi.spt
Cr ρ(r): rhoCr.spt
Fe φ(r): pFeFe.spt
Ni φ(r): pNiNi.spt
Cr φ(r): pCrCr.spt
Fe-Ni φ(r): pFeNi.spt
Fe-Cr φ(r): pFeCr.spt
Ni-Cr φ(r): pNiCr.spt

LAMMPS pair_style eam/alloy (2011--Bonny-G--Fe-Ni-Cr--LAMMPS--ipr2)
See Computed Properties
Notes: This is a modification to the previous LAMMPS version and was posted by Lucas Hale on May 26, 2021. To make the file compatible with LAMMPS versions after 29 Oct 2020, INF values at r=0 for the elemental r*phi tables were replaced by values computed using the parameters listed in the paper.
File(s):
Date Created: October 5, 2010 | Last updated: November 14, 2023