Next Article in Journal
Energy Supply System Modeling Tools Integrating Sustainable Livelihoods Approach—Contribution to Sustainable Development in Remote Communities: A Review
Next Article in Special Issue
Comparison of Deep Reinforcement Learning and PID Controllers for Automatic Cold Shutdown Operation
Previous Article in Journal
Sectoral Analysis of the Fundamental Criteria for the Evaluation of the Viability of Wave Energy Generation Facilities in Ports—Application of the Delphi Methodology
Previous Article in Special Issue
Operational Resilience of Nuclear-Renewable Integrated-Energy Microgrids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Momentum-Dependent Cosmic Ray Muon Computed Tomography Using a Fieldable Muon Spectrometer

by
Junghyun Bae
and
Stylianos Chatzidakis
*
School of Nuclear Engineering, Purdue University, West Lafayette, IN 47906, USA
*
Author to whom correspondence should be addressed.
Energies 2022, 15(7), 2666; https://doi.org/10.3390/en15072666
Submission received: 18 February 2022 / Revised: 15 March 2022 / Accepted: 2 April 2022 / Published: 5 April 2022
(This article belongs to the Special Issue Latest Advances in Nuclear Energy Systems)

Abstract

:
Cosmic ray muon tomography has been recently explored as a non-destructive technique for monitoring or imaging dense well-shielded objects, classically not achievable with traditional tomographic methods. As a recent example of technology transition from high-energy physics to real-world engineering applications, cosmic ray muon tomography has been used with various levels of success in nuclear nonproliferation. However, present muon detection systems have no momentum measurement capabilities and recently developed muon-based radiographic techniques rely only on muon tracking. This unavoidably reduces resolution and requires longer measurement times thus limiting the widespread use of cosmic ray muon tomography. Measurement of cosmic ray muon momenta has the potential to significantly improve the efficiency and resolution of cosmic ray muon tomography. In this paper, we propose and explore the use of momentum-dependent cosmic ray muon tomography using multi-layer gas Cherenkov radiators, a new concept for measuring muon momentum in the field. The muon momentum measurements are coupled with a momentum-dependent imaging algorithm (mPoCA) and image reconstructions are presented to demonstrate the benefits of measuring momentum in cosmic ray muon tomography.

1. Introduction

Monitoring large and dense objects, e.g., cargo containers or spent nuclear fuel dry casks, or large facilities, e.g., nuclear reactors, using cosmic ray muons could assess the stored content and provide information on whether it agrees with the declared content, maximize non-proliferation, and minimize potential security threats [1,2,3,4,5]. Conventional radiographic methods are limited, e.g., X-rays and gamma rays, and result in less than a few centimeters of penetration in lead [1]. Increasing the photon energy does not increase penetrating distance due to an increase in production of electron–positron pairs. Similarly, electrons cannot penetrate far into matter due to large radiation emission (Bremsstrahlung). For example, 1 GeV electrons can penetrate lead only within a few centimeters and commercially available portable electron accelerators can generate electrons with energies less than 10 MeV. Protons have been shown to penetrate several meters in dense materials and their potential has been demonstrated for large object imaging [2]. However, production of such subatomic particles requires large and immovable facilities, e.g., accelerators. Lastly, neutrinos can penetrate several kilometers in matter but their detection is impractical outside the laboratory, requiring kiloton size detectors.
A non-conventional class of ionizing particles, cosmic ray muons, presents certain advantages over the abovementioned charged and uncharged radiation. Cosmic ray muons are positively or negatively charged particles with ~200 times the mass of electrons, generated naturally in the upper layers of the atmosphere, and reaching the earth with a mean energy of 3-4 GeV at a rate of ~10,000 particles m−2 min−1 [6]. Muons can penetrate large and dense materials and their range increases with increasing energy. Figure 1 shows the range of muons as a function of atomic number and energy. Muons have the largest penetrating capabilities even in dense materials such as uranium and lead. In addition, muons are freely available and no radiological sources are required, resulting in a zero artificial radiation dose. Although the existence of cosmic muons and the theory which describes their transport and Multiple Coulomb Scattering (MCS) in matter is well understood, it is only recently that the applicability of muon scattering for imaging has been investigated [1,7]. Muon tomography is based on MCS and connects the expected muon scattering angle with the material atomic number and muon momentum. According to MCS, the distribution of muon scattering angles can be approximated using a Gaussian distribution and the characteristic parameters of the distribution can be determined from muon momentum and material properties.
Currently, cosmic ray muon tomography is solely based on scattering angle measurement without the capability for muon momentum measurement. Recent efforts at INFN have developed a large prototype based on MCS coupled with muon momentum reconstruction algorithms [8]. Although these efforts show promise, no portable detectors exist that can measure muon momentum in the field. This requires longer measurement times to compensate for the reduction in image resolution. To address this limitation, a Cherenkov-based muon spectrometer is currently under development for measuring muon momentum in the field [9,10]. Using pressurized gas Cherenkov radiators and depending on the pressure level, we can create a set of increasing muon momentum thresholds. Each radiator is triggered only when the incoming muon momentum exceeds the muon momentum threshold level. As a result, the signals from each radiator depend on the incoming muon momentum as not all Cherenkov radiators can generate a signal. By analyzing which radiator is triggered, we can then estimate the actual muon momentum. In addition, we extend the widely used reconstruction algorithm, Point of Closest Approach (PoCA) to include momentum into a new momentum dependent PoCA (mPoCA) algorithm. In the end, results comparing the two image reconstruction algorithms, PoCA and mPoCA, are presented to show the benefits of momentum measurement in cosmic ray muon applications.

