Computational Theory and Modeling Relevant to Fusion Reactor Materials Roger E. Stoller in collaboration with: Yuri Osetskiy, Haixuan Xu, Stanislav Golubov, Alexander Barashev Materials Science and Technology Division Oak Ridge National Laboratory

ICFRM-17 Tutorial Session: Fusion Reactor Materials

Aachen, Germany 11 October 2015

ornl Why do we care about modeling in structural materials? • Although irradiation experiments cannot be replaced by modeling alone, a purely experimental approach to understanding the effects of irradiation is also not practicable . costs for design and execution of reactor irradiations . costs of post-irradiation examination of radioactive materials . declining facilities for both irradiation and examination . combinatorial problem: broad range of materials, phenomena, and irradiation conditions - coolants, temperature, loading conditions, dose rate, i.e. because we have to … . dose because we have to • Recent advances in computational hardware, computational materials science, ... make it more feasible than ever to aggressively pursue the contribution of modeling see relevant chapters in: Comprehensive Nuclear Materials, R. J. M. Konings, T. R. Allen, R. E. Stoller, and S. Yamanaka, Editors, Elsevier (2012) ornl Components of Primary Radiation Damage Source Term

• incident particles of different types and energies produce different displacement production types of primary damage from elastic collisions – , heavy charged particles, , • produce differences in secondary damage accumulation • need to know both energy spectrum and absolute flux level

nuclear transmutation nuclear transmutation from from thermal neutrons, high-energy neutrons, E<1 eV threshold reactions, E>2 MeV

ornl What is unique about modeling the fusion environment? • the DT fusion reaction source term, unique aspect is 14.1 MeV neutrons higher energy atomic recoils for displacement damage production • higher levels of He and H from transmutation via threshold reactions, e.g. ~12 appm He/dpa and 50 appm H/dpa in AISI-316 for typical first note: many lower energy wall spectrum neutrons as well

ornl Early mechanistic models used to anticipate role of transmutant helium on swelling: potential fission-fusion difference • Helium production higher in fusion than LMFBR (FFTF) environment – experiments with isotopic tailoring of Ni in austentitic steel used to obtain intermediate and high helium in water moderated reactors (ORR and HFIR) – experiments confirmed model predictions of possible non- monotonic behavior Stoller and Odette, 1982

ornl Relevant phenomena with related computational and experimental methods

ornl Provide brief examples from various aspects of the multiscale modeling scheme • ab initio (VASP), development of interatomic potentials, He-Fe potential to account for He defects in iron • primary damage formation, and atomistic kinetic Monte Carlo • mesoscale (reaction rate theory) model of microstructural evolution – illustrate loosely coupled multiscale modeling (parameter passing) – 2D cluster dynamics – compare alternate kinetic models • molecular dynamics (atomistic) simulation of - defect interactions • displacement rate effects and vs. irradiation presentation focuses on bulk radiation effects, other unique aspects for plasma-facing materials ornl Application of ab initio calculations: Development of He-Fe potential • Wilson’s He-Fe pair potential developed in 1960s, Defects Fitted

still extensively used Unrelaxed Heint oct • A new He-Fe inter-atomic potential has been Unrelaxed Heint tet developed at ORNL Unrelaxed He • Extensive fitting to first-principles calculations of sub He octahedral point defects and clusters int

• First-principles calculations in VASP Heint tetrahedral

• Fitted to both relaxed and unrelaxed defects Heint midpoint oct-tet – Both forces and energies fitted Hesub • Must predict tetrahedral site as most stable site for He-He-V cluster He – Not possible with pair or EAM-style functional forms He-He-He-V cluster

– Hence 3-body potential needed Heint-Heint pair

• Large clusters e.g. He3V – Up to 128 in VASP

ornl ORNL three-body He-Fe potential

• Considerable effort invested to properly understand and implement 3-body potential: complex dependence

Energy of Fe-He-Fe triplet as a function of angle

Qjik and rik where rij is equal to 1.8 Angstrom

2 Y(rij ,rik ,Qijk )  cos Qijk  0.44fcut rij fcut rik  j,kFe

