Cross-section measurements for^58,60,61Ni(n,α)^55,57,58Fe reactions in the 4.50 – 5.50 MeV neutron energy region

Nickel is widely used for making stainless steel and other corrosion-resistant alloys which play an important role in the construction of nuclear reactors and accelerators. The⁵⁸Ni,⁶⁰Ni, and⁶¹Ni contents in natural nickel are 68.0769%, 26.2231%, and 1.1399%, respectively [1]. Research into the (n,α) reactions is important in nuclear engineering applications, because neutron-induced helium production would lead to helium accumulation and cause serious radiation damage to the materials. Besides, measurements of these cross sections could enhance our understanding of nuclear structure and nuclear reaction mechanisms. For example, the cross section of the⁵⁸Ni(n,α)⁵⁵Fe reaction would enable one to derive the level structure of the residual nucleus,⁵⁵Fe [2].

For the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe (Q= 2.899 MeV) reaction, existing measurement results are abundant because the activation method is available for this reaction. In the 4.50 – 5.50 MeV neutron energy region, six measurements [2-7] could be found in the EXFOR library [8], but there are noticeable discrepancies among these measurements. For example, the cross section of the⁵⁸Ni(n,α)⁵⁵Fe reaction measured by Gledenov (1997, 5.00 – 7.00 MeV) [6] is ~ 2 times higher than that obtained by Goverdovskiy (1992, 5.12 MeV) [4]. Besides, the coefficient of variation among different evaluations, including the ENDF/B-VIII.0 [9], ENDF/B-VII.1 [10], JENDL-4.0u+ [11], JEFF-3.3 [12], ROSFOND-2010 [13], and CENDL-3.1 [14] libraries, is 14.92% in the 4.50 MeV ≤E_n≤ 5.50 MeV region [15].

For the cross sections of the⁶⁰Ni(n,α)⁵⁷Fe (Q= 1.355 MeV) and⁶¹Ni(n,α)⁵⁸Fe (Q= 3.580 MeV) reactions, because the residual nuclei (⁵⁷Fe and⁵⁸Fe) of these two reactions are stable, the activation method is unavailable for measurement. For the cross sections of the⁶⁰Ni(n,α)⁵⁷Fe reactions, there is only one measurement (Khromyleva (2018, 6.00 – 7.15 MeV) [16]) in the EXFOR library [8]. For the⁶¹Ni(n,α)⁵⁸Fe reaction, there is no data in the entire neutron-energy region except for several results aroundE_n= 0.0253 eV [8]. In the 5.00 MeV ≤E_n≤ 5.50 MeV region, the coefficients of variation among different evaluations are 36.43% and 46.85% for the⁶⁰Ni(n,α)⁵⁷Fe and⁶¹Ni(n,α)⁵⁸Fe reactions, respectively [9-14].

Accurate measurements of the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions are required. In the present work, the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe reaction were measured at theE_n= 4.50, 4.75, 5.00, 5.25, and 5.50 MeV energy points and those of the⁶⁰Ni(n,α)⁵⁷Fe and⁶¹Ni(n,α)⁵⁸Fe reactions were obtained at theE_n= 5.00 and 5.50 MeV energy points. The present data for the⁶⁰Ni(n,α)⁵⁷Fe reaction are the first results of measurements below 6.0 MeV and those of the⁶¹Ni(n,α)⁵⁸Fe reactions are the first results of measurement in the MeV region. The details of the experiments are illustrated in Sect. 2, the data processing and the results are presented in Sect. 3 and Sect. 4, respectively, and the conclusion is drawn in Sect. 5.

The experiments were performed based on the 4.5-MV Van de Graaff accelerator at Peking University, China. As shown inFig. 1, the experimental apparatus consists of three main parts: the neutron source, the gridded twin ionization chamber (GIC) as the charged particle detector (with samples inside), and the scintillator detector for neutron energy spectrum measurement. The symmetric double sections of the GIC were defined as side 01 and side 02 , respectively.

Figure 1.(color online) Schematic drawing of the experimental apparatus.

Measurements were performed based on the 4.5-MV Van de Graaff accelerator at Peking University. The quasi mono-energetic neutrons were generated by the²H(d,n)³He reaction using the energetic deuteron beam from the accelerator to bombard a deuterium gas target 2.0 cm in length under a pressure of 3.0 atm. The deuterium gas target was separated from the vacuum tube of the accelerator by a 5.0 μm-thick molybdenum foil sample. The energy range of the incident deuteron was set to 2.091 – 2.868 MeV for the neutrons to be generated with the kinetic energy range of 4.50 – 5.50 MeV and an energy spread (1σ) of 0.14 – 0.19 MeV [17]. The deuteron beam current was ~ 2.0 μA throughout the measurement.

