Optimising Foil Selection for Neutron Activation Systems

Neutron spectrum unfolding using activation foils is currently the primary technique planned for measuring the neutron energy spectrum at the first wall of power-generating fusion reactors. Room for improvement in the effectiveness of current foil selection was identified, and a program produced to select foils procedurally in order to maximise the accuracy of the unfolding procedure. Using Kullback–Leibler Divergence to quantify the accuracy, the spectrum unfolded by the procedurally selected set of foils is found to be more accurate than the spectrum unfolded by a set of foils used in the literature.


Introduction
As power-generating nuclear fusion devices such as Spherical Tokamak for Energy Production (STEP) and DEMOnstration power plant (DEMO) are being designed, one neutronics challenge looms: our knowledge of the neutron spectrum at the first wall of the fusion reactor is severely limited. This is because the harsh conditions at the first wall forbid the use of traditional direct measurement techniques that can otherwise provide high-resolution neutron-spectrum measurements.
The designs of STEP and DEMO have not been finalized yet, but in general we can expect their neutron flux and temperature to be comparable to, if not exceeding, those of ITER. Even ITER's neutron fluxes (up to 10 14 cm À2 s À1 [1]) and temperature (up to 800 C [2]) will quickly damage most neutron detectors [3,4], and the strong magnetic field will interfere with the operations of electronic components [1].
However knowledge of the first wall spectrum is essential for monitoring the power output of the fusion reactor [5], predicting its tritium breeding ratio [6], inferring the radiation damage to the structural components, validating the radiation transport models for shielding components [7], as well as supplementing the plasma parameter measurements [8]. Many of these activities are critical to the safety, sustainability, and viability of fusion as a commercial energy source.
Neutron activation systems (NAS), which rely on neutron-spectrum unfolding using activation foils, remain the most robust technique currently available for measuring the first wall neutron spectrum, as the count rate can easily be controlled by the gamma-ray counting geometry and foil volume. However, the technique of neutron-spectrum unfolding using activation foils requires optimisation for measuring fusion neutron spectra. Previous investigations have only selected foils on an ad hoc basis (see [5,6,9] for example). In this paper, a programmatic approach to foil selection is shown to increase the precision provided by neutron activation systems.

Theory
By irradiating a set of foils with known composition and counting the radionuclides created as a result, one can infer the neutron spectrum. The number of radionuclides detected in the set of foils forms a vector N with length equal to the number of types of radionuclides produced. N is usually measured using gamma-ray spectrosocpy after extracting the foils from the irradiating environment, and the effectiveness of the gamma-ray spectrosocpy depends on the halflives of the products. The number of each type of radionuclides generated per unit fluence is a product of the number of reactant atoms present in the foil and the microscopic cross-section, the latter of which is dependent on the incident neutron energy. This product is termed the response matrix R. The neutron spectrum is denoted with /. The number of radionuclides present after irradiation can be modelled using the following equation, which when expressed as a summation equation, gives the number of the k th type of radionuclide created by the neutron-induced reactions as: where / i is the fluence in energy bin i, and n is the total number of bins in the spectrum. Each row of R should be sensitive to a different part of the spectrum, such that one can reconstruct the neutron spectrum / from N. This is achieved using an unfolding algorithm, which deduces a possible neutron spectrum / sol. that satisfies Eq. 1. To aid this process, a guessed spectrum, known as the a priori / 0 is typically provided to the unfolding algorithm to start the unfolding. This supplies the algorithm with more (and in the case of underdetermined unfolding, essential) information to increase the accuracy of the unfolded spectrum. Different unfolding algorithms yield different / sol. . Two such algorithms that will be used in this paper are maximum-entropy deconvolution (MAXED [10]) and pseudo-inverse. They are chosen because maximum-entropy and linear least-squares algorithms are the two most widely-used algorithms in NAS (as shown in [11]). MAXED is a well-known algorithm employed in many fusion neutron spectra unfolding investigations, including [12][13][14]. It maximises the the entropy S as a function of / sol. , (where the first term is equal to the negative Kullback-Leibler Divergence ÀD KL ð/ sol. k/ 0 Þ,) while keeping the v 2 of the unfolded spectra at a user-specified value [10].
The pseudo-inverse algorithm is classified as a 'linear least-squares method' [15] and is a special case of STAY'SL [16], where uncertainties on the a priori spectrum and the response matrix are assumed to be negligible, as is done for MAXED by default. The pseudo-inverse algorithm begins at / 0 in the / space and takes a step D/ towards the minimum v 2 manifold: By selecting a step-size \1, the solution's v 2 value can be controlled precisely. Regardless of the unfolding algorithm applied, their ability to correct the a priori will still depend on the quality of the information collected, which is dictated by the response matrix. The response matrix is controlled by the composition of the foils. Therefore we can optimise the response matrix and maximise the effectiveness of unfolding using a foilselector program that chooses the optimal combination of foils. The gist of this foil-selector program is as follows: 1. The user provides a list of candidate materials (a list of size M) that may be safely inserted into the neutron field for activation; • This would decide the isotopic composition of each candidate foil.
2. The user specifies the number of foils (L) that should be present in the foil set recommended by the foilselector; taking into consideration that • increasing the number of foils would improve the unfolding accuracy, • but the foil set would be costly and unwieldy to measure if L is too big.
3. The user specifies the gamma-ray spectrosocpy setup, specifying the • absolute efficiency of the gamma-ray counting setup, • irradiation duration in the neutron field, • transport duration from the neutron field to the detector, • and measurement duration on the gamma-ray detector.
4. The foil-selector generates all 1 possible foil sets by building foil sets of size L out of the M candidate materials.
• It decides the optimal foil geometry (area and thickness) for each chosen material to be formed into to maximize the number of c's detected.
• It then simulates the irradiation, measurement and unfolding of each foil set using the parameters provided by the user in the steps above. • It then evaluates the uncertainties of the unfolded solution spectrum. • It ranks these foil sets by these unfolded spectrum uncertainties.
5. Finally, the user picks the highest-ranking foil set(s) from the program output and purchases foils according to the recommended list of materials and their foil geometry to begin their unfolding experiment.
More details will be discussed in future publications.