ornl Interstitial He diffusion

6x10-8 • Fe Potentials:

– AM: Ackland-Mendelev /s) – A97: Ackland-Bacon 2 – FS: Finnis-Sinclair 10-8 • He-Fe Potentials

– OR: ORNL 3-body – JN: Juslin-Nordlund – W: Wilson Fe Fe-He E (eV) D (10-8 m2/s) m 0 • Results with ORNL Fe-He AM JN 0.062 9.64 -9 FS ORS 0.065 3.15 potential lies between those of 10 A97 ORS 0.073 3.27

other Fe-He potentials Hediffusion coeficient (m AM ORS 0.043 2.23 AM W 0.107 4.69 5x10-10 • ORNL potential shows 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 consistent results with different T /T (K-1) Fe potentials m

ornl Formation of Interstitial clusters

• Yellow atoms are He, BCC Fe lattice not shown. • He-He in Fe has strong binding – This leads to formation of interstitial clusters • Mobility decreases with size • When cluster is large enough, it creates Frenkel pair, emitting SIA – new He-vacancy cluster immobile – He and He-v clusters are new source of Frenkel pair – SIA formation leads to formation of small loops

ornl Comparison of potentials: He clustering • He diffusion and clustering properties • Behavior of different He-Fe potentials – Wilson potential creates FP (SIA emission) more quickly than ORNL potential, higher He binding energy. – Juslin-Nordlund potential does not form He clusters or create FP above ~400°C. • Provide parameters for mesoscale microstructural models

ornl He transport with ORNL potential

• Single interstitial helium atoms are very mobile – but easily captured by vacancies

• SIA can recombine with HexV if x small – creates mobile interstitial helium • Interstitial helium tends to form clusters • Small He clusters are mobile – The smaller, the more mobile • Large He clusters create Frenkel pair to increase local volume, reduce pressure – transformation to He-vacancy cluster, nascent bubble – He-vacancy cluster is immobile – ejected interstitials trapped by the He-V cluster and combine to form dislocation loops

ornl Molecular dynamics simulation of primary damage • MD simulations provide opportunity to investigate displacement cascade evolution, e.g. effects of lattice, PKA energy, T, etc. • Classical molecular dynamics, typical implementations: – many millions of atoms, solve Newton’s equation of motion – constant pressure or volume, periodic boundary condition – system may or may not be thermostated to prevent PKA from heating system – no electronic losses or -phonon coupling, energy of cascade simulation, hence for:, • to compare with standard NRT displacement model:

nNRT = 0.8Tdam /(2 Ed), Tdam= kinetic energy lost in elastic collisions

– EMD ~ Tdam (NRT) < EPKA

– e. g. for Tdam = 100 keV, EPKA = 175.8 keV and nNRT = 1000 (see appendix on use of SRIM for dpa calculations)

ornl Comparison of representative neutron and corresponding PKA energy spectra

• differences in neutron flux level lead to different atomic displacement rates • neutron energy spectrum differences lead to different PKA energy spectra – different coolants, water for HFIR and PWR vs. sodium for FFTF alter neutron energy spectrum, primarily influence lower energy – high energy influenced by neutron source, c.f. all fission with ITER fusion

ornl MD cascade database for iron • an extensive cascade database has been developed that covers a broad range of cascade energy and temperature – up to 200 keV at 100, 600, and 100 keV at 900K • this database includes a sufficient number of simulations at each condition to provide a good statistical measure of average cascade behavior – total number of point defects produced – in-cascade clustering fractions for both interstitials and vacancies – in-cascade cluster size distributions – Note: of cascade event leads to better statistics for defect production than for clustering • common use of Finnis-Sinclair potential provides basis for comparison other investigations: – pre-existing damage (relevant to cascade overlap) using 10 keV cascades – free surfaces (relevant to in situ experiments) using 10 and 20 keV cascades – nanograined iron using 10 and 20 keV cascades

ornl Time dependence of defect evolution in atomic displacement cascades

Movies 10keV

50keV

ornl Cascade-induced shock wave

• Three distinct phases in cascade evolution: – Supersonic – Transonic – Sonic Calder, et al. (Phil. Mag. 90, 2010) ornl The role of shock waves in defect cluster formation

February 27-March 3, 2011 – San Diego, ornl California Angular dependence of displacement threshold energy: effect of crystalline lattice on defect formation classical MD

Fe V

Zepeda-Ruiz, PRB 67, 2003 Bacon, et al., JNM 205, 1993 ab initio MD, P. Olsson, C. Domain, EDF R&D, unpublished

ornl

ornl If normalized to NRT displacements, MD results show reduced defect survival as cascade energy increases

- some of curve structure is significant, related to cascade morphology and subcascade formation - note small standard errors, measure of mean behavior - effect of temperature, 100 to 900K, is systematic but not strong ornl Many of surviving defects are in clusters formed during the cascade event

20 keV, 600K

Note: poorer statistics, larger standard errors, than for total defect survival • significant in nucleation of extended defects • not accounted for by dpa (or f >x MeV)

ornl Illustration of subcascade structure at peak damage condition for cascades at 100K

100 keV 10 keV

5 keV

• high energy cascades look like multiple lower energy events, leads to asymptotic behavior with energy • low energy events between subcascades have higher efficiency

ornl Influence of energy and temperature on SIA clustering

Stoller, JNM 276, 2000

ornl Vacancy clustering in iron • MD simulations reveal little vacancy clustering in the nearest-neighbor sense - but, vacancies are spatially correlated Typical uncollapsed vacancy cluster, 50 keV cascade at 100 K

arrows indicate 4nn • kMC aging of residual defects indicate that such loose or nascent clusters do tend to collapse into void-like configurations

ornl Influence of energy and temperature on vacancy clustering

Stoller, JNM 276, 2000

ornl MD Results: Application to Neutron Spectrum Effects • weighting of PKA spectrum using MD defect production • little effect of neutron energy spectrum observed in spectrum- averaged defect production (similar for other primary damage parameters) • yields effective source term for mesoscale models Point defect survival, per NRT

ornl MD results: averaged over neutron (pka) energy spectrum (from SPECTER/SPECOMP) T=100K Finnis-Sinclair ~proportional to potential absorbed energy only see Stoller and Greenwood, JNM 271 & 272 (1999) and ASTM STP Significance to 1366, 2000. damage accumulation?

ornl Overall Summary of MD Cascade Simulations • Cascade energy or PKA energy dependence of stable displacement production is more complex than standard NRT model • Subcascade formation plays a dominant role in controlling cascade morphology and stable damage production • Clusters formed directly in the cascade account for a substantial fraction of interstitials, many of which are in glissile configurations – interstitial clustering fraction and cluster sizes increase with both cascade energy and temperature • Significantly less in-cascade vacancy clustering appears, nascent clusters coalesce at longer times, e.g. in KMC simulations of cascade debris – vacancy clustering increases with cascade energy and decreases with temperature • Nearby free surfaces (relevant to in situ experiments), pre-existing damage (cascade overlap) and nanograined microstructure alter defect survival and clustering behavior Stable defect formation is the result of a complex process involving energetic and coupled many-body reactions

ornl See slides added at the end that include additional information on damage formation

• more details on defect clustering and cascade statistics • comparison of Fe with other metals • effect of pre-existing damage on cascade defect production • effect of free surfaces on cascade defect production • differences in cascade defect production in nanograined material • Appendix on recommended approach for using SRIM to compute dpa

ornl Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC)

