Microstructure of a heavily irradiated metal exposed to a spectrum of atomic recoils

At temperatures below the onset of vacancy migration, metals exposed to energetic ions develop dynamically fluctuating steady-state microstructures. Statistical properties of these microstructures in the asymptotic high exposure limit are not universal and vary depending on the energy and mass of the incident ions. We develop a model for the microstructure of an ion-irradiated metal under athermal conditions, where internal stress fluctuations dominate the kinetics of structural evolution. The balance between defect production and recombination depends sensitively not only on the total exposure to irradiation, defined by the fluence, but also on the energy of the incident particles. The model predicts the defect content in the high dose limit as an integral of the spectrum of primary knock-on atom energies, with the finding that low energy ions produce a significantly higher amount of damage than high energy ions at comparable levels of exposure to radiation.


INTRODUCTION
Metals exposed to irradiation by highly energetic particles develop nanoscale structural distortions through the ballistic displacement of atoms out of their crystal lattice sites.In the limit of high exposure, these radiation defects become numerous enough to significantly alter the microstructure, leading to detrimental changes in materials properties, such as reduction in thermal and electrical conductivity, volumetric swelling and dimensional changes, hardening, and embrittlement.The development of a quantitative simulation algorithm for modelling evolving fluctuating microstructures of metals and alloys in a radiation environment has recently acquired prominence in the context of virtual reactor design [1][2][3][4] , as failure resulting from the accumulation of radiation damage limits the service life-time of reactor components -both in the presently operating fission reactors 5 and in conceptual fusion reactor designs.In particular, it is deemed essential to be able to identify the difference between microstructures produced by exposure to various types of energetic particles, for example neutrons with various energy spectra 6 or energetic ions 7 used as cost-effective surrogates for neutron irradiation experiments.
Significant progress has been made recently in the theoretical development and experimental validation of quantitative models predicting the number of defects generated in a perfect crystalline matrix by an individual recoil atom resulting from a collision with an incident high energy particle [8][9][10] .However, these models do not describe the non-linear microstructural changes resulting from the accumulation of defects occurring beyond the low dose limit of exposure of 0.01 displacements per atom (dpa), where the overall density of defects exceeds approximately ∼ 0.01 at%.In the low dose limit, the clustering of defects is well described by the power law statis-tics, discovered in simulations 8 and confirmed by electron microscope observations 11 .These power laws are similar to those found in observations of fragmentation of solid projectiles occurring on impact 12 .In a heavily irradiated metal, the fact that the defects cluster directly in the cascade events appears less significant, since in a material already containing structural distortions, clustering also occurs due to elastic interaction between defects, irrespective of whether they accumulate sequentially or form simultaneously in collision cascades [13][14][15] .Eventually at doses above approximately 0.1 dpa, clusters of defects start coalescing, and this gives rise to the formation of an extended system-spanning dislocation network that, together with the population of isolated defects continuously produced by irradiation, forms a dynamically fluctuating driven steady state of the material [16][17][18] .
In this study, we develop a quantitative model for predicting the defect content in a heavily irradiated metal, informed by simulations of consecutive cascade impacts extending into the high dose range ∼1 dpa.Whilst in general, microstructure is expected to depend on temperature and radiation histories in an intricate and tangled manner, we find that in the athermal limit, where microstructure is driven towards a steady state by defect generation and stress relaxation and not by thermal diffusion, it is possible to formulate a predictive model for the defect content in a heavily irradiated material, analogous to but at the same time fundamentally distinct from models used for predicting defect generation in a perfect crystal lattice.
We find that radiation drives metals to different dynamic steady states containing different amounts of defects, depending on the spectrum of recoil energies, and that the defect content can be reversibly changed by altering the recoil spectrum.We predict that the recoil spectra involving a larger fraction of high energy particles produce microstructures with lower vacancy content, in qualitative and quantitative agreement with experimental assessments of vacancy concentration in tungsten irradiated by high energy ions.

