Next-Generation Neutron Detection Using a 6Li Glass Scintillator Composite

doi:10.21203/rs.3.rs-4060006/v1

Download PDF

Article

Next-Generation Neutron Detection Using a ⁶Li Glass Scintillator Composite

https://doi.org/10.21203/rs.3.rs-4060006/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Neutron detectors play a crucial role in fundamental science, nuclear security, safeguards, and civil applications. Traditionally, ³He-filled gas proportional counters are the gold standard for thermal neutron detection, prized for their efficiency, neutron/gamma discrimination, and stability—however, the scarcity of ³He has prompted the search for alternatives. Here, we introduce an innovative neutron detector design based on a scintillating composite consisting of ⁶Li-enriched glass scintillator particles dispersed in an organic matrix. A detector consisting of this scintillating composite, photomultiplier tubes (PMTs) for optical detection, and electronics for reading out the PMT signal in both pulse and current modes was prototyped and characterized using various neutron and gamma sources. The prototype (69.7 mm diameter and 260 mm length) achieves a detection efficiency of 6.7%, a die-away time of 10.5 µs, a negligible gamma misidentification probability (for ¹³⁷Cs and ⁶⁰Co gamma sources), and response linearity up to at least 3.7×10⁶ incident neutrons/s established via an innovative cross-calibration technique. This novel detector holds the potential to outperform traditional ³He-gas-based neutron systems, offering a viable alternative amidst the ongoing ³He shortage and promising advancements in neutron detection technology.

LA-UR-24-20569

Physical sciences/Physics/Applied physics

Physical sciences/Materials science/Materials for devices/Sensors and biosensors

Physical sciences/Optics and photonics/Applied optics/Optical sensors

Detecting neutrons amidst ubiquitous gamma-ray backgrounds is crucial in various scientific fields, including nuclear physics, neutron scattering science, fusion research, nuclear safeguards, security, nuclear engineering, and industrial well-logging applications [1–6]. Because of their neutral charge, neutrons do not produce ionization when passing through a material and instead interact primarily with atomic nuclei through scattering or nuclear reactions. As a result, direct detection via their induced charge, as is the case for alpha or beta particles, is not feasible for neutrons. Instead, neutron detection hinges on ionization by secondary particles, such as recoiling nuclei (protons) or charged reaction products. The principal mechanism for neutron detection also depends on the energy of the incident neutrons. “Fast” neutrons, which have energies on the keV to MeV scale, interact primarily through scattering on atomic nuclei, and the ionization induced by the recoiling nucleus can be used as the means for detection. Materials containing hydrogen are preferred for this application. Neutrons produced in nuclear fission are an example of fast neutrons. For lower energy “slow” neutrons (less than a few eV), most detection methods rely on nuclear reactions that become more probable at low energy, such as neutron capture. Some key isotopes, such as ³He, ⁶Li, and ¹⁰B, have very high neutron capture cross-sections for thermal neutrons (~ 25 meV). These reactions release a substantial amount of energy (hundreds of keV to a few MeV) in the form of charged particles, making them relatively easy to detect [2–5].

Neutron detectors based on ³He gas proportional counters have traditionally been considered the gold standard for thermal neutron detection. In such detectors, an incident neutron, usually first thermalized in a moderator material such as polyethylene, interacts with the ³He gas via the ³He(n,p)³H neutron capture reaction, which has a large thermal neutron cross-section (5,330 barns [4]). The charge impulses induced by the p and ³H interactions into the gas are then converted to logic pulses. Gamma rays, on the other hand, interact only weakly with the ³He gas. When the voltage pulses are viewed in a pulse-height spectrum, neutron capture produces significantly greater pulse heights than gamma interactions, and neutrons can thus be detected reliably as events having a pulse height greater than a pre-set threshold. These detectors offer excellent performance, such as high neutron detection efficiency, effective gamma-ray/neutron discrimination, and long-term stability, making them the most widely deployed type of neutron detector [4, 5]. According to A. Cho [7] and F. Sacchetti et al. [8], 84.5% of the usage of ³He gas in neutron detectors is dedicated to security applications, followed by neutron scattering and scientific applications (10%).

The potential for missing special nuclear material (SNM) in the civil nuclear fuel cycle is a significant concern for the international community. To address this concern, portal monitors, portable neutron detectors, and neutron well counters are commonly used detector systems in nuclear security, non-proliferation, and safeguards to detect potential nuclear material diversion for clandestine purposes. Another common application of ³He detectors is neutron multiplicity measurements for nuclear material control and accounting (NMC&A). Conventional passive neutron multiplicity counting (PNMC), using a thermal neutron well-counter, is a well-established method for assaying special nuclear materials like plutonium. In this technique, an item is placed inside a well-counter comprised of numerous (dozens to more than a hundred) ³He proportional counters embedded in high-density polyethylene, with ³He gas pressure ranging from 4 to 10 atm. The charge impulses collected from the ³He proportional detectors are converted into logic pulses and summed together to create the pulse train. In PNMC, the first three factorial moments of the neutron detected triggered count-rate histogram are extracted from the pulse train. They are called Singles, Doubles (correlated pairs), and Triples (correlated triplets); they are related to the properties of the special nuclear material under assay. The conventional ³He-based design has enabled the construction of thermal-neutron well counters with high neutron efficiency (up to about 50–60% for fission neutrons) for an item located in the center of the cavity. It is worth noting that in Doubles counting, the efficiency increases with the square of the efficiency of the detector, and in Triples counting, it increases with the cube of the efficiency, making the efficiency a demanding property of a multiplicity detector system. The ³He-based systems offer a relatively short detector time response, so-called die-away (some 10s of µs), as a short die-away time allows to minimize accidental coincidences [5, 9, 10]. Neutron-gamma separation, efficiency, short die-away, and long-term stability are critical requirements for this application. In recent decades, the limited supply and high cost of ³He have spurred extensive research into alternative neutron detection technologies. Exploring alternative techniques for neutron multiplicity counting – one of the more demanding neutron detection techniques described above – presents unique challenges [11–18].

