Estimation of double-differential cross-sections of 9 Be(p,xn) reaction for new nuclear data library JENDL-5

Abstract

Double-differential cross-sections of the ⁹Be(p,xn) reaction are newly evaluated based on the Wakabayashi’s function and neutronics analysis up to 12 MeV for a new nuclear data library, JENDL-5. We devoted our efforts especially to the re-optimization of the absolute cross-sections and the interpolation of the neutron energy spectra. Through the comparisons between the thick target yield measurements and Monte-Carlo simulations at different proton energies and neutron emission angles, we conclude that JENDL-5 gives the best evaluation in the world.

Keywords

Nuclear data Double-differential cross-section Wakabayashi’s function Interpolation Neutronics simulation Thick target yield JENDL-5

1. Introduction

The cross-section of the ⁹Be(p,xn) reaction is highly important for the design of neutron sources with small accelerators. However, its evaluated values in the world’s nuclear data libraries have large uncertainties since available experimental data-sets of the cross-sections are very scarce, and some data does not match each other. It is also because the study of the nuclear theories/modelings, which is often applied to the prediction of cross-sections or interpolation of experimental values, is still on a scientific challenge for the few-body systems, viz., light nuclei such as ⁹Be.

Looking at the status of absolute cross-section measurements which are available in the experimental database EXFOR [15], several data-sets are certainly found below ∼15 MeV. However, most of them were measured before the 70s, in which the data are given in a narrow energy region and discrepant each other by more than $20 %$ . Among the existing data, the experimental data of Bair et al. [2] gives the excitation function in a relatively wider energy range (4–14 MeV) than the others, but the measurement is a so-called relative one which means the absolute information is totally absent. According to the EXFOR database, there are similar situations also for the measured differential cross-sections.

In the nuclear data evaluation, the statistical theory [8] which is combined with nuclear structure models is often used for the interpolation of available experimental data. It is also applied to the prediction/estimation of differential and double-differential cross-sections (DDX). Indeed, our nuclear reaction model code CCONE [12] is known as one of the powerful tools for the nuclear data evaluation, with a minimal optimization of the model parameters therein, as it has been applied to the development of our nuclear data libraries, e.g., JENDL-4.0 [17]. However, it is difficult to apply such a kind of code to the light nuclei ( $A ≲ 20$ ) because they are few-body systems which may be out of the scope of the statistical theory. Also, understanding/modelling the structure of light nuclei is at the stage of a scientific study.

One of our special purpose nuclear data libraries, JENDL-4.0/HE [13], provides the evaluated cross-sections data up to 200 MeV to be used for the design of accelerator-based applications such as the neutron sources, accelerator-driven system for the transmutation of radioactive wastes, and so forth. Indeed, this library is implemented in the Particle and Heavy Ion Transport Code System (PHITS) code [16] for an optional use. This library certainly includes data of ⁹Be for incident protons, but we managed to apply the CCONE code to the data evaluation for ⁹Be despite the inaccuracy of the measurements and nuclear models mentioned above, which means the data for p+⁹Be in JENDL-4.0/HE could have a large uncertainty. Quite recently, one of the authors, Wakabayashi et al. [21], proposed a function which describes the cross section and DDX of ⁹Be(p,xn) where a number of parameters are optimized with selected experimental data-sets. The function is expected to be more suitable to the practical nuclear data evaluation since it is more flexible than the nuclear models.

The objectives of the present study are:

evaluation of the cross section for the ⁹Be(p,xn) reaction based on the Wakabayashi’s function,

generation of DDX with a physically reasonable interpolation scheme for the ⁹Be(p,xn) reaction based on the function

for our new nuclear data library JENDL-5 [11] which intends to provide evaluated nuclear data not only for the nuclear energy but also for the accelerator-based applications.

2. Data evaluation for ⁹Be(p,xn)

2.1. The Wakabayashi’s function

