Cross-section measurements of (n, 2n) and (n,p) reactions on^{124,126,128,130,131,132}Xe in the 14 MeV region and theoretical calculations of their excitation functions

Xenon is a natural element with the atomic numberZ= 54. It has nine stable isotopes with Z numbers following that of tin. Its isotopes are¹²⁶Xe(0.0890%),¹²⁴Xe(0.0952%),¹²⁸Xe(1.9102%),¹³⁰Xe(4.0710%),¹³⁶Xe(8.8573%),¹³⁴Xe(10.4357%),¹³¹Xe(21.232%),¹²⁹Xe(26.4006%), and¹³²Xe(26.9086%) [1]. These isotopes are prominent in the study of nuclear energy and nuclear physics. First, xenon is ideal for studying plasma conditions at high temperatures and pressure owing to its remarkable physical properties. In the design of a fusion reactor, adding an appropriate xenon gas and harnessing its pressure resistance and thermal neutron absorption characteristics can effectively prevent nuclear leakages, thereby ensuring the safety of nuclear reactors such as the international thermonuclear experimental reactor (ITER) and accelerator-driven subcritical systems (ADSs) [2]. The³H(d,n)⁴He fusion could provide output energy up to 14 MeV. In this case, the neutron yields would be (5~6)×10¹⁴n/s [3]. Therefore, the xenon gas activation inside a fusion reactor should be considered. Second, xenon isotopes are ascribed to the transitional region below the closed neutron shell withN= 82 [4]. The analysis of the transitional nucleus neutron cross-section is a fundamental interest of nuclear chemists and physicists. The possibility of comparing these competing approaches raises an additional interest from the experimentalists to obtain data on even–even Xe nuclei. Particularly important are the non-yrast level properties that provide strict tests for various nuclear models [5]. Nuclear reactions can be inducedvianeutrons with an energy of approximately 14 MeV (e.g., (n, 2n), (n,p), (n,α), (n,t), and (n,d) reactions [6,7]). Such reactions are essential for nuclear physics to elucidate the inertial confinement fusion (ICF) plasma [8,9]. The neutron′s energy distributed in the deuterium-tritium (DT) ICF plasma relies on the power of stopping for the corresponding tritons and deuterons. For this purpose, the cross-section values of xenon isotope nuclear reactions induced by 14 MeV neutrons are required. However, the preparation of appropriate rare-gas targets for irradiation is difficult, which explains the lack of data on the corresponding cross-sections. For the¹²⁴Xe isotopes, the cross sections of the¹²⁴Xe(n,2n)¹²³Xe reaction at approximately 14 MeV have been determined by a limited number of laboratories [10-13]. The cross-sections of the¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe,¹³²Xe(n, 2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I,¹³¹Xe(n,p)¹³¹I, and¹³²Xe(n,p)¹³²I reactions at ~14 MeV were obtained experimentally by only two groups [10,11]. They adopted sodium perxenate (Na₄XeO₆·2H₂O) and quinol-clathrates as the target sample, and the corresponding result was reported at only one energy point. In addition, significant differences have been identified among previous literature results [10-14]. These differences are mainly triggered by two factors:

1) Decay data deficiencies (characteristic γ-ray intensity). For the¹²⁴Xe(n,2n)¹²³Xe reaction, the cross-sections were measured by Sigg and Kuroda (1976) [10] and Kondaiahet al. (1968) [11] with the characteristic γ-ray intensity at 63.9% and 47.3% for the same characteristic γ-ray of 148.9 keV, respectively (refer toTable 1). Ref. [12] provides no information on decay. The datum for 48.9% was recently measured [1]. For the¹²⁷Xe product radionuclide, the cross-sections have also been determined by Sigg and Kuroda (1976) [10], with the characteristic γ-ray intensity at 60.8% for the same characteristic γ-ray of 202.8 keV (refer toTable 1). The data for 68.7% have been recently measured [1]. For the^129mXe product radionuclide, the cross-sections have been determined by Kondaiahet al. (1968) [11] using the characteristic γ-ray intensity to be 5.8% for the same characteristic γ-ray 196.56 keV. The data for 4.59% have been recently measured [1].

