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

Potential of Depth-of-Interaction-Based Detection Time Correction in Cherenkov Emitter Crystals for TOF-PET

Cherenkov light can improve the timing resolution of Positron Emission Tomography (PET) radiation detectors, thanks to its prompt emission. Coincidence time resolutions (CTR) of ~30 ps were recently reported when using 3.2 mm-thick Cherenkov emitters. However, sufficient detection efficiency requires thicker crystals, causing the timing resolution to be degraded by the optical propagation inside the crystal. We report on depth-of-interaction (DOI) correction to mitigate the time-jitter due to the photon time spread in Cherenkov-based radiation detectors. We simulated the Cherenkov and scintillation light generation and propagation in 3 × 3 mm2 lead fluoride, lutetium oxyorthosilicate, bismuth germanate, thallium chloride, and thallium bromide. Crystal thicknesses varied from 9 to 18 mm with a 3-mm step. A DOI-based time correction showed a 2-to-2.5-fold reduction of the photon time spread across all materials and thicknesses. Results showed that highly refractive crystals, though producing more Cherenkov photons, were limited by an experimentally obtained high-cutoff wavelength and refractive index, restricting the propagation and extraction of Cherenkov photons mainly emitted at shorter wavelengths. Correcting the detection time using DOI information shows a high potential to mitigate the photon time spread. These simulations highlight the complexity of Cherenkov-based detectors and the competing factors in improving timing resolution.

Effect of crystal-photodetector interface extraction efficiency on Cerenkov photons' detection time

Using Cerenkov photons to improve detector timing resolution in time-of-flight positron emission tomography scanners is promising since they constitute most of the signal rising edge. The main challenge in using Cerenkov light is its low yield per photoelectric interaction, which requires optimizing their complex optical transport in the detector. Monte Carlo simulations unlock information unavailable through benchtop measurements and help better understand the Cerenkov photon behavior. Although the first Cerenkov photons are emitted forward, part of the early triggering signal is lost due to poor light extraction at the crystal-photodetector interface. In addition, the electron path in the crystal, that determines the Cerenkov photon direction, is tortuous due to multiple scattering, causing the Cerenkov photons emitted after a few scatters to no longer be forward-directed. In this context, the transit time spread in the crystal, highly dependent on the detector geometry, plays a crucial role in the photon detection time. In this work, we performed optical simulations in bismuth germanium oxide using 511 keV gamma with GATE to investigate the optical photons extraction when modifying the index of refraction at the crystal-photodetector interface and the crystal aspect ratio. The mean detection time of the first, second, and third detected optical and Cerenkov photon separately was studied as a function of the total number of Cerenkov detected per event. For each configuration, we calculated the expected mean detection time using the probability of detection. Thinner crystals led to lower expected detection times due to the reduced transit time in the crystal. Reducing the refractive index discontinuity at the crystal-photodetector interface decreased all configurations expected mean detection time values. We showed that it not only improves the optical photons (scintillation and Cerenkov) detection efficiency at the photodetector face but directly ameliorates the probability of detecting the fastest one, reducing the effect of thicker materials and of losing the first detected photon information, both crucial to reduce the detector timing resolution. Thanks to their prompt emission and directionality at emission, Cerenkov photons represent the first detected optical photon in most configurations but increasing their detection efficiency is crucial to detect the fastest one.

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.