Deuterium depth profiling with D( 3 He, p)α nuclear reaction

Transcription

1 Seminar 1 b - 1.letnik, II.stopnja Deuterium depth profiling with D( He, p)α nuclear reaction Author: Mitja Kelemen Advisor: doc. dr. Primož Pelicon Ljubljana, september 201 Abstract The seminar is intended to describe the method of nuclear reactions D( He, p)α to detect deuterium in materials. Following the presentation of the nuclear reaction used for the analysis I concentrate on the description of the concepts needed to calculate the proton yield. I present related cases of measurements at the tandem accelerator laboratory of Jožef Stefan Institute. The described examples of analysis include deuterium depth profiling in tungsten exposed to atoms of deuterium and He microbeam analysis. 1

2 Contents 1 Introduction 2 Methods of deuterium detection and depth profiling Thermal Desorption Spectroscopy (TDS) Secondary ion mass spectroscopy (SIMS) Elastic Recoil Detection Analysis (ERDA) Neutron-ERDA (NERDA) Nuclear reactions resonance method (NRA) Nuclear reactions 5 4 Reaction kinematics 5 5 Scattering cross-section 7 6 Energy loss and stopping power 7 7 Proton yield 9 8 Experimental measurements 10 9 Some measurements and results from Jožef Stefan Institute Upgrade of classical NRA at Jožef Stefan Institute 12 2

3 1 Introduction One of most critical issues in the construction of thermonuclear reactor based on magnetic confinement in the tokamak configuration is the construction of the inest wall. During plasma operation of such device, the wall is subject to severe physical conditions, a combination of neutron and charged particle bombardment, heavy and uneven thermal load photon irradiation and atomic hydrogen exposure (Figure 1). This leads to erosion, deposition, gas capture and lattice damage of the wall. The build up of deposited layers also leads to fuel retention in the vessel walls as the deuterium or tritium fuel may be co-deposited along with the eroded materials. This is particularly true for carbon, which has a chemical affinity for hydrogen isotopes such as deuterium or tritium. While fuel retention is not a problem for present day tokamaks which are using deuterium plasmas, it will become an issue for next step devices, where the allowable vessel tritium inventory is limited, especially if carbon would be used as the dominant plasma-facing material. This is the main driver for developing the tungsten plasma-facing materials in view of the deuterium-tritium operation of International Thermonuclear Experimental Reactor (ITER). ITER walls will be made of Be, divertor zone of W. The most critical plasma striking zones of divertor wall will be made from carbon composites[1]. The use of carbon fiber composites (CFC), molybdenum and tungsten as plasma-facing materials in nuclear fusion research resulted in an increasing interest in quantitative profiling of deuterium in plasma-facing materials. Deuterium can penetrate these materials up to depths of several tens to several hundreds of µm, where it is trapped at low concentrations [2]. The assumption is that this deep penetration is an important process for the deuterium inventory in the tokamaks and other fusions reactors[]. Study on temperature dependence of diffusion speed are plan to take place in the end of 201 by Jožef Stefan Institute and Max-Planck-Institut für Plasmaphysik collaboration. To keep within the ITER regulatory limit, the tritium inventory build up in ITER must be controlled. The buildup is extrapolated from present day machines where no serious inventory issues are foreseen if tungsten plasma-facing covers are used. Techniques for tritium removal and tritium inventory measurements are under active development since they could be needed in the case carbon plasma-facing covers are used during deuterium-tritium experiments. The carbon materials penetration is related to the porosity of the material, while in tungsten the diffusion process is driven by the stress-field induced by implanted deuterium [2]. Poor understanding of these processes, and their investigation requires methods for quantitative depth profiling of deuterium up to depths of at least 10 (or better few tens) µm. Figure 1: Plasma wall interaction scheme[1]

