Biological effectiveness of He-3 and He-4 ion beams for cancer hadrontherapy: a study based on the BIANCA biophysical model.

While cancer therapy with protons and C-ions is continuously spreading, in the near future patients will be also treated with He-ions which, in comparison to photons, combine the higher precision of protons with the higher relative biological effectiveness (RBE) of C-ions. Similarly to C-ions, also for He-ions the RBE variation along the beam must be known as precisely as possible, especially for active beam delivery systems. In this framework the BIANCA biophysical model, which has already been applied to calculate the RBE along proton and C-ion beams, was extended to4He-ions and, following interface with the FLUKA code, was benchmarked against cell survival data on CHO normal cells and Renca tumour cells irradiated at different positions along therapeutic-like4He-ion beams at the Heidelberg Ion-beam Therapy centre, where the first He-ion patient will be treated soon. Very good agreement between simulations and data was obtained, showing that BIANCA can now be used to predict RBE following irradiation with all ion types that are currently used, or will be used soon, for hadrontherapy. Thanks to the development of a reference simulation database describing V79 cell survival for ion and photon irradiation, these predictions can be cell-type specific because analogous databases can be produced, in principle, for any cell line. Furthermore, survival data on CHO cells irradiated by a He-3 beam were reproduced to compare the biophysical properties of He-4 and He-3 beams, which is currently an open question. This comparison showed that, at the same depth, He-4 beams tend to have a higher RBE with respect to He-3 beams, and that this difference is also modulated by the considered physical dose, as well as the cell radiosensitivity. However, at least for the considered cases, no significant difference was found for the ratio between the RBE-weighted dose in the SOBP and that in the entrance plateau.

Rapid effective dose calculation for raster-scanning 4He ion therapy with the modified microdosimetric kinetic model (mMKM).

PURPOSE: To develop and verify effective dose (DRBE) calculation in 4He ion beam therapy based on the modified microdosimetric kinetic model (mMKM) and evaluate the bio-sensitivity of mMKM-based plans to clinical parameters using a fast analytical dose engine. METHODS: Mixed radiation field particle spectra (MRFS) databases have been generated with Monte-Carlo (MC) simulations for 4He-ion beams. Relative biological effectiveness (RBE) and DRBE calculation using MRFS were established within a fast analytical engine. Spread-out Bragg-Peaks (SOBPs) in water were optimized for two dose levels and two tissue types with photon linear-quadratic model parameters αph, ßph, and (α/ß)ph to verify MRFS-derived database implementation against computations with MC-generated mixed-field α and ß databases. Bio-sensitivity of the SOBPs was investigated by varying absolute values of ßph, while keeping (α/ß)ph constant. Additionally, dose, dose-averaged linear energy transfer, and bio-sensitivity were investigated for two patient cases. RESULTS: Using MRFS-derived databases, dose differences ≲2% in the plateau and SOBP are observed compared to computations with MC-generated databases. Bio-sensitivity studies show larger deviations when altering the absolute ßph value, with maximum D50% changes of ~5%, with similar results for patient cases. Bio-sensitivity analysis indicates a greater impact on DRBE varying (α/ß)ph than ßph in mMKM. CONCLUSIONS: The MRSF approach yielded negligible differences in the target and small differences in the plateau compared to MC-generated databases. The presented analyses provide guidance for proper implementation of RBE-weighted 4He ion dose prescription and planning with mMKM. The MRFS-DRBE calculation approach using mMKM will be implemented in a clinical treatment planning system.

First benchmarking of the BIANCA model for cell survival prediction in a clinical hadron therapy scenario.

In the framework of RBE modelling for hadron therapy, the BIANCA biophysical model was extended to O-ions and was used to construct a radiobiological database describing the survival of V79 cells as a function of ion type (1 ⩽ Z ⩽ 8) and energy. This database allowed performing RBE predictions in very good agreement with experimental data. A method was then developed to construct analogous databases for different cell lines, starting from the V79 database as a reference. Following interface to the FLUKA Monte Carlo radiation transport code, BIANCA was then applied for the first time to predict cell survival in a typical patient treatment scenario, consisting of two opposing fields of range-equivalent protons or C-ions. The model predictions were found to be in good agreement with CHO cell survival data obtained at the Heidelberg ion-beam therapy (HIT) centre, as well as predictions performed by the local effect model (version LEM IV). This work shows that BIANCA can be used to predict cell survival and RBE not only for V79 and AG01522 cells, as shown previously, but also, in principle, for any cell line of interest. Furthermore, following interface to a transport code like FLUKA, BIANCA can provide predictions of 3D biological dose distributions for hadron therapy treatments, thus laying the foundations for future applications in clinics.