H.X. Xu, Y.N. Osetsky, R.E. Stoller, Physical Review B, Brief Reports, 84, 132103 (2011) SEAKMC is particularly powerful for large systems with complex defects; accurately includes defect diffusion, defect interactions naturally occur ornl Active Volumes (AVs) – Spatial Localization

Small defects are localized, this can be exploited to speed calculations

Saddle point searches are carried out only within the AVs, different defects have different AV size and shape

Dumbbell Vacancy

ornl Short-term defect evolution using new self-evolving atomistic kinetic Monte Carlo (SEAKMC) method* The previously (for ~50 years) unexplained mechanism of <100> interstitial loop formation in iron has been determined for the first time using the SEAKMC method Results provide a direct link between experiments and atomistic simulations and new insights into defect evolution. − Formation process of <100> loops involves a distinctly atomistic interaction between two ½ <111> loops. − Reaction does not follow the conventional assumption of Burgers vector conservation between the reactants and the reaction product; process is different from all previously proposed mechanisms. − Results observed in SEAKMC using multiple interatomic potentials and simulation cell sizes; predictions at higher temperatures could be confirmed using molecular dynamics with much greater computational effort.

H. Xu, et al., Phys. Rev. Let. 110 (2013) ornl Application of mesoscale models

• mesoscale models are relevant to many phenomena in materials science and radiation effects – grain growth – dislocation evolution, by thermo-mechanical or radiation-induced processes – void swelling – precipitation of additional phases, and solute segregation – stress corrosion cracking, and irradiation-assisted SCC • size scale permits direct comparison with experiments such as TEM and mechanical property measurements • dependent on fundamental atomistic processes, and controls macroscopic observables such as strength, ductility, creep, ... • a primary application is the investigation of point defect and solute kinetics and microstructural evolution, includes both mean field models and Monte Carlo simulations based on reaction rate theory