4 2 Methods of deuterium detection and depth profiling In many cases, samples are from the walls of experimental fusion reactors. Consequently, the applied depth profiling method should be as insensitive to surface roughness as possible, prohibiting grazing incident or exit angles for the ion beam analysis methods, such as ERDA. 2.1 Thermal Desorption Spectroscopy (TDS) In the Thermal Desorption Spectroscopy (TDS), the specimen is heated at a constant heating rate in a vacuum. The change in hydrogen partial pressure due to hydrogen release from the specimen is monitored as a function of temperature or time. By carrying out the TDS, desorption peaks corresponding to the trapping sites are visible in the obtained TDS spectrum, partial pressure vs. temperature curve. After acquiring several TDS spectrum at different heating rates, desorption energy for each trapping state can be deterrent.from this it is possible to deduce total inventory and binding energy of hydrogen isotopes[4]. However, no information on depth concentration profiles is obtained. 2.2 Secondary ion mass spectroscopy (SIMS) This method allows us to do depth profiling of hydrogen isotopes up to a depth of about 100 µm [5]. However, it has some disadvantages as it is time consuming and not quantitative by itself. A quantification can be obtained by comparison with ion beam analysis results for the surface layer. For good accuracy the ion-beam analysed layer should be as thick as possible, which again raises the question of quantitative hydrogen isotopes analysis up to larger depths[5]. 2. Elastic Recoil Detection Analysis (ERDA) With ERDA (Elastic Recoil Detection Analysis) one detects lighter particles from a target sample that are recoiled by heavier incident ions from the probing beam. When a MeV He ion beam incident on a sample with H and D, the H and D in the sample will be scattered in the forward direction. Spectrometry of these forward recoiled atoms gives rise to the quantitative depth profiling of these species in the sample[6]. Main limitation of this method is that it can be used only in thin films and not to thick samples. This results also in high sensitivity on the surface of the samples if proper geometry is used Otherwise, it produces quantitative depth concentration profiling of deuterium. 2.4 Neutron-ERDA (NERDA) NERDA with incident 14 MeV neutrons beam from the sample back side allows depth profiling of hydrogen isotopes up to depths of hundreds of µm [2].The principle of measurements and analyses is similar to normal ERDA, except the use of neutron instead of heavy ions. However, the depth resolution of this method is only a few tens to about 100 µm [2] and counting statistics is poor due to the low incident neutron fluxes, resulting in low sensitivity. 2.5 Nuclear reactions resonance method (NRA) The D(d, p)t and the D( He, p)α reactions can be used for quantitative depth profiling of D up to large depths. The ( He, p) resonance method, exploits the broad resonance of this reaction, Figure by measuring the total proton yield as function of incident energy and a calculation of deuterium depth profile becomes possible [2][7]. The resonance method uses only the total number of detected protons and can be utilized with detectors having large solid angles. However, due to the broad resonance with a FWHM of 70 kev[2] and the non-negligible contribution of protons originating from outside the resonance region this method requires measurements at 8-10 different energies, thus being time-consuming; it requires good counting statistics, which largely compensates the advantage of a large detector; and it is often difficult to show that the derived depth profile is unambiguous. The paper concentrates on the NRA method. It is used by several groups to create quantitative depth profiles of deuterium, some also up to very large depths ( 100µm). And it is also in use at Jožef Stefan Institute. 4

5 Nuclear reactions During the measurement, a beam of He ion (1+ or 2+) is directed to the sample containing deuterium. The goal is to determine how much deuterium is accumulated in the sample and its depth distribution. Incident helium ions are elastically scattered in the sample, as well as induce nuclear reactions on deuterium nuclei 2He D 4 2 He p 2He D 1 1 p n + 2 He 2He D 2 1 D D p 2He D 1 T p p Q = 18, 5 MeV Q = 2, 7 MeV Q = 6, 52 MeV Q = 1, 95 MeV (1) In addition, helium ions can also interagate with carbon atoms in sample, if sample is based on CFC, which can be described by the nuclear reaction in which protons are produced, 2He C 14 7 N p Q =, 96 MeV 2He C 1 1 p N Q = 9, 85 MeV All these reactions produce various products, which are then detected in the detection set-up. Usually high-energy protons are detected, which are best suited for measurements, as they are charged and lose small amount of energy when passing through the sample material. From the measured data and properties of these nuclear reactions (kinematic relations and scatter cross-section) it is possible to say something about the depth profile of deuterium in the analysed layer. It should be noted that nuclear reactions (1) and (2) are not the only possible mode of interaction of ions with the helium atoms in the sample, since the incident helium ions are elasticity scattered on nucleus in target 2He D 2 He D (2) 2He C 12 6 C + 2 He () 2He C 1 6 C + 2 He These reactions in terms of depth profiling of deuterium are forming a low-energy background in the measurement. Problem is solved by placing protective film (usually Mylar or Clapton) in front of the detector. Fast protons pass, while scattered helium ions are stopped in the foil. For practical experiments and measurements is mostly used first of the given nuclear reaction D( He, p)α. Consequently, we will focus on this reaction and its application in practice. From now on if we write He it means He. Due to practical reasons, we detect fast protons from the first of given D( He, p)α reactions, as they are fast and could be easy separated from other detected particles. 4 Reaction kinematics For the beginning we must calculate amount of energy Q released in reaction D( He, p)α, which is defined as the difference of kinetic energy initial and final stats: Q = T f T i = c 2 (M i M f ) = c 2 (m He + m D (m α + m p + 2m e )) = = ( )MeV = MeV (4) From here we see that the reaction is exothermic (Q > 0), which means that the initial state will always get the selected final state, no matter what kind of energy the incident helium ion. After the formation of new particles, this energy is converted into kinetic energy of the reaction products (in this case the proton and α particle) in center-of-mass system. 5

