Comparison of optimized single and multifield irradiation plans of antiproton, proton and carbon ion beams

Robust Beam Selection Based on Water Equivalent Thickness Analysis in Passive Scattering Carbon-Ion Radiotherapy for Pancreatic Cancer

Carbon-ion radiotherapy (CIRT) is one of the most effective radiotherapeutic modalities. This study aimed to select robust-beam configurations (BC) by water equivalent thickness (WET) analysis in passive CIRT for pancreatic cancer. The study analyzed 110 computed tomography (CT) images and 600 dose distributions of eight patients with pancreatic cancer. The robustness in the beam range was evaluated using both planning and daily CT images, and two robust BCs for the rotating gantry and fixed port were selected. The planned, daily, and accumulated doses were calculated and compared after bone matching (BM) and tumor matching (TM). The dose-volume parameters for the target and organs at risk (OARs) were evaluated. Posterior oblique beams (120-240°) in the supine position and anteroposterior beams (0° and 180°) in the prone position were the most robust to WET changes. The mean CTV V95% reductions with TM were -3.8% and -5.2% with the BC for gantry and the BC for fixed ports, respectively. Despite ensuring robustness, the dose to the OARs increased slightly with WET-based BCs but remained below the dose constraint. The robustness of dose distribution can be improved by BCs that are robust to ΔWET. Robust BC with TM improves the accuracy of passive CIRT for pancreatic cancer.

Monte Carlo-Based Radiobiological Investigation of the Most Optimal Ion Beam Forming SOBP for Particle Therapy

Proton (p) and carbon (C) ion beams are in clinical use for cancer treatment, although other particles such as He, Be, and B ions have more recently gained attention. Identification of the most optimal ion beam for radiotherapy is a challenging task involving, among others, radiobiological characterization of a beam, which is depth-, energy-, and cell type- dependent. This study uses the FLUKA and MCDS Monte Carlo codes in order to estimate the relative biological effectiveness (RBE) for several ions of potential clinical interest such as p, ⁴He, ⁷Li, ¹⁰Be, ¹⁰B, and ¹²C forming a spread-out Bragg peak (SOBP). More specifically, an energy spectrum of the projectiles corresponding to a 5-cm SOBP at a depth of 8 cm was used. All secondary particles produced by the projectiles were considered and RBE was determined based on radiation-induced Double Strand Breaks (DSBs), as calculated by MCDS. In an attempt to identify the most optimal ion beam, using the latter data, biological optimization was performed and the obtained depth-dose distributions were inter-compared. The results showed that ¹²C ions are more effective inside the SOBP region, which comes at the expense of higher dose values at the tail (i.e., after the SOBP). In contrast, p beams exhibit a higher DSOPB/DEntrance ratio, if physical doses are considered. By performing a biological optimization in order to obtain a homogeneous biological dose (i.e., dose × RBE) in the SOBP, the corresponding advantages of p and ¹²C ions are moderated. ⁷Li ions conveniently combine a considerably lower dose tail and a DSOPB/DEntrance ratio similar to ¹²C. This work contributes towards identification of the most optimal ion beam for cancer therapy. The overall results of this work suggest that ⁷Li ions are of potential interest, although more studies are needed to demonstrate the relevant advantages. Future work will focus on studying more complex beam configurations.

Comparison of 126 MeV antiproton and proton—a FLUKA-based microdosimetric approach

Objective. This study aims at comparing dosimetric parameters of 126 MeV antiprotons and protons using microdosimetric approach. Approach. Microdosimetric distributions of 126 MeV proton and antiproton beams at 1 μm site size are calculated using the Monte Carlo-based FLUKA code. The distributions are calculated at various depths along the central axis in water phantom as well as at different off-axis locations. The study also includes calculations of secondary radiations produced by antiprotons and protons. Mean quality factor, Q ¯ is calculated using the ICRP 60 and ICRU 40 recommendations. The Relative Biological Effectiveness (RBE) of HSG tumour cell at 10% survival level is calculated based on Microdosimetric Kinetic Model. Main results. Q ¯ I C R P , Q ¯ ICRU and RBE for antiprotons are higher by a factor of about 3.60, 3.41 and 1.24, respectively, at Bragg-peak and higher by a factor of about 1.41, 1.76 and 1.05, respectively, at plateau region of depth-dose profile when compared to protons. At 15 cm depth along the central axis, Q ¯ ICRP , Q ¯ ICRU and RBE for protons are higher by a factor of about 1.42, 1.66 and 1.26, respectively, when compared to antiprotons. At the off-axis distance (Ld ) of 6 cm (at 11.5 cm depth in water), Q ¯ ICRP and Q ¯ ICRU of protons are higher than that of antiproton whereas the trend is opposite at off-axis distance of 4 cm. At Ld = 4 cm (at 11.5 cm depth in water), RBE of antiprotons is higher by about 4% than protons whereas at Ld = 6 cm, RBE of protons is higher by about 13% than antiprotons. Significance. The study shows that antiproton radiotherapy is advantageous as compared to protons considering enhancements in the absorbed dose and RBE-weighed dose values at the Bragg-peak.

Robust Angle Selection in Particle Therapy

The popularity of particle radiotherapy has grown exponentially over recent years owing to the marked advantage of the depth–dose curve and its unique biological property. However, particle therapy is sensitive to changes in anatomical structure, and the dose distribution may deteriorate. In particle therapy, robust beam angle selection plays a crucial role in mitigating inter- and intrafractional variation, including daily patient setup uncertainties and tumor motion. With the development of a rotating gantry, angle optimization has gained increasing attention. Currently, several studies use the variation in the water equivalent thickness to quantify anatomical changes during treatment. This method seems helpful in determining better beam angles and improving the robustness of planning. Therefore, this review will discuss and summarize the robust beam angles at different tumor sites in particle radiotherapy.

Monte Carlo simulations and dose measurements of 2D range-modulators for scanned particle therapy

This paper introduces the concept of a 2D range-modulator as a static device for generating spread-out Bragg peaks at very small distances to the target. The 2D range-modulator has some distinct advantages that can be highly useful for different research projects in particle therapy facilities. Most importantly, it creates an instantaneous, quasi-static irradiation field with only one energy, thus decreasing irradiation time tremendously. In addition, it can be manufactured fast and cost efficiently and its SOBP width and shape can be adjusted easily for the specific purpose/experiment. As the modulator is a static element, there is no need for rotation (e.g. like in a modulation wheel) or lateral oscillation and due to the small base structure period it can be positioned close to the target. Two different rapid prototyping manufacturing techniques were utilized. The modulation properties of one polymer and one steel modulator were investigated with both simulations and measurements. For this purpose, a sophisticated water phantom system (WERNER), that can perform fast, completely automated and high resolution dose measurements, was developed. Using WERNER, the dose distribution of a modulator can be verified quickly and reliably, both during experiments, as well as in a time constrained clinical environment. The maximum deviation between the Monte Carlo simulations and dose measurements in the spread-out Bragg peak region was 1.4% and 4% for the polymer and steel modulator respectively. They were able to create spread-out Bragg peaks with a high degree of dose homogeneity, thus validating the whole process chain, from the mathematical optimization and modulator development, to manufacturing, MC simulations and dose measurements. Combining the convenience, flexibility and cost-effectiveness of rapid prototyping with the advantages of highly customizable modulators, that can be adapted for different experiments, the 2D range-modulator is considered a very useful tool for a variety of research objectives. Moreover, we have successfully shown that the manufacturing of 2D modulators with high quality and high degree of homogeneity is possible, paving the way for the further development of the more complex 3D range-modulators, which are considered a viable option for the very fast treatment of moving targets and/or FLASH irradiation.

Comparison between bone matching and marker matching for evaluation of intra- and inter-fractional changes in accumulated dose of carbon ion radiotherapy for hepatocellular carcinoma

