SingleVolume Neutron Scatter Camera for HighEfficiency Neutron Imaging and Spectroscopy
Abstract
Neutron detection provides an effective method to detect, locate, and characterize sources of interest to nuclear security applications. Current neutron imaging systems based on doublescatter kinematic reconstruction provide good signal vs. background discrimination and spectral capability, but suffer from poor sensitivity due to geometrical constraints. This weakness can be overcome if both neutronproton scattering interactions are detected and resolved within one large contiguous active detector volume. We describe here a maximum likelihood approach to event reconstruction in a singlevolume system with no optical segmentation and sensitivity to individual optical photons on the surfaces of the scintillator. We present results from a Geant4based simulation establishing the feasibility of this singlevolume neutron scatter camera concept given notional performance of existing photodetector and readout technologies.
keywords:
fast neutron imaging, neutron scatter camera, nuclear security applications, neutron detection1 Introduction
Fissionenergy neutrons are a sensitive and specific signature of special nuclear material, due to their low and stable natural backgrounds, their penetrating nature, and the scarcity of benign neutronemitting materials. As such, neutron imaging has the potential to be an important tool in a range of nuclear security applications, including arms control treaty verification, emergency response, cargo screening, and standoff search. In particular, neutron scatter cameras (NSC) have become an established technology for fissionenergy neutron imaging and spectroscopy Ryan_1993 ; Vanier_2006 ; Marleau_2007 ; goldsmith_2016 ; PoitrassonRiviere_2014 . Based on the kinematic reconstruction of neutron double scatters, NSC systems provide excellent eventbyevent directional information for signaltobackground discrimination, reasonable imaging resolution, and good energy resolution when compared to competing imaging technologies PoitrassonRiviere_2014 ; Mascarenhas_2009 ; Brennan_2011 ; Weinfurther_2017 .
The operational principle of an NSC is based on the neutron elastically scattering on hydrogen at least twice in the detector’s active material, typically an organic scintillator. The incoming kinetic energy is obtained from
(1) 
as the sum of the energy deposited in the first scatter, which is assumed to generate a recoiling proton, plus the scattered neutron energy . This in turn is calculated from
(2) 
using the measured distance and time between the first two scatters. These quantities also allow for a calculation of the scattering angle defining a cone containing the incoming neutron direction:
(3) 
Most existing NSC systems consist of arrays of spatiallyseparated organic scintillator cells in order to isolate in two distinct cells the neutron’s two elastic scatters Ryan_1993 ; Vanier_2006 ; Marleau_2007 ; goldsmith_2016 ; PoitrassonRiviere_2014 . This requirement, however, has poor geometrical efficiency, constituting the primary drawback of current NSC systems.
The design concept of a singlevolume neutron scatter camera (SVSC), in which both neutron scatters occur in the same large active volume, addresses this drawback. Given that the interaction length of fissionenergy neutrons in organic materials is a few cm, the detector’s principle relies on the ability to resolve two proton recoils in a contiguous scintillator volume at spatial and temporal separations of order and respectively. This reconstruction of the neutron interaction locations must be based on the arrival time and position of the scintillation photons at the boundaries of the active volume. While traditional photomultiplier tubes (PMTs) cannot provide the necessary spatial and temporal resolution, recent advances in photodetector (PD) technology have made the approach possible Krizan_2014 . Photodetectors based on microchannel plate (MCP) electron multipliers Wiza_1979 inherently provide the spatial and temporal photon detection resolution that makes them attractive for singlevolume event reconstruction. These MCPPMTs have found application in a wide variety of areas ranging from medical imaging Kim_2012 to neutrino detection Li_2016a ; Back_2017 . Commercially available photodetectors have twodimensional pixels and subns timing planacon_datasheet , while models under development will potentially provide large area photocathode coverage at relatively low cost Adams_2015 . In this work, we investigate the feasibility of the SVSC concept using Monte Carlo simulation studies of a system with PD performance motivated by available MCPPMTs.
Throughout this paper, the term “singlevolume” means that the active detection volume is compact and essentially contiguous, rather than spread out spatially in discrete well separated volumes (cells, planes, etc.). We refer to an individual neutron scatter, most usefully neutronproton elastic scattering, as an “interaction”, and to a neutron history, which may include multiple interactions, as an “event”. In an experimental context, an event can mean the data collected by a single trigger, which could include information from zero (noise), one, or more (pileup) particles. “Event reconstruction” is an algorithm that uses the observable optical photon information to determine the locations and times of neutron interactions in the event. Many such reconstructed events can then be combined to reconstruct a neutron image and energy spectrum.
This article is structured as follows. We begin in Sec. 2 with the general motivation and designagnostic arguments for a singlevolume neutron scatter camera, as well as a survey of several possible design classes. In Sec. 3, we describe the “direct reconstruction” SVSC concept as simulated in this paper, emphasizing how the design choices are enabled by existing technology. The data processing and reconstruction steps required by this concept are also detailed in Sec. 3. The rest of the paper is dedicated to a demonstration of the feasibility of the SVSC concept based on Monte Carlo simulation data. Details of the simulated model and the corresponding physical behavior of the system directly obtained from the simulation are presented in Sec. 4. The results of the data processing stage applied to the simulated data, which will cover the event reconstruction and the source spectral and image reconstruction, are presented in Sec. 5 together with an analysis of the instrument’s limitations.
2 General considerations for a singlevolume neutron camera
We begin with a discussion of the motivation for and approaches to single volume neutron imaging, without reference to a specific concrete detector system design.
2.1 Advantages of singlevolume neutron imagers
The motivation for research into a singlevolume doublescatter neutron imager rests on the potential for significant improvements over cellbased scatter camera designs in two important aspects: imaging efficiency and compactness.
A singlevolume camera addresses two of the important efficiency limitations of the NSC: the requirement of no more than one neutron scatter in the first detector element (in order for a onetoone relationship between fractional energy deposited and scattering angle to be valid), and the geometrical efficiency for the scattered neutron to pass through a second detector element. Assuming the ability to resolve scattering events at separation in the single volume, an estimate of the improvement in efficiency for doublescatter neutron detection, relative to the NSC of Brennan_2011 , is about an order of magnitude.
To obtain this estimate, an MCNP Pozzi_2003 simulation was performed: a fission neutron pencil beam was fired at the center of one of the middle detectors in the conventional NSC front plane, and at the center of a volume of scintillator. The NSC was simulated with a standard plane separation of . The fraction of events satisfying simple requirements were tallied: first, that two interactions occur—as a function of the distance between the first two neutron scatter interactions for the singlevolume system, and with a front and rear plane requirement for the current NSC; second, that both interactions were hydrogen scatters (carbon recoils are typically too quenched to observe); and third, that both interactions deposited at least , an optimistic threshold. In Fig. 1, the results are shown for the singlevolume system as a function of the distance between the first two interactions. The corresponding fractions for the cellbased NSC are 0.026 for the first two interactions in the front and rear planes; 0.013 for both scatters producing proton recoils; and 0.008 for both protons having at least . The SVSC yields a significant improvement over the conventional NSC in potential efficiency when the minimum separation between interactions is low. For a minimum interaction spacing, the number of potentially detectable events is slightly over an order of magnitude higher in the SVSC than in the current NSC. Even at minimum interaction spacing, the fraction of usable neutrons is comparable to the current NSC. Further study will be needed to understand the relative power of events with small separation in a SVSC, since the energy and imaging resolutions will degrade with smaller separation distances.
Additionally, the overall detector footprint is more compact and lighter: about , with expected weight about for the same scintillator volume as the 32element SNL NSC, which takes up a volume of approximately and weighs more than . Easier transport and deployment is a key benefit for some applications. Moreover, a more compact detector allows to get closer to the source, which in turn increases detection rate according to the inverse square of the distance to the source, and improves spatial resolution at the object for a given angular resolution.
2.1.1 Four singlevolume imager concepts
Singlevolume neutron imaging is predicated on the ability to use the information encoded in the optical photons emitted by the scintillator to reconstruct the scattering history of the neutron in the active volume. Even if it is granted that from an information theory perspective there is in fact sufficient information retained by the optical photons to do this, it is still one of the key technical challenges of this research to develop detector designs and algorithms that accurately (and tractably) do it. Four general methods that have been considered for building a singlevolume neutron imager are as follows.
 Direct reconstruction