A sample changer with five sample positions was set at the common cathode of the GIC, and back-to-back double samples were placed at each of them as presented inTable 1. The sample changer could be rotated with a rotating a knob below the GIC without opening it. The samples used in the present measurements were as follows: a) the back-to-back compound α sources at sample position No.1 (²³⁴U, 4.775 MeV;²³⁹Pu, 5.155 MeV;²³⁸Pu, 5.499 MeV;²⁴⁴Cm, 5.805 MeV) were used to calibrate the detection system [18]; b) to monitor the neutron flux, three highly enriched (99.999%)²³⁸U₃O₈samples were prepared, and their nucleus number and unevenness were determined using theαspectrum measured by the GIC. Details of the method can be found in Ref. [19]. A²³⁸U₃O₈sample 43.0 mm in diameter, 62.1% in unevenness, and 600.8 μg/cm²in average thickness at the No.2 sample position was used to determine the absolute neutron flux. The other two²³⁸U₃O₈samples were glued onto the fission cathodes of the 01 and 02 sides to determine the relative neutron flux, respectively. Their diameters were 45.0 mm and 43.0 mm, unevennesses were 61.6% and 63.3%, and average thicknesses were 604.6 μg/cm²and 557.5 μg/cm², respectively. c) Two⁵⁸Ni samples at the No.3 sample position, and one⁶⁰Ni sample and one⁶¹Ni sample at the No.4 sample position were used to measure the foreground events, the data of which are listed inTable 2. The number of Ni atoms in the samples was determined by weighing the samples. Each Ni sample was prepared on a 0.1 mm thick tantalum backing with a diameter of 48.0 mm, as shown inFig. 2; d) the back-to-back 0.1 mm thick tantalum backings at the No.5 position, which were 48.0 mm in diameter, were used to measure the background events.

Figure 2.(color online) Pictures of the Ni samples.

The nucleus-number uncertainty was 1.0% for the⁵⁸Ni#I and⁶¹Ni samples, and 5.0% for the⁵⁸Ni#II and⁶⁰Ni samples. Although two⁵⁸Ni samples were used for the measurements, the⁵⁸Ni#II sample had a fairly large nucleus-number uncertainty because a small section was broken during operations. Therefore, the data obtained from the two⁵⁸Ni samples were processed individually and the cross sections obtained for the⁵⁸Ni#II sample were used to check those obtained for the⁵⁸Ni#I sample. The corrections for impurity interference are detailed in Sect. 3.5.

Comparing with our group's previous work [20], a new GIC with symmetric double sections was made and installed whose structure and electronics are shown inFig. 3. The old GIC has been in use for ~ 20 years, and has several shortcomings, such as gas leakage at the valve and signal instability when high voltages are applied to the electrodes. The structure of the new GIC is similar to that of the old one but it has better gas tightness and low signal noise at high voltages. Besides, the new GIC has two fission cathodes (each with the capacity for a²³⁸U sample) to monitor the relative neutron flux, while the old GIC had only one.

Figure 3.(color online) The structure and the electronics of the GIC: FC, fission cathode; A, anode; G, grid; C, cathode; FIFO, fan in–fan out (CAEN N625); PA, charge sensitive preamplifier (MESYTEC MPR-1); TCU, trigger control unit; SA, signal amplifier (ORTEC 572A); DA, delay amplifier (ORTEC 427A); LG, linear gate (ORTEC 542); SC, single channel analyzer (ORTEC 551); PDA, Signatec PDA14 waveform digitizer; GaGe, GaGe OVE-832-007 high resolution PCIe digitizer.

The GIC has seven electrodes: a common-cathode, two grids, two anodes, and two shields. In the present work, to monitor the relative neutron flux, two²³⁸U₃O₈samples were glued onto each shield, respectively. As such, the two shields would be used as the fission cathodes of the fission chambers. The distance was 6.1 cm from the cathode to the grid, 1.4 cm from the grid to the anode, 1.0 cm from the anode to the fission cathode, and 15.4 cm from the cathode to the front surface of the neutron source. Because there is only one foil each for the⁶⁰Ni and⁶¹Ni samples, the GIC would be rotated by 180° during the experiment to measure the forward (0° – 180°) and the backward (90°– 180°) cross sections of the (n,α) reaction. To obtain higher fission tallies, only the signal of the fission cathode near the neutron source was used. For sides 01 and 02, the signals from the cathode and anode were recorded by a PDA14 waveform digitizer that was activated by the trigger control unit. The trigger control unit would produce the external trigger, which was generated by the coincidence of the cathode and the anode signals. The trigger control unit effectively improved the anti-noise and anti-background performance of the GIC. The signal from the fission cathode was recorded by a GaGe OVE-832-007 high-resolution PCIe digitizer.