Determination of ion recombination and polarity effect correction factors for a plane-parallel ionization Bragg peak chamber under proton and carbon ion pencil beams.

Within the dosimetric characterization of particle beams, laterally-integrated depth-dose-distributions (IDDs) are measured and provided to the treatment planning system (TPS) for beam modeling or used as a benchmark for Monte Carlo (MC) simulations. The purpose of this work is the evaluation, in terms of ion recombination and polarity effect, of the dosimetric correction to be applied to proton and carbon ion curves as a function of linear energy transfer (LET). LET was calculated with a MC code for selected IDDs. Several regions of Bragg peak (BP) curve were investigated. The charge was measured with the plane-parallel BP-ionization chamber mounted in the Peakfinder as a field detector, by delivering a fixed number of particles at the maximum flux. The dose rate dependence was evaluated for different flux levels. The chamber was connected to an electrometer and exposed to un-scanned pencil beams. For each measurement the chamber was supplied with {±400, +200, +100} V. Recombination and polarity correction factors were then calculated as a function of depth and LET in water. Three energies representative of the clinical range were investigated for both particle types. The corrected IDDs (IDD k s) were then compared against MC. Recombination correction factors were LET and energy dependent, ranging from 1.000 to 1.040 (±0.5%) for carbon ions, while nearly negligible for protons. Moreover, no corrections need to be applied due to polarity effect being <0.5% along the whole IDDs for both particle types. IDD k s showed a better agreement than uncorrected curves when compared to MC, with a reduction of the mean absolute variation from 1.2% to 0.9%. The aforementioned correction factors were estimated and applied along the IDDs, showing an improved agreement against MC. Results confirmed that corrections are not negligible for carbon ions, particularly around the BP region.

Investigation of single carbon ion fragmentation in water and PMMA for hadron therapy.

Carbon ion radiotherapy is an attractive alternative to conventional radiotherapy, especially in case of deep-seated and radio-resistant tumors. As a consequence of inelastic nuclear reactions between primary particles and patient's tissues, the primary carbon ions may undergo nuclear fragmentation. The resulting decrease of primary ions and production of secondary fragments have to be carefully considered for accurate dose calculations in the treatment planning systems. The experimental data currently available provide only general information on carbon ion fragmentation and are not sufficient to cover the entire range of beam energies, target configurations and compositions relevant for radiotherapy. Therefore, new investigations were carried out to analyse the outcomes of the inelastic nuclear reaction processes on a single-ion-based approach. Measurements were performed at HIT, using 430 MeV/u carbon ion beams crossing water and PMMA targets. Unique in this method is the possibility of measuring number and type of fragments produced from each single carbon ion, provided that they are within the acceptance of the experimental apparatus. Concerning the amount of residual carbon ions behind water and PMMA targets with the same water equivalent thickness (WET), no significant differences were found. The experimental attenuation curve was well reproduced by the simulations. However, in the experiments, differences were observed regarding the amount of secondary fragments produced in water and in PMMA targets with the same WET. Differences were also found between experiments and simulations. These findings should be considered when dosimetric measurements are performed with PMMA instead of water phantoms. The found differences between experiments and simulations may contribute to improve the nuclear interaction and fragmentation models in Monte Carlo codes.

Sensitivity study of the microdosimetric kinetic model parameters for carbon ion radiotherapy.

