SIMULATING THE PRODUCTION OF FREE DEFECTS IN IRRADIATED METALS

Under cascade-producing irradiation by high energy neutrons or charged particles, only a small fraction of the initially displaced atoms contribute to the population of free defects that are available to migrate throughout the metal and cause microstructural changes. Although, in principle, computer simulations of free defect production could best be done using molecular dynamics, in practice, the wide ranges of time and distance scales involved can be done only by a combination of atomistic models that employ various levels of approximation. An atomic-scale, multi-model approach has been developed that combines molecular dynamics, binary collision models and stochastic annealing simulation. The annealing simulation is utilized in calibrating binary collision simulations to the results of molecular dynamics calculations, as well as to model the subsequent migration of the defects on more macroscopic time and size scales. The annealing simulation and the method of calibrating the multi-model approach are discussed, and the results of simulations of cascades in copper are presented. The temperature dependence of free defect production following simulated annealing of isolated cascades in copper shows a differential in the fractions of free vacancies and interstitial defects escaping from the cascade above stage V. This differential, a consequence of the direct formation of interstitial clusters in cascades and the relative thermal stability of vacancy and interstitial clusters during subsequent annealing, is the basis for the production bias mechanism of void swelling.