Molecular dynamics simulations of the defect evolution in tungsten on successive collision cascades

Utkarsh Bhardwaj,Manoj Warrier
2024-08-26
Abstract:Molecular dynamics (MD) simulations of successive collision cascades within the same simulation domain were performed using two different inter-atomic potentials (IAP) in tungsten, one EAM based and the other a `quantum accurate' machine learning potential, SNAP. The micro-structural changes are analyzed as a function of displacements per atom (dpa) for primary knock-on atom (PKA) energies of 20 keV and 50 keV, reaching up-to irradiation dose of 0.1 and 0.2 dpa, respectively. Five sample simulations are carried out for each case for observing stochastic differences in the evolution of damage. A detailed defect analysis is carried out to observe changes in different parameters such as the number of surface defects, defect density, defect morphology and size distribution etc., as a function of dpa. We explore the properties that are sensitive to the IAP used and those that are sensitive to the PKA energy and note their similarities with experimental results at various dpa values. The SNAP potential shows better agreement with the experiments for swelling and number of surface defects. However, it also predicts presence of high number of small sessile defects which may have definite affect on the processes of microstructural evolution and material properties.
Materials Science
What problem does this paper attempt to address?
The paper attempts to address the issue of using two different interatomic potentials (IAP) for molecular dynamics (MD) simulations of sequential collision cascades (SCC) in tungsten at different primary knock-on atom (PKA) energies, and to analyze the impact of these collisions on the microstructural evolution of the material. Specifically: 1. **Microstructural Evolution**: By conducting molecular dynamics simulations of sequential collision cascades within the same simulation domain, the study investigates the different effects of two interatomic potentials (one based on the Embedded Atom Method (EAM) and the other based on a machine learning quantum-accurate potential (SNAP)) on defect evolution in tungsten. 2. **Defect Characteristics Analysis**: The study analyzes the changes in parameters such as the number of defects, defect density, defect morphology, and size distribution at different displacements per atom (dpa), particularly at PKA energies of 20 keV and 50 keV, reaching irradiation doses of 0.1 and 0.2 dpa, respectively. 3. **Comparison with Experimental Results**: The results obtained under different IAP and PKA energy conditions are compared with experimental data, especially regarding the number of swelling and surface defects. 4. **Microstructural Stability**: The stability of defect clusters and their transport and interaction at different dpa values are explored to better understand the microstructural evolution process due to irradiation. 5. **Importance of Surface Defects**: The study shows that although there are significant differences in the defect morphology and size predicted by the two potentials, they exhibit similar trends in surface point defect density and the number of TEM-visible defect loops, with the SNAP potential showing a higher degree of agreement with experimental results. Through these studies, the paper aims to provide a more accurate theoretical basis for the microstructural evolution of tungsten materials under high-energy particle irradiation and to validate the performance of different interatomic potentials in simulations.