The energy dependence of the absolute cross-section is rather smooth except for the lower energy region ( $E_{n} ≲ 5 MeV$ ) where several resonances are observed in the available measurements. The neutron energy spectrum from the ⁹Be(p,xn) reaction is characterized by a continuum part and peaks which correspond with the transition to the ground and low-lying excited states of the residual nuclei $^{9} B^{(⋆)}$ . These characters are so complicated that it is difficult to predict them consistently only by a combination of nuclear models. Recently, one of the authors, Wakabayashi et al. [21], proposed a function which was designed to describe all the characters up to $E_{n} = 12 MeV$ , in which the parameters therein were optimized by selected experimental data-sets. According to the neutronics analysis with GEANT-4 [1], the function gives a reasonable estimation of thick target yield (TTY) better than those with existing nuclear data libraries. Unfortunately, MCNP [22] or PHITS [16] cannot use the function directly since they are designed to use nuclear data libraries or built-in models (e.g., Ref [4]) for high-energy nuclear reactions.

2.2. Incorporation of the function to nuclear data file

In the present study, we incorporated the Wakabayashi’s function to a nuclear data file, which finally enables us to perform neutronics analyses with both MCNP and PHITS. In a first step, the energy dependent cross-section of the ⁹Be(p,xn) reaction was calculated with the function, taking care of the selection of the energy grids to fully describe the resonant shape in the lower energy region. In the next step, we produced data table of the neutron energy spectra from the function at different proton energies and emission angles. In this case, much attention was paid also to the determination of neutron energy nodes especially around the peak positions to minimize the data size. Such point-wise data-sets were compiled in the ENDF-6 format [20] file.1 Since a normalized angle-energy distribution was required to be compiled for each incident energy in the ENDF-6 format, we normalized DDX with respect to the neutron energy $E_{n}$ and emission cosine angle μ for each proton energy $E_{p}$ as, $\begin{matrix} (1) & \int d E_{n} \int d μ f (μ, E_{p}, E_{n}) = 1 \end{matrix}$ where $f (μ, E_{p}, E_{n})$ is the normalized DDX to be given in the data file.

Fig. 1.

Evaluated and measured cross-sections of the ⁹Be(p,xn) reaction.

In compilation process, an interpolation scheme must be given to the data-tables because the data we generated from the function were forced to be discretized in principle. In the present study, the linear-linear interpolation was assumed for excitation function and also for each neutron energy spectrum at an emission angle. On the other hand, the two-dimensional unit-base [20] was assumed for the interpolation between the neutron spectra of two adjacent incident proton energies to minimize pile-ups of peaks which appear if the normal one-dimensional interpolation is used. Finally, we replaced the ⁹Be(p,xn) data in JENDL-4.0/HE with the present data tables from the threshold energy of the reaction to 12 MeV. Such a data merging is necessary because a complete data set of cross-sections (which includes those for the other types of reactions) is expected in the neutronics simulation.

Figure 1 shows the cross-section obtained from the function together with available experimental data [2,7,19] and evaluations. It is clear that the original function is designed to reproduce the measured data by Gibbons et al. [7] and the normalized data by Bair et al. [2]. Figure 2 illustrates an example of DDX extracted from the present data file with experimental data [3]. The two peaks observed in this spectrum correspond to the ground and excited states of ⁹B, with a full width at half maximum of 0.07 MeV [21]. As shown by one of the curves (INT = 2), it is difficult to produce sharp peaks correctly if we use the normal linear-linear interpolation. This issue encountered with the numerical data-table is solved in the present study with the unit-base interpolation, where we kept the ratio of the peak energy to maximum energy constant in the unit base, for the neutron spectrum interpolation as illustrated by the other curve (INT = 22). It is concluded that the function gives a reasonable estimation of the cross-section and DDX, and the applied interpolation schemes are reasonable.

Fig. 2.

Example of DDXs which are extracted from the present data file with the interpolation. The red curve (INT = 22) shows data with the unit-base interpolation while the blue curve (INT = 2) shows data with the normal linear-linear interpolation. Note that those curves are taken from the final data file (JENDL-5) where the absolute cross section is reduced by $15 %$ from the original function.

2.3. Neutronics analysis