This method is conceptually the simplest and yet technically the most challenging. The scintillator is a single monolithic volume, and optical photons are detected at the boundaries of the volume. The neutron interactions are reconstructed based on the information in the detected photon positions and times, accounting for the isotropic photon emission from the scintillator, the scintillator pulse shape, photon speed and attenuation in the scintillator, photodetector efficiency, etc. This paper documents studies of this concept.
 Optical segmentation

The “single volume” of scintillator is formed from multiple pillars of scintillator, separated by reflective boundary layers. The goal is to trap the optical photons from an interaction in one pillar, reading them out from both ends. For each interaction, the location along the long dimension of the pillar is determined by the amplitude ratio and/or time difference of the light observed at each end, while the transverse location is determined by which pillar the light is observed in. When two interactions occur in different pillars, their locations and times can be independently determined, constituting an event reconstruction. Simulation studies toward a realization of this concept are detailed in Weinfurther_2017 .
 Optical coded aperture

This method, like direct reconstruction, uses a single uninterrupted scintillator volume, but instead of directly detecting and reading out the optical photons on the boundary of the volume, they pass through a coded aperture on each side. This has the advantage of introducing a highfrequency spatial component to the detected photon distribution, which eases the spatial reconstruction of interaction locations. Disadvantages include the loss of a significant fraction (typically 50%) of the photons to the mask, and the presence of light guide material, which is necessary to allow the photons to spread before and after the coded mask, but also results in an inactive scattering/attenuating layer that neutrons must penetrate before reaching the active volume. A similar approach has been studied for gamma Compton imaging Ziock_2015 .
 Optical lattice