BACKGROUND AND PURPOSE

To determine whether bone matching (BM) or marker matching (MM) is the better positioning technique for carbon ion radiotherapy (CIRT) of primary hepatocellular carcinoma (HCC), we prospectively evaluated accumulated dose distributions with respect to intra- and inter-fractional anatomical changes.

MATERIALS AND METHODS

The accumulated doses in ten patients with HCC were evaluated, with the doses being calculated with respect to inter-fractional changes (InterDose) on treatment-room CT images on day 1 or day 2 of therapy (RefCT). This was accomplished by warping 3-day CT dose distributions to the RefCT through deformable registration. The accumulated doses were also calculated with respect to intra-fractional change (IntraDose) calculated by warping dose distributions for three 4DCT phases to the RefCT. Each dose was evaluated using dose-volume parameters for the clinical target volume (CTV) percentages receiving greater than 95% of the prescription dose (V₉₅).

RESULTS

The InterDose CTV V₉₅ values (mean [range]) were BM: 98.74% (95.62-100%), MM: 99.79% (98.55-100%), and the IntraDose values were BM: 99.46% (98.10-100%), MM: 99.74% (98.91-100%). Although all cases were acceptable with either matching method, MM provided better values than BM.

CONCLUSION

MM is a better positioning technique than BM for ensuring the target dose during and between fractions of CIRT. However, further analysis is required as our study included only a low number of cases.

Validation of new 2D ripple filters in proton treatments of spherical geometries and non-small cell lung carcinoma cases

A ripple filter (RiFi) is a passive energy modulator used in scanned particle therapy to broaden the Bragg peak, thus lowering the number of accelerator energies required for homogeneous target coverage, which significantly reduces the irradiation time. As we have previously shown, a new 6 mm thick RiFi with 2D groove shapes produced with 3D printing can be used in carbon ion treatments with a similar target coverage and only a marginally worse planning conformity compared to treatments with in-use 3 mm thick RiFis of an older 1D design. Where RiFis are normally not used with protons due to larger scattering and straggling effects, this new design would be beneficial in proton therapy too. Measurements of proton Bragg curves and lateral beam profiles were carried out for different RiFi designs and thicknesses as well as for no RiFi at the Heidelberg Ionenstrahl-Therapiezentrum. Base data for proton treatment planning were generated with the Monte Carlo code SHIELD-HIT12A with and without the 2D 6 mm RiFi. Plans on spherical targets in water were calculated with TRiP98 for a systematic RiFi performance analysis and for comparisons with carbon ion plans for the same respective energy depth step sizes. Plans for 9 stage I static non small cell lung cancer patients were calculated with Eclipse 13.7.15. Dose-volume-histograms, spatial dose distributions and dosimetric indexes were used for plan evaluation. Measurements confirm the functionality of the new 2D RiFi design, which reduces the beam spot size compared to 1D RiFis of the same thickness. Planning studies show that a 6 mm thick 2D RiFi could be used in proton therapy to lower the irradiation time. Although slightly worse planning conformity and dose homogeneity were found for plans with the RiFi compared to plans without, satisfactory results within the planning objective were obtained for all cases.

Technical Note: Predicting dose distribution with replacing stopping power ratio for inter-fractional motion and intra-fractional motion during carbon ion radiotherapy with passive irradiation method for stage I lung cancer

PURPOSE

We designed and evaluated a simple method for predicting the effects of intra-fractional and/or inter-fractional motion on dose distribution during carbon ion radiotherapy (CIRT) for solitary-lesion stage I lung cancer.

METHODS

