Positronium-based imaging requires realistic modelling of positronium (Ps) decay in matter. We introduce a modular Ps decay model implemented in GATE 9.4 and GATE 10, enabling the definition of an arbitrary number of decay channels characterised by lifetime, branching fraction, annihilation multiplicity (2g/3g), and optional prompt photon emission. The model is validated through analytical and numerical benchmarks, including lifetime distributions, branching fraction consistency, photon kinematics, and prompt photon emission. Its practical applicability is demonstrated using simulations of mixed annihilation scenarios and the NEMA IEC phantom with a large field-of-view PET system. The proposed model accurately reproduces input lifetime distributions as weighted sums of exponential components and correctly samples decay channel fractions. Simulated two- and three-photon annihilation kinematics are consistent with theoretical expectations. Complex mixtures of decay channels, including varying 3g-to-2g ratios and multi-component ortho-positronium lifetimes, are correctly modelled, with observable signatures reflected in both temporal and energy distributions. Phantom simulations demonstrate the capability to generate realistic positronium-sensitive datasets. This work provides the first general-purpose, multi-channel positronium decay model integrated into GATE, enabling realistic simulations of positronium behaviour in complex media. The model supports the development and optimisation of positronium-based imaging techniques, including PLI and multi-photon PET, and applies to medical imaging, industrial tomography, and fundamental physics studies. Its public availability and compatibility with standard GATE workflows make it a valuable tool for the broader research community.

Objective. Positronium-based imaging, including positronium lifetime imaging (PLI) and multi-photon positron emission tomography (PET), requires realistic modelling of positronium (Ps) decay in matter. Existing Monte Carlo (MC) implementations in Geant4 Application for Tomographic Emission (GATE) are limited to simplified single- or two-channel descriptions, preventing accurate representation of the multi-component lifetime structure observed in biological and material systems. The objective of this work is to develop and validate a flexible MC Ps decay source model supporting multiple annihilation channels within the GATE framework.Approach. We introduce a modular Ps decay model implemented in GATE 9.4 and GATE 10, enabling the definition of an arbitrary number of decay channels characterised by lifetime, branching fraction, annihilation multiplicity (/), and optional prompt photon emission. The model is validated through analytical and numerical benchmarks, including lifetime distributions, branching fraction consistency, photon kinematics, and prompt photon emission. Its practical applicability is demonstrated using simulations of mixed annihilation scenarios and the NEMA IEC phantom with a large field-of-view PET system.Results. The proposed model accurately reproduces input lifetime distributions as weighted sums of exponential components and correctly samples decay channel fractions. Simulated two- and three-photon annihilation kinematics are consistent with theoretical expectations. Complex mixtures of decay channels, including varying-to-ratios and multi-component ortho-positronium (o-Ps) lifetimes, are correctly modelled, with observable signatures reflected in both temporal and energy distributions. Phantom simulations demonstrate the capability to generate realistic positronium-sensitive datasets.Significance. This work provides the first general-purpose, multi-channel positronium decay model integrated into GATE, enabling realistic simulations of positronium behaviour in complex media. The model supports the development and optimisation of positronium-based imaging techniques, including PLI and multi-photon PET, and applies to medical imaging, industrial tomography, and fundamental physics studies. Its public availability and compatibility with standard GATE workflows make it a valuable tool for the broader research community.

Total-Body PET is one of the most promising medical diagnostics modalities. The high sensitivity provided by Total-Body technology can be advantageous for novel tomography methods like positronium imaging. Several efforts are ongoing to lower the price of the TB-PET systems. Among the alternatives, the Jagiellonian PET (J-PET) technology, based on plastic scintillator strips, offers a low-cost alternative. The work aimed to compare five Total-Body J-PET geometries as a possible next generation J-PET scanner design. We present comparative studies of performance characteristics of the cost-effective Total-Body PET scanners using J-PET technology. We investigated in silico five Total-Body scanner geometries. Monte Carlo simulations of the XCAT phantom, the 2-meter sensitivity line source and positronium sensitivity phantoms were performed. We compared the sensitivity profiles for 2-gamma and 3-gamma tomography, relative cost of the setups and performed quantitative analysis of the reconstructed images. The analysis of the reconstructed XCAT images reveals the superiority of the seven-ring scanners over the three-ring setups. However, the three-ring scanners would be approximately 2-3 times cheaper. The peak sensitivity values for two-gamma vary from 20 to 34 cps/kBq. The sensitivity curves for the positronium tomography have a similar shape to the two-gamma sensitivity profiles. The peak values are lower compared to the two-gamma cases, from about 20-28 times, with a maximum of 1.66 cps/kBq. The results show the feasibility of multi-organ imaging of all the systems to be considered for the next generation of TB J-PET designs. The relative cost for all the scanners is about 10-4 times lower compared to the cost of the uExplorer. These properties coupled together with J-PET cost-effectiveness, make the J-PET technology an attractive solution for broad application in clinics.

