Performance of LGAD strip detectors for particle counting of therapeutic proton beams.

Objective. The performance of silicon detectors with moderate internal gain, named low-gain avalanche diodes (LGADs), was studied to investigate their capability to discriminate and count single beam particles at high fluxes, in view of future applications for beam characterization and on-line beam monitoring in proton therapy.Approach. Dedicated LGAD detectors with an active thickness of 55µm and segmented in 2 mm2strips were characterized at two Italian proton-therapy facilities, CNAO in Pavia and the Proton Therapy Center of Trento, with proton beams provided by a synchrotron and a cyclotron, respectively. Signals from single beam particles were discriminated against a threshold and counted. The number of proton pulses for fixed energies and different particle fluxes was compared with the charge collected by a compact ionization chamber, to infer the input particle rates.Main results. The counting inefficiency due to the overlap of nearby signals was less than 1% up to particle rates in one strip of 1 MHz, corresponding to a mean fluence rate on the strip of about 5 × 107p/(cm2·s). Count-loss correction algorithms based on the logic combination of signals from two neighboring strips allow to extend the maximum counting rate by one order of magnitude. The same algorithms give additional information on the fine time structure of the beam.Significance. The direct counting of the number of beam protons with segmented silicon detectors allows to overcome some limitations of gas detectors typically employed for beam characterization and beam monitoring in particle therapy, providing faster response times, higher sensitivity, and independence of the counts from the particle energy.

Calibration method and performance of a time-of-flight detector to measure absolute beam energy in proton therapy.

BACKGROUND: The beam energy is one of the most significant parameters in particle therapy since it is directly correlated to the particles' penetration depth inside the patient. Nowadays, the range accuracy is guaranteed by offline routine quality control checks mainly performed with water phantoms, 2D detectors with PMMA wedges, or multi-layer ionization chambers. The latter feature low sensitivity, slow collection time, and response dependent on external parameters, which represent limiting factors for the quality controls of beams delivered with fast energy switching modalities, as foreseen in future treatments. In this context, a device based on solid-state detectors technology, able to perform a direct and absolute beam energy measurement, is proposed as a viable alternative for quality assurance measurements and beam commissioning, paving the way for online range monitoring and treatment verification. PURPOSE: This work follows the proof of concept of an energy monitoring system for clinical proton beams, based on Ultra Fast Silicon Detectors (featuring tenths of ps time resolution in 50 µm active thickness, and single particle detection capability) and time-of-flight techniques. An upgrade of such a system is presented here, together with the description of a dedicated self-calibration method, proving that this second prototype is able to assess the mean particles energy of a monoenergetic beam without any constraint on the beam temporal structure, neither any a priori knowledge of the beam energy for the calibration of the system. METHODS: A new detector geometry, consisting of sensors segmented in strips, has been designed and implemented in order to enhance the statistics of coincident protons, thus improving the accuracy of the measured time differences. The prototype was tested on the cyclotron proton beam of the Trento Protontherapy Center (TPC). In addition, a dedicated self-calibration method, exploiting the measurement of monoenergetic beams crossing the two telescope sensors for different flight distances, was introduced to remove the systematic uncertainties independently from any external reference. RESULTS: The novel calibration strategy was applied to the experimental data collected at TPC (Trento) and CNAO (Pavia). Deviations between measured and reference beam energies in the order of a few hundreds of keV with a maximum uncertainty of 0.5 MeV were found, in compliance with the clinically required water range accuracy of 1 mm. CONCLUSIONS: The presented version of the telescope system, minimally perturbative of the beam, relies on a few seconds of acquisition time to achieve the required clinical accuracy and therefore represents a feasible solution for beam commission, quality assurance checks, and online beam energy monitoring.

The INFN proton computed tomography system for relative stopping power measurements: calibration and verification.