In carbon ion therapy treatment planning, the relative biological effectiveness (RBE) is accounted for by optimization of the RBE-weighted dose (biological dose). The RBE calculation methods currently applied clinically in carbon ion therapy are derived from the microdosimetric kinetic model (MKM) in Japan and the local effect model (LEM) in Europe. The input parameters of these models are based on fit to experimental data subjected to uncertainties. We therefore performed a sensitivity study of the MKM input parameters, i.e. the domain radius (r d ), the nucleus radius (R n ) and the parameters of the linear quadratic (LQ) model (α x and ß). The study was performed with the FLUKA Monte Carlo code, using spread out Bragg peak (SOBP) scenarios in water and a biological dose distribution in a clinical patient case. Comparisons were done between biological doses estimated applying the MKM with parameters based on HSG cells, and with HSG parameters varied separately by ±{5, 25, 50}%. Comparisons were also done between parameter sets from different cell lines (HSG, V79, CHO and T1), as well as versions of the LEM. Of the parameters, r d had the largest impact on the biological dose distribution, especially on the absolute dose values. Increasing this parameter by 25% decreased the biological dose level at the center of a 3 Gy(RBE) SOBP by 14%. Variations in R n only influenced the biological dose distribution towards the particle range, and variations in α x resulted in minor changes in the biological dose, with an increasing impact towards the particle range. ß had the overall smallest influence on the SOBPs, but the impact could become more pronounced if alternative (LET dependent) implementations are used. The resulting percentage change in the SOBPs was generally less than the percentage change in the parameters. The patient case showed similar effects as with the SOBPs in water, and parameter variations had similar impact on the biological dose when using the clinical MKM and the general MKM. The clinical LEM calculated the highest biological doses to both tumor and surrounding healthy tissues, with a median target dose (D 50%) of 40.5 Gy(RBE), while the MKM with HSG and V79 parameters resulted in a D 50% of 34.2 and 36.9 Gy(RBE), respectively. In all, the observed change in biological dose distribution due to parameter variations demonstrates the importance of accurate input parameters when applying the MKM in treatment planning.

Monte Carlo simulation tool for online treatment monitoring in hadrontherapy with in-beam PET: A patient study.

Hadrontherapy is a method for treating cancer with very targeted dose distributions and enhanced radiobiological effects. To fully exploit these advantages, in vivo range monitoring systems are required. These devices measure, preferably during the treatment, the secondary radiation generated by the beam-tissue interactions. However, since correlation of the secondary radiation distribution with the dose is not straightforward, Monte Carlo (MC) simulations are very important for treatment quality assessment. The INSIDE project constructed an in-beam PET scanner to detect signals generated by the positron-emitting isotopes resulting from projectile-target fragmentation. In addition, a FLUKA-based simulation tool was developed to predict the corresponding reference PET images using a detailed scanner model. The INSIDE in-beam PET was used to monitor two consecutive proton treatment sessions on a patient at the Italian Center for Oncological Hadrontherapy (CNAO). The reconstructed PET images were updated every 10 s providing a near real-time quality assessment. By half-way through the treatment, the statistics of the measured PET images were already significant enough to be compared with the simulations with average differences in the activity range less than 2.5 mm along the beam direction. Without taking into account any preferential direction, differences within 1 mm were found. In this paper, the INSIDE MC simulation tool is described and the results of the first in vivo agreement evaluation are reported. These results have justified a clinical trial, in which the MC simulation tool will be used on a daily basis to study the compliance tolerances between the measured and simulated PET images.

MICRODOSIMETRIC SIMULATIONS OF CARBON IONS USING THE MONTE CARLO CODE FLUKA.

Therapeutic carbon ion beams produce a complex and variable radiation field that changes along the penetration depth due to the high density of energy loss along the particle track together with the secondary particles produced by nuclear fragmentation reactions. An accurate physical characterisation of such complex mixed-radiation fields can be performed by measuring microdosimetric spectra with mini tissue-equivalent proportional counters (mini-TEPCs), which are one of the most accurate devices used in experimental microdosimetry. Numerical calculations with Monte Carlo codes such as FLUKA can be used to supplement experimental microdosimetric measurements performed with TEPCs, but the nuclear cross sections and fragmentation models need to be benchmarked with experimental data for different energies and scenarios. The aim of this work is to compare experimental carbon microdosimetric data measured with the mini TEPC with calculated microdosimetry spectra obtained with FLUKA for 12C ions of 189.5 MeV/u in the Bragg peak region.

Dosimetric characterization of carbon fiber stabilization devices for post-operative particle therapy.