6 We are interested in how the kinetic energy of the proton depends on the energy of the incident He ion. this end, we solve the kinematic equations for a given nuclear reaction in the laboratory system[7]: To E He + E D = E p + E α p He + 0 = p p cos θ + p α cos φ 0 = p p sin θ p α sin φ where θ and φ are angles in the laboratory system, under which proton and a particle α fly after reaction is finished. The p x stands for the size of momentum three vectors of particles that appear in the reaction. When these equations are solved for the outgoing proton momentum[7], the result can be written as p p = E phe cos θ ± E 12 E 12 4m 2 p (E12 2 p2 He cos2 θ) 2(E12 2 (6) p2 He cos2 θ) wher E 12 and E 12 indicates E 12 = m D + p 2 He + m2 He E12 2 = E 12 p 2 He + m2 p m 2 α The expression (6) has two solutions in general, but it turns out that the only physical one, where we choose positive sign before the second article in the numerator. The kinetic energy of the proton in lab system is then equal to[7] T p = p 2 p + m 2 p m p (8) (5) (7) Figure 2: The kinetic energy of protons, resulting from a nuclear reaction D( He, p)α as a function of the scatter angle θ and energy of incident helium ions[7]. How T p depends on the kinetic energy of incident helium ions and the scatter angle θ is shown in the Figure 2. We can see that protons have incurred much higher energy then incident helium ions. This is due to the high Q, which is released in this reaction. Since the α particle is much heavier than the proton, proton takes the most of the energy. We obtain slow α particles and fast protons. From the resulting dependence of proton kinetic energy is also seen that the maximum proton kinetic energy when the proton is emitted in forward direction, and the lowest when it is emitted backwards. Our detector is positioned at 15. Higher incident energy of incident He result in lower detected energy of proton from ( He, p) reaction. 6

7 5 Scattering cross-section Because of the nature of nuclear reactions, i.e. the complex structure of matter and the forces which are in play in this reaction, the derivation of the scatter cross section is very complicated. Therefore, we rather lean to the experimental data. These were measured by numerous authors V.Kh. Alimov [8], W. Möller [9] and T.W. Bonner [10]. Which are otherwise given for different systems but all can be transformed in centre of mass system[7]. Figure : Differential scattering cross section for the reaction D( He, p)α the center of mass coordinate system. The points on these graphs are shown the measured scattering cross sections, which are measured by Alimov, Möller, Bonner. In addition to the full line marked phenomenological approximation of Möller and Besenbacher [7]. It turns out that the scattering cross-section shape is resonant, as shown in Figure. At this point it is also worth to point out that the bell-shaped ( resonance ) scatter cross-section shape is essential in profiling deuterium in the matter, because, as we shall see below, allows as to differentiate particles which are scattered at different depths in sample. 6 Energy loss and stopping power When scattering projectiles (in our case He ions) does not takes place on the surface of the sample, but rather deeper in the target, we have to consider that a particle loses a fraction of the energy on the path to the target nucleus with which will interact. Energy loss of the incident particle is due to particle Coulomb interaction with electrons in electron clouds on its path. Loss of energy per unit of travelled path of the particle is described by the stopping power S(E), which is defined as, S(E) = 1 de (9) ρ dz Energy dependence of the stopping power in the first approximation is described by Bethe-Bloch function [6], in general, is its dependence on energy is much more complex and we can not accurately calculate. In practice we use the average measured values by any combination of the projectile and the target, they can be found collected in the tables. 7

8 Figure 4: He ions energy at different depths and initial energies in a target of 100% W, calculated with SRIM 2008 code [11]. The result of stopping power is that He ions have only finate range in the selected target, display the Figure 4. Figure 5: The graph shows maximum range of He ions in selected target(tungsten with some proportions of deuterium) at different energy, calculated by SRIM 2008 code [11]. 8