The GIC's working gas was Xe + 8.5% H₂and the pressure was 0.855 atm so that theα-particles could be stopped before they reached the grids. The high voltages applied to the cathodes and anodes were -1200 and 600 V (the grid electrodes were grounded), which allowed for all the electrons from the ionization tracks to be collected.

With the GIC, information on the energy and angular distribution of the detected charged particle is obtained using the cathode and anode signal amplitudes. If a charged particle is emitted from the sample on the cathode and is stopped by the work gas before it reaches the gird, the cathode signalV_cand the anode signalV_acan be represented by [21]

and

whereG_candG_arespectively are the ratio constants of the cathode and anode that could be determined via energy calibration;Eis the energy of the emitted charged particle;dis the distance between the cathode and the grid; $ \bar{X} $ is the distance from the beginning to the center of gravity of the ionization trace; and $ \delta $ is the grid-inefficiency of the GIC. According to the theoretical calculations, in this work, the grid-inefficiency was 0.0121, which is quite low, meaning that the influence of the grid-inefficiency in the present work was negligible. The charged particle can be identified by determining the valid area in the cathode-anode two-dimensional spectrum, because according to Eqs. (1) and (2), different charged particles with the same energy and emitting angle would induce anode and cathode signals with different amplitudes due to their different ionization traces.

An EJ-309 liquid scintillator detector was used to obtain the neutron spectrum by unfolding the measured pulse-height spectra [22]. The axis of the scintillator detector was along the normal line of the electrodes of the ionization chamber and the 0° direction of the deuteron beam line. The distance from the entrance surface of the scintillator detector to the front surface of the neutron source was 2.60 m.

For each neutron energy point, measurements were performed with the sample sequentially at the No.1, No.2, No.3, No.4, and No.5 sample positions then the GIC was rotated by 180° and measurements were performed at the No.5, No.4, No.3, No.2, and No.1 sample positions. The total neutron beam duration was ~ 70 h.

During the experiment, the fission tallies obtained from the²³⁸U₃O₈sample glued onto the fission cathode near the neutron source were used to monitor the relative neutron flux, and the ratio of the fission count during the foreground measurement to that during the background measurement was approximately 1:1 – 2:1. The scintillator detector would run for ~ 15.0 min every hour to obtain the neutron spectrum.

First, the compoundαsources were used for energy calibration and to determine the valid area ofαevents. Second, the fission tallies from the²³⁸U₃O₈samples were calculated to determine the absolute and relative neutron flux. Third, the netαevents from the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe and⁶¹Ni(n,α)⁵⁸Fe reactions were obtained in the valid area after the subtraction of the background, and the detection efficiency ofαevents was calculated. Fourth, theαand fission events induced by the low-energy neutrons were corrected. Fifth, theαevents from the impurity isotopes were corrected. Finally, the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions were obtained, and the spread of neutron energy was unfolded into a specific energy using the iterative method.

The energy response of the detection system was calibrated using the compoundαsources. The energy of each channel could be calibrated using the fourαenergy groups, as shown inFig. 4. According to Eqs. (1) and (2), for theαparticles, if we let the emitting angles equal 0° and 90°, respectively, the 0° curve and 90° line can be drawn in the cathode-anode two-dimensional spectrum as shown inFig. 4, by which the valid area ofαevents can be determined. Taking the fluctuation in energy resolution into account, the valid area is a little lager than the region between the 0° curve and 90° line as shown by the two red curves inFig. 4.

Figure 4.(color online) The measured cathode-anode two-dimensional spectrum of the compoundαsources in the forward direction.

The total fission tallies from the²³⁸U₃O₈sample were used to determine the absolute neutron flux. A typical anode spectrum of the fission fragments is shown inFig. 5as an example.

Figure 5.(color online) The anode spectrum of the fission fragments for the absolute neutron flux measurement atE_n= 5.50 MeV.

