Open-Source Optimization of Hybrid Monte-Carlo Methods for Fast Response Modeling of NaI(Tl) and HPGe Gamma Detectors

TL;DR

This paper presents an open-source MATLAB-based Hybrid Monte-Carlo tool for fast, accurate modeling of gamma detector responses, making nuclear data analysis more accessible and efficient.

Contribution

It introduces a user-friendly, open-source MATLAB code utilizing Hybrid Monte-Carlo methods for rapid gamma detector response simulation, validated against experimental data.

Findings

The code achieves response calculations comparable to industry standards.

Simulation times are significantly reduced, e.g., 0.26 seconds on 16 cores.

Open-source format enhances accessibility for academic and research use.

Abstract

Modeling the response of gamma detectors has long been a challenge within the nuclear community. Significant research has been conducted to digitally replicate instruments that can cost over $100,000 and are difficult to operate outside a laboratory setting. Subsequently, there have been multiple attempts to create codes that replicate the response of sodium-iodide and high purity germanium detectors for the purpose of deriving data related to gamma ray interaction with matter. While robust programs do exist, they are often subject to export controls and/or they are not intuitive to use. Through the use of the Hybrid Monte-Carlo methods, MATLAB can be used to produce a fast first-order response of various gamma ray detectors. The combination of a graphics user interface with a numerical based script allows for an open-source and intuitive code. When benchmarked with experimental data from Co-60, Cs-137, and Na-22 the code can numerically calculate a response comparable to experimental and industry standard response codes. Through this code, it is shown that a savings in computational requirements and the inclusion of an intuitive user experience does not heavily compromise data when compared to other standard codes or experimental results. When the application is installed on a computer with 16 cores, the average time to simulate the benchmarked isotopes is 0.26 seconds and 1.63 seconds on a four-core machine. The results indicate that simple gamma detectors can be modeled in an open-source format. The anticipation for the MATLAB application is to be a tool that can be easily accessible and provide datasets for use in an academic setting requiring the gamma ray detectors. Ultimately, providing evidence that Hybrid Monte-Carlo codes in an open-source format can benefit the nuclear community.