The proposed method uses computed tomography (CT) images from treatment planning and intra-tumoral and/or inter-tumoral displacement. The predicted dose distribution (PDD) was calculated by replacing the current tumor region with the stopping power ratio (SPR) of the lung and replacing the moved tumor region with the SPR of the tumor. The actual dose distribution (ADD) was calculated without the replacement. Ten patients with solitary-lesion stage I lung cancer were retrospectively studied to evaluate the prediction method's accuracy. Four PDDs for intra-fractional motion (gate-in, exhalation, gate-out, inhalation phases during four-dimensional CT) and two PDDs for inter-fractional motion (CT images acquired 1-2 days before treatment) with bone- and tumor-matching methods were compared with each of six ADDs on each CT scan. Percentages of the planning/clinical target volumes (PTV/CTV) receiving >95% of the prescribed dose (V₉₅ ) and of minimum doses covering 95% of the PTV/CTV (D₉₅ ) were compared with dose volume histogram parameters.

RESULTS

The maximum tumor displacements occurred in the superior-inferior direction, with intra-fractional motion values of 3.75 and 8.97 mm for the superior and inferior directions, respectively, and inter-fractional values of 9.61 and 4.10 mm. The maximum average error for PTV V₉₅ regarding intra-fractional motion was -0.43% for the gate-out phase and -0.63% for the inhalation phase. There were no significant differences for these parameters (P = 0.541, P = 0.571). Average errors for PTV and CTV V₉₅ with inter-fractional motion with bone matching were 2.2% and 2.9%, respectively, with no significant differences (P = 0.387, P = 0.155).

CONCLUSIONS

The accuracy of the proposed method was good. Hence, it is feasible to use the proposed method during CIRT to predict dose distribution with respect to intra-fractional motion and/or inter-fractional motion of the tumor in patients with solitary-lesion stage I lung cancer.

Evaluation of the accuracy and clinical practicality of a calculation system for patient positional displacement in carbon ion radiotherapy at five sites

PURPOSE

We developed a system for calculating patient positional displacement between digital radiography images (DRs) and digitally reconstructed radiography images (DRRs) to reduce patient radiation exposure, minimize individual differences between radiological technologists in patient positioning, and decrease positioning time. The accuracy of this system at five sites was evaluated with clinical data from cancer patients. The dependence of calculation accuracy on the size of the region of interest (ROI) and initial position was evaluated for clinical use.

METHODS

For a preliminary verification, treatment planning and positioning data from eight setup patterns using a head and neck phantom were evaluated. Following this, data from 50 patients with prostate, lung, head and neck, liver, or pancreatic cancer (n = 10 each) were evaluated. Root mean square errors (RMSEs) between the results calculated by our system and the reference positions were assessed. The reference positions were manually determined by two radiological technologists to best-matching positions with orthogonal DRs and DRRs in six axial directions. The ROI size dependence was evaluated by comparing RMSEs for three different ROI sizes. Additionally, dependence on initial position parameters was evaluated by comparing RMSEs for four position patterns.

RESULTS

For the phantom study, the average (± standard deviation) translation error was 0.17 ± 0.05, rotation error was 0.17 ± 0.07, and ΔD was 0.14 ± 0.05. Using the optimal ROI size for each patient site, all cases of prostate, lung, and head and neck cancer with initial position parameters of 10 mm or under were acceptable in our tolerance. However, only four liver cancer cases and three pancreatic cancer cases were acceptable, because of low-reproducibility regions in the ROIs.

CONCLUSION

Our system has clinical practicality for prostate, lung, and head and neck cancer cases. Additionally, our findings suggest ROI size dependence in some cases.

3D range-modulator for scanned particle therapy: development, Monte Carlo simulations and experimental evaluation