2. Overview of Muon Tomography Detection Capabilities

Cosmic ray muon tomography is based on the principle of multiple Coulomb scattering (MCS). The estimated muon scattering angle distribution can be approximated using a Gaussian distribution with zero mean and standard deviation [11]:
f θ | 0 , σ θ 2 = 1 2 π σ θ exp 1 2 θ 2 σ θ 2
σ θ = 13.6   MeV β μ p μ c X X 0 1 + 0.088 log X X 0
where θ is the muon scattering angle, σθ is the width of Gaussian scattering angle distribution, X is the thickness of the scattering medium, and X0 is the radiation length. βμ is the ratio of the muon’s speed to the speed of light (c) and pμ is the muon momentum.
Early efforts have mainly focused on muon detectors for cargo scanning applications and proof of principle studies [12,13,14,15] and several detector types have been proposed for muon tomography. Existing detectors measure the locations and directions of incoming and outgoing muons and include drift tubes, scintillators, resistive plate chambers, and solid state detectors. Muon detectors are comprised of multiple parallel planes of position-sensitive chambers although cylindrical geometries have also been proposed [16,17,18,19]. A typical arrangement is shown in Figure 2. Cox et al. [20] summarized the requirements for muon detector development. Spatial resolution in the order of sub-mm, coincidence timing in the order of nanoseconds, and picosecond timing resolution would be required to differentiate between sub-GeV and GeV muon energies for future applications. They noted that efficient detectors with a large area and sensitivity over hundreds of GeV would be needed to optimize statistics and reduce the time taken to collect an adequate number of muons.
Muon events are recorded only when muons pass through all detector planes within a time window of a few nanoseconds. Only muon trajectories that fall within the superimposed grid boundaries are useful for reconstruction. Density maps can be generated with millimeter to centimeter spatial resolution and different materials can be discriminated depending on their atomic number using reconstruction algorithms. The discrimination capability and image quality depend on muon flux and are usually limited by incomplete information about muon path and momentum [21]. Minimizing the number of muons provides substantial time savings but can worsen the image quality, whereas a higher number of measured muons can improve image contrast and resolution but also increases measurement time. Recent efforts have shown that muons in the order of millions are required for imaging large-scale objects such as shipping containers [17,22,23] or dry casks [19,24,25]. This translates to measurement time in the order of hours or even days, depending on the detector characteristics, i.e., size and orientation. Table 1 provides a summary of existing muon detectors.
The PoCA image reconstruction algorithm has been widely used in cosmic ray muon tomography with applications in cargo scanning, spent fuel cask imaging, and reactor monitoring. The algorithm assumes that muons are scattered at a single point. This ignores multiple Coulomb scattering resulting in resolution limitations and image blur. The PoCA algorithm, due to its minimal computing requirements, is typically applied for reconstructing large objects with thousands of voxels, within a reasonable amount of time. PoCA assigns a point at the minimum distance between the incoming and outgoing trajectories. The intersection of the two trajectories is described mathematically as follows:
d i s t L 1 , L 2 = min d i s t P , Q
where p, Q are points on the incoming trajectory L1 and the outgoing trajectory L2, respectively. In 3D space, where trajectories may not intersect, the algorithm calculates the middle point of the shortest line between trajectories and the scattering angle, θ, is assigned to a voxel on a 3D grid. At the estimated point of scatter the following value is assigned:
s i = 1 2 θ x o u t θ x i n 2 + θ y o u t θ y i n 2
where θx and θy are the projected angles onto an x-y plane with respect to a positive x- and y-axis, respectively (Figure 2). As more muons pass through the object volume, each pixel (or voxel) is filled with scattering angles. The final value of each voxel is the average scattering angle divided by the voxel size L:
V o x e l   v a l u e = 1 N L i = 1 N s i
where N is the number of scattering angle values in a voxel. Dividing by the voxel size is needed to normalize the voxel value for geometries where the voxel size is not the same everywhere. Materials with a higher atomic number will generate higher voxel values as a result of increased scattering in this location. The computational complexity of the algorithm depends on the number of muons passing through the object volume and scales with O(N) [24].
Other reconstruction algorithms that have been explored for use in cosmic ray muon tomography belong to the class of filtered back-projection algorithms (FBP). When a large number of evenly spaced projections is not available, a necessary condition for FBP, or significant scattering takes place along the ray path, a solution can be provided using algebraic reconstruction techniques (ART). ART is computationally expensive but conceptually simpler than FBP. In this algebraic approach, a grid is superimposed on top of the object. Each grid pixel contains the unknown object function and a set of algebraic equations is set up for the determination of the unknown function. A necessary condition for ART is the correct estimation of the ray paths between the source and detector. Simulation and modeling in cosmic ray muon tomography is typically undertaken using GEANT4 (GEometry ANd Tracking) [30,31,32,33], a Monte Carlo code for transport of high energy particles in matter developed at CERN. As schematically shown in Figure 3, GEANT4 requires parameters related to object geometry, muon energy, and angular distributions.