2) Interfering reactions. For the¹³²Xe(n,p)¹³²I (β^-(100%),T_1/2=2.295 h) reaction, in previous research [10,11], the 773 keV gamma-ray emitted in the¹³²I decay was used to obtain the values of the¹³²Xe(n,p)¹³²I reaction cross section. However, the 773 keV gamma-ray did not only emerge from the¹³²I product radionuclide, but also from the¹³²Xe(n,p)^132mI (half-lifeT_1/2= 1.387 h,γ-ray energyE_γ= 772.6 keV, and intensityI_γ= 14.0%) reaction, and from the¹³⁴Xe(n,α)^131m1Te (half-lifeT_1/2= 33.25 h,γ-ray energyE_γ= 773.67 keV, and intensityI_γ= 36.8%) reaction (refer toTable 1, the decay data of the daughter nucleus, the detector used, and the monitor reaction are also given).

Owing to the lack of data and differences between these data, the libraries of evaluated nuclear data rely heavily on nuclear model calculations in the 13–15 MeV energy regions; hence, large differences exist among the various evaluations. For instance, for the¹²⁴Xe(n, 2n)¹²³Xe reaction, the evaluation result of database JENDL-4.0 [15] was 1.5 times that of ENDF/B-VIII.0 [16] (CENDL-3 [17]). For the¹³⁰Xe(n,p)¹³⁰I reaction, the result of database JENDL-4.0 [15] was 3 times that of ENDF/B-VIII.0 [16] (JEFF-3.3 [18], RUSFOND [19]). For the¹³¹Xe(n,p)¹³¹I reaction, the result of database JENDL-4.0 [15] was 4 times that of ENDF/B-VIII.0 [16], JEFF-3.3 [18], RUSFOND [19], and CENDL-3 [17]. For the¹³²Xe(n,p)¹³²I reaction, the result of database JENDL-4.0 [15] was 4 times that of RUSFOND [19]. Consequently, new measurements for the Xe(n,x) cross section are required in the 13–15 MeV neutron energy range.

This paper reports the¹²⁴Xe(n, 2n)¹²³Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe,¹²⁶Xe(n, 2n)¹²⁵Xe,¹³⁰Xe(n,p)¹³⁰I,¹³²Xe(n, 2n)^131mXe, and¹³¹Xe(n,p)¹³¹I reaction cross-section measurements at 13–15 MeV. We adopted natural high-purity pressured Xe gas as a target, and a high-resolution HPGe detector was utilized to measure the gamma activities of the product nuclei. During the irradiation process, each sample was wrapped in pure cadmium foils to circumvent the effects of the¹²⁴Xe(n,γ)¹²⁵Xe,¹²⁶Xe(n,γ)¹²⁷Xe,¹²⁸Xe(n,γ)^129mXe, and¹³⁰Xe(n,γ)^131mXe reactions induced by thermal neutrons on the¹²⁶Xe(n,2n)¹²⁵Xe,¹²⁸Xe(n,2n)¹²⁷Xe,¹³⁰Xe(n,2n)^129mXe, and¹³²Xe(n,2n)^131mXe reactions, respectively. The obtained results were compared and analyzed with previously reported results, the evaluation data in JENDL-4.0 [15], ENDF/B-VIII.0 [16], CENDL-3 [17], JEFF-3.3 [18], and RUSFOND [19] libraries, and the theoretical excitation curves acquired using the TALYS-1.95 [20] code.

The procedure for radioactive product identification and cross-section measurements is outlined in a different paper [21-23]. Only the most prominent characteristics associated with the measurements conducted in this study are presented.

