Patent Application
20240159921
The present invention relates to the field of radiation detection, more specifically the parallel detection of neutrons and gamma rays.
Neutrons emitted from spontaneous or induced fission provide a strong signature for the presence of Special Nuclear Material (SNM). The detection of SNM is of particular interest for homeland security applications — portal monitors, handheld instruments, backpack monitors, etc. Further, the simultaneous detection of neutrons and gamma rays as such is of particular interest for the monitoring of nuclear (power) systems and particle accelerators as well as for geophysical borehole investigations.
An efficient simultaneous detection of (slow) neutrons and gamma rays by just one detector instead of two types of detectors thereby would be a major advantage for the detection of SNM—not only in the context of homeland security applications.
Neutrons are neutral particles carrying a mass and a magnetic moment. Due to their nature, the interaction of neutrons with matter (disregarding their magnetic moment) is bound to an interaction with nuclei, and therefore occurs with low probabilities for most materials.
The (kinetic) energies of free neutrons span a range of many orders of magnitude. A large fraction of neutrons emitted in spontaneous or induced fission reaction have energies of up to 10 MeV, classified into three neutron energy regimes: slow neutrons (<1 keV), intermediate neutrons (1 keV to 0.5 MeV) and fast neutrons (0.5 MeV to 10 MeV).
The interaction of neutrons with matter—as neutron capture or neutron scattering—strongly depends on the neutron energy. Slow neutrons, in particular thermal neutrons, i.e. neutrons whose kinetic energy distribution corresponds to the Boltzmann distribution at room temperature, are more likely to be captured than fast neutrons. The neutron cross section defines the likelihood of an interaction between a neutron and a nucleus.
Most of the conventional neutron detection methods are based on the 3He(n,p) 3H, the 6Li(n,α) 3H, or the 10(n,α)7 Li reactions. These reactions emit ionizing particles, which are easy to detect and well distinguishable from gamma rays. Further, 3He, 6Li and 10B are distinguished by a large capture cross section for thermal neutrons.
The emitted charged particles generate short but dense ionization tracks with ranges of a few micrometers in solids, and up to some millimeters in gases of appropriate pressure. Consequently, the neutron-capturing component must be either a constituent of the detection medium (as, e.g., in 3He and BF3 tubes, scintillators comprising 6Li), or it must be provided as a thin layer in close contact with the detector medium (as, e.g., in 10B-coated proportional counters). The neutron-capturing component should consist of isotope-separated nuclides to avoid competing neutron-capture reactions reducing the detection efficiency. Therefore, corresponding detectors are often expensive. In addition, 3He has become very expensive recently due to a serious shortage.
Previously, scintillator crystals comprising natural Li or isotopically enriched or separated 6Li for a parallel detection of neutrons and gamma rays have been developed, amongst them CLYC (Glodo et al., IEEE Trans. Nucl. Sci. 55 (2008) 1206), CLLB (Shirwadkar et al., Nucl. Instr. Meth. In Phys. Res. A 652 (2011) 268), or NAIL (Yang et al., IEEE TNS 64 (2017) 2406). These advanced scintillator materials allow to distinguish neutron and gamma-ray interactions by analyzing the pulse shape of the detector signals. Neutron interactions produce tritons and alpha particles, while gamma-ray interactions produce photoelectrons, Compton electrons, or electron-positron pairs. Hence, scintillation detectors providing a reasonable efficiency for neutron detection as well as good spectroscopic performance for gamma ray detectors may be constructed using such scintillator crystals. However, these scintillator crystals are more complex and therefore expensive, rendering the price of corresponding detectors to numbers exceeding the price of common NaI(Tl) detectors of equal size by a factor of at least 2 to 5. As a result, such detectors may even become more expensive than a combination of a common NaI(Tl) gamma detector and a separate neutron detector (e.g., a 3He tube).
An alternative and cheaper solution for neutron detection in applications that are utilizing gamma ray detectors anyway is based on neutron capture reactions with subsequent emission of gamma rays, so-called Neutron Capture Detectors (NCD). An NCD detects neutrons by measuring the gamma radiation following neutron captures in a neutron-gamma converter arranged around or comprised in the gamma-ray detector, as disclosed in EP 24 600 32 A1 and WO 2011/012155 A1. Naturally, a major challenge of any NCD is to discriminate between the signal of gamma rays stemming from neutron capture and signal of ubiquitous gamma rays of external sources.
In neutron capture reactions a neutron is trapped by a nucleus AX resulting in a daughter nucleus A+X with an excitation corresponding to the neutron separation energy Sn(A+1X) of the daughter nucleus. Subsequently, this energy is released in form of a gamma-ray cascade. The cascade may involve continuum states and distinct gamma transitions. The discrete lines represent fingerprints of the capturing nuclide or a mixture of nuclides, e.g. of chemical elements in their natural isotope composition. The gamma rays emitted upon thermal neutron capture, in particular the so-called “prompt gamma rays from slow neutron capture” are used for elemental analysis by neutron activation, e.g. in homeland security applications and borehole investigations. The energies and intensities of prompt gamma rays from neutron capture are continually being characterized and tabulated, e.g. in the “Database of prompt gamma rays from slow neutron capture for elemental analysis”, IAEA (2007).
The detection of gamma rays generally is based on the measurement of electric currents. The electrons and positrons produced in gamma-ray interactions with the detector medium either directly produce electron-hole pairs in semiconductor materials or are converted into electric current by means of a scintillator material and a photodetector with internal or external amplification.
EP 17 170 025.5 relates to a method to detect both gamma rays and fast neutrons using a scintillator crystal providing different light pulse shapes for recoil electrons and recoil nuclei, such as Tl-doped NaI, for homeland security applications. This method is restricted, however, to the detection of fast neutrons.
EP 24 600 32 A1 relates to an NCD comprising a gamma ray detector supplemented with neutron-gamma converters for an as much as possible calorimetric detection of the gamma-ray cascades following neutron capture reactions. The limitation of this approach is due to the sensitivity of an NCD for ambient gamma radiation.
WO 2011/012155 A1 relates to an NCD consisting of one or more gamma ray detectors comprising a neutron-gamma converting isotope of a chemical element used as constituent or dopant of the active detector material for an as much as possible calorimetric detection of the gamma-ray cascades following neutron capture reactions. The limitation of this approach is due to the sensitivity of an NCD for ambient gamma radiation.
So far, the only way to distinguish between neutron capture gamma rays and gamma rays of other sources in an NCD is the rejection of detector signal corresponding to energy depositions below approximately 3 MeV, the NCD threshold, for being considered as neutron signals, since common radioactive nuclides do not emit noticeable fractions of gamma rays with energies above 2615 keV.
The major drawback of signal rejection below the NCD threshold is that neutron-capture gamma-ray cascades leading to energy dispositions in the gamma-ray detector below the NCD threshold are also rejected. Such cases would not occur if the NCD warranted a truly calorimetric detection of neutron-capture gamma-ray cascades. In practice the calorimetric regime could only be approached by an appropriate construction and by using a reasonably large and dense gamma-ray detector. Hence, the NCD threshold reduces the probability of detecting neutrons, i.e. the sensitivity of the detector. The sensitivity loss scales with the deviation from a truly calorimetric regime. Therefore, large sizes and corresponding high weights of NCD are required to counteract on the delimiting conditions of an NCD imposed by the NCD threshold.
Yakushev et al proposed in Nucl.Inst.Meth. A 848 (2017) 162 a method for thermal neutron detection in scintillators comprising iodine, as common Tl-doped NaI scintillators. Natural iodine only consists of the stable 127I isotope which provides a relatively large thermal neutron capture cross section of about 6.2 barn. This cross section is large enough to absorb a considerable fraction of thermal neutrons interacting with a Tl-doped NaI scintillator of usual dimensions (2-3″ diameter, 2-3″ height). The de-excitation cascade following neutron capture on 127I often involves the 137.8 keV state in 128I having a half-life (T1/2) of 845 ns. The de-excitation cascade of the 137.8 keV state in 128I involves low-energy gamma rays and converted transitions of which a large fraction is detected in the NaI detector itself—with a certain delay following the prompt part of the de-excitation cascade which may also be at least partially absorbed and then generate a signal in the NaI detector. Consequently, the resulting electric current signal of a gamma-ray cascade upon thermal neutron capture sometimes exhibits a double pulse within a time window of a few microseconds, wherein the delayed pulse fits to an energy deposition of approximately 138 keV. Hence, neutrons are detected if the distribution of time differences of double pulses exhibits the 845 ns decay component, whereby the strength of the 845 ns decay component reflects the thermal neutron flux the NCD is exposed to.
Using the technique of Yakushev et al, thermal neutrons can be discriminated from a gamma radiation background stemming from other (external) sources as long as the background is low enough to not significantly contaminate the double-pulse time-difference distribution by random pulse pile-ups. According to Yakushev et al, however, the analysis of the double pulses resulting from delayed gamma rays upon neutron capture is only possible because the time constant of the delayed gamma emission, 845 ns, is large in comparison to the light decay time of the Tl-doped NaI scintillator used in the demonstration which is approximately 250 ns.
While Yakushev et al's method is based on common NaI(Tl) scintillators, which are relatively cheap and widely used for gamma-ray spectroscopy in numerous applications, a number of disadvantages results from said method:
Said method requires the delayed pulses to not pile-up with the prompt signal in order to measure their energy deposited in the scintillator crystal and the time delay of distinct signals following each other.
Sakharov et al. (Nucl. Phys. A528 (1991) 317) investigated the level scheme of 128I. They state that the long-lived 137.85 keV state in 128I is often fed by another long-lived state, the 167.3 keV state having a half-life of 175 ns. The corresponding 29.5 keV transition is almost completely converted. According to their research, the 167.3 keV state is mainly fed by a 2-step transition, comprising a 142 keV gamma transition from a 376 keV state to a 234 keV state, slightly converted, followed by a strongly converted 67 keV transition to the 167.3 keV state. They found that this cascade of gamma rays and conversion electrons is often involved in the decay chain following neutron capture in 127I.
Further, pulse pile-up is a well-known and persisting problem in the field of radiation detection. Nuclear reactions and accompanying de-excitation cascades occur randomly. On the one hand, the time difference between consecutive pulses recorded in gamma ray detectors follows the probability density function of an exponential distribution. Hence, small time differences between consecutive pulses are more likely than large time differences. On the other hand, the pulse duration is controlled by the decay time of the scintillator which defines the time scale for resolving distinct signals. In case the time difference between consecutive pulses becomes shorter than the pulse duration, the signals pile up. The energy determination of a corresponding event, meaning an integration of the signal over time to measure the charge content, fails in case of a pulse pile-up, i.e. it cannot provide correct results for the two piled-up signals.
Pulse pile-ups can be detected, e.g. by comparing the ratio of two shifted integrals over each signal pulse. While single pulses will always return the same ratio of the shifted integrals, the ratio of shifted integrals of pulse pile-ups deviates.
Based on the identification of pile-up events, the signals of piled-up pulses often are rejected. The major disadvantage of pulse pile-up rejection, however, is the loss of potentially valuable signal.
In a publication by Scoullar et al, AIP Conf. Proc. 1412 (2011) 270, a method for real-time decoding of pulse pile-up events for a range of detectors, including scintillation based radiation detectors, is suggested. This method is model-based and characterizes the number, time-of-arrival and energy of all events in the detector output. In doing so, composite events are decoded and energy and time-of arrival of multi-pulse pile-ups are reconstructed, with a pulse pair resolution of down to 50 ns. An important constraint of this approach is the fixation of the pulse shape expected from the detector.
Pulse pile-up decoding or reconstruction techniques have so far been of particular advantage in high-flux pulsed beam experiments with low duty cycle and relatively slow detectors, such as Free Electron Lasers for which pulse pile-ups are unavoidable and signal collection is expensive. In usual spectroscopic measurements higher count rates lead to a more severe disturbance of the spectrum due to pulse pile-up. Sufficiently high input rates, however, usually allow the partial loss of signals, rendering the pulse pile-up rejection the method of choice.
Hence, an object of the present disclosure is to provide a system and a corresponding method allowing the detection of slow neutrons with common iodine-comprising scintillation detectors in a way that minimizes the above mentioned drawbacks at least partially thereby improving the sensitivity of common iodine-comprising scintillation detectors for neutron detection and the feasibility of their large-scale deployment for that purpose.
This problem may be solved by a radiation detection system and a method for detecting neutrons using delayed triple-coincidence capture gamma-ray detection. Preferred embodiments are described herein.
More specifically, the disclosure provides a radiation detection system for a parallel detection of gamma rays and neutrons, comprising a gamma-ray detector, a digitizer and an analyzer. The gamma-ray detector comprises a scintillator crystal comprising 127I and a photodetector with an amplifier, wherein said scintillator crystal is adapted to convert energy deposited by gamma rays or conversion electrons to optical photons. The photodetector is adapted to convert the optical photons to an electrical signal. The amplifier is adapted to amplify the electrical signal, wherein the electrical signal is in a known relationship with the energy deposited by the detected gamma rays or conversion electrons in said scintillator crystal. The digitizer comprises sampling Analog to Digital Converters, ADC, wherein said digitizer is adapted to sample the electrical signal of said gamma-ray detector with a predetermined frequency of at least 20 mega samples per second to generate digitized time series of the electrical signal. The analyzer is coupled operatively to said digitizer. Said digitizer is adapted to transmit the digitized time series to the analyzer, and the analyzer is adapted to analyze the digitized time series in order to identify signal components in the digitized time series with consecutive time delays between each other of at least 20 nanoseconds and at most 10 microseconds, wherein each signal component is due to an energy deposited in the scintillator crystal.
The analyzer is further adapted to identify a primary signal component in the digitized time series corresponding to an energy deposition E0 in a predefined range. The analyzer is adapted to search for a first delayed signal component in the digitized time series, the first delayed signal component corresponding to an energy deposition E1 of about 30 keV, preferably between 15 keV and 50 keV, and following the primary signal component in time, and to search for a second delayed signal component in the digitized time series, the second delayed signal component corresponding to an energy deposition E2 of about 138 keV, preferably between 100 keV and 200 keV, and following the first delayed signal component in time. Moreover, the analyzer is adapted to count the number of digitized time series comprising at least the first delayed signal component and the second delayed signal component as neutron events, thereby providing a measure for a neutron flux the scintillator crystal is exposed to.
The system has the advantage that no trigger holdoff is required that prevents detecting delayed signal components with a time difference of less than 1.8 μs to the primary signal component. Thus, also the signal component resulting from de-excitation of the long-lived excited state at about 167 keV having a half-life of just about 175 ns can be exploited to sharpen the search criteria for gamma rays upon neutron captures in the scintillator crystal. As a result, the radiation detection system manifests itself with a much stronger background suppression compared to other methods of neutron detection via neutron-capture gamma rays which renders the radiation detection system applicable in much stronger gamma radiation fields.
Preferably, the digitizer is even adapted to sample the electrical signal of said gamma-ray detector with a predetermined frequency of at least 100 mega samples per second.
Preferably, the analyzer is further adapted to quantify at least one event parameter of a group of event parameters for each of said digitized time series, the group of event parameters comprising:
Moreover, the analyzer is adapted to evaluate for each of said digitized time series whether at least one event parameter of the group of event parameters fulfills a predefined criteria to classify said digitized time series as neutron event, and to count the number of digitized time series classified as neutron event, thereby providing a measure for the neutron flux the scintillator crystal is exposed to.
A preferred embodiment of the system has the further advantage that even single neutron counts, in particular thermal neutron counts may be registered even in environments with an intense ambient gamma radiation. The relatively precise quantification of at least one event parameter is sufficient to suppress the background due to random triple pulse events such that individual digitized time series representing delayed triple coincidences can be attributed to neutron event, i.e. neutron captures in 127I of the scintillator crystal.
Preferably, the photodetector of the gamma-ray detector is a photomultiplier tube, a silicon photomultiplier (SiPM), or an avalanche photodiode.
Further preferably, the photomultiplier tube comprises a super-bialkali photokathode or an ultra-bialkali photokathode.
Preferably, said scintillator crystal is a NaI crystal, preferably with Tl doping, or a CsI crystal, preferably with Na or Tl doping, or a NaI crystal, preferably with Tl doping, comprising Li or B as co-dopant.
The radiation detection system therefore has the advantage that it can be equipped with common iodine-comprising crystals which are relatively cheap and already widely applied in gamma-ray spectroscopy.
Preferably, the gamma-ray detector, the digitizer and the analyzer of the radiation detection system are small enough to be arranged in a handheld device.
Preferably, the gamma-ray detector, the digitizer and the analyzer of the radiation detection system are small enough to be arranged in a backpack.
The radiation detection system therefore is particularly well adapted to an application in the field of homeland security.
Further, embodiments of the invention provide a method to detect neutrons and gamma rays, utilizing the radiation detection system. The method comprises a number of steps according to which neutrons interact with 127I of the scintillator crystal generating 128I upon neutron capture, wherein 128I de-excites both under a prompt emission of gamma radiation and at least partially via at least one long-lived excited state feeding another long-lived excited state which at least sometimes leads to two delayed energy depositions in the scintillator crystal corresponding to the two delayed de-excitation steps producing either conversion electrons or gamma radiation interacting with the scintillator crystal, following a primary signal component which is due to the prompt emission. The gamma-ray detector generates an electrical signal in succession of an interaction between gamma rays or conversion electrons and said scintillator crystal, wherein said electrical signal is in a known relationship with the energy deposited by the detected gamma rays or conversion electrons in said scintillator crystal. The digitizer samples the electrical signal of said gamma-ray detector with a predetermined frequency of at least 20 mega samples per second, preferably at least 100 mega samples per second, to generate a digitized time series of the detected gamma rays, and transmits the digitized time series to the analyzer. The analyzer identifies a primary signal component in the digitized time series corresponding to an energy deposition E0 in a predefined range, searches for a first delayed signal component in a digitized time series, the first delayed signal component corresponding to an energy deposition of about 30 keV, preferably between 15 keV and 50 keV, and following the primary signal in time, searches for a second delayed signal component, the second delayed signal component corresponding to an energy deposition of about 138 keV, preferably between 100 keV and 200 keV, and following the first delayed signal component in time. Moreover, the analyzer counts the number of digitized time series comprising at least the first delayed signal component and the second delayed signal component as neutron events, thereby providing a measure for the neutron flux the scintillator crystal is exposed to.
Preferably, the method further comprises the steps according to which the analyzer quantifies at least one event parameter of a group of event parameters for each of said digitized time series, the group of event parameters comprising:
In a further step, the analyzer evaluates for each of said digitized time series whether at least one event parameter of the group of event parameters fulfills a predefined criteria to classify said digitized time series as neutron event, and eventually counts the number of digitized time series classified as neutron event, thereby providing a measure for the neutron flux the scintillator crystal is exposed to.
The method provides the advantage that it facilitates simultaneous gamma-ray spectroscopy and neutron detection with common, off-the-shelf scintillation detectors using iodine-comprising scintillators, e.g., NaI (Tl), CsI (Na) or CsI (Tl). Therefore, the method enables neutron detection capabilities in common digital gamma-ray spectrometers with iodine-comprising scintillation detectors without requiring additional hardware.
By searching for a first delayed signal component and a second delayed signal component in the digitized time series immediately after the primary signal, i.e., without a trigger holdoff as long as 1.8 μs of the digitizer, the neutron capture events can be detected much more efficiently due to a much stronger background suppression compared to other methods of neutron detection using iodine-comprising scintillator crystals, e.g., Yakushev's method. It were the inventors who first realized that the de-excitation of the long-lived excited state at 167 keV in 128I to the long-lived excited state at 138 keV in 128I provides an effective handle to discriminate the gamma ray de-excitation cascade following a neutron-capture event from random pulse pile-ups in a strong gamma-ray radiation environment. In particular, depending on the radiation background it may often not be sufficient to rely on the double-pulse structure/delayed double coincidences as seen by Yakushev et al, according to which the delayed pulses correspond to an energy deposition of about 138 keV, since random pulses with energies which are equal within the resolution of the scintillator crystal may occur. In contrast, random pulses fulfilling the search criteria of the method are much less likely, rendering said method to be much more robust in strong gamma ray radiation backgrounds which are typical for homeland security applications.
Preferably, the analyzer finds and identifies at least one delayed signal component in said digitized time series by means of pulse pile-up reconstruction techniques, utilizing a method to decompose the digitized time series comprising piled-up signal components into these components, thereby quantifying at least one of the group of event parameters.
The use of pulse pile-up reconstruction techniques has the advantage over the sequential pulse by pulse detection proposed by Yakushev that also signal components at a time delay of less than 1.8 μs from the primary signal component may be distinguished and detected. This increases the neutron detection efficiency with respect to prior art by a factor of 2 to 4 and also improves the effect-to-background ratio even when only searching for double pulses, i.e. using delayed double-coincidence events. The gain basically depends on the minimum time distance between primary and delayed signal component that is required by the pile-up reconstruction algorithm for an unambiguous signal decoding.
Most preferably, the analyzer further identifies another primary signal component in the digitized time series corresponding to an energy deposition in a predefined range, searches for a delayed signal component in the digitized time series, the delayed signal component corresponding to an energy deposition of around 138 keV, preferably between 100 keV and 200 keV, and following said primary signal component in time, quantifies at least one event parameter of the group of event parameters for each of said digitized time series, the group of event parameters further comprising:
evaluates for each of said digitized time series whether at least one event parameter of the group of event parameters fulfills a predefined criteria to classify said digitized time series as thermal neutron event, and counts the number of digitized time series classified as neutron event, thereby providing a measure for the neutron flux the scintillator crystal is exposed to.
The combination of searching for delayed triple coincidences and delayed double coincidences using pulse pile-up reconstruction techniques eventually yields even better effect-to-background ratios allowing the parallel detection of neutrons and gamma rays even in environments with an intense ambient gamma radiation.
In an embodiment of the invention, the method according to the invention is applied in a radiation detection system which is small enough to be arranged in a backpack.
In a further embodiment of the invention, the method according to the invention is applied in a radiation detection system which is small enough to be arranged in a handheld device.
Some embodiments of the invention are explained in the following with reference to the figures. The described examples are provided for illustrative purposes and are not intended to limit the scope of the invention.
Shown is a radiation detection system 1 comprising a gamma-ray detector 10, a digitizer 40 and an analyzer 50. The gamma-ray detector 10 comprises a scintillator crystal 20 comprising 127 J and a photodetector 30 with an amplifier, wherein said scintillator crystal 20 is adapted to convert energy deposited by gamma rays or conversion electrons to optical photons. The photodetector 30 is adapted to convert the optical photons to an electrical signal and the amplifier is adapted to amplify the electrical signal, wherein the electrical signal is in a known relationship with the energy deposited by the detected gamma rays or conversion electrons in said scintillator crystal 20.
Preferably, the photodetector 30 is a photomultiplier tube with intrinsic signal amplification. The photomultiplier may as well be an avalanche photodiode or a silicon photomultiplier (SiPM).
Most preferably, the photomultiplier tube comprises a super-bialkali photokathode or an ultra-bialkali photokathode.
The scintillator crystal 20 comprises 127I intrinsically, acting as neutron-gamma converter. A neutron-gamma converter is adapted to emit gamma radiation upon capturing neutrons.
More specifically, said scintillator crystal 20 is a NaI crystal, preferably with Tl doping, or a CsI crystal, preferably with Na or Tl doping, or a NaI crystal, preferably with Tl doping, comprising Li or B as co-dopant.
Needless to say, scintillator crystal 20 and the photodetector 30 may be wrapped in light reflecting materials and arranged in a housing to protect them from external light and humidity (not shown).
The embodiment of the radiation detection system of
The embodiment of the radiation detection system of
Preferably, the analyzer 50 is adapted to analyze the digitized time series and to determine a measure for the neutron flux, particularly the slow neutron flux, the detector 10 is exposed to in real-time.
The analyzer 50 is further adapted to identify a primary signal component in the digitized time series corresponding to an energy deposition E0 in a predefined range. The analyzer 50 is adapted to search for a first delayed signal component in the digitized time series, the first delayed signal component corresponding to an energy deposition E1 of about 30 keV and following the primary signal component in time, and to search for a second delayed signal component in the digitized time series, the second delayed signal component corresponding to an energy deposition E2 of about 138 keV and following the first delayed signal component in time. The analyzer 50 is further adapted to count the number of digitized time series comprising at least the first delayed signal component and the second delayed signal component as neutron events, thereby providing a measure for the neutron flux the scintillator crystal is exposed to.
In practice, it may be sufficient to determine the energy depositions of the delayed signal components coarsely within time and energy windows, the energy windows e.g. ranging from 15 keV to 50 keV for the first delayed signal component and from 100 keV to 200 keV for the second delayed signal component. For single neutron counts, in particular in environments with very intense ambient gamma radiation (count rates far beyond 1000 counts per second), it has been proven useful to further suppress the background of random triple pulse events by quantifying at least one event parameter of a group of event parameters more precisely.
Preferably, the group of event parameters comprises the time difference between the primary signal component and the first delayed signal component, the time difference between the primary signal component and the second delayed signal component, the time difference between the first delayed signal component and the second delayed signal component, and the corresponding energy depositions of each signal component in the scintillator crystal 20, and combinations thereof.
The inventors were the first who realized that digitized time series representing triple coincidence events, wherein each pulse of the triple coincidence events fulfills relatively coarse grained time and energy windows, may be attributed with high accuracy to neutron events if at least one of the event parameters fulfills a predefined fine-grained criteria. Examples for such a predefined criteria are a time difference between the first delayed signal component and the second delayed signal component of less than 4 μs, an energy deposition corresponding to the primary signal component of less than 300 keV, or the energy depositions corresponding to the first delayed signal component and second delayed signal component of 30 keV and 138 keV, respectively, within the resolution of the scintillator crystal 20.
More specifically, the analyzer 50 is preferably adapted to carry out pulse pile-up reconstruction techniques, wherein the pulse pile-up reconstruction techniques comprise a method to decompose the digitized time series into their signal components. Each signal component may have a pulse shape comprising an onset, the onset being characterized by a sharp rise in finite time, and an exponential decay, the exponential decay being essentially characterized by the light decay time of the scintillator crystal 20. Based on the reconstructed signal components, the analyzer 50 may quantify at least one of the group of event parameters. Time differences may be computed between the onsets of two consecutive signal components and energy depositions in the scintillator crystal 20 of a corresponding signal component may be determined by integrating the corresponding signal component.
The analyzer 50 may preferably carry out the pulse pile-up reconstruction techniques and subsequent computation of event parameters in real-time to enable an online detection of neutron events, e.g. for homeland security applications with moving sources.
Preferably, the gamma-ray detector 10 and corresponding electronics 40 and 50 are small enough to be arranged in a handheld device. Equally preferably, the gamma-ray detector 10 and corresponding electronics 40 and 50 are small enough to be arranged in a backpack. These form factors of the radiation detection system according to embodiments of the invention have been proven particularly useful for homeland security applications.
The digitized time series shown in
Digitized time series representing delayed triple coincidences comprise three signal components, namely a primary signal component recorded at a time t1, a first delayed signal component recorded at a time t2, and a second delayed signal component recorded at a time t3. In
Digitized time series representing delayed double coincidences comprise only two signal components, namely a primary signal component recorded at a time t1 and a delayed signal component recorded at a time t2. Delayed double coincidences as mentioned by Yakushev et al. are characterized by a delayed signal component corresponding to an energy deposition of around 138 keV in the scintillator crystal 20.
The three panels of
Based on the sampling rate of the digital spectrometer 40 used for this demonstration, 1.8 μs correspond to 450 samples in the example of
The examples of delayed triple coincidences shown in
In a first step of the method, neutrons interact with 127I of the scintillator crystal generating 128I up on neutron capture. 128I de-excites both under a prompt emission of gamma radiation and at least partially via at least one long-lived excited state feeding another long-lived excited state, which at least sometimes leads to two delayed energy depositions in the scintillator crystal 20 corresponding to the two delayed de-excitation steps producing either conversion electrons or gamma radiation interacting with the scintillator crystal 20, following a primary signal component which is due to the prompt emission.
In a next step of the method, the gamma-ray detector 10 generates an electrical signal in succession of an interaction between gamma rays or conversion electrons and said scintillator crystal 20. The electrical signal is in a known relationship with the energy deposited by the detected gamma rays or conversion electrons in said scintillator crystal 20.
Next, the digitizer samples the electrical signal of said gamma-ray detector with a predetermined frequency of at least 20 mega samples per second, preferably at least 100 mega samples per second, to generate digitized time series of the detected gamma rays, and transmits the digitized time series to the analyzer 50.
In a further step of the method the analyzer 50 identifies a primary signal component in a digitized time series, denoted as Pulse 1 in
It may be sufficient to simply count the number of digitizes time series in which a primary signal component, a first delayed signal component and a second delayed signal component are identified to determine a measure the neutron flux the scintillator crystal is exposed to. Preferably, however, to further suppress random background events, the analyzer 50 quantifies at least one of a group of event parameters for each of the digitized time series, and evaluates for each of said digitized time series whether at least one event parameter of the group of event parameters fulfills a predefined criteria to classify said digitized time series as neutron event. With counting the number of digitized time series classified as neutron event the analyzer 50 provides a measure for a thermal neutron flux the scintillator crystal 20 is exposed to.
As demonstrated by the examples shown in
The event parameters quantified by the analyzer 50 of the delayed triple coincidences shown in
The analyzer 50 is preferably adapted to compute a measure for the neutron flux, particularly the slow neutron flux, the detector 10 is exposed to, from the event parameters, as explained with respect to
The data was measured with the detector 10, as described with respect to
Additional detector load was generated by complementary Th ore and 137Cs sources arranged at various distances from the scintillator crystal 20. The distributions of time differences measured at different detector loads are denoted by different symbols.
The t32 distributions were fitted with the fit function f(t)=a·e−λ·t+b, with λ=0.693/0.845 μs representing the decay constant corresponding to a half-life of 845 ns, and then normalized with respect to the fit parameter a. The fits curves are also shown in the figure.
The fit parameter a measures the strength of the second delayed signal component (half-life of 845 ns) in the t32 distribution characterizing the second delayed transition, and thus the number of neutrons detected in a corresponding measurement. This means that all distributions shown in
divided by the background contribution in the same time interval
N
bck=∫t1t2b dt=b·(t2t1)
for an integration time ranging from t1=0.16 μs, which is the minimum detectable time difference between piled-up signal components with the pulse pile-up reconstruction algorithm used for analyzing the data, to t2=4.0 μs, which is a reasonable limit covering most neutron signals but not too much background. Note that the choice of the integration time range affects the E/B ratio. While t1 should always be set to the minimum detectable pulse delay to maximize the number of neutron counts considered, the optimum t2 value is affected by the background contribution, which varies with the actual detector load. The E/B ratios given here thus represent no general result, but exemplary values illustrating the trend: Time distributions as well as the corresponding E/B ratios reflect the rising background contribution towards higher detector load. However, the neutron contribution, distinguished by an exponential decay with a decay constant corresponding to 845 ns half-life, is still visible and can be quantified even at a detector load around 15 kcps in terms of the fit parameter a.
Preferably, the analyzer 50 computes a measure for the thermal neutron flux the scintillator crystal 20 is exposed to by determining a distribution of time differences t32 between the second delayed signal component and the first delayed signal component, fitting said distribution using a function comprising an exponential decay component corresponding to the half-life of the second delayed signal component and a scaling factor, wherein the scaling factor represents the number of (thermal) neutrons detected. Similarly, the number of neutrons may be obtained from the distribution of time differences between the first delayed signal component and the primary signal component, and even from the distribution between the second delayed signal component and the primary signal component, as long as the functions are adapted to comprise the corresponding exponential decay components.
Also preferably, the analyzer 50 computes another event parameter, the time difference t21 between the first delayed signal component and the primary signal. By selecting only such delayed triple-coincidence events as potential neutron counts being considered for further analysis, e.g. for accumulating the t32 distribution, where t21 does not exceed a predefined limit, preferably 2 μs, more preferably 1.5 μs, the effect-to-background ratio is increased at the expense of the effect count rate.
Also preferably, the analyzer 50 computes another event parameter, the energy deposition in the scintillator crystal 20 corresponding to the primary signal, E0. By selecting only such delayed triple-coincidence events as potential neutron counts being considered for further analysis, e.g. for accumulating the t32 distribution, where E_0 corresponds, within the energy resolution of the gamma detector, to the energy of a strong transition leading to the 167.3 keV state in 128I, preferably 142 keV or 67 keV, or to the sum energy of such transitions, the effect-to-background ratio is increased at the expense of the effect count rate.
Similarly preferably, the analyzer 50 computes other event parameters, the energy depositions in the scintillator crystal 20 corresponding to the first and/or the second delayed signal components, E1 and/or E2, respectively. By selecting only such delayed triple-coincidence events as potential neutron counts being considered for further analysis, e.g. for accumulating the t32 distribution, where E1 corresponds, within the energy resolution of the gamma detector, to the 30 keV transition feeding the 137.8 keV state in 128I, and/or E2 corresponds, within the energy resolution of the gamma detector, to the 138 keV energy released during de-excitation of the 137.8 keV state in 128I, the effect-to-background ratio is increased at the expense of the effect count rate.
In general, said predefined criteria applied to one or more event parameters, eventually combined to a set of criteria, can be used to further increase the effect-to-background ration on the expense of the effect (i.e., neutron) count rate. This allows dynamically adapting the sensitivity of the radiation detection system 1 with respect to neutron detection to the ambient gamma-ray background by using sets of more or less restrictive criteria, depending on said ambient background. In case of low ambient gamma background, weak or even no criteria applied to the event parameters maximize the neutron count rate at still acceptable background in the neutron counts. In case of a strong ambient gamma background, leading to an increased number of random delayed triple coincidences, more restrictive event selection criteria reduce the neutron sensitivity but may still allow quantifying the neutron flux the scintillator crystal 20 is exposed to.
The event parameter t21, representing the time difference (delay) between the delayed signal component corresponding to an energy deposition of around 138 keV in the scintillator crystal 20 and the primary signal component, was quantified by the analyzer 50. The time distributions shown in
The t21 distributions were fitted with the function f(t)=a·e−λ·t+b, with λ=0.693/0.845 μs, but (for better precision) normalized with the neutron contributions obtained from the delayed triple coincidence data shown in
It is obvious that the delayed triple coincidences provide a much better—in these exemplary cases by about an order of magnitude better—effect-to-background ratio for neutron counting than the delayed double coincidences, while the neutron contribution measured with delayed double coincidences is only about twice the neutron contribution measured with delayed triple coincidences. Nevertheless the neutron contribution measured with delayed double coincidences provides an additional and valuable measure for the neutron flux the scintillator crystal 20 is exposed to.
In a preferred embodiment of the invention the analyzer 50 measures the neutron flux the scintillator crystal 20 is exposed to from the neutron contribution derived from delayed double coincidences and from the neutron contribution derived from delayed triple coincidences, both contributions obtained from a pile-up re-construction algorithm that allows decomposing piled-up pulses and a subsequent event parameter analysis, thus increasing the accuracy of the neutron flux measurement.
In contrast to
The analysis of delayed triple coincidence events obviously provides a much better, by more than an order of magnitude better, effect-to-background ratio than prior art, while the effect rates are rather comparable. Moreover, if delayed double-coincidence are not excluded but also considered separately as shown in
Hence,
1 Radiation detection system
10 Gamma-ray detector
20 Scintillator crystal
30 Photodetector with an amplifier
E/B Effect-to-background ratio
ND Number of detectible neutron counts based on delayed double coincidences
NT Number of detectible neutron counts based on delayed triple coincidences
NDY Number of detectible neutron counts using prior art
This application is a continuation of International Application No. PCT/EP2021/061467, filed Apr. 30, 2021, the contents of which are incorporated by reference herein.
| Number | Date | Country | |
|---|---|---|---|
| Parent | PCT/EP2021/061467 | Apr 2021 | US |
| Child | 18541836 | US |
