The nature of high-energy radiation damage in iron (original) (raw)
Related papers
The nature of high-energy radiation damage in iron: Modeling results
Understanding and predicting a material's performance in response to high-energy radiation damage, as well as designing future materials to be used in intense radiation environments, requires the knowledge of the structure, morphology and amount of radiation-induced structural changes 1-5 . We report the results of molecular dynamics simulations of high-energy radiation damage in iron in the range 0.2-0.5 MeV. We analyze and quantify the nature of collision cascades both at the global and local scale. We find that the structure of high-energy collision cascades becomes increasingly continuous as opposed to showing sub-cascade branching reported previously. At the local length scale, we find large defect clusters and novel small vacancy and interstitial clusters. These features form the basis for physical models aimed at understanding the effects of high energy radiation damage in structural materials.
Electronic effects in high-energy radiation damage in iron
Journal of Physics: Condensed Matter, 2014
Electronic effects are believed to be important in high-energy radiation damage processes where high electronic temperature is expected, yet their effects are not currently understood. Here, we perform molecular dynamics simulations of high-energy collision cascades in α-iron using the coupled two-temperature molecular dynamics (2T-MD) model that incorporates both effects of electronic stopping and electron-ion interaction. We subsequently compare it with the model employing the electronic stopping only, and find several interesting novel insights. The 2T-MD results in both decreased damage production in the thermal spike and faster relaxation of the damage at short times. Notably, the 2T-MD model gives a similar amount of the final damage at longer times, which we interpret to be the result of two competing effects: smaller amount of short-time damage and shorter time available for damage recovery.
Molecular Dynamics Simulations of High Energy Cascades in Iron
MRS Proceedings, 1994
A series of high-energy, up to 20 keV, displacement cascades in iron have been investigated for times up to 200 ps at 100 K using the method of molecular dynamics simulation. Thesimulations were carried out using the MOLDY code and a modified version of the many-bodyinteratomic potential developed by Finnis and Sinclair. The paper focuses on those results obtained at the highest energies, 10 and 20 keV. The results indicate that the fraction of the Frenkel pairs surviving in-cascade recombination remains fairly high in iron and that the fraction of the surviving point defects that cluster is lower than in materials such as copper. In particular, vacancy clustering appears to be inhibited in iron. Some of the interstitial clusters were observed to exhibit an unexpectedly complex, three-dimensional morphology. The observations are discussed in terms of their relevance to microstructural evolution and mechanical property changes in irradiated iron-based alloys.
The primary damage in Fe revisited by Molecular Dynamics and its binary collision approximation
MRS Proceedings, 2000
Molecular Dynamics (MD) is a very powerful tool for studying displacement cascades initiated by the neutrons when they interact with matter and thus evaluate the primary damage. The mean number of point defects created can be obtained with a fair standard error with a reasonable number of cascade simulations (10 to 20 [1]), however other cascades characteristics (spatial distribution, size and amount of defect clusters ...) display a huge variability. Therefore, they may need to be studied using faster methods such as the Binary Collision Approximation (BCA) which is several order of magnitude less time consuming. We have investigated the point defect distributions subsequent to atomic collision cascades by both MD (using EAM potentials for Fe) and its BCA. MD and its BCA lead to comparable point defect predictions. The significant similarities and differences are discussed.
Journal of Nuclear Materials, 1996
The evolution of high-energy displacement cascades in iron has been investigated for times up to 200 ps using the method of molecular dynamics simulation. The simulations were carried out using the MOLDY code and a modified version of the many-body interatomic potential developed by Finnis and Sinclair. Previously reported results have been supplemented by a series of 10 keV simulations at 900 K and 20 keV simulations at 100 K. The results indicate that the fraction of the Frenkel pairs escaping in-cascade recombination is somewhat higher and the fraction of the surviving point defects that cluster is lower in iron than in materials such as copper. In particular, vacancy clustering appears to be inhibited in iron. Many of the larger interstitial clusters were observed to exhibit a complex, three-dimensional morphology. The apparent mobility of the <111> crowdion and clusters of 411> crowdions was very high.
Journal of Nuclear Materials, 2011
ABSTRACT Atomic displacement cascades in solids are complex phenomena, the outcome of which can be statistically characterised by properties such as their spatial extent, morphology and the spatial correlation of defects. Some properties scale in a simple way with parameters such as the cascade energy, others have limited variability with energy, for example point defect cluster size distributions. Taking advantage of the latter invariance, we use object kinetic Monte Carlo simulations to demonstrate that most properties of displacement cascade play no significant role in the evolution of point defect cluster size distributions after long enough time. It is suggested that reliable long-term predictions are possible, when using only the self-interstitial and vacancy cluster size distributions from low energy displacement cascades as building blocks to represent the complete spectrum of cascade energies obtained under neutron irradiation conditions. This is shown on the basis of recursive properties of displacement cascades evidenced for the first time and taking only approximately into account the average volumes in which vacancies and self-interstitial atoms are confined.The model has been successfully used to simulate the evolution of point defect clusters in iron for displacement rates in the range of 10−6 dpa/s and doses of the order of 0.1 dpa. The applicability beyond this range and to more complex materials is discussed.
Philosophical Magazine, 2005
We describe a series of kinetic Monte Carlo simulations of post-cascade radiation damage evolution in -iron that illustrates the part played by elastic interaction between defects. Elastic interactions are included as a bias to the diffusion of mobile point defects and defect clusters. The simulations show that recombination fractions are reduced, and vacancy clustering is enhanced. The sensitivity of these effects to temperature, cascade energy, and geometric description of vacancy clusters is also investigated.
Length-scale Effects in Cascade Damage Production in Iron
MRS Proceedings, 2008
ABSTRACTMolecular dynamics simulations provide an atomistic description of the processes that control primary radiation damage formation in atomic displacement cascades. An extensive database of simulations describing cascade damage production in single crystal iron has been compiled using a modified version of the interatomic potential developed by Finnis and Sinclair. This same potential has been used to investigate primary damage formation in nanocrystalline iron in order to have a direct comparison with the single crystal results. A statistically significant number of simulations were carried out at cascade energies of 10 keV and 20 keV and temperatures of 100 and 600K to make this comparison. The results demonstrate a significant influence of nearby grain boundaries as a sink for mobile defects during the cascade cooling phase. This alters the residual primary damage that survives the cascade event. Compared to single crystal, substantially fewer interstitials survive in the na...
Journal of Nuclear Materials, 2006
Groups of displacement cascades calculated independently with different simulation models and computer codes are compared on a statistical basis. The parameters used for this comparison are the number of Frenkel pairs (FP) produced, the percentages of vacancies and self-interstitial atoms (SIAs) in clusters, the spatial extent and the aspect ratio of the vacancies and the SIAs formed in each cascade. One group of cascades was generated in the binary collision approximation (BCA) and all others by full molecular dynamics (MD). The MD results differ primarily due to the empirical interatomic potentials used and, to some extent, in code strategies. Cascades were generated in simulation boxes at different initial equilibrium temperatures. Only modest differences in the predicted numbers of FP are observed, but the other cascade parameters may differ by more than 100%. The consequences of these differences on long-term cluster growth in a radiation environment are examined by means of object kinetic Monte Carlo (OKMC) simulations. These were repeated with three different parameterizations of SIA and SIA cluster mobility. The differences encompassed low to high mobility, one-and three-dimensional migration of clusters, and complete immobility of large clusters. The OKMC evolution was followed until 0.1 dpa was reached. With the range of OKMC parameters used, cluster populations after 0.1 dpa differ by orders of magnitude. Using the groups of cascades from different sources induced no difference larger than a factor of 2 in the OKMC results. No correlation could be identified between the cascade parameters considered and the number densities of vacancies and SIAs predicted by OKMC to cluster in the long term. However, use of random point defect distributions instead of those obtained for displacement cascades as input for the OKMC modeling led to significantly different results.
The use of atomic level stress to characterize the structure of irradiated iron
Journal of Physics: Conference Series, 2012
The behaviour of irradiated material near a primary knock on atom immediately after impact is of great importance for designing reactor materials. Currently, molecular dynamics simulations with classical force fields provide the foundation for understanding the resulting cascade. However, modern density functional calculations can now treat large enough numbers of atoms that they can provide additional details of the magnetic and electronic nature of irradiated samples. In this paper we calculate from first principles the atomic level stresses for an instantaneous configuration following the initiation of a low energy cascade in iron.