Linear Energy Transfer (LET) is widely used to express the radiation quality of ion beams, when characterizing the biological effectiveness. However, averaged LET may be defined in multiple ways, and the chosen definition may impact the resulting reported value. We review averaged LET definitions found in the literature, and quantify which impact using these various definitions have for different reference setups. We recorded the averaged LET definitions used in 354 publications quantifying the relative biological effectiveness (RBE) of hadronic beams, and investigated how these various definitions impact the reported averaged LET using a Monte Carlo particle transport code. We find that the kind of averaged LET being applied is, generally, poorly defined. Some definitions of averaged LET may influence the reported averaged LET values up to an order of magnitude. For publications involving protons, most applied dose averaged LET when reporting RBE. The absence of what target medium is used and what secondary particles are included further contributes to an ill-defined averaged LET. We also found evidence of inconsistent usage of averaged LET definitions when deriving LET-based RBE models. To conclude, due to commonly ill-defined averaged LET and to the inherent problems of LET-based RBE models, averaged LET may only be used as a coarse indicator of radiation quality. We propose a more rigorous way of reporting LET values, and suggest that ideally the entire particle fluence spectra should be recorded and provided for future RBE studies, from which any type of averaged LET (or other quantities) may be inferred.

Ion beams passing through lung tissue show more pronounced energy straggling than expected for solid materials of the same thickness. Energy straggling in active scanning particle therapy can enlarge the pencil beam Bragg peaks in-depth as well as displace them, deteriorating the dose coverage of a target within the lung. While this is not yet considered in any known treatment planning system, we implement a mathematical model to be used for treatment planning, using TRiP98, which relies on a deterministic pencil beam algorithm. Through a randomization process based on a continuous Poisson probability distribution, the HU values of lung voxels are replaced with a modified value in successive iterations. The beam-modulation effect of the lung can thus be taken into account in treatment planning by recalculating the dose n times for n randomized CTs using the raster scan file of a plan that was optimized on the nonmodulated CT. The evaluation follows by averaging the resulting n dose distributions and comparing to the corresponding nonmodulated dose distribution, attending at dosimetric indices and dose-volume histograms. In this work, the functionality of these routines was tested for proton and carbon ion plans for two selected lung cancer patient cases with deep-seated tumors, showing that, with existing standard tools, it is possible to calculate the beam-modulation effect of the lung in TRiP98 in a transparent way. Variable model parameters, such as modulation power, voxel size and density voxel selection range, were evaluated. Furthermore, a systematic study for spherical geometries in a lung tissue CT cube is presented to investigate general trends.

Purpose: The dose response of Gafchromic EBT3 films exposed to proton beams depends on the dose, and additionally on the beam quality, which is often quantified with the linear energy transfer (LET) and, hence, also referred to as LET quenching. Fundamentally different methods to determine correction factors for this LET quenching effect have been reported in literature and a new method using the local proton fluence distribution differential in LET is presented. This method was exploited to investigate whether a more practical correction based on the dose- or fluence-averaged LET is feasible in a variety of clinically possible beam arrangements.

Methods: The relative effectiveness (RE) was characterized within a high LET spread-out Bragg peak (SOBP) in water made up by the six lowest available energies (62.4-67.5 MeV, configuration b(1)) resulting in one of the highest clinically feasible dose-averaged LET distributions. Additionally, two beams were measured where a low LET proton beam (252.7 MeV) was superimposed on b1, which contributed either 50% of the initial particle fluence or 50% of the dose in the SOBP, referred to as configuration b(2) and b(3), respectively. The proton LET spectrum was simulated with GATE/Geant4 at all measurement positions. The net optical density change differential in LET was integrated over the local proton spectrum to calculate the net optical density and therefrom the beam quality correction factor. The LET dependence of the film response was accounted for by an LET dependence of one of the three parameters in the calibration function and was determined from inverse optimization using measurement b(1). This method was then validated on the measurements of b(2) and b(3) and subsequently used to calculate the RE at 900 positions in nine clinically relevant beams. The extrapolated RE set was used to derive a simple linear correction function based on dose-averaged LET (L-d) and verify the validity in all points of the comprehensive RE set.

