Dose averaged linear energy transfer optimization for large sacral chordomas in carbon ion therapy.

BACKGROUND: Carbon ion beams are well accepted as densely ionizing radiation with a high linear energy transfer (LET). However, the current clinical practice does not fully exploit the highest possible dose-averaged LET (LETd) and, consequently, the biological potential in the target. This aspect becomes worse in larger tumors for which inferior clinical outcomes and corresponding lower LETd was reported. PURPOSE: The vicinity to critical organs in general and the inferior overall survival reported for larger sacral chordomas treated with carbon ion radiotherapy (CIRT), makes the treatment of such tumors challenging. In this work it was aimed to increase the LETd in large volume tumors while maintaining the relative biological effectiveness (RBE)-weighted dose, utilizing the LETd optimization functions of a commercial treatment planning system (TPS). METHODS: Ten reference sequential boost carbon ion treatment plans, designed to mimic clinical plans for large sacral chordoma tumors, were generated. High dose clinical target volumes (CTV-HD) larger than 250 cm 3 $250 \,{\rm cm}^{3}$ were considered as large targets. The total RBE-weighted median dose prescription with the local effect model (LEM) was D RBE , 50 % = 73.6 Gy $\textrm {D}_{\rm RBE, 50\%}=73.6 \,{\rm Gy}$ in 16 fractions (nine to low dose and seven to high dose planning target volume). No LETd optimization was performed in the reference plans, while LETd optimized plans used the minimum LETd (Lmin) optimization function in RayStation 2023B. Three different Lmin values were investigated and specified for the seven boost fractions: L min = 60 keV / µ m $\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ , L min = 80 keV / µ m $\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ and L min = 100 keV / µ m $\textrm {L}_{\rm min}=100 \,{\rm keV}/{\umu }{\rm m}$ . To compare the LETd optimized against reference plans, LETd and RBE-weighted dose based goals similar to and less strict than clinical ones were specified for the target. The goals for the organs at risk (OAR) remained unchanged. Robustness evaluation was studied for eight scenarios ( ± 3.5 % $\pm 3.5\%$ range uncertainty and ± 3 mm $\pm 3 \,{\rm mm}$ setup uncertainty along the main three axes). RESULTS: The optimization method with L min = 60 keV / µ m $\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ resulted in an optimal LETd distribution with an average increase of LET d , 98 % ${\rm {LET}}_{{\rm {d,}}98\%}$ (and LET d , 50 % ${\rm {LET}}_{{\rm {d,}}50\%}$ ) in the CTV-HD by 8.9 ± 1.5 keV / µ m $8.9\pm 1.5 \,{\rm keV}/{\umu }{\rm m}$ ( 27 % $27\%$ ) (and 6.9 ± 1.3 keV / µ m $6.9\pm 1.3 \,{\rm keV}/{\umu }{\rm m}$ ( 17 % $17\%$ )), without significant difference in the RBE-weighted dose. By allowing ± 5 % $\pm 5\%$ over- and under-dosage in the target, the LET d , 98 % ${\rm {LET}}_{{\rm {d,}}98\%}$ (and LET d , 50 % ${\rm {LET}}_{{\rm {d,}}50\%}$ ) can be increased by 11.3 ± 1.2 keV / µ m $11.3\pm 1.2 \,{\rm keV}/{\umu }{\rm m}$ ( 34 % $34\%$ ) (and 11.7 ± 3.4 keV / µ m $11.7\pm 3.4 \,{\rm keV}/{\umu }{\rm m}$ ( 29 % $29\%$ )), using the optimization parameters L min = 80 keV / µ m $\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ . The pass rate for the OAR goals in the LETd optimized plans was in the same level as the reference plans. LETd optimization lead to less robust plans compared to reference plans. CONCLUSIONS: Compared to conventionally optimized treatment plans, the LETd in the target was increased while maintaining the RBE-weighted dose using TPS LETd optimization functionalities. Regularly assessing RBE-weighted dose robustness and acquiring more in-room images remain crucial and inevitable aspects during treatment.

Sacral-Nerve-Sparing Planning Strategy in Pelvic Sarcomas/Chordomas Treated with Carbon-Ion Radiotherapy.

