Calculation update! New properties have been added to the website for dislocation monopole core structures, dynamic relaxes of both crystal and liquid phases, and melting temperatures! Currently, the results for these properties predominately focus on EAM-style potentials, but the results will be updated for other potentials as the associated calculations finish. Feel free to give us feedback on the new properties so we can improve their representations as needed.
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: H. Sharifi, and C.D. Wick (2025), "Developing interatomic potentials for complex concentrated alloys of Cu, Ti, Ni, Cr, Co, Al, Fe, and Mn", Computational Materials Science248, 113595. DOI: 10.1016/j.commatsci.2024.113595.
Abstract: Complex concentrated alloys (CCAs) are a new generation of metallic alloys composed of three or more principal elements with physical and mechanical properties that can be tuned by adjusting their compositions. The extensive compositional workspace of CCAs makes it impractical to perform a comprehensive search for a specific material property using experimental measurements. The use of computational methods can rapidly narrow down the search span, improving the efficiency of the design process. We carried out a high-throughput parameterization of modified embedded atom method (MEAM) interatomic potentials for combinations of Cu, Ti, Ni, Cr, Co, Al, Fe, and Mn using a genetic algorithm. Unary systems were parameterized based on DFT calculations and experimental results. MEAM potentials for 28 binary and 56 ternary combinations of the elements were parameterized to DFT results that were carried out with semi-automated frameworks. Specific attention was made to reproduce properties that impact compositional segregation, material strength, and mechanics.
Notes: This is a ternary listing for the 2025--Sharifi-H-Wick-C-D--Fe-Mn-Ni-Ti-Cu-Cr-Co-Al potential. This potential focuses on the structural analysis of alloys including shear strength and elastic constants, dislocation dynamics and their impact on alloy strength, and the analysis of defect effects, such as voids, on material properties. However, the potential was not optimized for temperature-dependent properties and was not fit to density, thermal expansion coefficients, or thermal conductivity data.
Citation: D.R. Tramontina, O.R. Deluigi, R. Pinzón, J. Rojas-Nunez, F.J. Valencia, R.C. Pasianot, S.E. Baltazar, R.I. Gonzalez, and E.M. Bringa (2023), "Probing radiation resistance in simulated metallic core–shell nanoparticles", Computational Materials Science227, 112304. DOI: 10.1016/j.commatsci.2023.112304.
Abstract: We present molecular dynamics (MD) simulations of radiation damage in Fe nanoparticles (NP) and bimetallic FeCu core–shell nanoparticles (CSNP). The CSNP includes a perfect body-centered cubic (bcc) Fe core coated with a face-centered cubic (fcc) Cu shell. Irradiation with Fe Primary Knock-on Atoms (PKA) with energies between 1 and 7 keV leads to point defects, without clustering beyond divacancies and very few slightly larger vacancy clusters, and without interstitial clusters, unlike what happens in bulk at the same PKA energies. The Fe-Cu interface and shell can act as a defect sink, absorbing radiation-induced damage and, therefore, the final number of defects in the Fe core is significantly lower than in the Fe NP. In addition, the Cu shell substantially diminishes the number of sputtered Fe atoms, acting as a barrier for recoil ejection. Structurally, the Cu shell responds to the stress generated by the collision cascade by creating and destroying stacking faults across the shell width, which could also accommodate further irradiation defects. We compare our MD results to Monte Carlo Binary Collision Approximation (BCA) simulations using the SRIM code, for the irradiation of an amorphous 3-layer thin film with a thickness equal to the CSNP diameter. BCA does not include defect recombination, so the number of Frenkel pairs is significantly higher than in MD, as expected. Sputtering yield (Y) is underestimated by BCA, which is also expected since the simulation is for a thin film at normal incidence. We also compare MD defect production to bulk predictions of the analytic Athermal Recombination Corrected Displacements Per Atom (arc-dpa) model. The number of vacancies in the Fe core is only slightly lower than arc-dpa predictions, but the number of interstitials is reduced by about one order of magnitude compared to vacancies, at 5 keV. According to the radiation resistance found for FeCu CSNP in our simulations, this class of nanomaterial could be suitable for developing new radiation-resistant coatings, nanostructured components, and shields for use in extreme environments, for instance, in nuclear energy and astrophysical applications.
Notes: The current interatomic potentials are a modified version of 2009--Bonny-G-Pasianot-R-C-Castin-N-Malerba-L--Fe-Cu-Ni, that include the ZBL correction at short distances, making them suitable for collision cascade simulations. Also, the Ni embedding function is currently modified for densities beyond 1.5 times the equilibrium value, in order to obtain a smooth equation of state behavior. The changes do not impact any of the previously published results.
Citation: G. Bonny, R.C. Pasianot, N. Castin, and L. Malerba (2009), "Ternary Fe-Cu-Ni many-body potential to model reactor pressure vessel steels: First validation by simulated thermal annealing", Philosophical Magazine89(34-36), 3531-3546. DOI: 10.1080/14786430903299824.
Abstract: In recent years, the development of atomistic models dealing with microstructure evolution and subsequent mechanical property change in reactor pressure vessel steels has been recognised as an important complement to experiments. In this framework, a literature study has shown the necessity of many-body interatomic potentials for multi-component alloys. In this paper, we develop a ternary many-body Fe–Cu–Ni potential for this purpose. As a first validation, we used it to perform a simulated thermal annealing study of the Fe–Cu and Fe–Cu–Ni alloys. Good qualitative agreement with experiments is found, although fully quantitative comparison proved impossible, due to limitations in the used simulation techniques. These limitations are also briefly discussed.
Notes: Notes from Giovanni Bonny: The references for the elements and binary potentials used in Fe-Cu-Ni are
Fe: 'potential 2' from M.I. Mendelev, A. Han, D.J. Srolovitz, G.J. Ackland, D.Y. Sun and M. Asta, Phil. Mag. A 83 (2003) 3977.
Cu: 'EAM 1' from Y. Mishin, M.J. Mehl, D.A. Papaconstantopoulos, A.F. Voter, J.D. Kress, Phys. Rev. B 63 (2001) 224106.
Ni: A.F. Voter and S.P. Chen, Mater. Res. Soc. Symp. Proc. 82 (1987) 175.
FeCu: R.C. Pasianot and L. Malerba, J. Nucl. Mater. 360 (2007) 118.
FeNi: G. Bonny, R.C. Pasianot and L. Malerba, Model. Simul. Mater. Sci. Eng. 17 (2009) 025010.
F_Ni.spt was modified for densities past 4.8 because of a discontinuity. Unless for cascade conditions (for which the potential was not stiffened), the properties should stay exactly the same (equilibrium density is around 1).
See Computed Properties Notes: This file was provided by Giovanni Bonny (Nuclear Materials Science Institute of SCK-CEN, Belgium) on 8 Feb. 2010. File(s):
See Computed Properties Notes: Listing found at https://openkim.org. This KIM potential is based on the files from 2009--Bonny-G--Fe-Cu-Ni--LAMMPS--ipr1. Link(s):