Dynamic evolution of mixed materials bombarded with multiple

Dynamic evolution of mixed materials bombarded with multiple ions beams and impurities Tatyana Sizyuk Ahmed Hassanein School of Nuclear Engineering Center for Materials Under e. Xtreme Environment Purdue University PFC community meeting at UCLA August 4 -6, 2010, Los Angeles CA

Outline n ITMC-DYN Computer Simulation Package Collision-induced reactions + Diffusion, Segregation, Chemical erosion, Surface recombination, and Desorption n Benchmarking for ELM averaged regimes n Benchmarking for TEXTOR discharges n PWI and material issues for ITER performance n Future code extension: (1) Chemical reaction rates + detailed processes of erosion chemical (2) 3 D target structure for nonideal surface erosion n Summary 2

ITMC-DYN Integrated Models n A. Hassanein, “Surface effects on sputtered atoms and their angular and energy dependence”, Fusion Technology 8 (1985) 1735. n A. Hassanein and D. L. Smith, “Elastic and inelastic surface effects on ion penetration and the resulting sputtering and backscattering”, Nuclear Instruments and Methods in Physics Research B 13 (1986) 225. n T. Sizyuk and A. Hassanein, "Dynamic analysis and evolution of mixed materials bombarded with multiple ions beams", J. Nucl. Mater. , available online 1 July 2010 n T. Sizyuk and A. Hassanein “Dynamic analysis of mixed ion beams/materials effects on the performance of ITER-like devices“, to be published J. Nucl. Mat. (2010) 3

ITMC-DYN Integrated Modeling q Collision models are integrated with detail models of time-dependent processes such as atom diffusion and segregation. Therefore, we consider actual fluences and exact irradiation times of laboratory experiments. q Time step is determined from processes of diffusion, molecular recombination, and surface segregation. q Additional surface boundary conditions are used to specify the recombination of a diatomic hydrogen isotope molecules and their release from the surface in modeling of both laboratory and reactor conditions. q Surface segregation models are implemented and helped to explain recent experimental results. q Chemical erosion of carbon is calculated in each layer near the surface, depending on carbon concentration. Erosion yield is compared with that predicted by well-defined and benchmarked formula. 4

Simulation Results of ITMC-DYN q Influence of binding energy on surface build-up and recession 2. 0 e. V binding energy of Cs Effect of potential and fluence on surface growth/erosion of Cs (on Si surface) with 0 e. V binding energy of Cs in modeling 5

Testing Dynamic Update of Materials Composition q Nano-scale multilayers formation – 90 e. V Fe+ and 50 e. V C+ -> Si High-resolution TEM micrograph of Fe/C multilayers with 4 bilayer periods: 0. 75 nm, 1. 25 nm, 2. 5 nm and 5. 0 nm in order from the Si substrate, respectively. ITMC-Dyn – 4 bilayer periods: ~1. 1 nm, ~ 1. 8 nm, ~ 2 nm, ~ 4. 3 nm. (each from 2 bilayers) 6

Testing Dynamic Update of Materials Composition q Nano-scale multilayers formation – 90 e. V Fe+ and 50 e. V C+ -> Si Sung-Yong Chun, “TEM-examination and computer simulation of nanoscale multilayers by pulsed cathodic arc deposition”, Journal of Ceramic Processing Research. Vol. 4, No. 3, pp. 115~117 (2003) ITMC-Dyn – 4 bilayer periods: ~1. 1 nm, ~ 1. 8 nm, ~ 2 nm, ~ 4. 3 nm. (each from 2 bilayers) 7

Benchmarking Multiple Beams on Targets q Beams of hydrogen ions containing carbon impurity ions were impinged on pure tungsten samples heated from 653 to 1050 K. q Ion beams consisted of 70% of 333 e. V H+, 10% of 500 e. V H+, 20% of 1000 e. V H+ and from 0. 1% to 1% of C ions with 1000 e. V energy. q Carbon concentration on the surface of tungsten was measured and conditions for blisters formation were investigated in dependence on the total fluence, on % of carbon concentration in ion beams, and on samples temperatures. 1) Y. Ueda, T. Shimada, and M. Nishikawa, "Impacts of carbon impurities in hydrogen plasmas on tungsten blistering", Nucl. Fusion 44 (2004) 62 2) T. Shimada, T. Funabiki, R. Kawakami, Y. Ueda, M. Nishikawa, "Carbon behavior on tungsten surface after carbon and hydrogen mixed beam irradiation", J. Nucl. Mater. , 329– 333 (2004) 747– 751 8

