The behavior of HfB2 under neutron irradiation

The behavior of HfB2 under neutron irradiation 1 corresp first-author Corresponding author. E-mail address:mnasseri@aeoi.org.ir Corresponding author. E-mail address:mnasseri@aeoi.org.ir mnasseri@aeoi.org.ir School of Plasma Physics and Nuclear Fusion, Institute of Nuclear Science and Technology (AEOI), Tehran, 141551339, Iran School of Plasma Physics and Nuclear Fusion, Institute of Nuclear Science and Technology (AEOI), Tehran, 141551339, Iran Due to the interesting properties of Hafnium diboride (HfB2) as a ceramic, it has drawn considerable attention from several researchers. To understand the radiation properties of HfB2 that may be used in the nuclear industry, the interaction of this composite under neutron irradiation was studied. It is obvious that, because HfB2 has Boron, this composite has the potential to absorb neutrons particularly in the range of thermal energy. The 10 B (n, α) 7 Li interaction is dominant among other interactions that produce alpha particles. The results revealed that alpha particles have a high fluctuation diagram following a saw-tooth spectrum, for thermal neutron energy which has to be carefully analyzed. If HfB2 intended to be used in fusion facilities, for example; in fusion reactors, its interaction with high neutrons (14 MeV) should be studied. The results for this case showed that, in the wide range of alpha energy, the amount of alpha particle is almost constant so this continuum spectrum is almost flat and there is just a small peak at 2.31 MeV that is belong to the very famous interaction ( 10 B (n, α) 7 Li).. HfB2 Geant4 Thermal Neutron High Energy Neutron Alpha particle. 1. Introduction 1 High efficiency materials in radiation environments and safety remarks in the nuclear industry has been the concern of Nuclear Material Scientists and necessitates them to introduce new materials. Each nuclear facility, even each part of it, requires particular materials to operate properly. One of those materials that was recently introduced is a kind of Ultra Advance Ceramic (UAC) called HfB2. HfB2 elements, Hafnium (Hf) and Boron (B), have a high capacity to absorb neutrons [ 1 1 - 3 3 1-3 ]. Neutrons thermal cross sections of stable Hf and B are 14 and 3980 Barns, respectively. Some HfB2 properties are listed in table (1) table (1) . 1 1 Properties 1 1 Value 1 1 Molar mass (gr) 1 1 200.112 1 1 Melting temperature (℃) 1 1 3380 1 1 Density at 25℃ (g/cm 3 ) 1 1 11.2 1 1 Water solubility 1 1 Insoluble 1 1 Crystal structure 1 1 Hexagonal 1 1 Young's modulus (GPa) 1 1 480 1 1 Hardness (GPa) 1 1 28 1 1 Coefficient of thermal expansion(℃ -1 ) 1 1 6.3 × 10 -6 1 1 Heat capacity at 25℃(J mol -1 ℃ -1 ) 1 1 49.5 1 1 Thermal conductivity (W m -1 ℃ -1 ) 1 1 104 These excellent thermal, mechanical, and nuclear properties are the reasons for using HfB2 as a control rod in fission power plant reactors (e.g. Pressurized Water Reactors (PWR)) [ 4 4 ]. Neut ron calculation and measurement in fusion facilities like tokomak are very important [ 5 5 , 6 6 ] and HfB2 also may play an important role as the first wall of tokomaks. Therefore, HfB2 ceramic with all these properties would be a good candidate for such purposes. Despite of these interesting properties, there are some difficulties in using HfB2 under irradiation of neutrons. This issue originates from the nuclear interactions especially 10 B (n, α) 7 Li. For example through 10 B (n, α) 7 Li interaction, the 10 B isotope transforms into an 7 Li isotope plus an alpha particle. There are two modes of interaction as follow: n + 10B➛7Li*+ α, Q= 2.31 MeV + γ (0.48 MeV)(93%) n + 10B ➛7Li+ α, Q=2.78 MeV (7%) Although these interactions are dominant, other interactions like: 11 B (n, α) 8 Li (β - ) → 8 Be * (2α) or 178 Hf (n, α) 175Yb or 180 Hf (n, α) 177Yb, may also happen. The emitted particle which is basically 4He will be absorbed by the material and will cause defects and dislocation in HfB2 which collapses its structure. According to the literature, thermal neutrons produce alpha particles with an average energy of 1.48MeV and aLithium ion with an average energy of 0.83MeV [ 7 7 ]. The interaction of neutron with HfB2 also, produces prompt γ rays which need more precautions when shielding is intended with this compound. Several experiments have been carried out to determine the cross section of the 10B (n, α) 7Li interaction [ 8 8 , 9 9 ]. Some investigators have compared the practical results with databases such as; the Evaluated Nuclear Data Format (ENDF) or the Japanese Evaluated Nuclear Data Library (JENDL), in a limited energy range [ 10 10 , 11 11 ]. In this study the well known Geant4 has been used to simulate the behavior of HfB2 interaction with neutrons in a wide range of neutron energy. The high Precision of Geant4 (HP) model covers the spectrum of 0-20MeV neutron energy. In this range, the Geant4 Nuclear Data Library (G4NDL) matches the ENDF and JENDL. The model considers neutron interaction with matter, including Capture, Elastic, Inelastic { (n,γ), (n,α), (n,2n), (n,nα), (n,np }, and also different models for the ion interactions [ 12 12 ]. Geant4 has been used for several simulation purposes and shows its validity [ 13 13 , 14 14 , 15 15 , 16 16 , and 17 17 ]. 2. Virtual Experiment 1 A rectangular bulk sample with surface area of 10cm×10cm and thicknesses of 2mm, 10mm, 20mm, 50mm and 100mm has been selected for the simulation. The bulk sample was bombarded by a square neutron beam with the same surface area in a fixed distance of 10cm from the surface of the bulk sample, as shown in figure (1) figure (1) . The sample surface was then hit by 10^6 neutrons during each virtual experiment. Various simulations have been carried out with mono energy neutrons of 0.025 eV, 1 eV, 1 keV, 100 keV, 1 MeV, 10 MeV and14 MeV to determine the production alpha particles resulting from10B(n, α)7Li or other (n, α) interactions. Figure (2) Figure (2) depicts all interactions from a 100 mm thickness of HfB2 slab, except 14 MeV neutrons. It was observed that, in the range of 0.025eV –100keV of neutron energies, the amount of alpha particle spectrums had a drastic fluctuation reaching 1.8 MeV. The compound did not produce any more alpha particles. In figure (2) figure (2) , the fluctuation is due to the alpha particle collisions. On the other hand, the fluctuation is due to counting all alpha particles either produced or collided ones. As we know when alpha particles collide with matter it loses its energy in specific amounts which means we have particles in discrete energy. So if we count all these particles then we have such fluctuations in the spectrum. This spectrum shows the alpha particles that exist in the material at any moment of time. 10 MeV energy has been divided into 256 channels; therefore with the discrete energy of alpha particles we must have more or less amounts in each channel. This fluctuation would be seen for high counting especially in the logarithmic scale. If we notice the spectrum for 1MeV or 10MeV neutron energy, we can see that all the alpha particles are not enough so we do not have such fluctuation in those spectrums. It was also recognized that, at 1MeV neutrons, the alpha particles have an energy range of 0- 4.5MeV. 4.5MeV alpha particle energy is not because of B (n, α) Li. These alphas could be created via other inelastic interactions of neutron with matter. At 10 MeVof neutron energy, the amount of produced alpha particles changes to around 4%, which indicates that only 4% of the 10 MeV neutron flux has been involved in the (n, α) interaction to produce alpha particles. This amount then gradually decreases and ends before 8MeV. The general shape of the alpha particles spectrum is the same for different thicknesses of the HfB2 slab. Figure (3) Figure (3) is plotted for 14MeV neutrons interactions with different thicknesses of HfB2 slab. It indicates that thicker slabs produce more alpha particles. Figure (3) Figure (3) also reveals that, the amount of lower energy alpha particles is higher with descending trends towards a 1.8 MeV alpha particle energy. At 1.8 MeV, a flat peak appears, then followed by a gradually continuing slope and it suddenly drops at 7.4 MeV. After that, the alpha particles’ energy ended at 8 MeV. The interaction of neutrons with HfB2 produces some amounts of energy which would be absorbed by the compound. The deposited energy, as part of a single neutron, was analyzed by different thicknesses of HfB2 slab and the results are plotted in figure (4) figure (4) . As reported by some researchers [ 7 7 ], in 10B (n, α)7Li interaction, the sum of the 4He and 7Li kinetic energies (0.83MeV+1.48MeV) is 2.31MeV. The starting point in figure (4) figure (4) shows that the 0.025eV thermal neutrons produce 2.48 MeV or average, which closely agrees with reported value. Departing from the starting point, the deposited energy within the compound reduces and reaches to its minimum value of around 0.2-0.5MeV, followed by a rapid increase thereafter. The increase of deposited energy could be due to the multiplication of elastic interactions that may produce low energy neutrons in conjunction with the repeating interaction of 10B (n, α) 7Li. This phenomenon has been analyzed and depicted separately in figure (5) figure (5) . The figure shows that neutrons’ elastic interactions starts from 10KeV onward and boosts swiftly as the energy of the neutrons becomes higher. 3. Conclusion 1 The behavior of HfB2 interaction with a wide energy range of neutrons was studied. Boron has a major influence on HfB2 performance and its large thermal neutron cross section plays the main role. Although, Boron can absorb more neutrons, the 10B (n, α) 7Li interaction may introduce some defects in HfB2 structure. Based on this concept, the amount of produced alpha particles and deposited energy within the HfB2 compound interacting with neutrons has been simulated. This virtual experiment approach is very helpful in inserting alpha detectors as well as neutron detectors in between the HfB2 composite, which is very challenging and even impossible. Neutron irradiation to the HfB2 composite affects the amount of Boron during this time. So, all computation of defective effects is time dependent for future studies. References M. A. Misheva and G. Kh. Toumbev , RADIAT EFF DEFECTS , 1985 , 90 : 57 - 66 . DOI: 10.1080/00337578508222516 http://doi.org/10.1080/00337578508222516 RADIAT EFF DEFECTS LI Zhefu ; XUE Xiangxin ; LIU Sulan ; LI Yong ; DUAN Peining , Nucl Sci Tech , 2012 . 23  ( 6 ):  344 - 348 . Nucl Sci Tech YIN Xiangbiao , WEI Yuezhou , ZU Jianhua , Nucl Sci Tech , 2013 . 24  ( 4 ): 040203 Nucl Sci Tech Pascale Cheminant , Xavier Deschanels , Loic Boulanger and Alain Thorel , KEY ENG MAT , 1997 , 132 : 643 - 646 . DOI: 10.4028/www.scientific.net/KEM.132-136.643 http://doi.org/10.4028/www.scientific.net/KEM.132-136.643. . KEY ENG MAT JIANG Xiao-Fei , CAO Jing , JIANG Chun-Yu , CAO Hong-Rui , SONG Xian-Ying , YIN Ze-Jie , Nucl Sci Tech , 2014 , 25  ( 4 ). DOI: 10.13538/j.1001-8042/nst.25.040401 http://doi.org/10.13538/j.1001-8042/nst.25.040401 Nucl Sci Tech LI Shiping , XU Xiufeng , CAO Hongrui , YANG Qingwei , YIN Zejie , Nucl Sci Tech , 2013  Vol. 24  ( 4 ). 040402 Nucl Sci Tech D. Simeone , O. Hablot , V. Micalet , P. Bellon and Y. Serruys , J NUCL MATER , 1997 , 246 : 206 - 214 , DOI: 10.1016/S0022-3115(97)00092-5 http://doi.org/10.1016/S0022-3115(97)00092-5 J NUCL MATER Georgios Giorginis , Vitali Khriatchkov , NUCL INSTRUM METH A , 2005 , 538 : 550 - 558 , DOI: 10.1016/j.nima.2004.08.106 http://doi.org/10.1016/j.nima.2004.08.106 NUCL INSTRUM METH A Simeone , X. Deschanels , D. Gosset , J.P. Bonal and E. Berthoumieux , J NUCL MATER , 2001 , 297 : 244 - 250 , DOI: 10.1016/S0022-3115(01)00644-4 http://doi.org/10.1016/S0022-3115(01)00644-4 J NUCL MATER Georgios Giorginis and Vitali Khryachkov , NUCL INSTRUM METH A , 2006 , 562 : 737 - 740 , DOI: 10.1016/j.nima.2006.02.035 http://doi.org/10.1016/j.nima.2006.02.035 NUCL INSTRUM METH A Guohui Zhang , Li’an Guo , Rongtai Cao , Jiaguo Zhang , Jinxiang Chen , APPL RADIAT ISOTOPES , 2008 , 66 : 1427 - 1430 , DOI: 10.1016/j.apradiso.2007.07.035 http://doi.org/10.1016/j.apradiso.2007.07.035 APPL RADIAT ISOTOPES Online , Available: http://geant4.cern.ch/ http://geant4.cern.ch/ XU Wenzhen , LIU Yanfen , TAN Zongquan , XIAO Ran , KONG Wei , YE Bangjiao , Nucl Sci Tech 2013  Vol. 24  ( 4 ): 040401 Nucl Sci Tech QIN Xue , ZHOU Rong , HAN Ji-Feng , YANG Chao-Wen , Nucl Sci Tech 2015 . 26  ( 1 ), DOI: 10.13538/j.1001-8042/nst.26.010501 http://doi.org/10.13538/j.1001-8042/nst.26.010501 Nucl Sci Tech Mohammad Mehdi Nasseri , MA Qing-Li , YIN Ze-Jie , WU Xiao-YI , Nucl Sci Tech 2004 . 15 ( 5 ), 304 - 307 Nucl Sci Tech Mohammad Mehdi Nasseri , MA Qing-Li , YIN Ze-Jie , WU Xiao-YI , Nucl Sci Tech 2004 . 15 ( 5 ), 308 - 311 Nucl Sci Tech Mohammad Mehdi Nasseri , MA Qing-Li , YIN Ze-Jie , Nucl Sci Tech 2004 . 15 ( 5 ), 361 - 364 Nucl Sci Tech This work was supported by the School of Plasma Physics and Nuclear Fusion, Institute of Nuclear Science and Technology (AEOI). The authors would like to thank Dr.Mohammad Mehdi Farhoudi for his useful comments. 10.1007/s41365-016-0032-7 2015-04-12 2015-07-28 2015-08-11 2016-04-06 2016-04-20 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2016 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-016-0032-7 http://dx.doi.org/10.1007/s41365-016-0032-7 . 2016 Nuclear Science and Techniques 02 27