Although the measured cross-sections including DDX are very scarce, there exists experimental data of TTY which were published in the current century, at different incident energies and emission angles, fortunately. Therefore, those measurements are quite useful to validate the present nuclear data file generated as described above. For example, the measured data by Howard et al. [9] are quite useful as they give neutron energy spectra at $E_{p} = 4, 5$ and 6 MeV from the most forward to backward emission angles. Before carrying out the neutronics analysis, the ACE format [6] file was generated for the Monte-Carlo simulation. We used the NJOY2016 code [14], which was modified for JENDL-5, for the data processing. Note that the two-dimensional unit-base interpolation was dealt as the linear-linear interpolation because it was not supported in NJOY and Monte-Carlo codes, unfortunately.

Neutronics simulations were performed by the MCNP code [22] with the ACE files of the present data file (original function), JENDL-4.0/HE [13] and ENDF/B-VIII.0 [5]. Figure 3 illustrates an example of the TTY simulation results together with the experimental data by Howard et. al. [9]. The result obtained with ENDF/B-VIII.0 is not consistent with the shape of the measured data. That with JENDL-4.0/HE seems to give a better estimation, but it largely underestimates the experimental data in the lower energy region. It is obvious that the result (original function) with the present file based on the original Wakabayashi’s function simulates the measured neutron spectrum better than the existing libraries. However, the result with the present data overestimates the measured spectrum by $10 \sim 20 %$ . Such a discrepancy is also systematically confirmed in other proton energies and emission angles.

Fig. 3.

Example of simulation results for the thick target yield for $E_{p}$ = 3 and 4 MeV at emission angle of 0 deg.

2.4. Estimation of cross-sections for JENDL-5

Even though we have confirmed that the function already gave the best estimation of TTY among the data libraries in the world, we noticed that the present data overestimated measured TTY systematically by $10 \sim 20 %$ . As we already mentioned above, the measured data-sets of the absolute cross-sections do not match each other. Besides, the experimental excitation function by Bair et. al. is a relative measurement, which means the absolute value was not known. Therefore, the present cross-section should also have a uncertainty of the same order, which may result in a systematic overestimation of the measured TTY.

For the new nuclear data library JENDL-5, we finally decided to reduce the absolute cross-section by $15 %$ from the original function as illustrated in Fig. 1. Figure 3 demonstrates that the prediction ability of the neutronics simulation with the re-optimized data file (JENDL-5) is much better than that based on the original function. From similar comparisons at other proton energies and neutron emission angles, we conclude that the present data file (JENDL-5) gives the best estimation worldwide.

2.5. JENDL-5 and remarks on unsolved issue of $p +^{9} Be$ evaluation

The new nuclear data library JENDL-5, which was compiled in the ENDF-6 format, was released in December 2021.2 This is the first major update of the previous version JENDL-4.0 which dates from 2010. The features of JENDL-5 relevant to the small accelerator might be summarised as follows:

revision of a large part of JENDL-4.0 such as neutron reaction data, thermal scattering law data

neutron activation cross section data

the first originally evaluated thermal scattering law data for 16 materials such as light and heavy water, benzene, and ethanol

reaction data induced by light charged particles (proton, deuteron and alpha-particle) and photon.

The evaluation for ⁹Be is included in the proton-induced sub-library where the present data are compiled up to 12 MeV. Above 12 MeV, the data are taken from JENDL-4.0/HE, where a discontinuity of the ⁹Be(p,xn) reaction cross section at the boundary energy does not appear between JENDL-5 and JENDL-4.0/HE as illustrated in Fig. 1. This is probably not a coincidence because the characteristics of nucleus tend to fade out as an excitation energy of compound nuclei (incident proton energy) becomes higher.

JENDL-5 application libraries including ACE format files will also be soon available. Touching upon the PHITS code, the current version of code package includes JENDL-4.0 in the ACE format. Similarly, JENDL-4.0/HE is also included for a number of nuclei relevant with accelerator applications. In the future, instead of JENDL-4 and JENDL-4/HE, JENDL-5 for neutrons below 20 MeV and for neutrons and protons up to 200 MeV for specific nuclides will be bundled to the PHITS package.