The purpose of this work was to design and manufacture a 3D range-modulator for scanned particle therapy. The modulator is intended to create a highly conformal dose distribution with only one fixed energy, simultaneously reducing considerably the treatment time. As a proof of concept, a 3D range-modulator was developed for a spherical target volume with a diameter of 5 cm, placed at a depth of 25 cm in a water phantom. It consists of a large number of thin pins with a well-defined shape and different lengths to modulate the necessary shift of the Bragg peak. The 3D range-modulator was manufactured with a rapid prototyping technique. The FLUKA Monte Carlo package was used to simulate the modulating effect of the 3D range-modulator and the resulting dose distribution. For that purpose, a special user routine was implemented to handle its complex geometrical contour. Additionally, FLUKA was extended with the capability of intensity modulated scanning. To validate the simulation results, dose measurements were carried out at the Heidelberg Ion Beam Therapy Center with a 400.41 MeV/u ¹²C beam. The high resolution dosimetric measurements show a good agreement between simulated and measured dose distributions. Irradiation of the monoenergetic raster plan took 3 s, which is approximately 20 times shorter than a comparable plan with 16 different energies. The combination of only one energy and a 3D range-modulator leads to a tremendous decrease in irradiation time. 'Interplay effects', typical for moving targets and pencil beam scanning, can be immensely reduced or disappear completely, making the delivery of a homogeneous dose to moving targets more reliable. Combining high dose conformity, very good homogeneity and extremely short irradiation times, the 3D range-modulator is considered to become a clinically applicable method for very fast treatment of lung tumours.

Changes in Rectal Dose Due to Alterations in Beam Angles for Setup Uncertainty and Range Uncertainty in Carbon-Ion Radiotherapy for Prostate Cancer

Background and Purpose

Carbon-ion radiotherapy of prostate cancer is challenging in patients with metal implants in one or both hips. Problems can be circumvented by using fields at oblique angles. To evaluate the influence of setup and range uncertainties accompanying oblique field angles, we calculated rectal dose changes with oblique orthogonal field angles, using a device with fixed fields at 0° and 90° and a rotating patient couch.

Material and Methods

Dose distributions were calculated at the standard angles of 0° and 90°, and then at 30° and 60°. Setup uncertainty was simulated with changes from −2 mm to +2 mm for fields in the anterior-posterior, left-right, and cranial-caudal directions, and dose changes from range uncertainty were calculated with a 1 mm water-equivalent path length added to the target isocenter in each angle. The dose distributions regarding the passive irradiation method were calculated using the K2 dose algorithm.

Results

The rectal volumes with 0°, 30°, 60°, and 90° field angles at 95% of the prescription dose were 3.4±0.9 cm³, 2.8±1.1 cm³, 2.2±0.8 cm³, and 3.8±1.1 cm³, respectively. As compared with 90° fields, 30° and 60° fields had significant advantages regarding setup uncertainty and significant disadvantages regarding range uncertainty, but were not significantly different from the 90° field setup and range uncertainties.

Conclusions

The setup and range uncertainties calculated at 30° and 60° field angles were not associated with a significant change in rectal dose relative to those at 90°.

Antiproton induced DNA damage: proton like in flight, carbon-ion like near rest

Biological validation of new radiotherapy modalities is essential to understand their therapeutic potential. Antiprotons have been proposed for cancer therapy due to enhanced dose deposition provided by antiproton-nucleon annihilation. We assessed cellular DNA damage and relative biological effectiveness (RBE) of a clinically relevant antiproton beam. Despite a modest LET (~19 keV/μm), antiproton spread out Bragg peak (SOBP) irradiation caused significant residual γ-H2AX foci compared to X-ray, proton and antiproton plateau irradiation. RBE of ~1.48 in the SOBP and ~1 in the plateau were measured and used for a qualitative effective dose curve comparison with proton and carbon-ions. Foci in the antiproton SOBP were larger and more structured compared to X-rays, protons and carbon-ions. This is likely due to overlapping particle tracks near the annihilation vertex, creating spatially correlated DNA lesions. No biological effects were observed at 28–42 mm away from the primary beam suggesting minimal risk from long-range secondary particles.

Monte Carlo simulations of new 2D ripple filters for particle therapy facilities

At particle therapy facilities with pencil beam scanning, the implementation of a ripple filter (RiFi) broadens the Bragg peak (BP), which leads to fewer energy steps from the accelerator required to obtain an homogeneous dose coverage of the planned target volume (PTV). At the Universitätsklinikum Gießen und Marburg, Germany, a new second generation RiFi has been developed with two-dimensional groove structures. In this work we evaluate this new RiFi design.

