arxiv: 2604.21173 · v2 · submitted 2026-04-23 · ⚛️ physics.med-ph

Recognition: unknown

Explicit positronium source modeling for Geant4 PET pipelines: controlled 2-gamma and 3-gamma generation and validation

Jason G Parker (1) ((1) Department of Radiology , Imaging Sciences , Indiana University School of Medicine)

Authors on Pith no claims yet

Pith reviewed 2026-05-08 13:23 UTC · model grok-4.3

classification ⚛️ physics.med-ph

keywords positronium modelingGeant4 simulationPET pipelines2-gamma annihilation3-gamma annihilationpositronium decaysimulation validation

0 comments

The pith

An explicit positronium source model in Geant4 gives users direct control over 2-gamma and 3-gamma annihilation channels, delays, and photon topology for PET simulations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper describes the addition of a controllable positronium source to Geant4 that lets researchers set the fractions of direct annihilation, para-positronium, and ortho-positronium events. It also allows choosing whether ortho-positronium decays into two or three photons, sampling delays, and adding prompt photons or positron range. Tests confirmed that the generated event fractions and lifetimes match the user inputs closely, with only modest extra computing cost. This matters for PET imaging research because accurate modeling of positronium behavior in tissue can improve how simulations predict real detector signals and help validate new imaging methods.

Core claim

The model was built to support direct annihilation, para-positronium, and ortho-positronium branches with user-defined fractions, explicit routing of ortho-positronium to two-photon or three-photon decay, exponential or fixed delay sampling, optional prompt-photon emission, and optional positron-range displacement. Event-level truth metadata were retained. In validation studies, the 2-gamma and 3-gamma fractions followed the requested parameters linearly, mean delays matched the lifetime settings, and computational cost increased linearly with event count.

What carries the argument

explicit positronium source model that routes events to different annihilation channels and decay modes while preserving metadata for validation

If this is right

2-gamma and 3-gamma fractions match requested values linearly
mean delays reproduce the set lifetimes with near one-to-one agreement
computational cost scales linearly with only modest overhead
downstream validation confirms correct handling of pure and mixed event datasets

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Researchers could use this to test how different positronium fractions affect image reconstruction in PET scanners
This framework might serve as a starting point for adding more detailed models of positronium formation in biological tissues
Simulations using this model could help design experiments to measure ortho-positronium lifetimes in different materials

Load-bearing premise

The new model integrates with Geant4's existing physics processes without creating inconsistencies in how photons are tracked or events are generated.

What would settle it

Running a large number of events with a specific ortho-positronium fraction set to produce three-photon decays and checking if the observed fraction deviates from the input value by more than a small margin would test the claim.

Figures

