× Updated! Potentials that share interactions are now listed as related models.
 
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 Science, 227, 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.

See Computed Properties
Notes: This file was provided by Roberto Pasianot on 8 June 2023.
File(s):
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 Magazine, 89(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).

LAMMPS pair_style eam/alloy (2009--Bonny-G--Fe-Cu-Ni--LAMMPS--ipr1)
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):
EAM tabulated functions (2009--Bonny-G--Fe-Cu-Ni--table--ipr1)
Notes: These files were provided by Giovanni Bonny on 8 Feb. 2010.
File(s):
Fe F(ρ): F_Fe.spt
Ni F(ρ): F_Ni.spt
Cu F(ρ): F_Cu.spt
Fe ρ(r): rhoFe.spt
Ni ρ(r): rhoNi.spt
Cu ρ(r): rhoCu.spt
Fe φ(r): pFeFe.spt
Ni φ(r): pNiNi.spt
Cu φ(r): pCuCu.spt
Fe-Ni φ(r): pFeNi.spt
Fe-Cu φ(r): pFeCu.spt
Cu-Ni φ(r): pCuNi.spt

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):
Date Created: October 5, 2010 | Last updated: October 31, 2023