ornl Example from reaction rate theory modeling

• starting point is continuity equations describing point defect (vacancy, Cv and interstitial, Ci) populations (analogous equations for solutes): D C C v v T v D C + ------U + G – C C – D C S = v v kT v v i v v v v t

DiCi T Ci D C + ------U + G – C C – D C S = i i kT i i i v i i i t where the  denote spatial derivatives.

• first term on the LHS describes point defect drift to discrete sinks, the Ui,v are interaction energies between the point defects and discrete sinks em • G iv  =  G dpa + G iv  is the total point defect generation rate, including thermal emission from sinks, and the D iv  are the point defect diffusivities T • S are the total sink strengths for continuum sinks (e.g. cavities, iv dislocation, grain boundaries, etc.), recombination rate coefficient is given in terms of an effective recombination radius, r : r  = 4rrDv + Di

ornl • typical assumptions/simplifications in mean-field models – the material is treated as a spatially-homogeneous effective medium with embedded effective sinks and sources for point defects – spatially-averaged point defect generation rates – these assumptions have been relaxed in particular cases, e.g. to investigate cascade-induced fluctuations in point defect concentrations – the models are formulated as a series of differential equations describing the production and fate of point defects and the corresponding evolution of the microstructure • With these approximations, and assuming that the irradiation produces only monomers, the time-dependent or steady state point defect concentrations can be obtained as a solution to the following equations:

dCiv em T ------= G + G – C C – D C S dt dpa iv i v iv iv iv

ornl • analogous equations can be written to describe an evolving point defect cluster population, for helium generation and distribution, and for the other microstructural components, e.g.: – if only the monomers are assumed to be mobile, an equation describing the di-interstitial population, C2i, can be written:

dC2i i 2 3i 2i 2i ------=  C +  C C – C  C +  C dt i i v v 3i 2ii i v v

– this equation can be generalized to formulate a master equation for interstitial (and vacancy) clusters of arbitrary size. However, this also generates an arbitrarily large number of equations - smallest visible defect clusters ~100 point defects – various grouping schemes have been devised to minimize number of equations, introduce potential source of error, most accurate is by Golubov, et al., Phil. Mag. A, 2001. – modern computers reduce driving force for limiting number of equations

ornl • rate equations can be written to describe the evolution of extended defects, e.g. cavity sink strength and growth rate: v v T 0.5 S = Z 4rN 1.+ r S iv iv v viv drv 1 v v v ------= -----Z D C – C – Z D C v vv v i i i dt rv • analogous expressions available for grain boundaries, , dislocation loops, etc. • greater or lesser detail can be built in as needed to simulate given phenomena e.g. defect nucleation vs. growth regimes, dislocation evolution, effects of solute segregation, etc.

ornl • solutions obtained from simultaneous integration of equations included in a given model and, when well calibrated with experimental data, such models have some predictive capability