This is similar to the optical segmentation approach, but aims to create virtual “pillars” in all three dimensions. Small cubes of scintillator are separated by very thin air gaps in a threedimensional array. The scintillatorair boundaries create total internal reflection (TIR) for photons hitting them at larger than the critical angle. Some photons experience TIR in each of the three dimensions of the lattice, and are therefore confined to a single “pillar” or row of cubes. Light therefore arrives as a bright spot on all six sides of the array, not just two. This approach has been studied for neutrino detection Lane_2015 .
It should be noted that the reconstruction uncertainties, and consequently the directional and spectral resolving power, of each event in singlevolume designs depend significantly on the design details, component performance (scintillator time profile, photodetector efficiency and resolution, etc.), and reconstruction algorithms used. On the one hand, a singlevolume design will be able to determine the location and time of neutron interactions with better resolution than the NSC. On the other hand, we expect reduced precision in the measurements of deposited energy when disentangling the light output from multiple interactions, and the short lever arm between the two neutron interactions also decreases imaging resolution with respect to the cellbased design. In the remainder of this paper, we focus on further exploring the direct reconstruction technique.
3 Singlevolume scatter camera design
We now proceed to specify a concrete detector system design and reconstruction algorithm for full simulation studies of a singlevolume imager.
3.1 SVSC system as simulated
The direct reconstruction SVSC design concept, schematically shown in Fig. 2, consists of an organic scintillator volume with each side optically coupled to and fully covered by fasttiming highgain MCPPMTs. These photodetectors play the central role of registering the arrival time and position of the isotropically emitted photons produced by particle interactions in the scintillator volume. As already mentioned, current MCPPMT technology, either available or under development, can provide photon timing resolution on the order of and mmscale spatial photon detection resolution. The simulation studies presented in Sec. 4 show the impact on the event reconstruction due to such finite MCPPMT resolution values by comparing results with idealized photodetectors able to provide exact time and spatial photon coordinates, and demonstrate the detector’s spectral and image reconstruction capabilities even with existing MCPPMT technology. Another enabling consideration of the SVSC concept is ability of high speed waveform sampling required to take advantage of the intrinsic time resolution of MCPPMTs, and especially to disentangle multiple photons arriving in one photodetector pixel within a timescale of nanoseconds. This is possible using switched capacitor array sampling, such as with the Domino Ring Sampler chips, which provide up to sampling frequency Ritt_2010 ; Bitossi_2016 ; Oberla_2014
When selecting the active material type, fast scintillation decay time represents the main driver, even more important than the ability to discriminate neutrons and gammas via pulse shape discrimination (PSD). Narrow scintillation pulses allow for better event timing reconstruction, especially when reconstructing events with several neutron interactions. Although PSD is desirable for this system, it may be superfluous since we already require the ability to reconstruct the time of flight of the detected particle, which provides neutron/gamma discrimination. High light output is also desired, as photon statistics is one limiting factor in the event reconstruction described below, although we note that the readout and signal processing are simplified if the photon occupancy per MCPPMT pixel is low. Here we simulate EJ232Q ej232q_datasheet , a quenched plastic scintillator with pulse width of less than but only 19% anthracene light output, a choice that emphasizes the pulse width over photostatistics. The scintillator size is set to a cube in order to guarantee a high probability of double neutron scatters, in accordance with the interaction length for fissionenergy neutrons. In general, the choice of active volume size will be dictated by a variety of factors, like detection rate and efficiency requirements, hardware constraints and applicationspecific tradeoffs, which go beyond the scope of this paper.
3.2 Reconstruction algorithms
Data processing and analysis are essential elements in the design and operation of the SVSC. They are accomplished in three separable analysis stages: singlephoton isolation in each individual pixel, neutron doublescatter event reconstruction, and the final image and spectrum reconstruction.
Singlephoton isolation is the process of assigning an arrival time and location to each detected photon using the digitized raw traces from each MCPPMT channel (or pixel). The specifics of this part of the analysis are highly coupled to the experimental behavior of the photodetector and electronics. For this simulation study, since the goal is to determine feasibility of the concept, we make the optimistic assumption that all detected photons can be individually identified.
For event reconstruction, we employ the direct reconstruction of particle interactions using the arrival position and time of isotropically emitted photons at the scintillator surface. Multiple methods have been explored; here we use the unbinned maximum likelihood (ML) technique. The likelihood function is maximized for a given event when the fixed set of observations (the list of detected photon coordinates and times) has the highest probability of occurring as the event parameters (the list of neutron interaction locations and times) are varied. Put another way, the ML algorithm finds the neutron trajectory that would be most likely to produce the observed data. We use the extended maximum likelihood Barlow_1990 ,
(4) 
where is the number of photons detected from neutron interaction and , is the number of detected photons, is the number of neutron interactions assumed, represents the detected position and time of photon , and is the probability to observe photon from interaction . Each photon could have come from any of the neutron interactions, so the probability to observe each photon is summed across interactions; each photon is an independent observation, so the photon probabilities multiply. Looking more closely at the probability function, we have
(5) 
The first term accounts for solid angle between the location of emission (interaction ) and the location of detection for photon . The angle between the line connecting the two locations and the normal of the photodetector surface is , and is the distance between the two locations. The second term accounts for attenuation in the scintillator material, where is the attenuation length. In the third term, is the pulse shape of the scintillator, where the time difference between the photon emission and detection has been corrected for the flight time of the photon ( is the speed of light in the medium). The likelihood maximization itself is performed using the SIMPLEX and MIGRAD algorithms as implemented by the MINUIT package in the ROOT software James:1994vla .
This likelihood maximization results in an estimate of the position, time, and intensity of each neutron interaction, and an associated variancecovariance matrix (via inverting the Hessian matrix at the loglikelihood maximum). From these fitted quantities, we then derive , , and the cone axis direction, as shown in Fig. 2, and their uncertainties. Finally, the kinematic equations Eqs. 1–3 are used to determine the ultimate quantities of interest, namely and the cone opening angle . In all cases, the fit errors are propagated to the derived results via the standard linearized approximation, including correlations.
Finally, image reconstruction is performed from many such neutron events using a Maximum Likelihood Estimation Maximization (MLEM) algorithm Shepp_1982 . We use a listmode implementation of MLEM, in which a simplistic system response is calculated for each observed event, under the assumption that the reconstructed cone axis and opening angle are accurate within uncertainties.
4 Simulation
A Monte Carlo model of the SVSC detector system was constructed using Geant v4.10.01.p02 Agostinelli_2003 . The scintillator material as described in Sec. 3.1 is centered at the coordinate origin. A thick acrylic light guide representing an optical coupling encloses the scintillator, which is then completely surrounded by a quartz layer of very short optical absorption length representing the MCPPMT. An isotropic pointlike neutron source is placed at from the origin along the axis.
The simulation uses the physical attributes listed in Table 1. Neutron interactions are modeled by the Geant4 NeutronHP crosssection package, which includes high precision hadronic models for energies below , appropriate for fissionenergy neutrons. Scintillation photons from electron recoils are generated according to the manufacturerreported scintillation efficiency for EJ232Q with 0.5% benzophenone concentration. Due to lack of experimental data on proton quenching in this scintillator, we use the energydependent proton light output measured for a similar plastic scintillator in Pozzi_2004 , and we model the carbon recoil light output according to Batchelor_1961 . The photon generation time profile follows the manufacturerreported exponential rise and decays times, shown in Table 1, with a reported scintillation pulse width of . The optical interfaces between the scintillator, the light guide and the photodetector are assumed to be perfectly smooth, and photons are either reflected or transmitted but are not absorbed in the materials interface. A wavelengthdependent refraction index is used for the acrylic light guide Kasarova_2007 . Optical photon detections per event are recorded over a window. In order to account for the photocathode quantum efficiency, photons absorbed in the photodetector volume are randomly saved with 25% efficiency.
Attribute  Value 

