## An Associated Particle Imaging System for Soil-Carbon Measurements

Mauricio Ayllon Unzueta a, b , Eoin Brodie , Craig Brown , Cristina Castanha , Charles Gary , c d c d Caitlin Hicks Pries , William Larsen , Bernhard Ludewigt , Andrew Rosenstrom , and Arun e a a a Persaud a

a

Accelerator Technology &amp; Applied Physics, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, US b Nuclear Engineering, University of California, Berkeley, CA 94709, US c Earth and Environmental Sciences, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, US d Adelphi Technology Inc., 2003 E. Bayshore Road, Redwood City, CA 94063, US e Department of Biological Sciences, Dartmouth College, Hanover, NH 03755, US

## ABSTRACT

We present first results from experimental data showing the capabilities of an Associated Particle Imaging system to measure carbon in soil and other elements. Specifically, we present results from a pre-mixed soil sample containing pure sand (SiO ) and 4 % carbon by 2 weight. Because the main isotopes of all those three elements emit characteristic high-energy gamma rays following inelastic neutron scattering, it is possible to measure their distribution with our instrument. A 3D resolution of several centimeters in all dimensions has been demonstrated.

detector, lanthanum bromide (LaBr ) and sodium io3 dine (NaI) gamma-ray detectors, and a fast digital data acquisition (DAQ) system, that allows for high-spatial resolution of several centimeters in all three dimensions. The basic API method is depicted in Figure 1.

Keywords: Associated Particle Imaging, 14 MeV neutrons, DT neutron generator, soil-carbon

## 1. INTRODUCTION

Associated Particle Imaging (API) has been used in the past for detection of illicit drugs, 1 explosives, 2 and special nuclear material (SNM) 3 among other applications. At Lawrence Berkeley National Laboratory (LBNL), we developed an API instrument for quantifying carbon distributions in soil. Non-destructive carbon-insoil measurement methods are important for understanding and quantifying soil-based carbon sequestration techniques that offer sequestration potential on a large scale since soil is the largest storage pool of terrestrial carbon. 4 Monitoring carbon in soil also supports improvements in soil health and crop yield. 5 Towards these goals, our API system employs state-ofthe-art components, including a custom neutron generator equipped with a position sensitive alpha-particle

Further author information: (Send correspondence to A. Persaud)

A. Persaud: E-mail apersaud@lbl.gov

Figure 1. Schematic illustration of API. A 14 1 MeV neu-. tron and an alpha particle are created at a known location. The neutron undergoes inelastic scattering off an isotope in the target, for example a 12 C nucleus, which emits a 4 4 MeV gamma ray. . The location of the interaction can be determined in 3D from the time difference between the detection of the alpha particle and the gamma ray and from the xy -measurement of the position of the alpha particle in the detector.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00950v1-figure-20592379-3cdf-4377-a40e-2281daa02262.jpg)

Neutrons are generated in deuterium(D) - tritium(T) fusion reactions (Equation 1) at beam energies of up to 100 keV. The two DT fusion reaction products, the alpha particle and the neutron, are emitted at 180 ◦ from each other in the center-of-mass reference frame. Therefore, by detecting the alpha particle in a positionsensitive detector, it is possible to calculate the direction of the associated neutron. The distance the neutron has traveled can be estimated from the recorded time difference between the detection of the alpha par- ticle and the associated neutron-induced gamma ray. Taken together one can reconstruct a full 3D distribution of the irradiated volume.

$$D + T \to \alpha \left ( 3. 5 \, \text{MeV} \right ) + n \, \left ( 1 4. 1 \, \text{MeV} \right ) \quad \left ( 1 \right ) \stackrel { \, \searrow \, \dots } { \text{Crv} t u } }$$

This technique offers an opportunity to nondestructively measure carbon distributions in heterogeneous soils, including changes with depth, and providing much larger representative samples than what is currently possible with standard core sampling techniques. Changes over time can be observed by repeatedly measuring the same location.

## 2. SYSTEM COMPONENTS

The main components of the API system include the neutron generator, the alpha detector, two gamma detectors (LaBr 3 and NaI), and the DAQ system.

## 2.1 Neutron generator

The neutron generator was designed and built by Adelphi Inc 6 based on their DT108-API generator. It is a sealed-type, compact neutron generator that consists of a microwave-driven ion source floated at a positive electrical potential, a high-voltage power supply with a maximum output of +100kV, a hydrogen getter that controls the gas pressure inside the chamber, and a titanium target explosion-bonded to a copper backing at ground potential. A sapphire vacuum window located 6 cm from the target enables light transport from the scintillator detecting the alpha particle to the photomultiplier tube (PMT). The neutron generator together with a schematic drawing is shown in Figure 2. The neutron generator combines a high output of 2 × 10 8 neutrons / s for high counting rates with a very small beam spot (2 mm diameter) on target for excellent spatial resolution.