Objective. This paper describes the procedure to calibrate the three-dimensional (3D) proton stopping power relative to water (SPR) maps measured by the proton computed tomography (pCT) apparatus of the Istituto Nazionale di Fisica Nucleare (INFN, Italy). Measurements performed on water phantoms are used to validate the method. The calibration allowed for achieving measurement accuracy and reproducibility to levels below 1%.Approach. The INFN pCT system is made of a silicon tracker for proton trajectory determination followed by a YAG:Ce calorimeter for energy measurement. To perform the calibration, the apparatus has been exposed to protons of energies ranging from 83 to 210 MeV. Using the tracker, a position-dependent calibration has been implemented to keep the energy response uniform across the calorimeter. Moreover, correction algorithms have been developed to reconstruct the proton energy when this is shared in more than one crystal and to consider the energy loss in the non-uniform apparatus material. To verify the calibration and its reproducibility, water phantoms have been imaged with the pCT system during two data-taking sessions.Main results. The energy resolution of the pCT calorimeter resulted to beσEE≅0.9%at 196.5 MeV. The average values of the water SPR in fiducial volumes of the control phantoms have been calculated to be 0.995±0.002. The image non-uniformities were below 1%. No appreciable variation of the SPR and uniformity values between the two data-taking sessions could be identified.Significance. This work demonstrates the accuracy and reproducibility of the calibration of the INFN pCT system at a level below 1%. Moreover, the uniformity of the energy response keeps the image artifacts at a low level even in the presence of calorimeter segmentation and tracker material non-uniformities. The implemented calibration technique allows the INFN-pCT system to face applications where the precision of the SPR 3D maps is of paramount importance.

Characterization of the INFN proton CT scanner for cross-calibration of x-ray CT.

Objective. The goal of this study was to assess the imaging performances of the pCT system developed in the framework of INFN-funded (Italian National Institute of Nuclear Physics) research projects. The spatial resolution, noise power spectrum (NPS) and RSP accuracy has been investigated, as a preliminary step to implement a new cross-calibration method for x-ray CT (xCT).Approach. The INFN pCT apparatus, made of four planes of silicon micro-strip detectors and a YAG:Ce scintillating calorimeter, reconstructs 3D RSP maps by a filtered-back projection algorithm. The imaging performances (i.e. spatial resolution, NPS and RSP accuracy) of the pCT system were assessed on a custom-made phantom, made of plastic materials with different densities ((0.66, 2.18) g cm-3). For comparison, the same phantom was acquired with a clinical xCT system.Main results. The spatial resolution analysis revealed the nonlinearity of the imaging system, showing different imaging responses in air or water phantom background. Applying the Hann filter in the pCT reconstruction, it was possible to investigate the imaging potential of the system. Matching the spatial resolution value of the xCT (0.54 lp mm-1) and acquiring both with the same dose level (11.6 mGy), the pCT appeared to be less noisy than xCT, with an RSP standard deviation of 0.0063. Concerning the RSP accuracy, the measured mean absolute percentage errors were (0.23+-0.09)% in air and (0.21+-0.07)% in water.Significance. The obtained performances confirm that the INFN pCT system provides a very accurate RSP estimation, appearing to be a feasible clinical tool for verification and correction of xCT calibration in proton treatment planning.

LIDAL, a Time-of-Flight Radiation Detector for the International Space Station: Description and Ground Calibration.

LIDAL (Light Ion Detector for ALTEA, Anomalous Long-Term Effects on Astronauts) is a radiation detector designed to measure the flux, the energy spectra and, for the first time, the time-of-flight of ions in a space habitat. It features a combination of striped silicon sensors for the measurement of deposited energy (using the ALTEA device, which operated from 2006 to 2012 in the International Space Station) and fast scintillators for the time-of-flight measurement. LIDAL was tested and calibrated using the proton beam line at TIFPA (Trento Institute for Fundamental Physics Application) and the carbon beam line at CNAO (National Center for Oncology Hadron-therapy) in 2019. The performance of the time-of-flight system featured a time resolution (sigma) less than 100 ps. Here, we describe the detector and the results of these tests, providing ground calibration curves along with the methodology established for processing the detector's data. LIDAL was uploaded in the International Space Station in November 2019 and it has been operative in the Columbus module since January 2020.

Does variable RBE affect toxicity risks for mediastinal lymphoma patients? NTCP-based evaluation after proton therapy treatment.

