2020年度核データ研究会
日本原子力学会核データ部会主催、
日本原子力学会「シグマ」調査専門委員会、理化学研究所仁科加速器科学研究センター、日本原子力研究開発機構原子力基礎工学研究センター、東京大学原子核科学研究センター、東京大学、KEK和光原子核科学センター共催、理研シンポジウムのサポートを得ています。
COVID-19 状況下での現地開催となります。対策については、COVID19 issue のページに記載がありますので、一読した上でご参加いただけると幸いです。
Nuclear Data Activities in Nishina Center
Nuclear data study for Accelerator Driven System at J-PARC
Nuclear spectroscopy at KISS
Production and Applications of Radioisotopes at RIKEN RI Beam Factory - Search for New Elements through Diagnosis and Therapy of Cancer -
Measurements of production cross sections of medical radioisotopes via charged-particle induced reactions
Development of Radioisotopes Production Method by Accelerator-based Neutron: Activity at Kyushu University 2020
Microscopic mean-field model is one of strong methods for providing and improving fission-related nuclear data.
It needs appropriate effective interaction, but there is no effective interaction designed for fission path.
In order to tackle this problem, we calculate the U-236 potential energy surface with respect to the elongation of a nucleus and the mass asymmetry with existing Skyrme effective interactions.
We report the energy characteristics of potential energy surface and important parts for correcting the fission barrier.
Mass-angle distribution (MAD) measurement of heavy and superheavy element fragmentation reactions is one of the powerful tools for investigating the mechanism of fission and fusion process. MAD shows a strong correlation between mass and angle when the quasi-fission event is dominant. It has characteristic that appears diagonal correlation as long as the quasi-fission event is dominant. This diagonal correlation could not be reproduced in previous our model before the introduction of the parameters.
In this study, we systematically evaluate the uncertainty model parameters contained in our model and clarify those model parameters to reproduce the diagonal correlation that appears in MAD. Using a dynamical model based on the fluctuation diffraction theorem that employs Langevin equations, we calculate the mass angle distribution and mass distribution of the four reaction systems 48Ti + 186W, 34S + 232Th, 48Ti + 208Pb, and 28S + 238U, which are dominated by quasi-fission. We were able to clarify the effects of uncertain model parameters on the mass angle distribution and mass distribution. In addition, we identified the values of model parameters that can reproduce the correlation between mass and angle. As a result, it found that the balance of tangential friction and moment of inertia values is important for the correlation between mass and angle.
It has been shown that fission has multiple modes, characterized by mass asymmetric fission and mass symmetric fission[1]. In neutron-rich heavy element region,it is argued that several fission modes coexist, with a significant change of their yields in accordance with the number of neutrons contained in the fissionig nucleus. A typical example is found in the isotope dependence of fission for fermium isotopes.
For Fm, the dominant mode transitions from the asymmetric splitting for ${}^\text{257}$Fm to the symmetric for ${}^\text{258}$Fm[2].
This transitions was interpreted as due to the lowering of the fission barrier for symmetric fission, and the becoming energy advantage fission path of symmetric fission then asymmetric fission.
It’s important to know of the potential energy surface structure and nuclear’s deformation process to understand fission mechanism in neutron-rich heavy element region[3].
${}^\text{258}$Md, which is the target of this work, is located near the boundary line where the transition from mass asymmetric fission to mass symmetric fission is expected to occur, and in recent years, the Japan Atomic Energy Agency has obtained the world's first fission data. As a result of data analysis, it shown that the mode of mass symmetric fission (superlong-mode) and the mode of asymmetric fission (standard-mode) coexist.
In this work, we compared the calculation using the fluctuation dissipation model(Langevin calculation)[4] with the experimental data, and considered the characteristics of the fission mode shown by the experimental data.
Reference
[1] U. Brosa, S. Grossmann, and A. Muller, Phys. Rep. 197,167 (1990).
[2] D.C. Hoffman et al., Phys. Rev. C, 21, 1980 (637)
[3] Y. Miyamoto et al., Phys. Rev. C, 99, 051601(R) (2019).
[4] S. Tanaka, Y. Aritomo, Y. Miyamoto, K. Hirose, and K. Nishio PRC 100, 064605 (2019).
Through joint research by the Japan Atomic Energy Agency (JAEA) and Kindai University, it has become clear that the yield distribution of fission products (fission fragments) changes significantly depending on the neutrons emitted from the compound nucleus. In the so-called multichance fission (MCF) concept, fission takes place after emitting several neutrons. This revives the shell structure of a nucleus responsible for mass-asymmetric fission, thus change the fission-fragment mass distribution. The effect of MCF is particularly important to treat high energy fissions, such as ADS system which transmute long lived minor actinide nucleus by fission. Until now, the calculation was performed by combining the fission model calculation (Langevin equation) and a statistical model using a code such as GEF [1,2]. However, this method does not introduce neutron emission during the fission process.
In the present work, we have introduced the neutron evaporation during fission process in the Langevin model. For this, a change of potential energy in each neutron evaporation step is treated. Fission fragment mass distribution of ${^{236-238}}$U,${^{238-240}}$Np, and ${^{240-242}}$Pu were calculated in the initial excitation energy range of E*=15-55MeV. The results show that the double-peak structure is maintained even at the highest excitation energies, and successfully reproduced the experimental data taken at the JAEA tandem facility [3-5].
Reference
[1] K, H, Schmidt, B. Jurado, C. Amouroux, and C. Schmitt, Nucl.Data Sheets 131, 107(2016).
[2] S.Tanaka, Y.Aritomo et al., Phys. Rev. C 100, 064605(2019).
[3] R. Leguillon, K. Nishio et al., Phys. Lett. B 761, 125(2016).
[4] K. Hirose, K. Nishio et al., Phys. Rev. Lett. 119, 222501(2017).
[5] M. Vermeulen, K. Nishio et al., Phys. Rev. C, in print.
Nuclear fission plays an essential role in nucleosynthesis by the rapid-neutron-capture process (r-process), which is a cosmic origin of heavy elements beyond iron. For very neutron-rich environments in neutron star mergers, the strong r-process can be achieved, and the nucleosynthesis path goes into the trans-uranium region. In such conditions, fission is important to shape the r-process abundances due to fission recycling, which determined the termination of the r-process in the heavy nuclei region. Besides abundance prediction, fission is also a key role as the main heating source of kilonovae, which are electromagnetic transients of neutron star mergers. A sign of fission heating may have been observed in the light curve of the kilonova associated with the gravitational wave, GW170817. The precise understanding of fission becomes much crucial in the era of gravitational astronomy.
In this study, we calculate the fission-fragment mass distributions of very neutron-rich nuclei, which are important for the nucleosynthesis calculations of the r-process, but experimental nuclear data is not available. We adopt the Langevin method [1], widely adopted in the study of low-energy fission in the past few years. We found that the calculated mass distributions for uranium, of which Z distribution is calculated with UCD (unchanged charge distribution assumption), well reproduce experimental data in JENDL (${}^{232}{\rm U}$ to ${}^{238}{\rm U}$) [2]. We also found that the fission distribution changes from the two peak feature (asymmetric fission) to the one-peak (symmetric fission) as the neutron number increases. The confirmation by future experiments would be desirable for these theoretical predictions to develop a complete theory set of fission distributions applicable to r-process nucleosynthesis simulations.
[1] S. Tanaka, Y. Aritomo, Y. Miyamoto, K. Hirose, and K. Nishio PRC 100, 064605 (2019).
[2] K. Shibata, O. Iwamoto et al.: "JENDL-4.0: A New Library for Nuclear Science and Engineering," J. Nucl. Sci. Technol. 48(1), 1-30 (2011).
Inelastic scattering is useful tool to explore the nuclear structure in the excited states. In particular, the inelastic excitation to the continuum energy states above the particle decay threshold, which is often called breakup reaction, is very important because we can pin down a specific nuclear structure by controlling the exit channels, which are the combination of the emitted fragments.
A typical and good example about the inelastic scattering to the continuum can be seen in the breakup reaction of $^{12}$Be into the $\alpha$ + $^8$He and $^6$He + $^6$He channels. In this experiment, the careful multi-pole decomposition analysis (MDA) was performed, and the MDA analysis elucidates that many resonant states with a sharp width exist in the spins from $J^\pi=$ 0$^+$ to 8$^+$. The 0$^+$ resonances in the $\alpha$ + $^8$He channel appear in a close energy spacing of 0.5 MeV in the lower energy region below $E_{e.x.} \leq 15$ MeV, which is quite consistent to the energy scheme expected from the cluster resonances.
Basic and important quantities in the analysis of the resonance enhancement embedded in the continuum strength are the resonance parameters, such as the resonance energy and the decay width. In determining the resonance parameters, the evaluation of the non-resonant background strength is indispensable because the resonant enhancement, which has the strong energy dependence, are embedded in the non-resonant background contribution with a broad structure. Since the background strength is structure-less and must have the weak energy dependence, the shape of the non-resonant background strength is often assumed by the simple analytic function or evaluated from the simple reaction mechanism, such as the direct breakup without the final state interaction between the decaying fragments.
In this report, we investigate the structure of the non-resonant background continuum, which is generated by the binary breakup reaction, and explore the prescription to evaluate the background contribution by extending the Migdal-Watson formula for the s-wave breakup in the charge neutral system. In the calculation of the strength function for the binary breakup, we employ the complex scaling method (CSM), which is a powerful tool to describe the few-body continuum states. We handle the breakup reaction of $^{20}$Ne into $\alpha$ + $^{16}$O and $^{12}$Be into $\alpha$ + $^8$He. From the CSM calculation and the Migdal-Watoson theory, we propose the analytic function, which is appropriate to evaluate the background contribution for the binary breakup.
Cobalt (Co) is one of the structural materials in nuclear and accelerator facilities. It is contained in carbon steel and concrete as well as SUS304. $^{59}$Co is only stable isotope of Cobalt. The nuclear data of $^{59}$Co are considered to be important specifically for radioactivity estimation of $^{58,60}$Co related to decommissioning of the facilities. JENDL-4.0 includes the nuclear data of $^{59}$Co, which based the evaluation in 1988. Major revision was carried out at the JENDL-3.3 evaluation in 2001, followed by the covariance estimation in 2002. After the release of JENDL-3.3, many measured data for capture, (n,2n), (n,p), and (n,$\alpha$) reactions have been published. Therefore, the reconsideration of nuclear data is required for JENDL-5.
The evaluation of $^{59}$Co was divided into three energy regions: resolved resonance region, unresolved resonance region, and fast neutron energy region. In the resolved resonance region, the resonance parameters and scattering radius were taken from de Saussure et al. (1992). In the unresolved resonance region, the data of thick sample of de Saussure et al. were adopted, supplemented with the data of thin sample for large resonances. In the fast neutron energy region, the nuclear reaction model code CCONE was used to calculate cross sections, angular distributions and double differential cross sections. The evaluation was performed based on many types of measured data. The obtained results are in good agreement with the measured data and will be shown in the poster presentation.
Oak Ridge National Laboratory released the SCALE6.2 code [1] in 2016 (the latest version is SCALE6.2.4). The ORIGEN code [1] in SCALE6.2 is completely different from the ORIGEN-S code [2] until SCALE6.0 [2].
1) ORIGEN uses one group cross section data generated from a specified neutron spectrum and a multigroup activation library with the COUPLE code [1], not three group cross section data with a typical neutron spectrum.
2) The input format of ORIGEN is easy to use and understand.
3) It is expected that the calculation accuracy improves because ORIGEN uses one group cross section data generated from neutron spectra in all calculation points.
4) The calculation time of ORIGEN including COUPLE is at most about twice of that of ORIGEN-S even for 200 groups.
We expect that ORIGEN in SCALE6.2 will be mainly used for activation calculations in nuclear facility decommissioning. Thus we produced a SCALE6.2 ORIGEN library from JENDL Activation Cross Section File for Nuclear Decommissioning 2017 (JENDL/AD-2017) [3] with the AMPX-6 [4] in order to popularize JENDL/AD-2017 widely. The processing conditions are as follows.
- Temperature : 300 K
- Group structure : 200 groups (the same as one of libraries attached in SCALE6.2)
- Weight function : Maxwell+1/E+Fission spectrum + 1/E (above 10 MeV)
- Infinite dilution
We tested the SCALE6.2 ORIGEN library of JENDL/AD-2017 with the JPDR decommissioning data [5], which demonstrated the library had no problem.
References
[1] ORNL, “SCALE: A Modular Code System for Performing Standardized Computer Analyses for Licensing Evaluation,” ORNL/TM-2005/39 Version 6, Oak Ridge National Laboratory (2009).
[2] (Ed.) W.A. Wieselquist, R.A. Lefebvre, M.A. Jessee, “SCALE Code System,” ORNL/TM-2005/39 Version 6.2.4, Oak Ridge National Laboratory (2020).
[3] https://wwwndc.jaea.go.jp/ftpnd/jendl/jendl-ad-2017.html
[4] D. Wiarda, M.E. Dunn, N.M. Greene, M.L. Williams, C. Celik, L.M. Petrie, “AMPX-6: A Modular Code System for Processing ENDF/B,” ORNL/TM-2016/43, Oak Ridge National Laboratory (2016).
[5] T. Sukegawa, N. Sasamoto, K. Fujiki, “ACCURACY VERIFICATION FOR CALCULATION OF INVENTORY IN JPDR DUE TO NEUTRON ACTIVATION,” INDC(JPN)-164, International Atomic Energy Agency (1993).
To predict the operating lifetime of materials in high-energy radiation environments at accelerator facilities, Monte Carlo codes such as PHITS, MARS, and FLUKA are used to calculate the number of displacements per atom (dpa) related to the number of Frenkel pairs. The Norgertt–Robinson–Torrens (NRT) model has been widely used to predict the number of “initial” Frenkel pairs (NRT-dpa). For more accurate estimation of the actual damage production, athermal-recombination-corrected displacement damage (arc-dpa) was proposed, recently. For the validation of codes, it is necessary to measure displacement cross-sections of metals in relation to changes in electrical resistivity at cryogenic temperature (around 4 K) where the recombination of Frenkel pairs by thermal motion is well suppressed. The comparison between the experimental data and calculated results for proton irradiation with energies from 0.1 to 30 GeV indicates that the arc-dpa results are good agreements with the experimental data and the NRT-dpa results are larger than the data by a factor of around three.
In this presentation, we introduce our experimental plan for displacement cross sections with 120-GeV protons at Fermilab Test Beam Facility (FTBF) in Fermi National Accelerator Laboratory (FNAL). Experiments will be performed at the M03 beam line high rate tracking area in FTBF for the US fiscal year 2022 (October 2021 – September 2022). For the preparation of experiments, we developed the sample assembly with four wire sample of Al, Cu, Nb and W with 250-$\mu$m diameter and 4-cm length. Annealing under vacuum condition (~10$^{-4}$ Pa) was performed by heating an Al sample at 840 K for 30 minutes, a Cu sample at 1289 K for 30 minutes, a Nb sample at 1923 K for 15 minutes, and a W sample at 2473 K for 15 minutes, respectively. The sample assembly will be maintained at around 4 K by using a Gifford–McMahon (GM) cryocooler in a vacuum chamber. Then, changes in the electrical resistivity of samples will be obtained under 120-GeV proton irradiation. Recovery of the accumulated defects through isochronal annealing, which is related to the defect concentration in the sample, will also be measured after the cryogenic irradiation.
This work was supported by JSPS KAKENHI Grant Number JP19H02652.
We require reliable nuclear data that can appropriately evaluate the radiation characteristics of fuel debris for the purpose such as development of new sensors, non-destructive assay technologies and optimization of radiation shielding. In the past, even if different results were obtained depending on calculation codes, it was difficult to clarify what caused the differences. To overcome it, we have developed a new reliable code to calculate radiation decay and radioactive source spectra that can accurately treats with large amounts of nuclides and all decay modes in the decay data file.
As the first step, we compared the photon spectra of fuel debris by using the recent decay data files: JENDL/DDF-2015, decay sub-libraries of ENDF/B-VIII.0 and JEFF-3.3. As shown in Fig.1, the result of JENDL/DDF-2015 is smaller than those of ENDF/B-VIII.0 and JEFF-3.3. This is mainly caused by the following reasons:
(1) X-ray data of 137mBa (T1/2 =2.6 min.) in JENDL/DDF-2015 is missing. The 137mBa is generated from β- decay from large amount of 137Cs (T1/2 = 30 years) and it will remain for a long time by radiation equilibrium.
(2) Gamma ray data of 241Am in 60 keV is missing in JENDL/DDF-2015.
(3) Gamma ray data of 106Rh (T1/2 = 2 hour) is missing in JENDL/DDF-2015 in the energy range from 3.0 to 3.4 MeV. The 106Rh is in the radiation equilibrium with 106Ru (T1/2=1.0 year)
In the presentation, we will report requests for the modifications on the decay schemes and branching ratios of decay modes for the next JENDL decay data file.
For accurate prediction of neutronic characteristics for accelerator-driven system (ADS) and a source term of spallation neutrons for reactor physics experiments for the ADS at Kyoto University Critical Assembly (KUCA), we have launched an experimental program to measure nuclear data on ADS using the Fixed Field Alternating Gradient (FFAG) accelerator at Kyoto University (Period: October 2019 – March 2023). This program is composed of two subprograms, focusing on two nuclear reaction mechanisms, (1) spallation reactions and (2) high-energy fission, for incident proton energies from several tens of MeV to 100 MeV. In the first subprogram, we will measure neutron energy spectra of double-differential cross-sections (DDXs) and thick-target neutron-yields (TTNYs) for several targets (i.e. Pb, Bi, Fe, etc.); in the second subprogram, fission fragment mass number induced from heavy targets (i.e. Pb, Bi) will be measured. In this poster session, the present status of the first subprogram will be presented.
The Accurate Neutron-Nucleus Reaction Measurement Instrument (ANNRI) beamline in the Materials and Life Science (MLF) experimental facility of the Japan Proton Accelerator Research Complex (J-PARC) provides the most intense neutron beam available in the world and was carefully designed to precisely measure neutron-induced reactions using the time-of-flight (TOF) method. Currently, the J-PARC accelerator is operated in double-bunch mode in which two 0.1 μs wide proton
bunches are shot into a spallation target with a time difference of 0.6 μs. Because of this, events detected with a specific time-of-flight (TOF) have two different energies as they could have been originated from each of the two different proton pulses. This is particularly important in the continuum
region (keV region) where the cross section can be expressed as a smooth averaged function. In this region, it is impossible to separate the contribution from each proton pulse and, hence, this mode introduces serious ambiguities into the cross-section measurements.
A neutron filtering system has been designed in order to bypass the double-bunched structure of the neutron beam as part of the “Study on accuracy improvement of fast-neutron capture reaction data of long-lived MAs for development of nuclear transmutation systems” project. Filter materials were introduced into the ANNRI beamline in order to produce quasi-monoenergetic neutron filtered beams. The materials suitable to be used as filters present sharp minima in the total cross-section due to the interference between the potential and s-wave resonance scattering. Neutrons having that energy can be transmitted through the filters and, therefore, produce a quasi-monoenergetic beam. Filter assemblies consisting of Fe with a thickness of 20 cm, and Si with thicknesses of 20 cm and 30 cm of Si were used separately to produce filtered neutron peaks with energies of 24 keV (Fe) and of 54 and 144 keV (Si).
In this study, the characteristics and performance of the neutron filtering system at ANNRI using Fe and Si determined from both measurements and simulations are presented. The incident neutron flux was analyzed by means of transmission experiments using Li-glass detectors and capture experiments using a boron sample which was measured with a NaI(Tl) spectrometer. Moreover, simulations using the PHITS code were performed in order to determine the energy distribution of the integrated filtered peaks and assess the reliability of experimental results. Finally, preliminary results of the capture cross section of $^{197}$Au are presented using the NaI(Tl) spectrometer alongside the neutron filtering system.
Precise nuclear data for neutron-induced reactions are necessary for the design of nuclear transmutation system. Nevertheless, current uncertainties of nuclear data for minor actinide (MA) does not achieve requirements for the design of transmutation facilities. Measurements of the neutron capture cross section are ongoing at the Accurate Neutron Nucleus Reaction measurement Instrument (ANNRI) in the Materials and Life science experiment Facility (MLF) of the Japan Proton Accelerator Research Complex (J-PARC). The determination of an incident neutron flux for measurements of neutron capture cross section is one of the main causes that affect the final uncertainty of the cross section results.
In the present work, we suggest a new method to reduce systematic uncertainties of capture cross section measurements. The method employs change of the self-shielding effect with sample rotation angle. In the new technique, a sample area density of a boron sample which is used for measurements of the incident neutron spectrum. In capture cross section measurements in ANNRI, a boron sample is placed to determine the incident neutron spectrum by counting 478 keV ${\gamma}$-ray from the ${{}^{10}}$B(n,${\alpha}{\gamma}$)${}^{7}$Li reaction. The uncertainty of the boron sample area density that is usually calculated from the mass and the area introduces the uncertainty of the incident neutron spectrum. In this method, the boron sample is tilted with respect to the neutron beam direction, thereby changing the effective area. The neutron self-shielding effect increases with the effective area density. This results in change of the shapes of time-of-flight(TOF) spectrum of 478 keV ${\gamma}$-ray counts form the ${{}^{10}}$B(n,${\alpha}{\gamma}$)${}^{7}$Li reaction with the tilted angle. Comparing the difference of the TOF spectra at different angles and assuming the 1/v energy dependence of cross section of the ${{}^{10}}$B(n,${\alpha}{\gamma}$)${}^{7}$Li reaction, the area density of the boron sample can be determined without using the sample mass and area.
Theoretical and experimental studies on the new method are ongoing. Calculation using Monte Carlo simulation code PHITS were carried out to study the feasibility of the present method. Test experiments using a sample rotation system at ANNRI were also performed. Preliminary results will be given in this poster session.
Highly precise neutron nuclear data is required in nuclear transmutation research of long-lived minor actinides (MA) in nuclear waste. In neutron capture cross section measurement, monitoring the number of the incident neutrons is necessary. However, in measurement with J-PARC/ANNRI, direct neutron monitoring system has not been employed. To make measurement with ANNRI more robust, an additional neutron beam monitor is required. Conventional neutron detectors cannot be used as a beam monitor at ANNRI because of two reasons, high counting rate environment and gamma-flash. The neutron flux at ANNRI is one of the highest in the world. Gamma-flash, an intense gamma-ray burst produced when the proton beam pulse bombards the spallation target, can paralyze a detector generally used in nuclear data measurement. In general, a semiconductor detector or an inorganic scintillator, which is adopted for a neutron detector, has relatively longer response time and is unsuitable for beam monitoring at ANNRI.
Therefore, a combination of a thin plastic scintillator and a $^6$LiF foil was selected as a detection system, whose fast response enabled detecting neutrons at a high counting rate. Low gamma ray sensitivity of a thin plastic scintillator allows measuring fast TOF region without count loss or detector paralysis. The geometry of the $^6$LiF foil, the plastic scintillator, and photomultiplier tube (PMT) was designed. The optimal thickness of the $^6$LiF foil was determined with simulation codes, SRIM and PHITS. A $^6$LiF foil was made by vacuum deposition method. A test detector system was built to study the feasibility of the method.
The detector system was tested under the high neutron irradiation condition at J-PARC /ANNRI. A neutron TOF spectrum was successfully measured without significant count loss or detector paralysis. A neutron energy spectrum was driven from difference of TOF spectrum with and without $^6$LiF. The neutron spectrum was compared with a past neutron spectrum and good agreement was obtained. Statistic error was 0.68 % at 6.0 meV even though measurement times in this study was pretty short (~11 min).
BNCT is a promising cancer therapy which kills tumor cells while suppressing exposure dose to normal tissues. Normally, the neutron field of BNCT, which is produced by a nuclear reactor or an accelerator-based neutron source, has an energy distribution spreading within thermal, epi-thermal and fast neutron regions. Because epi-thermal neutrons are generally used for BNCT, we must measure the epi-thermal neutron flux intensity to evaluate the therapeutic effect and patient’s exposure dose. In addition, we also have to measure the fast neutron flux intensity to evaluate the exposure dose that may be harmful to the human body. However, it is quite difficult to measure such intensities directly and accurately because there is no suitable neutron spectrometer and no activation material covering epi-thermal or fast neutrons separately. The objective of this work is hence to develop new detectors to precisely measure the absolute integral flux intensities of epi-thermal and fast neutrons.
An epi-thermal neutron detector we develop controls its sensitivity by using cadmium and polyethylene. A fast neutron detector controls by using cadmium, B4C and polyethylene. The shape of the epi-thermal neutron detector is a cube, each side of which is 5.52 cm covered with a cadmium sheet. However, the epi-thermal neutron detector is a little sensitive to fast neutrons. To clarify the fast neutron contribution, we develop the fast neutron detector. To extract only fast neutrons, the fast neutron detector consists of two sub-detectors, and the fast neutron intensity is estimated by making difference of the two sub-detectors. The shape of one of them is a cube covered with polyethylene with a side of 4.4 cm and that of the other is a cube covered with B4C with a side of 4.6 cm. Moderated neutrons are measured by activation reaction of 71Ga (n, γ) 72Ga of a GaN foil positioned at the center of the detector. Design calculations were carried out by MCNP5.
After fabricating the detectors, in order to test the performance of the epi-thermal and fast neutron detectors, verification experiments were conducted at KUR, Kyoto University and FNL facility, Tohoku University, respectively.
As the result, the epi-thermal neutron flux intensity could be measured with an error of 3.9 % by correcting the high energy neutron contribution with the calculated value. The fast neutron flux intensity could be measured accurately, that is, the experimental and calculated values agree well within the error range.
<span> The elastic scattering reaction cross section data commonly show smaller in backward angles compared to those of forward angles when the energy of the incident neutron is high. However, in high neutron flux field, such as fusion reactor, the back-scattering reaction cross section is becoming not negligible on the calculation result. Until now, there were differences reported between experimental and calculated values of neutron benchmark experiments using a DT neutron source, which focused on back-scattering phenomena like a gap streaming experiment. For this problem, the author’s group developed a benchmark method for large-angle scattering cross sections and has carried out experiments with an iron sample for the last few years. The benchmark method was successfully established based on the activation of Nb foil having a large activation cross section at around 14 MeV.
<span> We are now planning to carry out benchmark experiments for other fusion structural materials such as tungsten, lead, F82H and so on. And in the next step, we aim to consider benchmark experiments for lighter materials like Li, Be, B, C, N and O. In this case, the energy of neutrons generated by backscattering is low. Especially for Li, being one of the most important materials in fusion reactor, back-scattering neutrons cannot be captured by Nb foil due to the high threshold energy of 93Nb(n,2n) reaction.
<span> In this study, to solve this problem, we examined possible nuclides having a low activation reaction threshold energy, so that it can react with low energy neutrons generated by the backscattering of Li, and simultaneously having not too low threshold energy, so that the influence of room-return neutrons can be eliminated properly. The optimization was achieved by calculating and comparing the number of counts for all the possible reactions of all the existing stable nuclides considering appropriate irradiation and measuring times. The activation reaction cross section data were taken from JENDL/AD-2017.
<span> As a result, we have found that 181Ta(n,2n) was the most suitable reaction giving us the largest number of counts in an acceptable short experimental time. Then experiments were carried out to confirm whether 181Ta(n,2n) cross section was consistent with the nuclear data.
A new nondestructive measurement technique has been developed to evaluate the amount of water in concrete. A concrete wall is irradiated with fast neutrons to activate a gold foil set on the concrete. By evaluating in advance the relation of the gold activity and water content by calculations, we can determine the water content in the concrete, the water content of which is not known. In this study, to validate the present technique experiments were performed with concrete samples having different water contents, which were made from only cement and water. It was confirmed from the experiments that water content could be estimated by the present nondestructive measurement technique though the system is still simple with cement and water. Now we are examining the validity for concretes made from cement, water and sand.
The double differential cross-section (DDX) of the photoneutron is an important quantity for radiation shielding and shielding calculation of the electron accelerator design. Shielding calculation is usually carried out by Monte Carlo simulation codes, which use the nuclear data library in the calculation. We measured the DDXs of the $\rm(\gamma,xn)$ reaction using 16.6 MeV polarized photon on $\rm^{nat.}Pb$, $\rm^{197}Au$, $\rm^{nat.}Sn$, $\rm^{nat.}Cu$, $\rm^{nat.}Fe$, and $\rm^{nat.}Ti$ targets [1], and showed the neutron spectra including the evaporation and direct components. In this presentation, we compared the DDXs from the photonuclear data JENDL/PD-2016.1 and the experiment to check their consistency. The DDXs were extracted from the JENDL/PD-2016.1 library by our python-based software. The abundances of each target's isotopes were considered in calculating the DDXs from the JENDL/PD-2016.1 library. This comparison showed the differences in the photoneutron’s energy distributions, and the differences were mostly in the direct component. The evaporation components were found in both JENDL/PD-2016.1 and the experimental data; however, they were not totally consistent. This result can be the first comparison of the DDXs from JENDL/PD-2016.1 and the experiment. The quantitative comparison and discussion will be presented at the symposium.
Keywords: differential double cross-section photoneutron, 16.6 MeV polarized photon, JENDL/PS-2016.1.
[1] T.K. Tuyet et al., “Double differential cross section of the $\rm(\gamma,xn)$ reaction on medium-heavy nuclei for 16.6 MeV polarized photons. 2020 Fall meeting of AESJ, Sep. 16th – 18th.
To optimize disposal of low-level radioactive waste originating from decommissioning of nuclear facilities, required are 1) reliable assessment of radioactivity level by calculation and measurement and 2) a good estimate of the uncertainty of those results for the classification of radioactive waste. In order to improve the reliability of the calculations in clearance verification, we established a procedure of estimating the uncertainty of radioactivity concentration due to that of nuclear data. For that, we estimated covariance of neutron cross sections of important nuclides that account for over 90 % in ΣD/C of concrete material and carbon steel by employing a propagation of uncertainties in the resonance parameters and statistical model parameters with nuclear data code group T6. Here, D stands for radioactivity concentration, and C stands for clearance level. Then, we developed a new method to calculate uncertainty of radioactivity with Total Monte Carlo method by connecting randomly perturbed endf-format files generated by the T6 calculation to a cross section processing code NJOY and an activation calculation code ORIGEN2 using ORLIBJ40, a set of cross section library based on JENDL-4.0. It was concluded that the uncertainty of the radioactivity due to that of nuclear data for nuclides which dominate the ΣD/C is sufficiently small, and the main factor of uncertainty of radioactivity comes from that of the neutron flux.
RIKEN Accelerator driven compact Neutron Source-II (RANS-II) based on the $^7$Li(p, n)$^7$Be reaction for neutron production with 2.49 MeV proton beam, has been under beam commissioning to demonstrate specific performance of the system. RIKEN has a prospect of realizing novel non-destructive neutron inspection for infrastructures with the use of RANS. As prominent characteristics, RANS-II has the maximum neutron energy of 0.8 MeV, which is lower than that of 5 MeV at RANS based on the $^9$Be(p, n)$^9$B reaction with 7 MeV proton injection, and gives extremely forward favored angular distribution with respect to the proton beam direction. Also, it should be emphasized that RANS-II system is installed in a relatively small space isolated by concrete shield with boron containment. Accordingly, there should be quite large differences in neutronic performances between RANS-II and RANS in terms of neutron spectrum and angular distributions. In preparation of experiments at RANS-II, the simulation of radiation fields for neutron and $\gamma$-ray in RANS-II experimental hall plays a critical important role for designing experimental set-up in low background.
Then, we have performed simulations to characterize radiation fields of RANS-II The cross section libraries implemented in PHITS are utilized in neutron and $\gamma$-ray transportations. Several important conditions of RANS-II modeling are as follows:
・ The lithium (Li) target is made by depositing thin Li layer of about 100 $\mu$m on a 5 mm thick Cu substrate cooled by water in the target station.
・ The target station with about 90 cm side cubic shape, configures five layers; polyethylene, lead, borated polyethylene, lead and iron, to reduce the radiation leakage.
・ There is a hole with a 15 $\times$ 15 cm$^2$ cross section in the forward direction.
・ The experimental hall has dimensions of 14 $\times$ 5.5 $\times$ 3.0 m$^3$ surrounded by floor and wall made of concrete (partly the borated concrete) and polyethylene ceiling.
As a result, the neutron and γ-ray distribution spreads widely in the experimental hall due to the wide openings in the target station. The scattered radiation could be the major contributor to the background of experiments. On the other hand, primary neutrons produced at the target are shielded reasonably by the target station. To design effective collimators for high quality beam extraction, we have calculated neutron beam profiles with parameters of collimator diameters and materials. It is shown that there is an optimized collimator configuration to extract suitable beam.
Reliable assessment of radioactivity in target and structural materials for high-energy accelerator facilities such as accelerator-driven systems and spallation neutron sources requires detailed information on nuclide production cross sections by spallation reactions. To obtain the systematic cross section data for nuclide productions by spallation reactions, we have conducted irradiation experiments at Japan Proton Accelerator Research Complex (J-PARC). So far, we have measured nuclide production cross sections for light to medium-heavy target elements (Z≤47) with proton energies ranging from 0.4 to 3.0 GeV. To investigate heavier target elements, we conducted an experiment for target elements with atomic number around Z=70, including $^{\rm nat}{\rm Lu}$ (Z=71) target.
Four sets of Ho ($90 \rm mg/cm^2$), Lu ($100 \rm mg/cm^2$), and Re ($210 \rm mg/cm^2$) foils were packed in aluminum containers together with 0.1-mm-thick aluminum catchers to avoid recoil contamination. Each set of targets was irradiated with 0.4-, 1.3-, 2.2-, and 3.0-GeV protons accelerated by 3-GeV Rapid Cycling Synchrotron (RCS). The beam current was monitored by a current transformer installed in front of the irradiation position. After the irradiation, gamma-rays emitted from the samples were detected by two high-purity Germanium detectors (relative efficiency 20%, Canberra Co., Ltd.).
The measured cross sections were compared with theoretical predictions by Particle and Heavy Ion Transport code System (PHITS) [1], and INCL++/ABLA[2,3]. The figure shows experimental and calculated $^{\rm nat}{\rm Lu}(p,X)^{\rm nat}{\rm Be}$ reaction cross sections. INCL/GEM model implemented in PHITS underestimated the experimental cross sections by a factor of about 2.
In the presentation, we will report our experimental results for the natLu target, and more detailed discussion on reaction mechanics will be given.
References
[1] T. Sato, et al., J. Nucl. Sci. Technol. 55, pp.684 (2018).
[2] D. Mancusi et al., Phys. Rev. C 91:034602 (2015).
[3] A. Kelic et al., Proceedings of Joint ICTP-IAEA Advanced Workshop on Model Codes for Spallation Reactions, ICTP Trieste, Italy pp.181 (2008).
COherent Muon to Electron Transition (COMET) is an experiment at J-PARC, which will search for coherent neutrino-less conversion of a muon to an electron in muonic atoms. The experiment will be carried out in two steps: Phase-I and Phase-II. In the Phase-I experiment, 3.2 kW 8 GeV proton beam irradiates a 70-cm long graphite target to produce negative pions. The negative pions are captured in the magnetic field and delivered to pion-decay and muon-transport sections. The Phase-I experiment aims to detect the $\mu^-e$ conversion events and measure the beam-related background events for the Phase-II experiment. Now, it has been planned that the maintenance by radiation workers would be conducted after the 150-day operation and the following 180-day cooling. It is necessary to evaluate the residual radiation dose for the safety of the workers during the maintenance. On this study, we calculated fluxes of neutron, photon, proton and other charged particles in the beam room during the beam operation and the residual activity after the cooling time by using Monte Carlo simulation code PHITS version 3.22 and DCHAIN-PHITS version 3.21. The calculation results show that the design of components around the target and beam dump needs to be improved to reduce the radioactivities after the cooling time.
Recently, a small modular reactor (SMR) with inherent and passive safety has been receiving attention all over the world. In Japan, a very small modular reactor, namely, MoveluXTM(Mobile-Very-Small reactor for Local Utility in X-mark) has been developing by Toshiba Energy Systems＆Solutions Corporation. MoveluXTM is a thermal reactor that uses a calcium hydride as a neutron moderator. The use of a Sn-Pb alloy as an in-core heat transport medium is being considered. The Sn-Pb alloy is in a solid state when the reactor is started, and becomes liquid since the core temperature reaches 660℃ during operation. Therefore, the total cross section data of the Sn-Pb alloy is important for evaluating the effect of the change in the total cross section depending on the state of Sn-Pb alloy on the reactor characteristics. However, there are no reports on experimental data for total cross section of Sn-Pb alloys in both solid and liquid states in spite of the fact that it is important data for nuclear engineering. In the present study, the neutron total cross section was obtained from neutron transmission measurements by the time-of-flight (TOF) method using the Kyoto University Institute for Institute for Integrated Radiation and Nuclear Science – Linear Accelerator (KURNS-LINAC). The sample temperature was changed from room temperature (solid) up to 300℃ (liquid). The total cross sections of solid and liquid states were compared and the change in Bragg edge due to the difference of crystal structure was observed in the energy region below 0.01 eV. Comparing the total cross sections of the solid and the solid resolidificated after melting, it was confirmed that some Bragg edges, which are thought to be due to the crystal structure of Pb, disappeared by the resolidification. At the poster presentation, the detail of the total cross section measurement experiment and the results obtained so far will be discussed.
Processing of spent fuel from nuclear power plants is a worldwide problem. The high-level radioactive waste is the product after the reprocessing of spent fuel, which includes minor actinides and fission products of radioactive waste. Especially, $^{90}$Sr (T$_{1/2}$ = 28.8 years) is the highest radiotoxic nuclide in the fission products. It is highly desired to develop nuclear transmutation technology using accelerator facilities to reduce these harmful nuclides. The simplest way can be to irradiate a neutron beam on the radioactive waste. However, it is not well known that $^{90}$Sr is transmuted into how much and which nuclide in this reaction. Therefore, it is essential to study, in advance, the reaction-cross-sections to each nuclide from $^{90}$Sr. From this point of view, the inverse kinematics, i.e. including the $^{90}$Sr beam incident on light-particle targets, is an effective method the reaction products can be identified at the forward directions.
To realize this purpose, we have planned the proton- and deuteron-induced reaction-cross-section measurements in inverse kinematics and performed the experiment using the BigRIPS separator [1] and the ZeroDegree spectrometer [1] at the RIKEN Radioactive Isotope Beam Factory. The radioactive $^{90}$Sr beam with 104 MeV/u, produced and separated in the BigRIPS, incident on the C, CH$_{2}$, and CD$_{2}$ targets. The reaction products in the forward directions were transferred to the ZeroDegree and identified using the detectors at the focal plane. The reaction-cross-sections were obtained from the measured yields of each reaction channel. At this time, the contributions from carbon and beam-line materials were subtracted as a background. The obtained reaction-cross-sections were compared to the PHITS calculation [2] and the data with different energy of 185 MeV/u [3].
[1] T. Kubo, et al., Progr. Theor. Exp. Phys. 2012, 03C003 (2012).
[2] T. Sato, et al., J. Nucl. Sci. Technol. 50, 913 (2013).
[3] H. Wang, et al., Phys. Lett. B 754, 104 (2016).
A long-lived isotope of Hf, 175Hf (T1/2 = 70 d), is useful for basic studies for rutherfordium (Rf, Z = 104). This isotope is producible in no-carrier-added form in the proton- and deuteron-induced reactions on natLu. However, excitation functions of these nuclear reactions have been scarcely studied. In this work, we measured the excitation functions of the natLu(p,xn)175Hf and natLu(d,xn)175Hf reactions up to 18-MeV proton and 24-MeV deuteron energies using a stack-foil technique and a γ-ray spectrometry. We performed these experiments at RIKEN and Institute for Nuclear Research (ATOMKI). The target stacks of Ta/Lu/Ti and Lu/Ti foils were irradiated for 2 h with proton or deuteron beams of approximately 180‒240 nA. After the irradiation, each foil was subjected to γ-ray spectrometry with Ge detectors. We noticed that the half-life of 173Hf is slightly longer than that adopted in the current nuclear database. Therefore, we measured a precision half-life of 173Hf in a separate experiment. In this work, we could measure the excitation functions of the natLu(p,xn)173,175Hf and natLu(d,x)173,175Hf, 173,174m,174g,176m,177m,177gLu reactions. Thick-target yields of 175Hf were also deduced from the measured excitation functions. The yields are 0.47 MBq/µA·h at 17.2-MeV proton and 2.0 MBq/µA·h at 24.0 MeV deuteron. We determined the half-life of 173Hf to be 24.176 ± 0.012 h which is 0.58 ± 0.10 h longer than that in the database.
Quantification of radioactivity of fission products (FP) is very important for assessment of decay heat after shutdown of a core, etc. For such assessments, comprehensive data sets of fission yield and decay chain, such as JENDL/FPY&FPD-2011, have been developed. However, validation of each nuclide in such data sets has still been cumbersome. In this work, two detection techniques of FPs are studied to give data for such validation.
In order to characterize reactions occurred in nuclear fuel, gamma ray spectroscopy was conducted at Kyoto university critical facility assembly (KUCA). At KUCA, uranium (U)fuel of 93 wt% 235U enrichment was loaded in C-core. They were moderated and shielded by light water. The core power during the critical operation was 4.6 mW. Outside the tank of the core, a HP-Ge detector of 30 % relative efficiency was set and the gamma ray was measured. As the results, peak spectra of fission products such as 90, 95, 97Y, 90,90m,91Rb, 87,88Br, 136Te, etc. were detected although they were overwrapped by prompt gamma ray components. Due to the prompt components, the relative statistical accuracy was from 2 to 20 %. Thanks to the measurements during the critical operation, gamma rays of half-life shorter than 4 s was achieved.
Contrarily, 238U(n,g) gamma ray spectroscopy was conducted with the same HP-Ge for neutrons of thermal and resonance energy at the KU-LINAC-pulsed neutron source facility. The time of flight (TOF) of neutron was measured associated with beam pulse to identify the incident neutron energy. The repetition rate of the pulse was 50 Hz. In the TOF spectrum after the so called “thermal neutron peak”, time-background region was identified. The gamma ray in the region out of phase of the beam pulse was considered emitted by decay of radioactive material of which half-life is longer than 20 ms. The measured peak structure was found fairly resemble to that of measured at KUCA.
The detection efficiency of the gamma rays at KUCA was calculated with MCNP-5. That at LINAC was experimentally determined. With the measured count rate and the efficiency, the gamma ray emission rate was deduced and compared against that calculated based on JENDL/FPY&FPD-2011. The ratio of the measured to the calculated value against each gamma ray by the two experiments show fairly resemble trend. That indicates the both experiments are promising to give reference data for validation of FP yield and decay data sets such as JENDL/FPY&FPD-2011.
Isotope production in spallation reaction of 93Zr and 93Nb induced by proton and deuteron
Development of energy-degraded RI beam and expansion of nuclear reaction studies
Theoretical analysis of deuteron-induced reactions and development of deuteron nuclear database
Roles and current status of reactor physics experiment in research reactors
Deep Learning for Basic Science
Data-driven approaches for nuclear shell-model calculations
Nuclear data generation using machine learning
Exploration of automated data processing for mass production of nuclear data at RIBF