3. Fieldable Muon Spectrometer and Momentum-Dependent PoCA

Recent simulations showed that a significant number of muons (>106 muons) would be needed to adequately image large objects, e.g., identify missing spent nuclear fuel in dry casks [19,24,25]. This was confirmed by a recent experimental effort by LANL [26] which resulted in ~200 h of measurement time in order to register only 1.62 × 105 muons, about 100 times lower than theoretically possible. The simulated results did not account for any additional measurement time needed to address noise and/or detector measurement uncertainties. Minimizing time measurement requirements would require maximizing the number of muons used for imaging to muons measured, i.e., imaging efficiency. This motivates the use of muon momentum as an additional feature to improve cosmic ray muon tomography beyond what is currently feasible.

3.1. Fieldable Muon Spectrometer

When the velocity of a muon surpasses the speed of light in an optically transparent medium, or a radiator, Cherenkov light is emitted (Figure 4). Therefore, there is a threshold muon momentum level for muons to induce Cherenkov radiation in a given radiator and it is derived as follows [34]:
p t h c = m μ c 2 n 2 1
where pth is the threshold muon momentum, n is the refractive index of the medium, and mμc2 is the muon rest mass energy. In gas radiators, the refractive index depends on pressure and temperature and it can be approximated by the Lorentz–Lorenz equation:
n 1 + 3 A m p R T
where Am is the molar refractivity, p is the gas pressure, T is the absolute temperature, and R is the universal gas constant. Therefore, the threshold muon momentum is given by:
p t h = m μ c R 3 A m T p
Equation (8) shows the relation between the Cherenkov threshold momentum for muons and gas pressure. We choose CO2 because it can be pressurized (up to 5.7 MPa) without condensation. Continuous threshold momentum levels with an interval of 1.0 GeV/c are designed by pressurizing CO2 gas radiators, pth = 1.0, 2.0, …, 5.0 GeV/c (Figure 4). We used a SiO2 solid radiator to achieve the lowest threshold momentum level (pth = 0.1 GeV/c) because it is not possible to obtain with pressurized CO2 at room temperature. As a result, a signal will be generated only when the momentum threshold is less than the actual muon momentum. This applies even in the case where a muon interacts with all radiators. The actual muon momentum can be estimated by measuring the Cherenkov signals generated or not generated in each radiator.

3.2. Momentum-Dependent PoCA

To include momentum in a PoCA algorithm, the measured muon momentum is additionally tagged in a PoCA point along with a scattering angle and this momentum-dependent PoCA imaging algorithm is named mPoCA. Therefore, each muon event includes five values: 3D Cartesian components of a PoCA point (Px, Py, and Pz), the reconstructed spatial muon scattering angle, θspace, and the measured muon momentum, pμ,
X i = P x ,   P y ,   P z , θ space ,   p μ i i = 1 ,   ,   N μ
where Nμ is the total number of muon events. All measured muon momenta are classified by a finite number of momentum levels depending on the number of radiators in the muon spectrometer. Therefore, Xi (Equation (9)) can be decomposed with respect to momentum groups:
X T = X 1 T , X 2 T ,   ,   X N p 1 T ,   X N p T
X j = V l m n 1 θ 1   V l m n 2 θ 2 V l m n L X θ L j
where LX is the length of Xj, Np is the number of muon momentum groups, and Vlmn is the PoCA voxel when l, m, and n are integers. Because the expected muon scattering angle distribution depends on muon momentum, each XjT (j = 1, 2, …, Np) will have a unique scattering angle distribution. This decomposition process was not possible without muon momentum measurement.

4. GEANT4 Simulations

The feasibility of the proposed fieldable muon spectrometer was evaluated using GEANT4. To generate muons from the actual measured muon spectrum, a “Muon Event Generator” is used based on a phenomenological model that captures the main characteristics of the experimentally measured spectrum coupled with a set of statistical algorithms [36,37]. The muons generated have zenith angles in the range 0–90° and energies in the range 1–100 GeV. The “Muon Event Generator” generates the necessary user-defined histograms for use with the GEANT4 G4GeneralParticleSource macro file. Hence, a stochastic muon transport simulation using GEANT4 was performed to evaluate the characteristics of the proposed muon spectrometer.
To measure muon momentum, we modeled a Cherenkov muon spectrometer in GEANT4 which consisted of one glass (1 cm), five CO2 gas radiators (10 cm), and thin photon-absorbing black-aluminum foils (0.1 cm) to minimize photon losses to the environment. To measure the produced optical photons, an array of silicon photomultiplier (SiPM) was placed in each radiator. The overall length was 51.7 cm and the surface area 0.2 × 0.2 m2. SiO2 with 1.45 refractive index was used for the glass radiator. Pressurized CO2 gases were used as gas radiators to provide sequentially increasing refractive indices and increasing momentum threshold levels. Some CO2 gas radiators needed to be depressurized to achieve pth = 4 and 5 GeV/c. Any optical photon events were terminated once they arrived either at the photon absorbers or outside of “world” volume during the simulations. All major physics are included in the GEANT4 reference physics list, QGSP_BERT [33].