INTRODUCTION: Mediastinal lymphoma (ML) is a solid malignancy affecting young patients. Modern combined treatments allow obtaining good survival probability, together with a long life expectancy, and therefore with the need to minimize treatment-related toxicities. We quantified the expected toxicity risk for different organs and endpoints in ML patients treated with intensity-modulated proton therapy (IMPT) at our centre, accounting also for uncertainties related to variable RBE. METHODS: Treatment plans for ten ML patients were recalculated with a TOPAS-based Monte Carlo code, thus retrieving information on LET and allowing the estimation of variable RBE. Published NTCP models were adopted to calculate the toxicity risk for hypothyroidism, heart valve defects, coronary heart disease and lung fibrosis. NTCP was calculated assuming both constant (i.e. 1.1) and variable RBE. The uncertainty associated with individual radiosensitivity was estimated by random sampling α/ß values before RBE evaluation. RESULTS: Variable RBE had a minor impact on hypothyroidism risk for 7 patients, while it led to significant increase for the remaining three (+24% risk maximum increase). Lung fibrosis was slightly affected by variable RBE, with a maximum increase of â‰… 1%. This was similar for heart valve dysfunction, with the exception of one patient showing an about 10% risk increase, which could be explained by means of large heart volume and D1 increase. DISCUSSION: The use of NTCP models allows for identifying those patients associated with a higher toxicity risk. For those patients, it might be worth including variable RBE in plan evaluation.

Investigation of In-Field and Out-of-Field Radiation Quality With Microdosimetry and Its Impact on Relative Biological Effectiveness in Proton Therapy.

PURPOSE: Using microdosimetry, this study investigated the relative biological effectiveness (RBE) and quality factor (Q¯) variations in field and out of field as a function of radiation quality for clinical protons. METHODS AND MATERIALS: A water phantom with a spread-out Bragg peak (SOBP) was irradiated to acquire microdosimetric spectra at several distal and lateral depths with a tissue equivalent proportional counter. The measurements were used as inputs to microdosimetric kinetic and Loncol models to determine the RBE spatial distribution and compare it with predictions from the dose-averaged linear energy transfer-based McNamara model. Q¯ values and biological and dose equivalent values were also calculated. RESULTS: The data demonstrated that radiation quality changed more rapidly with depth than lateral distance from the SOBP. In beam, yD ranged from approximately 4 keV/µm at the entrance to 8 keV/µm at the SOBP far end, reaching approximately 15 keV/µm at the penumbra. Out of field, the overall highest value of 23 ± 2 keV/µm was observed at the beam-edge penumbra. Radiation quality changes caused RBE deviations from the clinical value of 1.1, whose extent depends on the approach used for assessing radiation quality as well as on the radiobiological model. For RBE10, microdosimetry-based models appeared to better reproduce the radiobiological data than the dose-averaged linear energy transfer model. Out of field, both the RBE and Q¯ values appeared to have limitations in describing the radiation biological effectiveness. This research also presents a first comprehensive benchmark of TOPAS code against in-field and out-of-field microdosimetric spectra of therapeutic protons. CONCLUSIONS: Further investigation will be necessary to evaluate the quantitative effects of RBE variations on treatment planning and assess the clinical consequences in terms of both tumor control and normal-tissue toxicity. The achievement of this goal calls for accurate radiobiological data to validate the RBE models.

Modelling the HPRT-gene mutation induction of particle beams: systematicin vitrodata collection, analysis and microdosimetric kinetic model implementation.

Objective. Since the early years, particle therapy treatments have been associated with concerns for late toxicities, especially secondary cancer risk (SCR). Nowadays, this concern is related to patients for whom long-term survival is expected (e.g. breast cancer, lymphoma, paediatrics). In the aim to contribute to this research, we present a dedicated statistical and modelling analysis aiming at improving our understanding of the RBE for mutation induction (RBEM˜) for different particle species.Approach. We built a new database based on a systematic collection of RBE data for mutation assays of the gene encoding for the purine salvage enzyme hypoxanthine-guanine phosphoribosyltransferase from literature (105 entries, distributed among 3 cell lines and 16 particle species). The data were employed to perform statistical and modelling analysis. For the latter, we adapted the microdosimetric kinetic model (MKM) to describe the mutagenesis in analogy to lethal lesion induction.Main results. Correlation analysis between RBE for survival (RBES) andRBEM˜reveals significant correlation between these two quantities (ρ= 0.86,p< 0.05). The correlation gets stronger when looking at subsets of data based on cell line and particle species. We also show that the MKM can be successfully employed to describeRBEM˜,obtaining comparably good agreement with the experimental data. Remarkably, to improve the agreement with experimental data the MKM requires, consistently in all the analysed cases, a reduced domain size for the description of mutation induction compared to that adopted for survival.Significance. We were able to show that RBESandRBEM˜are strongly related quantities. We also showed for the first time that the MKM could be successfully applied to the description of mutation induction, representing an endpoint different from the more traditional cell killing. In analogy to the RBES,RBEM˜can be implemented into treatment planning system evaluations.