9 7 Proton yield So far we have developed and describe all the mechanisms and processes that occur in scattering of He nuclei on the nuclei of deuterium. In the experiment a beam of helium ions with energy E 0 is hitting thick layered sample under incident angle α, scattered particles at an angle ψ (see Figure 6). Scattering particles angle is θ = π (α + ψ). If we divide the sample into a series of thin layers and we can ask, how many particles, that are generated in a given layer, enter the detector. The number of detected particles from the selected layer is donated by target yield dy and depends on the cross-section dσi dω of the selected nuclear reaction in the core layer selected, the solid angle Ω of the detector, the particles flux j p and the chemical composition of the individual layers (density, atomic mass and mass elements in layers). Since the treated layer consists of several types of atoms (carbon and deuterium), measured yield of each element in the layers separately[7]: dy i = dσ i(e(e 0, z)) Ωj 0 t ρn A C i (z)adz (10) dω M i where we denoted by C i concentration of selected element in the layers. And i runs through all the elements in the layer at which we are interested in. If you now introduce Q He i = j 0 t ρn A M i Ain considering that dz = dx cos α we can write yield of thin layer as: He dy i = Q i cos α C i(x) dσ i(e(e 0, z)) dω det dx (11) dω wher dω det is the fraction of the solid angle of the detector, which overlooks the particles flying from the sample at an angle ψ. Figure 6: Scattering geometry of the experiment. Incident projectile energy is E 0, the energy just before the nuclear reaction E, energy just after nuclear reaction E, and the output energy from target is E. With α denoted the angle of the incident particle beam, ψ is out going angle and θ is scattering angle of a particle. δω is a small part of the solid angle, which dislodged particles are scattered at an angle θ[7] Total yield per unit of time in selected sample under investigation is computed by integrating expression (11) 9

10 over the entire thickness of the layers and the entire solid angle of the detector. At the end this is resulting in He Y i (E 0, L) = Q i cos α Ω det L 0 C i (X) dσ i(e(e 0, z)) dx (12) dω It is important to point out that in the derivation of the proton yield we assumed to have a homogeneous mixture of deuterium and carrier material (CFC or W) that is without crystal structure. The reason for this lies in the fact that programs which are used for analysis of the measured data, in their Monte Carlo calculation do not take the crystal structure as parameter [11] [12]. 8 Experimental measurements In the experiment, the beam of incident helium ion is focused on the given sample and we are measuring the energy distribution of the proton yield as a function of proton energy. This relationship is obtained by dividing dyp de p (E p ) by energy dependence scatter cross section, we can express from the equation (11)[7]: 1 Y 0 dy D de p (E p ) 1 dσ(e p) dω = C D (E p ) (1) Since we are interested in concentration as a function of position x in the layer and not from the energies E p, we must overwrite (1) in the appropriate dependency. Therefore, we must take into account the proper stopping power to make a connection between the energy E p and the position of the particle in the layer. Since the stopping power is dependent on the concentration of all the elements in the layer, it is not easy problem. Therefore, the concentration profile is determined by an iterative procedure. To be able to deduce the concentration of deuterium from the measured Y p, we must de convolute expression (12).We may do this by a suitable program (i.e. simnra[12]) for simulated concentration profiles, to seek out appropriate energy dependence yields and trying to adjust them to the measured points. Figure 7: The energy distribution yields dy de as a function of the kinetic energy of protons detected at different energies of incident He ions for d = 10 µm thick layer of carbon containing 7% of deuterium, which is evenly distributed over the layers[7]. 10

11 9 Some measurements and results from Jožef Stefan Institute Measurements are executed in a vacuum chamber in one of the experimental lines in the Jožef Stefan Institute tandem accelerator. Measurements are carried out as follows: on each sample, 5 measurements at different energies of the He beam are executed: 4.2 MeV,. MeV, 2.5 MeV, 1.5 MeV and 0.74 MeV. We measure the energy distribution of protons emitted at each of 5 incident energies. In this case we obtain spectrum seen on figure 8 for one energy. Figure 8: The spectrum of fast protons, obtained by He beam at energy of 1.5 Mev, measured at Jožef Stefan Institute. The resulting spectra are then processed with the program simnra. The calculated yields of a model build target are compared with the results of measurements and a depth distribution of deuterium is determined. The result is seen on figure 9. 11