Properties of the steady state
The generation of radiation defects is simulated by explicitly propagating atomic trajectories from the moment an incident high energy particle collides with an atom in a crystal lattice, often called the recoil or primary knockon atom.If sufficient kinetic energy E R is transferred to the recoil atom, it is displaced from its lattice site, leading to a cascade that proceeds through three distinct stages of evolution, see Fig. 1.In the ballistic phase lasting approximately ∼ 0.1 ps, the recoil atom initiates more recoils, which themselves may initiate more recoils, and so on, thereby causing a cascade of self-similar displacements, until no atoms remain with kinetic energy sufficient to create more recoils.This is followed by the heat spike phase of duration ∼ 1 ps, caused by the redistribution of the remaining kinetic energy of displaced atoms to the surrounding crystal lattice, which produces a localised region with average atomic energies well over the melting point.In the final cooling phase extending to ∼ 10 ps, the ions cool down and recrystallise.The recrystallisation process concludes on a time scale comparable to that of atomic motion, and therefore crystal defects may remain in the volume of the former heat spikethese are the aforementioned radiation defects.
The state of the art modelling of high dose microstructure at atomic resolution is the sequential simulation of collision cascades 13 .Whereas a single collision cascade in a conventional size simulation cell introduces the dose of ∼ 10 −4 dpa 17 , higher levels of exposure can be reached by repeating thousands of such cascade simulations successively in the same evolving microstructure.It is computationally feasible to reach high doses of φ ∼ 1 dpa using this algorithm 18,19 , however, the total simulated time is unlikely to extend far beyond 1 µs for system sizes of interest due to the computational limitations of molecular dynamics.The simulated dose rates of dφ/dt ∼ 1 dpa/µs are many orders of magnitude above even the comparatively high dose rates of dφ/dt ∼ 1 dpa/hour characteristic of ion irradiation experiments.
Under such extreme dose rates, thermally activated processes that contribute to damage recovery in a simulated microstructure are partly or wholly suppressed because the defects involved are statistically more likely to be reconstituted by a collision cascade before they get a chance to undergo any thermally activated transformation.Reversing the argument, one may also conclude that for a specific thermally-activated process and a given dose rate, there is a temperature below which the process is effectively frozen out and does not occur, thereby not contributing to microstructural evolution.It is in this athermal regime that consecutive cascade simulations offer a valid description of a heavily irradiated material.
It is possible to define the athermal regime more pre- thermal factor thermal factor Estimating the phase space of the athermal regime.a Simulation and model for the mean number of defects and molten atoms generated by a cascade with given damage energy in tungsten.The standard error is smaller than the marker size.b Thermal factor (5) for dose rates dφ/dt representative of experiment and of cascade simulations.With decreasing temperature, thermally-activated processes with higher energy barriers become frozen out (ξ 1), as any given region in the crystal becomes more likely to be impacted and molten by a cascade before a thermally activated transformation occurs.For pure tungsten, E melt = 2.70 eV and the transformation barriers are: for vacancy migration E V m = 1.52 eV, for self-climb 20 2E V m = 3.04 eV, and for vacancy-mediated climb 21 cisely for a thermally-activated process characterised by a given activation energy E A , for example vacancy migration.The mean time between the events is given by Arrhenius' law where ν 0 is the attempt frequency approximately equal to the Debye frequency divided by 2π.In a system containing N atoms, assume that a collision cascade occurs on average after every time interval t c .While the heat spike phase of the cascade does not produce a true equilibrium liquid phase, there is significant mixing possible during this time interval 22 , and it is a convenient simplification to consider the heat spike as the melting and recrystallisation of N melt atoms.In this scenario, the mean time between the successive cascade impacts at any atomic site is given by According to the formal definition of dose, a single cascade increments it by the ratio of the number of atoms displaced in the ballistic phase N d to the system size: ∆φ = N d /N .Eliminating t c using the expression for dose rate φ = N d /(N t c ), we define the thermal factor ξ as the ratio of mean times (2) to (1): Next, we use the estimate 23 , where E d ∼ 10-100 eV is the threshold displacement energy, defined as the average kinetic energy required to displace an atom out of its crystal lattice site averaged over crystallographic directions 24 , and the damage energy T d < E R is the remaining kinetic energy after accounting for the energy lost to electronic excitation 25 .Finally, we parameterise the size of the molten cascade region by the phenomenological expression ) where E melt is the effective energy per atom required to melt the crystal 26 , as obtained from fits to single cascade simulations similarly to defect production models 9 , see Fig. 2a, and E min d is the threshold displacement energy minimum with respect to crystallographic directions.For the damage energy region of interest T d > E d , the thermal factor can be simplified to The thermal factor is the mean number of times a defect takes part in a thermal migration event before being enveloped by the next cascade.In the athermal regime ξ 1, a defect is much more likely to be reconstituted by a cascade than undergo a thermal hop, while in the thermal regime ξ 1, a defect is able to migrate many times, transferring mass through the system or coalescing with other defects, before being embroiled in another cascade.Examples of thermal factors for various processes in pure tungsten are shown in Fig. 2. For vacancies migrating in pure tungsten, equation (5) shows that ξ 1 at room temperature and dose rates representative of ion irradiation experiments, suggesting that consecutive cascade simulation indeed offers a valid description of microstructural evolution under these conditions.

