Objective. The GPU-based Ultra-fast Monte Carlo positron emission tomography simulator (UMC-PET) incorporates the physics of the emission, transport and detection of radiation in PET scanners. It includes positron range, non-colinearity, scatter and attenuation, as well as detector response. The objective of this work is to present and validate UMC-PET as a a multi-purpose, accurate, fast and flexible PET simulator. Approach. We compared UMC-PET against PeneloPET, a well-validated MC PET simulator, both in preclinical and clinical scenarios. Different phantoms for scatter fraction (SF) assessment following NEMA protocols were simulated in a 6R-SuperArgus and a Biograph mMR scanner, comparing energy histograms, NEMA SF, and sensitivity for different energy windows. A comparison with real data reported in the literature on the Biograph scanner is also shown. Main results. NEMA SF and sensitivity estimated by UMC-PET where within few percent of PeneloPET predictions. The discrepancies can be attributed to small differences in the physics modeling. Running in a 11 GB GeForce RTX 2080 Ti GPU, UMC-PET is ∼1500 to ∼2000 times faster than PeneloPET executing in a single core Intel(R) Xeon(R) CPU W-2155 @ 3.30 GHz. Significance. UMC-PET employs a voxelized scheme for the scanner, patient adjacent objects (such as shieldings or the patient bed), and the activity distribution. This makes UMC-PET extremely flexible. Its high simulation speed allows applications such as MC scatter correction, faster SRM estimation for complex scanners, or even MC iterative image reconstruction.

Despite its success in many biomedical applications, Positron Emission Tomography (PET) has the drawback of typically having lower spatial resolution and higher noise respect to other medical imaging techniques. The best achievable spatial resolution in PET scanners is limited by factors such as the positron range, non-collinearity and the size of the detector crystals. In this work, we present a novel method that uses series of image reconstructions (super-iterations) to go beyond the expected resolution-noise limits for a given PET acquisition. The image quality improvement is achieved using the projections of the previous image reconstruction to redistribute the measured counts of each line-of-response (LOR) into several subLORs, from which a new activity distribution with better quality is reconstructed. The method was evaluated with data from the preclinical scanner 4R-SuperArgus PET/CT, using the NEMA NU4-2008 image quality phantom, a cold Derenzo phantom, and an in-vivo FDG cardiac study on a rat. Resolution and recovery coefficient (RC) improvement of ~10% was achieved while keeping the same noise level. Qualitative results from the in-vivo study also confirm this improvement in image quality. The proposed method is able to achieve significantly better images at the expense of a modest increase of the computational time, and it could be also applied to other modalities, such as SPECT and Compton Cameras.

Real-time positron emission tomography (PET) may provide information from first-shot images, enable PET-guided biopsies, and allow awake animal studies. Fully-3D iterative reconstructions yield the best images in PET, but they are too slow for real-time imaging. Analytical methods such as Fourier back projection (FBP) are very fast, but yield images of poor quality with artifacts due to noise or data incompleteness. In this work, an image reconstruction based on the pseudoinverse of the system response matrix (SRM) is presented. w. To implement the pseudoinverse method, the reconstruction problem is separated into two stages. First, the axial part of the SRM is pseudo-inverted (PINV) to rebin the 3D data into 2D datasets. Then, the resulting 2D slices can be reconstructed with analytical methods or by applying the pseudoinverse algorithm again. The proposed two-step PINV reconstruction yielded good-quality images at a rate of several frames per second, compatible with real time applications. Furthermore, extremely fast direct PINV reconstruction of projections of the 3D image collapsed along specific directions can be implemented.

Among other factors such as random, attenuation and scatter corrections, uniform spatial resolution is key to performing accurate quantitative studies in Positron emission tomography (PET). Particularly in preclinical PET studies involving simultaneous acquisition of multiple animals, the degradation of image resolution due to the depth of interaction (DOI) effect far from the center of the Field of View (FOV) becomes a significant concern. In this work, we incorporated a spatially-variant resolution model into a real time iterative reconstruction code to obtain accurate images of multi-animal acquisition. We estimated the spatially variant point spread function (SV-PSF) across the FOV using measurements and Monte Carlo (MC) simulations. The SV-PSF obtained was implemented in a GPU-based Ordered subset expectation maximization (OSEM) reconstruction code, which includes scatter, attenuation and random corrections. The method was evaluated with acquisitions from two preclinical PET/CT scanners of the SEDECAL Argus family: a Derenzo phantom placed 2 cm off center in the 4R-SuperArgus, and a multi-animal study with 4 mice in the 6R-SuperArgus. The SV-PSF reconstructions showed uniform spatial resolution without significant increase in reconstruction time, with superior image quality compared to the uniform PSF model.

Introduction: The positron range (PR) effect is a significant factor limiting spatial resolution in positron emission tomography (PET), particularly for high-resolution systems and non-standard isotopes. Objective: This study introduces a novel analytical model to accurately and rapidly describe positron range distributions (PRd) for various PET radioisotopes. Approach: The proposed model explicitly incorporates the Coulomb repulsion effect, the multi-branch nature of certain β⁺ emitters, and the scaling of PR with electronic density. To minimise bias, we used a histogram-free statistical method to derive the cumulative PRd from Monte Carlo (MC) simulated annihilation datasets, avoiding arbitrary histogram binning. A comparative analysis of PR estimates was conducted across three major MC radiation transport algorithm packages: PENELOPE (via PenEasy/PeneloPET), GEANT4 (via GATE), and EGS5 (via PHITS), revealing notable discrepancies between codes, versions, and input configurations, especially at short distances from the source. Main results: The new analytical model demonstrated an excellent reproduction of the simulated data for isotopes including C-11, N-13, O-15, F-18, Cu-64, Ga-68, Rb-82, and I-124, achieving in general coefficients of determination (R²) greater than 0.995 and mean absolute percentage errors (MAPE) < 20%. Compared to previous methods, our model provides a more accurate description of PRd at short and long distances and offers improved R² values. Significance: This work provides a robust framework for generating accurate annihilation point spread function (aPSF) kernels, facilitating improved PR correction in quantitative Nuclear Medical imaging and supporting research with diverse radioisotopes.

Preclinical PET animal studies require immobilization of the animal, typically accomplished through the administration of anesthesia, which may affect the radiotracer biodistribution. The use of ${}^{18}$F point sources attached to the rat head is one of the most promising methods for motion compensation in awake rat PET studies. However, the presence of radioactive markers may degrade image quality. In this study, we aimed to investigate the most favorable conditions for preclinical PET studies using awake rats with attached point sources. Firstly, we investigate the optimal activity conditions for the markers and rat-injected tracer using Monte Carlo simulations to determine the parameters of maximum detectability without compromising image quality. Additionally, we scrutinize the impact of delayed window correction for random events on marker detectability and overall image quality within these studies. Secondly, we present a method designed to mitigate the influence of rapid rat movements, which resulted in a medium loss of events of around 30%, primarily observed during the initial phase of the data acquisition. We validated our study with PET acquisitions from an awake rat within the acceptable conditions of activity and motion compensation parameters. This acquisition revealed an 8% reduction in resolution compared to a sedated animal, along with a 6% decrease in signal-to-noise ratio (SNR). These outcomes affirm the viability of our method for conducting awake preclinical brain studies.