Results: The uncorrected film dose deviated up to 26% from the reference dose, whereas the corrected film dose agreed within 3% in all three beams in water (b(1), b(2) and b(3)). The LET dependence of the calibration function started to strongly increase around 5 keV/lm and flatten out around 30 keV/mu m. All REs calculated from the proton fluence in the nine simulated beams could be approximated with a linear function of dose-averaged LET (RE = 1.0258-0.0211 lm/keV L-d). However, no functional relationship of RE- and fluence-averaged LET could be found encompassing all beam energies and modulations.

Conclusions: The film quenching was found to be nonlinear as a function of proton LET as well as of the dose-averaged LET. However, the linear relation of RE on dose-averaged LET was a good approximation in all cases. In contrast to dose-averaged LET, fluence-averaged LET could not describe the RE when multiple beams were applied.

The ion recombination is examined in parallel-plate ionization chambers in scanning proton beams at the Danish Centre for Particle Therapy and the Skandion Clinic. The recombination correction factor k(s) is investigated for clinically relevant energies between 70 MeV and 244 MeV for dose rates below 400 Gy min(-1) in air. The Boutillon formalism is used to separate the initial and general recombination. The general recombination is compared to predictions from the numerical recombination code IonTracks and the initial recombination to the Jaffe theory. k(s) is furthermore calculated with the two-voltage method (TVM) and extrapolation approaches, in particular the recently proposed three-voltage (3VL) method. The TVM is in agreement with the Boutillon method and IonTracks for dose rates above 100 Gy min(-1). However, the TVM calculated k(s) is closer related to the Jaffe theory for initial recombination for lower dose rate, indicating a limited application in scanning light ion beams. The 3VL is in turn found to generally be in agreement with Boutillon's method. The recombination is mapped as a function of the dose rate and proton energy at the two centres using the Boutillon formalism: the initial recombination parameter was found to be A = (0.10 +/- 0.01) V at DCPT and A = (0.22 +/- 0.13) V at Skandion, which is in better agreement with the Jaffe theory for initial recombination than previously reported values. The general recombination parameter was estimated to m2=(4.7 +/- 0.1).103V2nA-1cm-1m2=(7.2 +/- 0.1).103V2nA-1cm-1<i. Furthermore, the numerical algorithm IonTracks is demonstrated to correctly predict the initial recombination at low dose rates and the general recombination at high dose rates.

Purpose: Conventional X-ray radiotherapy induces a pro-inflammatory response mediated by altered expression of inflammation-regulating cytokines. Proton scanning and X-ray irradiation produce distinct changes to cytokine gene expression in vitro suggesting that proton beam therapy may induce an inflammatory response dissimilar to that of X-ray radiation. The purpose of the present study was to determine whether proton scanning beam radiation and conventional X-ray photon radiation would induce differential regulation of circulating cytokines in vivo.

Materials and methods: Female CDF1 mice were irradiated locally at the right hind leg using proton pencil beam scanning or X-ray photons. Blood samples were obtained from two separate mice groups. Samples from one group were drawn by retro-orbital puncture 16 months post irradiation, while samples from the other group were drawn 5 and 30 days post irradiation. Concentration of the cytokines IL-6, IL-1 beta, IL-10, IL-17A, IFN-gamma, and TNF alpha was measured in plasma using bead-based immunoassays.

Results: The cytokines IL-6, IL-1 beta, IL-10, IFN-gamma, and TNF alpha were expressed at lower levels in plasma samples from proton-irradiated mice compared with X-ray-irradiated mice 16 months post irradiation. The same cytokines were downregulated in proton-irradiated mice 5 days post irradiation when compared to controls, while at day 30 expression had increased to the same level or higher. X-ray radiation did not markedly change expression levels at days 5 and 30.

Conclusions: The inflammatory response to proton and X-ray irradiation seem to be distinct as the principal pro-inflammatory cytokines are differentially regulated short- and long-term following irradiation. Both the development of normal tissue damage and efficacy of immunotherapy could be influenced by an altered inflammatory response to irradiation.

