Over the past decades, digital cameras have replaced their analog counterparts and despite the advancement in smaller pixel-pitch, faster read-out and lower read-out noise, the principles of how digital cameras record images have not changed much. For that matter, there are two main types of digital image sensors, namely the charge-coupled device (CCD) [1] and the complementary metal–oxide–semiconductor (CMOS) sensor [2], with the main difference between the two being the read-out. For this work, the sequential (or individual) pixel read-out of most CMOS sensors is crucial, because it enables the recently developed Timepix3 chip [3] to register the Time-of-Arrival (ToA) and the intensity of light in each pixel (through the Time-over-Threshold (ToT) method) for each event simultaneously and sparsely with up to 80 Mhits/s. Here, a “hit” or “event” is defined as any one of the 256×256 pixels composing the chip (with a pitch of 55 μm and a nominal time-resolution of 1.5625 ns) being activated, i.e. having some charge over a certain threshold being deposited on it at a specific ToA for a specific ToT. While the Timepix3 chip has a low sensitivity to visible light, Amsterdam Scientific Instruments (ASI) has developed the TPX3Cam [4–6], a light sensitivity-enhanced silicon sensor based Timepix3 camera system that can image and time-stamp light flashes of more than 103 photons with high quantum efficiency (> 90%) in the 400 to 1000 nm wavelength-range [7]. It is suited for application to imaging of electrons, ions or even single photons. For the latter case, the threshold of 103 photons is certainly too high for the TPX3Cam to detect a single photon directly. To overcome this limit, an image intensifier is needed to amplify light that would otherwise not reach the threshold. This applies in particular for example to light coming off a neutron scintillator screen, as it is the case for this work. With high enough amplification (>> 103), even single photons can be detected [5,6]. Utilizing this capability, we have developed a new approach for particle counting using scintillators, changing the principles of how the radiographs are recorded. This is achieved with parameters such as spatial and time resolution being defined after the measurement and not directly defined ab-initio by the pixel-size or the timing resolution of the system, but more so by the capability of detecting clusters of photons emitted by the scintillator after a particle interaction via the coincidence of many pixels activated over time.
Using simple Center of Mass (CoM) algorithms for single particle interactions and detection of the consequently emitted photons, the position of the particle’s interaction can be determined more accurately. This bears analogy to a similar digital imaging principle to that of the already applied principle in the determination of the position of stars in astronomy since 1980 [8]. In a sense, each interaction on the scintillator material creates a “star-like” signature. Thus, the idea to use CoM for improving position accuracy in digital imaging is basically as old as digital imaging itself. Nonetheless, applying a CoM to a single object, versus running such an analysis on single photon events for imaging purposes efficiently requires a detector technology that has shifted from a frame-based to an event-based system, such as the Timepix3 chip, as well as the necessary computational power to perform the analysis in real-time.
First studies on this concept applied to neutron imaging were reported on a frame-based imaging system [9]. However, the reported measurements faced limitations in acquisition dead-time and data processing time due to the slow read-out and noisy data inherent to the camera system in use for that study. Moreover, results were produced from post processed data only. While real-time single-event data-acquisition for neutron imaging was demonstrated on a Timepix quad detector combined with a neutron sensitive Microchannel-Plate (MCP) [10] with successful application to energy-resolved neutron imaging [11], the detector system in these applications has a fixed FoV based on the active area of 28×28mm2. For the work presented here, it was crucial to find a solution that would allow for a variable FoV as is essential state-of-the-art in neutron imaging today. Therefore, using a scintillator-based light collection system has the advantage compared to an MCP detector [10] to use optics for adjustment of the FoV, which is crucial for imaging of different length scales and objects of different sizes. A requirement that applies to essentially all neutron imaging beam-lines.
Aside from the FoV, additional focus of this study was to achieve as high of a spatial resolution as possible by taking advantage of the event-based data-acquisition. A recent example of using a TPX3Cam optical camera [5] showed that it is possible to improve the optical resolution of the camera beyond the pixel size of the sensor via CoM photon detection. An effective resolution equal to ⅕ of the pixel size of the sensor with a temporal resolution of better than 5 ns was demonstrated. In comparison to the described setup in that work, here the target for the camera concept is neutron imaging and, thus, essentially detecting the light a scintillator screen emits after a neutron is absorbed. Hence, the setup is very similar in principle with the caveat of light emission from the scintillator being time-dependent and generally based on several exponential decay functions [12–14]. Therefore, in this work, we go beyond the identification of single photons for optical resolution improvement to detecting clusters of photons that can be attributed to single neutron events, as will be shown. The added complexity in turn provides a solution to the afterglow problem for scintillators [13] and potentially enables discrimination of events in routine imaging measurements, where e.g. neutron and gamma events could be distinguished based on the light-yield decay properties of the interaction in the scintillator. This principle is generally also applied to photo-multiplier tube (PMT) detectors to reduce the gamma background for neutron transmission measurements [15–17].
The first measurements with neutrons for this work were conducted at the Beamline for neutron Optics and other Approaches (BOA) [18] at the Paul Scherrer Institute (PSI) in Switzerland, with data collected in “regular” white beam neutron imaging mode and Time-of-Flight ToF transmission mode enabled by a disk chopper installed in the beam to determine the scintillator response and to characterize the system. Furthermore, Bragg-Edge imaging [19,20] and dynamic response measurements of different neutron scintillator materials were performed, but will be reported at a later stage elsewhere. This first report focuses on the concept of the event-based image reconstruction. In particular, we demonstrate a two-step process in the analysis, that is applying a CoM on individual pixels for photon identification in a first step and a CoM algorithm on the individual photon events in a second step to calculate a more accurate position of neutrons interacting with the scintillator, making this detector concept a very powerful tool for radiography using, but not limited to, neutrons.
1.1. Event-mode imaging setup for neutrons
The components of the system, listed sequentially viewed from the source (see fig. 1), are a ZnS:LiF scintillator for thermal neutrons positioned in the direct beam, a mirror positioned at 45° relative to the scintillator to reflect photons coming off the scintillator out of the direct beam path, a focusing lens to adjust the focus onto the scintillator surface, an MCP image intensifier (with P47 phosphor, 106 amplification, 6 μm MCP pore size, 8 μm pitch, a collimation ratio L/D = 50:1, two stack MCP mounted in chevron configuration) including a relay lens, and the ASI optical camera TPX3Cam based on a Timepix3 readout chip. It should be noted that the principle experimental setup for the measurements is identical to that of “regular” neutron imaging detectors using scintillators [21], apart from the camera system. Moreover, for the measurements, the BOA beamline was essentially operated in its standard configuration with only the camera replaced by the intensified TPX3Cam camera. A photograph of the setup at the BOA beamline and a schematic of the setup are shown in fig. 1A and B, respectively.
As it is crucial for measurements using lenses that they are properly focused, the relay lens of the intensifier was focused using a low flux of photons hitting the light intensifier, such that the average size for a single photon on the detector resulted in the activation of ~3×3 pixels (average per photon event). The second focusing lens (between mirror and light intensifier, see fig. 1B) was adjusted using short acquisitions while observing and minimizing the radius of clusters of photons on the scintillator (neutron events) as illustrated in fig. 1C. Inspecting the events, the clusters of photons can be identified by the naked eye. Additionally, single photons without a direct correlation in position or time with other photon events were identified and as will be shown in the following section can be attributed to the scintillator afterglow. This result was as expected since each neutron interaction excites the scintillator with the consequential emission of a large number of photons (>> 10) within a short time-frame (~1 ms). In order to characterize the spatial and temporal parameters of the events, ToF measurements were performed to correlate the spatial information of the photons with respect to their timing for each neutron interaction.
1.2. Neutron Time-of-Flight measurements of the detector response
In a first step to entangle the spatial and temporal information of the photons, measurements were performed at the BOA beamline with the instrument operating in ToF mode [18]. Data with a large FoV (15×15cm2) were acquired using a single disc chopper. The chopper disc to detector distance was set to 6.45 m. The wavelength resolution, ,in this configuration resulted in ~1% (at 4 Å), calculated using the chopper disc radius of 200 mm, the slit of 1.8 mm with a length of 30 mm from the outer diameter of the disc and the 25 Hz rotation frequency.
The number of interactions on the scintillator was measured in this configuration in photon and neutron event-mode for a total acquisition time of 20.4 h and is shown in fig. 2A. It should be noted that the term “flux” in fig. 2 does not correspond to the actual flux of photons or neutrons at BOA in the configuration as shown in fig. 1, but is correlated to the number of detected events in photon or neutron event-mode. For that matter, for each detected neutron, about ten photons were detected on average by the TPX3Cam.
Two different areas on the detector with lower and higher incident intensities, denominated in fig. 2A by “dark” and “bright” regions, respectively, were selected for a direct comparison of the detector response under different neutron fluxes. The lower and higher intensity regions result from the inhomogeneity of the incident neutron flux and will be addressed in the following section (see fig. 3B). For both regions, it can be observed that at a ToF of t ~5 ms (or ~4.5 Å) the flux reaches its maximum, whereby t > 30 ms (corresponding to wavelengths > 25 Å) the profile flattens. It is assumed that in the flat region mostly photons from the long-lived fluorescent lifetimes are emitted, also known as scintillator afterglow. Therefore, as a good approximation, single photons t > 30 ms are emitted randomly from the scintillator surface, not contributing to the photon clusters as discussed earlier. As a result, with the here presented neutron event-mode, the afterglow can be essentially eliminated via rejection of these events. This can be observed by the significantly lower intensities of neutrons at t > 30 ms relative to the peak of the profile at ~5 ms in fig. 2A. More specifically, fig. 2A shows that independent of the intensity of the neutron event-based data for bright or dark regions on the detector area (localized higher or lower neutron flux), the measured flux drops to the same value of ~30 neutrons/cm2/s at t > 30 ms and can be attributed to a constant background. This background is typical for neutron ToF measurements [22] due to the experimental cave made from concrete filling up with essentially a "neutron gas”, which to a good approximation can be assumed constant throughout the entire wavelength range. The photon event-based data in contrast shows a significantly different “stationary” value at t > 30 ms for bright and dark detector regions, resulting from the fluorescent afterglow effect that is depended on the localized exposure of the scintillator by neutrons, inducing the consequently emitted photons.
Assuming that for t → ∞ the signal on the detector is purely noise, the signal-to-noise (S/N) ratio for the measured photon and neutron event-based fluxes averaged over the entire detector area are shown in fig. 2 B). Comparing the neutron and photon profiles within the brightest regions on the detector, the S/N improved by up to a factor of 7.5, with a S/N of up to 300.
1.3. Large field-of-view imaging (15×15 cm2 active area)
Using the setup described in the previous section, several objects were imaged with a 15×15 cm2 FoV at 6.45 m distance from the chopper to the scintillator and a total acquisition time of 9.27 h. A photograph of the objects, in particular a cruciform steel sample for biaxial tensile load studies, positioned at 10 cm distance from the scintillator is shown in fig. 3A, while fig. 3B shows an open beam image taken without sample. Fig. 3C - E show normalized neutron radiographs of a sample constructed from different event-modes. Inspecting the radiographs, the photon (fig. 3D) and neutron (fig. 3E) event-based images show incrementally improved image qualities when going from raw (fig. 3C) to photon image reconstruction in a first step and from photon to neutron event image reconstruction in a second step. The improved image quality becomes particularly apparent in the enlarged sections of the images as shown in fig. 3C - E below the full FoV.
Vertical stripes were observed in all radiographs and can be attributed to a very non-uniform incident beam with flux depleted regions that remain as artefacts despite open beam normalization. The non-uniformity is caused by a multi-channel bender neutron optic upstream of the beam exit into the instrument. Through the utilized beam limiting slit at the chopper in the instrument the bender channel structure is imaged, causing a horizontal variation in intensities by more than an order of magnitude within the measured FoV as can be observed in fig. 3B. Nonetheless, despite the large variations in intensity, most of the detector area was reasonably well normalized.
The modulation transfer function (MTF) using the slanted edge method [23] applied to the slanted edges provided by the steel cruciform (fig. 3) for the different event-modes was computed and is shown in fig. 4. Using the raw event-data, the resolution at 50% MTF equates to 0.508 lp/mm or 984 μm. Applying the CoM for individual pixels to detect photon events, the resolution at 50% MTF improved to 0.703 lp/mm or 711 μm using ½ pixel-pitch in the image reconstruction. By furthermore applying the CoM for individual photon events for neutron detection as described earlier, the resolution at 50% MTF improved to 1.289 lp/mm or 388 μm using ½ pixel-pitch in the image reconstruction. While the resolution increased significantly going from raw to photon to neutron event-mode, it can be shown that by furthermore reducing the effective sub-pixel size to ¼ pixel-pitch the resolution continues to increase, with an improved resolution at 50% MTF of 1.466 lp/mm or 341 μm using the neutron event-based data. However, by splitting the events into smaller and smaller bins, the standard deviation for each sub-pixel increases, limiting in this case the practical sub-pixel size to ~¼ pixel-pitch at which the MTF can still be still be computed.
1.4. Small field-of-view imaging (4×4 cm2 active area)
For high spatial resolution measurements, the chopper was removed and a white-beam radiograph of a test pattern (Siemens star made from Gd) with an active FoV of 4×4 cm2, as shown in fig. 5, was recorded at the BOA beamline. The total acquisition time was 2.2 h. Using the integrated raw event-data, a radiograph at intrinsic resolution of the sensor is shown in fig. 5A. Inspecting the line-pairs of the test pattern, a resolution close to the intrinsic resolution of the sensor for that FoV with ~3 lp/mm or ~166.7 μm was observed. For CoM event-based neutron image reconstruction the resolution improved to ~10 lp/mm or ~50 μm. This yields an increase by a factor of 3.3 in spatial resolution comparing raw event-data with neutron event-based image reconstruction. Inspecting the regions around the scintillator (rounded square region of 3×3cm2 with one edge broken-off), it is noticeable that essentially no counts were measured in the active FoV that was not covered by the scintillator material. For the computation of transmission values in these regions, this resulted in a division by zero. However, to visualize this, the values were set to zero, as can be observed by the black background in fig. 5B surrounding the scintillator. In contrast, fig. 5A shows transmission values around 0.5 in these regions, resulting from the division of the photon background (reflected photons within the scintillator box, see fig. 5A) with and without the sample in the beam. As discussed previously, this background is similar to the afterglow scenario in the sense of randomized single-photons that do not contribute to photon clusters, but still make their way to the detector. In the data processing, these events are discarded. The result is that light reflections are essentially eliminated and do not impact the computed transmission values.
While the acquisition time of 2.2 h seems relatively long for a radiograph that can resolve features down to 50 μm at a 4×4 cm2 FoV, the aim was to collect sufficient data for an in-depth analysis of the events on the scintillator. It should be noted that to avoid any saturation of the MCP or TPX3Cam, the incoming neutron beam was collimated down to 1×1 cm2 using Cd-slits to reduce the hit-rate on the detector to be well below the saturation limit to avoid loss of photons, whereby the setup would have allowed to open the slits to 4×4 cm2, potentially increasing the flux by a factor of 16. At present, the Timepix4 chip [24] is already being developed and will potentially offer an increase in hit-rate roughly by an order of magnitude compared to Timepix3 [3], with each new generation of Timepix chip presumably increasing the maximum hit-rate further.