ornl However, the success of these models in fitting data can be deceiving (the devil is in the details ...) • data fitting with incomplete models leads to use of “effective” parameter values, use of parameter-rich models limits confidence in model extrapolation • for example, if point defect absorption dominated by dislocations

with sink strength Sd, swelling rate is proportional to product of  d and net dislocation bias (Zi -1):

dV Gdpa d ------= ------Z – 1 d i dt S v • In this case, data fitting can not be used to obtain unique set of model parameters. MD cascade simulations provide independent estimate of  and thereby permit better estimate of d Zi

ornl MD cascade simulations provide independent estimate d of  and thereby permit better estimate of Zi

• more complex models may be ‘stiffer’ with respect to arbitrary parameter choices, but more complex models introduce additional parameters

ornl Cluster dynamics simulation of bubble evolution under irradiation: two-dimensional master equation df(,,) x m t Jx(  1,) mJxm  (,)  Jxm (,  1)  Jxm (,) dt x x m m Mechanisms include: • Radiation-induced generation of point • Migration, clustering, and absorption of defects and He atoms point defects and He atoms • Frenkel pair creation by over- • Radiation-induced resolutioning of He pressurized He-vacancy clusters atoms • Evaporation of He atoms and f(x, m+1) m+1 vacancies

P (x) Q (m+1) He He • Discrete equations up to x=m=10

P (x-1) P (x) v v • Grouping scheme of Golubov, et al. m f(x-1, m) f(x+1, m) Q (x) + Q (x+1) + i i (Phil. Mag. A81 (2001) 643) for larger Qv(x,m) Qv(x+1,m) sizes PHe(x) QHe(m) m-1 f(x, m-1)

x-1 x x+1

ornl Movie to illustrate 2D calculations

ornl Bubble size distributions for nominal input parameters

He bubble size distributions obtained at 100 dpa for irradiation temperatures of 200°C and 330°C

ornl Comparison of kinetic models • Mean field reaction rate theory (MFRT) and object kinetic Monte Carlo (OkMC) have a common basis in reaction rate theory • MFRT has a long history in radiation effects modeling; OkMC has seen increasing use since recent advances in computational power have expanded the accessible domain for these models • Since MFRT and OkMC models are similar kinetic models – can be used simulate the same phenomena – many details are handled differently in the two approaches • A direct comparison of the MFRT and OkMC models – point defect cluster dynamics modeling, relevant to nucleation and evolution of radiation-induced defect structures – illustrated relative strengths and weaknesses of the two approaches (see Stoller, et al. J. Nucl. Mater. 382, 2008)

ornl MFRT-OkMC: Inherent differences parameter or MFRT OkMC mechanism solution method deterministic stochastic inferred from processes and time explicit variable reaction rates smeared, effective space medium, possible full spatial dependence multi-region MFRT time and space- averaged, but c.f. defect production discrete in time and space Mansur’s cascade diffusion model sink strength, e.g. explicit mathematical inferred from fate of point dislocations expression defects limited (computationally) by defect or sink essentially unlimited simulation cell size, i.e. N ≥ density 1/(x·y·z) ornl Implications for comparison of simulation results: Making sure you are simulating the same problem

• deterministic vs. stochastic: not possible to do exact point-to- point comparisons • spatial correlations when simulating cascade defect production: lost in MFRT, how to handle in OkMC – electron irradiation with only FP production most straight forward – spatial relationship of point defects and clusters if use cascade results as input in OkMC • minimum achievable defect/sink density in OkMC

– e.g. for bcc iron, x=y=z=300ao, V=6.42x10-22 m-3, 1/V=1.57x1021 m3 18 24 -3 – c.f. typical defect densities in irradiated materials = 10 to 10 m – imposes limitation on range of irradiation temperature (low) and displacement rate (high)

