Atomistic-Object Kinetic Monte Carlo simulations of irradiation damage in tungsten
arXiv:1904.04089 · doi:10.1088/1361-651X/ab1a1e
Abstract
We describe the development of a new object kinetic Monte Carlo code where the elementary defect objects are off-lattice atomistic configurations. Atomic-level transitions are used to transform and translate objects, to split objects and to merge them together. This gradually constructs a database of atomic configurations -- a set of relevant defect objects and their possible events generated on-the-fly. Elastic interactions are handled within objects with empirical potentials at short distances, and between spatially distinct objects using the dipole tensor formalism. The model is shown to evolve mobile interstitial clusters in tungsten faster than an equivalent molecular dynamics simulation, even at elevated temperatures. We apply the model to the evolution of complex defects generated using molecular dynamics simulations of primary radiation damage in tungsten. We show that we can evolve defect structures formed in cascade simulations to experimentally observable timescales of seconds while retaining atomistic detail. We conclude that the first few nanoseconds of simulation following cascade initiation would be better performed using molecular dynamics, as this will capture some of the near-temperature-independent evolution of small highly-mobile interstitial clusters. We also conclude that, for the 20keV PKA cascades annealing simulations considered here, internal relaxations of sessile objects difficult to capture using conventional object KMC with idealised object geometries establish the conditions for long timescale evolution.
References in corpus (14)
- Ab initio Random Structure Searching
- High energy collision cascades in tungsten: dislocation loops structure and clustering scaling laws
- Point defect modeling in materials: coupling ab initio and elasticity approaches
- The Kinetic Activation-Relaxation Technique: A Powerful Off-lattice On-the-fly Kinetic Monte Carlo Algorithm
- Kinetic Activation Relaxation Technique
- Monte Carlo with Absorbing Markov Chains: Fast Local Algorithms for Slow Dynamics
- Elastic trapping of dislocation loops in cascades in ion-irradiated tungsten foils
- Displacement cascades and defects annealing in tungsten, Part I: defect database from molecular dynamics simulations
- MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion
- Direct observation of size scaling and elastic interaction between nano-scale defects in collision cascades
- Evolution of small defect clusters in ion-irradiated 3C-SiC: combined cluster dynamics modeling and experimental study
- Displacement cascades and defects annealing in tungsten, Part II: Object kinetic Monte Carlo Simulation of Tungsten Cascade Aging
- Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification
- Self-Learning Kinetic Monte Carlo Simulations of Al Diffusion in Mg
Cited by in corpus (3)
- Parameter-free quantitative simulation of high dose microstructure and hydrogen retention in ion-irradiated tungsten
- Graph Theory Based Approach to Characterize Self Interstitial Defect Morphology
- Neural-Parareal: Dynamically Training Neural Operators as Coarse Solvers for Time-Parallelisation of Fusion MHD Simulations