12 Figure 9: The result of the depth distribution of deuterium in a tungsten sample, measured at Jožef Stefan Institute. 10 Upgrade of classical NRA at Jožef Stefan Institute In the laboratories of the Jožef Stefan Institute (JSI) beside of common broadbeem NRA, a so-called µnra is developed. In this method, the beam of He ions is focused through quadrupole lenses on the approximate size of 4 4µm 2, which at energy of 4.5MeV penetrate in to tungsten target approximately 4µm. This enables a good lateral resolution, which provides information on the lateral distribution of deuterium in the walls of fusion reactors. Due to experimental set-up, namely Al foil in front of the detector and measuring at one beam energy we loses all information on the depth distribution of deuterium. During the measurements, in total four detectors are use simultaneously: an NRA PIPS detector, an RBS PIPS detector, an HPGe X-ray detector[1] for -helium-induced X-ray emission (HIXE) and additionally a PIPS detector for beam dose normalization that detects ions scattered from the beam chopper[1], shown on Figure 10. Each detected event in the set of detectors is record and save in a list mode together with information on the beam position at the moment of the event. On the basis of these measurements were obtained the following results shown in Figure 9. Where shown part of elemental inventory of one element of limiter from Tore Supra tokamak. Conclusion In the seminar I tried to demonstrate the use of nuclear reaction D( He, p)α for determining the concentration of deuterium in the sample. At the same time I presented importance of this method for the determination of deuterium deposition in the walls of fusion reactors. Allowing us to better predict the operation of a new generation of fusion reactors, for example ITER. I demonstrated the use of traditional methods and its upgrade mnra. And present some measurement results to Jožef Stefan Institute. 12

13 Figure 10: NRA detection set-up for deuterium analysis with a He focused beam at the JSI microprobe. The three detection systems, NRA PIPS detector, RBS PIPS detector and HPGe X-ray detector for -helium-induced X-ray emission (HIXE) were all positioned at an angle of 15 with respect to the beam direction[14]. Figure 11: Qualitative mapping of the castellation gap of the limiter tile in the erosion zone (Tore Supra DITS experiment, Gap F26T9 QL1) with micro-nra (D) and HIXE (Fe, Cr, Ti, Cu). To cover the total depth of the gap of 6 mm, four frames of µm 2 size were measured, each consecutively shifted 1500µm along the gap depth axis, and merged into a single gap image [14]. References [1]Primož Pelicon, Primož Vavpetič, Zdravko Rupnik, Iztok Čadež, Sabina Markelj, Nataša Grlj, Mirko Ribič, Zvone Grabnar. Influence of the first wall topography on the deuterium retention, 21 st International Conference Nuclear Energy for New Europe, Ljubljana, 2012 [2] Mayer et al. Quantitaaive depth profiling of deuterium up to very large depths, Nucl.Inst. and Meth.in Phys. Rese. B 267 (2009) 506 1

14 []Masashi Shimada, G. Cao, Y. Hatano, T. Oda, Y. Oya, M. Hara and P. Calderoni1. Deuterium Depth Profile in Neutron-Irradiated Tungsten Exposed to Plasma, 2011,1th International Workshop on Plasma-Facing Materials and Components for Fusion Applications [4]Takahiro Izumi, Goroh Itoh. Thermal Desorption Spectroscopy Study on the Hydrogen Trapping States in a Pure Aluminum, Materials Transactions, Vol. 52, No. 2 (2011) pp. 10 to 14 [5]Shane E. Harton, Fred A. Stevie, Dieter P. Griffis, Harald Ade. SIMS depth profiling of deuterium labeled polymers in polymer multilayers, Applied Surface Science 252 (2006) [6] Žiga Šmit. Spektroskopske metode s pospešenimi ioni, 200. DMFA-založništvo [7]Miha Mihovilovič. Določanje globinske porazdelitve devterija z jedrsko reakcijo D( He, p)α, 2008, Univerza v Ljubljani, podiplomski seminar [8] V.Kh. Alimov et al., Nucl. Instr. and Meth. B 24 (2005) 169 [9] W. Möller, F. Besenbacher, Nucl. Instr. and Meth. 168 (1980) 111 [10] T.W. Bonner, Phys. Rev. 88 (1952) 47 [11]J. F. Ziegler, SRIM-The Stopping and Range of Ions in Matter, [12] SIMNRA Home Page by Matej Mayer, mam/, Max-Planck-Institut fã 1 4 r Plasmaphysik, [1] Mitja Kelemen. Spektrometrija mikro-pixe, 2012, Univerza v Ljubljani, seminar [14]P. Pelicon, P. Vavpetič, N. Grlj, I.Čadež,S. Markelj, S. Brezinsek, A. Kreter, T. Dittmar, E. Tsitrone, B. Pégourié, P. Languille, M. Rubel, T. Schwarz-Selinger. Fuel retention study in fusion reactor walls by micro-nra deuterium mapping, Nuclear Instruments and Methods in Physics Research B 269 (2011)