To minimize radiation-induced lumbosacral neuropathy (RILSN), we employed sacral-nerve-sparing optimized carbon-ion therapy strategy (SNSo-CIRT) in treating 35 patients with pelvic sarcomas/chordomas. Plans were optimized using Local Effect Model-I (LEM-I), prescribed DRBE|LEM-I|D50% (median dose to HD-PTV) = 73.6 (70.4-76.8) Gy (RBE)/16 fractions. Sacral nerves were contoured between L5-S3 levels. DRBE|LEM-I to 5% of sacral nerves-to-spare (outside HD-CTV) (DRBE|LEM-I|D5%) were restricted to <69 Gy (RBE). The median follow-up was 25 months (range of 2-53). Three patients (9%) developed late RILSN (≥G3) after an average period of 8 months post-CIRT. The RILSN-free survival at 2 years was 91% (CI, 81-100). With SNSo-CIRT, DRBE|LEM-I|D5% for sacral nerves-to-spare = 66.9 ± 1.9 Gy (RBE), maintaining DRBE|LEM-I to 98% of HD-CTV (DRBE|LEM-I|D98%) = 70 ± 3.6 Gy (RBE). Two-year OS and LC were 100% and 93% (CI, 84-100), respectively. LETd and DRBE with modified-microdosimetric kinetic model (mMKM) were recomputed retrospectively. DRBE|LEM-I and DRBE|mMKM were similar, but DRBE-filtered-LETd was higher in sacral nerves-to-spare in patients with RILSN than those without. At DRBE|LEM-I cutoff = 64 Gy (RBE), 2-year RILSN-free survival was 100% in patients with <12% of sacral nerves-to-spare voxels receiving LETd > 55 keV/µm than 75% (CI, 54-100) in those with ≥12% of voxels (p < 0.05). DRBE-filtered-LETd holds promise for the SNSo-CIRT strategy but requires longer follow-up for validation.

The sensitivity of radiobiological models in carbon ion radiotherapy (CIRT) and its consequences on the clinical treatment plan: Differences between LEM and MKM models.

PURPOSE: Carbon ion radiotherapy (CIRT) relies on relative biological effectiveness (RBE)-weighted dose calculations. Japanese clinics predominantly use the microdosimetric kinetic model (MKM), while European centers utilize the local effect model (LEM). Despite both models estimating RBE-distributions in tissue, their physical and mathematical assumptions differ, leading to significant disparities in RBE-weighted doses. Several European clinics adopted Japanese treatment schedules, necessitating adjustments in dose prescriptions and organ at risk (OAR) constraints. In the context of these two clinically used standards for RBE-weighted dose estimation, the objective of this study was to highlight specific scenarios for which the translations between models diverge, as shortcomings between them can influence clinical decisions. METHODS: Our aim was to discuss planning strategies minimizing those discrepancies, ultimately striving for more accurate and robust treatments. Evaluations were conducted in a virtual water phantom and patient CT-geometry, optimizing LEM RBE-weighted dose first and recomputing MKM thereafter. Dose-averaged linear energy transfer (LETd) distributions were also assessed. RESULTS: Results demonstrate how various parameters influence LEM/MKM translation. Similar LEM-dose distributions lead to markedly different MKM-dose distributions and variations in LETd. Generally, a homogeneous LEM RBE-weighted dose aligns with lower MKM values in most of the target volume. Nevertheless, paradoxical MKM hotspots may emerge (at the end of the range), potentially influencing clinical outcomes. Therefore, translation between models requires great caution. CONCLUSIONS: Understanding the relationship between these two clinical standards enables combining European and Japanese based experiences. The implementation of optimal planning strategies ensures the safety and acceptability of the clinical plan for both models and therefore enhances plan robustness from the RBE-weighted dose and LETd distribution point of view. This study emphasizes the importance of optimal planning strategies and the need for comprehensive CIRT plan quality assessment tools. In situations where simultaneous LEM and MKM computation capabilities are lacking, it can provide guidance in plan design, ultimately contributing to enhanced CIRT outcomes.

Commissioning and clinical implementation of an independent dose calculation system for scanned proton beams.

PURPOSE: Experimental patient-specific QA (PSQA) is a time and resource-intensive process, with a poor sensitivity in detecting errors. Radiation therapy facilities aim to substitute it by means of independent dose calculation (IDC) in combination with a comprehensive beam delivery QA program. This paper reports on the commissioning of the IDC software tool myQA iON (IBA Dosimetry) for proton therapy and its clinical implementation at the MedAustron Ion Therapy Center. METHODS: The IDC commissioning work included the validation of the beam model, the implementation and validation of clinical CT protocols, and the evaluation of patient treatment data. Dose difference maps, gamma index distributions, and pass rates (GPR) have been reviewed. The performance of the IDC tool has been assessed and clinical workflows, simulation settings, and GPR tolerances have been defined. RESULTS: Beam model validation showed agreement of ranges within ± 0.2 mm, Bragg-Peak widths within ± 0.1 mm, and spot sizes at various air gaps within ± 5% compared to physical measurements. Simulated dose in 2D reference fields deviated by -0.3% ± 0.5%, while 3D dose distributions differed by 1.8% on average to measurements. Validation of the CT calibration resulted in systematic differences of 2.0% between IDC and experimental data for tissue like samples. GPRs of 99.4 ± 0.6% were found for head, head and neck, and pediatric CT protocols on a 2%/2 mm gamma criterion. GPRs for the adult abdomen protocol were at 98.9% on average with 3%/3 mm. Root causes of GPR outliers, for example, implants were identified and evaluated. CONCLUSION: IDC has been successfully commissioned and integrated into the MedAustron clinical workflow for protons in 2021. IDC has been stepwise and safely substituting experimental PSQA since February 2021. The initial reduction of proton experimental PSQA was about 25% and reached up to 90% after 1 year.