Xenon gas (99.999% purity with natural isotopic composition) was stored in a stainless-steel sphere with a 20 mm inner diameter and 1 mm wall thickness. The xenon gas weight was obtained by subtracting the container weight from the total weight. The gas weight and pressure were 2.3–5.2 g and 120–250 atm, respectively. Five of such samples of xenon gas and 10 monitors (20 mm diameter), comprising 99.99% pure and 0.12 mm thick Nb foil and a 99.999% pure and 0.3 mm thick Al foil, were prepared for irradiation. A ZrNbAl–Xe–AlNbZr stack was prepared from these five samples and mounted at 0°, 45°, 90°, 110°, and 135° angles relative to the beam. The end of the neutron production target was 5 cm apart from the sample sphere center (seeFig. 1). The neutron fluxes for the low-threshold of the¹³⁰Xe(n,p)¹³⁰I,¹³¹Xe(n,p)¹³¹I, and¹³²Xe(n,p)¹³²I reactons were monitored by implementing the²⁷Al(n, α)²⁴Na reaction (E_th= 3.249 MeV). The neutron fluxes for the high-threshold of the¹²⁴Xe(n, 2n)¹²³Xe,¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe, and¹³²Xe(n, 2n)^131mXe reactions were monitored by implementing the⁹³Nb(n, 2n)^92mNb reaction cross-section (E_th= 8.792 MeV). The cross-sections of the²⁷Al(n, α)²⁴Na and⁹³Nb(n, 2n)^92mNb monitor reactions were adopted from a different paper [24]. Subsequently, a cadmium foil (99.95% pure and 1 mm thick) was wrapped around the samples to minimize the influence of the¹²⁴Xe(n, γ)¹²⁵Xe,¹²⁶Xe(n, γ)¹²⁷Xe,¹²⁸Xe(n, γ)^129mXe, and¹³⁰Xe(n,γ)^131mXe reactions activatedviathermal neutrons on the¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe, and¹³²Xe(n, 2n)^131mXe reactions, respectively.

Figure 1.(color online) Sample position with the target assembly.

The K-400 neutron generator at the Institute of Nuclear Physics and Chemistry of the China Academy of Engineering Physics was utilized to irradiate the samples for 2 h. The neutron yields were ~(4–5)×10¹⁰n/s. The³H(d,n)⁴He reaction produced 14 MeV neutrons when excited by a 134 keV deuteron beam energy and a 200 μA beam current. The generator described above uses a solid³H–Ti target that is 2.65 mg/cm²thick. To monitor the neutron flux within the irradiation, the accompanyingα-particles were utilized. An Au–Si surface barrier detector was set at an angle of 135° relative to the beam direction. The obtained data were adopted for the correction of the neutron flux deviation.

The average energies of neutron during the irradiating tests at the 0°, 45°, 90°,110°, and 135° incident angles were established via two approaches. The first approach adopted the cross-section ratios of the⁹³Nb(n, 2n)^92mNb and⁹⁰Zr(n, 2n)^89m+gZr reactions [25]. Another approach utilized the expressions obtained from Ref. [23]. The results obtained by the two approaches are consistent within their uncertainties. These neutron energies were 14.8, 14.4, 14.1, 13.8, and 13.5 MeV, respectively. The neutron energy uncertainty at 5 cm was assessed to be equal to 0.2 MeV, following the consideration of the specimen and the d⁺beam sizes (~4 mm) [25].

After irradiation, the samples decay for different amounts of time (ranging from 55 min to 5 days depending on the half-life of the daughter nucleus), and the gamma-ray activities of the¹²³Xe,¹²⁵Xe,¹²⁷Xe,^129mXe,^131mXe,¹³⁰I,¹³¹I,²⁴Na,^92mNb, and¹³²I nuclei were measured by the calibrated detector of GEM-60P HPGe (with the crystal length and diameter equal to 72.3 and 70.1 mm, respectively). The energy resolution at 1.332 MeV and a relative efficiency were equal to 1.69 keV and about 68%, respectively. The detector was calibrated using standard sources ofγ-ray. Monte Carlo simulations were also performed to correct the geometry difference between theγ-ray calibration sources used for determining the efficiency of the HPGe detectors and the Xe sphere. The partial spectra ofγ-ray for Xe samples collected after 4.5 h and 32.7 h irradiations are presented inFigs. 2and3, respectively. In addition to the strong lines of interest at 148.9 and 536.07 keV inFig. 2and 188.4, 202.86, 163.93, 196.56, and 364.49 keV inFig. 3, in the former figure, the weaker¹³²Iγ-ray transitions at 522.65 keV are also indicated. The ORTEC® GammaVision® Gamma Spectrum Analysis computer program examined the peak area [26] using the emulation package of ORTEC® MAESTRO® MCA for data collection and analysis [26].

Figure 2.(color online) γ-ray spectrum of xenon obtained with a 68% relative efficient HPGe detector atE_n=13.5 MeV, 4.5 h after the neutron activation of xenon isotopes and an acquisition time of 4.9 h. The lines of interest are 148.9, 522.6, and 536.07 keV. In the energy ranges of 120–180 keV and 500–600 keV, the enlarged portion of the spectrum is given.

Figure 3.(color online) γ-ray spectrum of xenon obtained with a 68% relative efficient HPGe detector atE_n=14.8 MeV, 32.7 h after the neutron activation of xenon isotopes, and an acquisition time of 10.2 h. The lines of interest are 188.4, 202.86, 163.93, 196.56, and 364.49 keV. In the energy ranges of 350–550 keV, the enlarged portion of the spectrum is given.

Table 2presents the above mentioned reactions, the characteristics of the corresponding reaction product radioactive decay, and their abundance.

The reaction cross-sections were calculated as follows [27-31]:

where the subscripts0andxdenote the standard monitor and measured reactions values, respectively, whileFrepresents an activity correction factor:

wheref_c,f_s,f_g, andf_Ωdenote the correcting factors representing the coincidence summing impact of the cascadeγ-rays created by the nuclides, the self-absorption of the γ-rays, the samples counting geometries, and the solid angle of neutron flux, respectively. The calculations of the coincidence summing correction factors were based on the approach explained in previous studies [32,33].

Some of these factors were calculated by the following weighted average formula:

whereμdenotes the coefficient of linear attenuation at photon energyE(seeTable 2),Lis the distance from the specimen to the HPGe surface, while Ω₁, Ω₂, Ω₃represent the solid angle of the front monitor disk, sample and back monitor disks relative to the center of the target, respectively. The coefficient of mass attenuation (μ/ρ) was obtained using the γ-energy taken from a previous study [34]. The linear attenuation coefficient of natural xenon gas was determined as $ \mu =\rho (\mu /\rho ) $ , whereρ(in g/cm³) represents the Xe gas density in stainless-steel spheres. The linear attenuation coefficient of stainless steel was determined as $\mu = 7.86(\mu /\rho )$ , where the density of iron is 7.86 g/cm³. The comprehensive self-absorption correction factors are presented inTable 3, and the results are distributed between 1.072 and 2.014. The value of the correction factorf_Ωis 1.0756.

In the process of calculating the cross sections of the¹³⁰Xe(n,2n)^129mXe,¹³²Xe(n,2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I, and¹³¹Xe(n,p)¹³¹I reactions, the contributions of the interfering¹²⁹Xe(n, n′)^129mXe,¹³¹Xe(n, n′)^131mXe,¹³¹Xe(n, d)¹³⁰I, and¹³²Xe(n, d)¹³¹I reactions were respectively subtracted using the following equations [7,35]:

The interfering¹³¹Xe(n, d)¹³⁰I and¹³²Xe(n, d)¹³¹I reaction cross-sections were 0.09, 0.15, 0.22, 0.32, and 0.50 mb and 0.30, 0.41, 0.56, 0.80, and 1.11 mb, respectively, at 13.5, 13.8, 14.1, 14.4, and 14.8 MeV, respectively. These values were obtained from the interpolations performed using JEFF-3.3 [16]. The minor contributions to the activities of the product nuclei¹³⁰I and¹³²I from the interfering¹³²Xe(n, t)¹³⁰I and¹³⁴Xe(n, t)¹³²I reactions were neglected because their cross-sections were below 0.04 mb in the 13–15 MeV range [16]). Unfortunately, neither the experimental measurement nor evaluated data were available for the¹²⁹Xe(n, n′)^129mXe and¹³¹Xe(n, n′)^131mXe reactions. Hence, theoretical calculations were conducted to infer these data.

Our primary sources of uncertainty were the counting statistics (0.2%–22%), detector efficiency (3%), cross section of the monitor reaction (0.4%–1.5%), sample weight (0.05%–0.9%), measurements and decay times (0.1%), and the self-absorption of γ-ray (1%). Other uncertainties originated from the experimental and standard nuclear parameters, including the product nuclear half-lives (0.01%–1.1%), specificγ-ray branching ratios (0.2%– 3.1%), and target isotope abundance (0.02%–0.3%). We adopted a quadratic sum rule for the uncertainty analysis [36].

The cross sections of the¹²⁴Xe(n, 2n)¹²³Xe,¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe,¹³²Xe(n, 2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I,¹³¹Xe(n,p)¹³¹I, and¹³²Xe(n,p)¹³²I reactions were assessed by the TALYS-1.95 software [20] up to 20 MeV. TALYS-1.95 is a nuclear model code applied to replicate nuclear reactions comprising tritons, neutrons, photons, protons, deuterons,³He, andα-particles aimed at the target nuclei. It can be adopted for the reactions with incident energies up to 200 MeV [37-41]. The different parameters in the TALYS-1.95 code were adjusted according to our measured data and previous experiments conducted by other researchers for the¹²⁴Xe(n,2n)¹²³Xe,¹²⁶Xe(n,2n)¹²⁵Xe,¹³⁰Xe(n,2n)^129mXe,¹³²Xe(n,2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I,¹³¹Xe(n,p)¹³¹I, and¹³²Xe(n,p)¹³²I reactions.

The offlineγ-ray spectroscopic measuring technique was adopted to obtain cross-sections for the¹²⁴Xe(n, 2n)¹²³Xe,¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe,¹³²Xe(n, 2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I,¹³¹Xe(n,p)¹³¹I, and¹³²Xe(n,p)¹³²I reactions. The obtained results were then examined with those from the JENDL-4.0 [15], ENDF/B-VIII.0 [16], CENDL-3 [17], JEFF-3.3 [18], and RUSFOND [19] data libraries and TALYS-1.95 model code [20]. For the first time, this paper reports the cross-sections of several reactions (excluding the¹²⁴Xe(n, 2n)¹²³Xe and¹³²Xe(n,p)¹³²I) at 13.5, 13.8, and 14.1 MeV (seeTable4).

For the¹²⁴Xe(n, 2n)¹²³Xe reaction, the gamma-ray with an intensity ofI_γ=48.9% and 148.9 keV emitted in the decay of¹²³Xe was utilized to obtain the cross-section of the¹²⁴Xe(n, 2n)¹²³Xe reaction. Sigg and Kuroda [10] usedI_γ=63.9% while Kondaiahet al. [11] usedI_γ=47.3% for the same ray. However, studies demonstrate that the obtained data forI_γ=48.9% exhibits higher measuring accurately. Accordingly, the cross-sections are modified as follows:

where σ_literatureis the literature cross-section,I_γ(old)represents the corresponding intensity of the 148.9 keV gamma-ray in the literature, andI_γ(new)is 48.9%. Consequently, the divergence of the four evaluation curves given by databases JENDL-4.0 [15], ENDF/B-VIII.0 [16], CENDL-3 [17], JEFF-3.3 [18], and RUSFOND [19] is very obvious (ENDF/B-VIII.0 [16] is the same as CENDL-3 [17]). The evaluation result of database JENDL-4.0 [15] was 1.5 times that of ENDF/B-VIII.0 [16] (CENDL-3 [17]). The evaluation curves of the ENDF/B-VIII.0 [16] and CENDL-3 [17] data libraries were approximately similar to the simulated curves of excitation processed using TALYS-1.95 [20] (refer toFig. 4). The results obtained from our measurements increase upon increasing the neutron energy in the 14 MeV range, similar to the theoretical and evaluation data. However, some differences were still observed between the results of earlier experiments. Only four laboratories have reported neutron-induced experimental cross-sections regarding the reaction of¹²⁴Xe(n, 2n)¹²³Xe [10–13]. In the 13–15 MeV neutron region of energy, the results obtained in our study are consistent (within the experimental uncertainties) with those obtained by the TALYS-1.95 code (ldmodel 2: The back-shifted Fermi-gas model), with those obtained from the RUSFOND [19] assessments, and with the values reported by Sigg and Kuroda [10] and Kondaiahet al. [11]. In addition, the evaluated data obtained from JENDL-4.0 [15] and Bazan [12] exhibited larger cross-sections than those obtained from the TALYS-1.95 code (ldmodels 2 and 4), ENDF/B-VIII.0 [16] and CENDL-3 [17] assessments, and the study by Bhikeet al. [13].

Figure 4.(color online) Cross-section of¹²⁴Xe(n,2n)¹²³Xe reaction.

For the¹²⁶Xe(n, 2n)¹²⁵Xe reaction, only two laboratories reported the cross section data in the previous period, and both provided only one energy point [10,11]. In this study, the¹²⁶Xe(n, 2n)¹²⁵Xe reaction cross-sections were determined via the⁹³Nb(n, 2n)^92mNb monitor reaction with a high threshold. The samples were wrapped in a Cd foil to minimize the influence of the low-energy neutrons derived from²H agglomeration within the tritium target over time and background neutrons. Furthermore, the OMP parametersr_vanda_vwere adjusted.Figure 5presents the cross-sections of the¹²⁶Xe(n, 2n)¹²⁵Xe reaction together with calculations obtained from the TALYS-1.95 code and ENDF/B-VIII.0 [16] (RUSFOND [19]), JEFF-3.3 [18], and JENDL-4.0 [15] libraries (marked as solid lines). At 13.5 and 13.8 MeV, the data measured in this study are consistent with the assessed outcomes from the JENDL-4.0 [15] libraries, and within the limit of experimental uncertainties. In the 13-15 MeV regions, our data are higher than the literature values [10,11] (refer toTable 4andFig. 5), as well as the assessed outcomes obtained from the ENDF/B-VIII.0 [16] (RUSFOND [19]) and JEFF-3.3 [18] libraries, and the results of the TALYS-1.95 calculation. ENDF/B-VIII.0 [16] and RUSFOND [19] databases provide the same results. At neutron energies of 13.5, 13.8, and 14.1 MeV, our results are reported for the first time.

Figure 5.(color online) Cross-section of the¹²⁶Xe(n,2n)¹²⁵Xe reaction.

To date, only two laboratories have experimentally determined and reported the cross-sections of the¹²⁸Xe(n, 2n)¹²⁷Xe reaction at one energy point. These values could be found in the corresponding nuclear reaction databases at 14 MeV. Early measurement [10] was carried out using the 202.8 keV (I_γ=60.8%) gamma-ray emitted in the¹²⁷Xe. However, for the same ray, data for theI_γ=68.7% has been measured accurately in recent years. Hence, these literature cross sections are corrected by Eq, (10). These results provide a solid foundation to verify the reliability of the experimental results and the validity of the theoretical calculations presented here. In the theoretical calculation, the models for ther_vanda_vparameters of the optical model potential (OMP) were adjusted. The excitation curves in the ENDF/B-VIII.0 [16] (which is similar to RUSFOND [19]) and JEFF-3.3 [18] (JENDL-4.0 [15]) databases are very similar to the theoretical curves of the excitation attainedviaTALYS-1.95 [20] within the range of neutron energies with a 15 MeV threshold (refer toFig. 6). Only minor differences among these curves can be observed. Within the limits of experimental uncertainties, the experimental data obtained here agree with the evaluation curves obtained from JEFF-3.3 [18] (JENDL-4.0 [15]) and the corrected cross-sections of Sigg and Kuroda [10], including the experiment values reported by Kondaiahet al. [11].

Figure 6.(color online) Cross-section of the¹²⁸Xe(n,2n)¹²⁷Xe reaction.

The IAEA database does not contain any evaluation for the cross sections of the¹³⁰Xe(n, 2n)^129mXe. Only two experimental cross sections were reported for the 13.5–14.8 MeV range [10,11], and these two reports provide the cross-sections at either 14.6 or 14.4 MeV. The OMP parametersr_vanda_vfor the¹³⁰Xe(n, 2n)^129mXe reaction, as well as the model for the level densities, were adjusted. These cross-sections are illustrated inFig. 7. Our data for the 13.5–14.8 MeV range disagreed with the results of experimental papers [10,11] even when considering the experimental uncertainties. FromFig. 7, it can also be observed that our values are lower than those presented in Refs. [10,11]; however, in the 13–15 MeV region, the values obtained here are higher than those obtained from the TALYS-1.95 code. Furthermore, at the energy of the 14.6 MeV, the values reported by Sigg and Kuroda [10] are approximately 1.6 times greater than the values obtained in this research. This discrepancy was very likely due to the interference of the¹²⁹Xe(n,n′)^129mXe reaction. To date, neither experimental measurements nor evaluation results have been reported for the¹²⁹Xe(n,n′)^129mXe reaction. However, according to the theoretical calculation findings of the TALYS-1.95 code (utilizing default parameters), the cross-section values of this reaction at 13–15 MeV are between 143 and 166 mb.

Figure 7.(color online) Cross-section of¹³⁰Xe(n,2n)^129mXe reaction.

For the¹³²Xe(n, 2n)^131mXe reaction, only two works [10,11] have reported cross-sections, both of which provide the cross-section data at one energy value. Although this reaction is also affected by the interfering¹³¹Xe(n,n′)^131mXe reaction, its impact is not as large as that of the interfering¹²⁹Xe(n,n′)^129mXe reaction. The effect of the cross-section value of the¹³¹Xe(n,n′)^131mXe reaction at 13–15 MeV is 106–121 mb (theoretical results of the TALYS-1.95 with default parameters). According to Eq. (6), this value was considered in this work. Similarly,Fig. 8presents the cross section of the¹³²Xe(n, 2n)^131mXe reaction. In addition, for the energy range between 13.3 and 15 MeV, our experimental results are in excellent agreement with the theoretical calculation (ld-model 5: comprising microscopic level densities (Skyrme force) from Hilaire’s combinatorial tables [41]) and the results obtained by Sigg and Kuroda [10] within their uncertainty limits (seeFig. 8); however, these results were larger than the experiment values reported by Kondaiah et al. [11].

Figure 8.(color online) Cross-section of¹³²Xe(n,2n)^131mXe reaction.

Figure 9presents the cross sections of the¹³⁰Xe(n,p)¹³⁰I reaction based on the results reported in [10,11]. InFig. 9, the values obtained using the TALYS-1.95 calculations and the data obtained from the ENDF/B-VIII.0 [16] (JEFF-3.3 [18] and RUSFOND [19]) and JENDL-4.0 [15] libraries are represented in continuous lines. As illustrated inFig. 9, there are significant dissimilarities among the results obtained from TALYS-1.95 [20], including the evaluation values acquired using JENDL-4.0 [15] at a neutron energy of approximately 14 MeV. The values evaluated by the JENDL-4.0 [15] for the¹³⁰Xe(n,p)¹³⁰I reaction are 3 times higher those assessed by ENDF/B-VIII.0 [16], JEFF-3.3 [18], and RUSFOND [19] at a neutron energy of approximately 14 MeV. Above the 15 MeV neutron energy, in addition to becoming significantly larger, the difference between the two evaluation curves also changes in the opposite direction. The main reason for this phenomenon is the lack of necessary experimental data. Furthermore, the present result is consistent with the results reported in [11] within the limits of experimental uncertainties. At the 14.8 MeV energy point, the results obtained by Sigg and Kuroda [10] were 40% higher than those measured in this study.

Figure 9.(color online) Cross-section of¹³⁰Xe(n,p)¹³⁰I reaction.

For the¹³¹Xe(n,p)¹³¹I reaction, the following conditions were adopted for the measurements: gamma-rays at 364.49 keV (Iγ = 81.5%) emitted via¹³¹I (half-life,T_1/2= 8.0252 d). The measurements reported in the literature adopt the same decay data as this work.Figure 10compares our data with those of the and literature [10,11], together with the data obtained by TALYS-1.95 and from the JENDL-4.0 [15], ENDF/B-VIII.0 [16], CENDL-3 [17], JEFF-3.3 [18], and RUSFOND [19] libraries. The agreement between these data is the best at the 13–15 MeV range (seeFig. 10), despite the lack of experimental data at an early stage. Within the neutron energy rage of 13.5–14.8 MeV, our values agree well with the experimental data obtained in [10] and with the evaluated data obtained from the CENDL-3 [17] libraries, within the limits of experimental uncertainty. However, all of our data are below those of JENDL-4.0 [15] (refer toFig. 10), and the result of the JENDL-4.0 database [15] was 4 times that of ENDF/B-VIII.0 [16], JEFF-3.3 [18], RUSFOND [19], and CENDL-3 [17]. For the theoretical calculation of the¹³¹Xe(n,p)¹³¹I reaction, although the OMP parametersr_vanda_v, as well as the level densities, were adjusted to 2, 2, and ld-model 3, respectively, the results remain slightly higher than those obtained in the present study and the evaluated data obtained from the ENDF/B-VIII.0 [16], CENDL-3 [17], JEFF-3.3 [18], and RUSFOND [19] libraries.

Figure 10.(color online) Cross-section of¹³¹Xe(n,p)¹³¹I reaction.

In this study, the latest decay data (half-lifeT_1/2= 2.295 h, characteristicγ-ray energyE_γ= 522.65 keV, andγ-ray intensityI_γ= 16.0 %) of the¹³²I product were used to deduce the cross-section value of the¹³²Xe(n,p)¹³²I reaction. The characteristic gamma-ray energy of 773 keV was selected in [10,11] to determine the activity of the¹³²I product. However, this line originates from both¹³²I and also from^132mI (half-lifeT_1/2= 1.387 h,γ-ray energyE_γ= 772.6 keV, and intensityI_γ= 14.0 %) of the¹³²Xe(n,p)^132mI reaction, including the^131m1Te (half-lifeT_1/2= 33.25 h,γ-ray energyE_γ= 773.67 keV, and intensityI_γ= 36.8 %) of the¹³⁴Xe(n, α)^131m1Te reaction. The results obtained in this study are presented inFig. 11, together with all other data from previous reports [10,11], theoretical calculation results form the TALYS-1.95 (ldmodel 2) code, and the evaluation data from the ENDF/B-VIII.0 [16], CENDL-3 [17], JEFF-3.3 [18], and RUSFOND [19] libraries. As can be observed fromFig. 11, the divergence of the three group evaluation data curves and theoretical curve is relatively large, and this difference gradually becomes larger with an increase in neutron energy. At around 14 MeV, the result of the JENDL-4.0 database [15] was 4 times that of RUSFOND [19]. In this case, the experimental results play an important role in eliminating the differences between the evaluation data and optimizing the parameters of the nuclear theoretical model.

Figure 11.(color online) Cross-section of¹³²Xe(n,p)¹³²I reaction.

The cross-sections of multiple reactions including¹²⁴Xe(n, 2n)¹²³Xe,¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe,¹³²Xe(n, 2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I,¹³¹Xe(n,p)¹³¹I, and¹³²Xe(n,p)¹³²I reactions were measured for the 13–15 MeV neutron energies via the activation approach. The excitation functions for these reactions up to a neutron energy of 20 MeV were calculated by the Talys-1.95 program using adjusted parameters. The obtained results were compared with previous experiments, evaluated data, and theoretical results based on the Talys-1.95 code. In general, within the experimental uncertainties, the present cross-section values were consistent with those of previous experiments and the theoretical excitation curves at corresponding energies. Our novel achievements are expected to facilitate the quality of the neutron cross-section databases and assist possible novel evaluations of the cross sections of xenon isotopes in the incident neutron energy range from the threshold energy to 20 MeV. Moreover, reliable cross section data can play a key role in the verification of nuclear reaction models and the determination of model parameters, especially nuclear reactions for which experimental data are scarce. Evidently, for the first time, this study reports the experimental values of cross-sections at 13.5, 13.8, and 14.1 MeV for the¹²⁶Xe(n, 2n)¹²⁵Xe,¹²⁸Xe(n, 2n)¹²⁷Xe,¹³⁰Xe(n, 2n)^129mXe,¹³²Xe(n, 2n)^131mXe,¹³⁰Xe(n,p)¹³⁰I, and¹³¹Xe(n,p)¹³¹I reactions, as well as at 13.5 MeV for the¹³²Xe(n,p)¹³²I reaction. Hence, this study addresses the limitations of the single energy of some cross section data and also broadens the energy range of neutron cross section data.

We would like to thank the Intense Neutron Generator group at the China Academy of Engineering Physics for performing the irradiations.