ornl ornl ornl Summary: MFRT-OkMC Comparison • Common issues, primarily related to parameterization and basic model formulation – e.g. binding and migration energies for small point defect clusters, nature of defect migration (1-D vs 3D) – are all (or all relevant) mechanisms and reactions accounted for? • Both approaches produce similar results on well-posed problems, similar dependence on key irradiation and material variables • Issues for MFRT – account for effects of discrete time and spatial events, e.g. primary damage • Issues for OkMC: primarily computational – limited to relatively low doses – system size limitatation constrains available regime of irradiation temperature and dose rate • MFRT-OkMC differences remaining (see Stoller, et al. JNM 382, 2008) – cluster densities when statistical nucleation difficult – damage accumulation under cascade damage conditions – uncorrelated vs correlated cascade debris ornl MD-Based Dislocation Dynamics: Objectives - Understand atomic-scale details of dislocation interactions in complicated environment (point defects, impurities, large obstacles, temperature, variable loading conditions, etc.) Particular issues : - mechanical property changes (strengthening, embrittlement, etc.) - mechanisms of hardening due to localized obstacles - mechanisms of microstructure modification by moving dislocations, e.g. cleared channels

ornl MD-based dislocation dynamics

Y X - [111] Y - [110] Z X Z - [112]

ornl Direct comparison of MD dislocation dynamics simulation with in situ TEM experiment

73nm

Simulation and experiment reveal dynamic interaction of a gliding dislocation with stacking fault tetrahedron: critical role of free surface in experiment

Importance of simulation accurately simulating the experiment

ornl

160 2nm at T=300K Results for dislocation – 140 bubble interactions Strong • obstacle strength is a non- monontonic function of radius

120

(MPA) MAX  100 Intermediate

80 0.0 0.5 1.0 1.5 2.0 Weak

He/Vac 24 ) void 20

<110> He/Vac=0.25 a 16 He/Vac=0.5 He/Vac=0.75 12 He/Vac=1.0 8 He/Vac=2.0

