Abstract:Purpose: To compare the accuracy with which different hadronic inelastic physics models across ten Geant4 Monte Carlo simulation toolkit versions can predict positron-emitting fragments produced along the beam path during carbon and oxygen ion therapy. Materials and Methods: Phantoms of polyethylene, gelatin or poly(methyl methacrylate) were irradiated with monoenergetic carbon and oxygen ion beams. Post-irradiation, 4D PET images were acquired and parent 11 C, 10 C and 15 O radionuclides contributions in each voxel were determined from the extracted time activity curves. Next, the experimental configurations were simulated in Geant4 Monte Carlo versions 10.0 to 11.1, with three different fragmentation models - binary ion cascade (BIC), quantum molecular dynamics (QMD) and the Liege intranuclear cascade (INCL++) - 30 model-version combinations. Total positron annihilation and parent isotope production yields predicted by each simulation were compared between simulations and experiments using normalised mean squared error and Pearson cross-correlation coefficient. Finally, we compared the depth of maximum positron annihilation yield and the distal point at which positron yield decreases to 50% of peak between each model and the experimental results. Results: Performance varied considerably across versions and models, with no one version/model combination providing the best prediction of all positron-emitting fragments in all evaluated target materials and irradiation conidiations. BIC in Geant4 10.2 provided the best overall agreement with experimental results in the largest number of test cases. QMD consistently provided the best estimates of both the depth of peak positron yield (10.4 and 10.6) and the distal 50%-of-peak point (10.2), while BIC also performed well and INCL generally performed the worst across most Geant4 versions. Conclusions: Best spatial prediction of annihilation yield and positron-emitting fragment production during carbon and oxygen ion therapy was found to be 10.2.p03 with BIC or QMD. These version/model combinations are recommended for future heavy ion therapy research.

Simulation of Bragg Curves Produced by Protons, Alpha-Particles, and Carbon Ions in Water

Bragg-Curve Simulation of Carbon-Ion Beams for Particle-therapy Applications: a study with the GEANT4 toolkit

Proton and Carbon-Ion Minibeam Therapy: From Modeling to Treatment

Analysis of Bragg Curve Parameters and Lateral Straggle for Proton and Carbon Beams

Simulation of Proton and Carbon-12 Ion Beam for Tumor/Cancer Treatment

Proton relative biological effectiveness for the induction of DNA double strand breaks based on Geant4

SU‐GG‐T‐346: Accurate Simulation for Therapeutic Protons Range‐Energy and Range‐Density Tables

Comparative study of RBE and cell survival fractions for $^{1}$H, $^{4}$He, $^{12}$C and $^{16}$O beams using Geant4 and Microdosimetric Kinetic model

Modeling the Biophysical Effects in a Carbon Beam Delivery Line using Monte Carlo Simulation

SU‐FF‐T‐411: Boundary Study of Bragg Peak Shift and Bragg Peak Degradation in Proton Dose Calculation

Development of a more accurate Geant4 quantum molecular dynamics model for hadron therapy

Benchmarking geant4 nuclear models for hadron therapy with 95 MeV/nucleon carbon ions

MCHIT - Monte Carlo model for proton and heavy-ion therapy

A comparison of proton ranges in complex media using GATE/Geant4, MCNP6 and FLUKA

Ionization cross sections of water clusters as an alternative to liquid phase water: A Geant4-DNA simulation study

Development and verification of a Geant4 model of the electron beam mode in a clinical linear accelerator

In vitro modified microdosimetric kinetic model–based predictions for B14-150 cells survival in 450 MeV/u carbon ion beam with aluminum ridge filter for biologically optimized spread-out Bragg peak

Jump stochastic differential equations for the characterisation of the Bragg peak in proton beam radiotherapy

Parameterization of dose profiles of therapeutic minibeams of protons, $^{4}$He, $^{12}$C, and $^{16}$O

A quantitative assessment of Geant4 for predicting the yield and distribution of positron-emitting fragments in ion beam therapy