METHODS

The Monte Carlo (MC) code SHIELD-HIT12A is used to determine the RiFi-induced inhomogeneities in the dose distribution for various ion types, initial particle energies and distances from the RiFi to the phantom surface as well as in the depth of the phantom. The beam delivery and monitor system (BAMS) used at Marburg, the Heidelberg Ionentherapiezentrum (HIT), Universitätsklinikum Heidelberg, Germany and the GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany is modeled and simulated. To evaluate the PTV dose coverage performance of the new RiFi design, the heavy ion treatment planning system TRiP98 is used for dose optimization. SHIELD-HIT12A is used to prepare the facility-specific physical dose kernels needed by TRiP, and for recalculating the physical dose distribution after TRiP optimization.

RESULTS

At short distances from the RiFi to the phantom surface fine structures in the dose distribution are observed. For various RiFis, ion types and initial particle energies the distance dmax at which maximum dose inhomogeneity occurs is found and an expression for dmax is deduced. The distance d0.01 at which the dose inhomogeneity is less than 1% is estimated and used as a threshold distance at which dose distributions are considered homogeneous. The MC data are found to agree with analytical expressions for dmax and d0.01; both are inversely related to the angular distribution. Increasing scatter from the beam delivery and monitoring system results in reduced dmax and d0.01. Furthermore, dmax and d0.01 are found to be proportional to the RiFi period λ.

CONCLUSION

Our findings clearly indicate that the dose inhomogeneity induced by RiFis does not add uncertainties to the dose distribution in the clinical setting. The new RiFi design can be used in treatments to obtain homogeneous PTV dose coverage with fewer energy steps while improving lateral penumbra, thereby reducing the required treatment time.

Dose and dose averaged LET comparison of ¹H, ⁴He, ⁶Li, ⁸Be, ¹⁰B, ¹²C, ¹⁴N, and ¹⁶O ion beams forming a spread-out Bragg peak

PURPOSE

Modern clinical accelerators are capable of producing ion beams from protons up to neon. This work compares the depth dose distribution and corresponding dose averaged linear energy transfer (LET) distribution, which is related to the biological effectiveness, for different ion beams (¹H, ⁴He, ⁶Li, ⁸Be, ¹⁰B, ¹²C, ¹⁴N, and ¹⁶O) using multi-energetic spectra in order to configure spread-out Bragg peaks (SOBP).

METHODS

Monte Carlo simulations were performed in order to configure a 5 cm SOBP at 8 cm depth in water for all the different ion beams. Physical dose and dose averaged LET distributions as a function of depth were then calculated and compared. The superposition of dose distribution of all ions is also presented for a two opposing fields configuration. Additional simulations were performed for (12)C beams to investigate the dependence of dose and dose averaged LET distributions on target depth and size, as well as beam configuration. These included simulations for a 3 cm SOBP at 7, 10, and 13 cm depth in water, a 6 cm SOBP at 7 depth in water, and two opposing fields of 6 cm SOBP.

RESULTS

Alpha particles and protons present superior physical depth dose distributions relative to the rest of the beams studied. Dose averaged LET distributions results suggest higher biological effectiveness in the target volume for carbon, nitrogen and oxygen ions. This is coupled, however, with relatively high LET values-especially for the last two ion species-outside the SOBP where healthy tissue would be located. Dose averaged LET distributions for ⁸Be and ¹⁰B beams show that they could be attractive alternatives to ¹²C for the treatment of small, not deeply seated lesions. The potential therapeutic effect of different ion beams studied in this work depends on target volume and position, as well as the number of beams used.

CONCLUSIONS

The optimization of beam modality for specific tumor cites remains an open question that warrants further investigation and clinically relevant results.

A semi-empirical model for the therapeutic range shift estimation caused by inhomogeneities in proton beam therapy