Although we have obtained the nuclear data evaluation for $p +^{9} Be$ which gives the best estimation of the neutron spectra in the world, there still exists an issue to be solved in the future. Indeed, we always observed overestimations of the measured TTY of Howard et al. at backward angles as illustrated in Fig. 4. The other libraries also give shapes which are different from the experimental data. There should be a fundamental issue on description for the angular distribution both in the Wakabayashi’s function and nuclear models used in the previous libraries. In near future, we plan to update the function with incorporation of the coupled-channel theory [18] only for the direct reaction to overcome such a remaining issue. Also, it may be meaningful to look into a recent progress on the machine-learning approach to the data evaluation [10].

Fig. 4.

Example of simulation results for the thick target yield at the backward angle.

3. Conclusion

We have evaluated the cross section and DDX of the ⁹Be(p,xn) reaction based on the function proposed by Wakabayashi et. al. up to 12 MeV. Through the neutronics analysis, we finally decided to reduce the absolute cross-section by $15 %$ from the original one to provide a better agreement with the TTY measurements. The evaluated cross section and DDX data are compiled in the ENDF-6 format with physically reasonable interpolation schemes. This data is merged in our new nuclear data library JENDL-5 released last December. We conclude that JENDL-5 gives the best prediction worldwide of neutron spectra from a beryllium target bombarded with proton.

Footnotes

We adopted the laboratory angle-energy law (MF = 6, LAW = 7) in the ENDF-6 format to compile the normalized spectra.

https://wwwndc.jaea.go.jp/jendl/j5/j5.html

Acknowledgements

The authors would like to thank members of RIKEN/RANS for a useful discussion on this study.

References

Allison

Amako

Apostolakis

Arce

Asai

Aso

Bagli

Bagulya

Banerjee

Barrand

Beck

B.R.

Bogdanov

A.G.

Brandt

Brown

J.M.C.

Burkhardt

Canal

Cano-Ott

Chauvie

Cho

Cirrone

G.A.P.

Cooperman

Cortés-Giraldo

M.A.

Cosmo

Cuttone

Depaola

Desorgher

Dong

Dotti

Elvira

V.D.

Folger

Francis

Galoyan

Garnier

Gayer

Genser

K.L.

Grichine

V.M.

Guatelli

Guèye

Gumplinger

Howard

A.S.

Hřivnáčová

Hwang

Incerti

Ivanchenko

V.N.

Jones

F.W.

Jun

S.Y.

Kaitaniemi

Karakatsanis

Karamitros

Kelsey

Kimura

Koi

Kurashige

Lechner

Lee

S.B.

Longo

Maire

Mancusi

Mantero

Mendoza

Morgan

Murakami

Nikitina

Pandola

Paprocki

Perl

Petrović

Pia

M.G.

Pokorski

Quesada

J.M.

Raine

Reis

M.A.

Ribon

Ristić Fira

Romano

Russo

Santin

Sasaki

Sawkey

Shin

J.I.

Strakovsky

I.I.

Taborda

Tanaka

Tomé

Toshito

Tran

H.N.

Truscott

P.R.

Urban

Uzhinsky

Verbeke

J.M.

Verderi

Wendt

B.L.

Wenzel

Wright

D.H.

Wright

D.M.

Yamashita

Yarba

and Yoshida

, Recent developments in Geant4, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment835 (2016), 186–225, https://www-sciencedirect-com-443.web.bisu.edu.cn/science/article/pii/S0168900216306957 . doi:10.1016/j.nima.2016.06.125.

Bair

J.K.

Jones

C.M.

and Willard

H.B.

, Neutrons from the proton bombardment of ⁶Li, ⁷Li, ⁹Be, ¹¹B and ¹⁸O, Nuclear Physics53 (1964), 209. doi:10.1016/0029-5582(64)90598-X.

Bauer

R.W.

Anderson

J.D.

and Wong

, A search for an excited state of ⁹B near 1.7 MeV, Nuclear Physics56 (1964), 117. doi:10.1016/0029-5582(64)90458-4.