5. Results

A proposed overview implementation in nuclear safeguards is shown in Figure 5. The spectrometer is installed between objects and muon trackers to measure momentum for all muon events. Because the performance of a muon spectrometer depends on the number of radiators, we assessed the performance by increasing the number of radiators to 10 and 100 radiators and extending the momentum range from 0.1 to 10.0 GeV/c. We performed four simulations with a different number of radiators using 104 muon samples in each case: (i) 102 radiators (fine measurement resolution, ±0.05 GeV/c) without noise, (ii) 102 radiators with noise (i.e., including scintillation and transition photon emission), (iii) 10 radiators (coarse resolution, ±0.5 GeV/c) without noise, and (iv) 10 radiators with noise. When the instrumental uncertainty is not considered, the cosmic ray muon spectrum was successfully reconstructed using 102 radiators (Figure 6). Although the increased number of radiators improved resolution, it negatively impacted the signal-to-noise ratio due to the decreased Cherenkov signals. When noise was added, the spectrum was slightly shifted to the right because noise (additional photon signal) increases the probability of false classification. We could improve the accuracy by suppressing the noise level because the momentum overestimation pattern was predictable. For the simulation results with 10 radiators, the cosmic ray muon spectrum was adequately reconstructed even if there were fewer radiators.
To statistically evaluate our muon spectrometer, we calculated the classification rate (CR), which is the probability that the system correctly identifies the actual muon momentum. The computed classification rates using 104 mono-energetic muons for various discrimination levels to remove noise, i.e., 0, 1, and 2, are shown in Figure 6. These linear discriminators were used to remove any noise signal from the total signal. When muon momentum was less than 1.0 GeV/c, the overall CR was less than 60% because the scintillation photons dominated in the glass radiator. When 1.0 < pth < 4.0 GeV/c, all gas radiators showed high CRs because the predictable noise level (mostly from the scintillation) was discriminated out from the total signal. This situation was further improved with a linear discriminator. However, when 4.0 < pth < 6.0 GeV/c, CRs began to decrease. This is attributed to low Cherenkov photon yields due to the depressurized CO2 (~a half atmosphere pressure) gas. When pth > 6.0 GeV/c, on the other hand, CR rebounded because Cherenkov photon yields gradually increase as muon momentum increases. However, it is not efficient to use any linear discriminator in high pth (>5.0 GeV/c).
To explore the benefits of muon momentum measurement in nuclear safeguards, e.g., special nuclear material (SNM) monitoring, three spherical special nuclear material samples, i.e., significant amounts (SQ) of high-enriched uranium (HEU), low-enriched uranium (LEU), and plutonium (Pu), were studied using Monte Carlo simulations where they were surrounded by various thicknesses of lead shielding. 5 × 103 muons were used to build scattering angle variance distributions [38,39]. In Figure 7, the mean scattering angle variance of distributions for HEU, LEU, and Pu are plotted as a function of thickness of lead shielding with 2σ errors when the muon momentum follows a cosmic ray muon spectrum with pμ = 0.2–100.0 GeV/c (left column) and pμ = 3 GeV/c (right column). Without measuring muon momentum, the scattering angle variance distributions for the three types of SNMs begin to overlap as the thickness of lead shielding increases. On the other hand, they are clearly distinguishable with muon momentum knowledge even when the thickness of lead shielding is 30 cm. The following simulation results show the reconstructed images of spherical LEU samples surrounded by 5 cm thick lead shielding (Figure 8). We used 105 muon samples in GEANT4 to simulate Monte Carlo muon transports and muon tomographic images were reconstructed using a (i) PoCA algorithm with a cosmic ray muon spectrum when pμ = 0.2–100.0 GeV/c, (ii) PoCA algorithm with mono-energetic muon, pμ = 3 GeV/c, and (iii) mPoCA algorithm where materials were classified in three groups, M1, M2, and M3, for SNMs, shielding materials, and others, respectively. In the mPoCA algorithm, the most probable material group was recorded instead of scattering angle values (Equation (5)) in each voxel. The results show that the imaging resolution is significantly improved and boundaries between SNM (LEU) and lead shielding (5 cm) are visually distinguished when the number of noise voxels in the air gap decreases.

6. Conclusions

Cosmic ray muons are relativistic particles with energies in the order of GeV that can play an important role in non-proliferation and nuclear material monitoring applications because they present unique advantages over existing methods. Research on muon detectors and imaging approaches for real-time monitoring and imaging of large objects using muon scattering, displacement, and transmission are constantly improving and expanding the cosmic ray muon applications in various fields. In the present work, we evaluated a Cherenkov muon spectrometer using multi-layer pressurized gas Cherenkov radiators. Six sequential threshold momentum levels were set using a solid radiator (SiO2) and pressurized CO2 gases, pth = 0.1, 1.0, 2.0, 3.0, 4.0, and 5.0 GeV/c. Then, we performed detailed GEANT4 simulations to evaluate the functionality of the proposed fieldable muon spectrometer. We were able to measure muon momentum with a resolution of ±0.5 GeV/c (a mean CR is 87%) within a momentum range of 0.1 to 10.0 GeV/c. We showed the potential of measuring muon momentum in nuclear nonproliferation, e.g., special nuclear material monitoring, by comparing muon scattering distributions with and without muon momentum measurement. Finally, we extended the widely used reconstruction algorithm, PoCA, to include momentum into a new momentum-dependent PoCA (mPoCA) algorithm. With the mPoCA, imaging resolution improved and noise decreased by a factor of 3 to 4.