BACKGROUND:Total-body (TB) Positron Emission Tomography (PET) is one of the most promising medical diagnostics modalities, opening new perspectives for personalized medicine, low-dose imaging, multi-organ dynamic imaging or kinetic modeling. The high sensitivity provided by total-body technology can be advantageous for novel tomography methods like positronium imaging, demanding the registration of triple coincidences. Currently, state-of-the-art PET scanners use inorganic scintillators. However, the high acquisition cost reduces the accessibility of TB PET technology. Several efforts are ongoing to mitigate this problem. Among the alternatives, the Jagiellonian PET (J-PET) technology, based on axially arranged plastic scintillator strips, offers a low-cost alternative solution for TB PET. PURPOSE:The work aimed to compare five total-body J-PET geometries with plastic scintillators suitable for multi-organ and positronium tomography as a possible next-generation J-PET scanner design. METHODS:We present comparative studies of performance characteristics of the cost-effective total-body PET scanners using J-PET technology. We investigated in silico five TB scanner geometries, varying the number of rings, scanner radii, and other parameters. Monte Carlo simulations of the anthropomorphic XCAT phantom, the extended 2-m sensitivity line source and positronium sensitivity phantoms were used to assess the performance of the geometries. Two hot spheres were placed in the lungs and in the liver of the XCAT phantom to mimic the pathological changes. We compared the sensitivity profiles and performed quantitative analysis of the reconstructed images by using quality metrics such as contrast recovery coefficient, background variability and root mean squared error. The studies are complemented by the determination of sensitivity for the positronium lifetime tomography and the relative cost analysis of the studied setups. RESULTS:The analysis of the reconstructed XCAT images reveals the superiority of the seven-ring scanners over the three-ring setups. However, the three-ring scanners would be approximately 2-3 times cheaper. The peak sensitivity values for two-gamma vary from 20 to 34 cps/kBq and are dominated by the differences in geometrical acceptance of the scanners. The sensitivity curves for the positronium tomography have a similar shape to the two-gamma sensitivity profiles. The peak values are lower compared to the two-gamma cases, from about 20-28 times, with a maximum value of 1.66 cps/kBq. This can be contrasted with the 50-cm one-layer J-PET modular scanner used to perform the first in-vivo positronium imaging with a sensitivity of 0.06 cps/kBq. CONCLUSIONS:The results show the feasibility of multi-organ imaging of all the systems to be considered for the next generation of TB J-PET designs. Among the scanner parameters, the most important ones are related to the axial field-of-view coverage. The two-gamma sensitivity and XCAT image reconstruction analyzes show the advantage of seven-ring scanners. However, the cost of the scintillator materials and SiPMs is more than two times higher for the longer modalities compared to the three-ring solutions. Nevertheless, the relative cost for all the scanners is about 10-4 times lower compared to the cost of the uExplorer. These properties coupled together with J-PET cost-effectiveness and triggerless acquisition mode enabling three-gamma positronium imaging, make the J-PET technology an attractive solution for broad application in clinics.

Discrete symmetries play an important role in particle physics with violation of CP connected to the matter-antimatter imbalance in the Universe. We report the most precise test of P, T and CP invariance in decays of ortho-positronium, performed with methodology involving polarization of photons from these decays. Positronium, the simplest bound state of an electron and positron, is of recent interest with discrepancies reported between measured hyperfine energy structure and theory at the level of 10-4 signaling a need for better understanding of the positronium system at this level. We test discrete symmetries using photon polarizations determined via Compton scattering in the dedicated J-PET tomograph on an event-by-event basis and without the need to control the spin of the positronium with an external magnetic field, in contrast to previous experiments. Our result is consistent with QED expectations at the level of 0.0007 and one standard deviation.

Objective: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J-PET) scanner for intra-treatment proton beam range monitoring. Approach: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries (three dual-heads and three cylindrical). We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread-Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J-PET scanner prototype dedicated to the proton beam range assessment. Main results: The investigations indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10$^{-5}$ coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple-layer dual-head geometry. Significance: We performed simulation studies demonstrating that the feasibility of the J-PET detector for PET-based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre-clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double-layer cylindrical and triple-layer dual-head J-PET geometry configurations seem promising for the future clinical application. Experimental tests are needed to confirm these findings.