The present study focuses on developing a novel neutron detector, offering a viable alternative to ³He-based neutron detectors. The innovative detector presented here is based on scintillating composites, wherein ⁶Li-enriched glass scintillator particles serve as neutron detection centers, effectively addressing the sensitivity to gamma rays exhibited by the ⁶Li-enriched glass scintillator. We focus is primarily on nuclear safeguards and security applications, such as neutron multiplicity counting. Similar to ³He, ⁶Li has a sizeable thermal neutron capture cross-section (940 barns) via the ⁶Li(n,α)³H reaction (reaction Q-value = 4.78 MeV, ³H at 2.73 MeV, and α at 2.05 MeV) [19]. In contrast to ³He, the ⁶Li⁺ ion can be a high atom density constituent of solid scintillators, allowing sizeable neutron capture cross-sections to be realized. Furthermore, ⁶Li is much more readily sourced than ³He and can thus be a potential alternative for efficient neutron capture. The neutron detector developed in this work utilizes scintillating composites consisting of small ⁶Li glass scintillator cubes (GS20^™, Scintacor, Cambridge, UK; 95% ⁶Li enriched, 6.6 wt% ⁶Li [20]) of 1.5 mm × 1.5 mm × 1.5 mm in size that are dispersed in a non-scintillating optically transparent mineral oil matrix at a density of 2.8 vol% [17]. Incident fast neutrons, such as fission neutrons, are moderated to thermal energies within the organic matrix and subsequently captured by ⁶Li in the ⁶Li glass scintillator cubes, where they convert to charged particles through the ⁶Li(n,α)³H reaction. These reaction products deposit their energy in the ⁶Li glass matrix, a portion of which excites the Ce³⁺ dopant, which subsequently decays to the ground state via the 5d→4f optical transition around 395 nm (57 ns lifetime). Through this process, ~ 6,000 optical photons are emitted isotropically per neutron capture event [17, 20]. This scintillation light then propagates through the optically transparent mineral oil matrix to photodetectors that are optically coupled to the composite. In contrast, the secondary electrons ejected by gamma-ray interactions within the ⁶Li glass scintillator and the mineral oil only deposit relatively small amounts of energy in the small scintillator particles, most of these electrons passing through the scintillator cube, producing significantly fewer photons than the scintillation from a single neutron capture event. As we will show, this enables excellent neutron-gamma discrimination via simple pulse-height discrimination akin to ³He gas proportional counters (see Fig. 1) [16–18].

In the following sections, we present a comprehensive analysis of the neutron detection properties of the ⁶Li glass composite scintillator detector described above. This novel detector system showcases excellent capabilities: excellent neutron-gamma discrimination, a short die-away time, and a detection efficiency that is comparable or superior to traditional ³He detectors embedded in high-density polyethylene moderators are demonstrated. These properties are particularly favorable for the future construction of thermal neutron multiplicity counters with superior performance and attractive cost compared to the current ³He-based neutron multiplicity counting systems. Additionally, we assess the performance of this ⁶Li glass composite scintillator detector under exposure to a high-intensity ²⁵²Cf spontaneous fission neutron source. The results indicate that the ⁶Li glass composite scintillator detector significantly surpasses the performance of the ³He-based detector system, presenting an attractive alternative for a range of neutron detection applications in which high neutron fluxes (continuous or pulsed) are present.

The structure of the paper is as follows: design and fabrication of the detector, results of performance evaluations, followed by discussion of the applicability to multiplicity neutron detection. Finally, the paper concludes by reporting on the future directions of the research.

Design and Detector Fabrication

The fabrication process of the scintillating composite neutron detector is described in the Materials and Methods section [17]. The detector consisted of 8200 cubes of GS20™ (Scintacor, Cambridge, UK) ⁶Li glass scintillator, each 1.5 mm × 1.5 mm × 1.5 mm in size and attached to fused silica rods. These assemblies were then immersed in mineral oil inside a 69.7 mm diameter by 260 mm long fused silica tube. Hamamatsu R6233 photomultiplier tubes were glued to each end of the fused silica tube. The fused silica tube was wrapped with an optical reflector and finally with a 0.5 mm thick sheet of cadmium. Figure 2 shows a picture of the detector.

Detector Performance Evaluation

Pulse-Height Spectra

The detector was exposed to a 10.5 µCi ²⁵²Cf fast fission neutron source, as well as to ¹³⁷Cs, ⁶⁰Co, and ²³²Th gamma sources to determine the gamma/neutron separation in the pulse-height spectrum and the associated gamma misidentification probability. Figure 3a shows the pulse-height spectrum for exposure to the ²⁵²Cf source. The pulse-height distribution comprises a peak centered around channel 1560 that is due to neutron capture events on ⁶Li followed by Ce³⁺ scintillation emission as well as a broad distribution at lower channel numbers that is due to gamma-ray interactions and some neutron capture events that did not deposit their full energy because the reaction products escaped the scintillator particle. A Gaussian fit to the neutron capture peak (solid orange curve in Fig. 3a) indicates an energy resolution (defined as full width at half max (FWHM) relative to the peak centroid) of 15.3% with a centroid at channel 1560. The intersection of the Gaussian fit to the neutron capture peak and the gamma distribution was determined to be at channel 1310. As a result, the region of interest (ROI), defining neutron capture events, was set from channel 1310 to 2000 (the highest channel recorded in the multi-channel analyzer) for the measurements.

Figure 3b shows the pulse-height spectra due to pure gamma-ray sources in comparison with the neutron distribution measured above. The emission energies of the gamma source used here are 662 keV (¹³⁷Cs), 1173 and 1333 keV (⁶⁰Co), and 2614 keV (²³²Th). Notably, ¹³⁷Cs is a significant contributor to the gamma emission from irradiated nuclear fuel, given its half-life of approximately 30 years [5, 21].