The purpose of this study was to devise a simple semi‐empirical model to estimate the range shift in clinical practices with high‐Z inhomogeneity in proton beam. A semi‐empirical model utilizing the logarithmic dependence on Z in stopping power from Bohr's classical approach has been developed to calculate the range shift due to the presence of inhomogeneity. Range shift from metallic plates of atomic number Z of various thicknesses were measured in water using a parallel plate ionization chamber and calculated with the FLUKA Monte Carlo code. The proton range shifts for bone and polymethyl methacrylate (PMMA) were estimated using the semi‐empirical model and compared with Monte Carlo calculation. The semi‐empirical equation to determine range shift and water equivalent thickness is presented. The model predicts a shift of the therapeutic range to within 2.5% accuracy for initial proton energies of 50 to 250 MeV and atomic numbers from 3.3 (effective Z for water) to 82. This equation is independent of beam energy, and thus provides range shift from high‐Z materials without the knowledge of proton energy. The proposed method of calculating the therapeutic range shift accurately requires only knowledge of the effective or actual atomic number of the inhomogeneity and the thickness of the inhomogeneity along the beam direction. The model generalizes the range shift calculation for any material based on its effective atomic number, and permits reliable prediction of the range shift for material combinations where no data is currently available. The proposed model can be readily implemented in routine clinical practice for proton range shift estimation and quality assurance on the treatment planning.

PACS numbers: 87.53.‐j, 87.55.‐x, 87.53.Bn, 87.55.D‐, 87.55.Qr, 87.55.K‐

Neutron fluence in antiproton radiotherapy, measurements and simulations

INTRODUCTION

A significant part of the secondary particle spectrum from antiproton annihilation consists of fast neutrons, which may contribute to a significant dose background found outside the primary beam.

MATERIALS AND METHODS

Using a polystyrene phantom as a moderator, we have performed absolute fluence measurements of the thermalized part of the fast neutron spectrum using Lithium-6 and -7 Fluoride TLD pairs. The results were compared with the Monte Carlo particle transport code FLUKA.

RESULTS

The experimental results are found to be in good agreement with simulations. The thermal neutron kerma resulting from the measured thermal neutron fluence is insignificant compared to the contribution from fast neutrons.

DISCUSSION

The secondary neutron fluences encountered in antiproton therapy are found to be similar to values calculated for pion treatment, however exact modeling under more realistic treatment scenarios is still required to quantitatively compare these treatment modalities.

Dose- and LET-painting with particle therapy

Tumour hypoxia is one of the limiting factors in obtaining tumour control in radiotherapy. The high-LET region of a beam of heavy charged particles such as carbon ions is located in the distal part of the Bragg peak. A modulated or spread out Bragg peak (SOBP) is a weighted function of several Bragg peaks at various energies, which however results in a dilution of the dose-average LET in the target volume. Here, we investigate the possibility to redistribute the LET by dedicated treatment plan optimisation, in order to maximise LET in the target volume. This may be a strategy to potentially overcome hypoxia along with dose escalation or dose painting. The high-LET region can be shaped in very different ways, while maintaining the distribution of the absorbed dose or biological effective dose. Treatment plans involving only carbon ion beams, show very different LET distributions depending on how the fields are arranged. Alternatively, a LET boost can be applied in multi-modal treatment planning, such as combining carbon ions with protons and/or photons. For such mixed radiation modalities, significant "LET boosts" can be achieved at nearly arbitrary positions within the target volume. Following the general understanding of the relationship between hypoxia, LET and the oxygen enhancement ratio (OER), we conclude, that an additional therapeutic advantage can be achieved by confining the high-LET part of the radiation in hypoxic compartments of the tumour, and applying low-LET radiation to the normoxic tissue. We also anticipate that additional advantages may be achieved by deliberate sparing of normal tissue from high LET regions. Consequently, treatment planning based on simultaneous dose and LET optimisation has a potential to achieve higher tumour control and/or reduced normal tissue control probability (NTCP).