Diagnosis of Langerhans cell histiocytosis on cytological examination of cerebrospinal fluid: Report of the first case.

Langerhans cell histiocytosis (LCH) is a disease of unknown etiology characterized by a proliferation of histiocytic cells resembling dendritic Langerhans cells. LCH can be unifocal or multifocal, with one- or many-organ involvement. The serous fluids are rarely involved. Cytological diagnosis of LCH is possible and relies on recognition of the typical cytomorphological features and subsequent immunocytochemical confirmation. Given the possibility of multisystem involvement, after diagnosing LCH it is necessary to carry out staging exams such as a bone survey, abdominal ultrasound, complete blood count, screening for diabetes insipidus and pulmonary function tests. We present the first case of LCH where the diagnosis was reached on cytological material from the cerebrospinal fluid. To the best of our knowledge, this is the first such case reported in the international literature to date. The morphological and immunocytochemical characteristics of our case are described, and the relevant literature is reviewed.

Quantification of biological range uncertainties in patients treated at the Krakow proton therapy centre.

BACKGROUND: Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation. METHODS AND MATERIALS: The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach. RESULTS: The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm. DISCUSSION: Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area. CONCLUSIONS: Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

Biological Impact of Target Fragments on Proton Treatment Plans: An Analysis Based on the Current Cross-Section Data and a Full Mixed Field Approach.

Clinical routine in proton therapy currently neglects the radiobiological impact of nuclear target fragments generated by proton beams. This is partially due to the difficult characterization of the irradiation field. The detection of low energetic fragments, secondary protons and fragments, is in fact challenging due to their very short range. However, considering their low residual energy and therefore high LET, the possible contribution of such heavy particles to the overall biological effect could be not negligible. In this context, we performed a systematic analysis aimed at an explicit assessment of the RBE (relative biological effectiveness, i.e., the ratio of photon to proton physical dose needed to achieve the same biological effect) contribution of target fragments in the biological dose calculations of proton fields. The TOPAS Monte Carlo code has been used to characterize the radiation field, i.e., for the scoring of primary protons and fragments in an exemplary water target. TRiP98, in combination with LEM IV RBE tables, was then employed to evaluate the RBE with a mixed field approach accounting for fragments' contributions. The results were compared with that obtained by considering only primary protons for the pristine beam and spread out Bragg peak (SOBP) irradiations, in order to estimate the relative weight of target fragments to the overall RBE. A sensitivity analysis of the secondary particles production cross-sections to the biological dose has been also carried out in this study. Finally, our modeling approach was applied to the analysis of a selection of cell survival and RBE data extracted from published in vitro studies. Our results indicate that, for high energy proton beams, the main contribution to the biological effect due to the secondary particles can be attributed to secondary protons, while the contribution of heavier fragments is mainly due to helium. The impact of target fragments on the biological dose is maximized in the entrance channels and for small α/ß values. When applied to the description of survival data, model predictions including all fragments allowed better agreement to experimental data at high energies, while a minor effect was observed in the peak region. An improved description was also obtained when including the fragments' contribution to describe RBE data. Overall, this analysis indicates that a minor contribution can be expected to the overall RBE resulting from target fragments. However, considering the fragmentation effects can improve the agreement with experimental data for high energy proton beams.

Study of relationship between dose, LET and the risk of brain necrosis after proton therapy for skull base tumors.

PURPOSE: We investigated the relationship between RBE-weighted dose (DRBE) calculated with constant (cRBE) and variable RBE (vRBE), dose-averaged linear energy transfer (LETd) and the risk of radiographic changes in skull base patients treated with protons. METHODS: Clinical treatment plans of 45 patients were recalculated with Monte Carlo tool FRED. Radiographic changes (i.e. edema and/or necrosis) were identified by MRI. Dosimetric parameters for cRBE and vRBE were computed. Biological margin extension and voxel-based analysis were employed looking for association of DRBE(vRBE) and LETd with brain edema and/or necrosis. RESULTS: When using vRBE, Dmax in the brain was above the highest dose limits for 38% of patients, while such limit was never exceeded assuming cRBE. Similar values of Dmax were observed in necrotic regions, brain and temporal lobes. Most of the brain necrosis was in proximity to the PTV. The voxel-based analysis did not show evidence of an association with high LETd values. CONCLUSIONS: When looking at standard dosimetric parameters, the higher dose associated with vRBE seems to be responsible for an enhanced risk of radiographic changes. However, as revealed by a voxel-based analysis, the large inter-patient variability hinders the identification of a clear effect for high LETd.