Objective. The Jagiellonian positron emission tomography (J-PET) technology, based on plastic scintillators, has been proposed as a cost effective tool for detecting range deviations during proton therapy. This study investigates the feasibility of using J-PET for range monitoring by means of a detailed Monte Carlo simulation study of 95 patients who underwent proton therapy at the Cyclotron Centre Bronowice (CCB) in Krakow, Poland. Approach. Discrepancies between prescribed and delivered treatments were artificially introduced in the simulations by means of shifts in patient positioning and in the Hounsfield unit to the relative proton stopping power calibration curve. A dual-layer, cylindrical J-PET geometry was simulated in an in-room monitoring scenario and a triple-layer, dual-head geometry in an in-beam protocol. The distribution of range shifts in reconstructed PET activity was visualized in the beam's eye view. Linear prediction models were constructed from all patients in the cohort, using the mean shift in reconstructed PET activity as a predictor of the mean proton range deviation. Main results. Maps of deviations in the range of reconstructed PET distributions showed agreement with those of deviations in dose range in most patients. The linear prediction model showed a good fit, with coefficient of determination r (2) = 0.84 (in-room) and 0.75 (in-beam). Residual standard error was below 1 mm: 0.33 mm (in-room) and 0.23 mm (in-beam). Significance. The precision of the proposed prediction models shows the sensitivity of the proposed J-PET scanners to shifts in proton range for a wide range of clinical treatment plans. Furthermore, it motivates the use of such models as a tool for predicting proton range deviations and opens up new prospects for investigations into the use of intra-treatment PET images for predicting clinical metrics that aid in the assessment of the quality of delivered treatment.

In positron emission tomography (PET) studies, convolutional neural networks (CNNs) may be applied directly to the reconstructed distribution of radioactive tracers injected into the patient's body, as a pattern recognition tool. Nonetheless, unprocessed PET coincidence data exist in tabular format. This paper develops the transformation of tabular data into n-dimensional matrices, as a preparation stage for classification based on CNNs. This method explicitly introduces a nonlinear transformation at the feature engineering stage and then uses principal component analysis to create the images. We apply the proposed methodology to the classification of simulated PET coincidence events originating from NEMA IEC and anthropomorphic XCAT phantom. Comparative studies of neural network architectures, including multilayer perceptron and convolutional networks, were conducted. The developed method increased the initial number of features from 6 to 209 and gave the best precision results (79.8) for all tested neural network architectures; it also showed the smallest decrease when changing the test data to another phantom.

Total-Body PET imaging is one of the most promising newly introduced modalities in the medical diagnostics. State-of-the-art PET scanners use inorganic scintillators such as L(Y)SO or BGO, however, those technologies are very expensive, prohibitng the broad total-body PET applications. We present the comparative studies of performance characteristics of the cost-effective Total-Body PET scanners using Jagiellonian PET (J-PET) technology that is based on plastic scintillators. Here, we investigated in silico five realistic Total-Body scanner geometries, varying the number of rings, scanner radius, and distance between the neighbouring rings. Monte Carlo simulations of two NEMA phantoms (2-meter sensitivity line source and image quality) and the anthropomorphic XCAT phantom, were used to assess the performance of the tested geometries. We compared the sensitivity profiles and we performed the quantitative analysis of the reconstructed images by using the quality metrics such as contrast recovery coefficient, background variability and root mean squared error. The optimal scanner design was selected for the first Total-Body J-PET scanner configuration.

Objective. This paper reports on the implementation and shows examples of the use of the ProTheRaMon framework for simulating the delivery of proton therapy treatment plans and range monitoring using positron emission tomography (PET). ProTheRaMon offers complete processing of proton therapy treatment plans, patient CT geometries, and intra-treatment PET imaging, taking into account therapy and imaging coordinate systems and activity decay during the PET imaging protocol specific to a given proton therapy facility. We present the ProTheRaMon framework and illustrate its potential use case and data processing steps for a patient treated at the Cyclotron Centre Bronowice (CCB) proton therapy center in Krakow, Poland. Approach. The ProTheRaMon framework is based on GATE Monte Carlo software, the CASToR reconstruction package and in-house developed Python and bash scripts. The framework consists of five separated simulation and data processing steps, that can be further optimized according to the user’s needs and specific settings of a given proton therapy facility and PET scanner design. Main results. ProTheRaMon is presented using example data from a patient treated at CCB and the J-PET scanner to demonstrate the application of the framework for proton therapy range monitoring. The output of each simulation and data processing stage is described and visualized. Significance. We demonstrate that the ProTheRaMon simulation platform is a high-performance tool, capable of running on a computational cluster and suitable for multi-parameter studies, with databases consisting of large number of patients, as well as different PET scanner geometries and settings for range monitoring in a clinical environment. Due to its modular structure, the ProTheRaMon framework can be adjusted for different proton therapy centers and/or different PET detector geometries. It is available to the community via github (Borys et al 2022).