The detection efficiency (ε_f) of the fission fragments for the absolute neutron flux measurement is (75.55 – 84.93)% in the present work. A Monte Carlo simulation was used for the threshold and self-absorption corrections. The simulation code was written using Matlab-2019a [23]. The stopping power of the fission fragments in the samples was calculated by SRIM-2013 [24], the mass distribution of the fission fragments was calculated by the GEF code [25], and the angular distributions of the fission fragments were obtained from Ref. [26]. Details of the simulation can be found in Ref. [19]. The black curve inFig. 5shows the simulation result for the fission fragments, which agrees well with the measurement spectrum. The shape of the simulation spectrum was almost identical at different neutron energy points, because the variation ofE_nwas negligible compared to the total kinetic energy (TKE) of²³⁸U fission fragments (~ 170 MeV). The change ofε_fin the present results was due to the different thresholds. The threshold position was adjusted at different energy points to obtain higher fission tallies while avoiding the background. The total fission tallies could be determined by the fission tallies within the thresholds (N_f) divided byε_f.

As described in Sect. 2.5, the absolute neutron flux was measured as an individual experimental step, while the relative neutron flux was continuously measured for different experimental steps. The relative neutron flux was obtained using the relative fission tallies within the thresholds of the²³⁸U₃O₈sample glued onto the fission cathode for the different experimental steps. For the two fission cathodes, the thresholds were fixed, so the detection efficiency of the fission fragments can be regarded as invariant and be eliminated when calculating the relative tallies.