The excellent separation between neutron and gamma events observed in Fig. 3 is a direct result of the structure of the ⁶Li glass composite scintillator. The ranges of the ³H and α produced by a neutron capture reaction on ⁶Li are small compared to the ⁶Li glass scintillator cube dimensions. On average, ³H and α deposit nearly all their energy into the glass host, where part of it is converted to scintillation photons. In contrast, gamma rays only deposit a small fraction of their energy into the ⁶Li glass scintillator cubes, resulting in a significantly smaller number of scintillation photons and, thus, a smaller detected pulse height when compared to neutron capture events. Neutrons could, therefore, be detected with high confidence by setting a pulse height threshold at channel 1310 (see Fig. 3a) and counting events with pulse heights greater than this threshold.

To evaluate the probability of misidentifying a gamma event as a neutron event, we define the gamma misidentification probability as the fraction of the gamma distribution present in the ROI chosen for the neutron capture peak (channels 1310–2000), which is indicated by the dashed blue vertical line in Fig. 3b. Although the pulse-height spectra were corrected for background, the statistics achieved during the acquisition time in this measurement was insufficient for the ¹³⁷Cs and ⁶⁰Co spectra to converge to their actual values in the ROI. Instead, the Compton tails of the ¹³⁷Cs, ⁶⁰Co, and ²³²Th spectra were fitted by double-exponential functions to obtain better representations of the respective gamma contributions in the ROI. From these fits, we find a gamma misidentification probability of 9.41×10^− 3 for ²³²Th. In contrast, for ¹³⁷Cs and ⁶⁰Co, the Compton edges drop off well below the ROI, and the respective gamma misidentification probabilities are zero.

Note that, given the position and width of the neutron capture peak relative to the gamma distribution in the pulse-height spectrum, there is ample room to further reduce the gamma misidentification probability by increasing the pulse-height threshold, i.e. increasing the lower channel bound of the ROI. This would be particularly effective for ²³²Th, which emits several gamma rays with high energy. The tradeoff of increasing the pulse-height threshold is that some neutron capture events will then not be counted. The choice depends on the needs of the application. The gamma misidentification probabilities given above should, therefore, be viewed as upper limits.

Neutron/Gamma Discrimination Stability to Variation of High Voltages

Measurements were conducted to evaluate the sensitivity of the neutron counts to variations in the high voltages (HV) applied to the PMTs. For this purpose, the variation of the combined ROI count rate was measured for different PMT A and PMT B high voltage settings with the detector exposed uniformly to a ²⁵²Cf neutron source. The individual ROI count rates were initially matched by setting PMT A to 1000 V and adjusting PMT B to 983 V, and then the resulting combined ROI count rate served as the reference (shown in bold in Table 1). The PMT voltages were then varied ± 20 V from the gain-matched condition to assess the resulting variation of the combined ROI count rate relative to the reference.

Table 1

Sensitivity of the detector count rate to the variation of the high voltages applied to the PMTs, where 1000 V (PMT A) and 983 V (PMT B) were the nominal operating voltages.
		PMT B
		963 V	983 V	1003 V
PMT A	980 V	96.96%	99.12%	99.86%
	1000 V	99.59%	100.00%	100.49%
	1020 V	99.44%	100.34%	100.17%

The results in Table 1 show that the neutron count rate varies by only a fraction of a percent when adjusting the PMT high voltages by ± 20 V. This stability of the neutron count rate in response to high-voltage changes demonstrates the reliability of the ⁶Li glass composite scintillator detector’s neutron counting performance. This characteristic is crucial to ensure accurate measurements over extended durations. Furthermore, this behavior resembles that of ³He proportional neutron detectors operating in the so-called high voltage plateau, where the neutron count rate remains practically independent of the HV variation (within a limited specific range of operating voltage) [5].

Cross Calibration and Dynamic Range

Figures 4a and 4b show the combined count rate and combined current from both PMTs of the ⁶Li glass composite scintillator detector, respectively, measured at various neutron fluxes. The variation in the neutron flux was achieved by positioning the detector at various distances ranging from 50 to 230 cm from the ²⁵²Cf source (emitting 5.2×10⁸ n/s) located at the neutron-free-in-air (NFIA) facility at the Central Health Physics Calibration Facility (CHPCF at Los Alamos National Laboratory, TA-36). The neutron count rate was found to depend linearly on the neutron flux up to ~ 2.6×10³ n/(s⋅cm²), corresponding to a count rate of ~ 3×10⁴ cps (see Fig. 4a). Deviation from linear behavior was observed for greater neutron flux, as the multi-channel analyzer (Ortec 927 ASPEC) became increasingly saturated for dead time fractions greater than 25%. In contrast, the PMT current depended linearly on all neutron flux measured in this experiment (see Fig. 4b).

The range of neutron fluxes < 2.6×10³ n/(s⋅cm²) defines a cross-calibration regime in which the neutron count rate is not yet saturated (Fig. 4a) and a PMT current can be measured (Fig. 4b). Using this data, Fig. 5a plots the PMT current vs the neutron count rate in this cross-calibration regime. The linear fit represents an approximation for the charge per neutron capture event produced by the PMTs, allowing the current to be cross-calibrated to the count rate. A cross-calibration constant of 78.5 pC per neutron capture was found, a value characteristic of this particular detector configuration. Note that this approach tends to overestimate the charge per neutron capture event because the PMT current results from all interactions that generate scintillation light (neutron and gamma), while the count rate in the ROI is almost exclusively due to neutron capture events. Being cognizant of this limitation, the cross-calibration constant was then used to convert the current measurements into neutron count rates, revealing the wide dynamic range of the ⁶Li glass composite scintillator detector. As shown in Fig. 5b, the detector response in current mode was found to be linear over a dynamic range extending from single event neutron counting at background levels up to the maximum total incident neutron flux of 3.7×10⁶ n/s that was achievable at the NFIA facility.