Effects of nuclear interaction corrections and trichrome fragment spectra modelling on dose and linear energy transfer distributions in carbon ion radiotherapy.

Background and Purpose: Nuclear interaction correction (NIC) and trichrome fragment spectra modelling improve relative biological effectiveness-weighted dose (DRBE) and dose-averaged linear energy transfer (LETd) calculation for carbon ions. The effect of those novel approaches on the clinical dose and LET distributions was investigated. Materials and Methods: The effect of the NIC and trichrome algorithm was assessed, creating single beam plans for a virtual water phantom with standard settings and NIC + trichrome corrections. Reference DRBE and LETd distributions were simulated using FLUKA version 2021.2.9. Thirty clinically applied scanned carbon ion treatment plans were recalculated applying NIC, trichrome and NIC + trichrome corrections, using the LEM low dose approximation and compared to clinical plans (base RS). Four treatment sites were analysed: six prostate adenocarcinoma, ten head and neck, nine locally advanced pancreatic adenocarcinoma and five sacral chordoma. The FLUKA and clinical plans were compared in terms of DRBE deviations for D98%, D50%, D2% for the clinical target volume (CTV) and D50% in ring-like dose regions retrieved from isodose curves in base RS plans. Additionally, region-based median LETd deviations and global gamma parameters were evaluated. Results: Dose deviations comparing base RS and evaluation plans were within ± 1% supported by γ-pass rates over 97% for all cases. No significant LETd deviations were reported in the CTV, but significant median LETd deviations were up to 80% for very low dose regions. Conclusion: Our results showed improved accuracy of the predicted DRBE and LETd. Considering clinically relevant constraints, no significant modifications of clinical protocols are expected with the introduction of NIC + trichrome.

Validation of a deep-learning segmentation model for adult and pediatric head and neck radiotherapy in different patient positions.

Background and purpose: Autocontouring for radiotherapy has the potential to significantly save time and reduce interobserver variability. We aimed to assess the performance of a commercial autocontouring model for head and neck (H&N) patients in eight orientations relevant to particle therapy with fixed beam lines, focusing on validation and implementation for routine clinical use. Materials and methods: Autocontouring was performed on sixteen organs at risk (OARs) for 98 adult and pediatric patients with 137 H&N CT scans in eight orientations. A geometric comparison of the autocontours and manual segmentations was performed using the Hausdorff Distance 95th percentile, Dice Similarity Coefficient (DSC) and surface DSC and compared to interobserver variability where available. Additional qualitative scoring and dose-volume-histogram (DVH) parameters analyses were performed for twenty patients in two positions, consisting of scoring on a 0-3 scale based on clinical usability and comparing the mean (Dmean) and near-maximum (D2%) dose, respectively. Results: For the geometric analysis, the model performance in head-first-supine straight and hyperextended orientations was in the same range as the interobserver variability. HD95, DSC and surface DSC was heterogeneous in other orientations. No significant geometric differences were found between pediatric and adult autocontours. The qualitative scoring yielded a median score of ≥ 2 for 13/16 OARs while 7/32 DVH parameters were significantly different. Conclusions: For head-first-supine straight and hyperextended scans, we found that 13/16 OAR autocontours were suited for use in daily clinical practice and subsequently implemented. Further development is needed for other patient orientations before implementation.

Investigation on the physical dose filtered by linear energy transfer for treatment plan evaluation in carbon ion therapy.