Technical Note: CT calibration for proton treatment planning by cross-calibration with proton CT data.

PURPOSE: This study explores the possibility of a new method for x-ray computed tomography (CT) calibration by means of cross-calibration with proton CT (pCT) data. The proposed method aims at a more accurate conversion of CT Hounsfield Units (HU) into proton stopping power ratio (SPR) relative to water to be used in proton-therapy treatment planning. METHODS: X-ray CT scan was acquired on a synthetic anthropomorphic phantom, composed of different tissue equivalent materials (TEMs). A pCT apparatus was instead adopted to obtain a reference three-dimensional distribution of the phantom's SPR values. After rigid registration, the x-ray CT was artificially blurred to the same resolution of pCT. Then a scatter plot showing voxel-by-voxel SPR values as a function of HU was employed to link the two measurements and thus obtaining a cross-calibrated x-ray CT calibration curve. The cross-calibration was tested at treatment planning system and then compared with a conventional calibration based on exactly the same TEMs constituting the anthropomorphic phantom. RESULTS: Cross-calibration provided an accurate SPR mapping, better than by conventional TEMs calibration. The dose distribution of single beams optimized on the reference SPR map was recomputed on cross-calibrated CT, showing, with respect to conventional calibration, minor deviation at the dose fall-off (lower than 1%). CONCLUSIONS: The presented data demonstrated that, by means of reference pCT data, a heterogeneous phantom can be used for CT calibration, paving the way to the use of biological samples, with their accurate description of patients' tissues. This overcomes the limitations of conventional CT calibration requiring homogenous samples, only available by synthetic TEMs, which fail in accurately mimicking the properties of biological tissues. Once a heterogeneous biological sample is provided with its corresponding reference SPR maps, a cross-calibration procedure could be adopted by other PT centers, even when not equipped with a pCT system.

Relative stopping power measurements and prosthesis artifacts reduction in proton CT.

We present a set-up for proton computed tomography (pCT), composed of a microstrip silicon tracker and a YAG:Ce calorimeter, able to directly measure the relative stopping power (RSP) maps to be used in hadron therapy. The system, tested with an electron density phantom at the Trento proton Therapy Center, is able to correlate measured and expected RSP with discrepancies of the order of 1% or less. Furthermore, pCT tomographies of an anthropomorphous head phantom taken with our device, when compared with x-ray CT images of the same object, evidence a significant reduction of artifacts induced by titanium spinal bone prosthesis and tungsten dental filling.

Proton pencil beam scanning reduces secondary cancer risk in breast cancer patients with internal mammary chain involvement compared to photon radiotherapy.

PURPOSE: Proton pencil beam scanning (PBS) represents an interesting option for the treatment of breast cancer (BC) patients with nodal involvement. Here we compare tangential 3D-CRT and VMAT to PBS proton therapy (PT) in terms of secondary cancer risk (SCR) for the lungs and for contralateral breast. METHODS: Five BC patients including supraclavicular (SVC) nodes in the target (Group 1) and five including SVC plus internal-mammary-nodes (IMNs, Group 2) were considered. The Group 1 patients were planned by PT versus tangential 3D-CRT in free-breathing (FB). The Group 2 patients were planned by PT versus VMAT considering both FB and deep-inspiration breath hold (DIBH) irradiation. The prescription dose to the target volume was 50 Gy (2 Gy/fraction). A constant RBE = 1.1 was assumed for PT. The SCR was evaluated with the excess absolute risk (EAR) formalism, considering also the age dependence. A cumulative EAR was finally computed. RESULTS: According to the linear, linear-exponential and linear-plateau dose response model, the cumulative EAR for Group 1 patients after PT was equal to 45 ± 10, 17 ± 3 and 15 ± 3, respectively. The corresponding relative increase for tangential 3D-CRT was equal to a factor 2.1 ± 0.5, 2.1 ± 0.4 and 2.3 ± 0.4. Group 2 patients showed a cumulative EAR after PT in FB equal to 65 ± 3, 21 ± 1 and 20 ± 1, according to the different models; the relative risk obtained with VMAT increased by a factor 3.5 ± 0.2, 5.2 ± 0.3 and 5.1 ± 0.3. Similar values emerge from DIBH plans. Contrary to photon radiotherapy, PT appears to be not sensitive to the age dependence due to the very low delivered dose. CONCLUSIONS: PBS PT is associated to significant SCR reduction in BC patients compared to photon radiotherapy. The benefits are maximized for young patients with both SVC and IMNs involvement. When combined with the improved sparing of the heart, this might contribute to the establishment of effective patient-selection criteria for proton BC treatments.