SC cube size  
SC Efficiency  
SC pulse rise time  
SC pulse decay time  
SC absorption length  
SC refraction index  1.58 
LG refraction index  1.49–1.59 
PD quantum efficiency  0.25 
PD absorption length  
PD refraction index  1.57 
PD time spread  
PD pixel size 
The left plot of Fig. 3 shows the histograms of the number of neutron elastic collisions within the scintillator volume, normalized by the number of primary neutrons that interact in the detector. When a threshold is imposed on the neutron deposited energy—which corresponds to about of proton light output in EJ232Q or, equivalently, to only about 30 detected photons—the number of abovethreshold neutron collisions per event is reduced drastically. This indicates that the transfer of most of the neutron’s energy to the scintillator is usually done in few elastic scatters.
The middle and right histograms of Figure 3 show the histograms of the time and distance between the first two neutron collisions. These quantities, along with the energy transferred in the first collision, contain the information from which the neutron’s incoming energy and direction can be extracted. For zero threshold, the mean distance between the first two neutron collisions is about , increasing to for collisions above , and the mean time remains in both threshold cases.
The arrival time and spatial coordinates of the individual optical photons from all neutron collisions, regardless of deposited energy, constitute the input data for the event reconstruction algorithm described above. In order to include the effect of finite timing resolution in the MCPPMT and electronics, the photon arrival time was Gaussian smeared with standard deviation equal to , which is informed by the transit time spread in existing detectors and bandwidth of available electronics planacon_datasheet ; Grigoryev_2016 . The photon twodimensional coordinates along each photodetector plane—e.g., the coordinates of the photodetector at —were also Gaussian smeared. For the coordinates smearing, we use a standard deviation equal to , which matches the standard deviation of a rectangular distribution with width representing the pixel size of the Planacon XP85012.
In the following two ways, the current study takes advantage of information that is known in simulated events but not in real data, and therefore potentially overestimates the system performance. First, we constrain the total number of neutron collisions, , to be the simulated number of neutron scatters depositing more than . In real data, it should be possible to estimate the number of energetic collisions in the scintillator volume from the spatially integrated photon detection time profile combined with the timesliced spatial distributions of photon detection locations. We note that there is no obviously “correct” value for : counting real but lowenergy collisions into adds irrelevant degrees of freedom to the MINUIT minimization, increasing the probability of incorrect convergence of Eq. (4). A process of iterative minimizations with ranging over most likely values according to the observed signals could be part of the event reconstruction algorithm when working with real data.
Second, in this study, small deviations from the interactions’ actual location coordinates (), time (), and number of emitted photons (50) are used to set the initial guesses in the minimization algorithm. (Events for which the minimization algorithm does not move from these initial parameter values are rejected as failed fits.) In the limit of a globally convex likelihood function, this procedure does not produce a bias, as the fit will converge to the same global minimum regardless of the initial position. However, in the case that the function has multiple local minima, this procedure could bias the results by overestimating the accuracy of the initial guess that could be achieved. Again, in real data, initial guesses could be extracted from rough analysis of the photon detections, but we leave the development and demonstration of such an algorithm to future work.
Idealizations in our simulated model could play a more relevant role in the reconstruction algorithm success. Complete wavelength dependence of optical processes is not included. The optical interfaces between the scintillator, the light guide and the photocathode are assumed to be perfectly polished surfaces with nearly matching indexes of refraction, and thus, most simulated photons are transmitted and absorbed in the photocathode. However, a realistic detector will likely suffer from increased photon scattering in the optical surfaces due to imperfect optical coupling and surface roughness. Any scattered photons that are later detected by the photodetectors will not satisfy the reconstruction’s assumption of being directly emitted from interaction locations, and thus, will affect the event reconstruction success. Additionally, the detected photon count of a real detector is expected to be smaller than that of our simulated model due to incomplete photocathode coverage and other unmodeled optical absorption effects. Furthermore, we assume in simulation that each photon can be individually resolved and accurate arrival time can be assigned, but the overlapping of their corresponding singlephotoelectron waveforms, already discussed in Sec. 3, will affect the photon count and timing per photodetector pixel. Finally, here we simulate and reconstruct only neutron events. In reality, some gamma interactions will be misidentified as neutrons and contribute erroneously to a neutron image.
Due to the assumptions and limitations of the simulation model described above, the event reconstruction and the image reconstruction results presented in the next section should be understood as an upper limit on the performance of the SVSC detector concept. Put another way, we aim here to establish feasibility of the SVSC design, rather than to accurately predict performance of a system.
5 Simulation results
In order to quantify the event reconstruction performance, we define the “true” variables in a given event as the positions, times, and intensities of the first two np scatters in the simulation, as well as the higherlevel quantities (e.g. or ) calculated from those. Note that in some classes of events, such as those with an initial carbon scatter before two hydrogen scatters, the “true” reconstructed cone may not include the incident direction since the kinematic assumptions are not valid. But the goal of the event reconstruction is to determine the attributes of the proton recoils, so we evaluate it on that basis.
For each variable, we construct resolution plots, i.e. histograms of the difference between the reconstructed and the true variables; these are plotted in Figs. 4–7 and will be discussed below in turn. The differences are denoted with the symbol followed by the variable symbols defined in Fig. 2. In these plots, the photon arrival time and coordinates that entered the reconstruction algorithms have been smeared as described in the previous section, and only events for which the likelihood maximization converged are included. Each histogram is plotted with a fitted Gaussian function, whose mean () and width () are given in the plot as well as in Table 2.
Smeared  Exact  