BACKGROUND: Large tumor size has been reported as a predicting factor for inferior clinical outcome in carbon ion radiotherapy (CIRT). Besides the clinical factors accompanied with such tumors, larger tumors receive typically more low linear energy transfer (LET) contributions than small ones which may be the underlying physical cause. Although dose averaged LET is often used as a single parameter descriptor to quantify the beam quality, there is no evidence that this parameter is the optimal clinical predictor for the complex mixed radiation fields in CIRT. PURPOSE: Purpose of this study was to investigate on a novel dosimetric quantity, namely high-LET-dose ( D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the physical dose filtered based on an LET threshold) as a single parameter estimator to differentiate between carbon ion treatment plans (cTP) with a small and large tumor volume. METHODS: Ten cTPs with a planning target volume, PTV ≥ 500 cm 3 $\mathrm{PTV}\ge {500}\,{{\rm cm}^{3}}$ (large) and nine with a PTV < 500 cm 3 $\mathrm{PTV}<{500}\,{{\rm cm}^{3}}$ (small) were selected for this study. To find a reasonable LET threshold ( L thr $\textrm {L}_{\textrm {thr}}$ ) that results in a significant difference in terms of D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the voxel based normalized high-LET-dose ( D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) distribution in the clinical target volume (CTV) was studied on a subset (12 out of 19 cTPs) for 18 LET thresholds, using standard distribution descriptors (mean, variance and skewness). The classical dose volume histogram concept was used to evaluate the D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ distributions within the target of all 19 cTPs at the before determined L thr $\textrm {L}_{\textrm {thr}}$ . Statistical significance of the difference between the two groups in terms of mean D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ volume histogram parameters was evaluated by means of (two-sided) t-test or Mann-Whitney-U-test. In addition, the minimum target coverage at the above determined L thr $\textrm {L}_{\textrm {thr}}$ was compared and validated against three other thresholds to verify its potential in differentiation between small and large volume tumors. RESULTS: An L thr $\textrm {L}_{\textrm {thr}}$ of approximately 30 keV / µ m ${30}\,{\rm keV/}\umu {\rm m}$ was found to be a reasonable threshold to classify the two groups. At this threshold, the D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ were significantly larger ( p < 0.05 $p<0.05$ ) in small CTVs. For the small tumor group, the near-minimum and median D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ (and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) in the CTV were in average 9.3 ± 1.5 Gy $9.3\pm {1.5}\,{\rm Gy}$ (0.31 ± 0.08) and 13.6 ± 1.6 Gy $13.6\pm {1.6}\,{\rm Gy}$ (0.46 ± 0.06), respectively. For the large tumors, these parameters were 6.6 ± 0.2 Gy $6.6\pm {0.2}\,{\rm Gy}$ (0.20 ± 0.01) and 8.6 ± 0.4 Gy $8.6\pm {0.4}\,{\rm Gy}$ (0.28 ± 0.02). The difference between the two groups in terms of mean near-minimum and median D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ ( D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) was 2.7 Gy (11%) and 5.0 Gy (18%), respectively. CONCLUSIONS: The feasibility of high-LET-dose based evaluation was shown in this study where a lower D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ was found in cTPs with a large tumor size. Further investigation is needed to draw clinical conclusions. The proposed methodology in this work can be utilized for future high-LET-dose based studies.

Planning Strategy to Optimize the Dose-Averaged LET Distribution in Large Pelvic Sarcomas/Chordomas Treated with Carbon-Ion Radiotherapy.

To improve outcomes in large sarcomas/chordomas treated with CIRT, there has been recent interest in LET optimization. We evaluated 22 pelvic sarcoma/chordoma patients treated with CIRT [large: HD-CTV ≥ 250 cm3 (n = 9), small: HD-CTV < 250 cm3 (n = 13)], DRBE|LEM-I = 73.6 (70.4-73.6) Gy (RBE)/16 fractions, using the local effect model-I (LEM-I) optimization and modified-microdosimetric kinetic model (mMKM) recomputation. We observed that to improve high-LETd distribution in large tumors, at least 27 cm3 (low-LETd region) of HD-CTV should receive LETd of ≥33 keV/µm (p < 0.05). Hence, LETd optimization using 'distal patching' was explored in a treatment planning setting (not implemented clinically yet). Distal-patching structures were created to stop beams 1-2 cm beyond the HD-PTV-midplane. These plans were reoptimized and DRBE|LEM-I, DRBE|mMKM, and LETd were recomputed. Distal patching increased (a) LETd50% in HD-CTV (from 38 ± 3.4 keV/µm to 47 ± 8.1 keV/µm), (b) LETdmin in low-LETd regions of the HD-CTV (from 32 ± 2.3 keV/µm to 36.2 ± 3.6 keV/µm), (c) the GTV fraction receiving LETd of ≥50 keV/µm, (from <10% to >50%) and (d) the high-LETd component in the central region of the GTV, without significant compromise in DRBE distribution. However, distal patching is sensitive to setup/range uncertainties, and efforts to ascertain robustness are underway, before routine clinical implementation.

Proton or Carbon Ion Therapy for Skull Base Chordoma: Rationale and First Analysis of a Mono-Institutional Experience.

BACKGROUND: Skull base chordomas are radio-resistant tumors that require high-dose, high-precision radiotherapy, as can be delivered by particle therapy (protons and carbon ions). We performed a first clinical outcome analysis of particle therapy based on the initial 4-years of operation. METHODS: Between August 2017 and October 2021, 44 patients were treated with proton (89%) or carbon ion therapy (11%). Prior gross total resection had been performed in 21% of lesions, subtotal resection in 57%, biopsy in 12% and decompression in 10%. The average prescription dose was 75.2 Gy RBE in 37 fractions for protons and 66 Gy RBE in 22 fractions for carbon ions. RESULTS: At a median follow-up of 34.3 months (range: 1-55), 2-, and 3-year actuarial local control rates were 95.5% and 90.9%, respectively. The 2-, and 3-year overall and progression-free survival rates were 97.7%, 93.2%, 95.5% and 90.9%, respectively. The tumor volume at the time of particle therapy was highly predictive of local failure (p < 0.01), and currently, there is 100% local control in patients with tumors < 49 cc. No grade ≥3 toxicities were observed. There was no significant difference in outcome or side effect profile seen for proton versus carbon ion therapy. Five patients (11.4%) experienced transient grade ≤2 radiation-induced brain changes. CONCLUSIONS: The first analysis suggests the safety and efficacy of proton and carbon ion therapy at our center. The excellent control of small to mid-size chordomas underlines the effectiveness of particle therapy and importance of upfront maximum debulking of large lesions.