Figure 2. Left: image of the DT108-API neutron generator tube and right: conceptual illustration of neutron generation and alpha detection.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00950v1-figure-fb091a9f-4927-4003-a736-3ab3e4f46201.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00950v1-figure-e2e3f178-d7c6-4ef7-8d1b-4865b9646ad4.jpg)

## 2.2 Alpha detector

The alpha detector consists of several components: an inorganic scintillator crystal, Yttrium Aluminum Perovskite (YAP), with dimension 50 mm × 50 mm × 1 mm (Crytor ) which is placed approximately 6 cm from the 7 beam target inside the neutron generator. The photons generated inside the YAP are transported through a 3 mm vacuum sapphire window (MPF Products ) to 8 the photocathode of a PMT (Hamamatsu H13700-03 ) 9 located outside the vacuum chamber. An aluminum coating on the beam-facing side of the YAP acts as a mirror and doubles the amount of photons arriving at the PMT. The coating also blocks back-scattered deuterium and tritium ions from hitting the scintillator. The PMT is currently read out with a four-corner readout board enabling accurate measurement of the position ( &lt; 1 mm on the alpha-detector, corresponding to &lt; 5 cm at 60 cm from the neutron source). For more details on the alpha detector please see our previous publication. 10 However, the current board design also limits the maximum allowable alpha detection rate to about 2 5 . × 10 5 alphas / s (corresponding to 5 × 10 6 neutrons / s taking the solid angle of the alpha detector into account) due to longer RC delays of the measured signals and the fact that each alpha particle creates a signal on all four channels. Both facts resulting in pile-up issues at higher rates. Future experiments will utilize a new board design to allow operation at higher rates.

## 2.3 Lanthanum bromide and sodium iodide gamma detectors

The gamma detectors, both inorganic scintillators, were chosen because of their fast scintillation light decay time and relatively good energy resolution. The 3 inch LaBr 3 detector (Saint Gobain B390S) excels at both of these requirements since its scintillation light has a 25 ns decay time, and its energy resolution at 4 4 MeV is . approximately 1 %. Additionally, we utilize a larger (5 inch) NaI detector (Alpha Spectra 20I20/5(9823)BN). However, LaBr 3 is significantly more expensive than the NaI scintillators. Scintillation light form NaI has a decay time of 120 ns and a resolution of 3 % at 4 4 MeV. . The NaI detector performance will be compared to that of the LaBr 3 detector, and it will be used to boost the count rate from deeper soil layers where we might be able to trade-off the timing precision for total measurement time.

## 2.4 Data acquisition system

For data acquisition and pre-processing, we use the PIXIE-16, 11 which is a 16-channel digital DAQ. It is capable of complex coincidence logic between multiple

Figure 3. API measurements of soil, graphite, and sand samples. Left: experimental setup; right: x y -distribution of detected events. The measurement also shows pickup of events from the wall and some neutron shielding on the left side.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00950v1-figure-ae540ce5-0407-4079-a1fe-6ca492db5666.jpg)

channels, implements an energy filter with on-the-fly pileup detection, as well as a constant fraction discrimination (CFD) algorithm for time-of-arrival detection. It consists of 16 analog-digital converters (500 MHz) that feed into four field-programmable gate arrays (FPGAs) and a single digital signal processor (DSP) for the online calculation. In list mode, events consisting of timestamped energies are buffered locally and then read out to a computer and written to disk.

The time of arrival is obtained with a CFD algorithm, which produces an output signal a fixed time after the leading edge of the pulse has reached a constant fraction of the pulse amplitude. This allows for the timing signal to be amplitude independent over a wide energy range [12, p. 662]. Measurements (in a forthcoming publication) using a 1 inch thick dry sand sample where used to obtain the full width half maximum (FWHM) of the time difference between the alpha detection (YAP scintillator) and each of the gamma signals using a customized firmware for the PIXIE-16 system with optimized parameters for the CFD algorithm. Time resolutions of 1 3 ns for the LaBr . 3 and 1 6 ns for . the NaI detector were measured this way. The measurement still includes contributions to the FWHM from the horizontal extend of the sample as well as a component due to the measurement in the alpha detector (estimated to be on the same order as the LaBr 3 detector). Further optimization of the parameters used in the CFD algorithm could potentially improve the time resolution to approximately 1 ns for the YAP - LaBr 3 detector combination. The currently achieved FWHM of 1 3 ns corresponds to a depth resolution of 7 cm. .

## 3. EXPERIMENTAL SETUP AND RESULTS

Sand (SiO ) 2 was mixed with 4 % carbon by weight to simulate real soil and placed in several aluminum boxes (each 3 inch × 4 inch × 5 inch). Well graded