(mm)  
(mm)  
(mm)  
(ns)  
(mm)  
(mm)  
(mm)  
(ns)  
(mm)  
(ns)  
()  
()  
We observe first that the distributions are not well described by a single Gaussian. We attribute this to two effects. First, each distribution includes a large diversity of different types of events, with variability in the interactions’ distance to the photodetectors, in the distance between interactions, in the energy allocation among the interactions, in the total number of photonemitting interactions, etc. As such, even if event reconstruction resulted in accurate and normal uncertainties for each type of event, these histograms would not necessarily be Gaussian in shape since they are effectively composed of multiple distributions, each with its own standard deviation. Second, some events are reconstructed with an erroneous topology. For example, if the first two np scatters are very close in space and time, they may be reconstructed as a single interaction, and the third scatter reconstructed as the second. This second type of error will lead to very long nonGaussian tails in the distributions.
Still, the widths of the histograms, expressed in terms of the standard deviations of the corresponding Gaussian fits, provide a measure of the resolution achievable with the event reconstruction algorithm, while the Gaussian means resulting in values significantly shifted from zero represent biases in the reconstruction. The first pair of columns of Table 2 list the values of and for all plotted variables where the smeared time and coordinates were input to the algorithm, while the second pair gives the results when the exact and , before the Gaussian smearing, were used. With a similar column breakdown, Table 3 contains efficiency results, including the fraction of incident neutrons that contain a given number of proton recoils above a threshold, and the fraction of those for which the reconstruction is “successful”, defined as the fraction of events for which the likelihood maximization converges, with the added restriction that the reconstructed locations and times of the first two interactions are within and of the true values, respectively. The latter condition removes some, but not all, of the events reconstructed with an erroneous topology. As in Table 2, these values are given for the reconstructions with and without resolution smearing.
# proton recoils  Fraction of incident  Reconstruction success rate (%)  

fission neutrons (%)  Smeared  Exact  
2  
3  
4 
Figures 4–5 show the histograms for the locations, times, and detected photon count of the first and second abovethreshold proton recoils: (, , , , ) and (, , , , ), respectively. These histograms represent the resolutions of the individual interaction variables directly obtained from the event reconstruction algorithm. The reconstruction resolution for the first interaction is about in each spatial dimension and in time. The reconstruction of the second interaction is somehwat less accurate, as indicated by larger widths of the corresponding histograms. One possible cause is that belowthreshold collisions tend to occur later in the neutron trajectory as its kinetic energy decreases, and photons emitted from those belowthreshold collisions will be more likely assigned to the second reconstructed interaction. Added to this is the fact that the energy deposited in the second interaction is on average smaller than the energy deposited in the first interaction, producing a smaller second photon pulse which provides relatively less information to the reconstruction algorithm.
The and plots show a more significant bias for these quantities, on the order of their resolution. The reconstructed quantities, defined as the number of photons detected from interaction , are simply the from the maximized likelihood (Eq. (4)). The true quantities are the expected number of photons generated in that interaction (i.e. before Poisson fluctuation) multiplied by the photodetector quantum efficiency. Therefore, this definition of the true values leaves out both Poisson fluctuations and optical attenuation in the scintillator; the latter effect is consistent with the sign of the observed bias. Note that this bias propagates to the proton energy calculated from this number of photons, and ultimately therefore to the neutron energy and the scattering angle . To address the bias, we plan to recast the likelihood to depend on the emitted rather than the detected number of photons from each interaction, which more directly relates to the deposited energy, but requires implementing a nontrivial efficiency integral.
The histograms of the distance and time between the first and the second interaction, plotted in Fig. 6, are more relevant quantities for the event reconstruction, as the individual positions and times do not appear in the kinematic equations. The cores of these distributions are thought to represent well reconstructed events and show Gaussian standard deviation values of and . These quantities are derived respectively from combining the coordinates and times of the two interactions, which are not independent variables but instead are correlated by the reconstruction algorithm. Therefore, the asymmetric nonGaussian tails of the and histograms reveal correlations in the misreconstruction of the individual interaction variables, and their impact on the reconstruction algorithm limits is discussed below. Also included in Fig. 6 is the histogram of the energy deposited by the first recoiling proton, , calculated from . The quantities , , and constitute the inputs used to calculate the incoming neutron energy and scattering angle according to Eqs. 1–3.
In order to evaluate the effect of the event reconstruction in the final reconstruction of the source spectrum and image, the histograms of and are plotted in Fig. 7. We also plot the histogram of the angular difference between the reconstructed and the true vector defining the axis of the incoming neutron cone, since errors in the cone axis direction will also effect the image reconstruction along with errors in the cone opening angle . Note that is unique among these resolution plots in being an unsigned angular difference, and consequently we do not expect a symmetric distribution peaked at zero; therefore we do not include a Gaussian fit for . As for a Rayleigh distribution, we use the mode of the distribution as an estimate of the underlying angular uncertainty in each of two independent orthogonal components; this is the value given in Table 2. However, it should be noted that the distribution has a much longer tail than the corresponding Rayleigh distribution, so as in the case of the Gaussian widths, this value at best represents the “core” of well reconstructed events.
The most salient feature of these histograms is the broad peak for negative values of representing reconstructed scattering angles smaller than the true angles. These misreconstructed events correlate with the extended high tail of values, as well as other tails in the previous plots of reconstructed quantities. Further investigation indicated that these misreconstructed events are primarily observed when the distance between the first two neutron interactions in the scintillator is small, i.e. less than about . Figure 8 shows as well as plotted against the true value, illustrating that result. These events also of course tend to have small time difference between the first two interactions. One likely explanation for this observation, as mentioned above, is that when the first two interactions are close in space and time, they are reconstructed as a single interaction; the reconstruction algorithm then identifies a subsequent interaction or group of interactions as the second one.
Figure 9 provides a graphical comparison of the effect that smearing the photon time () and position coordinates have in the reconstruction quality, also showing the intermediate cases in which only time or only position resolutions were applied. The figure demonstrates that the time resolution is the main contributor to the error in event reconstruction; it also reduces the reconstruction success rate by 20%–30%. In contrast, the smearing in photon detection position within the pixel size of currently available MCP PMT technology has little impact on the event reconstruction. However, the pixel size drives not only position uncertainty but also photon occupancy per pixel, which does not affect these results but will be an important consideration for a hardware implementation.
Furthermore, when exact photon times are used, the negative tail in the histogram decreases in size, and events with the first two interactions as close as and can be successfully reconstructed. While this indicates that reduction in the photon detection time resolution beyond can potentially improve the SVSC performance further, it also establishes an intrinsic limit in event reconstruction resolution within this algorithm framework.
Finally, we consider the highestlevel reconstructed distributions, the energy spectrum and image. A reconstructed source energy spectrum is presented in Fig. 10, where it is compared with the “true” spectrum calculated using the known true values of the interactions’ location, time and deposited energy, as well as with the actual incident energy spectrum of the same set of events. The latter two distributions are quite similar albeit with some differences at low , but the reconstructed energy spectrum has two qualitative issues. First, for the accurately reconstructed events in the core of the spectrum, is biased low due to the absence of a correction for optical attenuation as discussed above. Second, there is a long tail of misreconstructed events with high ; in fact, fully 20% of the reconstructed events overflow this plot ().
The reconstructed image of Fig. 11 shows good agreement with the simulated source location, indicated by a black cross in the lower zoomedin view. In the case of the image, misreconstructed events tend to contribute randomly, which reduces their importance in determining a source centroid. This image was reconstructed with listmode MLEM under the assumption of a simple system response, in which for example the two reconstructed interactions are always taken to be the first two neutron interactions. This method produces qualitatively better results than simple backprojection. However, we expect that further improvements could be obtained in the future with a more complete system response function, including better treatment of experimental resolutions and alternate event topologies, such as, for example, including the possility of an initial carbon scatter. Quantification of the SVSC image quality for point and extended sources will be done in future work.
6 Conclusions
We have described the theoretical advantages of singlevolume doublescatter neutron imaging, namely a high imaging efficiency and a compact form factor. We also demonstrated the feasibility of the direct reconstruction SVSC concept via Monte Carlo simulation and reconstruction of a neutron point source.
For the simulated detector system, the results indicate about 13.6% of incident fissionenergy neutrons result in at least two proton scattering interactions above , with about half of those events successfully reconstructed using the described event likelihood maximization algorithm. A core resolution of () in each spatial dimension and () in time is observed for the first (second) interaction, while a limitation of the current likelihood expression results in a bias in the reconstructed number of photons for each interaction. Resolutions of higherlevel quantities, such as the distance between the two interactions, have complicated distributions, indicating the presence of nontrivial correlations among the lowlevel quantitites.
Most of the poorly reconstructed events have small distance (and time) between the two first interactions; this is expected both because these events are more difficult to reconstruct and because, even for a given spatial/temporal resolution, they contain less neutron energy and direction information due to the reduced lever arm.
The response of the system to gammas has not been addressed here, but is an important part of ongoing simulation studies. In particular we would like to understand the ability of the system to discriminate neutron from gamma events based on interaction timing in the absence of pulse shape discrimination.
Physical construction, calibration, and operation of the simulated system presents several technical challenges, and the next desired step is demonstration of the SVSC concept in a laboratory prototype system.
Acknowledgements
We thank Glenn Jocher and Kurtis Nishimura for sharing knowledge and wisdom from their experience with a closely related detector system. We thank Marc Ruch for his software implementation of listmode MLEM for doublescatter imaging, as well as for various fruitful discussions on compact neutron imaging.
Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DENA0003525.
References
References