Technical note: In silico benchmarking of the linear energy transfer-based functionalities for carbon ion beams in a commercial treatment planning system.

BACKGROUND: The increasing number of studies dealing with linear energy transfer (LET)-based evaluation and optimization in the field of carbon ion radiotherapy (CIRT) indicates the rising demand for LET implementation in commercial treatment planning systems (TPS). Benchmarking studies could play a key role in detecting (and thus preventing) computation errors prior implementing such functionalities in a TPS. PURPOSE: This in silico study was conducted to benchmark the following two LET-related functionalities in a commercial TPS against Monte Carlo simulations: (1) dose averaged LET (LETd ) scoring and (2) physical dose filtration based on LET for future LET-based treatment plan evaluation and optimization studies. METHODS: The LETd scoring and LET-based dose filtering (in which the deposited dose can be separated into the dose below and above the user specified LET threshold) functionalities for carbon ions in the research version RayStation (RS) 9A-IonPG TPS (RaySearch Laboratories, Sweden) were benchmarked against GATE/Geant4 simulations. Pristine Bragg peaks (BPs) and cuboid targets, positioned at different depths in a homogeneous water phantom and a setup with heterogeneity were used for this study. RESULTS: For all setups (homogeneous and heterogeneous), the mean absolute (and relative) LETd difference was less than 1 keV/ µ $\umu$ m (3.5%) in the plateau and target and less than 2 keV/ µ $\umu$ m (8.3%) in the fragmentation tail. The maximum local differences were 4 and 6 keV/ µ $\umu$ m, respectively. The mean absolute (and relative) physical dose differences for both low- and high-LET doses were less than 1 cGy (1.5%) in the plateau, target and tail with a maximum absolute difference of 2 cGy. CONCLUSIONS: No computation error was found in the tested functionalities except for LETd in lateral direction outside the target, showing the limitation of the implemented monochrome model in the tested TPS version.

Organs at risk dose constraints in carbon ion radiotherapy at MedAustron: Translations between LEM and MKM RBE models and preliminary clinical results.

BACKGROUND: Carbon ion radiotherapy (CIRT) treatment planning is based on relative biological effectiveness (RBE) weighted dose calculations. A large amount of clinical evidence for CIRT was collected in Japan with RBE estimated by the modified microdosimetric kinetic model (MKM) while all European centres apply the first version of the local effect model (LEM). Japanese schedules have been used in Europe with adapted prescription dose and organs at risk (OAR) dose constraints. Recently, less conservative adapted LEM constraints have been implemented in clinical practice. The aim of this study was to analyse the new set of LEM dose constraints for brain parenchyma, brainstem and optic system considering both RBE models and evaluating early clinical data. MATERIAL AND METHODS: 31 patients receiving CIRT at MedAustron were analysed using the RayStation v9A planning system by recalculating clinical LEM-based plans in MKM. Dose statistics (D1cm3, D5cm3, D0.1cm3, D0.7cm3, D10%, D20%) were extracted for relevant critical OARs. Curve fitting for those values was performed, resulting in linear quadratic translation models. Clinical and radiological toxicity was evaluated. RESULTS: Based on derived fits, currently applied LEM constraints matched recommended MKM constraints with deviations between -8% and +3.9%. For particular cases, data did not follow the expected LEM vs MKM trends resulting in outliers. Radiological (asymptomatic) toxicity was detected in two outlier cases. CONCLUSION: Respecting LEM constraints does not automatically ensure that MKM constraints are met. Constraints for both RBE models need to be fulfilled for future CIRT patients at MedAustron. Careful selection of planning strategies is essential.

Novel Carbon Ion and Proton Partial Irradiation of Recurrent Unresectable Bulky Tumors (Particle-PATHY): Early Indication of Effectiveness and Safety.