PURPOSE: The aim of this study was to evaluate the dosimetric impact caused by recently introduced carbon fiber reinforced polyetheretherketone (CF/PEEK) stabilization devices, in comparison with conventional titanium (Ti) implants, for post-operative particle therapy (PT). METHODS: As a first step, protons and carbon ions Spread-Out Bragg Peaks (SOBPs) were delivered to CF/PEEK and Ti screws. Transversal dose profiles were acquired with EBT3 films to evaluate beam perturbation. Effects on image quality and reconstruction artifacts were then investigated. CT scans of CF/PEEK and Ti implants were acquired according to our clinical protocol and Hounsfield Unit (HU) mean values were evaluated in three regions of interest. Implants and artifacts were then contoured in the sample CT scans, together with a target volume to simulate a spine tumor. Dose calculation accuracy was assessed by comparing optimized dose distributions with Monte Carlo simulations. In the end, the treatment plans of nine real patients (seven with CF/PEEK and two with Ti stabilization devices) were retrospectively analyzed to evaluate the dosimetric impact potentially occurring if improper management of the spine implant was carried out. RESULTS: As expected, CF/PEEK screw caused a very slight beam perturbation in comparison with Ti ones, leading to a lower degree of dose degradation in case of contouring and/or set-up uncertainties. Furthermore, CF/PEEK devices did not determine appreciable HU artifacts on CT images thus improving image quality and, as a final result, dose calculation accuracy. CONCLUSIONS: CF/PEEK spinal fixation devices resulted dosimetrically more suitable than commonly-used Ti implants for post-operative PT.

Fred: a GPU-accelerated fast-Monte Carlo code for rapid treatment plan recalculation in ion beam therapy.

Ion beam therapy is a rapidly growing technique for tumor radiation therapy. Ions allow for a high dose deposition in the tumor region, while sparing the surrounding healthy tissue. For this reason, the highest possible accuracy in the calculation of dose and its spatial distribution is required in treatment planning. On one hand, commonly used treatment planning software solutions adopt a simplified beam-body interaction model by remapping pre-calculated dose distributions into a 3D water-equivalent representation of the patient morphology. On the other hand, Monte Carlo (MC) simulations, which explicitly take into account all the details in the interaction of particles with human tissues, are considered to be the most reliable tool to address the complexity of mixed field irradiation in a heterogeneous environment. However, full MC calculations are not routinely used in clinical practice because they typically demand substantial computational resources. Therefore MC simulations are usually only used to check treatment plans for a restricted number of difficult cases. The advent of general-purpose programming GPU cards prompted the development of trimmed-down MC-based dose engines which can significantly reduce the time needed to recalculate a treatment plan with respect to standard MC codes in CPU hardware. In this work, we report on the development of fred, a new MC simulation platform for treatment planning in ion beam therapy. The code can transport particles through a 3D voxel grid using a class II MC algorithm. Both primary and secondary particles are tracked and their energy deposition is scored along the trajectory. Effective models for particle-medium interaction have been implemented, balancing accuracy in dose deposition with computational cost. Currently, the most refined module is the transport of proton beams in water: single pencil beam dose-depth distributions obtained with fred agree with those produced by standard MC codes within 1-2% of the Bragg peak in the therapeutic energy range. A comparison with measurements taken at the CNAO treatment center shows that the lateral dose tails are reproduced within 2% in the field size factor test up to 20 cm. The tracing kernel can run on GPU hardware, achieving 10 million primary [Formula: see text] on a single card. This performance allows one to recalculate a proton treatment plan at 1% of the total particles in just a few minutes.

Investigation of mixed ion fields in the forward direction for 220.5 MeV/u helium ion beams: comparison between water and PMMA targets.

Currently there is a rising interest in helium ion beams for radiotherapy. For benchmarking of the physical beam models used in treatment planning, there is a need for experimental data on the composition and spatial distribution of mixed ion fields. Of particular interest are the attenuation of the primary helium ion fluence and the build-up of secondary hydrogen ions due to nuclear interactions. The aim of this work was to provide such data with an enhanced precision. Moreover, the validity and limits of the mixed ion field equivalence between water and PMMA targets were investigated. Experiments with a 220.5 MeV/u helium ion pencil beam were performed at the Heidelberg Ion-Beam Therapy Center in Germany. The compact detection system used for ion tracking and identification was solely based on Timepix position-sensitive semiconductor detectors. In comparison to standard techniques, this system is two orders of magnitude smaller, and provides higher precision and flexibility. The numbers of outgoing helium and hydrogen ions per primary helium ion as well as the lateral particle distributions were quantitatively investigated in the forward direction behind water and PMMA targets with 5.2-18 cm water equivalent thickness (WET). Comparing water and PMMA targets with the same WET, we found that significant differences in the amount of outgoing helium and hydrogen ions and in the lateral particle distributions arise for target thicknesses above 10 cm WET. The experimental results concerning hydrogen ions emerging from the targets were reproduced reasonably well by Monte Carlo simulations using the FLUKA code. Concerning the amount of outgoing helium ions, significant differences of 3-15% were found between experiments and simulations. We conclude that if PMMA is used in place of water in dosimetry, differences in the dose distributions could arise close to the edges of the field, in particular for deep seated targets.