Defect content in a driven steady state
In the athermal limit, where all the thermally-activated processes are frozen out, irradiated microstructure develops through the balance between defect generation by cascades and defect recombination during the recrystallisation of cascade heat spikes.Considering the fact that irradiation by highly energetic particles generates a broad spectrum of recoil energies, one poses a question if the Pristine tungsten is damaged by 100 keV cascades until saturation.c-d Switching to 100 eV cascades drives the system towards a steady state with higher vacancy content.e-f Switching back to 100 keV cascades returns the system to the previous steady state with lower vacancy content through formation and coalescence of vacancy loops.Green (pink) lines show the position of 1/2 111 ( 100 ) dislocation lines detected using the DXA method 27 .
damage produced by the cascades with a given recoil energy can be reversed by cascades with a different damage energy, or whether the crystal structure retains memory of its irradiation history.
Figure 3 illustrates microstructural evolution of the initially pristine tungsten that is successively damaged to saturation first by 100 keV energy cascades, followed by 100 eV energy cascades, and finally again by 100 keV energy cascades.The vacancy content, shown in figure 3g, is found to saturate to different levels depending on the recoil energy: switching from 100 keV to 100 eV cascades drives the vacancy content from 0.32 at.% to ∼1 at.%, after which switching back to 100 keV cascades drives the vacancy content back to 0.32 at.%.The final vacancy content is such as if the 100 eV cascade exposure had not occurred at all.
We observe that the vacancy content saturation to a steady state is entirely determined by the recoil energy, irrespective of the irradiation history.While the migration barrier for some interstitial-type defects is sufficiently small to allow for their athermal migration 28 , driven by the fluctuating stress fields developing in the irradiated microstructure 14 , vacancies and vacancy clusters are effectively immobile in the athermal limit of microstructural evolution, remaining stationary between the successive cascade impacts.In the high dose limit, every part of the initially crystalline material has molten and recrystallised at least once, with the average cascade impacts generating a fixed number of vacancies and a fixed mean number of recrystallised atoms, depending on the damage energy.
Based on the above reasoning, we propose the following expression for the average vacancy concentration c sat V in the athermal steady state driven by collision cascades with damage energy T d : where N sat V (T d ) is the number of vacancies remaining after recrystallisation of N sat melt (T d ) atoms in the steadystate microstructure.These quantities are essentially unknown, and their determination would require extensive analysis of collision cascades in the already saturated microstructures.However, under the assertion that the microstructure of the athermal steady state is qualitatively well described as a crystal with highly elevated vacancy content, see figure 3d, the models for cascades in a perfect crystal, surprisingly, might also be able to describe the saturated steady state.Using the above expression for the size of the cascade melt in an initially pristine crystal (4), we also make use of the athermal recombination corrected displacement model (arc-dpa) introduced by Nordlund et al. 9 , in the form modified by Yang and Olsson 10 , to obtain an estimate for the number of vacancies generated by a single cascade: where function b + c is a phenomenological expression accounting for the sub-linear defect production at intermediate damage energies, with parameters b and c fitted such that N arc d accurately matches the number of defects obtained from molecular dynamics simulations for a given material 9,10,29 .
The central outcome of this analysis is an approximation for the formally unknown number of vacancies and recrystallised atoms in the saturated regime by their counterparts in a perfect crystal: N sat V ≈ N arc d and N sat melt ≈ N melt .Hence, we arrive at a simple expression for predicting the saturated vacancy concentration c sat V in a high-dose microstructure where N melt and N arc d are given by Eqs. ( 4) and ( 7), respectively.Expression ( 8) is equivalent to stating that in a heavily irradiated material in the athermal limit, the history of microstructural evolution is completely irrelevant for the vacancy type defects.We do not discuss here the self-interstitials, because while formally the vacancy and interstitial counts must be equal, interstitials exhibit a strong drive towards coalescence 30 , and the eventual formation of extended dislocation networks 14 or even complete crystal planes 15,18 , and as such may not be readily recognisable as point defects.
We have validated the vacancy concentration model ( 8) by running cascade simulations in tungsten up to the dose of 1 dpa over a broad range of recoil energies.We find that vacancy concentration saturates well below the dose of 1 dpa, see Fig. 4a, in agreement with independent simulations 13,14,16 and experimental observations 31,32 .This can be rationalised by recognising that the number of molten atoms per cascade is much higher than the number of ballistically displaced atoms, as suggested by their respective characteristic energy scales T d /E melt T d /E d .In Fig. 4b we compare the saturated vacancy concentration with the analytical prediction (8), and find that the prediction is in qualitative and quantitative agreement with the simulation results.Most notably, we find that impacts with smaller recoil energies that are characteristic of irradiation by light particles, produce a significantly higher vacancy content than impacts with high recoil energies, characteristic of irradiation by heavy particles.
Predictions of the athermal saturated vacancy content can be made for a variety of materials, as parameters for N arc d are available in literature for a broad selection of metals and alloys 9,10,29 , and as E melt can either be obtained from cascade simulations or estimated.We have tested the transferability of the model using face-centered cubic (fcc) copper and hexagonal close-packed (hcp) zirconium for which we performed cascade simulations using 100 eV, 1 keV, and 10 keV recoil energies up to 1 dpa, see table I.The predictions are in qualitative agreement with data from high dose cascade simulations.In principle, a few data points from high dose cascade simulations could be used to adjust the model parameters for better quantitative agreement, as shown in Fig. 4b Vacancy content as a function of recoil energy.a The vacancy concentration in tungsten irradiated up to a dose of 1 dpa, as obtained from cascade simulations, is found to decrease with increasing recoil energy.Shaded areas indicate the standard error obtained from multiple simulations.b Comparison of the vacancy concentration from simulation in the ∼1 dpa high dose limit (points) with the analytic prediction (line) given by the number density of vacancies generated in a cascade in an initially pristine crystal.Fitting the model by improves quantitative agreements (dashed), but removes the predictive aspect of the model.Mean values are averaged over 0.9-1.0dpa, with error bars showing the standard deviation, or in the case of multiple simulations, the mean standard error, over the same interval.
we emphasise that high dose simulations require several orders of magnitude more computational resources than the single-cascade simulations used for parameterising the predictive model; generating the single-cascade data shown in Fig. 2 required only ∼5,000 core-hours, compared to ∼1,000,000 core-hours required for accumulating the high dose data shown in Fig. 4.
In the general case where there is a distribution of damage energies p(T d ), the probability for an atom in the crystal to become part of a heat spike of an incident cascade is equal to p hit (T d ) = p(T d )N sat melt (T d )/N , where N is the total number of atoms.The mean saturated vacancy concentration, averaged across a volume region of interest, is therefore obtained as the expectation value of the vacancy concentration c sat V (T d ) over the distribution after some simplification using Eq. ( 6).