BACKGROUND: We present the early results of a novel partial bulky-tumor irradiation using particles for patients with recurrent unresectable bulky tumors who failed previous state-of-the-art treatments. METHODS: First, eleven consecutive patients were treated from March 2020 until December 2021. The targeted Bystander Tumor Volume (BTV) was created by subtracting 1 cm from Gross Tumor Volume (GTV) surface. It reflected approximately 30% of the central GTV volume and was irradiated with 30-45 Gy RBE (Relative Biological Effectiveness) in three consecutive fractions. The Peritumoral Immune Microenvironment (PIM) surrounding the GTV, containing nearby tissues, blood-lymphatic vessels and lymph nodes, was considered an organ at risk (OAR) and protected by highly conservative constraints. RESULTS: With the median follow up of 6.3 months, overall survival was 64% with a median survival of 6.7 months; 46% of patients were progression-free. The average tumor volume regression was 61% from the initial size. The symptom control rate was 91%, with an average increase of the Karnofsky Index of 20%. The abscopal effect has been observed in 60% of patients. CONCLUSIONS: Partial bulky-tumor irradiation is an effective, safe and well tolerated treatment for patients with unresectable recurrent bulky disease. Abscopal effects elucidate an immunogenic pathway contribution. Extensive tumor shrinkage in some patients might permit definitive treatment-otherwise previously impossible.

Technical note: Impact of beamline-specific particle energy spectra on clinical plans in carbon ion beam therapy.

PURPOSE: The Local Effect Model version one (LEM I) is applied clinically across Europe to quantify the relative biological effectiveness (RBE) of carbon ion beams. It requires the full particle fluence spectrum differential in energy in each voxel as input parameter. Treatment planning systems (TPSs) use beamline-specific look-up tables generated with Monte Carlo (MC) codes. In this study, the changes in RBE weighted dose were quantified using different levels of details in the simulation or different MC codes. METHODS: The particle fluence differential in energy was simulated with FLUKA and Geant4 at 500 depths in water in 1-mm steps for 58 initial carbon ion energies (between 120.0 and 402.8 MeV/u). A dedicated beam model was applied, including the full description of the Nozzle using GATE-RTionV1.0 (Geant4.10.03p03). In addition, two tables generated with FLUKA were compared. The starting points of the FLUKA simulations were phase space (PhS) files from, firstly, the Geant4 nozzle simulations, and secondly, a clinical beam model where an analytic approach was used to mimic the beamline. Treatment plans (TPs) were generated with RayStation 8B (RaySearch Laboratories AB, Sweden) for cubic targets in water and 10 clinical patient cases using the clinical beam model. Subsequently, the RBE weighted dose was re-computed using the two other fluence tables (FLUKA PhS or Geant4). RESULTS: The fluence spectra of the primary and secondary particles simulated with Geant4 and FLUKA generally agreed well for the primary particles. Differences were mainly observed for the secondary particles. Interchanging the two energy spectra (FLUKA vs. GEANT4) to calculate the RBE weighted dose distributions resulted in average deviations of less than 1% in the entrance up to the end of the target region, with a maximum local deviation at the distal edge of the target. In the fragment tail, larger discrepancies of up to 5% on average were found for deep-seated targets. The patient and water phantom cases demonstrated similar results. CONCLUSION: RBE weighted doses agreed well within all tested setups, confirming the clinical beam model provided by the TPS vendor. Furthermore, the results showed that the open source and generally available MC code Geant4 (in particular using GATE or GATE-RTion) can also be used to generate basic beam data required for RBE calculation in carbon ion therapy.

Carbon Ion Dose Constraints in the Head and Neck and Skull Base: Review of MedAustron Institutional Protocols.

BACKGROUND: Dose constraints are of paramount importance for the outcome of any radiotherapy treatment. In this article, we report dose-volume constraints as well as currently used fractionation schedules for carbon ion radiotherapy as applied in MedAustron (Wiener Neustadt, Austria). MATERIALS AND METHODS: For fractionation schedules, both German and Japanese regimes were used. From the clinical experience of National Institute of Radiological Sciences (Chiba, Japan) and Heidelberg Ion Therapy (Heidelberg, Germany; formerly GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany) and the work by colleagues in Centro Nazionale Adroterapia Oncologica (Pavia, Italy) recalculating the dose from the microdosimetric kinetic model to the local effect model, we have set the dose constraints for critical organs of the head and neck area. Where no clinical data was available, an educated guess was made, based on data available from photon and proton series. RESULTS: We report the constraints for the optic nerve and chiasm, brainstem, spinal cord, cochlea, brain parenchyma, salivary gland, eye and adnexa, and mandibular/maxillary bone; constraints are grouped based on a fractionation scheme (German versus Japanese) and the risk of toxicity (safe, low to middle, and middle to high). CONCLUSION: We think validation of dose constraints should present a relevant part of the activity of any carbon ion radiotherapy facility, and we anticipate future multicentric, joint evaluations.

Results of an independent dosimetry audit for scanned proton beam therapy facilities.

