Data-driven RBE parameterization for helium ion beams

Helium ion beams are expected to be available again in the near future for clinical use. A suitable formalism to obtain relative biological effectiveness (RBE) values for treatment planning (TP) studies is needed. In this work we developed a data-driven RBE parameterization based on published in vitro experimental values. The RBE parameterization has been developed within the framework of the linear-quadratic (LQ) model as a function of the helium linear energy transfer (LET), dose and the tissue specific parameter (α/β)ph of the LQ model for the reference radiation. Analytic expressions are provided, derived from the collected database, describing the RBEα = αHe/αph and Rβ = βHe/βph ratios as a function of LET. Calculated RBE values at 2 Gy photon dose and at 10% survival (RBE10) are compared with the experimental ones. Pearson's correlation coefficients were, respectively, 0.85 and 0.84 confirming the soundness of the introduced approach. Moreover, due to the lack of experimental data at low LET, clonogenic experiments have been performed irradiating A549 cell line with (α/β)ph = 5.4 Gy at the entrance of a 56.4 MeV u(-1)He beam at the Heidelberg Ion Beam Therapy Center. The proposed parameterization reproduces the measured cell survival within the experimental uncertainties. A RBE formula, which depends only on dose, LET and (α/β)ph as input parameters is proposed, allowing a straightforward implementation in a TP system.

Ionoacoustic characterization of the proton Bragg peak with submillimeter accuracy

PURPOSE

Range verification in ion beam therapy relies to date on nuclear imaging techniques which require complex and costly detector systems. A different approach is the detection of thermoacoustic signals that are generated due to localized energy loss of ion beams in tissue (ionoacoustics). Aim of this work was to study experimentally the achievable position resolution of ionoacoustics under idealized conditions using high frequency ultrasonic transducers and a specifically selected probing beam.

METHODS

A water phantom was irradiated by a pulsed 20 MeV proton beam with varying pulse intensity and length. The acoustic signal of single proton pulses was measured by different PZT-based ultrasound detectors (3.5 and 10 MHz central frequencies). The proton dose distribution in water was calculated by Geant4 and used as input for simulation of the generated acoustic wave by the matlab toolbox k-WAVE.

RESULTS

In measurements from this study, a clear signal of the Bragg peak was observed for an energy deposition as low as 10(12) eV. The signal amplitude showed a linear increase with particle number per pulse and thus, dose. Bragg peak position measurements were reproducible within ±30 μm and agreed with Geant4 simulations to better than 100 μm. The ionoacoustic signal pattern allowed for a detailed analysis of the Bragg peak and could be well reproduced by k-WAVE simulations.

CONCLUSIONS

The authors have studied the ionoacoustic signal of the Bragg peak in experiments using a 20 MeV proton beam with its correspondingly localized energy deposition, demonstrating submillimeter position resolution and providing a deep insight in the correlation between the acoustic signal and Bragg peak shape. These results, together with earlier experiments and new simulations (including the results in this study) at higher energies, suggest ionoacoustics as a technique for range verification in particle therapy at locations, where the tumor can be localized by ultrasound imaging. This acoustic range verification approach could offer the possibility of combining anatomical ultrasound and Bragg peak imaging, but further studies are required for translation of these findings to clinical application.

Monte Carlo comparison of x-ray and proton CT for range calculations of proton therapy beams

Proton computed tomography (CT) has been described as a solution for imaging the proton stopping power of patient tissues, therefore reducing the uncertainty of the conversion of x-ray CT images to relative stopping power (RSP) maps and its associated margins. This study aimed to investigate this assertion under the assumption of ideal detection systems. We have developed a Monte Carlo framework to assess proton CT performances for the main steps of a proton therapy treatment planning, i.e. proton or x-ray CT imaging, conversion to RSP maps based on the calibration of a tissue phantom, and proton dose simulations. Irradiations of a computational phantom with pencil beams were simulated on various anatomical sites and the proton range was assessed on the reference, the proton CT-based and the x-ray CT-based material maps. Errors on the tissue's RSP reconstructed from proton CT were found to be significantly smaller and less dependent on the tissue distribution. The imaging dose was also found to be much more uniform and conformal to the primary beam. The mean absolute deviation for range calculations based on x-ray CT varies from 0.18 to 2.01 mm depending on the localization, while it is smaller than 0.1 mm for proton CT. Under the assumption of a perfect detection system, proton range predictions based on proton CT are therefore both more accurate and more uniform than those based on x-ray CT.

On the parametrization of lateral dose profiles in proton radiation therapy

PURPOSE

The accurate evaluation of the lateral dose profile is an important issue in the field of proton radiation therapy. The beam spread, due to Multiple Coulomb Scattering (MCS), is described by the Molière's theory. To take into account also the contribution of nuclear interactions, modern Treatment Planning Systems (TPSs) generally approximate the dose profiles by a sum of Gaussian functions. In this paper we have compared different parametrizations for the lateral dose profile of protons in water for therapeutical energies: the goal is to improve the performances of the actual treatment planning.

METHODS