Boudard

Cugnon

David

J.-C.

Leray

and Mancusi

, New potentialities of the Liege intranuclear cascade model for reactions induced by nucleons and light charged particles, Phys. Rev. C87 (2013), 014606. doi:10.1103/PhysRevC.87.014606.

Brown

D.A.

Chadwick

M.B.

Capote

Kahler

A.C.

Trkov

Herman

M.W.

Sonzogni

A.A.

Danon

Carlson

A.D.

Dunn

Smith

D.L.

Hale

G.M.

Arbanas

Arcilla

Bates

C.R.

Beck

Becker

Brown

Casperson

R.J.

Conlin

Cullen

D.E.

Descalle

M.-A.

Firestone

Gaines

Guber

K.H.

Hawari

A.I.

Holmes

Johnson

T.D.

Kawano

Kiedrowski

B.C.

Koning

A.J.

Kopecky

Leal

Lestone

J.P.

Lubitz

Damián

J.I.M.

Mattoon

C.M.

McCutchan

E.A.

Mughabghab

Navratil

Neudecker

Nobre

G.P.A.

Noguere

Paris

Pigni

M.T.

Plompen

A.J.

Pritychenko

Pronyaev

V.G.

Roubtsov

Rochman

Romano

Schillebeeckx

Simakov

Sin

Sirakov

Sleaford

Sobes

Soukhovitskii

E.S.

Stetcu

Talou

Thompson

van der Marck

Welser-Sherrill

Wiarda

White

Wormald

J.L.

Wright

R.Q.

Zerkle

Žerovnik

and Zhu

, ENDF/B-VIII.0: The 8th major release of the nuclear reaction data library with CIELO-project cross sections, new standards and thermal scattering data, Nuclear Data Sheets148 (2018), 1–142. Special Issue on Nuclear Reaction Data.

Conlin

J.L.

and Romano

, A Compact ENDF (ACE) Format Specification, Technical Report, LA-UR-19-29016, U.S. Department of Energy, Office of Scientific and Technical Information, 2019. doi:10.2172/1561065.

Gibbons

J.H.

and Macklin

R.L.

, Total neutron yields from light elements under proton and alpha bombardment, Physical Review114 (1959), 571. doi:10.1103/PhysRev.114.571.

Hauser

and Feshbach

, The inelastic scattering of neutrons, Phys. Rev.87 (1952), 366. doi:10.1103/PhysRev.87.366.

Howard

W.B.

Grimes

S.M.

Massey

T.N.

Al-Quraishi

S.I.

Jacobs

D.K.

Brient

C.E.

and Yanch

J.C.

, Measurement of the thick-target ⁹Be(p,n) neutron energy spectra, Nuclear Science and Engineering138 (2001), 145. doi:10.13182/NSE01-A2206.

10.

Iwamoto

and Kunieda

, G-HyND: A hybrid nuclear data estimator with Gaussian processes, Journal of Nuclear Science and Technology59(3) (2022), 334–344. doi:10.1080/00223131.2021.1971120.

11.

Iwamoto

et al., Japanese Evaluated Nuclear Data Library version 5: JENDL-5, 2023, J. Nucl. Sci. Technol., to be published (accepted).

12.

Iwamoto

Kunieda

Minato

and Shibata

, The CCONE code system and its application to nuclear data evaluation for fission and other reactions, Nucl. Data Sheets131 (2016), 259. doi:10.1016/j.nds.2015.12.004.

13.

Kunieda

Iwamoto

Minato

Okamoto

et al., Overview of JENDL-4.0/HE and benchmark calculations, in: Proceedings of the 2016 Symposium on Nuclear Data: JAEA-Conf 2016-004, 2017, pp. 41–46.

14.

MacFarlane

R.E.

and Kahler

A.C.

, Methods for processing ENDF/B-VII with NJOY, Nuclear Data Sheets111(12) (2010), 2739–2890. doi:10.1016/j.nds.2010.11.001.

15.

Otuka

Dupont

Semkova

Pritychenko

Blokhin