Experimental validation
As performing a direct quantitative measurement of vacancy content in an irradiated material is highly nontrivial 40 , we opted for an indirect comparison to depthresolved deuterium concentration measurements.We consider deuterium retention in irradiated tungsten as, due to the expected use of tungsten for plasma-facing components in fusion reactor designs, the data derived from such experiments are available in the literature over a broad variety of irradiation conditions.In the deuterium retention experiments considered here, the material is first damaged by ion irradiation and then exposed to a deuterium plasma at elevated temperatures T ∼ 400 K over the course of several days.As the deuterium particle energies are much smaller than the threshold displacement energy, the deuterium is only implanted into a shallow surface layer from where it subsequently diffuses into the irradiation-damaged layer of the material whilst binding to trapping sites.First-principles simulations found that the deuterium binds preferentially to void surfaces, with a single tungsten vacancy trapping up to 5 deuterium atoms with trapping energies above 1.1 eV 41 .We compiled nuclear reaction analysis measurements of peak deuterium concentrations in tungsten irradiated to doses φ 0.3 dpa under athermal condition T < 500 K for a broad range of irradiation energies from a number of studies 32,[36][37][38][39] .For each particle energy represented in the experimental data, we simulated the recoil spectra, see Fig. 5, and used them to compute the expected saturated vacancy concentration using Eq. ( 9).In Fig. 5b we compare the predicted deuterium concentrations computed as c sat D = 5 c sat V 16 with the experimentally measured deuterium concentration in the saturation regime.The predictions closely follow the trends found in experiment.We find that ions that have a higher probability of generating smaller cascades also generate microstructures characterised by higher deuterium retention, which affirms the prediction that low energy cascades correspond to higher large-dose-microstructure vacancy concentrations, see Fig. 4. We also computed the deuterium concentrations using the fitted model shown in Fig. 4b, and found that the results differ only by a mean relative error of 7 %.The fitted model does offer better accuracy, see the Supplemental Material, but also requires data from computationally intensive high dose simulations, which on balance appear unwarranted since there are other comparatively less controlled sources of uncertainty present here.
Although the experiments 32,[36][37][38][39] were conducted at temperatures below which vacancies become mobile, referred to as stage III recovery 42 , some lower-barrier interstitial migration processes may still contribute towards the interstitial-vacancy recombination.We do not expect this effect to play a significant role in tungsten, where interstitial migration barriers are so low that they are by the fluctuating stress fields present in the microstructure under irradiation.An additional source of systematic error is our assumption that the entire vacancy-type content produced in a cascade is represented solely in the form of mono-vacancies, ignoring the possibility of vacancy clustering.From the comparison with our high dose cascade overlap simulations, we find that even with the uncertainty introduced by using an empirical interatomic potential, we can still estimate the deuterium content to around 10% of the experimentally measured value.
Deuterium retention experiments for neutronirradiated tungsten under athermal conditions are comparatively scarce.Shimada et al. 43 report a value of 0.8 % retention for tungsten irradiated to 0.3 dpa at 70 °C in the High Flux Isotope Reactor (HFIR) reactor, with deuterium subsequently implanted at 200 °C.Using a recoil spectrum representative of HFIR 44 , we predict a deuterium concentration of 1.5 %.It is encouraging that the prediction overestimates the retention, given that the experiment was conducted at a temperature where some extent of defect recombination through vacancy migration is expected.