Figure 4. API measurements of the same setup as shown in Figure 3, showing the y z -projection. The neutron generator is located at a z -location of 0 cm. The API cone is clearly visible from some neutron shielding at the sides. Due to the contrast, the graphite and sand boxes are not clearly visible in this view.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00950v1-figure-95e2c9df-0111-4e24-875b-3a0322099b79.jpg)

sand was sieved to remove fine pebbles then mixed with agricultural-grade worm castings, both oven-dried at 175 ◦ C to remove excess moisture. Worm castings have a relatively high carbon content (38 3 % C by weight as . measured by loss on ignition method) and were thoroughly combined with sand to produce homogeneous, artificial soil mixtures of 4 % carbon by weight. Boxes containing this mixture where stacked in a 2 × 3 × 3 array centered below the neutron source with the top surface of the boxes located 70 cm from the neutron source, and irradiated at a neutron generator output of 5 × 10 6 neutrons / s (acceleration voltage of 50 kV). Additionally, we placed a graphite brick of dimensions 22 cm × 13 cm × 6 5 cm . and another aluminum box filled with pure sand right next to the test soil to act as standards, as shown in Figure 3. Data was acquired in coincidence mode for 9 hours (future improvement will greatly reduce this measurement time, see discussion in Section 4). Figure 3 shows the experimental setup and the reconstructed x y -image and Figure 4 shows the y z -projection of the same measurement. Spectra of the different regions in space according to the three different samples are shown in Figure 5.

Note that the API technique allows for the measurement of gamma-ray spectra from selected regions in space, greatly reducing counts from unwanted events that would show up in other measurement techniques. This has proven to be useful in other areas such as the one reported by Litvak et al. 13 where spectroscopic analysis is rather difficult due to high background levels coming from structural materials. The gamma spectra in Figure 5 shows several peaks which are the result of inelastic neutron scattering (INS) on the three main elements (Si, C, O), several transitions between excited states of each nuclei, and interactions of gamma rays within the detector crystal. For instance, the graphite brick spectrum shows mostly the result of the transition

Figure 5. LaBr 3 gamma spectra taken with of soil, sand, and graphite samples. See Figure 3 for details on the experimental setup. Gamma rays with energies &lt; 3 MeV are suppressed due the CFD threshold being set too high in the PIXIE-16 DAQ. Also note that the resolution at 4 4 MeV is . about 3 %, which is higher than previously measured resolution for the same detector. This is most likely due to gain drift during the long measurement.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00950v1-figure-46fbc8e1-ab12-454f-810d-9a4c0499bc6d.jpg)

of the 4 4 MeV excited state in . 12 C to its ground state. The peak at 3 9 MeV is the single escape peak due to . pair production within the crystal. Also notice the interfering gamma line from 28 Si at 4 5 MeV in the soil . and sand spectra that needs to be taken into account for an accurate measurement of carbon distribution in soil. The unique ability to obtain gamma spectra from specific regions of space isolated from structural materials and unwanted background is clearly exemplified in the gamma spectrum of the graphite brick above. Table 1 shows the main gamma peaks of interest identified in our soil spectrum.

Table 1. Main inelastic gamma lines associated with the most common elements in soil. One important note is that the 28 Si peak is close to the 12 C peak which it hard to resolve, especially for the NaI detector.

| Element   |   Energy [MeV] |
|-----------|----------------|
| 28 Si     |          1.78  |
| 27 Al     |          2.21  |
| 16 O      |          2.74  |
| 28 Si     |          2.84  |
| 16 O      |          3.68  |
| 16 O      |          3.85  |
| 12 C      |          4.44  |
| 28 Si     |          4.497 |
| 16 O      |          6.13  |

## 4. CONCLUSION AND OUTLOOK

First experimental results and characterization of capabilities where obtained for an API system developed at LBNL. We showed that by performing 3D reconstruction of the location of the neutron scattering centers using a position sensitive alpha detector together with time-of-flight analysis, we can achieve a &lt; 5 cm lateral resolution at 60 cm from the neutron source and a depth resolution of &lt; 7 cm. Additionally, we demonstrated the ability to acquire energy spectra from specific regions of interest thereby greatly reducing the unwanted events resulting from gammas produced outside the volume of interest. All three detectors show the expected energy and time resolution.

We are currently developing a new alpha detector readout board that is capable of handling higher rates (approximately × 5). Further planned changes to the readout board will allow us to operate at the maximum neutron generator output of 2 × 10 8 neutrons / s. These changes will allow us to cut the measurement time for the shown measurement down to less than 15 minutes. Additional detectors could be employed to reduce the measurement time even further. We are also in the process of developing the analysis software to extract accurate carbon density information from the gamma spectra and neutron flux measurements. Finally, we designed together with Adelphi Inc. a new neutron generator with several improvements that will have an impact on the overall resolution of the system, data quality, and portability.

