Monte Carlo (MC) simulations of beam interaction and transport in matter are increasingly considered as essential tools to support several aspects of radiation therapy. Despite the vast application of MC to photon therapy and scattered proton therapy, clinical experience in scanned ion beam therapy is still scarce. This is especially the case for ions heavier than protons, which pose additional issues like nuclear fragmentation and varying biological effectiveness. In this work, we present the evaluation of a dedicated framework which has been developed at the Heidelberg Ion Beam Therapy Center to provide automated FLUKA MC simulations of clinical patient treatments with scanned proton and carbon ion beams. Investigations on the number of transported primaries and the dimension of the geometry and scoring grids have been performed for a representative class of patient cases in order to provide recommendations on the simulation settings, showing that recommendations derived from the experience in proton therapy cannot be directly translated to the case of carbon ion beams. The MC results with the optimized settings have been compared to the calculations of the analytical treatment planning system (TPS), showing that regardless of the consistency of the two systems (in terms of beam model in water and range calculation in different materials) relevant differences can be found in dosimetric quantities and range, especially in the case of heterogeneous and deep seated treatment sites depending on the ion beam species and energies, homogeneity of the traversed tissue and size of the treated volume. The analysis of typical TPS speed-up approximations highlighted effects which deserve accurate treatment, in contrast to adequate beam model simplifications for scanned ion beam therapy. In terms of biological dose calculations, the investigation of the mixed field components in realistic anatomical situations confirmed the findings of previous groups so far reported only in homogenous water targets. This work can thus be useful to other centers commencing clinical experience in scanned ion beam therapy.
BackgroundIn proton radiation therapy a constant relative biological effectiveness (RBE) of 1.1 is usually assumed. However, biological experiments have evidenced RBE dependencies on dose level, proton linear energy transfer (LET) and tissue type. This work compares the predictions of three of the main radio-biological models proposed in the literature by Carabe-Fernandez, Wedenberg, Scholz and coworkers.MethodsUsing the chosen models, a spread-out Bragg peak (SOBP) as well as two exemplary clinical cases (single field and two fields) for cranial proton irradiation, all delivered with state-of-the-art pencil-beam scanning, have been analyzed in terms of absorbed dose, dose-averaged LET (LETD), RBE-weighted dose (DRBE) and biological range shift distributions.ResultsIn the systematic comparison of RBE predictions by the three models we could show different levels of agreement depending on (α/β)x and LET values. The SOBP study emphasizes the variation of LETD and RBE not only as a function of depth but also of lateral distance from the central beam axis. Application to clinical-like scenario shows consistent discrepancies from the values obtained for a constant RBE of 1.1, when using a variable RBE scheme for proton irradiation in tissues with low (α/β)x, regardless of the model. Biological range shifts of 0.6– 2.4 mm (for high (α/β)x) and 3.0 – 5.4 mm (for low (α/β)x) were found from the fall-off analysis of individual profiles of RBE-weighted fraction dose along the beam penetration depth.ConclusionsAlthough more experimental evidence is needed to validate the accuracy of the investigated models and their input parameters, their consistent trend suggests that their main RBE dependencies (dose, LET and (α/β)x) should be included in treatment planning systems. In particular, our results suggest that simpler models based on the linear-quadratic formalism and LETD might already be sufficient to reproduce important RBE dependencies for re-evaluation of plans optimized with the current RBE = 1.1 approximation. This approach would be a first step forward to consider RBE variations in proton therapy, thus enabling a more robust choice of biological dose delivery. The latter could in turn impact clinical outcome, especially in terms of reduced toxicities for tumors adjacent to organs at risk.
We introduce the automation of the range difference calculation deduced from particle-irradiation induced β(+)-activity distributions with the so-called most-likely-shift approach, and evaluate its reliability via the monitoring of algorithm- and patient-specific uncertainty factors. The calculation of the range deviation is based on the minimization of the absolute profile differences in the distal part of two activity depth profiles shifted against each other. Depending on the workflow of positron emission tomography (PET)-based range verification, the two profiles under evaluation can correspond to measured and simulated distributions, or only measured data from different treatment sessions. In comparison to previous work, the proposed approach includes an automated identification of the distal region of interest for each pair of PET depth profiles and under consideration of the planned dose distribution, resulting in the optimal shift distance. Moreover, it introduces an estimate of uncertainty associated to the identified shift, which is then used as weighting factor to 'red flag' problematic large range differences. Furthermore, additional patient-specific uncertainty factors are calculated using available computed tomography (CT) data to support the range analysis. The performance of the new method for in-vivo treatment verification in the clinical routine is investigated with in-room PET images for proton therapy as well as with offline PET images for proton and carbon ion therapy. The comparison between measured PET activity distributions and predictions obtained by Monte Carlo simulations or measurements from previous treatment fractions is performed. For this purpose, a total of 15 patient datasets were analyzed, which were acquired at Massachusetts General Hospital and Heidelberg Ion-Beam Therapy Center with in-room PET and offline PET/CT scanners, respectively. Calculated range differences between the compared activity distributions are reported in a 2D map in beam-eye-view. In comparison to previously proposed approaches, the new most-likely-shift method shows more robust results for assessing in-vivo the range from strongly varying PET distributions caused by differing patient geometry, ion beam species, beam delivery techniques, PET imaging concepts and counting statistics. The additional visualization of the uncertainties and the dedicated weighting strategy contribute to the understanding of the reliability of observed range differences and the complexity in the prediction of activity distributions. The proposed method promises to offer a feasible technique for clinical routine of PET-based range verification.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.