PURPOSE: An independent dosimetry audit based on end-to-end testing of the entire chain of radiation therapy delivery is highly recommended to ensure consistent treatments among proton therapy centers. This study presents an auditing methodology developed by the MedAustron Ion Beam Therapy Center (Austria) in collaboration with the National Physical Laboratory (UK) and audit results for five scanned proton beam therapy facilities in Europe. METHODS: The audit procedure used a homogeneous and an anthropomorphic head phantom. The phantoms were loaded either with an ionization chamber or with alanine pellets and radiochromic films. Homogeneously planned doses of 10Gy were delivered to a box-like target volume in the homogeneous phantom and to two clinical scenarios with increasing complexity in the head phantom. RESULTS: For all tests the mean of the local differences of the absolute dose to water determined with the alanine pellets compared to the predicted dose from the treatment planning system installed at the audited institution was determined. The mean value taken over all tests performed was -0.1±1.0%. The measurements carried out with the ionization chamber were consistent with the dose determined by the alanine pellets with a mean deviation of -0.5±0.6%. CONCLUSION: The developed dosimetry audit method was successfully applied at five proton centers including various "turn-key" Cyclotron solutions by IBA, Varian and Mevion. This independent audit with extension to other tumour sites and use of the correspondent anthropomorphic phantoms may be proposed as part of a credentialing procedure for future clinical trials in proton beam therapy.

Beam monitor calibration of a synchrotron-based scanned light-ion beam delivery system.

PURPOSE: This paper presents the implementation and comparison of two independent methods of beam monitor calibration in terms of number of particles for scanned proton and carbon ion beams. METHODS: In the first method, called the single-layer method, dose-area-product to water (DAPw) is derived from the absorbed dose to water determined using a Roos-type plane-parallel ionization chamber in single-energy scanned beams. This is considered the reference method for the beam monitor calibration in the clinically relevant proton and carbon energy ranges. In the second method, called the single-spot method, DAPw of a single central spot is determined using a Bragg-peak (BP) type large-area plane-parallel ionization chamber. Emphasis is given to the detailed characterization of the ionization chambers used for the beam monitor calibration. For both methods a detailed uncertainty budget on the DAPw determination is provided as well as on the derivation of the number of particles. RESULTS: Both calibration methods agreed on average within 1.1% for protons and within 2.6% for carbon ions. The uncertainty on DAPw using single-layer beams is 2.1% for protons and 3.1% for carbon ions with major contributions from the available values of kQ and the average spot spacing in both lateral directions. The uncertainty using the single-spot method is 2.2% for protons and 3.2% for carbon ions with major contributions from the available values of kQ and the non-uniformity of the BP chamber response, which can lead to a correction of up-to 3.2%. For the number of particles, an additional dominant uncertainty component for the mean stopping power per incident proton (or the CEMA) needs to be added. CONCLUSION: The agreement between both methods enhances confidence in the beam monitor calibration and the estimated uncertainty. The single-layer method can be used as a reference and the single-spot method is an alternative that, when more accumulated knowledge and data on the method becomes available, can be used as a redundant dose monitor calibration method. This work, together with the overview of information from the literature provided here, is a first step towards comprehensive information on the single-spot method.

A GATE/Geant4 beam model for the MedAustron non-isocentric proton treatment plans quality assurance.

PURPOSE: To present a reference Monte Carlo (MC) beam model developed in GATE/Geant4 for the MedAustron fixed beam line. The proposed model includes an absolute dose calibration in Dose-Area-Product (DAP) and it has been validated within clinical tolerances for non-isocentric treatments as routinely performed at MedAustron. MATERIAL AND METHODS: The proton beam model was parametrized at the nozzle entrance considering optic and energy properties of the pencil beam. The calibration in terms of absorbed dose to water was performed exploiting the relationship between number of particles and DAP by mean of a recent formalism. Typical longitudinal dose distribution parameters (range, distal penumbra and modulation) and transverse dose distribution parameters (spot sizes, field sizes and lateral penumbra) were evaluated. The model was validated in water, considering regular-shaped dose distribution as well as clinical plans delivered in non-isocentric conditions. RESULTS: Simulated parameters agree with measurements within the clinical requirements at different air gaps. The agreement of distal and longitudinal dose distribution parameters is mostly better than 1 mm. The dose difference in reference conditions and for 3D dose delivery in water is within 0.5% and 1.2%, respectively. Clinical plans were reproduced within 3%. CONCLUSION: A full nozzle beam model for active scanning proton pencil beam is described using GATE/Geant4. Absolute dose calibration based on DAP formalism was implemented. The beam model is fully validated in water over a wide range of clinical scenarios and will be inserted as a reference tool for research and for independent dose calculation in the clinical routine.

Clinical implementation and commissioning of the MedAustron Particle Therapy Accelerator for non-isocentric scanned proton beam treatments.