Author Contributions

Conceptualization, S.C. and J.B.; methodology, S.C. and J.B.; software, S.C. and J.B.; validation, S.C. and J.B.; formal analysis, S.C. and J.B.; investigation, S.C. and J.B.; resources, S.C. and J.B.; data curation, S.C. and J.B.; writing—original draft preparation, S.C. and J.B.; writing—review and editing, S.C. and J.B.; visualization, S.C. and J.B.; supervision, S.C.; project administration, S.C.; funding acquisition, S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC was funded by the Purdue Research Foundation grant number [60000029].

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank Hany Abdel-Khalik at Purdue University for invitation and support in writing this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Borozdin, K.N.; Hogan, G.E.; Morris, C.; William, C.; Saunders, A.; Schultz, L.J. Radiographic imaging with cosmic-ray muons. Nature 2003, 422, 20–22. [Google Scholar] [CrossRef] [PubMed]
  2. Gustafsson, J. Tomography of Canisters for Spent Nuclear Fuel Using Cosmic-Ray Muons; Uppsala University: Uppsala, Sweden, 2005. [Google Scholar]
  3. Liu, Z.; Chatzidakis, S.; Scaglione, J.M.; Liao, C.; Yang, H.; Hayward, J.P. Muon Tracing and Image Reconstruction Algorithms for Cosmic Ray Muon Computed Tomography. IEEE Trans. Image Process. 2019, 28, 426–435. [Google Scholar] [CrossRef] [PubMed]
  4. Chatzidakis, S.; Choi, C.K.; Tsoukalas, L.H. Interaction of cosmic ray muons with spent nuclear fuel dry casks and determination of lower detection limit. Nucl. Inst. Methods Phys. Res. A 2016, 828, 37–45. [Google Scholar] [CrossRef] [Green Version]
  5. Chatzidakis, S. Classification and Imaging of Spent Nuclear Fuel Dry Casks Using Cosmic Ray Muons; Oak Ridge National Laboratory (ORNL): Oak Ridge, TN, USA, 2017; pp. 237–245. [Google Scholar]
  6. Hagiwara, K.; Hikasa, K.; Nakamura, K.; Tanabashi, M.; Aguilar-Benitez, M. Review of Particle Properties (Particle Data Group). Phys. Rev. D 2002, 66, 01001. [Google Scholar] [CrossRef] [Green Version]
  7. Jenneson, P.M. Large vessel imaging using cosmic-ray muons. Nucl. Inst. Methods Phys. Res. A 2004, 525, 346–351. [Google Scholar] [CrossRef]
  8. Checchia, P.; Benettoni, M.; Bettella, G.; Conti, E.; Cossutta, L.; Furlan, M.; Gonella, F.; Klinger, J.; Montecassiano, F.; Nebbia, G.; et al. INFN muon tomography demonstrator: Past and recent results with an eye to near-future activities. Philos. Trans. R. Soc. A 2019, 377, 20180065. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Bae, J.; Chatzidakis, S. A Cosmic Ray Muon Spectrometer Using Pressurized Gaseous Cherenkov Radiators. arXiv 2021, arXiv:2110.15061. [Google Scholar]
  10. Bae, J.; Chatzidakis, S. Fieldable Muon Momentum Measurement using Coupled Pressurized Gaseous Cherenkov Detectors. arXiv 2022, arXiv:2201.02591. [Google Scholar]
  11. Lynch, G.R.; Dahl, O.I. Approximations to multiple Coulomb scattering. Nucl. Inst. Methods Phys. Res. 1991, 58, 6–10. [Google Scholar] [CrossRef]
  12. Österlund, M.; Wallin, L.; Mermod, P.; Flodin, A.; Nilsson, L.; Gustafsson, J.; Hayashi, M.; Pomp, S.; Ohrn, A.; Donnard, J.; et al. Tomography of canisters for spent nuclear fuel. Proc. Sci. 2006, 25, 030. [Google Scholar] [CrossRef] [Green Version]
  13. Jonkmans, G.; Anghel, V.N.P.; Jewett, C.; Thompson, M. Nuclear waste imaging and spent fuel verification by muon tomography. Ann. Nucl. Energy 2013, 53, 267–273. [Google Scholar] [CrossRef] [Green Version]
  14. Clarkson, A.; Hamilton, D.; Hoek, M.; Ireland, D.; Johnstone, J.; Kaiser, R.; Keri, T.; Lumsden, S.; Mahon, D.; McKinnon, B.; et al. The design and performance of a scintillating-fibre tracker for the cosmic-ray muon tomography of legacy nuclear waste containers. Nucl. Inst. Methods Phys. Res. A 2014, 745, 138–149. [Google Scholar] [CrossRef] [Green Version]
  15. Thomay, C.; Velthuis, J.; Poffley, T.; Baesso, P.; Cussans, D.; Frazão, L. Passive 3D imaging of nuclear waste containers with Muon Scattering Tomography. J. Instrum. 2016, 11, P03008. [Google Scholar] [CrossRef] [Green Version]
  16. Anghel, V.; Armitage, J.; Baig, F.; Boniface, K.; Boudjemline, K.; Bueno, J.; Charles, E.; Drouin, P.L.; Erlandson, A.; Gallant, G.; et al. A plastic scintillator-based muon tomography system with an integrated muon spectrometer. Nucl. Inst. Methods Phys. Res. A 2015, 798, 12–23. [Google Scholar] [CrossRef] [Green Version]
  17. Schultz, L.J. Cosmic Ray Muon Radiography. Ph.D. Thesis, Portland State University, Portland, Oregon, 2003. [Google Scholar]
  18. Hoch, R. Advances in Cosmic Ray Muon Tomography Reconstruction Algorithms. Master’s Thesis, Florida Institute of Technology, Melbourne, FL, USA, 2009. [Google Scholar]
  19. Poulson, D.; Durham, J.M.; Guardincerri, E.; Morris, C.L.; Bacon, J.D.; Plaud-Ramos, K.; Morley, D.; Hecht, A.A. Cosmic ray muon computed tomography of spent nuclear fuel in dry storage casks. Nucl. Inst. Methods Phys. Res. A 2017, 842, 48–53. [Google Scholar] [CrossRef] [Green Version]
  20. Cox, L.; Adsley, P.; O’Malley, J.; Quillin, S.; Steer, C.; Clemett, C. Detector requirements for a cosmic ray muon scattering tomography system. In Proceedings of the 2008 IEEE Nuclear Science Symposium Conference Record, Dresden, Germany, 19–25 October 2008; pp. 706–710. [Google Scholar] [CrossRef]
  21. Chatzidakis, S.; Hausladen, P.A.; Croft, S.; Chapman, J.A.; Jarrell, J.J.; Scaglione, J.M.; Choi, C.K.; Tsoukalas, L.H. Exploring the Use of Muon Momentum for Detection of Nuclear Material Within Shielded Spent Nuclear Fuel Dry Casks. Am. Nucl. Soc. Annu. Meet. 2017, 116, 190–193. [Google Scholar]
  22. Tripathy, S.; Datta, J.; Jash, A.; Bouri, S.; Majumdar, N.; Mukhopadhyay, S.; Sarkar, S. Material Identification with Cosmic Ray Muons using RPC. J. Instrum. 2019, 14, C07007. [Google Scholar] [CrossRef] [Green Version]
  23. Riggi, S.; Antonuccio-Delogu, V.; Bandieramonte, M.; Becciani, U.; Costa, A. Muon tomography imaging algorithms for nuclear threat detection inside large volume containers with the Muon Portal detector. Nucl. Inst. Methods Phys. Res. A 2013, 728, 59–68. [Google Scholar] [CrossRef] [Green Version]
  24. Chatzidakis, S.; Choi, C.K.; Tsoukalas, L.H. Analysis of Spent Nuclear Fuel Imaging Using Multiple Coulomb Scattering of Cosmic Muons. IEEE Trans. Nucl. Sci. 2016, 63, 2866–2874. [Google Scholar] [CrossRef]
  25. Liu, Z.; Hayward, J.; Liao, C.; Yang, H. Detection of Missing Assemblies and Estimation of the Scattering Densities in a VSC-24 Dry Storage Cask with Cosmic-Ray-Muon-Based Computed Tomography. arXiv 2017, arXiv:1706.07072. [Google Scholar]
  26. Durham, J.M.; Guardincerri, E.; Morris, C.L.; Poulson, D.; Bacon, J.D.; Chichester, D.; Fabritius, J.; Fellows, S.; Plaud-Ramos, K.O.; Morley, D.J.; et al. Cosmic Ray Muon Imaging of Spent Nuclear Fuel in Dry Storage Casks. J. Nucl. Mater. Manag. 2016, 44. [Google Scholar]
  27. Green, J.A.; Alexander, C.; Asaki, T.; Bacon, J.; Blanpied, G.; Borozdin, K.; Canabal-Rey, A.; Cannon, M.; Chartrand, R.; Clark, D.; et al. Optimizing the tracking efficiency for cosmic ray muon tomography. In Proceedings of the 2006 IEEE Nuclear Science Symposium Conference Record, San Diego, CA, USA, 29 October–1 November 2006; Volume 1, pp. 285–288. [Google Scholar] [CrossRef]
  28. Pesente, S.; Vanini, S.; Benettoni, M.; Bonomi, G.; Calvini, P.; Checchia, P.; Conti, E.; Gonella, F.; Nebbia, G.; Squarcia, S.; et al. First results on material identification and imaging with a large-volume muon tomography prototype. Nucl. Inst. Methods Phys. Res. A 2009, 604, 738–746. [Google Scholar] [CrossRef]
  29. Åström, E.; Bonomi, G.; Calliari, I.; Calvini, P.; Checchia, P.; Donzella, A.; Faraci, E.; Forsberg, F.; Gonella, F.; Hu, X.; et al. Precision measurements of linear scattering density using muon tomography. J. Instrum. 2016, 11, P07010. [Google Scholar] [CrossRef] [Green Version]
  30. Agostinelli, S.; Allison, J.; Amako, K.; Apostolakis, J.; Araujo, H.; Arce, P.; Asai, M.; Axen, D.; Banerjee, S.; Barrand, G.; et al. GEANT4—A simulation toolkit. Nucl. Inst. Methods Phys. Res. Sect. A 2003, 506, 250–303. [Google Scholar] [CrossRef] [Green Version]
  31. Allison, J.; Amako, K.; Apostolakis, J.E.A.; Araujo, H.A.A.H.; Dubois, P.A.; Asai, M.A.A.M.; Barrand, G.A.B.G.; Capra, R.A.C.R.; Chauvie, S.A.C.S.; Chytracek, R.A.C.R.; et al. Geant4 developments and applications. IEEE Trans. Nucl. Sci. 2006, 53, 270–278. [Google Scholar] [CrossRef] [Green Version]
  32. Bogdanov, A.G.; Burkhardt, H.; Ivanchenko, V.N.; Kelner, S.R.; Kokoulin, R.P.; Maire, M.; Rybin, A.M.; Urban, L. Geant4 Simulation of High Energy Muon Interactions. In Proceedings of the IEEE Symposium Conference Record Nuclear Science, Rome, Italy, 16–22 October 2004; pp. 2043–2047. [Google Scholar]
  33. Geant4 Collaboration. Physics Reference Manual (Release 10.4); CERN: Geneva, Switzerland, 2017. [Google Scholar]
  34. Jelley, J.V. Čerenkov radiation: Its origin, properties and applications. Contemp. Phys. 1961, 3, 45–57. [Google Scholar] [CrossRef]
  35. Bae, J.; Chatzidakis, S. Fieldable muon spectrometer using multi-layer pressurized gas Cherenkov radiators and its applications. Sci. Rep. 2022, 12, 2559. [Google Scholar] [CrossRef]
  36. Chatzidakis, S.; Chrysikopoulou, S.; Tsoukalas, L.H. Developing a cosmic ray muon sampling capability for muon tomography and monitoring applications. Nucl. Inst. Methods Phys. Res. A 2015, 804, 33–42. [Google Scholar] [CrossRef] [Green Version]
  37. Smith, J.A.; Duller, N.M. Effects of pi meson decay-absorption phenomena on the high-energy mu meson zenithal variation near sea level. J. Geophys. Res. 1959, 64, 2297–2305. [Google Scholar] [CrossRef]
  38. Bae, J.; Chatzidakis, S. A New Semi-Empirical Model for Cosmic Ray Muon Flux Estimation. Prog. Theor. Exp. Phys. 2022, ptac016. [Google Scholar] [CrossRef]
  39. Bae, J.; Chatzidakis, S. The Effect of Cosmic Ray Muon Momentum Measurement for Monitoring Shielded Special Nuclear Materials. arXiv 2021, arXiv:2109.02470. [Google Scholar]