Results
Each unfolded spectrum / sol. is an attempt at reconstructing the true spectrum / truth . The highest fidelity reconstruction of / truth is considered the most accurately unfolded spectrum. The following sections compare the accuracy of unfolded spectra when different response matrices are used.

Demonstrating the Importance of Foil Selection with a Toy Model
Using a coarsely binned, hypothetical spectrum, Fig. 1 shows the effect of using different response matrices on the accuracy of the unfolded solution. A neutron spectrum is divided into five equal-sized energy bins. The fluences in bins 1 À 5 are arbitrarily chosen as [10,5,20, 20, 10] a.u. respectively, and set as the 'truth' spectrum (solid blue line). Two sets of foils, set A and set B, are irradiated by this truth spectrum. The number of detectable radionuclides produced through each reaction in response to a neutron in bins 1 to 5 is shown in Table 1. Set A was chosen to contain the same number of types of reactions as set B (i.e., three types of reactions) to ensure a fair comparison, and both sets of foils produce a similar number of radionuclides upon irradiation. The irradiation of foil sets A and B were simulated by calculating the number of radionuclides produced using Eq. 1. The error on the detected number of radionuclides N k was set as rðN k Þ ¼ ffiffiffiffiffi ffi N k p , assuming the counting statistics is the dominant source of error. Using a naïve a priori (a flat spectrum, dash-dotted red line) and an arbitrary small number (0.01) as the target v 2 , the unfolding was performed using MAXED. On Fig. 1 the spectrum unfolded using the radionuclide measurements from set A is plotted in the orange dashed line, while that of set B is plotted in the dotted green line.
The unfolding was repeated using different values of target v 2 , using MAXED as well as pseudo-inverse. The accuracy of the unfolded solutions, as measured using Kullback-Leibler Divergence D KL [17], are shown in Table 2, and are dimensionless. Kullback-Leibler Divergence between two normalized distributions P and Q is defined as Fig. 1 The toy model spectrum unfolded with MAXED using a target v 2 ¼ 0:01 A higher D KL represents a larger deviation of the unfolded solution from the true spectrum. For reference, the a priori itself has D KL ð/ truth ; / 0 Þ ¼ 0:11086. One can conclude from Table 2 that the spectra unfolded using foil set A are more accurate than spectra unfolded using foil set B, irrespective of the target v 2 and the unfolding algorithm used.
The improvement in the unfolded spectra can be explained by the better coverage of the spectrum by set A's response matrix than set B's. Set B's foils are mostly insensitive to neutrons in bins 3-5, thus the deviation of the a priori spectrum from the truth spectrum in bins 3 and 4 are underestimated by the unfolding algorithms. Meanwhile, set A's foils are highly sensitive to bins 3-5, allowing set A to detect the deviation of the a priori spectrum in bins 3-4 and accurately reproduce the neutronspectrum profile in bins 3-5.