Neutron Detection Efficiency

The neutron detection efficiency is defined as the number of neutrons detected in the ROI per the number of neutrons incident upon the detector. With reference to fission (fast) neutrons, it was evaluated both experimentally (²⁵²Cf fission source) and compared with Monte Carlo simulations. The slope of the linear fit in Fig. 5b corresponds to the experimentally determined neutron detection efficiency, and a value of 6.7% was found.

Additionally, MCNP 6.2 Monte Carlo simulations, with ENDF VIII library, were performed, assuming the NIST-calibrated ²⁵²Cf neutron source was located 2.3 meters from the center of the detector, perpendicular to the detector axis [22, 23]. In the MCNP code, the ⁶Li glass particle / mineral oil composite prototype simulated with the same detector geometry as described in the Fabrication Section. The detector model was integrated into a high-fidelity validated model of the NFIA facility (see Methods section). Simulations of the experiment were performed using an isotropic ²⁵²Cf source and the full NFIA facility geometry to accurately capture the effects of scattering in the room. The detector efficiency was calculated to be 7.3 ± 0.2%.

Die-Away Time

The die-away time of the detector was measured using the Rossi-α method [5, 9], which captures the time distribution of neutron events following every detected neutron on the pulse train chosen as the initiating neutron event. The Rossi-α distribution is constructed by recording the time interval between subsequent neutron events within a specified time window. The process involves moving to the next neutron event trigger and repeating the acquisition of the arrival time distribution [5]. Figure 6a shows the detector die-away time response measured using the Rossi-α method. The constant background is due to uncorrelated neutron detection events. The transient exhibits a fast component at early times (t < 8 µs) that is attributed to neutrons still in the moderation process. At longer times (t > 8 µs), the behavior is dominated by thermal neutron diffusion, and the decay is single exponential. The dashed line is a single exponential fit to the data in the 8–60 µs range, yielding a decay time of 10.5 ± 0.1 µs.

As for the efficiency, MCNP Monte Carlo simulations were performed to estimate the neutron die-away time [23], and the results are presented in Fig. 6b. The time distribution was fitted with a simple exponential model within the 8–60 µs range to determine the detector’s die-away time, and a value of 10.21 ± 0.16 µs was found. This aligns closely with the experimental value of 10.5 ± 0.1 µs.

The die-away time is dominated by thermal neutron diffusion through the mineral oil moderator to the ⁶Li glass scintillator cubes. The high density of ⁶Li neutron capture centers throughout the moderator leads to a shorter average diffusion path compared to a detector consisting of ³He tubes encased in polyethylene. Thus, the die-away time in the composite detector is shorter than the several tens of µs die-away time that is typical of ³He proportional counter-based systems embedded in polyethylene. This characteristic is crucial as it minimizes the impact of accidental coincidence in neutron coincidence/multiplicity type detector systems (see Applicability to Neutron Multiplicity Counting section of the paper) [9, 10].

Position Sensitivity of the Neutron Flux Interaction

Scintillation light experiences some attenuation as it propagates through the composite. If a neutron capture event occurs in the middle between the two PMTs, the resulting scintillation light is attenuated equally as it propagates the same distance, on average, to each of the PMTs, and the pulse heights registered by the PMTs are, on average, equal. If, however, the neutron capture event occurs closer to PMT A than to PMT B, then the resulting scintillation light is attenuated less on its shorter path to PMT A, and the pulse height registered by PMT A is thus greater than that registered by PMT B. This effect can be utilized to detect the neutron capture event's location along the composite cylinder's long axis.

Figure 7 shows the change in the pulse-height spectra recorded by PMT A (red curves) and PMT B (green curves) for a ²⁵²Cf source placed at five different positions along the long axis of the composite cylinder, with $z=0,$ indicating the center of the composite cylinder. The pulse-height spectra recorded for uniform illumination of the detector are also shown for reference. Note that the pulse-height spectra shown here are for the individual PMTs, and the resolution is therefore significantly worse than that shown in Figs. 1 and 3, which used the combined PMT A and PMT B signal. As expected, the individual pulse-height spectra of PMT A and PMT B for the symmetrical case of $z=0$ were found to be identical as the scintillation light propagation losses from the center to PMT A and from the center to PMT B are identical. A symmetric case was also realized for uniform neutron illumination, resulting in identical pulse-height spectra for PMT A and PMT B. In contrast, the pulse-height spectra deviated as the source was moved away from the center of the composite cylinder, and the respective light propagation losses became different. As expected, the changes were found to be anti-symmetrical in that the pulse-height spectra for $z=\pm 5 \text{cm}$ and for $z=\pm 11 \text{cm}$ they were identical except for PMT A and PMT B being reversed.

One way to quantify the dissimilarity of the pulse-height spectra $A$ and $B$ is by defining the dimensionless quantity $\delta$ as

$$\delta =\text{sgn}\left[\text{max}\left(A\right)-\text{max}\left(B\right)\right]\times \frac{\sum _{i=1}^{N}{\left({A}_{i}-{B}_{i}\right)}^{2}}{\sum _{i=1}^{N}{\left({A}_{i}+{B}_{i}\right)}^{2}}$$

where $A$ and $B$ are the vectors representing the pulse-height spectra of PMT A and PMT B, respectively, $N$ is the number of elements in the pulse-height spectrum vectors, $\text{sgn}$ is the sign function