Helium ions at the heidelberg ion beam therapy center: comparisons between FLUKA Monte Carlo code predictions and dosimetric measurements.

In the field of particle therapy helium ion beams could offer an alternative for radiotherapy treatments, owing to their interesting physical and biological properties intermediate between protons and carbon ions. We present in this work the comparisons and validations of the Monte Carlo FLUKA code against in-depth dosimetric measurements acquired at the Heidelberg Ion Beam Therapy Center (HIT). Depth dose distributions in water with and without ripple filter, lateral profiles at different depths in water and a spread-out Bragg peak were investigated. After experimentally-driven tuning of the less known initial beam characteristics in vacuum (beam lateral size and momentum spread) and simulation parameters (water ionization potential), comparisons of depth dose distributions were performed between simulations and measurements, which showed overall good agreement with range differences below 0.1 mm and dose-weighted average dose-differences below 2.3% throughout the entire energy range. Comparisons of lateral dose profiles showed differences in full-width-half-maximum lower than 0.7 mm. Measurements of the spread-out Bragg peak indicated differences with simulations below 1% in the high dose regions and 3% in all other regions, with a range difference less than 0.5 mm. Despite the promising results, some discrepancies between simulations and measurements were observed, particularly at high energies. These differences were attributed to an underestimation of dose contributions from secondary particles at large angles, as seen in a triple Gaussian parametrization of the lateral profiles along the depth. However, the results allowed us to validate FLUKA simulations against measurements, confirming its suitability for 4He ion beam modeling in preparation of clinical establishment at HIT. Future activities building on this work will include treatment plan comparisons using validated biological models between proton and helium ions, either within a Monte Carlo treatment planning engine based on the same FLUKA code, or an independent analytical planning system fed with a validated database of inputs calculated with FLUKA.

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.

The introduction of 'new' ion species in particle therapy needs to be supported by a thorough assessment of their dosimetric properties and by treatment planning comparisons with clinically used proton and carbon ion beams. In addition to the latter two ions, helium and oxygen ion beams are foreseen at the Heidelberg Ion Beam Therapy Center (HIT) as potential assets for improving clinical outcomes in the near future. We present in this study a dosimetric validation of a FLUKA-based Monte Carlo treatment planning tool (MCTP) for protons, helium, carbon and oxygen ions for spread-out Bragg peaks in water. The comparisons between the ions show the dosimetric advantages of helium and heavier ion beams in terms of their distal and lateral fall-offs with respect to protons, reducing the lateral size of the region receiving 50% of the planned dose up to 12 mm. However, carbon and oxygen ions showed significant doses beyond the target due to the higher fragmentation tail compared to lighter ions (p and He), up to 25%. The Monte Carlo predictions were found to be in excellent geometrical agreement with the measurements, with deviations below 1 mm for all parameters investigated such as target and lateral size as well as distal fall-offs. Measured and simulated absolute dose values agreed within about 2.5% on the overall dose distributions. The MCTP tool, which supports the usage of multiple state-of-the-art relative biological effectiveness models, will provide a solid engine for treatment planning comparisons at HIT.

Ion recombination correction factor in scanned light-ion beams for absolute dose measurement using plane-parallel ionisation chambers.

Based on international reference dosimetry protocols for light-ion beams, a correction factor (k s) has to be applied to the response of a plane-parallel ionisation chamber, to account for recombination of negative and positive charges in its air cavity before these charges can be collected on the electrodes. In this work, k s for IBA PPC40 Roos-type chambers is investigated in four scanned light-ion beams (proton, helium, carbon and oxygen). To take into account the high dose-rates used with scanned beams and LET-values, experimental results are compared to a model combining two theories. One theory, developed by Jaffé, describes the variation of k s with the ionization density within the ion track (initial recombination) and the other theory, developed by Boag, describes the variation of k s with the dose rate (volume recombination). Excellent agreement is found between experimental and theoretical k s-values. All results confirm that k s cannot be neglected. The solution to minimise k s is to use the ionisation chamber at high voltage. However, one must be aware that charge multiplication may complicate the interpretation of the measurement. For the chamber tested, it was found that a voltage of 300 V can be used without further complication. As the initial recombination has a logarithmic variation as a function of 1/V, the two-voltage method is not applicable to these scanned beams.