Unfolding a Fusion Neutron Spectrum
An analytical method has been developed to select a set of foils that gives a better response matrix for unfolding experiments, such that they produce more accurate unfolded solutions. To demonstrate this in a fusion-relevant scenario, foils were selected for unfolding an expected D-T spectrum of ITER, retrieved from the reference input spectra section of the FISPACT-II Wiki [18] (solid blue line in Fig. 2). To generate a physically plausible a priori spectrum, an approximation is created by fitting a power-law background (/ bg ðEÞ / E À0:910 ) plus a Gaussian peak at 14.079 MeV with a standard deviation of 378.92 keV, corresponding to a purely Ohmic heated D-T plasma with temperature = 25.2 keV [19] (red dash-dotted line, Fig. 2). The target v 2 was again set to 0.01.
A set of foils (set C) was chosen using the foil-selector program and compared against another set of foils (set D) used in a fusion-neutron-spectrum unfolding experiment [12]. To ensure a fair comparison set C was chosen to have the same number of foils (3) as set D. The compositions and detectable products of both foil sets are listed in Tables 3 and 4. The response matrix R is constructed from the microscopic cross-section for the production of each radionuclide, times the number of reactants in the foil available to react with the neutrons. Therefore if multiple production pathways exists via the reaction of multiple reactants j to form product k, then the k th row of the response matrix is formed by Eq. 6: where N D j is the number density of the j th reactant, V is the volume of the foil, and r jk ðE i Þ is the microscopic crosssection of the production of radionuclide k from reactant j at the i th neutron group's mean energy. The microscopic cross-section values are retrieved from the TENDL-2017 nuclear data library [20] and re-binned into the VITAMIN-J 175 group structure to match the group structure of the ITER-DT spectrum. The irradiation of both sets of foils (C and D) were simulated using Eq. 1. Similar to Sect. 3.1 the number of detectable reactants N k has associated standard . In particular, we draw the reader's attention to the difference between set C and set D: set C was chosen such that it has 1. more types of reactions, and 2. more detectable radionuclide half-lives (on the order of 1 h) than set D's foils. Point 1 reduces the degree of underdetermination when unfolding, and point 2 allows a larger fraction of the radionuclides produced to be detected via gamma-ray spectroscopy while minimising the number of pile-up events. A product is considered detectable only if it emits at least one unique gamma-ray line or X-ray line resolvable from the rest of the gamma-ray spectrum using a HPGe detector. For ease of calculation, we assume that 100% of the radionuclides decayed are detected. As demonstrated by Fig. 2, the neutron spectrum unfolded using foil set C fits the underlying spectrum ('truth') better, especially around the D-T peak at % 14 MeV, where the orange dashed line overlaps with the solid blue line consistently. The unfolding procedure is repeated several times using various values of target v 2 and different algorithms (pseudo-inverse as well as MAXED), and the unfolded spectra accuracy is recorded in Table 5. For reference, the a priori spectrum used here has D KL ð/ truth ; / 0 Þ ¼ 0:0861. One can conclude from Table 5 that the spectrum unfolded using foil set C is more accurate than the spectrum unfolded using foil set D, even after controlling for the target v 2 and the unfolding algorithm used.
A pertinent question is why both sets of foils' unfolded spectra fit the high energy end of the truth spectrum well but fit the low energy end of the truth spectrum poorly. This can be explained by examining the response matrices for sets C and D, or by taking a closer look at Tables 3 and  4. Most neutron-induced reactions are either capture reactions (n,c) or threshold reactions. In fact, most of the reactions that form the response matrices are threshold reactions [(e.g., (n,p), (n,2n), etc.], mostly or completely insensitive to neutrons that enter the foil below its threshold energy. Most nuclear reactions included in these matrices have threshold energy [ 7:5 MeV; thus it is almost impossible to find any reaction with non-negligible microscopic cross-sections at \5 MeV in this response matrix. Therefore the only reactions that are sensitive to the low energy part of the spectrum in Fig. 2 are the seven capture reactions in set C and two capture reactions in set D. These few capture reactions take on the mammoth task of correcting the neutron fluence in over a hundred bins 2 . It is, therefore, no surprise that they fail to adjust the neutron spectrum satisfactorily at low energies. The limited variety of reactions is a difficult issue to solve if we continue to restrict to using very few (three) elements. On the other hand, the large number of threshold reactions in set C, in addition to the coarse binning of the VITAMIN-J 175 group structure at high energy ([ 5 MeV), reduces the effective degree of underdetermination on the high energy side of the spectrum, leading to excellent accuracy of set C's unfolded spectrum around the D-T peak. It is possible that a different choice of group structure will yield an even better selection of foils and lower unfolded D KL values, but group structure optimisation is beyond the scope of this paper, so the analysis above adheres to the VITAMIN-J 175 group structure commonly used in fusion neutronics [5,21,22].

Conclusion
Better unfolding accuracy was achieved using a set of activation foils selected using an automated foil selection procedure. This is demonstrated by the lower D KL values for the respective set of procedurally selected foils in Tables 1 and 5, irrespective of other unfolding parameters (target v 2 and choice of unfolding algorithms). The overall merit of the procedurally-selected set C in Sect. 3.2 is clear: it unfolds to a more accurate spectrum than if set D were used, especially around the D-T peak. This means it can give reliable information about how many neutrons are created from the plasma through the D-T reaction, which is useful for deducing the power output of the reactor. However, it cannot provide reliable information about how neutrons scatter to lower energies and thermalise within the tokamak. This issue may be solved by adding more foils with threshold reactions of lower threshold energies into the foil set in the future.
This research allows existing and future neutron activation systems to improve their neutron measurement capability by recommending the best combination of foil compositions. Data Availability All data generated or analysed during this study are included in this published article.

Conflict of interest None declared.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.