(1)
J. Ryan, K. Bennett, H. Debrunner, D. Forrest, J. Lockwood, M. Loomis,
M. McConnell, D. Morris, V. Schönfelder, B. Swanenburg, W. Webber,
Comptel
measurements of solar flare neutrons, Advances in Space Research 13 (9)
(1993) 255 – 258.
doi:10.1016/02731177(93)90487V.
URL http://www.sciencedirect.com/science/article/pii/027311779390487V 
(2)
P. E. Vanier, Analogies between
neutron imaging and gammaray imaging, in: Proc.SPIE, Vol. 6319, 2006, pp.
6319 – 6319 – 8.
doi:10.1117/12.682717.
URL http://dx.doi.org/10.1117/12.682717  (3) P. Marleau, J. Brennan, K. Krenz, N. Mascarenhas, S. Mrowka, Advances in imaging fission neutrons with a neutron scatter camera, in: 2007 IEEE Nuclear Science Symposium Conference Record, Institute of Electrical & Electronics Engineers (IEEE), 2007, pp. 170–172. doi:10.1109/nssmic.2007.4436310.

(4)
J. E. M. Goldsmith, M. D. Gerling, J. S. Brennan,
A compact neutron scatter camera
for field deployment, Review of Scientific Instruments 87 (8) (2016) 083307.
arXiv:http://dx.doi.org/10.1063/1.4961111, doi:10.1063/1.4961111.
URL http://dx.doi.org/10.1063/1.4961111 
(5)
A. PoitrassonRivière, M. C. Hamel, J. K. Polack, M. Flaska, S. D.
Clarke, S. A. Pozzi,
Dualparticle
imaging system based on simultaneous detection of photon and neutron
collision events, Nuclear Instruments and Methods in Physics Research
Section A: Accelerators, Spectrometers, Detectors and Associated Equipment
760 (Supplement C) (2014) 40–45.
doi:10.1016/j.nima.2014.05.056.
URL http://www.sciencedirect.com/science/article/pii/S0168900214005889 
(6)
N. Mascarenhas, J. Brennan, K. Krenz, P. Marleau, S. Mrowka,
Results
With the Neutron Scatter Camera, IEEE Transactions on Nuclear Science
56 (3) (2009) 1269–1273.
doi:10.1109/TNS.2009.2016659.
URL http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=5076056 
(7)
J. Brennan, E. Brubaker, R. Cooper, M. Gerling, C. Greenberg, P. Marleau,
N. Mascarenhas, S. Mrowka,
Measurement
of the fast neutron energy spectrum of an source using a
neutron scatter camera, IEEE Transactions on Nuclear Science 58 (5) (2011)
2426–2430.
doi:10.1109/TNS.2011.2163192.
URL http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=6012498  (8) K. Weinfurther, J. Mattingly, E. Brubaker, J. Steele, Modelbased design evaluation of a compact, highefficiency neutron scatter camera, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 883 (2018) 115–135. arXiv:1710.06480v1, doi:10.1016/j.nima.2017.11.025.