Optimizing the modified microdosimetric kinetic model input parameters for proton and 4He ion beam therapy application.

Models able to predict relative biological effectiveness (RBE) values are necessary for an accurate determination of the biological effect with proton and 4He ion beams. This is particularly important when including RBE calculations in treatment planning studies comparing biologically optimized proton and 4He ion beam plans. In this work, we have tailored the predictions of the modified microdosimetric kinetic model (MKM), which is clinically applied for carbon ion beam therapy in Japan, to reproduce RBE with proton and 4He ion beams. We have tuned the input parameters of the MKM, i.e. the domain and nucleus radii, reproducing an experimental database of initial RBE data for proton and He ion beams. The modified MKM, with the best fit parameters obtained, has been used to reproduce in vitro cell survival data in clinically-relevant scenarios. A satisfactory agreement has been found for the studied cell lines, A549 and RENCA, with the mean absolute survival variation between the data and predictions within 2% and 5% for proton and 4He ion beams, respectively. Moreover, a sensitivity study has been performed varying the domain and nucleus radii and the quadratic parameter of the photon response curve. The promising agreement found in this work for the studied clinical-like scenarios supports the usage of the modified MKM for treatment planning studies in proton and 4He ion beam therapy.

Experimental dosimetric comparison of 1H, 4He, 12C and 16O scanned ion beams.

At the Heidelberg Ion Beam Therapy Center, scanned helium and oxygen ion beams are available in addition to the clinically used protons and carbon ions for physical and biological experiments. In this work, a study of the basic dosimetric features of the different ions is performed in the entire therapeutic energy range. Depth dose distributions are investigated for pencil-like beam irradiation, with and without a modulating ripple filter, focusing on the extraction of key Bragg curve parameters, such as the range, the peak-width and the distal 80%-20% fall-off. Pencil-beam lateral profiles are measured at different depths in water, and parameterized with multiple Gaussian functions. A more complex situation of an extended treatment field is analyzed through a physically optimized spread-out Bragg peak, delivered with beam scanning. The experimental results of this physical beam characterization indicate that helium ions could afford a more conformal treatment and in turn, increased tumor control. This is mainly due to a smaller lateral scattering than with protons, leading to better lateral and distal fall-off, as well as a lower fragmentation tail compared to carbon and oxygen ions. Moreover, the dosimetric dataset can be used directly for comparison with results from analytical dose engines or Monte Carlo codes. Specifically, it was used at the Heidelberg Ion Beam Therapy Center to generate a new input database for a research analytical treatment planning system, as well as for validation of a general purpose Monte Carlo program, in order to lay the groundwork for biological experiments and further patient planning studies.

The FLUKA Monte Carlo code coupled with the NIRS approach for clinical dose calculations in carbon ion therapy.

Particle therapy facilities often require Monte Carlo (MC) simulations to overcome intrinsic limitations of analytical treatment planning systems (TPS) related to the description of the mixed radiation field and beam interaction with tissue inhomogeneities. Some of these uncertainties may affect the computation of effective dose distributions; therefore, particle therapy dedicated MC codes should provide both absorbed and biological doses. Two biophysical models are currently applied clinically in particle therapy: the local effect model (LEM) and the microdosimetric kinetic model (MKM). In this paper, we describe the coupling of the NIRS (National Institute for Radiological Sciences, Japan) clinical dose to the FLUKA MC code. We moved from the implementation of the model itself to its application in clinical cases, according to the NIRS approach, where a scaling factor is introduced to rescale the (carbon-equivalent) biological dose to a clinical dose level. A high level of agreement was found with published data by exploring a range of values for the MKM input parameters, while some differences were registered in forward recalculations of NIRS patient plans, mainly attributable to differences with the analytical TPS dose engine (taken as reference) in describing the mixed radiation field (lateral spread and fragmentation). We presented a tool which is being used at the Italian National Center for Oncological Hadrontherapy to support the comparison study between the NIRS clinical dose level and the LEM dose specification.

A phenomenological relative biological effectiveness approach for proton therapy based on an improved description of the mixed radiation field.