A.I.

Aikawa

Babykina

Bossant

Chen

Dunaeva

Forrest

R.A.

Fukahori

Furutachi

Ganesan

Gritzay

O.O.

Herman

Hlavač

Katō

Lalremruata

Lee

Y.O.

Makinaga

Matsumoto

Mikhaylyukova

Pikulina

Pronyaev

V.G.

Saxena

Schwerer

Simakov

S.P.

Soppera

Suzuki

Takács

Tao

Taova

Tárkányi

Varlamov

V.V.

Wang

Yang

S.C.

Zerkin

and Zhuang

, Towards a more complete and accurate experimental nuclear reaction data library (EXFOR): International collaboration between nuclear reaction data centres (NRDC), Nuclear Data Sheets120 (2014), 272–276. doi:10.1016/j.nds.2014.07.065.

16.

Sato

Iwamoto

Hashimoto

Ogawa

Furuta

Abe

S-i.

Kai

Tsai

P.-E.

Matsuda

Iwase

Shigyo

Sihver

and Niita

, Features of Particle and Heavy Ion Transport code System (PHITS) version 3.02, Journal of Nuclear Science and Technology55(6) (2018), 684–690. doi:10.1080/00223131.2017.1419890.

17.

Shibata

Iwamoto

Nakagawa

Iwamoto

Ichihara

Kunieda

Chiba

Furutaka

Otuka

Ohsawa

Murata

Matsunobu

Zukeran

Kamada

and Katakura

, A new library for nuclear science and engineering, J. Nucl. Sci. Technol. (Tokyo)48 (2011), 1–30, JENDL-4.0. doi:10.1080/18811248.2011.9711675.

18.

Soukhovitskiĩ

E.Sh.

Chiba

Capote

Quesada

J.M.

Kunieda

and Morogovskiĩ

G.B.

, Supplement to OPTMAN code, Manual Version 10 (2008), JAEA-Data/Code, 2008-025, Japan Atomic Energy Agency, 2008.

19.

Teranishi

and Furubayashi

, Level width of the ground state of ⁹B, Physics Letters9 (1964), 157. doi:10.1016/0031-9163(64)90124-6.

20.

Trkov

Herman

and Brown

D.A.

(eds), ENDF-6 Formats Mannual: Data Formats and Procedures for the Evaluated Nuclear Data Files ENDF/B-VI, ENDF/B-VII and ENDF/B-VIII, Written by the Members of the Cross Sections Evaluation Working Group, National Nuclear Data Center Brookhaven National Laboratory Upton, NY 11973-5000, USA, 2018.

21.

Wakabayashi

Taketani

Hashiguchi

Ikeda

Kobayashi

Wang

Yan

Harada

Ikeda

and Otake

, A function to provide neutron spectrum produced from the ⁹Be + p reaction with protons of energy below 12 MeV, Journal of Nuclear Science and Technology55(8) (2018), 859–867. doi:10.1080/00223131.2018.1445566.

22.

Werner

C.J.

Bull

J.S.

Solomon

C.J.

Brown

F.B.

McKinney

G.W.

Rising

M.E.

Dixon

D.A.

Martz

R.L.

Hughes

H.G.

Cox

L.J.

Zukaitis

A.J.

Armstrong

J.C.

Forster

R.A.

and Casswell

, MCNP Version 6.2 Release Notes, Technical Report, 2018-02-05, U.S. Department of Energy, Office of Scientific and Technical Information, 2018. doi:10.2172/1419730.

Estimation of double-differential cross-sections of 9 Be(p,xn) reaction for new nuclear data library JENDL-5

Abstract

Keywords

1. Introduction

2. Data evaluation for 9Be(p,xn)

2.1. The Wakabayashi’s function

2.2. Incorporation of the function to nuclear data file

2.5. JENDL-5 and remarks on unsolved issue of p + 9 Be evaluation

Footnotes

Acknowledgements

References

2. Data evaluation for ⁹Be(p,xn)

2.5. JENDL-5 and remarks on unsolved issue of $p +^{9} Be$ evaluation