Improving atomic displacement and replacement calculations with physically realistic damage models.

Atomic collision processes are fundamental to numerous advanced materials technologies such as electron microscopy, semiconductor processing and nuclear power generation. Extensive experimental and computer simulation studies over the past several decades provide the physical basis for understanding the atomic-scale processes occurring during primary displacement events. The current international standard for quantifying this energetic particle damage, the Norgett-Robinson-Torrens displacements per atom (NRT-dpa) model, has nowadays several well-known limitations. In particular, the number of radiation defects produced in energetic cascades in metals is only ~1/3 the NRT-dpa prediction, while the number of atoms involved in atomic mixing is about a factor of 30 larger than the dpa value. Here we propose two new complementary displacement production estimators (athermal recombination corrected dpa, arc-dpa) and atomic mixing (replacements per atom, rpa) functions that extend the NRT-dpa by providing more physically realistic descriptions of primary defect creation in materials and may become additional standard measures for radiation damage quantification.

Direct Observation of Defect Range and Evolution in Ion-Irradiated Single Crystalline Ni and Ni Binary Alloys.

Energetic ions have been widely used to evaluate the irradiation tolerance of structural materials for nuclear power applications and to modify material properties. It is important to understand the defect production, annihilation and migration mechanisms during and after collision cascades. In this study, single crystalline pure nickel metal and single-phase concentrated solid solution alloys of 50%Ni50%Co (NiCo) and 50%Ni50%Fe (NiFe) without apparent preexisting defect sinks were employed to study defect dynamics under ion irradiation. Both cross-sectional transmission electron microscopy characterization (TEM) and Rutherford backscattering spectrometry channeling (RBS-C) spectra show that the range of radiation-induced defect clusters far exceed the theoretically predicted depth in all materials after high-dose irradiation. Defects in nickel migrate faster than in NiCo and NiFe. Both vacancy-type stacking fault tetrahedra (SFT) and interstitial loops coexist in the same region, which is consistent with molecular dynamics simulations. Kinetic activation relaxation technique (k-ART) simulations for nickel showed that small vacancy clusters, such as di-vacancies and tri-vacancies, created by collision cascades are highly mobile, even at room temperature. The slower migration of defects in the alloy along with more localized energy dissipation of the displacement cascade may lead to enhanced radiation tolerance.

Influence of chemical disorder on energy dissipation and defect evolution in concentrated solid solution alloys.

A grand challenge in materials research is to understand complex electronic correlation and non-equilibrium atomic interactions, and how such intrinsic properties and dynamic processes affect energy transfer and defect evolution in irradiated materials. Here we report that chemical disorder, with an increasing number of principal elements and/or altered concentrations of specific elements, in single-phase concentrated solid solution alloys can lead to substantial reduction in electron mean free path and orders of magnitude decrease in electrical and thermal conductivity. The subsequently slow energy dissipation affects defect dynamics at the early stages, and consequentially may result in less deleterious defects. Suppressed damage accumulation with increasing chemical disorder from pure nickel to binary and to more complex quaternary solid solutions is observed. Understanding and controlling energy dissipation and defect dynamics by altering alloy complexity may pave the way for new design principles of radiation-tolerant structural alloys for energy applications.

Recombination radius of a Frenkel pair and capture radius of a self-interstitial atom by vacancy clusters in bcc Fe.

The recombination radius of a Frenkel pair is a fundamental parameter for the object kinetic Monte Carlo (OKMC) and mean field rate theory (RT) methods that are used to investigate irradiation damage accumulation in irradiated materials. The recombination radius in bcc Fe has been studied both experimentally and numerically, however there is no general consensus about its value. The detailed atomistic processes of recombination also remain uncertain. Values from 1.0a0 to 3.3a0 have been employed as a recombination radius in previous studies using OKMC and RT. The recombination process of a Frenkel pair is investigated at the atomic level using the self-evolved atomistic kinetic Monte Carlo (SEAKMC) method in this paper. SEAKMC calculations reveal that a self-interstitial atom recombines with a vacancy in a spontaneous reaction from several nearby sites following characteristic pathways. The recombination radius of a Frenkel pair is estimated to be 2.26a0 by taking the average of the recombination distances from 80 simulation cases. In addition, we apply these procedures to the capture radius of a self-interstitial atom by a vacancy cluster. The capture radius is found to gradually increase with the size of the vacancy cluster. The fitting curve for the capture radius is obtained as a function of the number of vacancies in the cluster.

Cascade morphology transition in bcc metals.

Energetic atom collisions in solids induce shockwaves with complex morphologies. In this paper, we establish the existence of a morphological transition in such cascades. The order parameter of the morphology is defined as the exponent, b, in the defect production curve as a function of cascade energy (N(F) ~ E(MD)(b)). Response of different bcc metals can be compared in a consistent energy domain when the energy is normalized by the transition energy, µ, between the high- and the low-energy regime. Using Cr, Fe, Mo and W data, an empirical formula of µ as a function of displacement threshold energy, E(d), is presented for bcc metals.

Solving the puzzle of <100> interstitial loop formation in bcc Iron.

The interstitial loop is a unique signature of radiation damage in structural materials for nuclear and other advanced energy systems. Unlike other bcc metals, two types of interstitial loops, 1/2<111> and <100>, are formed in bcc iron and its alloys. However, the mechanism by which <100> interstitial dislocation loops are formed has remained undetermined since they were first observed more than fifty years ago. We describe our atomistic simulations that have provided the first direct observation of <100> loop formation. The process was initially observed using our self-evolving atomistic kinetic Monte Carlo method, and subsequently confirmed using molecular dynamics simulations. Formation of <100> loops involves a distinctly atomistic interaction between two 1/2<111> loops, and does not follow the conventional assumption of dislocation theory, which is Burgers vector conservation between the reactants and the product. The process observed is different from all previously proposed mechanisms. Thus, our observations might provide a direct link between experiments and simulations and new insights into defect formation that may provide a basis to increase the radiation resistance of these strategic materials.

Self-evolving atomistic kinetic Monte Carlo: fundamentals and applications.

The fundamentals of the framework and the details of each component of the self-evolving atomistic kinetic Monte Carlo (SEAKMC) are presented. The strength of this new technique is the ability to simulate dynamic processes with atomistic fidelity that is comparable to molecular dynamics (MD) but on a much longer time scale. The observation that the dimer method preferentially finds the saddle point (SP) with the lowest energy is investigated and found to be true only for defects with high symmetry. In order to estimate the fidelity of dynamics and accuracy of the simulation time, a general criterion is proposed and applied to two representative problems. Applications of SEAKMC for investigating the diffusion of interstitials and vacancies in bcc iron are presented and compared directly with MD simulations, demonstrating that SEAKMC provides results that formerly could be obtained only through MD. The correlation factor for interstitial diffusion in the dumbbell configuration, which is extremely difficult to obtain using MD, is predicted using SEAKMC. The limitations of SEAKMC are also discussed. The paper presents a comprehensive picture of the SEAKMC method in both its unique predictive capabilities and technically important details.