$$\text{sgn}\left(x\right)≔ \left\{\begin{array}{c}-1 \text{if} x<0\\ 0 \text{if} x=0\\ 1 \text{if} x>0\end{array}\right.$$

and $\text{max}\left(x\right)$ returns the greatest element within vector $x$. That is, the second factor is the normalized difference in the pulse-height spectral shapes, and the first factor signifies the sign of the spectral difference. As shown in Fig. 8, a linear relationship between $\delta$ and the source position was found, and the respective slope was found to be 27.8 ± 0.8 cm for this experimental configuration. Given two pulse-height spectra $A$ and $B$, this calibration constant can now be used to determine the position of the source along the long axis of the composite cylinder. This position sensitivity may be attractive for applications requiring precise spatial information in neutron detection, such as nuclear security (e.g., verification at ports) [7] and neutron radiography. Note, however, that this method does not provide the location of individual neutron capture events but rather requires acquiring many neutron capture events to obtain the necessary pulse-height spectra for the analysis.

Applicability to Neutron Multiplicity Counting

The performance of neutron multiplicity counting relies on crucial factors, including neutron detection efficiency, die-away time, sensitivity to gamma rays, count rate capability, and counting stability during extended measurements. To assess the performance of the present ⁶Li-glass particle detector for a proposed use in neutron multiplicity counting, we conducted initial Monte Carlo simulations using the MCNP 6.2 code. The results of the previous sections demonstrate that MCNP simulations can effectively replicate the detector's performance.

The MCNP simulations involved filling a counter with ⁶Li-glass-composite material, mimicking the form factor of the widely used Active Well Coincidence Counter (AWCC) [5] – a typical system used in nuclear safeguards based on ³He counters embedded in high-density polyethylene. The modeling assumed uniformly distributed 1 mm³ lithium glass particles with a 4 mm pitch. While light transport was not considered, the calculation estimated the neutron reaction rate. Additionally, different from the original AWCC based on ³He, detection material was filled in the well's top and bottom end plugs, enabling a detector solid angle of approximately 4π. This is a possibility which was impracticable with the traditional ³He-based neutron counter.

The results of the Monte Carlo simulations revealed impressive performance improvements for the AWCC-⁶Li-composite-based detector. It exhibited approximately 2.5 times higher efficiency (ε) and around 3.6 times shorter die-away time (τ) compared to the ³He-based AWCC. This achievement brought the performance of the ⁶Li-glass-particle-based AWCC on par with the exemplary characteristics of the Epithermal Neutron Multiplicity Counter (ENMC) [24] but with a significantly reduced footprint. The performance metric, often expressed as Figure-of-Merit (FOM) for multiplicity neutron detectors, is represented by the equation FOM = ε³ / τ. A higher multiplicity FOM corresponds to better multiplicity detector performance [11, 14]. Results are presented in Table 2. For reference, the ENMC boasts a FOM of 12,500, while the simulated ⁶Li-glass-particle-based AWCC demonstrated an impressive theoretical FOM of 20,500. Even when considering the loss due to light transport, the simulation results suggest that our design outperforms ³He-based multiplicity systems. In addition, the performance calculated so far assumes a “point” item at the center of the cavity. Note that the position-sensitive characteristic of the ⁶Li composite detector may provide further information related to the distribution of the nuclear material in the item under analysis. These findings promise to advance neutron multiplicity counting capabilities in various scientific and nuclear security applications.

Table 2

Comparison of neutron multiplicity counter performances between ³He- based and ⁶Li particle-based AWCCCs.
Characteristic	AWCC (³He)	ENMC (³He)	AWCC (⁶Li particle-based 4-mm pitch)
Outer dimension (cm)	48 (diameter)	66 (per side, box geometry)	48 (diameter)
Weight (estimated) (kg)	84	280	122
Neutron detection efficiency (fission/fast neutrons) with item in the center of the cavity	26%	65%	66%
Die-away time (µs)	51	22	14
Multiplicity FOM (ε³/τ)	345	12,500	20,500

Future Directions

The detection of shielded nuclear material in nonproliferation and safeguards applications has seen advancements through active interrogation techniques, which involve using intense neutron pulses to interrogate the presence of nuclear material, particularly uranium items [25]. Recent developments in laser-driven pulse neutron sources have showcased the production of highly intense, pulse-forward neutron beams, opening the door to a new generation of active interrogation systems [26, 27]. Neutron yields on the order of more than 10¹⁰ neutrons per sr per pulse have been demonstrated with about 80 J laser energy, in a conversion two-step approach where laser light impinging on a suitable target produced a burst of ions that are converted in material such as Be or Li. A pulse of fast neutrons with a distribution of several MeV is produced in a few ns. These laser-driven neutron sources for applications in global security systems demand high-performance neutron detectors in terms of efficiency and dynamic range to measure induced neutrons [26–29]. In this context, the ⁶Li-glass particle detector offers advantages over traditional ³He detector systems for such applications.

The compactness and modularity of the present ⁶Li glass particle detection system are essential features. They are particularly valuable for nuclear security and emergency applications, where portability is crucial. In such cases, Silicon Photomultipliers (SiPMs) can be utilized instead of Photomultiplier Tubes (PMTs) to ensure even greater compactness with a weight reduction. In addition, the form factor (e.g., square size for the detectors) can be selected to facilitate modularity and portability, where several detector modules can be easily arranged to assemble a coincidence/multiplicity counter directly on the field.

Our novel neutron detection technology holds significant promise for enhancing active interrogation capabilities and addressing nuclear security challenges.

Detector Fabrication

A GS20™ ⁶Li glass scintillator plate (Scintacor, Cambridge, UK), with dimensions of 158 mm × 158 mm × 1.5 mm and polished to an optical finish on each side, was diced into 1.5 mm × 1.5 mm × 1.5 mm cubes using a diamond cutting wheel. The resulting cubes had two opposing polished faces and four as-diced rougher faces. Next, UV-curable NOA88 optical adhesive (Norland Products, Cranbury, NJ) was used to glue the ⁶Li glass scintillator cubes to the outer surfaces of fused silica tubes having a cross-section of 2.7 mm × 2.7 mm, a wall thickness of 0.35 mm, and a length of 256.8 mm (Technical Glass Products, Inc., Painesville Twp., OH). 50 scintillator cubes were affixed with a center-to-center spacing of 4.7 mm to each lateral side of the square tubes, using the as-diced sides of the cubes as a gluing interface. The linear array was offset by 2.35 mm on adjacent sides of the square tube. A total of 8200 GS20 scintillator cubes were used to fabricate 41 such scintillating tube assemblies. A hexagonal pattern (7.24 mm lattice constant) of 45 holes (4.0 mm diameter) was drilled into two 4.0 mm thick and 65 mm diameter disks of Solacryl® SUVT acrylic (Spartech Polycast, Maryland Heights, MO). These two disks and four-square fused silica tubes (without GS20 cubes) were used to create a structure to support the 41 scintillating tube assemblies that were inserted through the holes in the acrylic disks.

This assembly was then inserted into a fused silica tube of 69.7 mm internal diameter, 2.25 mm wall thickness, and 260 mm length (Technical Glass Products, Inc., Painesville Twp., OH), which also had two tubular ports at each end of the lateral surface. This assembly defines the active volume of the prototype built.

Next, the windows of two Hamamatsu R6233 photomultiplier tubes were glued to each open end of the fused silica tube containing the scintillator rod assembly. This produced a liquid-tight cylindrical container that was then filled with Drakeol® 9LT (Calumet Specialty Product Partners, Karns City, PA) mineral oil through the two tubular ports, thus completely submersing the scintillating tube assembly in mineral oil. After sealing the tubular ports, a 250-µm diameter silica optical fiber was spiral-wound (7 windings) around the length of the fused silica tube to create an air gap for the subsequent wrapping with Vikuiti™ Enhanced Specular Reflector (ESR) film (3M Company, Maplewood, MN). The air gap ensured that scintillation light emitted within the assembly and incident on the lateral surface of the fused silica tube with an angle greater than the angle of total internal reflectance (TIR) would undergo lossless TIR rather than lossy reflection at the ESR film. The latter only interacted with the fraction of light that did not undergo TIR and escaped from the fused silica tube assembly. Finally, a 0.5-mm thick cadmium sheet was wrapped around the outside of the assembly to prevent any thermal neutrons from entering the detector through the lateral surface. This comprehensive fabrication process successfully created the oil-based scintillating detector, a picture of which is shown in Fig. 2.

Pulse Height Spectra Measurement Setup

Pulse height spectra were recorded by summing the outputs from the anodes of the two PMTs and conditioning it with an in-house charge-sensitive pre-amplifier with a sensitivity of 10 mV/pC and a decay constant of 200 µs, followed by an Ortec 672 shaping amplifier set to a gain of 34 and a shaping time constant of 1.5 µs. The processed signal pulses were then analyzed using an Ortec 927 ASPEC multi-channel analyzer (MCA). Data acquisition was performed with the Ortec Maestro software. The voltages of PMT A and PMT B were adjusted to 1000 V and 983 V, respectively, to match the respective positions of the ⁶Li neutron capture peak in the pulse-height histograms, thus ensuring equal gain of the two PMT/pre-amplifier chains.

Die-Away Time Measurement Setup

To measure the time distribution of neutron detection events, the PMT outputs from both sides were fed into a chain consisting of a preamplifier, Ortec 672 shaper/amplifier, followed by an Ortec 551 timing Single Channel Analyzer (SCA). The SCA threshold was adjusted only to include the neutron peak from the system’s pulse height distribution. The arrival times of the pulses were then recorded using a pulse train recorder [30].

Position Sensitivity Measurement Setup

To assess the position sensitivity capability of the composite detector, we conducted experiments using a tungsten-aperture collimated ²⁵²Cf fission neutron source. The source was placed inside a 2-inch (5.08 cm) diameter and 4 inch (10.16 cm) long cylindrical tungsten collimator that had a 3/8 inch (0.9525 cm) opening through the center. The source was pushed to one side of the collimator and placed in direct contact with the detector. The source and collimator were placed at various positions along the lateral surface of the composite detector to measure its position sensitivity. The signals from the PMTs on each side were connected to independent preamplifiers and shaper channels (not summed) with the same parameters as described in the Pulse Height Measurement Setup section. The high voltages of both PMTs remained as described in the Pulse Height Measurement Setup configuration to preserve the gain matching. The pulses from each side were recorded in coincidence using a digital oscilloscope. To calculate a residual metric, 105 electrical pulses were accumulated at each spatial position.

Cross Calibration and Dynamic Range Measurements Setup

A series of experiments were conducted at the Central Health Physics Calibration Facility (CHPCF at Los Alamos National Laboratory, TA-36) in the neutron-free-in-air (NFIA) facility, which is designed to reduce undesired neutron scattering by the experimental hall. A ²⁵²Cf fission neutron source with a NIST (National Institute of Standards and Technology) traceable calibration was used. To obtain the total neutron flux [neutrons/(s⋅cm²)] at the various ⁶Li glass composite scintillator detector distances from the ²⁵²Cf source, both direct and scattered contributions were calculated using MCNP6.2 coupled with the ENDF VIII library (see section Reference measurements and neutron flux at the Central Health Physics Calibration Facility)

Reference Measurements and Neutron Flux at the Central Health Physics Calibration Facility

The Central Health Physics Calibration Facility (CHPCF) is a secondary standards laboratory at LANL that maintains several neutron, gamma, and beta reference fields for the routine calibration of health physics instruments. These fields are also made available for special studies by researchers both within and outside the Laboratory.

As a secondary standards lab, all reference fields at the CHPCF are traceable to a national or international standard. Of interest to this study is the neutron-free-in-air (NFIA) reference field located inside a large hall (12 m × 12 m × 7 m). The neutron sources (²⁵²Cf) are stored below grade and raised pneumatically from a control console outside the hall. The sources are raised about 4 m above the concrete floor to a low-scatter aluminum platform upon which sits the irradiator assembly. A linear positioning stage (LPS) located on the platform allows precise positioning of instruments with respect to the irradiator. Positioning of the LPS at distances ranging from 30–300 cm from the source is also done remotely from the control console. Instruments placed on the LPS are aligned vertically and horizontally using lasers. The wall directly behind (and nearest) the irradiator is made of drywall to minimize scattering. The remaining walls in the hall are constructed of concrete.

The NFIA sources are calibrated at NIST (National Institute of Standards and Technology) with respect to 4π emission rate using a Mn-bath technique with typical uncertainties on the order of 1.6% (k = 1). MCNP6.2, coupled with ENDF VIII cross-section library, was used for calculations of neutron fluence based on a comprehensive model of the NFIA facility (including the source and irradiator construction) normalized per source neutron to provide the neutron (spectral) flux. These MCNP-calculated quantities were then scaled based on the NIST-calibrated emission rate. The irradiator software at the control console automatically applies a decay-corrected dose rate for each source and distance daily.

For the purposes of the performance evaluation of the ⁶Li-composite-based detector, the direct and scattered contributions to neutron fluence at selected distances from a neutron source were calculated using MCNP6.2 as described before, but this time considered the source term to be comprised of not only the doubly encapsulated source but also the aluminum rabbit and the immediate environs of the irradiator assembly. The rest of the hall was modeled as a void, and point detector tallies (F5 type) at distances of interest were used to determine the “direct” fluence based on the total F5 tally score. Subtracting these results from those based on the earlier comprehensive model gave the combined contribution of room and air scatter to the total fluence and dose. The least squares fit of the “direct” fluence showed a 1/d² dependence for distances ranging from 30 to 300 cm, as expected for a point-like source as it is a ²⁵²Cf fission neutron source. When calculated in this fashion, the “direct” fluence was 10% higher in the horizontal (i.e. the instrument) plane than based on just the bare source 4π emission rate. This increase was partly due to the inherent anisotropy of the cylindrical source capsule (~ 2% effect) supplemented by scatter within the extended source term. Measurements using neutron spectrometers such as Bonner spheres and ROSPEC have validated the dose rates derived using the approach described above. A detailed uncertainty budget analysis for the NFIA indicated typical total uncertainties on the fluence of ± 6.0% for the NFIA ²⁵²Cf sources (k = 2). These uncertainties are well within those expected by DOE for a secondary standards laboratory, i.e. ±20% (k = 2).

Competing interests

The authors declare no competing interests.

Acknowledgments

Research presented in this article was supported by the Laboratory Directed Research and Development program of Los Alamos National Laboratory under project number 20210431ER. We acknowledge the support received from N. Wehmann, C. Roybal, J. Bland, J. Mendez, M. Olivas, D. Jones, R. Weinmann-Smith, and M. Swinhoe.

Byrne, J. (2013). Neutrons, Nuclei and Matter: an exploration of the physics of slow neutrons, Dover Publications.
Segre, E. (1978). Nuclei and Particles: an introduction to Nuclear and Subnuclear physics. Benjnamin-Cummings Pub Co. ISBN-10 0805386017.
Beckhurts, K.H., Wirtz, K. (1964). Neutron Physics. Springer. ISBN-10:3642876161.
Knoll, G.F. (2010). Radiation Detection and Measurement. John Wiley & Sons Inc. ISNB-10:9780470131480.
Reilly, D., Ensslin, N., Smith, H., Jr, Kreiner, S. (Eds) (March, 1991). Passive Nondestructive assay of nuclear materials. US Nuclear Regulatory Commission report NUREG/CR-5550. Los Alamos National Laboratory report LA-UR-90-732. ISBN 0-16-032724-5.
Neutrons for the Nation. Discovery and Applications while Minimizing the Risk of Nuclear Proliferation. A Report by the American Physical Society on Public Affairs. (July 2018). APS Physics. https://www.aps.org/policy/reports/popa-reports/upload/APSNeutronsfortheNation.pdf
Cho, O. (2009) Helium-3 shortage could put freeze on low-temperature research”, Science 326, 778-779.
Sacchetti, F., Colonna, N., Faccini, R., Guerard, B., Hall-Wilton, R., Murtas, F., Petrillo, C., Pietropaolo, A., Rhodes, N., Quintieri, L., Tardocchi, M., & Valente, P. (2015). "³He-free neutron detectors and their applications." Eur. Phys. J. Plus, 130, 53.
Croft, S., Favalli, A., Hauck, D. K., Henzlova, D., Santi, P. A. (2012). Feynman variance-to-mean in the context of passive neutron coincidence counting. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 686, 136-144.
Favalli, A., Croft, S., Santi, P. (2015). Point model equations for neutron correlation counting: Extension of Böhnel׳s equations to any order. Nuclear Instruments and Methods in Physics, 795, 370-375
Henzlova, D., Kouzes, R., McElroy, R., Peerani, P., Aspinall, M., Baird, K., Bakel, A., Borella, M., Bourne, M., Bourva, L., Cave, F., Chandra, R., Chernikova, D., Croft, S., Dermody, G., Dougan, A., Ely, J., Fanchini, E., Finocchiaro, P., Gavron, V., Kureta, M., Ianakiev, K. D., Ishiyama, K., Lee, T., Martin, C., McKinny, K., Menlove, H. O., Orton, C., Pappalardo, A., Pedersen, B., Peranteau, D., Plenteda, R., Pozzi, S., Schear, M., Seya, M., Siciliano, E., Stave, S., Sun, L., Swinhoe, M. T., Tagziria, H., Vaccaro, S., Takamine, J., Weber, A. L., Yamaguchi, T., & Zhu, H. (2015). "Current Status of Helium-3 Alternative Technologies for Nuclear Safeguards." U.S. DOE, Los Alamos National Laboratory report (LA-UR-15-21201).
Kouzes, R. T., Ely, J. H., Erikson, L. E., Kernan, W. J., Lintereur, A. T., Siciliano, E. R., Stephens, D. L., Stromswold, D. C., Van Ginhoven, R. M., & Woodring, M. L. (2010). Neutron detection alternatives to 3He for national security applications. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 623(3), 1035-1045.
Peerani, P., Tomanin, A., Pozzi, S., Dolan, J., Miller, E., Flaska, M., Battaglieri, M., De Vita, R., Ficini, L., Ottonello, G., Ricco, G., Dermody, G., & Giles, C. (2012). "Testing on novel neutron detectors as an alternative to ³He for security applications." Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 696, 110-120.
Lintereur, A.T., Ely, J.H., Kouzes, R.T., Siciliano, E.R., Swinhoe, M.T., Woodring, M.L.(2012) Alternatives to Helium-3 for neutron multiplicity counters, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC), Anaheim, CA, USA, 2012, pp. 547-553, doi: 10.1109/NSSMIC.2012.6551168.
Ianakiev, K.D., Swinhoe, M.T., Favalli, A., Chung, K., MacArthur, D.W., (2011). ⁶Li foil scintillation sandwich thermal neutron detector, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Volume 652, Issue 1.
Ianakiev, K. D., Hehlen, M. P., Swinhoe, M. T., Favalli, A., Iliev, M. L., Lin, T. C., Bennett, B. L., Barker, M. T. (2015). "Neutron detector based on Particles of ⁶Li glass scintillator dispersed in an organic lightguide matrix." Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 784, 189-193.
Hehlen, M. P., Wiggins, B. W., Favalli, A., Iliev, M., Ianakiev, K. D. (2018). "Light propagation in a neutron detector based on ⁶Li glass scintillator particles in an organic matrix." Journal of Applied Physics, 124, 124502. https://doi.org/10.1063/1.5047433.
Wiggins, B. W., Favalli, A., Iliev, M., Ianakiev, K. D., Hehlen, M. P. (2019). "Computational investigation of arranged scintillating particle composites for fast neutron detection." Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 915, 17-23.
McGregor, D. S., Klann, R. T., Gersch, H. K., Yang, Y. H. (2001). Thin-film-coated bulk GaAs detectors for thermal and fast neutron measurements. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 466(1), 126-141.
https://scintacor.com/products/6-lithium-glass/ last accessed October 14^th, 2023.
Favalli, A., Iliev, M., Ianakiev, K., Hunt, A. W., Ludewigt, B. (2018). Delayed gamma-ray spectroscopy with lanthanum bromide detector for non-destructive assay of nuclear material. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 877, 192-196.
Kulesza, J. A., Adams, T. R., Armstrong, J. C., Bolding, S. R., Brown, F. B., Bull, J. S., Burke, T. P., Clark, A. R., Forster III, R. A., Giron, J. F., Grieve, T. S., Josey, C. J., Martz, R. L., McKinney, G. W., Pearson, E. J., Rising, M. E., Solomon Jr., C. J., Swaminarayan, S., Trahan, T. J., Wilson, S. C., Zukaitis, A. J. MCNP® Code Version 6.3.0 Theory & User Manual. Los Alamos National Laboratory Tech. Rep. LA-UR-22-30006, Rev. 1. Los Alamos, NM, USA. September 2022.
Hendricks, J.S., Swinhoe, M.T., Favalli, A. Monte Carlo N-Particle Simulations for Nuclear Detection and Safeguards. (2022) Springer.
Menlove, H. O., Rael, C. D., Kroncke, K. E., & DeAguero, K. J. (2004). Manual for the Epithermal Neutron Multiplicity Detector (ENMC) for Measurement of Impure MOX and Plutonium Samples. Los Alamos National Laboratory, USA
Favalli, A., Mehner, H.-C., Crochemore, J.-M., Pedersen, B. (2009). Pulsed Neutron Facility for Research in Illicit Trafficking and Nuclear Safeguards. IEEE Transactions on Nuclear Science, vol. 56, no. 3, pp. 1292-1296, June 2009, doi: 10.1109/TNS.2009.2015297.
Favalli, A., Guler, N., Henzlova, D. et al. (2019). Characterizing laser-plasma ion accelerators driving an intense neutron beam via nuclear signatures. Sci Rep 9, 2004.
Favalli, A., Bridgwater, J., Croft, S., Falk, K., Fernandez, J. C., Gautier, D., ... & Swinhoe, M. (2016, September). Laser-driven neutron source for detection of nuclear material. In Proceeding of 2016 Advances in Nuclear Nonproliferation Technology and Policy Conference, Santa Fe, NM, USA.
Huang, C.K., D. P. Broughton, S. Palaniyappan, A. Junghans, M. Iliev, S. H. Batha, R. E. Reinovsky, A. Favalli. (2022) High-yield and high-angular-fluence neutron generation from deuterons accelerated by laser-driven collisionless shock. Appl. Phys. Lett.; 120 (2): 024102. https://doi.org/10.1063/5.0075960
Roth, M, et al. (2013). Bright laser-driven neutron sources based on relativistic transparency of solids. Phys. Rev. Lett, 110, 4.
https://www.iki.kfki.hu/radsec/research/PTR-32_ThrowAway2.pdf (last access, Dec.27, 2021).

There is NO Competing Interest.

Download PDF

Version 1

posted

You are reading this latest preprint version

Next-Generation Neutron Detection Using a ⁶Li Glass Scintillator Composite

Status:

Version 1

Abstract

Figures

Introduction

Design and Detector Fabrication

Results

Detector Performance Evaluation

Pulse-Height Spectra

Neutron/Gamma Discrimination Stability to Variation of High Voltages

Cross Calibration and Dynamic Range

Neutron Detection Efficiency

Die-Away Time

Position Sensitivity of the Neutron Flux Interaction

Discussion

Applicability to Neutron Multiplicity Counting

Future Directions

Materials and Methods

Detector Fabrication

Pulse Height Spectra Measurement Setup

Die-Away Time Measurement Setup

Position Sensitivity Measurement Setup

Cross Calibration and Dynamic Range Measurements Setup

Declarations

Competing interests

Acknowledgments

References

Additional Declarations

Status:

Version 1