CONCLUSION AND OUTLOOK
The exposure of metals to high dose irradiation at relatively low temperatures eventually leads to a saturation of defect content.The actual magnitude of this content depends on the primary recoil energy, with low energy cascades generating up to 5 times higher vacancy concentrations than high energy cascades.Our central result is that the number density of vacancies generated within the characteristic volume of the cascade can be used as a predictor for the saturated vacancy concentration, valid over many orders of magnitude in recoil energy spectra.The model is transferable across different metals, provided that the irradiation conditions are athermal, meaning that the dose rate is high relative to the W 20.3 MeV (1) W 12.3 MeV (2) Cu 9 MeV (3) Fe 9 MeV (4) W 5.5 MeV (5) Au 5 MeV (6) Si 7.5 MeV (7) He 1 MeV (8) H 3 MeV (9)   rate of thermal relaxation of defect configurations.The model prediction is validated by comparison to experimental data on the saturated deuterium concentration in highly irradiated tungsten over a broad range of ion energies and masses.
The saturated vacancy concentration acts as a single number quantifying radiation damage in simplest terms, and can be used to predict the outcome of irradiation with different particle energies and masses a priori.This is analogous to how total exposure to radiation is measured in units of displacements per atoms, transferable across different materials and computed from the entire recoil spectrum.While the prediction of the saturated damage state is formally valid in the particular regime of irradiation conditions where thermal defect migration is negligible, this prediction further serves as a rigorously defined baseline for characterising deviations from the steady-state, for instance when thermal diffusion is active.

