A Monte Carlo positronium decay source model with multiple annihilation channels in GATE

Authors:
Open Access
Publisher website
Google Scholar
Abstract
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 ( 2 γ / 3 γ ), 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 3 γ -to- 2 γ 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.
Funding Information
  • Ministère de l’Europe et des Affaires étrangères (PHC POLONIUM program (project number: 55194PB))
  • Fundacja na rzecz Nauki Polskiej (FIRST TEAM FENG.02.02-IP.05-0152/23)
  • Austrian Ministry for Women, Science and Research (WTZ PL 11/2026)
  • Narodowa Agencja Wymiany Akademickiej (BPN/BAT/2025/1/00009)
  • French Ministry for Higher Education, Research and Space (“PHC POLONIUM” program (project number: 55194)

This publication has 37 references indexed in Scilit: