Displacement cascades and defect annealing in tungsten, Part II: Object kinetic Monte Carlo simulation of tungsten cascade aging
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Two classes for the molecular dynamics program LAMMPS – one compute style and one dump style – are presented that are designed to identify, count, and output point defects in cascade damage and related molecular dynamics simulations. The calculations are done in parallel across multiple MPI processes, and can be done at the same time as the original simulation. This drastically reduces storage requirements by eliminating the need to post-process every atom in the system. Non-cubic lattices, free surfaces, and large voids can be eliminated from the output by suitable choices of the reference lattice. The classes are derived from LAMMPS’s Dump and Compute classes, and pose no additional overhead to LAMMPS if they are not used.
A detailed study has been carried out on recrystallised Zr-1.0Nb alloys corroded and irradiated under different conditions, including ex-autoclave and ex-reactor samples. After 540 days in reactor and damage around 5 dpa, the neutron irradiated sample shows no serious evidence for radiation enhanced corrosion and is still in pre-transition stage. The good corrosion resistance of the neutron irradiated Zr-1.0Nb can be related to the higher volume fraction of tetragonal phase and fewer interconnected nano porosity/cracks in the oxide. This indicates less tetragonal to monoclinic transition, leads to more protective oxide in the neutron irradiated sample, containing little evidence for short circuit paths for the penetration of oxygen or water towards the metal-oxide interface. These observations on a sample with a slow overall oxidation rate are consistent with the hypothesis that interconnected porosity can lead to early transitions and rapid oxidation. Tetragonal oxide can be either stabilised by irradiation, or stabilised by local release of impurity species from SPPs such as dissolution of Fe from Zr-Nb-Fe precipitates or radiation introduced precipitates (RIPs) which is likely to be small β-Nb clusters. The oxide consists of well-aligned columnar-equiaxed microstructure in the autoclave sample while a more complex oxide grain structure was observed in the neutron-irradiated sample. As oxide continues to grow, there are more < a > and < c > loops, dissolved Fe and RIPs in the metal matrix, however, the corrosion rate is low enough for the tetragonal oxide to stabilise and suboxide + Zr(Osat) phases exist for protectiveness, so there is no enhanced corrosion after radiation. In situ ion irradiation in the TEM revealed no visible defect clusters or voids in the oxide, suggesting that radiation damage to the metal matrix rather than oxide may have a stronger effect on corrosion mechanisms after neutron irradiation, however, cascade damages are not visible in this case. Neutron irradiation also seems to have little effect on promoting fast oxidation or dissolution of β-Nb precipitates into the surrounding oxide or metal during irradiation. Furthermore, these results are discussed in the light of the current mechanisms for corrosion of nuclear fuel cladding alloys.
With significant improvement in High Temperature Superconductors (HTS), several projects are adopting HTS technology for fusion power systems. Compact HTS tokamaks offer potential advantages including lower plant costs, enhanced plasma control, and ultimately lower cost of electricity. However, as compact reactors have a reduced radial build to accommodate shielding, HTS degradation due to radiation damage or heating is a significant and potentially design limiting issue. Shielding must mitigate threats to the superconducting coils: neutron cascade damage, heat deposition and potentially organic insulator damage due x-rays. Unfortunately, there are currently no hi-performance shielding materials to enable the potential performance enhancement offered by HTS. In this work, we present a manufacturing method to fabricate a new class of composite shields that are high performance, high operating temperature, and simultaneously neutron absorbing and neutron moderating. The composite design consists of an entrained metal-hydride phase within a radiation stable MgO ceramic host matrix. We discuss the fabrication, characterization, and thermophysical performance data for a series of down-selected composite materials inspired by future fusion core designs and their operational performance metrics. To our knowledge these materials represent the first ceramic composite shield materials containing significant metal hydrides.
Enhanced radiation tolerance of nanostructured metals is attributed to the high density of interfaces that can absorb radiation-induced defects. Here, cavity evolution mechanisms during cascade damage, helium implantation, and annealing of nanocrystalline nickel are characterized via in situ transmission electron microscopy (TEM). Films subjected to self-ion irradiation followed by helium implantation developed evenly distributed cavity structures, whereas films exposed in the reversed order developed cavities preferentially distributed along grain boundaries. Post-irradiation annealing and orientation mapping demonstrated uniform cavity growth in the nanocrystalline structure, and cavities spanning multiple grains. Furthermore, these mechanisms suggest limited ability to reduce swelling, despite the stability of the nanostructure.
This project was to use light ion beam induced charge (IBIC) to detect damage cascades generated by a single heavy ion, and thereby reveal details of the shape of the cascade and the physics of recombination of carriers that interact with the cluster. Further IBIC measurements using the hardware and software of this project will improve the accuracy of theoretical models used to predict electrical degradation in devices exposed to radiation environments. In addition, future use of light ion IBIC detection of single ion-induced damage could be used to locate single ion implantation sites in quantum computing applications. This project used Sandia's Pelletron and nanoImplanter (nI) to produce heavy ion-induced collision cascades in p-n diodes, simulating cascades made by primary knock-on atoms recoiled by neutrons. Si and Li beams from the nI were used to perform highly focused scans generating IBIC signal maps where regions of lower charge collection efficiency were observed without incurring further damage. The very first use of ion channeled beams for IBIC was explored to maximize ionization, improve contrast and provide very straight line trajectories to improve lateral resolution.
The effect of stress on irradiation responses of highly oriented pyrolytic graphite (HOPG) was studied by combing molecular dynamics (MD) simulation, proton irradiation, and Raman characterization. MD simulations of carbon knock-on at energies < 60 eV were used to obtain average threshold displacement energies (E¯d) as a function of strain ranging from 0 to 10%. Simulations at a higher irradiation energy of 2–5 keV were used to study the effect of strain on damage cascade evolution. With increasing tensile strain, E¯d was reduced from 35 eV at 0% strain to 31 eV at 10% strain. The strain-reduced E¯d led to a higher damage peak and more surviving defects (up to 1 ps). Furthermore, high strains induced local cleavage around the cavities, as one additional mechanism of damage enhancement. Experimentally, HOPG film was folded, and the folded region with the maximum tensile stress was irradiated by a 2 MeV proton beam. Raman characterization showed significantly enhanced D to G modes in comparison to the stress-free irradiation. Based on the strain dependence of E¯d and the Kinchin–Pease model, a formula for displacement estimation under different tensile strains is proposed. The stress effects need to be considered in graphite applications in a reactor’s harsh environment where both neutron damage and stress are present.
Modern semiconductor devices, such as gate-all-around nanosheet field-effect transistors (GAA NS FETs), are smaller than displacement damage cascades from fission neutrons. In this regime, device failure may occur through low-probability single events, rather than by parametric degradation previously seen in larger devices. Here, we present a statistical model that predicts the probability of a damage event in a small device and the probability distribution of the magnitude, i.e., number of displacements within the device, from each event. The model is developed first for neutron irradiation and then for energetic ion irradiation. The model is consistent with results from recent experiments in which lithium-ion irradiation produced stepwise increases in subthreshold current in GAA NS FETs.
The project is aimed at understanding the fundamentals of radiation response of low dimensional carbon systems and irradiation-induced mechanical property changes, with a focus on the unique phenomena caused by their geometry, boundary, and quantum size effects. In comparison with their bulk counterparts (graphite), radiation responses of low-dimensional carbon systems are substantially different. The differences exist at almost every stage of defect development, namely, displacement creation, damage cascade, thermal spike formation defect recombination, defect clustering, and structural reconstruction. Many traditional concepts in ion-solid interaction theory do not apply at the nanoscale or require substantial modification. Through this project, we systematically studied various aspects of radiation damage development in graphene, carbon nanotubes (CNTs), and graphite.
Ion implantation has emerged as a powerful tool for manipulating and investigating oxide quantum materials, where strong electron correlations give rise to exotic phases and functionalities. This Review examines how controlled ion bombardment enables targeted modulation of structural, electronic, and magnetic properties in complex oxides. Recent advances in ion beam techniques is analyzed that provide unprecedented control over defect profiles, from point defects to extended damage cascades. In application, ion irradiation has demonstrated the ability to tune quantum phenomena including metal‐insulator transitions, magnetism, and superconductivity through selective doping, defect creation, and lattice strain. The ability to write quantum states with nanoscale precision using focused ion beams is opening new paths for quantum device fabrication. Challenges in understanding radiation damage mechanisms in strongly correlated systems is discussed, and identify opportunities for designing quantum functionality through defect‐engineering approaches. The technique's demonstrated capability to write quantum states with nanoscale precision, combined with its compatibility with existing semiconductor processing infrastructure, suggests a clear pathway for translating fundamental discoveries in quantum materials into practical technological applications.
By means of molecular dynamics (MD) simulations, we have calculated Kaptiza resistance of UO2 with or without radiation damage. For coincident site lattice boundaries of different configurations, the boundary thermal resistance of unirradiated UO2 can be well described by a parameter-reduced formula by using boundary energies as variables. We extended the study to defect-loaded UO2 by introducing damage cascades in close vicinity to the boundaries. Following cascade annealing and defect migrations towards grain boundaries, the boundary energy increases and so does Kaptiza resistance. The correlations between these two still follow the same formula extracted from the unirradiated UO2. The finding will benefit multi-scale modeling of UO2 thermal properties under extreme radiation conditions by combining effects from boundary configurations and damage levels.
Abstract While radiation is known to degrade AlGaN/GaN high-electron-mobility transistors (HEMTs), the question remains on the extent of damage governed by the presence of an electrical field in the device. In this study, we induced displacement damage in HEMTs in both ON and OFF states by irradiating with 2.8 MeV Au 4+ ion to fluence levels ranging from 1.72 × 10 10 to 3.745 × 10 13 ions cm −2 , or 0.001–2 displacement per atom (dpa). Electrical measurement is done in situ , and high-resolution transmission electron microscopy (HRTEM), energy dispersive x-ray (EDX), geometrical phase analysis (GPA), and micro-Raman are performed on the highest fluence of Au 4+ irradiated devices. The selected heavy ion irradiation causes cascade damage in the passivation, AlGaN, and GaN layers and at all associated interfaces. After just 0.1 dpa, the current density in the ON-mode device deteriorates by two orders of magnitude, whereas the OFF-mode device totally ceases to operate. Moreover, six orders of magnitude increase in leakage current and loss of gate control over the 2-dimensional electron gas channel are observed. GPA and Raman analysis reveal strain relaxation after a 2 dpa damage level in devices. Significant defects and intermixing of atoms near AlGaN/GaN interfaces and GaN layer are found from HRTEM and EDX analyses, which can substantially alter device characteristics and result in complete failure.
Atomistic modeling of radiation damage through displacement cascades is deceptively non-trivial. Because of the high energy and stochastic nature of atomic collisions, individual primary knock-on atom (PKA) cascade simulations are computationally expensive and ill-suited for length and dose upscaling. Here, we propose a reduced-order atomistic cascade model capable of predicting and replicating radiation events in metals across a wide range of recoil energies. Our methodology approximates cascade and displacement damage production by modeling the cascade as a core-shell atomic structure composed of two damage production estimators, namely an athermal recombination corrected displacements per atom (arc-dpa) in the shell and a replacements per atom (rpa) representing atomic mixing in the core. These estimators are calibrated from explicit PKA simulations and a standard displacement damage model that incorporates cascade defect production efficiency and mixing effects. We highlight the predictability and accuracy of our reduced-order atomistic cascade method for the cases of copper and niobium by comparing its results with those from full PKA simulations in terms of defect production as well as the resulting cascade evolution and structure. We offer examples for simulating high energy cascade fragmentation and large dose ion-bombardment to demonstrate its possible applicability. Finally, we discuss the various practical considerations and challenges associated with this methodology especially when simulating subcascade formation and dose effects.
Formation of large volume fractions of Mn-Ni-Si precipitates (MNSPs) causes excess irradiation embrittlement of reactor pressure vessel (RPV) steels at high, extended-life fluences. Thus, a new and unique, semi-empirical cluster dynamics model was developed to study the evolution of MNSPs in low-Cu RPV steels. The model is based on CALPHAD thermodynamics and radiation enhanced diffusion kinetics. The thermodynamics dictates the compositional and temperature dependence of the free energy reductions that drive precipitation. The model treats both homogeneous and heterogeneous nucleation, where the latter occurs on cascade damage, like dislocation loops. The model has only four adjustable parameters that were fit to an atom probe tomography (APT) database. The model predictions are in semi-quantitative agreement with systematic Mn, Ni and Si composition variations in alloys characterized by APT, including a sensitivity to local tip-to-tip variations even in the same steel. The model predicts that heterogeneous nucleation plays a critical role in MNSP formation in lower alloy Ni contents. Single variable assessments of compositional effects show that Ni plays a dominant role, while even small variations in irradiation temperature can have a large effect on the MNSP evolution. Within typical RPV steel ranges, Mn and Si have smaller effects. Furthermore, the delayed but then rapid growth of MNSPs to large volume fractions at high fluence is well predicted by the model. For purposes of illustration, the effect of MNSPs on transition temperature shifts are presented based on well-established microstructure-property and property-property models.
Tungsten is the main candidate material for plasma-facing armour components in future fusion reactors. Bombardment with energetic fusion neutrons causes collision cascade damage and defect formation. Interaction of defects with helium, produced by transmutation and injected from the plasma, modifies defect retention and behaviour. Here we investigate the residual lattice strains caused by different doses of helium-ion-implantation into tungsten and tungsten-rhenium alloys. Energy and depth-resolved synchrotron X-ray micro-diffraction uniquely permits the measurement of lattice strain with sub-micron 3D spatial resolution and ~10 -4 strain sensitivity. Increase of helium dose from 300 appm to 3000 appm increases volumetric strain by only ~2.4 times, indicating that defect retention per injected helium is ~3 times higher at low helium doses. This suggests defect retention is not a simple function of implanted helium dose, but strongly depends on material composition and presence of impurities. Conversely, analysis of W-1 wt% Re alloy samples and of different crystal orientations shows that both the presence of rhenium, and crystal orientation, have a comparatively small effect on defect retention. These insights are key for the design of armour components in future reactors where it will be essential to account for irradiation-induced dimensional change when predicting component lifetime and performance.
Gold is a noble metal typically stable as a solid in a face-centered cubic (FCC) structure under ambient conditions; however, under particular circumstances aberrant allotropes have been synthesized. In this work, we document the phase transformation of 25 nm thick nanocrystalline (NC) free-standing gold thin-film via in situ ion irradiation studied using atomic-resolution transmission electron microscopy (TEM). Utilizing precession electron diffraction (PED) techniques, crystallographic orientation and the radiation-induced relative strains were measured and furthermore used to determine that a combination of surface and radiation-induced strains lead to an FCC to hexagonal close packed (HCP) crystallographic phase transformation upon a 10 dpa radiation dose of Au 4+ ions. Contrary to previous studies, HCP phase in nanostructures of gold was stabilized and did not transform back to FCC due to a combination of size effects and defects imparted by damage cascades.
Radiation-induced segregation (RIS) poses a significant challenge for ferritic Fe–Cr alloys under irradiation, as it can compromise mechanical integrity and increase susceptibility to intergranular corrosion. Yet, the mechanisms governing Cr segregation remain incompletely understood. Here, in this study, we present a physics-based rate-theory model parameterized using self-consistent mean field theory-based Onsager transport coefficients to investigate RIS at the grain boundary in dilute Fe-(0.1 at. %) Cr. Under equal production rates of vacancies and self-interstitial atoms (SIA), and their equal absorption rates by bulk dislocations, the model simulates the experimentally observed transition from Cr enrichment at low temperatures to depletion at higher temperatures. Under these unbiased conditions, systematic investigation reveals that while temperature-dependent transport properties dictate the segregation direction, dose rate, grain size, and dislocation density only influence the magnitude and spatial extent of Cr segregation. However, under more realistic conditions of preferential vacancy production within damage cascade and/or preferential SIA absorption by bulk dislocations, the enrichment-to-depletion transition shifts to lower temperatures. Our findings demonstrate that RIS predictions based solely on transport coefficients are valid only under symmetric point defect flux conditions, and that biases in defect production and absorption must be considered for accurate predictions. This work provides a mechanistic framework for understanding RIS in ferritic alloys and informs alloy design for advanced nuclear systems.
The nondestructive investigation of single vacancies and vacancy clusters in ion-irradiated samples requires a depth-resolved probe with atomic sensitivity to defects. The recent development of short-pulsed positron beams provides such a probe. Here, we combine depth-resolved Doppler broadening and positron annihilation lifetime spectroscopies to identify vacancy clusters in ion-irradiated Fe and measure their density as a function of depth. Despite large concentrations of dislocations and voids in the pristine samples, positron annihilation measurements uncovered the structure of vacancy clusters and the change in their size and density with irradiation dose. When combined with transmission electron microscopy measurements, the study demonstrates an association between the increase in the density of small vacancy clusters with irradiation and a remarkable reduction in the size of large voids. This, previously unknown, mechanism for the interaction of cascade damage with voids in ion-irradiated materials is a consequence of the high porosity of the initial microstructure.