Abstract:Background: Monte Carlo (MC) simulations are considered the gold-standard for accuracy in radiotherapy dose calculation; so far however, no commercial treatment planning system (TPS) provides a fast MC for supporting clinical practice in carbon ion therapy. Purpose: To extend and validate the in-house developed fast MC dose engine MonteRay for carbon ion therapy, including physical and biological dose calculation. Methods: MonteRay is a CPU MC dose calculation engine written in C++ that is capable of simulating therapeutic proton, helium and carbon ion beams. In this work, development steps taken to include carbon ions in MonteRay are presented. Dose distributions computed with MonteRay are evaluated using a comprehensive validation dataset, including various measurements (pristine Bragg peaks, spread out Bragg peaks in water and behind an anthropomorphic phantom) and simulations of a patient plan. The latter includes both physical and biological dose comparisons. Runtimes of MonteRay were evaluated against those of FLUKA MC on a standard benchmark problem. Results: Dosimetric comparisons between MonteRay and measurements demonstrated good agreement. In terms of pristine Bragg peaks, mean errors between simulated and measured integral depth dose distributions were between -2.3% and +2.7%. Comparing SOBPs at 5, 12.5 and 20 cm depth, mean absolute relative dose differences were 0.9%, 0.7% and 1.6% respectively. Comparison against measurements behind an anthropomorphic head phantom revealed mean absolute dose differences of 1.2 % ± 1.1 % $1.2\% \pm 1.1\;\%$ with global 3%/3 mm 3D-γ passing rates of 99.3%, comparable to those previously reached with FLUKA (98.9%). Comparisons against dose predictions computed with the clinical treatment planning tool RayStation 11B for a meningioma patient plan revealed excellent local 1%/1 mm 3D-γ passing rates of 98% for physical and 94% for biological dose. In terms of runtime, MonteRay achieved speedups against reference FLUKA simulations ranging from 14× to 72×, depending on the beam's energy and the step size chosen. Conclusions: Validations against clinical dosimetric measurements in homogeneous and heterogeneous scenarios and clinical TPS calculations have proven the validity of the physical models implemented in MonteRay. To conclude, MonteRay is viable as a fast secondary MC engine for supporting clinical practice in proton, helium and carbon ion radiotherapy.

Experimental validation of a FLUKA Monte Carlo simulation for carbon‐ion radiotherapy monitoring via secondary ion tracking

Experimental validation of a FLUKA Monte Carlo simulation for carbon-ion radiotherapy monitoring via secondary ion tracking

An in-vivo treatment monitoring system for ion-beam radiotherapy based on 28 Timepix3 detectors

Inhomogeneity detection within a head-sized phantom using tracking of charged nuclear fragments in ion beam therapy

Dosimetric verification in water of a Monte Carlo treatment planning tool for proton, helium, carbon and oxygen ion beams at the Heidelberg Ion Beam Therapy Center

Three dimensional reconstruction of therapeutic carbon ion beams in phantoms using single secondary ion tracks

Development and verification of an electron Monte Carlo engine for applications in intraoperative radiation therapy

First experimental verification of prompt gamma imaging with carbon ion irradiation

Monte Carlo investigation of sub-millimeter range verification in carbon ion radiation therapy using interaction vertex imaging

Development and validation of MonteRay, a fast Monte Carlo dose engine for carbon ion beam radiotherapy

Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study

Nanodosimetric investigation of the track structure of therapeutic carbon ion radiation. Part 1: Measurement of ionization cluster size distributions

Nanodosimetric investigation of the track structure of therapeutic carbon ion radiation. Part 2: Detailed radiation transport and track structure simulation

Evaluating an analytical prediction algorithm of positron emitter distributions in patient data for PET monitoring of carbon ion therapy: A simulation study

Mini-GRID radiotherapy on the CLEAR very-high-energy electron beamline: collimator optimization, film dosimetry, and Monte Carlo simulations

Analysis methods for in-beam PET images in proton therapy treatment verification: a comparison based on Monte Carlo simulations

Sub-mm Accuracy Results Measured from the First Prospective Clinical Trial of a Novel Real-Time IGRT System, Kilovoltage Intrafraction Monitoring (KIM)

Real-time delivered dose assessment in carbon ion therapy of moving targets

Intra- and Inter-Fraction Relative Range Verification in Heavy-Ion Therapy Using Filtered Interaction Vertex Imaging

Phase Space Generation for Proton and Carbon Ion Beams for External Users’ Applications at the Heidelberg Ion Therapy Center

Monte Carlo simulation study of an in vivo four-dimensional tracking system with a diverging collimator for monitoring radiation source (Ir-192) location during brachytherapy: proof of concept and feasibility