METHODS
Molecular dynamics cascade simulations.MD simulations were run using lammps 45 with interatomic interactions described by empirical potentials of the embedded atom model type: W (Ref: 33 ), Cu (Ref: 34 ), and Zr (Ref: 35 ).Simulation cells were initialised with perfect crystal structure, with periodic boundary conditions applied to all three directions and box dimensions chosen consistent to the 0 K equilibrium lattice constant.Starting with an initial dose of 0 dpa, irradiation damage is introduced following an iterative procedure: First, a small number of atoms N C are selected at random and their kinetic energy set to the recoil energy, with randomly-oriented velocities.If N C > 1, we avoided selecting random atoms too close to one another in order to avoid spurious effects caused by overlapping cascades.Next, the simulation is propagated in the NVE ensemble including some damping terms, see below, until the crystal cools below 100 K and at least 5 ps are simulated, after which velocities are zeroed out and atomic coordinates relaxed to a local energy minimum using the method of conjugate gradients.This procedure is repeated until a target dose is reached, with each repetition incrementing the dose by 0.8T d N C /(2E d N ), where N is the number of atoms in the simulation cell.We chose N C such that the dose is incremented by ∼0.0002 dpa per iteration as a compromise between computational efficiency and avoiding excessive heating of the crystal.System sizes varied between 65,000 and 2,000,000 atoms, with smaller cells and multiple repetitions used for lower energy recoils.The vacancy content was determined using the method described in Ref. 16 .Microstructural images were rendered with ovito 27 .
Energy loss to electronic excitation was included by adding a frictional force of magnitude mv/τ to all atoms with kinetic energy higher than 10 eV, where m is the atomic mass and τ is the time damping constant extracted from the low velocity selfion stopping limit of the srim code 46 .Electronic stopping time constants are 0.90 ps, 1.04 ps, and 2.37 ps for W, Cu, and Zr, respectively.For atoms with temperature below their melting temperature, another damping term was added to model energy loss due to electron-phonon coupling 47 with time constants of 16.0 ps, 26.4 ps, and 12.5 ps for W, Cu, and Zr, respectively.We refer to the Supplemental Material for more information on how these parameters were obtained.
Parameters for the arc-dpa and melt models for the interatomic potentials used here were obtained from single cascade simulations and are listed in the Supplemental Material.The size of the melt was determined using ovito 27 by filtering out atoms of the stable crystal structure type as determined by the common neighbour analysis method, after which atoms with less than 3/4 of the equilibrium coordination number are filtered out to remove the crystalline region surrounding the melt.The remaining atoms have a liquid-like radial distribution function.The melt size is then obtained as the maximum size of the melt over the course of a cascade.
Recoil spectra simulations.Recoil energies were obtained using the srim software 46 , which is based on the binary collision approximation method.Recoil energies were converted into damage energies using the Lindhard stopping formula 25 as described in Ref: 23 .The damage energies were binned, and bin counts divided by the bin width and by the total number of recoils, to obtain a discrete rep-resentation of the normalised probability distribution p(T d ) for a recoil to have the damage energy T d , shown in Fig. 5a.The recoil spectrum for the HFIR reactor is shown in the Supplemental Material.

FIG. 1 .
FIG.1.Displacement damage in the driven dynamic steady state of a material exposed to radiation.Damage evolution after initialisation of a cascade with 30 keV recoil energy in a-c pristine tungsten and e-g tungsten previously irradiated to 1 dpa.Shown are only the crystal defects identified using the Wigner-Seitz defect analysis as interstitials (red) and vacancies (blue).d,h Evolution of the defect content produced by cascade initiation.In pristine crystalline tungsten, the defect content increases linearly with each cascade event, while in the already irradiated tungsten, the defect content fluctuates around a mean value; no new defects are generated as the microstructure is in a steady-state with respect to irradiation.

FIG. 3 .
FIG. 3. Reversibility of driven high-dose microstructure.a-bPristine tungsten is damaged by 100 keV cascades until saturation.c-d Switching to 100 eV cascades drives the system towards a steady state with higher vacancy content.e-f Switching back to 100 keV cascades returns the system to the previous steady state with lower vacancy content through formation and coalescence of vacancy loops.Green (pink) lines show the position of 1/2 111 ( 100 ) dislocation lines detected using the DXA method27 .
FIG.5.Saturated deuterium content in the high dose limit.a Distribution of cascade damage energies for various irradiation ions.Ions with higher momentum are more likely to generate higher energy cascades.b Comparison of analytically predicted (9) and experimentally measured deuterium concentrations in the saturated limit.Irradiation with lower momentum ions leads to a higher deuterium content, consistent with the higher vacancy content produced by smaller cascades.Colors and labels are consistent with damage energy spectra from Fig.5a.Measured concentrations extracted from (6) 32 , (1,3,4,7,8) 36 , (5) 37 , (2) 38 , (9) 39 .

TABLE I .
Saturated vacancy content in at.% for three types of crystal structures damaged by cascades with given recoil energy.Comparison between prediction (8) and simulation at 1 dpa shows that the model qualitatively reproduces the drop in vacancy content as the recoil energy increases.
5. Saturated deuterium content in the high dose limit.aDistribution of cascade damage energies for various irradiation ions.Ions with higher momentum are more likely to generate higher energy cascades.bComparison of analytically predicted (9) and experimentally measured deuterium concentrations in the saturated limit.Irradiation with lower momentum ions leads to a higher deuterium content, consistent with the higher vacancy content produced by smaller cascades.Colors and labels are consistent with damage energy spectra from Fig.5a.Measured concentrations extracted from (6) 32 , (1,3,4,7,8) 36 , (5) 37 , (2) 38 , (9) 39 .