Proton or photon radiosurgery for cardiac ablation of ventricular tachycardia? Breath and ECG gated robust optimization.

PURPOSE: Ventricular tachycardia (VT) is a life-threatening heart disorder. The aim of this preliminary study is to assess the feasibility of stereotactic body radiation therapy (SBRT) photon and proton therapy (PT) plans for the treatment of VT, adopting robust optimization technique for both irradiation techniques. METHODS: ECG gated CT images (in breath hold) were acquired for one patient. Conventional planning target volume (PTV) and robust optimized plans (25GyE in single fraction) were simulated for both photon (IMRT, 5 and 9 beams) and proton (SFO, 2 beams) plans. Robust optimized plans were obtained both for protons and photons considering in the optimization setup errors (5 mm in the three orthogonal directions), range (±3.5%) and the clinical target volume (CTV) motion due to heartbeat and breath-hold variability. RESULTS: The photon robust optimization method, compared to PTV-based plans, showed a reduction in the average dose to the heart by about 25%; robust optimization allowed also reducing the mean dose to the left lung from 3.4. to 2.8 Gy for 9-beams configuration and from 4.1 to 2.9 Gy for 5-beams configuration. Robust optimization with protons, allowed further reducing the OAR doses: average dose to the heart and to the left lung decreased from 7.3 Gy to 5.2 GyE and from 2.9 Gy to 2.2 GyE, respectively. CONCLUSIONS: Our study demonstrates the importance of the optimization technique adopted in the treatment planning system for VT treatment. It has been shown that robust optimization can significantly reduce the dose to healthy cardiac tissues and that PT further increases this gain.

Evaluation of proton beam radiation-induced skin injury in a murine model using a clinical SOBP.

The purpose of this paper is to characterize the skin deterministic damage due to the effect of proton beam irradiation in mice occurred during a long-term observational experiment. This study was initially defined to evaluate the insurgence of myelopathy irradiating spinal cords with the distal part of a Spread-out Bragg peak (SOBP). To the best of our knowledge, no study has been conducted highlighting high grades of skin injury at the dose used in this paper. Nevertheless these effects occurred. In this regard, the experimental evidence of significant insurgence of skin injury induced by protons using a SOBP configuration will be shown. Skin damages were classified into six scores (from 0 to 5) according to the severity of the injuries and correlated to ED50 (i.e. the radiation dose at which 50% of animals show a specific score) at 40 days post-irradiation (d.p.i.). The effects of radiation on the overall animal wellbeing have been also monitored and the severity of radiation-induced skin injuries was observed and quantified up to 40 d.p.i.

Accurate proton treatment planning for pencil beam crossing titanium fixation implants.

PURPOSE: To present a planning strategy for proton pencil-beam scanning when titanium implants need to be crossed by the beam. METHODS: We addressed three issues: the implementation of a CT calibration curve to assign to titanium the correct stopping power; the effect of artefacts on CT images and their reduction by a dedicated algorithm; the differences in dose computation depending on the dose engine, pencil-beam vs Monte-Carlo algorithms. We performed measurement tests on a simple cylinder phantom and on a real implant. These phantoms were irradiated with three geometries (single spots, uniform mono-energetic layer and uniform box), measuring the exit dose either by radio-chromic film or multi-layer ionization chamber. The procedure was then applied on two patients treated for chordoma. RESULTS: We had to set in the calibration curve a mass density equal to 4.37 g/cm3 to saturated Hounsfield Units, in order to have the correct stopping power assigned to titanium in TPS. CT artefact reduction algorithm allowed a better reconstruction of the shape and size of the implant. Monte-Carlo resulted accurate in computing the dose distribution whereas the pencil-beam algorithm failed due to sharp density interfaces between titanium and the surrounding material. Finally, the treatment plans obtained on two patients showed the impact of the dose engine algorithm, with 10-20% differences between pencil-beam and Monte-Carlo in small regions distally to the titanium screws. CONCLUSION: The described combination of CT calibration, artefacts reduction and Monte-Carlo computation provides a reliable methodology to compute dose in patients with titanium implants.