## ACKNOWLEDGMENTS

The authors also would like to thank Takeshi Katayanagi for mechanical support. The information, data, or work presented herein was funded by the Advanced Research Projects Agency-Energy (ARPAE), U.S. Department of Energy, under Contract No. DEAC02-05CH11231.

## REFERENCES

- [1] C. L. Fontana, A. Carnera, M. Lunardon, F. Pino, C. Sada, F. Soramel, L. Stevanato, G. Nebbia, C. Carasco, B. Perot, A. Sardet, G. Sannie, A. Iovene, C. Tintori, K. Grodzicki, M. Moszyski, P. Sibczyski, L. Swiderski, and S. Moretto, 'Detection system of the first rapidly relocatable tagged neutron inspection system (rrtnis), developed in the framework of the european h2020 cbord project,' Physics Procedia 90 , pp. 279 - 284, 2017. Conference on the Application of Accelerators in Research and Industry, CAARI 2016, 30 October 4 November 2016, Ft. Worth, TX, USA.
- [2] C. Carasco, B. Perot, S. Bernard, A. Mariani, J.-L. Szabo, G. Sannie, T. Roll, V. Valkovic, D. Sudac, G. Viesti, M. Lunardon, C. Bottosso, D. Fabris, G. Nebbia, S. Pesente, S. Moretto,
- A. Zenoni, A. Donzella, M. Moszynski, M. Gierlik, T. Batsch, D. Wolski, W. Klamra, P. L. Tourneur, M. Lhuissier, A. Colonna, C. Tintori, P. Peerani, V. Sequeira, and M. Salvato, 'In-field tests of the euritrack tagged neutron inspection system,' Nucl. Instrum. Methods Phys. Res. A 588 (3), pp. 397 405, 2008.
- [3] Q. Ji, B. Ludewigt, J. Wallig, W. Waldron, and J. Tinsley, 'Development of a time-tagged neutron source for snm detection,' Physics Procedia 66 , pp. 105 - 110, 2015. The 23rd International Conference on the Application of Accelerators in Research and Industry - CAARI 2014.
- [4] J. Sanderman, T. Hengl, and G. J. Fiske, 'Soil carbon debt of 12,000 years of human land use,' Proc. Natl. Acad. Sci. U. S. A. 114 , pp. 95759580, Sept. 2017.
- [5] J. Huang, B. Minasny, A. B. McBratney, J. Padarian, and J. Triantafilis, 'The location- and scalespecific correlation between temperature and soil carbon sequestration across the globe,' Sci. Total Environ. 615 , pp. 540-548, Feb. 2018.
- [6] Adelphi Technology Inc. http://adelphitech.com/, 2018.
- [7] CRYTUR, spol. s r.o. https://crytur.cz, 2018.
- [8] MPF Products. https://mpfpi.com/, 2017.
- [9] Hamamatsu Photonics K.K. https://www.hamamatsu.com, 2018.
- [10] M. A. Unzueta, W. Mixter, Z. Croft, J. Joseph, B. Ludewigt, and A. Persaud, 'Position sensitive alpha detector for an associate particle imaging system,' arXiv.org , p. 1811.08591, Nov. 2018.
- [11] XIA, Pixie-16 User Manual , Aug. 2018.
- [12] G. Knoll, Radiation Detection and Measurement , Wiley, 3rd ed., 2010.
- [13] M. Litvak, Y. Barmakov, S. Belichenko, R. Bestaev, E. Bogolubov, A. Gavrychenkov, A. Kozyrev, I. Mitrofanov, A. Nosov, A. Sanin, V. Shvetsov, D. Yurkov, and V. Zverev, 'Associated particle imaging instrumentation for future planetary surface missions,' Nucl. Instrum. Methods Phys. Res. A 922 , pp. 19 - 27, 2019.
Mauricio Ayllon Unzueta Eoin Brodie Craig Brown Cristina Castanha Charles Gary Caitlin Hicks Pries William Larsen Bernhard Ludewigt Andrew Rosenstrom Arun Persaud
An Associated Particle Imaging System for Soil-Carbon Measurements
We present first results from experimental data showing the capabilities of
an Associated Particle Imaging system to measure carbon in soil and other
elements. Specifically, we present results from a pre-mixed soil sample
containing pure sand (SiO$_2$) and 4% carbon by weight. Because the main
isotopes of all those three elements emit characteristic high-energy gamma rays
following inelastic neutron scattering, it is possible to measure their
distribution with our instrument. A 3D resolution of several centimeters in all
dimensions has been demonstrated.