(9)
P. Križan,
Overview
of photon detectors for fast single photon detection, Journal of
Instrumentation 9 (10) (2014) C10010–C10010.
doi:10.1088/17480221/9/10/C10010.
URL http://stacks.iop.org/17480221/9/i=10/a=C10010?key=crossref.dd2d276bf48d2aec459f867b12392922 
(10)
J. L. Wiza, Microchannel
plate detectors, Nuclear Instruments and Methods 162 (13) (1979) 587–601.
doi:10.1016/0029554x(79)907341.
URL http://dx.doi.org/10.1016/0029554X(79)907341 
(11)
H. Kim, C.T. Chen, H. Frisch, F. Tang, C.M. Kao,
A prototype TOF PET
detector module using a microchannel plate photomultiplier tube with
waveform sampling, Nuclear Instruments and Methods in Physics Research
Section A: Accelerators, Spectrometers, Detectors and Associated Equipment
662 (1) (2012) 26–32.
doi:10.1016/j.nima.2011.09.059.
URL http://dx.doi.org/10.1016/j.nima.2011.09.059 
(12)
V. A. Li, R. Dorrill, M. J. Duvall, J. Koblanski, S. Negrashov, M. Sakai, S. A.
Wipperfurth, K. Engel, G. R. Jocher, J. G. Learned, L. Macchiarulo,
S. Matsuno, W. F. McDonough, H. P. Mumm, J. Murillo, K. Nishimura, M. Rosen,
S. M. Usman, G. S. Varner, Invited
article: miniTimeCube, Review of Scientific Instruments 87 (2) (2016)
021301.
doi:10.1063/1.4942243.
URL https://doi.org/10.1063/1.4942243  (13) A. R. Back, J. F. Beacom, M. Bergevin, E. CatanoMur, S. Dazeley, E. Drakopoulou, F. D. Lodovico, A. Elagin, J. Eisch, V. Fischer, S. Gardiner, R. Hatcher, J. He, R. Hill, T. Katori, F. Krennrich, R. Kreymer, M. Malek, C. L. McGivern, M. Needham, M. O’Flaherty, G. D. O. Gann, B. Richards, M. C. Sanchez, M. Smy, R. Svoboda, E. Tiras, M. Vagins, J. Wang, A. Weinstein, M. Wetstein, Accelerator Neutrino Neutron Interaction Experiment (ANNIE): Preliminary results and physics phase proposal, Tech. rep., Iowa State University (2017). arXiv:1707.08222v2.
 (14) Photonis, Planacon XP85012 Data Sheet (January 2013).

(15)
B. Adams, A. Elagin, H. Frisch, R. Obaid, E. Oberla, A. Vostrikov, R. Wagner,
J. Wang, M. Wetstein,
Timing characteristics of
large area picosecond photodetectors, Nuclear Instruments and Methods in
Physics Research Section A: Accelerators, Spectrometers, Detectors and
Associated Equipment 795 (2015) 1–11.
doi:10.1016/j.nima.2015.05.027.
URL http://dx.doi.org/10.1016/j.nima.2015.05.027 
(16)
S. A. Pozzi, E. Padovani, M. Marseguerra,
Mcnppolimi:
a montecarlo code for correlation measurements, Nuclear Instruments and
Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors
and Associated Equipment 513 (3) (2003) 550 – 558.
doi:10.1016/j.nima.2003.06.012.
URL http://www.sciencedirect.com/science/article/pii/S0168900203023027 
(17)
K. Ziock, J. Braverman, L. Fabris, M. Harrison, D. Hornback, J. Newby,
Event
localization in bulk scintillator crystals using coded apertures, Nuclear
Instruments and Methods in Physics Research Section A: Accelerators,
Spectrometers, Detectors and Associated Equipment 784 (Supplement C) (2015)
382 – 389, symposium on Radiation Measurements and Applications 2014 (SORMA
XV).
doi:https://doi.org/10.1016/j.nima.2015.01.013.
URL http://www.sciencedirect.com/science/article/pii/S016890021500039X  (18) C. Lane, S. M. Usman, J. Blackmon, C. Rasco, H. P. Mumm, D. Markoff, G. R. Jocher, R. Dorrill, M. Duvall, J. G. Learned, V. Li, J. Maricic, S. Matsuno, R. Milincic, S. Negrashov, M. Sakai, M. Rosen, G. Varner, P. Huber, M. L. Pitt, S. D. Rountree, R. B. Vogelaar, T. Wright, Z. Yokley, A new type of neutrino detector for sterile neutrino search at nuclear reactors and nuclear nonproliferation applications, ArxivarXiv:1501.06935v1.

(19)
S. Ritt, R. Dinapoli, U. Hartmann,
Application
of the DRS chip for fast waveform digitizing, Nuclear Instruments and
Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors
and Associated Equipment (2010) 1–3doi:10.1016/j.nima.2010.03.045.
URL http://linkinghub.elsevier.com/retrieve/pii/S0168900210006091 
(20)
M. Bitossi, R. Paoletti, D. Tescaro,
Ultrafast sampling and
data acquisition using the drs4 waveform digitizer, IEEE Transactions on
Nuclear Science 63 (4) (2016) 2309 – 2316.
doi:10.1109/TNS.2016.2578963.
URL http://dx.doi.org/10.1109/TNS.2016.2578963 
(21)
E. Oberla, J.F. Genat, H. Grabas, H. Frisch, K. Nishimura, G. Varner,
A 15GSa/s, 1.5GHz
bandwidth waveform digitizing ASIC, Nuclear Instruments and Methods in
Physics Research Section A: Accelerators, Spectrometers, Detectors and
Associated Equipment 735 (2014) 452–461.
arXiv:1309.4397v1,
doi:10.1016/j.nima.2013.09.042.
URL http://dx.doi.org/10.1016/j.nima.2013.09.042  (22) E. Technology, Ej232Q Data Sheet.

(23)
R. Barlow,
Extended
maximum likelihood, Nuclear Instruments and Methods in Physics Research
Section A: Accelerators, Spectrometers, Detectors and Associated Equipment
297 (3) (1990) 496 – 506.
doi:https://doi.org/10.1016/01689002(90)913348.
URL http://www.sciencedirect.com/science/article/pii/0168900290913348  (24) F. James, MINUIT Function Minimization and Error Analysis: Reference Manual Version 94.1, Tech. rep., CERN (1994).

(25)
L. A. Shepp, Y. Vardi,
Maximum likelihood
reconstruction for emission tomography., IEEE transactions on medical
imaging 1 (2) (1982) 113–22.
doi:10.1109/TMI.1982.4307558.
URL http://www.ncbi.nlm.nih.gov/pubmed/18238264 
(26)
S. Agostinelli, J. Allison, K. Amako, J. Apostolakis, H. Araujo, P. Arce,
M. Asai, D. Axen, S. Banerjee, G. Barrand, F. Behner, L. Bellagamba,
J. Boudreau, L. Broglia, A. Brunengo, H. Burkhardt, S. Chauvie, J. Chuma,
R. Chytracek, G. Cooperman, G. Cosmo, P. Degtyarenko,
A. Dell'Acqua, G. Depaola, D. Dietrich, R. Enami,
A. Feliciello, C. Ferguson, H. Fesefeldt, G. Folger, F. Foppiano, A. Forti,
S. Garelli, S. Giani, R. Giannitrapani, D. Gibin, J. G. Cadenas,
I. González, G. G. Abril, G. Greeniaus, W. Greiner, V. Grichine,
A. Grossheim, S. Guatelli, P. Gumplinger, R. Hamatsu, K. Hashimoto, H. Hasui,
A. Heikkinen, A. Howard, V. Ivanchenko, A. Johnson, F. Jones, J. Kallenbach,
N. Kanaya, M. Kawabata, Y. Kawabata, M. Kawaguti, S. Kelner, P. Kent,
A. Kimura, T. Kodama, R. Kokoulin, M. Kossov, H. Kurashige, E. Lamanna,
T. Lampén, V. Lara, V. Lefebure, F. Lei, M. Liendl, W. Lockman,
F. Longo, S. Magni, M. Maire, E. Medernach, K. Minamimoto, P. M. de Freitas,
Y. Morita, K. Murakami, M. Nagamatu, R. Nartallo, P. Nieminen, T. Nishimura,
K. Ohtsubo, M. Okamura, S. O'Neale, Y. Oohata, K. Paech,
J. Perl, A. Pfeiffer, M. Pia, F. Ranjard, A. Rybin, S. Sadilov, E. D. Salvo,
G. Santin, T. Sasaki, N. Savvas, Y. Sawada, S. Scherer, S. Sei, V. Sirotenko,
D. Smith, N. Starkov, H. Stoecker, J. Sulkimo, M. Takahata, S. Tanaka,
E. Tcherniaev, E. S. Tehrani, M. Tropeano, P. Truscott, H. Uno, L. Urban,
P. Urban, M. Verderi, A. Walkden, W. Wander, H. Weber, J. Wellisch,
T. Wenaus, D. Williams, D. Wright, T. Yamada, H. Yoshida, D. Zschiesche,
Geant4—a
simulation toolkit, Nuclear Instruments and Methods in Physics Research
Section A: Accelerators, Spectrometers, Detectors and Associated Equipment
506 (3) (2003) 250–303.
doi:10.1016/s01689002(03)013688.
URL http://dx.doi.org/10.1016/S01689002(03)013688 
(27)
S. A. Pozzi, J. A. Mullens, J. T. Mihalczo,
Analysis
of neutron and photon detection position for the calibration of plastic
(bc420) and liquid (bc501) scintillators, Nuclear Instruments and Methods
in Physics Research Section A: Accelerators, Spectrometers, Detectors and
Associated Equipment 524 (1) (2004) 92 – 101.
doi:10.1016/j.nima.2003.12.036.
URL http://www.sciencedirect.com/science/article/pii/S0168900204001342 
(28)
R. Batchelor, W. Gilboy, J. Parker, J. Towle,
The response of organic
scintillators to fast neutrons, Nuclear Instruments and Methods 13 (1961)
70–82.
doi:10.1016/0029554x(61)901719.
URL http://dx.doi.org/10.1016/0029554X(61)901719 
(29)
S. N. Kasarova, N. G. Sultanova, C. D. Ivanov, I. D. Nikolov,
Analysis
of the dispersion of optical plastic materials, Optical Materials 29 (11)
(2007) 1481 – 1490.
doi:10.1016/j.optmat.2006.07.010.
URL http://www.sciencedirect.com/science/article/pii/S0925346706002473 
(30)
V. A. Grigoryev, V. A. Kaplin, T. L. Karavicheva, A. B. Kurepin, E. F.
Maklyaev, Y. A. Melikyan, D. V. Serebryakov, W. H. Trzaska, E. M. Tykmanov,
Study of the planacon
xp85012 photomultiplier characteristics for its use in a cherenkov detector,
Journal of Physics: Conference Series 675 (4) (2016) 042016.
URL http://stacks.iop.org/17426596/675/i=4/a=042016