Citation: M.R. Fellinger, H. Park, and J.W. Wilkins (2010), "Force-matched embedded-atom method potential for niobium",
Physical Review B,
81(14), 144119. DOI:
10.1103/physrevb.81.144119.
Abstract: Large-scale simulations of plastic deformation and phase transformations in alloys require reliable classical interatomic potentials. We construct an embedded-atom method potential for niobium as the first step in alloy potential development. Optimization of the potential parameters to a well-converged set of density-functional theory (DFT) forces, energies, and stresses produces a reliable and transferable potential for molecular-dynamics simulations. The potential accurately describes properties related to the fitting data and also produces excellent results for quantities outside the fitting range. Structural and elastic properties, defect energetics, and thermal behavior compare well with DFT results and experimental data, e.g., DFT surface energies are reproduced with less than 4% error, generalized stacking-fault energies differ from DFT values by less than 15%, and the melting temperature is within 2% of the experimental value.
IMD option EAM (2010--Fellinger-M-R--Nb--IMD--ipr1)
Notes: These files were provided by Michael Fellinger, Hyoungki Park, and John Wilkins (The Ohio State University) and posted with their permission on 14 July 2010. Details of the fitting procedure and testing can be found in the reference listed above.
File(s):
See Computed PropertiesNotes: These files were provided by Michael Fellinger, Hyoungki Park, and John Wilkins (The Ohio State University) and posted with their permission on 14 July 2010. Mike Fellinger also provided the additional note: "The Nb.eam.alloy file is in the setfl format suitable for the
LAMMPS MD code. This format requires r*phi and rho to be tabulated from r = 0 to r = r_cut. The domain of phi and rho in the published potential is 1.738 ≤ r ≤ 4.75 A. For phi, we extend the cubic polynomial for 1.738 ≤ r ≤ 2.073 A to r = 0. For rho, we linearly extrapolate from r = 1.738 A to r = 0. The potential in the
IMD format is tabulated with 5,001 points for each function. The potential in the LAMMPS setfl format is tabulated with 10,001 points for each function. Comparisons of the two tabulations show very slight differences in some defect energies, probably due to the different numbers of tabulation points."
File(s):
See Computed PropertiesNotes: Listing found at https://openkim.org. This KIM potential is based on the files from 2010--Fellinger-M-R--Nb--LAMMPS--ipr1.
Link(s):