Figures reproduced from arXiv: 2604.21173 by Imaging Sciences, Indiana University School of Medicine), Jason G Parker (1) ((1) Department of Radiology.

**Figure 1.** Figure 1: shows that annihilation-channel control was recovered cleanly in both the native Geant4 reference sweep and the explicit positronium source model. In the native sweep, the observed 2γ fractions were 0.7537, 0.4984, and 0.2529 for requested control settings of 0.25, 0.50, and 0.75, with complementary 3γ fractions of 0.2463, 0.5016, and 0.7471, respectively. The end-member conditions were also recovered esse… view at source ↗

**Figure 3.** Figure 3: FIGURE 3 view at source ↗

read the original abstract

An explicit positronium (Ps) source model was implemented in Geant4 to provide direct event-level control over annihilation channel selection, decay timing, and photon emission topology. The implementation supports direct annihilation, para-positronium (p-Ps), and ortho-positronium (o-Ps) branches with user-defined fractions, explicit routing of o-Ps events to two-photon (2-gamma) or three-photon (3-gamma) decay, exponential or fixed delay sampling, optional prompt-photon emission, and optional positron-range displacement. Event-level truth metadata were retained to support downstream validation and analysis. The implementation was evaluated in controlled Geant4 studies using native reference configurations, explicit branch-fraction sweeps, lifetime sweeps, timing benchmarks, and a frozen point-source downstream test harness. Observed 2-gamma and 3-gamma fractions followed the requested control parameters with the expected linear behavior, and measured mean delays reproduced the prescribed lifetime settings with near one-to-one agreement. Computational cost scaled linearly with event count, with modest overhead relative to native Geant4 operation. A minimal downstream validation framework was used to verify branch-consistent handling of pure and mixed datasets, including expected method-source compatibility and recovery of valid events under unified 2-gamma and 3-gamma routing. These results establish a practical and internally consistent code framework for explicit positronium modeling in Geant4 and provide a simple pathway for PET researchers to incorporate controlled Ps generation into existing simulation pipelines.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

A narrow but practical Geant4 extension for explicit positronium control in PET simulations, with clean internal validation and no obvious interface bugs.

read the letter

The paper's main contribution is a Geant4 source model that lets users directly set fractions for direct annihilation, para-positronium, and ortho-positronium, then route the ortho events to either two or three photons while controlling delays and optional prompt photons. It also keeps per-event truth metadata so downstream code can verify what was generated. That combination of explicit routing and retained info is new for Geant4 PET pipelines, even if the underlying physics is standard.

Referee Report

0 major / 4 minor

Summary. The manuscript describes the implementation of an explicit positronium source model in Geant4 for PET pipelines. The model supports user-defined branching fractions for direct annihilation, para-positronium (p-Ps), and ortho-positronium (o-Ps), with explicit routing of o-Ps to 2-gamma or 3-gamma decays, exponential or fixed delay sampling, optional prompt-photon emission, and positron-range displacement. Event-level truth metadata is retained. Validation uses controlled studies (branch-fraction sweeps, lifetime sweeps, timing benchmarks, and a downstream mixed-dataset harness) showing that observed 2-gamma/3-gamma fractions follow input parameters linearly and mean delays match prescribed lifetimes with near one-to-one agreement, alongside linear computational scaling and pipeline compatibility.

Significance. If the reported internal consistency holds, this provides a practical extension to Geant4 for PET researchers requiring explicit control over positronium annihilation channels and timing, addressing gaps in native capabilities for multi-gamma and timing-sensitive simulations. The parameter-sweep validation approach directly tests interface behavior without fitted parameters, the retention of truth metadata supports downstream analysis, and the modest overhead plus downstream harness are strengths for reproducibility and integration into existing pipelines.

minor comments (4)

Abstract: the phrase 'native reference configurations' is used without definition or reference to specific Geant4 versions/physics lists; a brief clarification would aid readers in understanding the baseline comparison.
Results/Validation: while linear behavior and near one-to-one lifetime agreement are reported, the manuscript would benefit from a table or figure summarizing quantitative metrics (e.g., slope, intercept, or RMS deviation) across the parameter sweeps to make the agreement more precise.
The computational cost discussion notes linear scaling with modest overhead, but no specific benchmark values (e.g., CPU time per event or percentage overhead relative to native Geant4) are provided; including these would strengthen the practicality claim.
Consider adding a summary table of all user-controllable parameters (branch fractions, lifetimes, delay modes, etc.) with their allowed ranges and default behaviors for improved usability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of our manuscript on the explicit positronium source model in Geant4 and for recommending minor revision. The assessment correctly captures the implementation details, validation approach, and intended utility for PET pipelines. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity: direct input-output validation of source model

full rationale

The paper implements a controllable positronium source in Geant4 and validates it by sweeping user-specified inputs (branch fractions, lifetimes, delays) and confirming that simulated 2-gamma/3-gamma fractions and mean delays match those inputs with linear and near one-to-one agreement. This is standard implementation testing against prescribed parameters rather than any derivation, fitted prediction, or self-referential claim. No equations, self-citations, ansatzes, or uniqueness theorems appear in the load-bearing steps; the framework is self-contained as code plus direct consistency checks with native Geant4 behavior.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The work adds no new physical axioms or entities; it supplies control interfaces and validation harnesses over Geant4's existing positron annihilation and photon transport processes.

free parameters (2)

user-defined branch fractions
Explicitly set by the user rather than fitted by the paper.
lifetime and delay parameters
Prescribed by the user for each test configuration.

axioms (1)

domain assumption Geant4's native physics lists correctly model positronium formation, annihilation, and photon propagation
The new model routes events into these existing processes.

pith-pipeline@v0.9.0 · 5580 in / 1204 out tokens · 53092 ms · 2026-05-08T13:23:16.447841+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

20 extracted references

[1]

Introduction Positronium (Ps), the bound state of an electron and a positron, is a physically distinct annihilation intermediate whose decay topology depends on spin state. In the simplest picture relevant to PET- related modeling, para-positronium decays predominantly via two-photon (2γ) emission, whereas ortho-positronium decays predominantly via three-...
[2]

Methods 2.1 Positronium event generation A custom positronium (Ps) source model was implemented in Geant4 through an explicit event provider coupled to the primary generator. In this mode, each event was sampled from a three- component branching model comprising direct two-photon annihilation, para-positronium (p-Ps) two- photon decay, and ortho-positroni...
[3]

Results 3.1 Branching fraction validation Figure 1 shows that annihilation-channel control was recovered cleanly in both the native Geant4 reference sweep and the explicit positronium source model. In the native sweep, the observed 2γ fractions were 0.7537, 0.4984, and 0.2529 for requested control settings of 0.25, 0.50, and 0.75, with complementary 3γ fr...
[4]

The evaluations in this paper showed the requested controls were recovered quantitatively and reproducibly

Discussion A positronium-aware source model with explicit control of annihilation class, decay timing, and photon topology was implemented in Geant4 and validated against the expected branching and lifetime behavior (Figures 1–2). The evaluations in this paper showed the requested controls were recovered quantitatively and reproducibly. Reproducibility wa...
[5]

Downstream frozen-harness tests showed branch-consistent handling of explicit 2γ, explicit 3γ, and mixed datasets

Conclusion An explicit positronium source model with controllable branching, lifetime, prompt emission, and annihilation topology was implemented and validated within Geant4. Downstream frozen-harness tests showed branch-consistent handling of explicit 2γ, explicit 3γ, and mixed datasets. The framework provides a practical pathway for incorporating contro...
[6]

No human participants, human data, human tissue, or animals were involved

Declarations Ethical Approval and Consent to participate This study involved software development, simulation, and computational validation only. No human participants, human data, human tissue, or animals were involved. Ethical approval and consent to participate were therefore not required. Consent for publication The authors consent to publication of t...
[7]

Conti, M. and L. Eriksson, Physics of pure and non-pure positron emitters for PET: a review and a discussion. EJNMMI Phys, 2016. 3(1): p. 8

2016
[8]

Eur Phys J C Part Fields, 2016

Kaminska, D., et al., A feasibility study of ortho-positronium decays measurement with the J-PET scanner based on plastic scintillators. Eur Phys J C Part Fields, 2016. 76(8): p. 445

2016
[9]

Sci Adv, 2024

Moskal, P ., et al., Positronium image of the human brain in vivo. Sci Adv, 2024. 10(37): p. eadp2840

2024
[10]

Sci Adv, 2021

Moskal, P ., et al., Positronium imaging with the novel multiphoton PET scanner. Sci Adv, 2021. 7(42): p. eabh4394

2021
[11]

Phys Med Biol, 2019

Moskal, P ., et al., Feasibility study of the positronium imaging with the J-PET tomograph. Phys Med Biol, 2019. 64(5): p. 055017

2019
[12]

Roncali, E. and S.R. Cherry, Simulation of light transport in scintillators based on 3D characterization of crystal surfaces. Phys Med Biol, 2013. 58(7): p. 2185–98

2013
[13]

Stockhoff, and S.R

Roncali, E., M. Stockhoff, and S.R. Cherry, An integrated model of scintillator-reﬂector properties for advanced simulations of optical transport. Phys Med Biol, 2017. 62(12): p. 4811–4830

2017
[14]

Phys Med Biol, 2024

Salvadori, J., et al., PET digitization chain for Monte Carlo simulation in GATE. Phys Med Biol, 2024. 69(16)

2024
[15]

Phys Med Biol, 2022

Sarrut, D., et al., The OpenGATE ecosystem for Monte Carlo simulation in medical physics. Phys Med Biol, 2022. 67(18)

2022
[16]

Nuclear Instruments & Methods in Physics Research Section a-Accelerators Spectrometers Detectors and Associated Equipment, 2003

Agostinelli, S., et al., GEANT4-a simulation toolkit. Nuclear Instruments & Methods in Physics Research Section a-Accelerators Spectrometers Detectors and Associated Equipment, 2003. 506(3): p. 250–303

2003
[17]

Nuclear Instruments & Methods in Physics Research Section a-Accelerators Spectrometers Detectors and Associated Equipment, 2016

Allison, J., et al., Recent developments in GEANT4. Nuclear Instruments & Methods in Physics Research Section a-Accelerators Spectrometers Detectors and Associated Equipment, 2016. 835: p. 186–225

2016
[18]

Phys Med Biol, 2019

Efthimiou, N., et al., Implementation and validation of time-of-ﬂight PET image reconstruction module for listmode and sinogram projection data in the STIR library. Phys Med Biol, 2019. 64(3): p. 035004

2019
[19]

Phys Med Biol, 2012

Thielemans, K., et al., STIR: software for tomographic image reconstruction release 2. Phys Med Biol, 2012. 57(4): p. 867–83

2012
[20]

Stewart, C.A., Welch, V ., Plale, B., Fox, G, Pierce, M, Sterling, T, Indiana University Pervasive Technology Institute. 2017

2017