Benchmarking Geant4 hadronic models for prompt-gamma monitoring in carbon ion therapy

Purpose: The real-time monitoring of the spread-out Bragg peak would allow the planned dose delivered during treatment to be directly verified, but this poses a major challenge in modern ion beam therapy. A possible method to achieve this goal is to exploit the production of secondary particles by the nuclear reactions of the beam with the patient and correlate their emission profile to the planned target volume position. In this study, we present both the production rate and energy spectra of the prompt- produced by the interactions of the C-12 ion beam with a polymethyl methacrylate (PMMA) target. We also assess three different Monte Carlo models for prompt-simulation based on our experimental data. Methods: The experiment was carried out at the GSI Helmholtz Centre for Heavy Ion Research, Darmstadt, Germany with a 220MeV/u C-12 ions beam impinging on a 5x5x20cm(3) polymethyl methacrylate beam stopping target, with the prompt- being detected by a hexagonally-shaped barium fluoride scintillator with a circumscribed radius of 5.4cm and a length of 14cm, placed at 60 degrees and 90 degrees with respect to the beam direction. Monte Carlo simulations were carried out with three different hadronic models from the Geant4 code: binary ion cascade (BIC), quantum molecular dynamics (QMD), and Liege intranuclear cascade (INCL++). Results: An experimental prompt-yield of 1.06x10(-2)sr(-1) was measured at 90 degrees. A good agreement was observed between the shapes of the experimental and simulated energy spectra, especially with the INCL++ physics list. The prompt-gamma yield obtained with this physics list was compatible with the measurement within 2 sigma, with a relative difference of 26% on average. BIC and QMD physics lists proved to be less accurate than INCL++, with the difference between the measured and simulated yields exceeding 100%. The differences between the three physics lists were ascribed to important discrepancies between the models of the physical processes producing prompt- emissions. Conclusion: In conclusion, this study provides prompt- yield values in agreement with previously published results for different carbon ions energies. This work demonstrates that the INCL++ physics list from Geant4 is more accurate than BIC and QMD to reproduce prompt-gamma emission properties. (C) 2017 American Association of Physicists in Medicine