Blisters formation and dependence on target temperature q Blisters formation depends on mobility of implanted gas and on concentration of gas atoms in the target q Accumulation of H in tungsten target in the experiments can be due to lattice imperfections, ion-induced damage (that is not applicable for low energies of H ions), and carbon deposition q In these experiments carbon was deposited in the form of graphite or in binding with tungsten (WC or W 2 C) on the surface of the samples 9

Modeling of Ions/Target Interactions (ITMC-DYN) Beam of H with C Impurities on W surface Fluence 3 x 1024/m 2 0. 8 % of C T. Shimada, T. Funabiki, R. Kawakami, Y. Ueda, M. Nishikawa, "Carbon behavior on tungsten surface after carbon and hydrogen mixed beam irradiation", J. Nucl. Mater. , 329– 333 (2004) 747– 751 ITMC-Dyn modeling with diffusion coefficient of H in dependence on carbon concentration in W target and for temperature of target 653 K 10

Modeling of Ions/Target Interactions (ITMC-DYN) Beam of H with C Impurities on W surface Fluence 3 x 1024/m 2 0. 8 % of C It was reported in the experiment that at 453 K a high densities of blisters formed ITMC-Dyn modeling with diffusion coefficient of H in dependence on carbon concentration in W target and for temperature of target 453 K 11

Modeling of Ions/Target Interactions (ITMC-DYN) Beam of H with C Impurities on W surface Fluence 3 x 1023/m 2 0. 8 % of C It the experiment blisters appeared under this fluence and target temperature 653 K ITMC-Dyn modeling with diffusion coefficient of H in dependence on carbon concentration in W target and for temperature of target 653 K 12

Dependence of blisters formation on impurity concentrations The lower content of carbon of 0. 11% in comparison with 0. 84% resulted in the following differences in samples: q Peak carbon concentration changed from 40% to more than 60% q its location is shifted 15 -20 nm from the surface in the case of 0. 11% and 78 nm in the case of 0. 84% q target erosion is almost 10 times higher under the ion fluence containing 0. 84% C The difference in the target erosion explains the difference in peak location. Carbon atoms, redistributed inside the W target by incoming hydrogen ions, are then removed by surface erosion and this process is accelerated in the case of 0. 84% C concentration in beam. 13

Modeling of Ions/Target Interactions (ITMC-DYN) Beam of H with C Impurities on W surface q Conclusion: Low level of impurity contents in plasma can significantly affect erosion lifetime, largely increase hydrogen isotope retention, and enhance bubble/blister formation in candidate reactor materials in fusion environments. q Evidence from other experiments: “It is assumed that tungsten carbide formed on the W surface under exposure to the carbon-seeded D plasmas serves as a barrier layer for diffusion and prevents the outward transport of deuterium, thus increasing the D retention in the bulk of tungsten. ” V. Kh. Alimov, et al. , “Deuterium retention in tungsten exposed to low-energy, high-flux clean and carbon-seeded deuterium plasmas”, J. Nucl. Mater. 375 (2008) 192 14

Effects of segregation and etching on carbon concentration in target q Surface segregation moves carbon atoms from bulk to surface layer because of increasing chemical potential of carbon in the bulk of W/C compound. q Implementing surface segregation can shift peak of carbon concentration profile closer to the surface q However simulating this processes self-consistently with ions beam deposition, target atoms sputtering, and atom cascade redistribution have the effect that surface segregation loses its effectiveness 15

Influence of Post Irradiation Processes – Surface Segregation, Diffusion and Etching q Post-ion-deposition surface segregation, and subsequent processes of etching of samples q Etching is done using 1. 7 ke. V Ar ions in experiments for measurement of carbon concentration profile q This can explain the reason of the enhanced carbon accumulation on the surface 16

Test of Diffusion Modeling Carbon deposition, knock-out by H and diffusion in W Self-consistently 17

Influence of Post Irradiation Processes – Surface Segregation and Diffusion 18

Influence of Post Irradiation Processes – Etching by 1. 7 ke. V Ar ions 19

Benchmarking TEXTOR discharges n Textor limiter was exposed at a radial distance R = 48. 5 cm from plasma center and 2. 5 cm behind the last close flux surface (LCFS) during 125 plasma discharges with a total duration of 559 plasma seconds n Edge plasma parameters during the experiment were measured with fast probe and He-beam diagnostics. The ion fluence on the limiter estimated using He-beam data had a value of 4. 9 · 1019 ion/cm 2 averaged over the area of castellation n Surface temperature of the limiter was varied from 200 C to 400 C which agrees well with values expected for upper vertical targets of ITER divertor [A. Litnovsky, et al. , J. Nucl. Mater. 367– 370 (2007) 1481] 20

Benchmarking TEXTOR discharges At carbon concentration in beam of 3%, carbon layer covers tungsten surface and can prevent hydrogen diffusion to the bulk Parameters for modeling were used based on fast probe measurements: 120 e. V D ions, 300 e. V C ions Parameters for TEXTOR experiments: Ion flux - 2. 5 x 1017 D/sm 2/s; 1%, 2% and 3% of C in the beam; 600 K temperature in W 21

Blistering at similar to TEXTOR conditions 1 x 1017 D/cm 2/s flux ‘clean’ and carbonseeded D plasmas (ion fluence of about 2 x 1020 D/cm 2); 200 e. V/D Impurity concentration in the near surface layer: V. Kh. Alimov, et al. , “Deuterium retention in tungsten exposed to low-energy, high-flux clean and carbon-seeded deuterium plasmas”, J. Nucl. Mater. 375 (2008) 192 22

Integrated Modeling issues for PWI n Erosion, re-deposition Mixed ions beams with impurities, erosion of compounds, coupling with above surface plasma codes n D and T retention Hydrogen diffusion, molecular recombination and desorption, chemical erosion, saturation in C and Be, agglomeration and blistering in W n He retention He diffusion, bubbles and blisters formation, W surface modification n Surface morphology Physical sputtering and impurities deposition in dependence on surface conditions n Mixing and alloying Chemical reactions rates for compounds formation V. Philipps, 19 th PSI San Diego 2010 23

Code extension: reaction rates for compounds formation at first wall materials n Elementary reactions, equations for reaction fluxes and experimental reaction parameters (sample of Be and C compound): formation reaction fluxes 2 Be + C Be 2 C → 2 Be + C destruction reaction fluxes Solving rate equations: Be 2 C formation: ΔE=1. 8 e. V (exp. ), k 0 = 10 -29 m 4/s Be 2 C dissociation: ΔE=3. 0 e. V (from ΔHf), k 0 = 1013 s-1 Ch. Linsmeier, 19 th PSI San Diego 2010 24

Code extension: detailed processes of chemical erosion Three processes determine the chemical erosion of carbon under low-energy hydrogen bombardment: n Reaction of thermalized ions within implanted surface proceeds via the hydrogenation of carbon atoms at edges of graphitic planes to CH 3–C complexes n Thermal reaction is enhanced by radiation damage induced in material that provides open bonds for hydrogen attachment. n At low surface temperatures: kinetic hydrocarbon emission or surface effects, since hydrocarbon radicals are bound to surface with much less binding energy (≈1 e. V) than are carbon atoms in their regular lattice environment (7. 4 e. V) J. Roth, et al. , Nucl. Fusion 44 (2004) 25

Code extension: 3 D target structure for exploring nonideal surface erosion For exploring of: Sputtering and deposition in W fuzz Shin Kajita, 19 th PSI San Diego 2010 Surface roughness effects on impurities deposition Y. Ueda, J. Nucl. Mater. 390– 391 (2009) q This will be accomplished by significantly increasing mesh size and target structure with parallel implementation of code running at ANL Supercomputers 26

Summary & Conclusion q A new dynamic version of ITMC code is developed (ITMC-DYN) q This version is capable of dynamically simulating target composition changes and all physical processes as a function of irradiation time q New models and simulation of diffusion, segregation, chemical erosion, and etching processes are developed and implemented q Initial results of modeling were presented and compared quite well with recent laboratory experiments of mixed ion beams and the effect of blister formation in tungsten material q Future work ($$ ? ) will include more detailed reaction processes and simulating more complicated experiments, real reactor conditions, and to propose/recommend laboratory experiments to enhance our understanding of mixed materials behavior and effects 27

Thank You 28

Examples of Simulation Results q Influence of binding energy on surface build-up and recession W. Eckstein et al. , ”Ion-induced alkali-silicon interfaces: Atomistic simulations of collisional effects”, Nucl. Instr. Meth. B 119 (1996) 477486 ITMC-DYN – Effect of potential and fluence on surface growth/erosion of Cs (on Si surface) with 0 e. V binding energy of Cs in modeling 29

Rates of mass losses Fluence 3 x 1024/m 2; target temperature – 653 K; C/H = 0. 84% 30

n n n n PWI and Material Issues for ITER Performance Erosion, re-deposition T retention He retention Surface morphology Cracking Melting and melt stability Mixing and alloying V. Philipps, 19 th PSI San Diego 2010 31

Code extension: detailed processes of chemical erosion n Total methane erosion can be calculated as sum of temperature dependent erosion YT (Mech model) n And kinetic part Y 0 (Hopf model for temperature 300 K) We will implement also in model: n Inventory of different hybridization states of carbon atoms (sp 2, sp 3 and sp 3 H) n Control of upper limit of H atoms bounded to C (1/3) n Control of maximum sp 3 centers as function of temperature 32