Proton therapy treatment planning systems (TPSs) are based on the assumption of a constant relative biological effectiveness (RBE) of 1.1 without taking into account the found in vitro experimental variations of the RBE as a function of tissue type, linear energy transfer (LET) and dose. The phenomenological RBE models available in literature are based on the dose-averaged LET (LET D ) as an indicator of the physical properties of the proton radiation field. The LET D values are typically calculated taking into account primary and secondary protons, neglecting the biological effect of heavier secondaries. In this work, we have introduced a phenomenological RBE approach which considers the biological effect of primary protons, and of secondary protons, deuterons, tritons (Z = 1) and He fragments (3He and 4He, Z = 2). The calculation framework, coupled with a Monte Carlo (MC) code, has been successfully benchmarked against clonogenic in vitro data measured in this work for two cell lines and then applied to determine biological quantities for spread-out Bragg peaks and a prostate and a head case. The introduced RBE formalism, which depends on the mixed radiation field, the dose and the ratio of the linear-quadratic model parameters for the reference radiation [Formula: see text], predicts, when integrated in an MC code, higher RBE values in comparison to LET D -based parameterizations. This effect is particular enhanced in the entrance channel of the proton field and for low [Formula: see text] tissues. For the prostate and the head case, we found higher RBE-weighted dose values up to about 5% in the entrance channel when including or neglecting the Z = 2 secondaries in the RBE calculation. TPSs able to proper account for the mixed radiation field in proton therapy are thus recommended for an accurate determination of the RBE in the whole treatment field.

Characterization of a commercial scintillation detector for 2-D dosimetry in scanned proton and carbon ion beams.

INTRODUCTION: Pencil beam scanning technique used at CNAO requires beam characteristics to be carefully assessed and periodically checked to guarantee patient safety. This study aimed at characterizing the Lynx® detector (IBA Dosimetry) for commissioning and periodic quality assurance (QA) for proton and carbon ion beams, as compared to EBT3 films, currently used for QA checks. METHODS AND MATERIALS: The Lynx® is a 2-D high-resolution dosimetry system consisting of a scintillating screen coupled with a CCD camera, in a compact light-tight box. The scintillator was preliminarily characterized in terms of short-term stability, linearity with number of particles, image quality and response dependence on iris setting and beam current; Lynx® was then systematically tested against EBT3 films. The detector response dependence on radiation LET was also assessed. RESULTS: Preliminary results have shown that Lynx is suitable to be used for commissioning and QA checks for proton and carbon ion scanning beams; the cross-check with EBT3 films showed a good agreement between the two detectors, for both single spot and scanned field measurements. The strong LET dependence of the scintillator due to quenching effect makes Lynx® suitable only for relative 2-D dosimetry measurements. CONCLUSION: Lynx® appears as a promising tool for commissioning and periodic QA checks for both protons and carbon ion beams. This detector can be used as an alternative of EBT3 films, allowing real-time measurements and analysis, with a significant time sparing.

Biologically optimized helium ion plans: calculation approach and its in vitro validation.

Treatment planning studies on the biological effect of raster-scanned helium ion beams should be performed, together with their experimental verification, before their clinical application at the Heidelberg Ion Beam Therapy Center (HIT). For this purpose, we introduce a novel calculation approach based on integrating data-driven biological models in our Monte Carlo treatment planning (MCTP) tool. Dealing with a mixed radiation field, the biological effect of the primary (4)He ion beams, of the secondary (3)He and (4)He (Z = 2) fragments and of the produced protons, deuterons and tritons (Z = 1) has to be taken into account. A spread-out Bragg peak (SOBP) in water, representative of a clinically-relevant scenario, has been biologically optimized with the MCTP and then delivered at HIT. Predictions of cell survival and RBE for a tumor cell line, characterized by [Formula: see text] Gy, have been successfully compared against measured clonogenic survival data. The mean absolute survival variation ([Formula: see text]) between model predictions and experimental data was 5.3% ± 0.9%. A sensitivity study, i.e. quantifying the variation of the estimations for the studied plan as a function of the applied phenomenological modelling approach, has been performed. The feasibility of a simpler biological modelling based on dose-averaged LET (linear energy transfer) has been tested. Moreover, comparisons with biophysical models such as the local effect model (LEM) and the repair-misrepair-fixation (RMF) model were performed. [Formula: see text] values for the LEM and the RMF model were, respectively, 4.5% ± 0.8% and 5.8% ± 1.1%. The satisfactorily agreement found in this work for the studied SOBP, representative of clinically-relevant scenario, suggests that the introduced approach could be applied for an accurate estimation of the biological effect for helium ion radiotherapy.