Figure 1. Range (mass density units) of various radiation types in various materials (left). Range in (cm) uranium as a function of energy (right) [5].
Figure 1. Range (mass density units) of various radiation types in various materials (left). Range in (cm) uranium as a function of energy (right) [5].
Energies 15 02666 g001
Figure 2. Typical arrangement of muon detectors in muon tomography applications.
Figure 2. Typical arrangement of muon detectors in muon tomography applications.
Energies 15 02666 g002
Figure 3. GEANT4 model construction using object geometry, muon energy, and angular distributions (top). In addition, post data analysis for decision-making and schematics of the PoCA algorithm principle are shown (center). The 2D reconstruction of a spent nuclear fuel dry cask using two different imaging techniques, muon scattering and transmission (bottom) (adapted from [5]).
Figure 3. GEANT4 model construction using object geometry, muon energy, and angular distributions (top). In addition, post data analysis for decision-making and schematics of the PoCA algorithm principle are shown (center). The 2D reconstruction of a spent nuclear fuel dry cask using two different imaging techniques, muon scattering and transmission (bottom) (adapted from [5]).
Energies 15 02666 g003
Figure 4. Overview of Cherenkov muon spectrometer using one solid (SiO2) and five pressurized CO2 radiators (top) [35]. The expected yields of Cherenkov and scintillation photons of each radiator (italic numbers represent radiator IDs) as a function of muon energy (lower left) and their wavelength distribution (lower right) are shown [35].
Figure 4. Overview of Cherenkov muon spectrometer using one solid (SiO2) and five pressurized CO2 radiators (top) [35]. The expected yields of Cherenkov and scintillation photons of each radiator (italic numbers represent radiator IDs) as a function of muon energy (lower left) and their wavelength distribution (lower right) are shown [35].
Energies 15 02666 g004
Figure 5. Cherenkov muon spectrometer implementation in nuclear nonproliferation. The spectrometer is installed between objects and muon trackers.
Figure 5. Cherenkov muon spectrometer implementation in nuclear nonproliferation. The spectrometer is installed between objects and muon trackers.
Energies 15 02666 g005
Figure 6. Cosmic ray muon spectra reconstructed using the 10 (upper left) and 100 (upper right) radiators without (dotted line) and with noise (solid line). 104 muon samples were simulated. The histogram represents the classified actual cosmic ray muon spectrum using 10 and 102 momentum levels, respectively. The lower figure shows classification rates for various muon momenta from 0.1 to 10.0 GeV/c using 104 muons. CO2 radiators with uniform radiator size and linear threshold momentum levels were used.
Figure 6. Cosmic ray muon spectra reconstructed using the 10 (upper left) and 100 (upper right) radiators without (dotted line) and with noise (solid line). 104 muon samples were simulated. The histogram represents the classified actual cosmic ray muon spectrum using 10 and 102 momentum levels, respectively. The lower figure shows classification rates for various muon momenta from 0.1 to 10.0 GeV/c using 104 muons. CO2 radiators with uniform radiator size and linear threshold momentum levels were used.
Energies 15 02666 g006
Figure 7. The mean muon scattering angle variance for significant quantities LEU, HEU, and Pu as a function of lead shielding thickness with error bars (2σ) when pμ follows a cosmic ray muon spectrum, 0.2–100 GeV/c (top left) and pμ = 3 GeV/c (top right). The muon scattering angle variance distribution when a lead shielding thickness, LPb = 30 cm (bottom row). The number of muons is 5 × 103.
Figure 7. The mean muon scattering angle variance for significant quantities LEU, HEU, and Pu as a function of lead shielding thickness with error bars (2σ) when pμ follows a cosmic ray muon spectrum, 0.2–100 GeV/c (top left) and pμ = 3 GeV/c (top right). The muon scattering angle variance distribution when a lead shielding thickness, LPb = 30 cm (bottom row). The number of muons is 5 × 103.
Energies 15 02666 g007
Figure 8. The layout of a spherical LEU sample surrounded by 5 cm thick lead shielding (top). Three reconstructed images (cross-sectional images) using a (i) PoCA algorithm with a cosmic ray muon spectrum when pμ = 0.2–100.0 GeV/c (left), (ii) PoCA algorithm with mono-energetic muon momentum, pμ = 3 GeV/c (center), and (iii) mPoCA algorithm with a cosmic ray muon spectrum when pμ = 0.2–100.0 GeV/c (right) are shown. The number of muons is 5×103.
Figure 8. The layout of a spherical LEU sample surrounded by 5 cm thick lead shielding (top). Three reconstructed images (cross-sectional images) using a (i) PoCA algorithm with a cosmic ray muon spectrum when pμ = 0.2–100.0 GeV/c (left), (ii) PoCA algorithm with mono-energetic muon momentum, pμ = 3 GeV/c (center), and (iii) mPoCA algorithm with a cosmic ray muon spectrum when pμ = 0.2–100.0 GeV/c (right) are shown. The number of muons is 5×103.
Energies 15 02666 g008
Table 1. Main characteristics of various muon tomography detector prototypes.
Table 1. Main characteristics of various muon tomography detector prototypes.
Detector TypeActive AreaSampling VolumeResolutionEfficiency
LANL prototype [17]Delay line readout multiple drift wire chambers60 × 60 cm260 × 60 × 27 cm3400 μm>96%
Mini Muon Tracker [26]Drift wire chambers1.2 × 1.2 m2100 × 100 × 60 cm3400 μm (2 mrad)>90%
Large Muon Tracker [27]Drift wire chambers2.0 × 2.0 m22.0 × 2.0 × 1.5 m3400 μm (2 mrad)>96%
Muon Portal INFN [28]Plastic scintillator strips3 × 6 m23 × 6 × 2.8 m33.5 mm (3.5 mrad)>90%
INFN [29]Drift wire chambers3 × 2.5 m23 × 2.5 × 1.6 m3200 μm>90%
CRIPT [16]Plastic scintillator bars2.0 × 2.0 m21.6 × 2.0 × 2.0 m33 mm (6 mrad)>99.5%
RPC [22]Resistive plate chambers0.5 × 0.5 m20.09 × 0.5 × 0.5 m3300–900 μm>95%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bae, J.; Chatzidakis, S. Momentum-Dependent Cosmic Ray Muon Computed Tomography Using a Fieldable Muon Spectrometer. Energies 2022, 15, 2666. https://doi.org/10.3390/en15072666

AMA Style

Bae J, Chatzidakis S. Momentum-Dependent Cosmic Ray Muon Computed Tomography Using a Fieldable Muon Spectrometer. Energies. 2022; 15(7):2666. https://doi.org/10.3390/en15072666

Chicago/Turabian Style

Bae, Junghyun, and Stylianos Chatzidakis. 2022. "Momentum-Dependent Cosmic Ray Muon Computed Tomography Using a Fieldable Muon Spectrometer" Energies 15, no. 7: 2666. https://doi.org/10.3390/en15072666

APA Style

Bae, J., & Chatzidakis, S. (2022). Momentum-Dependent Cosmic Ray Muon Computed Tomography Using a Fieldable Muon Spectrometer. Energies, 15(7), 2666. https://doi.org/10.3390/en15072666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop