Orientation Effect in D + D Reaction Initiated by 20 keV Deuterons at Channeling in Textured Cvd-diamond Target
Aleksandr Bagulya, Oleg Dalkarov, Mikhail Negodaev, Yuriy Pivovarov, Aleksei Rusetskii, Timur Tukhfatullin
PP R E P R I N T
A. V. BAGULYA, O. D. DALKAROV, M. A. NEGODAEV, Yu. L. PIVOVAROV, A. S. RUSETSKII, T. A. TUKHFATULLIN O RIENTATION E FFECT IN d(d,n) Hе R
EACTION I NITIATED BY
20 keV D
EUTERONS AT C HANNELING IN T EXTURED
CVD-D
IAMOND T ARGET L P I — 2 0 1 6
ISSN 2410-4914
РЕПРИНТЫ
ФИАН ИМ . П . Н . Л ЕБЕДЕВА
ISSN 2410-4914
Главный редактор В . И . Ритус , зам . главного редактора А . А . Гиппиус , научный секретарь С . А . Богачев , ответственный секретарь Л . В . Селезнев
Редакционная коллегия : В . С . Бескин , А . А . Горбацевич , О . Д . Далькаров , Е . И . Демихов , И . Г . Зубарев , К . П . Зыбин , А . А . Ионин , Н . Н . Колачевский , Е . Р . Корешева , С . Ф . Лихачев , А . С . Насибов , И . Д . Новиков , В . Н . Очкин , Н . Г . Полухина , В . С . Лебедев , Н . Н . Сибельдин , Д . Р . Хохлов , С . А . Чайковский
И нформац и я
Препринты
ФИАН им . П . Н . Лебедева издаются с г . Издатель : Федеральное государственное бюджетное учреждение науки
Физический институт им . П . Н . Лебедева
Российской академии наук ( ФИАН ) Адрес редакции : Россия , 119991
Москва , Ленинский просп ., 53,
ФИАНТел .: +7 (499) 132-6137, +7 (499) 783-3640;E-mail: [email protected]
Страница сборника « Препринты
ФИАН им . П . Н . Лебедева » в интернете : http://preprints.lebedev.ru/ © Физический институт им . П . Н . Лебедева
РАН , 2016
A.V. Bagulya , O.D. Dalkarov , M.A. Negodaev , Yu.L. Pivovarov , A.S. Rusetskii , T.A.Tukhfatullin P.N. Lebedev Physical Institute, Russian Academy of Sciences, Moscow, Russia National Research Tomsk Polytechnic University, Tomsk, Russia
Orientation Effect in d(d,n) H е Reaction Initiated by 20 keV Deuterons at Channeling in Textured CVD – Diamond Target Corresponding author. E-mail: [email protected] Corresponding author. E-mail: [email protected]
LPI — 2016
ABSTRACT Orientation effect of increasing the enhancement factor of DD reaction in CVD-Diamond was investigated by simulation. It is obtained that the flux peaking effect up to 2.2 times increases the relative enhancement factor for a parallel beam and up to 1.2 times for the deuteron beam with angular divergence equals 3 critical channeling angles. Qualitative agreement with the experiment was obtained.
1. Introduction
The interactions of the deuterium beam with deuterium enriched fixed targets are investigated in [1–4] using HELIS accelerator facility at the P.N. Lebedev Physical Institute of the Russian Academy of Sciences (LPI). In Ref. [5] the authors investigated the neutron yield in the reaction: )MeV8.0(He)MeV45.2( ndd (1) using textured CVD-Diamond target and 20 keV deuterium beam from HELIS accelerator which delivers the beam with small angular and energy divergences. In [5] the authors suggested that the observed enhancement of neutron yield is connected both with the screening and channeling effects. To clarify the role of channeling in enhancement of neutron yield in d(d,n) He reaction in CVD-Diamond crystal target, we present here the results of computer simulations. The deuterons trajectories in crystal are simulated using the computer code Basic Channeling with Mathematica™ (BCM-1.0) [6], which allows calculate angular and spatial distribution of channeled particles in a thin crystals, see e.g. [7].
2. Experimental setup and results of experiment
Investigation of DD reaction yield was performed at HELIS accelerator facility. This multi-purpose accelerator facility operates continuous ion beams with currents up to 50 mA and energies up to 50 keV. A schematic diagram of the HELIS setup is shown in Fig. 1. Nuclear DD reaction in the interactions of the deuterium ion beam with deuterium-enriched fixed targets were conduct using the polycrystalline deuterium-enriched of 400 µm thick polycrystalline CVD diamond. The method of target formation was the follows. The film was grown on a 57 mm diameter silicon substrate using the microwave plasma-assisted chemical vapour deposition (MPCVD) system [8].
Fig. 1 . Schematic diagram of the HELIS setup: 1 – the ion source (duoplasmatron); 2 – electromagnetic lens; 3 – three-stage chamber of differential pumping; 4 – non contact current meter; 5 – auxiliary ion source; 6 and 10 – Faraday cages; 7 – targets chamber; 8 – the device for calorimetric measurement of the ion beam current; 9 – electrostatic analyzer; 11 – water cooled target holder; 12 – target; 13 – feeder of gas in vacuum chamber; 14 – substrate; 15 – substrate heater.
The black diamond film was obtained, with numerous structural defects in the crystallites, such as twins and amorphous carbon inclusions with a size up to 1 nm. Further the diamond film was separated from the substrate by etching of silicon in a mixture of hydrofluoric nitric and acetic acids, and cut by a Nd: YAG laser into discs of 18 mm diameter. The structure of polycrystalline diamond is anisotropic and not homogeneous. The crystallites are growing in the form of columns, perpendicular towards the surface. The transversal size of crystallites increases from ≈
1 µm in the layer close to the substrate to about 50 µm on the growth side. The growth surface of the sample shows a clear crystalline structure with (100) grain orientation. The neutron flux, produced in the DD reaction, was measured in the longitudinal and transverse direction with respect to the beam axis by using a multichannel neutron detector based of He counters. The relative yield of the DD reaction was determined as Y dd = n n /( S I d ), where n n is the neutron flux, S is the irradiated area of the target and I d denotes the ion beam current. The observed neutron yield measured in longitudinal and transverse directions with respect to the ion beam as a function of the angle between the beam direction and the norm to the target plane is shown in Fig. 2. The enhancement of the DD reaction yield with decreasing angle is clearly observed. Fig. 2 . The relative neutron yield obtained with the CVD-diamond sample as a function of the angle β between the beam and the target plane norm, measured in longitudinal ( ■ ) and transverse ( ♦ ) directions with respect to the ion beam. Ion beam with the energy of E d =20 keV and the current of 50 µA [5]. Fig. 3.
Diamond lattice for <100> direction. Small green circles are positions of deuterium atoms, and large red circles are positions of carbon atoms. Blue vertical line are the (400) channeling planes.
3. Simulations results and discussion
In the first calculation we suggest that deuterium atoms are arranged in perfect diamond lattice between the (400) planes (see in Fig. 3.). In this case, channeled deuterium ions move mainly between the channeling plane and passes near the deuterium atoms situated in the target, which should lead to the enhancement of the DD reaction. a b c d
Fig. 4.
Simulated trajectories of 20 keV deuterium ions in (400) C crystal and corresponding flux density: (a), (b) parallel beam, (c), (d) beam with angular spread Δ θ = θ c /5 ( θ c =1.4 o is the critical channeling angle). Simulations of the ions trajectories are carried out using the code “Basic Channeling with Mathematica” BCM–1.0 [6]. This code computes numerical solutions of the classical equations of motion using the model of continuum potential. Results of calculation of trajectories for 20 keV deuterium ions in a C crystal channeled along the (400) plane are depicted in Fig. 4. Because the potential for planar-channeled positive particles is close to harmonic, the trajectories are characterized by a single oscillation wavelength λ and the oscillation period weakly depend on initial transversal coordinate. So the bright flux peaking effect is clearly observed. For the crystal thickness equals to 0.15 µm we observe the 10 maxima on the flux density function for ideal parallel beam (Fig. 4b). Due to these maxima we expect the enhancement of the DD reaction yield. Computer code allows simulating the beam with difference angular spread; thereby we take into account the angular divergence of the real ion beam and the fact that in textured CVD-Diamond the crystallites can be oriented differently with respect to the target surface as well. Thus in the simulation for each entry points of deuterium ion into the crystal the five incidence angle was generated using random real generator with the normal distribution and specific standard deviation Δ θ . Results of simulation of trajectories for 20 keV deuterium ions with angular spread Δ θ = θ c /5 ( θ c =1.4 o is the critical channelling angle) are shown in the Fig. 4c, d. Because of the the oscillation period weakly depends on the initial transvers energy of deuterons the flux maxima on spatial distribution are not sharp as for the case of parallel beam, that should lead to the decrease of the enhancement factor (EF) of reaction. Probability of the DD reaction depends on the impact parameter of the deuterium ions, penetrating trough the crystal with respect to the deuterium atoms situated in the target. To estimate the relative EF of DD reaction it is necessary to estimate the number of trajectories passing close to the center of the channel where target deuterium atoms are located: L nn xx punchL nn xx pch dxzndxfdz dxzndxfdzR ,,)( , (2) where f ch ( x,z ) is the flux density of channeled ions in C crystal, f unch ( x,z ) is the flux density of ions moving in amorphous target, L is the crystal thickness, n is the number of channel, d p is the interplanar distance. Δ x =0.1 Å is the interval inside which we count the number of trajectories. A charged particle traversing an amorphous target is deflected by many small-angle scattering events. Most of this deflection is due to Coulomb scattering from nuclei and therefore is called multiple Coulomb scattering (MS). The MS distribution is Gaussian for small deflection angles. Thus f unch ( x,z ) is defined as a Gaussian function with standard deviation MSunch , where Δ θ is angular spread of Fig. 5 . Relative EF of the DD reaction versus angle of incidence D + ions into the C crystal for different angular spreed of the beam: ■ Δ θ = θ с / ● Δ θ = θ с / ▲ Δ θ = θ с , ♦ Δ θ =2 θ с , ▼ Δ θ =3 θ с ( θ c =1.39 o is the critical channeling angle). the ion beam and Δ θ MS is the angular divergence due to MS. According [9] for diamond crystal and deuterium ions with energy E =20 keV Δ θ MS = o . The calculation result of relative EF of DD reaction as function of incidence angle of D + ion with respect to channeling plane is shown on the Fig. 5. The relative EF increased with decreasing incident angle of deuterium ions with respect to channeling plane because in this case it is increasing the number of ions involved in channeling motion. The similar dependence of neutron yield was observed in the experiment [5]. Flux peaking effect up to 2.2 times increases the EF of DD reaction for beam with angular divergence Δθ = θ с /5, and it is reduces to 1.2 with increasing angular divergence of ion beam to Δ θ =3 θ с .
4. Conclusions
Orientation effect of increasing the EF of DD reaction in CVD-Diamond was investigated by simulation. It is obtained that the flux peaking effect up to 2.2 times increases the relative EF for a parallel beam and up to 1.2 times for the deuteron beam with angular divergence equals 3 critical channeling angles. Qualitative agreement with the experiment [5] was obtained.
Acknowledgements
The theoretical work was supported by National Research Tomsk Polytechnic University grant No VIU-NRiI-23/2016. 0
References [1]
A.V. Bagulya, O.D. Dalkarov, M.A. Negodaev, A.S. Rusetskii, A.P. Chubenko, Bull. Lebedev Phys. Inst. 39 (9) (2012) 247–253. [2]
A.V. Bagulya, O.D. Dalkarov, M.A. Negodaev, A.S. Rusetskii, A.P. Chubenko, Bull. Lebedev Phys. Inst. 39 (12) (2012) 325–329. [3]
A.V. Bagulya, O.D. Dalkarov, M.A. Negodaev, A.S. Rusetskii, A.P. Chubenko, A.L. Shchepetov, Bull. Lebedev Phys. Inst. 40 (10) (2013) 282–284. [4]
A.V. Bagulya, O.D. Dalkarov, M.A. Negodaev, A.S. Rusetskii, A.P. Chubenko, A.L. Shchepetov, Bull. Lebedev Phys. Inst. 40 (11) (2013) 305–309. [5]
A.V. Bagulya, O.D. Dalkarov, M.A. Negodaev, Nucl. Instr. and Meth. B 355 (2015) 340. [6]
O. V. Bogdanov, E. I. Fiks, K. B.Korotchenko, Yu. L.Pivovarov, T. A.Tukhfatullin, J. Phys.: Conf. Ser. 236 (2010) Article number 012029. [7]
Y.Takabayashi, Yu.L.Pivovarov, T.A. Tukhfatullin, Phys. Lett. B 751 (2015) 453. [8]
V.G. Ralchenko, E. Pleuler, F.X. Lu, D.N. Sovyk, A.P. Bolshakov, S.B. Guo, W.Z. Tang, I.V. Gontar, A.A. Khomich, E.V. Zavedeev, V.I. Konov, Diam. Relat. Mater. 23 (2012) 172–177. [9]
K.A. Olive, et al.,
Particle Data Group, Chin. Phys. C 38 (2014) 090001. ормат х П . л . 0,75. Тираж экз . Заказ № Отпечатано с оригинал - макета заказчика в типографии РИИС
ФИАН
Москва , Ленинский проспект Александр
Васильевич
БАГУЛЯ
Олег
Дмитриевич
ДАЛЬКАРОВ
Михаил
Александрович
НЕГОДАЕВ
Юрий
Леонидович
ПИВОВАРОВАлексей
Сергеевич
РУСЕЦКИЙТимур
Ахатович
ТУХВАТУЛЛИН О РИЕНТАЦИОННЫЙ
ЭФФЕКТ В РЕАКЦИИ d(d,n) He ВЫЗВАННЫЙ
КАНАЛИРОВАННИЕМ
ДЕЙТРОНОВ С ЭНЕРГИЕЙ КЭ В В ТЕКСТУРИРОВАННОЙ