We have simulated typical dose profiles at the CNAO (Centro Nazionale di Adroterapia Oncologica) beamline with the FLUKA code and validated them with data taken at CNAO considering different energies and depths. We then performed best fits of the lateral dose profiles for different functions using ROOT and MINUIT.

RESULTS

The accuracy of the best fits was analyzed by evaluating the reduced χ(2), the number of free parameters of the functions and the calculation time. The best results were obtained with the triple Gaussian and double Gaussian Lorentz-Cauchy functions which have 6 parameters, but good results were also obtained with the so called Gauss-Rutherford function which has only 4 parameters.

CONCLUSIONS

The comparison of the studied functions with accurate and validated Monte Carlo calculations and with experimental data from CNAO lead us to propose an original parametrization, the Gauss-Rutherford function, to describe the lateral dose profiles of proton beams.

Investigation of EBT2 and EBT3 films for proton dosimetry in the 4-20 MeV energy range

Radiochromic films such as Gafchromic EBT2 or EBT3 films are widely used for dose determination in radiation therapy because they offer a superior spatial resolution compared to any other digital dosimetric 2D detector array. The possibility to detect steep dose gradients is not only attractive for intensity-modulated radiation therapy with photons but also for intensity-modulated proton therapy. Their characteristic dose rate-independent response makes radiochromic films also attractive for dose determination in cell irradiation experiments using laser-driven ion accelerators, which are currently being investigated as future medical ion accelerators. However, when using these films in ion beams, the energy-dependent dose response in the vicinity of the Bragg peak has to be considered. In this work, the response of these films for low-energy protons is investigated. To allow for reproducible and background-free irradiation conditions, the films were exposed to mono-energetic protons from an electrostatic accelerator, in the 4-20 MeV energy range. For comparison, irradiation with clinical photons was also performed. It turned out that in general, EBT2 and EBT3 films show a comparable performance. For example, dose-response curves for photons and protons with energies as low as 11 MeV show almost no differences. However, corrections are required for proton energies below 11 MeV. Care has to be taken when correction factors are related to an average LET from depth-dose measurements, because only the dose-averaged LET yields similar results as obtained in mono-energetic measurements.

Comparing the biological washout of β+-activity induced in mice brain after 12C-ion and proton irradiation

In clinical ion beam therapy, protons as well as heavier ions such as carbon are used for treatment. For protons, β(+)-emitters are only induced by fragmentation reactions in the target (target fragmentation), whereas for heavy ions, they are additionally induced by fragmentations of the projectile (further referred to as autoactivation). An approach utilizing these processes for treatment verification, by comparing measured Positron Emission Tomography (PET) data to predictions from Monte Carlo simulations, has already been clinically implemented. For an accurate simulation, it is important to consider the biological washout of β(+)-emitters due to vital functions. To date, mathematical expressions for washout have mainly been determined by using radioactive beams of (10)C- and (11)C-ions, both β(+)-emitters, to enhance the counting statistics in the irradiated area. Still, the question of how the choice of projectile (autoactivating or non-autoactivating) influences the washout coefficients, has not been addressed. In this context, an experiment was carried out at the Heidelberg Ion Beam Therapy Center with the purpose of directly comparing irradiation-induced biological washout coefficients in mice for protons and (12)C-ions. To this aim, mice were irradiated in the brain region with protons and (12)C-ions and measured after irradiation with a PET/CT scanner (Siemens Biograph mCT). After an appropriate waiting time, the mice were sacrificed, then irradiated and measured again under similar conditions. The resulting data were processed and fitted numerically to deduce the main washout parameters. Despite the very low PET counting statistics, a consistent difference could be identified between (12)C-ion and proton irradiated mice, with the (12)C data being described best by a two component fit with a combined medium and slow washout fraction of 0.50 ± 0.05 and the proton mice data being described best by a one component fit with only one (slow) washout fraction of 0.73 ± 0.06.

Automation and uncertainty analysis of a method for in-vivo range verification in particle therapy

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.

A method to increase the nominal range resolution of a stack of parallel-plate ionization chambers

A detector prototype based on a stack of 61 parallel-plate ionisation chambers (PPIC) interleaved with absorber plates of polymethyl methacrylate (PMMA) was assembled for transmission imaging purposes in ion beam therapy. The thickness of the absorber sheets in the PPIC stack determines the nominal range resolution of the detector. In the current set-up, 3 mm PMMA slabs are used. The signal of the 61 active channels of the stack thereby provides a discrete approximation of the Bragg curve in the detector. In this work, a data processing method to increase the range resolution (MIRR) in a stack of ionization chambers is presented. In the MIRR the position of the maximum of the Bragg curve is deduced from the ratio of measured signals in adjacent PPIC channels. The method is developed based on Bragg curves obtained from Monte Carlo simulations and validated with experimental data of a wedge-shaped PMMA phantom acquired with the PPIC stack using carbon ion beams. The influence of the initial beam energy and of phantom inhomogeneities on the MIRR is quantitatively evaluated. Systematic errors as well as inaccuracies related to signal noise are discussed and quantified. It is shown that with the MIRR an increased range resolution of 0.7 mm PMMA equivalent or 0.8 mm water equivalent thickness is achieved for the considered experimental data.