Purpose: To identify differential cellular responses after proton and photon irradiation by comparing transcriptomes of primary fibroblasts irradiated with either radiation type. Methods and Materials: A panel of primary dermal fibroblast cultures was irradiated with low and higher linear energy transfer (LET) proton beams. Cobalt-60 photon irradiation was used as reference. Dose was delivered in 3 fractions of 3.5 Gy (relative biological effectiveness) using a relative biological effectiveness of 1.1 for proton doses. Cells were harvested 2 hours after the final fraction was delivered, and RNA was purified. RNA sequencing was performed using Illumina NextSeq 500 with high-output kit. The edgeR package in R was used for differential gene expression analysis. Results: Pairwise comparisons of the transcriptomes in the 3 treatment groups showed that there were 84 and 56 differentially expressed genes in the low LET group compared with the Cobalt-60 group and the higher LET group, respectively. The higher LET proton group and the Cobalt-60 group had the most distinct transcriptome profiles, with 725 differentially regulated genes. Differentially regulated canonical pathways and various regulatory factors involved in regulation of biological mechanisms such as inflammation, carcinogenesis, and cell cycle control were identified. Conclusions: Inflammatory regulators associated with the development of normal tissue complications and malignant transformation factors seem to be differentially regulated by higher LET proton and Cobalt-60 photon irradiation. The reported transcriptome differences could therefore influence the progression of adverse effects and the risk of developing secondary cancers.

Ionization quenching in ion beam dosimetry is often related to the fluence- or dose-averaged linear energy transfer (LET). Both quantities are however averaged over a wide LET range and a mixed field of primary and secondary ions. We propose a novel method to correct the quenched luminescence in scintillators exposed to ion beams. The method uses the energy spectrum of the primaries and accounts for the varying quenched luminescence in heavy, secondary ion tracks through amorphous track structure theory. The new method is assessed against more traditional approaches by correcting the quenched luminescence response from the BCF-12, BCF-60, and 81-0084 plastic scintillators exposed to a 100 MeV pristine proton beam in order to compare the effects of the averaged LET quantities and the secondary ions. Calculations and measurements show that primary protons constitute more than 92% of the energy deposition but account for more than 95% of the luminescence signal in the scintillators. The quenching corrected luminescence signal is in better agreement with the dose measurement when the secondary particles are taken into account. The Birks model provided the overall best quenching corrections, when the quenching corrected signal is adjusted for the number of free model parameters. The quenching parameter kB for the BCF-12 and BCF-60 scintillators is in agreement with literature values and was found to be kB = (10.6 +/- 0.1) x 10(-2) mu m keV(-1) for the 81-0084 scintillator. Finally, a fluence threshold for the 100 MeV proton beam was calculated to be of the order of 10(10) cm(-2), corresponding to 110 Gy, above which the quenching increases non-linearly and the Birks model no longer is applicable.

Radiotherapy beams of protons or heavier ions generate secondary particles through nuclear interactions over different patient tissues. The resulting particle spectra depend on the tissue composition and on charge and energy of the primary beam ions. In proton radiotherapy, predictive radiobiological models usually apply dose-averaged linear energy transfer (LET). Microdosimetry-based models for proton or heavier ion primary beams also rely on dose-averaged quantities, the values of which depend on whether the produced secondaries are included or excluded in the calculation. In turn, this will affect the results of calculations of the relative biological effectiveness (RBE) of these beams. In this brief note, we study quantitatively the influence of the secondary radiation spectra on the averaged expectation values of LET and their impact on predictions of RBE. It is noted that for microdosimetry-based quantities and for corresponding LET-based parameters the trends are similar and that fluence-averaged quantities should be studied more closely.

Gafchromic EBT3 films are applied in proton radiotherapy for 2D dose mapping because they demonstrate spatial resolution well below 1mm. However, the film response must be corrected in order to reach the accuracy of dose measurements required for the clinical use. The in-house developed AnalyseGafchromic software allows to analyze and correct the measured response using triple channel dose calibration, statistical scan-to-scan fluctuations as well as experimentally determined dose and LET dependence. Finally, the optimized protocol for evaluation of response of Gafchromic EBT3 films was applied to determine 30 x 40 cm(2) dose profiles of the scanning therapy unit at the Cyclotron Centre Bronowice, CCB in Krakow, Poland.

In this chapter, we describe two different methods, analytical (pencil beam) algorithms and Monte Carlo simulations, used to obtain the intended dose distributions in patients and evaluate their strengths and shortcomings. We discuss the difference between the prescribed physical dose and the biologically effective dose, the relative biological effectiveness (RBE) between ions and photons and the dependence of RBE on the linear energy transfer (LET). Lastly, we show how LET- or RBE-based optimization can be used to improve treatment plans and explore how the availability of multimodality ion beam facilities can be used to design a tumor-specific optimal treatment.