4 along<110> ( 0 5 10 15 20 25 30 35 40 along <112>, (nm) ornl (Partial) list of issues in understanding how to use charged particle irradiations • charged particle vs. neutron

– a few orders of magnitude in Gdpa available – effect of surface, greater for in situ electron irradiation – effect of gradients – PKA spectrum – implanted , chemistry and interstitials – lack of transmutatation products, can be “simulated” • heavy ion vs. electron – PKA spectrum, lack of cascades (clusters) in electron irradiation – surface influence in electron irradiation • neutron vs. neutron

– less range in Gdpa available – PKA spectrum effects, more modest that in ion/neutron comparison • neutron-neutron comparisons are least ambiguous 56 ornl Comparing ion and neutron irradiation: Packan and Farrell, JNM 78, 1978, pure nickel

57 ornl Simple model-based comparison • Mansur’s equation for temperature shift (JNM 78,1978):

• Old Stoller-Odette model, ASTM STP 955 75K • displacement rate varied in model • Mansur Eqn. gives 73K • shift of peak vs shift of cut-off • shift will be sensitive to ~175K m Ev • "shift" may be different at different doses 58 ornl Comparing ion and neutron irradiations • ... its complicated • when microstructure is sensitive to dose rate, it is reasonable to expect a temperature-shift in radiation-induced phenomena • relative complexity and stability of unirradiated microstructure is likely to influence T-shift, e.g. FM vs austenitic stainless steels • new models need to account for influence of mobile interstitial clusters under cascade-producing irradiation • point defect partitioning is the name of the game, whether dislocation bias, production bias, high sink densities, ... • when assessing models and their predictions vis a vis experimental results, need to think carefully about interpretation – model limitations, missing mechanisms – “effective” parameters – the art is in knowing what matters

59 ornl Overall Summary

• Substantial progress in understanding and predicting the behavior of materials has been provided by theory and computational modeling • Further progress requires additional research: – Advances in electronic structure theory to obtain intrinsic and defect properties in iron and its alloys, including the effects of He and H – Computationally-efficient, physically-robust interatomic potentials for multicomponent alloys, including effects of directional bonding and magnetism – Advanced atomistic and mesoscale models describing the defects and processes that interact in complex ways in multicomponent, multiphase alloys – “Properly” linked and multiscale (atomistic, mesoscale, and continuum) deformation and fracture models for predicting hardening, plastic instability, changes in ductile-to-brittle fracture, dimensional instability, and creep/creep rupture behavior under realistic loading conditions • It is critical to validate models* using data from special-purpose and engineering experiments, objective is a predictive capability to avoid technical surprises, e.g. new damage mechanisms at high doses and long times, how to use ion irradiation to understand neutron irradiation * need to think about model assumptions and what the simulation results mean ornl

ornl Additional slides on:

Cascade statistics, morphology, and other materials

Appendix on calculating dpa:

ornl

ornl Comparison of MD defect formation in a range of bcc, fcc, and hcp metals

Bacon, et al., JNM 276, 2000

• different values, but similar power-law depending in defect survival • energy dependence of SIA clustering similar, varies by a factor <2

ornl Comparison of 10 keV cascades in iron and copper based on molecular dynamics simulations

Fe Cu

• Similar embedded type interatomic potentials, edge length of

simulation cells is 50ao. • Note higher level of in-cascade clustering in more-compact copper cascade. • Defect survival (relative to NRT) is lower in copper than in iron. ornl Comparison of 20 keV cascades in iron and copper based on molecular dynamics simulations

Fe Cu

• Similar embedded atom type interatomic potentials, edge length of

simulation cells are 50ao. • Note higher level of in-cascade clustering in more-compact copper cascade. • Defect survival (relative to NRT) is lower in copper than in iron. ornl Influence of cascade structure on time signature

• systematic differences in the time dependence of defect formation are observed around the time of peak damage

ornl Influence of cascade structure on time signature

• subcascade formation reduces influence of pressure wave that creates many small-range, transient displacements • note that stable defect formation is reduced for structures with higher peak values ornl Formation of point defect clusters

Calder, et al. (Phil. Mag. 90, 2010) have carried out detailed analysis of high-energy cascades. By tracking individual atoms, and local temperature and atom density, they show: • the formation of large SIA clusters in Fe is related to the formation of hypersonic recoils (>10 speed of sound) near the sonic front of the primary cascade • site of the SIA cluster is determined within ~0.1 ps • vacancy cluster collapse is a result of prior SIA clustering

ornl

ornl

ornl Major difference between high-energy cascades: interstitial cluster size distribution

ornl

ornl

ornl Pre-existing Damage: Effect on Defect Survival

ornl

ornl

R. E. Stoller, “The Effect of Free Surfaces on Cascade Damage Production in Iron,” Journal of Nuclear Materials 307-311 (2002) 935-940. R. E. Stoller and S. G. Guiriec, “Secondary Factors Influencing Cascade Damage Formation,” Journal of Nuclear Materials 329-333 (2004) 1228-1232. For extended study: see Y.N. Osetsky, A.F. Calder and R.E. Stoller, "How do Energetic Ions Damage Metallic Surfaces?," Current Opinion in Solid State & Materials Science 19 (2015) 277-286.

ornl

ornl

ornl

ornl Investigation of potential length-scale effects in primary damage formation

• The MD database provides a good description of primary radiation damage formation from atomic displacement cascades, but most is in perfect, single crystals. • Large number of mobile defects produced in a displacement cascade, nearby grain boundaries could potentially reduce the residual damage from any given cascade. • Previous work (e.g. Samaras, et al. in Ni) indicated that the high volume fraction of grain boundaries in bulk nanocrystalline materials could provide efficient point defect sinks and/or recombination sites.

ornl Approach

• Parallel version of MOLDY MD code using OpenMP on shared memory platform • Voronoi technique used to create nanocrystalline system – “Nucleation site” distribution based on an fcc lattice, melted through Monte Carlo using a Lennard-Jones potential • Choose nanograin “nucleation sites” → Fill in grains → Remove overlapping atoms → Equilibrate structure – Nanocrystalline system equilibrated ~200 ps • Simulate cascade event → Perform nearest neighbor defect analysis → Differentiate grain boundary reconstruction from in-grain defects – Cascade simulation was run for ~15 ps, well into the region where the defect count stabilize

ornl Nanocrystalline System: investigate length scale effects on primary damage • MD simulation cell, 100 lattice parameters (~28.6 nm) on edge, periodic boundaries • System shown contains 32 grains • ~10 nm grain size • Nanocrystalline system equilibrated ~200 ps • Representative cascade size: 5 keV 20 keV

ornl Defect Visualization

• Grain boundaries (GBs) distinguished using spherical approximation – 1st and 2nd nearest neighbor analysis – Before and after visualization can show GB movement or reconstruction • Possible defects flagged when more than 0.3 of a lattice parameter from an original atom site • Possible in-grain defects flagged when more than 0.5 a lattice parameter from original GBs

ornl Cascade Simulations

• Small study of cascade energy and temperature, 8 simulations at each condition: – 10 keV, 100K – 20 keV, 100K – 20 keV, 600K • Results compared with “normal” cascade events to determine difference in residual damage between single crystal and nanocrystalline iron cascades • Previous cascade database for comparison: – 10 keV, 100K: 15 events – 20 keV, 100K: 10 events – 20 keV, 600K: 8 events

ornl Comparison of Final Defect State:10 keV, 100K

Single Crystal Iron Nanocrystalline Iron

Residual defects - Interstitial atoms are green, vacant sites are red, grain boundary atoms within 3 lattice parameters of any defect are black.

ornl Comparison with 100 K database

Mean values are shown with the standard error. • Higher vacancy survival in nanograined material • Much lower interstitial survival in nanograined material

ornl Results: Stable defect production

Note: Ni≠Nv

Mean values are shown with the standard error.

• Higher vacancy survival in nanograined material • Much lower interstitial survival in nanograined material ornl Results: In-cascade clustering

Much less in-cascade interstitial clustering in nanograined material

ornl Defect production in nanograined Fe • Strong influence of microstructural length scale (grain size) on primary damage production • Reduced interstitial survival and clustering will reduce formation of radiation-induced microstructural components • Excess vacancy production may lead to higher supersaturation, greater propensity for cavity formation • Impact of altered primary damage behavior needs to be evaluated over longer time scale, e.g. using mean field rate theory or Monte Carlo models • Further analysis of vacancy clustering and grain boundary motion is underway

ornl Appendix on calculating dpa: meaning and limitations, and relationship with damage correlation

ornl • observed radiation damage depends on some measure of radiation exposure, “damage flux,” f>[0.1, 0.5, 1.0] MeV, but ions, electrons, photons?? • need to account for differences in PKA energy spectra to be able to correlate data from different types of irradiation • secondary displacement model by Norgett, Robinson and Torrens, Nucl. Engr. and Des. 33 (1975); based on earlier work by Kinchin and Pease (1955), damage partitioning model by Lindhard, et al. (1963) [see ASTM standards E521 and E693]

• number of displacements, nNRT, proportional to fraction of PKA energy \ deposited in elastic collisions, Tdam: 0.8T n  dam NRT 2 E d • The best way to think of the NRT dpa is to focus on Tdam as an exposure parameter – a measure of absorbed energy – does not account for anything other than the potential total number of atomic displacements, e.g. no information on in-cascade formation of point defect clusters – does not account for transmutation production – does not account for any effects due to ionization number of dpa may or may not correlate with any particular observed effect of radiation exposure ornl Compare NRT, TRIM, and MD displacements

1000 800 500

see Stoller, et al., NIMB 310, 2013 for recommendations on the use of SRIM for computing dpa ornl Number of displacements from TRIM • varies with method chosen: – use integral of output file: vacancy.txt or “Total target vacancies” – compute ~damage energy by integrating output files – Kinchin-Pease vs. “full cascade” mode • if reporting dpa for comparison with neutron irradiation – use Kinchin-Pease – use standard displacement threshold energy, i.e. 40 eV for iron (see ASTM E521) – set lattice and surface binding energy to 0.0 PKA energy NRT damage NRT Average MD, TRIM K-P TRIM Full (keV) energy (keV) displacements 100K Cascade 1.0 0.81 8 --- 7.7 to 9.4 ---

78.7 50 500 168 533 to 540 566 to 572 529* 1052-1075* * upper row from alternate damage energy integrations, lower row from vacancy.txt

ornl