PURPOSE: This paper describes the clinical implementation and medical commissioning of the MedAustron Particle Therapy Accelerator (MAPTA) for non-isocentric scanned proton beam treatments. METHODS: Medical physics involvement during technical commissioning work is presented. Acceptance testing procedures, including advanced measurement methods of intra-spill beam variations, are defined. Beam monitor calibration using two independent methods based on a dose-area product formalism is described. Emphasis is given to the medical commissioning work and the specificities related to non-isocentric irradiation, since a key feature of MedAustron is the routine delivery of non-isocentric scanned proton beam treatments. RESULTS: Key commissioning results and beam stability trend lines for more than 2 yr of clinical operation have been provided. Intra-spill beam range, size, and position variations were within specifications of 0.3 mm, 15%, and 0.5 mm, respectively. The agreement between two independent beam monitor calibration methods was better than 1.0%. Non-isocentric treatment delivery allowed lateral penumbra reduction of up to about 30%. Daily QA measurements of the beam range, size, position, and dose were always within 1 mm, 10%, 1 mm, and 2% from the baseline data, respectively. CONCLUSIONS: Non-isocentric treatments have been successfully implemented at MedAustron for routine scanned proton beam therapy using horizontal and vertical fixed beamlines. Up to now every patient was treated in non-isocentric conditions. The presented methodology to implement a new Scanned Ion Beam Delivery (SIBD) system into clinical routine for proton therapy may serve as a guidance for other centers.

Evaluation of electromagnetic and nuclear scattering models in GATE/Geant4 for proton therapy.

PURPOSE: The dose core of a proton pencil beam (PB) is enveloped by a low dose area reaching several centimeters off the central axis and containing a considerable amount of the dose. Adequate modeling of the different components of the PB profile is, therefore, required for accurate dose calculation. In this study, we experimentally validated one electromagnetic and two nuclear scattering models in GATE/Geant4 for dose calculation of proton beams in the therapeutic energy window (62-252 MeV) with and without range shifter (RaShi). METHODS: The multiple Coulomb scattering (MCS) model was validated by lateral dose core profiles measured for five energies at up to four depths from beam plateau to Bragg peak region. Nuclear halo profiles of single PBs were evaluated for three (62.4, 148.2, and 252.7 MeV) and two (97.4 and 124.7 MeV) energies, without and with RaShi, respectively. The influence of the dose core and nuclear halo on field sizes varying from 2-20 cm was evaluated by means of output factors (OFs), namely frame factors (FFs) and field size factors (FSFs), to quantify the relative increase of dose when increasing the field size. RESULTS: The relative increase in the dose core width in the simulations deviated negligibly from measurements for depths until 80% of the beam range, but was overestimated by up to 0.2 mm in σ toward the end of range for all energies. The dose halo region of the lateral dose profile agreed well with measurements in the open beam configuration, but was notably overestimated in the deepest measurement plane of the highest energy or when the beam passed through the RaShi. The root-mean-square deviations (RMSDs) between the simulated and the measured FSFs were less than 1% at all depths, but were higher in the second half of the beam range as compared to the first half or when traversing the RaShi. The deviations in one of the two tested hadron physics lists originated mostly in elastic scattering. The RMSDs could be reduced by approximately a factor of two by exchanging the default elastic scattering cross sections for protons. CONCLUSIONS: GATE/Geant4 agreed satisfyingly with most measured quantities. MCS was systematically overestimated toward the end of the beam range. Contributions from nuclear scattering were overestimated when the beam traversed the RaShi or at the depths close to the end of the beam range without RaShi. Both, field size effects and calculation uncertainties, increased when the beam traversed the RaShi. Measured field size effects were almost negligible for beams up to medium energy and were highest for the highest energy beam without RaShi, but vice versa when traversing the RaShi.

Implementation of dosimetry equipment and phantoms at the MedAustron light ion beam therapy facility.

PURPOSE: To describe the implementation of dosimetry equipment and phantoms into clinical practice of light ion beam therapy facilities. This work covers not only standard dosimetry equipment such as computerized water scanners, films, 2D-array, thimble, and plane parallel ionization chambers, but also dosimetry equipment specifically devoted to the pencil beam scanning delivery technique such as water columns, scintillating screens or multilayer ionization chambers. METHOD: Advanced acceptance testing procedures developed at MedAustron and complementary to the standard acceptance procedures proposed by the manufacturer are presented. Detailed commissioning plans have been implemented for each piece of dosimetry equipment and include an estimate of the overall uncertainty budget for the range of clinical use of each device. Some standard dosimetry equipment used in many facilities was evaluated in detail: for instance, the recombination of a 2D-array or the potential use of a microdiamond detector to measure reference transverse dose profiles in water in the core of the primary pencil beams and in the low-dose nuclear halo (over four orders of magnitude in dose). RESULTS: The implementation of dosimetry equipment as described in this work allowed determining absolute spot sizes and spot positions with an uncertainty better than 0.3 mm. Absolute ranges are determined with an uncertainty comprised of 0.2-0.6 mm, depending on the measured range and were reproduced with a maximum difference of 0.3 mm over a period of 12 months using three different devices. CONCLUSION: The detailed evaluation procedures of dosimetry equipment and phantoms proposed in this work could serve as a guidance for other medical physicists in ion beam therapy facilities and also in conventional radiation therapy.