The Accuracy of Cerenkov Photons Simulation in Geant4/Gate Depends on the Parameterization of Primary Electron Propagation

Energetic electrons traveling in a dispersive medium can produce Cerenkov radiation. Cerenkov photons' prompt emission, combined with their predominantly forward emission direction with respect to the parent electron, makes them extremely promising to improve radiation detector timing resolution. Triggering gamma detections based on Cerenkov photons to achieve superior timing resolution is challenging due to the low number of photons produced per interaction. Monte Carlo simulations are fundamental to understanding their behavior and optimizing their pathway to detection. Therefore, accurately modeling the electron propagation and Cerenkov photons emission is crucial for reliable simulation results. In this work, we investigated the physics characteristics of the primary electrons (velocity, energy) and those of all emitted Cerenkov photons (spatial and timing distributions) generated by 511 keV photoelectric interactions in a bismuth germanate crystal using simulations with Geant4/GATE. Geant4 uses a stepwise particle tracking approach, and users can limit the electron velocity change per step. Without limiting it (default Geant4 settings), an electron mean step length of ~250 μm was obtained, providing only macroscopic modeling of electron transport, with all Cerenkov photons emitted in the forward direction with respect to the incident gamma direction. Limiting the electron velocity change per step reduced the electron mean step length (~0.200 μm), leading to a microscopic approach to its transport which more accurately modeled the electron physical properties in BGO at 511 keV. The electron and Cerenkov photons rapidly lost directionality, affecting Cerenkov photons' transport and, ultimately, their detection. Results suggested that a deep understanding of low energy physics is crucial to perform accurate optical Monte Carlo simulations and ultimately use them in TOF PET detectors.

Integration of polarization in the LUTDavis model for optical Monte Carlo simulation in radiation detectors

OBJECTIVE

Cerenkov photons have distinctive features from scintillation photons. Among them is their polarization: their electric field is always perpendicular to the direction of propagation of light and parallel to the plane of incidence. Scintillation photons are instead considered unpolarized.

APPROACH

This study aims at understanding and optimizing the reflectance of polarized Cerenkov photons for optical Monte Carlo simulation of scintillation detectors with Geant4/GATE. First, the Cerenkov emission spectrum and polarization were implemented in the previously developed look-up-table Davis model of crystal reflectance. Next, we modified Geant4/GATE source code to account for scintillation and Cerenkov photons LUTs simultaneously. Then, we performed optical Monte Carlo simulations in BGO using GATE to show the effect of Cerenkov features on the photons' momentum at the photodetector face, using two surface finishes, with and without reflector.

MAIN RESULTS

In this work, we describe the new features added to the algorithm and GATE. We showed that Cerenkov characteristics affect their probability to be reflected/refracted and thus their travel path within a material.

SIGNIFICANCE

We showed the importance of accounting for accurate Cerenkov photons reflectance while performing advanced optical Monte Carlo simulations.