The cathode-anode two-dimensional spectra of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions were analyzed.Fig. 6(a)shows the two-dimensional spectrum for the measurement of the⁵⁸Ni(n,α)⁵⁵Fe reaction in the forward direction (from the⁵⁸Ni#I sample) after subtracting the background events atE_n= 5.50 MeV. The background tallies were normalized using the fission tallies from the²³⁸U₃O₈sample glued onto the fission cathode. The effectiveαevents could be selected in the valid area described in Sect. 3.1. Then, the two-dimensional spectrum was projected onto the anode channel asFig. 6(b)shows, in which the thresholds were set to separate the effectiveαevents from the background that could not be subtracted. The tally of netαevents within the thresholds ( ${N_\alpha }$ ) can be determined by

Figure 6.(color online) (a) The cathode-anode two-dimensional spectrum after subtracting background events and (b) the anode projection spectrum of the foreground, background, and net events for the measurement of the⁵⁸Ni(n,α)⁵⁵Fe reaction in the forward direction (from the⁵⁸Ni#I sample) atE_n= 5.50 MeV.

where $N_\alpha ^{{\rm{fore}}}$ and $N_\alpha ^{{\rm{back}}}$ are the total tallies within the thresholds of the foreground and background measurements, respectively. ${C^{{\rm{f\_b}}}}$ is the background normalization coefficient, which is the ratio between the fission tallies from the²³⁸U₃O₈sample glued onto the fission cathode during the foreground measurement and those during the background measurement. The detection efficiency of theαparticles was determined from the simulation results of the anode spectrum. Detection efficiency was related to the neutron energy and the threshold of which the details will be described in Sect. 3.4.

The neutron spectrum at each energy point was obtained by unfolding the pulse height spectra measured by an EJ-309 liquid scintillator detector [22]. The neutron energy spectrum atE_n= 5.5 MeV is shown inFig. 7as an example (neutrons with the energy below 1 MeV were ignored because they hardly affect the results). The energy of the neutrons is divided into two parts: the main neutron region around the peak atE_n= 5.50 MeV, and the low-energy neutron region. The variation in the spectrum resulting from the difference between the positions of the scintillator detector and the GIC was corrected using the Monte Carlo method [27]. In the present work, the low-energy neutrons account for 9.29% – 16.50% of the total, and theαevents and fission events they induced should be corrected. $k_\alpha ^{{\rm{low}}}$ and $k_f^{{\rm{low}}}$ are defined as the proportion ofαevents and fission events within the thresholds induced by low-energy neutrons, respectively, and they would be calculated as follows.

Figure 7.(color online) Neutron spectrum atE_n= 5.50 MeV (neutrons with the energy under 1 MeV were ignored).

The angular and energy distributions of theαevents from the (n,α) reaction will change at different neutron energy points, thus changing the detection efficiency; therefore, the anode spectra induced by the neutrons with different energies should be simulated.

For the simulation, the relative intensities of the neutrons at different energy bins were obtained from the measured neutron spectrum, the stopping power ofαparticles in the samples was calculated using SRIM-2013 [24], the angular and energy distributions ofαparticles were calculated using TALYS-1.9 [28] and, for the initial condition, the cross sections were obtained from the ENDF/B-VIII.0 library [9]. Because the tallies of the simulated spectrum may be systematically higher or lower than those measured, the simulated spectrum was multiplied by a factor to fit the measurements, which means that only the variation trend in the cross sections obtained from the ENDF/B-VIII.0 library was used in the present results. The detection efficiencies ofαevents (ε_α) were obtained from the simulated spectrum. According to the simulation,ε_αis 86.73% – 92.64% for the⁵⁸Ni(n,α)⁵⁵Fe reaction in the 4.50 MeV ≤E_n≤ 5.50 MeV region, 69.47% – 87.00% for the⁶⁰Ni(n,α)⁵⁷Fe reaction, and 84.11% – 91.48% for the⁶¹Ni(n,α)⁵⁸Fe reaction in the 5.00 MeV ≤E_n≤ 5.50 MeV region.

Because the cross sections were not measured for low neutron energy in the present work, it is necessary to use the evaluation cross sections to calculate the final results. To estimate the uncertainty of the evaluation cross sections, the ENDF/B-VII.1 library was used, and the final result was recalculated. Then, the difference between the new result and the previous result was regarded as the uncertainty.

As an example,Fig. 8shows theαevents of the⁵⁸Ni(n,α)⁵⁵Fe reaction (in the forward direction atE_n= 5.50 MeV) induced by the total neutrons, main neutrons, low-energy neutrons, and the neutrons with 5.50 MeV of specific energy as “Simulation (total neutrons),” “Simulation (main neutrons),” “Simulation (low-energy neutrons),” and “Simulation ([5.50 MeV] neutrons),” respectively. Using the simulated spectra of “Simulation (low-energy neutrons)” and “Simulation (total neutrons),” the proportion ofαevents within the thresholds induced by the low-energy neutrons ( $k_\alpha ^{{\rm{low}}}$ ) could be determined. The values of $k_\alpha ^{{\rm{low}}}$ are 2.93% – 5.25%, 0.22% – 1.35%, and 2.96% – 4.44% for the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions, respectively.

Figure 8.(color online) The measured and simulated anode projection spectrum for the measurement of the⁵⁸Ni(n,α)⁵⁵Fe reaction in the forward direction (from the⁵⁸Ni#I sample) atE_n= 5.50 MeV. The simulatedαevents were induced by the neutrons in different energy regions.

As described in Sect. 3.2, the detection efficiency of the fission events could be regarded as invariant with the change of neutron energy when the thresholds were fixed; it can then be eliminated when the proportion of fission events induced by the neutrons with different energies is calculated. Therefore, the proportion of the fission events within the thresholds induced by the low-energy neutrons ( $k_f^{{\rm{low}}}$ ) could be determined using the standard cross section of the²³⁸U(n,f) reaction [29] and the measured neutron spectrum. According to calculation, $k_f^{{\rm{low}}}$ is 9.41% – 16.32% at different neutron energy points.

The correction coefficient of the events induced by low-energy neutrons (ρ^low) was introduced by

The values of ${\rho ^{{\rm low}}}$ are 1.064 – 1.132, 1.125 – 1.164, and 1.090 – 1.127 for the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions, respectively. The values of ${\rho ^{{\rm{low}}}}$ are greater than 1, because the low-energy neutrons would have a stronger influence on fission events than onαevents.

As shown inTable 2, in addition to the main isotopes in the samples there are some impurity isotopes, including⁵⁸Ni,⁶⁰Ni,⁶¹Ni,⁶²Ni, and⁶⁴Ni. Theαevents arising from the⁶²Ni and⁶⁴Ni impurity isotopes hardly affect the results because the cross sections of the (n, a) reactions of these two isotopes were smaller than those of other isotopes by two orders of magnitude [15]. Therefore, only the influences of the⁵⁸Ni,⁶⁰Ni, and⁶¹Ni isotopes were taken into consideration.

The interference arising from the impurities was negligible in the⁵⁸Ni(n,α)⁵⁵Fe reaction measurements because the cross sections of the (n,α) reactions of the⁶⁰Ni and⁶¹Ni impurity isotopes were fairly small compared with those in the⁵⁸Ni(n,α)⁵⁵Fe reaction, and the purity of the⁵⁸Ni sample was very high (99.84%). In the⁶⁰Ni(n,α)⁵⁷Fe and⁶¹Ni(n,α)⁵⁸Fe reaction measurements, the interference from the impurities in the samples was non-negligible because the cross sections of the two reactions were approximately 70% – 90% lower than those of the⁵⁸Ni(n,α)⁵⁵Fe reaction and⁵⁸Ni was the main impurity isotope in the⁶⁰Ni and⁶¹Ni samples.

To correct the interference of the impurity isotopes, the anode spectrum contributed by different isotopic compositions was simulated according to their proportions presented inTable 1. A typical example is shown inFig. 9. The conditions for the simulation were similar to those described in Sect. 3.3 but the cross sections of the (n,α) from the impurity isotopes were obtained from the present results: in the measurement of the⁵⁸Ni(n,α)⁵⁵Fe reaction, the interference from the impurities was negligible and consequently the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe reaction could be obtained without correcting for the impurity isotopes; then, for the measurement of the⁶⁰Ni(n,α)⁵⁷Fe reaction, the interference was mainly from the⁵⁸Ni(n,α)⁵⁵Fe reaction, the cross sections of which already had been obtained; lastly, for the measurement of the⁶¹Ni(n,α)⁵⁸Fe reaction, the interference was mainly from the⁵⁸Ni(n,α)⁵⁵Fe and⁶⁰Ni(n,α)⁵⁷Fe reactions, the cross sections of which had also been obtained already.

Figure 9.(color online) The measured and simulated anode projection spectrum for the measurement of the⁶¹Ni(n,α)⁵⁵Fe reaction in the forward direction atE_n= 5.50 MeV. The simulatedαevents were from all isotopes in the⁶¹Ni sample, including the⁶¹Ni isotope and the⁵⁸Ni and⁶⁰Ni impurity isotopes.

The tallies of theαevents from different isotopes could be determined from the simulated results. The proportion ofαevents arising from the impurity isotope (β^impurity) could be determined using the ratio between the simulated tallies within the thresholds of the impurity isotopes and those of all isotopes.β^impurityis ~ 0% for the⁵⁸Ni(n,α)⁵⁵Fe reaction, 0.89% – 1.07% for the⁶⁰Ni(n,α)⁵⁷Fe reaction, and 11.69% – 12.95% for the⁶¹Ni(n,α)⁵⁸Fe reaction.

Either the forward or backward cross sections (σ_α) of the (n,α) reactions can be calculated by

whereσ_fis the standard cross section of the²³⁸U(n, f) reaction, and the values ofσ_fin the present results are 0.5592, 0.5579, 0.5483, 0.5517, and 0.5482 b atE_n= 4.50, 4.75, 5.00, 5.25, and 5.50 MeV, respectively [29].N_UandN_Niare the numbers of²³⁸U and^{58, 60, 61}Ni nuclei in the samples, respectively.N_αandN_fare the tallies of the netαevents and fission events described in Sects. 3.3 and 3.2, respectively.ε_fandε_αare the detection efficiencies for the fission fragments andαparticles at the specific neutron energy described in Sect. 3.2 and 3.3, respectively.ρ^lowis the correction coefficient of the events induced by the low-energy neutrons described in Sect. 3.5.β^impurityis the proportion ofαevents arising from the impurity isotope described in Sect. 3.5.G= 0.9890 – 1.003 is the ratio between the average neutron flux in the area of the nickel sample and that of the²³⁸U₃O₈sample in the sample changer, andGis obtained by the Monte Carlo method (Gis introduced because there is a slight difference between the diameters of the two samples).C^f_fissionis the ratio between the fission tallies from the²³⁸U₃O₈sample glued onto the fission cathode during the foreground measurement and those obtained during the absolute neutron flux measurement.R^unfoldingis the unfolding coefficient that is explained as follows.

According to the measured neutron spectrum shown inFig. 7, in the main neutron region, the width of the neutron peak is non-negligible and the corresponding uncertainty ofE_nis 2.5% – 4.2% in the 4.50 MeV ≤E_n≤ 5.50 MeV region. The spread of neutron energy could be unfolded using the iterative method. $k_\alpha ^{{{\bar E}_{{n}}}}$ and $k_f^{{{\bar E}_{{n}}}}$ , which are the ratios between the tallies of theαand fission events within the thresholds induced by the neutrons with the specific energy ${\bar E_{{n}}}$ and those induced by the main neutrons that were determined for unfolding the neutron spectrum.

$k_\alpha ^{{{\bar E}_{{n}}}}$ can be determined using the simulation method. As shown inFig. 8, theαevents induced by main neutrons and the neutrons with the specific energy ${\bar E_{{n}}}$ can be determined using the simulated tallies of “Simulation (main neutrons)” and “Simulation ([5.50 MeV] neutrons).” The conditions for simulation were the same as those described in Sect. 3.3 except for the cross section of the (n,α) reaction. The forward or backward cross section of the (n,α) reaction was obtained from the last iteration using Eq. (3). In the main neutron region, the neutron energies are very close and the excitation functions of the (n,α) reactions are smooth; therefore, the cross section can be obtained from the measured results using the linear interpolation or extrapolation method. $k_f^{{{\bar E}_{{n}}}}$ could be determined using the standard cross section of the²³⁸U(n,f) reactions [29] and the measured neutron spectrum.

Next, the unfolding coefficientR^unfoldingwas calculated by

Then, the new cross sectionσ_αwould be calculated using Eq. (5). This deconvolution process was iterated until the variation of the cross section was less than 0.1% (usually 10 times).R^unfoldingwas set to 1 as the initial value. In the present work, the final value ofR^unfoldingis 0.9705–1.018, which would slightly correct the original results by -2.95%-1.80%. The value ofR^unfoldingis mainly affected by the average neutron energy in the main neutron region. If the average neutron energy is smaller than the specific energy ${\bar E_{{n}}}$ ,R^unfoldingwould be greater than 1; otherwise, it would be less than 1. The unfolding method could decrease the uncertainty ofE_nfrom 2.5% – 4.2% to 1.4% – 1.5%, and the final uncertainty ofE_nis contributed by the energy resolution of the scintillator detector.

The forward and backward cross sections of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions can be calculated using Eq. (5). The sources of uncertainty and their magnitudes are presented inTable 3. The total (n,α) cross section could be obtained by adding the forward cross section to the backward. Further, the forward/backward ratios in the laboratory reference system could be obtained by division. The results are shown inTables. 4–6andFigs. 10–12(results of the⁵⁸Ni(n,α)⁵⁵Fe reaction were obtained from the⁵⁸Ni#I sample).

Figure 10.(color online) The present⁵⁸Ni(n,α)⁵⁵Fe cross sections (obtained from the⁵⁸Ni#I sample) compared with existing measurements and evaluations and TALYS-1.9 calculations [8,15,28,30].

Figure 11.(color online) The present⁶⁰Ni(n,α)⁵⁷Fe cross sections compared with existing measurements, evaluations and TALYS-1.9 calculations [8,15,28,30].

Figure 12.(color online) The present⁶¹Ni(n,α)⁵⁸Fe cross sections compared with existing measurements and evaluations and TALYS-1.9 calculations [8,15,28,30].

As described in Sect. 2.2, the present cross sections of the⁵⁸Ni(n,α)⁵⁵Fe reaction were obtained from the⁵⁸Ni#I sample, and they were checked using the results obtained from the⁵⁸Ni#II sample, asFig. 13shows. The results obtained from the two samples agreed, thus verifying the reliability of the measurement results .

Figure 13.(color online) The present⁵⁸Ni(n,α)⁵⁵Fe cross sections obtained from the⁵⁸Ni#I and⁵⁸Ni#II samples compared with evaluations and TALYS-1.9 calculations [15,28,30].

TALYS-1.9 [28] was used to calculate the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions. To attain better agreement with the present results, several input parameters, including those of the optic, level density, and stripping models, were adjusted from the default input values, as listed inTable 7(parameters of the stripping model mainly affect the results forE_n> 10.0 MeV). The calculated cross sections and forward/backward ratios agree well with the measured results, as presented inTables 4–6andFigs. 10–12. With the adjusted parameters, the calculated cross sections of other major reaction channels, including the (n,tot), (n,el), (n,n’), and (n,p) reactions, and the angular distributions of elastic scattering have been checked. The calculated results of these reactions also agree well with most existing measurements and evaluations, thus verifying the reliability of the adjusted parameters.

The present cross sections were compared with existing measurement data obtained from the EXFOR library [8] and evaluations obtained from the ENDF/B-VIII.0 [9], ENDF/B-VII.1 [10], JENDL-4.0u+ [11], JEFF-3.3 [12], ROSFOND-2010 [13], CENDL-3.1 [14], and TENDL-2019 [30] libraries:

1) For the⁵⁸Ni(n,α)⁵⁵Fe reaction, the present cross sections in the 4.50 MeV ≤E_n≤ 5.50 MeV region agree well with the measurement data of Fessler (1999, 5.36 – 19.4 MeV) [7], T. Sanami (1998, 4.51 – 6.51 MeV) [2], S. M. Qaim (1984, 5.36 – 9.49 MeV) [3], and the ENDF/B-VII.1 evaluation [10]. Besides, the excitation function calculated using TALYS-1.9 with the adjusted parameters also accords with the ENDF/B-VII.1 evaluation [10].

The present cross sections are 21.6% lower than those in the measurement data of Gledenov (1997, 5.00 – 7.00 MeV) [6] and 79.1% higher than those of Goverdovskiy (1992, 5.12 MeV) [4] in the 4.50 MeV ≤E_n≤ 5.50 MeV region. Compared with other measurement data in this neutron energy region, the Gledenov [6] and Goverdovskiy [4] results respectively have the highest and the lowest values. Uncertainties in their results are greater than 10%, while those in the present results are less than 7%. The V. V. Ketlerov (1996, 3.55 – 6.83 MeV) [5] measurement data suggest that the excitation function has a “valley structure” aroundE_n= 4.92 MeV. However, the present results and none of the evaluations show this structure [9-14,30].

2) For the⁶⁰Ni(n,α)⁵⁷Fe reaction, there is only one measurement (Khromyleva (2018, 6.00 – 7.15 MeV) [16] in the MeV region, which consistent with our measurement results. The present cross sections agree well with the data from the JEFF-3.3 (= ROSFOND-2010) library [12,13]. The excitation function calculated using TALYS-1.9 with the adjusted parameters accords with both measurements and the JEFF-3.3 (= ROSFOND-2010) library [12,13].

3) For the⁶¹Ni(n,α)⁵⁸Fe reaction, there is no measurement in the MeV region and there is a significant deviation among different evaluations. Compared with various other evaluations, the present results better agree with the ENDF/B-VIII.0 library [9]. The excitation function calculated using TALYS-1.9 with the adjusted parameters is close to the ENDF/B-VIII.0 library forE_n< 6.0 MeV [9] and the TENDL-2019 library forE_n> 6.0 MeV [30].

As presented inTables 4,5,and6the nearly symmetrical distributions for the forward and backward directions of the outgoing alpha-particles mean that the compound mechanism is predominant for the⁵⁸Ni(n,α)⁵⁵Fe and⁶¹Ni(n,α)⁵⁸Fe reactions. However, in the case of the⁶⁰Ni(n,α)⁵⁷Fe reaction, a slight non-statistical effect may be available.

Above theE_n= 5.50 MeV region, the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions noticeably grow correspondingly with the increase in neutron energy. The cross sections of the three reactions in the 8.50 MeV ≤E_n≤ 10.50 MeV region are approximately 2 – 10 times larger than those in the 4.50 MeV ≤E_n≤ 5.50 MeV region [9]. However, in the 8.50 MeV ≤E_n≤ 10.50 MeV region, the existing measurement data for the⁵⁸Ni(n,α)⁵⁵Fe reactions are scarce, and there is no measurement for the⁶⁰Ni(n,α)⁵⁷Fe and⁶¹Ni(n,α)⁵⁸Fe reactions [8]. Further, there are significant discrepancies among different evaluations in this neutron energy region [9]. Taking these factors into consideration, measurements of the cross sections of the three reactions are planned in the 8.50 MeV ≤E_n≤ 10.50 MeV region based on the HI-13 tandem accelerator of the China Institute of Atomic Energy (CIAE).

In the present work, the cross sections of the⁵⁸Ni(n,α)⁵⁵Fe,⁶⁰Ni(n,α)⁵⁷Fe, and⁶¹Ni(n,α)⁵⁸Fe reactions were measured in the 4.50 – 5.50 MeV neutron-energy region using the 4.5-MV Van de Graaff accelerator, GIC detector, enriched nickel isotopic foil samples,²³⁸U₃O₈samples, and the EJ-309 liquid scintillator detector. The present data of the⁵⁸Ni(n,α)⁵⁵Fe reaction are consistent with the data of most measurements and the ENDF/B-VII.1 evaluation [8,10]. The present data of the⁶⁰Ni(n,α)⁵⁷Fe reaction are the first measurement results below 6.0 MeV, and those of the⁶¹Ni(n,α)⁵⁸Fe reactions are the first measurement results in the MeV region. The results calculated using TALYS-1.9 with the parameters adjusted agree well with the present data. These results would contribute significantly to clarifying the deviations and discrepancies among the measurements and evaluations to date.

The authors are indebted to the operations crew of the 4.5-MV Van de Graaff accelerator at Peking University. Dr. Qiwen Fan from the China Institute of Atomic Energy is acknowledged and thanked for preparing the⁵⁸Ni,⁶⁰Ni, and⁶¹Ni samples.