Since 18 of December 2019 conferences.iaea.org uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.

26–28 May 2025
IAEA
Europe/Vienna timezone
7th International Workshop on Models and Data for Plasma-Material Interactions in Fusion Devices (MoD-PMI)

An impasse in molecular dynamics simulations of primary damage in W: which interatomic potential yields the correct defect morphologies?

Not scheduled
20m
Press Room (IAEA)

Press Room

IAEA

Vienna International Centre Wagramer Strasse 5 1220 Vienna Austria
poster

Speaker

Manoj Warrier (Bhabha Atomic Research Center)

Description

Five widely used interatomic potentials (IAP) for tungsten were used to simulate collision cascades in crystal tungsten. Three of the IAP were embedded atom model (EAM) based [1,2,3] whilst the other two were machine learning (ML) based [4,5]. The molecular dynamics (MD) simulations were carried out for primary knock-on atoms (PKA) having energies 5, 10, 20, 50, 75, 100 and 150 keV. The PKA were directed along 100 random directions for sufficient statistics at each PKA energy, amounting to a total of 3500 collision cascade simulations. Several terabytes of position data of the atoms were then analyzed using CSaransh [6] to obtain the number of defects, the number of defects in clusters, the defect size distributions, the various classes of defects based on their morphology and orientation using a graph theory based method [7]. It is seen that whilst the number of defects predicted by the IAP compares well with each other and with the arcDPA predictions within statistical errors, other aspects of the defects like the number of defects in clusters, the defect size distributions, the number of defects in various classes, etc., show substantial variation for the various IAP. It is seen that there is no trend for these parameters with traditional properties of the IAP like stiffness/range or their ratio S/R. The main difference is the formation of C15 like rings in two of the IAP, whilst in the other three IAP <111> defect clusters dominate [8]. An experimental and a simulation approach to validate the IAP based on the morphology of the defect it produces will be discussed.

REFERENCES:
1. N. Juslin, B. Wirth, Interatomic potentials for simulation of he bubble formation in w, Journal of Nuclear Materials 432 (1) (2013) 61 – 66.
2. A. Sand, J. Dequeker, C. Becquart, C. Domain, K. Nordlund, Non-equilibrium properties of interatomic potentials in cascade simulations in tungsten, Journal of Nuclear Materials 470 (2016) 119–127.
3. C. Björkas, K. Nordlund, S. Dudarev, Modelling radiation effects using the ab-initio based tungsten and vanadium potentials, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 267 (18) (2009) 3204 – 3208
4. M. A. Wood, A. P. Thompson, Quantum-accurate molecular dynamics potential for tungsten (2017) arXiv:1702.07042.
5. J. Byggmästar, K. Nordlund, F. Djurabekova, Modeling refractory high-entropy alloys with efficient machine-learned interatomic potentials: Defects and segregation, Physical Review B 104 (10).
6. U. Bhardwaj https://github.com/haptork/csaransh
7. U. Bhradwaj, A. Sand, M. Warrier, Graph theory based approach to characterize self interstitial defect morphology, Camputational Materials Science, 195 (2021) 110474.
8. M. Warrier, U. Bhradwaj, Defect cluster morphology in W from collision cascades: Results comparing five interatomic potentials arxiv:2042.00359 (https://arxiv.org/abs/2402.00359

Primary author

Manoj Warrier (Bhabha Atomic Research Center)

Co-author

Utkarsh Bhardwaj (Bhabha Atomic Research Centre)

Presentation materials

There are no materials yet.