Modelling the risk of radiation induced alopecia in brain tumor patients treated with scanned proton beams.

PURPOSE: To develop normal tissue complication probability (NTCP) models for radiation-induced alopecia (RIA) in brain tumor patients treated with proton therapy (PT). METHODS AND MATERIALS: We analyzed 116 brain tumor adult patients undergoing scanning beam PT (median dose 54 GyRBE; range 36-72) for CTCAE v.4 grade 2 (G2) acute (≤90 days), late (>90 days) and permanent (>12 months) RIA. The relative dose-surface histogram (DSH) of the scalp was extracted and used for Lyman-Kutcher-Burman (LKB) modelling. Moreover, DSH metrics (Sx: the surface receiving ≥ X Gy, D2%: near maximum dose, Dmean: mean dose) and non-dosimetric variables were included in a multivariable logistic regression NTCP model. Model performances were evaluated by the cross-validated area under the receiver operator curve (ROC-AUC). RESULTS: Acute, late and permanent G2-RIA was observed in 52%, 35% and 19% of the patients, respectively. The LKB models showed a weak dose-surface effect (0.09 ≤ n ≤ 0.19) with relative steepness 0.29 ≤ m ≤ 0.56, and increasing tolerance dose values when moving from acute and late (22 and 24 GyRBE) to permanent RIA (44 GyRBE). Multivariable modelling selected S21Gy for acute and S25Gy, for late G2-RIA as the most predictive DSH factors. Younger age was selected as risk factor for acute G2-RIA while surgery as risk factor for late G2-RIA. D2% was the only variable selected for permanent G2-RIA. Both LKB and logistic models exhibited high predictive performances (ROC-AUCs range 0.86-0.90). CONCLUSION: We derived NTCP models to predict G2-RIA after PT, providing a comprehensive modelling framework for acute, late and permanent occurrences that, once externally validated, could be exploited for individualized scalp sparing treatment planning strategies in brain tumor patients.

Clinical implementation in proton therapy of multi-field optimization by a hybrid method combining conventional PTV with robust optimization.

To implement a robust multi-field optimization (MFO) technique compatible with the application of a Monte Carlo (MC) algorithm and to evaluate its robustness. Nine patients (three brain, five head-and-neck, one spine) underwent proton treatment generated by a novel robust MFO technique. A hybrid (hMFO) approach was implemented, planning dose coverage on isotropic PTV compensating for setup errors, whereas range calibration uncertainties are incorporated into PTV robust optimization process. hMFO was compared with single-field optimization (SFO) and full robust multi-field optimization (fMFO), both on the nominal plan and the worst-case scenarios assessed by robustness analysis. The SFO and the fMFO plans were normalized to hMFO on CTV to obtain iso-D95 coverage, and then the organs at risk (OARs) doses were compared. On the same OARs, in the normalized nominal plans the potential impact of variable relative biological effectiveness (RBE) was investigated. hMFO reduces the number of scenarios computed for robust optimization (from twenty-one in fMFO to three), making it practicable with the application of a MC algorithm. After normalizing on D95 CTV coverage, nominal hMFO plans were superior compared to SFO in terms of OARs sparing (p  < 0.01), without significant differences compared to fMFO. The improvement in OAR sparing with hMFO with respect to SFO was preserved in worst-case scenarios (p  < 0.01), confirming that hMFO is as robust as SFO to physical uncertainties, with no significant differences when compared to the worst case scenarios obtained by fMFO. The dose increase on OARs due to variable RBE was comparable to the increase due to physical uncertainties (i.e. 4-5 Gy(RBE)), but without significant differences between these techniques. hMFO allows improving plan quality with respect to SFO, with no significant differences with fMFO and without affecting robustness to setup, range and RBE uncertainties, making clinically feasible the application of MC-based robust optimization.