Experimental validation of a FLUKA Monte Carlo simulation for carbon-ion radiotherapy monitoring via secondary ion tracking.

BACKGROUND: In-vivo monitoring methods of carbon ion radiotherapy (CIRT) includes explorations of nuclear reaction products generated by carbon-ion beams interacting with patient tissues. Our research group focuses on in-vivo monitoring of CIRT using silicon pixel detectors. Currently, we are conducting a prospective clinical trial as part of the In-Vivo Monitoring project (InViMo) at the Heidelberg Ion Beam Therapy Center (HIT) in Germany. We are using an innovative, in-house developed, non-contact fragment tracking system with seven mini-trackers based on the Timepix3 technology developed at CERN. PURPOSE: This article focuses on the implementation of the mini-tracker in Monte Carlo (MC) based on FLUKA simulations to monitor secondary charged nuclear fragments in CIRT. The main objective is to systematically evaluate the simulation accuracy for the InViMo project. METHODS: The implementation involved integrating the mini-tracker geometry and the scoring mechanism into the FLUKA MC simulation, utilizing the finely tuned HIT beam line. The systematic investigation included varying mini-tracker angles (from 15 ∘ $15^\circ$ to 45 ∘ $45^\circ$ in 5 ∘ $5^\circ$ steps) during the irradiation of a head-sized phantom with therapeutic carbon-ion pencil beams. To evaluate our implemented FLUKA framework, a comparison was made between the experimental data and data obtained from MC simulations. To ensure the fidelity of our comparison, experiments were performed at the HIT using the parameters and setup established in the simulations. RESULTS: Our research demonstrates high accuracy in reproducing characteristic behaviors and dependencies of the monitoring method in terms of fragment distributions in the mini-tracker, track angles, emission profiles, and fragment numbers. Discrepancies in the number of detected fragments between the experimental data and the data obtained from MC simulations are less than 4% for the angles of interest in the InViMo detection system. CONCLUSIONS: Our study confirms the potential of our simulation framework to investigate the performance of monitoring inter-fractional anatomical changes in patients undergoing CIRT using secondary nuclear charged fragments escaping from the irradiated patient.

ESTRO-EPTN radiation dosimetry guidelines for the acquisition of proton pencil beam modelling data.

Proton therapy (PT) is an advancing radiotherapy modality increasingly integrated into clinical settings, transitioning from research facilities to hospital environments. A critical aspect of the commissioning of a proton pencil beam scanning delivery system is the acquisition of experimental beam data for accurate beam modelling within the treatment planning system (TPS). These guidelines describe in detail the acquisition of proton pencil beam modelling data. First, it outlines the intrinsic characteristics of a proton pencil beam-energy distribution, angular-spatial distribution and particle number. Then, it lists the input data typically requested by TPSs. Finally, it describes in detail the set of experimental measurements recommended for the acquisition of proton pencil beam modelling data-integrated depth-dose curves, spot maps in air, and reference dosimetry. The rigorous characterization of these beam parameters is essential for ensuring the safe and precise delivery of proton therapy treatments.

An in-vivo treatment monitoring system for ion-beam radiotherapy based on 28 Timepix3 detectors.

Ion-beam radiotherapy is an advanced cancer treatment modality offering steep dose gradients and a high biological effectiveness. These gradients make the therapy vulnerable to patient-setup and anatomical changes between treatment fractions, which may go unnoticed. Charged fragments from nuclear interactions of the ion beam with the patient tissue may carry information about the treatment quality. Currently, the fragments escape the patient undetected. Inter-fractional in-vivo treatment monitoring based on these charged nuclear fragments could make ion-beam therapy safer and more efficient. We developed an ion-beam monitoring system based on 28 hybrid silicon pixel detectors (Timepix3) to measure the distribution of fragment origins in three dimensions. The system design choices as well as the ion-beam monitoring performance measurements are presented in this manuscript. A spatial resolution of 4 mm along the beam axis was achieved for the measurement of individual fragment origins. Beam-range shifts of 1.5 mm were identified in a clinically realistic treatment scenario with an anthropomorphic head phantom. The monitoring system is currently being used in a prospective clinical trial at the Heidelberg Ion Beam Therapy Centre for head-and-neck as well as central nervous system cancer patients.

Compliance of volunteers in a fully-enclosed patient rotation system for MR-guided radiation therapy: a prospective study.

BACKGROUND: Particle therapy makes a noteworthy contribution in the treatment of tumor diseases. In order to be able to irradiate from different angles, usually expensive, complex and large gantries are used. Instead rotating the beam via a gantry, the patient itself might be rotated. Here we present tolerance and compliance of volunteers for a fully-enclosed patient rotation system in a clinical magnetic resonance (MR)-scanner for potential use in MR-guided radiotherapy, conducted within a prospective evaluation study. METHODS: A patient rotation system was used to simulate and perform magnetic resonance imaging (MRI)-examinations with 50 volunteers without an oncological question. For 20 participants, the MR-examination within the bore was simulated by introducing realistic MRI noise, whereas 30 participants received an examination with image acquisition. Initially, body parameters and claustrophobia were assessed. The subjects were then rotated to different angles for simulation (0°, 45°, 90°, 180°) and imaging (0°, 70°, 90°, 110°). At each angle, anxiety and motion sickness were assessed using a 6-item State-Trait-Anxiety-Inventory (STAI-6) and a modified Motion Sickness Assessment Questionnaire (MSAQ). In addition, general areas of discomfort were evaluated. RESULTS: Out of 50 subjects, three (6%) subjects terminated the study prematurely. One subject dropped out during simulation due to nausea while rotating to 45°. During imaging, further two subjects dropped out due to shoulder pain from positioning at 90° and 110°, respectively. The average result for claustrophobia (0 = no claustrophobia to 4 = extreme claustrophobia) was none to light claustrophobia (average score: simulation 0.64 ± 0.33, imaging 0.51 ± 0.39). The mean anxiety scores (0% = no anxiety to 100% = maximal anxiety) were 11.04% (simulation) and 15.82% (imaging). Mean motion sickness scores (0% = no motion sickness to 100% = maximal motion sickness) of 3.5% (simulation) and 6.76% (imaging) were obtained across all participants. CONCLUSION: Our study proves the feasibility of horizontal rotation in a fully-enclosed rotation system within an MR-scanner. Anxiety scores were low and motion sickness was only a minor influence. Both anxiety and motion sickness showed no angular dependency. Further optimizations with regard to immobilization in the rotation device may increase subject comfort.

Assessment of fluence- and dose-averaged linear energy transfer with passive luminescence detectors in clinical proton beams.

Objective.This work investigates the use of passive luminescence detectors to determine different types of averaged linear energy transfer (LET-) for the energies relevant to proton therapy. The experimental results are compared to reference values obtained from Monte Carlo simulations.Approach.Optically stimulated luminescence detectors (OSLDs), fluorescent nuclear track detectors (FNTDs), and two different groups of thermoluminescence detectors (TLDs) were irradiated at four different radiation qualities. For each irradiation, the fluence- (LET-f) and dose-averaged LET (LET-d) were determined. For both quantities, two sub-types of averages were calculated, either considering the contributions from primary and secondary protons or from all protons and heavier, charged particles. Both simulated and experimental data were used in combination with a phenomenological model to estimate the relative biological effectiveness (RBE).Main results.All types ofLET-could be assessed with the luminescence detectors. The experimental determination ofLET-fis in agreement with reference data obtained from simulations across all measurement techniques and types of averaging. On the other hand,LET-dcan present challenges as a radiation quality metric to describe the detector response in mixed particle fields. However, excluding secondaries heavier than protons from theLET-dcalculation, as their contribution to the luminescence is suppressed by ionization quenching, leads to equal accuracy betweenLET-fandLET-d. Assessment of RBE through the experimentally determinedLET-dvalues agrees with independently acquired reference values, indicating that the investigated detectors can determineLET-with sufficient accuracy for proton therapy.Significance.OSLDs, TLDs, and FNTDs can be used to determineLET-and RBE in proton therapy. With the capability to determine dose through ionization quenching corrections derived fromLET-, OSLDs and TLDs can simultaneously ascertain dose,LET-, and RBE. This makes passive detectors appealing for measurements in phantoms to facilitate validation of clinical treatment plans or experiments related to proton therapy.

The potential of mixed carbon-helium beams for online treatment verification: a simulation and treatment planning study.

Objective.Recently, a new and promising approach for range verification was proposed. This method requires the use of two different ion species. Due to their equal magnetic rigidity, fully ionized carbon and helium ions can be simultaneously accelerated in accelerators like synchrotrons. At sufficiently high treatment energies, helium ions can exit the patient distally, reaching approximately three times the range of carbon ions at an equal energy per nucleon. Therefore, the proposal involves adding a small helium fluence to the carbon ion beam and utilizing helium as an online range probe during radiation therapy. This work aims to develop a software framework for treatment planning and motion verification in range-guided radiation therapy using mixed carbon-helium beams.Approach.The developed framework is based on the open-source treatment planning toolkit matRad. Dose distributions and helium radiographs were simulated using the open-source Monte Carlo package TOPAS. Beam delivery system parameters were obtained from the Heidelberg Ion Therapy Center, and imaging detectors along with reconstruction were facilitated by ProtonVDA. Methods for reconstructing the most likely patient positioning error scenarios and the motion phase of 4DCT are presented for prostate and lung cancer sites.Main results.The developed framework provides the capability to calculate and optimize treatment plans for mixed carbon-helium ion therapy. It can simulate the treatment process and generate helium radiographs for simulated patient geometry, including small beam views. Furthermore, motion reconstruction based on these radiographs seems possible with preliminary validation.Significance.The developed framework can be applied for further experimental work with the promising mixed carbon-helium ion implementation of range-guided radiotherapy. It offers opportunities for adaptation in particle therapy, improving dose accumulation, and enabling patient anatomy reconstruction during radiotherapy.

Carbon ion radiotherapy of hepatocellular carcinoma provides excellent local control: The prospective phase I PROMETHEUS trial.

Background & Aims: Inoperable hepatocellular carcinoma (HCC) can be treated by stereotactic body radiotherapy. However, carbon ion radiotherapy (CIRT) is more effective for sparing non-tumorous liver. High linear energy transfer could promote therapy efficacy. Japanese and Chinese studies on hypofractionated CIRT have yielded excellent results. Because of different radiobiological models and the different etiological spectrum of HCC, applicability of these results to European cohorts and centers remains questionable. The aim of this prospective study was to assess safety and efficacy and to determine the optimal dose of CIRT with active raster scanning based on the local effect model (LEM) I. Methods: CIRT was performed every other day in four fractions with relative biological effectiveness (RBE)-weighted fraction doses of 8.1-10.5 Gy (total doses 32.4-42.0 Gy [RBE]). Dose escalation was performed in five dose levels with at least three patients each. The primary endpoint was acute toxicity after 4 weeks. Results: Twenty patients received CIRT (median age 74.7 years, n = 16 with liver cirrhosis, Child-Pugh scores [CP] A5 [n = 10], A6 [n = 4], B8 [n = 1], and B9 [n = 1]). Median follow up was 23 months. No dose-limiting toxicities and no toxicities exceeding grade II occurred, except one grade III gamma-glutamyltransferase elevation 12 months after CIRT, synchronous to out-of-field hepatic progression. During 12 months after CIRT, no CP elevation occurred. The highest dose level could be applied safely. No local recurrence developed during follow up. The objective response rate was 80%. Median overall survival was 30.8 months (1/2/3 years: 75%/64%/22%). Median progression-free survival was 20.9 months (1/2/3 years: 59%/43%/43%). Intrahepatic progression outside of the CIRT target volume was the most frequent pattern of progression. Conclusions: CIRT of HCC yields excellent local control without dose-limiting toxicity. Impact and implications: To date, safety and efficacy of carbon ion radiotherapy for hepatocellular carcinoma have only been evaluated prospectively in Japanese and Chinese studies. The optimal dose and fractionation when using the local effect model for radiotherapy planning are unknown. The results are of particular interest for European and American particle therapy centers, but also of relevance for all specialists involved in the treatment and care of patients with hepatocellular carcinoma, as we present the first prospective data on carbon ion radiotherapy in hepatocellular carcinoma outside of Asia. The excellent local control should encourage further use of carbon ion radiotherapy for hepatocellular carcinoma and design of randomized controlled trials. Clinical Trials Registration: The study is registered at ClinicalTrials.gov (NCT01167374).

Dosimetric study for breathing-induced motion effects in an abdominal pancreas phantom for carbon ion mini-beam radiotherapy.

BACKGROUND: Particle mini-beam therapy exhibits promise in sparing healthy tissue through spatial fractionation, particularly notable for heavy ions, further enhancing the already favorable differential biological effectiveness at both target and entrance regions. However, breathing-induced organ motion affects particle mini-beam irradiation schemes since the organ displacements exceed the mini-beam structure dimensions, decreasing the advantages of spatial fractionation. PURPOSE: In this study, the impact of breathing-induced organ motion on the dose distribution was examined at the target and organs at risk(OARs) during carbon ion mini-beam irradiation for pancreatic cancer. METHODS: As a first step, the carbon ion mini-beam pattern was characterized with Monte Carlo simulations. To analyze the impact of breathing-induced organ motion on the dose distribution of a virtual pancreas tumor as target and related OARs, the anthropomorphic Pancreas Phantom for Ion beam Therapy (PPIeT) was irradiated with carbon ions. A mini-beam collimator was used to deliver a spatially fractionated dose distribution. During irradiation, varying breathing motion amplitudes were induced, ranging from 5 to 15 mm. Post-irradiation, the 2D dose pattern was analyzed, focusing on the full width at half maximum (FWHM), center-to-center distance (ctc), and the peak-to-valley dose ratio (PVDR). RESULTS: The mini-beam pattern was visible within OARs, while in the virtual pancreas tumor a more homogeneous dose distribution was achieved. Applied motion affected the mini-beam pattern within the kidney, one of the OARs, reducing the PVDR from 3.78  ± $\pm$  0.12 to 1.478  ± $\pm$  0.070 for the 15 mm motion amplitude. In the immobile OARs including the spine and the skin at the back, the PVDR did not change within 3.4% comparing reference and motion conditions. CONCLUSIONS: This study provides an initial understanding of how breathing-induced organ motion affects spatial fractionation during carbon ion irradiation, using an anthropomorphic phantom. A decrease in the PVDR was observed in the right kidney when breathing-induced motion was applied, potentially increasing the risk of damage to OARs. Therefore, further studies are needed to explore the clinical viability of mini-beam radiotherapy with carbon ions when irradiating abdominal regions.

Clinical Workflow of Cone Beam Computer Tomography-Based Daily Online Adaptive Radiotherapy with Offline Magnetic Resonance Guidance: The Modular Adaptive Radiotherapy System (MARS).

PURPOSE: The Ethos (Varian Medical Systems) radiotherapy device combines semi-automated anatomy detection and plan generation for cone beam computer tomography (CBCT)-based daily online adaptive radiotherapy (oART). However, CBCT offers less soft tissue contrast than magnetic resonance imaging (MRI). This work aims to present the clinical workflow of CBCT-based oART with shuttle-based offline MR guidance. METHODS: From February to November 2023, 31 patients underwent radiotherapy on the Ethos (Varian, Palo Alto, CA, USA) system with machine learning (ML)-supported daily oART. Moreover, patients received weekly MRI in treatment position, which was utilized for daily plan adaptation, via a shuttle-based system. Initial and adapted treatment plans were generated using the Ethos treatment planning system. Patient clinical data, fractional session times (MRI + shuttle transport + positioning, adaptation, QA, RT delivery) and plan selection were assessed for all fractions in all patients. RESULTS: In total, 737 oART fractions were applied and 118 MRIs for offline MR guidance were acquired. Primary sites of tumors were prostate (n = 16), lung (n = 7), cervix (n = 5), bladder (n = 1) and endometrium (n = 2). The treatment was completed in all patients. The median MRI acquisition time including shuttle transport and positioning to initiation of the Ethos adaptive session was 53.6 min (IQR 46.5-63.4). The median total treatment time without MRI was 30.7 min (IQR 24.7-39.2). Separately, median adaptation, plan QA and RT times were 24.3 min (IQR 18.6-32.2), 0.4 min (IQR 0.3-1,0) and 5.3 min (IQR 4.5-6.7), respectively. The adapted plan was chosen over the scheduled plan in 97.7% of cases. CONCLUSION: This study describes the first workflow to date of a CBCT-based oART combined with a shuttle-based offline approach for MR guidance. The oART duration times reported resemble the range shown by previous publications for first clinical experiences with the Ethos system.

Linear Energy Transfer Measurements and Estimation of Relative Biological Effectiveness in Proton and Helium Ion Beams Using Fluorescent Nuclear Track Detectors.

PURPOSE: Our objective was to develop a methodology for assessing the linear energy transfer (LET) and relative biological effectiveness (RBE) in clinical proton and helium ion beams using fluorescent nuclear track detectors (FNTDs). METHODS AND MATERIALS: FNTDs were exposed behind solid water to proton and helium (4He) ion spread-out Bragg peaks. Detectors were imaged with a confocal microscope, and the LET spectra were derived from the fluorescence intensity. The track- and dose-averaged LET (LETF and LETD, respectively) were calculated from the LET spectra. LET measurements were used as input on RBE models to estimate the RBE. Human alveolar adenocarcinoma cells (A549) were exposed at the same positions as the FNTDs. The RBE was calculated from the resulting survival curves. All measurements were compared with Monte Carlo simulations. RESULTS: For protons, average relative differences between measurements and simulations were 6% and 19% for LETF and LETD, respectively. For helium ions, the same differences were 11% for both quantities. The position of the experimental LET spectra primary peaks agreed with the simulations within 9% and 14% for protons and helium ions, respectively. For the RBE models using LETD as input, FNTD-based RBE values ranged from 1.02 ± 0.01 to 1.25 ± 0.04 and from 1.08 ± 0.09 to 2.68 ± 1.26 for protons and helium ions, respectively. The average relative differences between these values and simulations were 2% and 4%. For A549 cells, the RBE ranged from 1.05 ± 0.07 to 1.47 ± 0.09 and from 0.89 ± 0.06 to 3.28 ± 0.20 for protons and helium ions, respectively. Regarding the RBE-weighted dose (2.0 Gy at the spread-out Bragg peak), the differences between simulations and measurements were below 0.10 Gy. CONCLUSIONS: This study demonstrates for the first time that FNTDs can be used to perform direct LET measurements and to estimate the RBE in clinical proton and helium ion beams.

Patterns of Temporal Lobe Reaction and Radiation Necrosis after Particle Radiotherapy in Patients with Skull Base Chordoma and Chondrosarcoma-A Single-Center Experience.

BACKGROUND: The current study aims to evaluate the occurrence of temporal lobe reactions and identify possible risk factors for patients who underwent particle therapy of the skull base. METHODS: 244 patients treated for skull base chordoma (n = 144) or chondrosarcoma (n = 100) at the Heidelberg Ion Beam Therapy Center (HIT) using a raster scan technique, were analyzed. Follow-up MRI-scans were matched with the initial planning images. Radiogenic reactions were contoured and analyzed based on volume and dose of treatment. RESULTS: 51 patients with chordoma (35.4%) and 30 patients (30%) with chondrosarcoma experienced at least one temporal lobe reaction within the follow-up period (median 49 months for chondrosarcoma, 62 months for chordoma). Age, irradiated volume, and dose values were significant risk factors for the development of temporal lobe reactions with the highest significance for the value of DMax-7 being defined as the dose maximum in the temporal lobe minus the 7cc with the highest dose (p = 0.000000000019; OR 1.087). CONCLUSION: Temporal lobe reactions are a common side effect after particle therapy of the skull base. We were able to develop a multivariate model, which predicted radiation reactions with a specificity of 99% and a sensitivity of 52.2%.

Daily AI-Based Treatment Adaptation under Weekly Offline MR Guidance in Chemoradiotherapy for Cervical Cancer 1: The AIM-C1 Trial.

(1) Background: External beam radiotherapy (EBRT) and concurrent chemotherapy, followed by brachytherapy (BT), offer a standard of care for patients with locally advanced cervical carcinoma. Conventionally, large safety margins are required to compensate for organ movement, potentially increasing toxicity. Lately, daily high-quality cone beam CT (CBCT)-guided adaptive radiotherapy, aided by artificial intelligence (AI), became clinically available. Thus, online treatment plans can be adapted to the current position of the tumor and the adjacent organs at risk (OAR), while the patient is lying on the treatment couch. We sought to evaluate the potential of this new technology, including a weekly shuttle-based 3T-MRI scan in various treatment positions for tumor evaluation and for decreasing treatment-related side effects. (2) Methods: This is a prospective one-armed phase-II trial consisting of 40 patients with cervical carcinoma (FIGO IB-IIIC1) with an age ≥ 18 years and a Karnofsky performance score ≥ 70%. EBRT (45-50.4 Gy in 25-28 fractions with 55.0-58.8 Gy simultaneous integrated boosts to lymph node metastases) will be accompanied by weekly shuttle-based MRIs. Concurrent platinum-based chemotherapy will be given, followed by 28 Gy of BT (four fractions). The primary endpoint will be the occurrence of overall early bowel and bladder toxicity CTCAE grade 2 or higher (CTCAE v5.0). Secondary outcomes include clinical feasibility, quality of life, and imaging-based response assessment.

Segmentation of 71 Anatomical Structures Necessary for the Evaluation of Guideline-Conforming Clinical Target Volumes in Head and Neck Cancers.

The delineation of the clinical target volumes (CTVs) for radiation therapy is time-consuming, requires intensive training and shows high inter-observer variability. Supervised deep-learning methods depend heavily on consistent training data; thus, State-of-the-Art research focuses on making CTV labels more homogeneous and strictly bounding them to current standards. International consensus expert guidelines standardize CTV delineation by conditioning the extension of the clinical target volume on the surrounding anatomical structures. Training strategies that directly follow the construction rules given in the expert guidelines or the possibility of quantifying the conformance of manually drawn contours to the guidelines are still missing. Seventy-one anatomical structures that are relevant to CTV delineation in head- and neck-cancer patients, according to the expert guidelines, were segmented on 104 computed tomography scans, to assess the possibility of automating their segmentation by State-of-the-Art deep learning methods. All 71 anatomical structures were subdivided into three subsets of non-overlapping structures, and a 3D nnU-Net model with five-fold cross-validation was trained for each subset, to automatically segment the structures on planning computed tomography scans. We report the DICE, Hausdorff distance and surface DICE for 71 + 5 anatomical structures, for most of which no previous segmentation accuracies have been reported. For those structures for which prediction values have been reported, our segmentation accuracy matched or exceeded the reported values. The predictions from our models were always better than those predicted by the TotalSegmentator. The sDICE with 2 mm margin was larger than 80% for almost all the structures. Individual structures with decreased segmentation accuracy are analyzed and discussed with respect to their impact on the CTV delineation following the expert guidelines. No deviation is expected to affect the rule-based automation of the CTV delineation.

Energy painting: helium-beam radiography with thin detectors and multiple beam energies.

Objective.Compact ion imaging systems based on thin detectors are a promising prospect for the clinical environment since they are easily integrated into the clinical workflow. Their measurement principle is based on energy deposition instead of the conventionally measured residual energy or range. Therefore, thin detectors are limited in the water-equivalent thickness range they can image with high precision. This article presents ourenergy paintingmethod, which has been developed to render high precision imaging with thin detectors feasible even for objects with larger, clinically relevant water-equivalent thickness (WET) ranges.Approach.A detection system exclusively based on pixelated silicon Timepix detectors was used at the Heidelberg ion-beam therapy center to track single helium ions and measure their energy deposition behind the imaged object. Calibration curves were established for five initial beam energies to relate the measured energy deposition to WET. They were evaluated regarding their accuracy, precision and temporal stability. Furthermore, a 60 mm × 12 mm region of a wedge phantom was imaged quantitatively exploiting the calibrated energies and five different mono-energetic images. These mono-energetic images were combined in a pixel-by-pixel manner by averaging the WET-data weighted according to their single-ion WET precision (SIWP) and the number of contributing ions.Main result.A quantitative helium-beam radiograph of the wedge phantom with an average SIWP of 1.82(5) % over the entire WET interval from 150 mm to 220 mm was obtained. Compared to the previously used methodology, the SIWP improved by a factor of 2.49 ± 0.16. The relative stopping power value of the wedge derived from the energy-painted image matches the result from range pullback measurements with a relative deviation of only 0.4 %.Significance.The proposed method overcomes the insufficient precision for wide WET ranges when employing detection systems with thin detectors. Applying this method is an important prerequisite for imaging of patients. Hence, it advances detection systems based on energy deposition measurements towards clinical implementation.

Essential parameters needed for a U-Net-based segmentation of individual bones on planning CT images in the head and neck region using limited datasets for radiotherapy application.

Objective.The field of radiotherapy is highly marked by the lack of datasets even with the availability of public datasets. Our study uses a very limited dataset to provide insights on essential parameters needed to automatically and accurately segment individual bones on planning CT images of head and neck cancer patients.Approach.The study was conducted using 30 planning CT images of real patients acquired from 5 different cohorts. 15 cases from 4 cohorts were randomly selected as training and validation datasets while the remaining were used as test datasets. Four experimental sets were formulated to explore parameters such as background patch reduction, class-dependent augmentation and incorporation of a weight map on the loss function.Main results.Our best experimental scenario resulted in a mean Dice score of 0.93 ± 0.06 for other bones (skull, mandible, scapulae, clavicles, humeri and hyoid), 0.93 ± 0.02 for ribs and 0.88 ± 0.03 for vertebrae on 7 test cases from the same cohorts as the training datasets. We compared our proposed solution approach to a retrained nnU-Net and obtained comparable results for vertebral bones while outperforming in the correct identification of the left and right instances of ribs, scapulae, humeri and clavicles. Furthermore, we evaluated the generalization capability of our proposed model on a new cohort and the mean Dice score yielded 0.96 ± 0.10 for other bones, 0.95 ± 0.07 for ribs and 0.81 ± 0.19 for vertebrae on 8 test cases.Significance.With these insights, we are challenging the utilization of an automatic and accurate bone segmentation tool into the clinical routine of radiotherapy despite the limited training datasets.

Comparative evaluation of a surface-based respiratory monitoring system against a pressure sensor for 4DCT image reconstruction in phantoms.

Four-dimensional computed tomography (4DCT), which relies on breathing-induced motion, requires realistic surrogate information of breathing variations to reconstruct the tumor trajectory and motion variability of normal tissues accurately. Therefore, the SimRT surface-guided respiratory monitoring system has been installed on a Siemens CT scanner. This work evaluated the temporal and spatial accuracy of SimRT versus our commonly used pressure sensor, AZ-733 V. A dynamic thorax phantom was used to reproduce regular and irregular breathing patterns acquired by SimRT and Anzai. Various parameters of the recorded breathing patterns, including mean absolute deviations (MAD), Pearson correlations (PC), and tagging precision, were investigated and compared to ground-truth. Furthermore, 4DCT reconstructions were analyzed to assess the volume discrepancy, shape deformation and tumor trajectory. Compared to the ground-truth, SimRT more precisely reproduced the breathing patterns with a MAD range of 0.37 ± 0.27 and 0.92 ± 1.02 mm versus Anzai with 1.75 ± 1.54 and 5.85 ± 3.61 mm for regular and irregular breathing patterns, respectively. Additionally, SimRT provided a more robust PC of 0.994 ± 0.009 and 0.936 ± 0.062 for all investigated breathing patterns. Further, the peak and valley recognition were found to be more accurate and stable using SimRT. The comparison of tumor trajectories revealed discrepancies up to 7.2 and 2.3 mm for Anzai and SimRT, respectively. Moreover, volume discrepancies up to 1.71 ± 1.62% and 1.24 ± 2.02% were found for both Anzai and SimRT, respectively. SimRT was validated across various breathing patterns and showed a more precise and stable breathing tracking, (i) independent of the amplitude and period, (ii) and without placing any physical devices on the patient's body. These findings resulted in a more accurate temporal and spatial accuracy, thus leading to a more realistic 4DCT reconstruction and breathing-adapted treatment planning.

A Monte Carlo study on the secondary neutron generation by oxygen ion beams for radiotherapy and its comparison to lighter ions.

Objective.To study the secondary neutrons generated by primary oxygen beams for cancer treatment and compare the results to those from primary protons, helium, and carbon ions. This information can provide useful insight into the positioning of neutron detectors in phantom for future experimental dose assessments.Approach.Mono-energetic oxygen beams and spread-out Bragg peaks were simulated using the Monte Carlo particle transport codesFLUktuierende KAskade, tool for particle simulation, and Monte Carlo N-Particle, with energies within the therapeutic range. The energy and angular distribution of the secondary neutrons were quantified.Main results.The secondary neutron spectra generated by primary oxygen beams present the same qualitative trend as for other primary ions. The energy distributions resemble continuous spectra with one peak in the thermal/epithermal region, and one other peak in the fast/relativistic region, with the most probable energy ranging from 94 up to 277 MeV and maximum energies exceeding 500 MeV. The angular distribution of the secondary neutrons is mainly downstream-directed for the fast/relativistic energies, whereas the thermal/epithermal neutrons present a more isotropic propagation. When comparing the four different primary ions, there is a significant increase in the most probable energy as well as the number of secondary neutrons per primary particle when increasing the mass of the primaries.Significance.Most previous studies have only presented results of secondary neutrons generated by primary proton beams. In this work, secondary neutrons generated by primary oxygen beams are presented, and the obtained energy and angular spectra are added as supplementary material. Furthermore, a comparison of the secondary neutron generation by the different primary ions is given, which can be used as the starting point for future studies on treatment plan comparison and secondary neutron dose optimisation. The distal penumbra after the maximum dose deposition appears to be a suitable location for in-phantom dose assessments.

Characteristics of breathing-adapted gating using surface guidance for use in particle therapy: A phantom-based end-to-end test from CT simulation to dose delivery.

To account for intra-fractional tumor motion during dose delivery in radiotherapy, various treatment strategies are clinically implemented such as breathing-adapted gating and irradiating the tumor during specific breathing phases. In this work, we present a comprehensive phantom-based end-to-end test of breathing-adapted gating utilizing surface guidance for use in particle therapy. A commercial dynamic thorax phantom was used to reproduce regular and irregular breathing patterns recorded by the GateRT respiratory monitoring system. The amplitudes and periods of recorded breathing patterns were analysed and compared to planned patterns (ground-truth). In addition, the mean absolute deviations (MAD) and Pearson correlation coefficients (PCC) between the measurements and ground-truth were assessed. Measurements of gated and non-gated irradiations were also analysed with respect to dosimetry and geometry, and compared to treatment planning system (TPS). Further, the latency time of beam on/off was evaluated. Compared to the ground-truth, measurements performed with GateRT showed amplitude differences between 0.03 ± 0.02 mm and 0.26 ± 0.03 mm for regular and irregular breathing patterns, whilst periods of both breathing patterns ranged with a standard deviation between 10 and 190 ms. Furthermore, the GateRT software precisely acquired breathing patterns with a maximum MAD of 0.30 ± 0.23 mm. The PCC constantly ranged between 0.998 and 1.000. Comparisons between TPS and measured dose profiles indicated absolute mean dose deviations within institutional tolerances of ±5%. Geometrical beam characteristics also varied within our institutional tolerances of 1.5 mm. The overall time delays were <60 ms and thus within both recommended tolerances published by ESTRO and AAPM of 200 and 100 ms, respectively. In this study, a non-invasive optical surface-guided workflow including image acquisition, treatment planning, patient positioning and gated irradiation at an ion-beam gantry was investigated, and shown to be clinically viable. Based on phantom measurements, our results show a clinically-appropriate spatial, temporal, and dosimetric accuracy when using surface guidance in the clinical setting, and the results comply with international and institutional guidelines and tolerances.

A phantom to simulate organ motion and its effect on dose distribution in carbon ion therapy for pancreatic cancer.

Objective. Carbon ion radiotherapy is a promising radiation technique for malignancies like pancreatic cancer. However, organs' motion imposes challenges for achieving homogeneous dose delivery. In this study, an anthropomorphicPancreasPhantom forIon-beamTherapy (PPIeT) was developed to simulate breathing and gastrointestinal motion during radiotherapy.Approach. The developed phantom contains a pancreas, two kidneys, a duodenum, a spine and a spinal cord. The shell of the organs was 3D printed and filled with agarose-based mixtures. Hounsfield Units (HU) of PPIeTs' organs were measured by CT. The pancreas motion amplitude in cranial-caudal (CC) direction was evaluated from patients' 4D CT data. Motions within the obtained range were simulated and analyzed in PPIeT using MRI. Additionally, GI motion was mimicked by changing the volume of the duodenum and quantified by MRI. A patient-like treatment plan was calculated for carbon ions, and the phantom was irradiated in a static and moving condition. Dose measurements in the organs were performed using an ionization chamber and dosimetric films.Main results. PPIeT presented tissue equivalent HU and reproducible breathing-induced CC displacements of the pancreas between (3.98 ± 0.36) mm and a maximum of (18.19 ± 0.44) mm. The observed maximum change in distance of (14.28 ± 0.12) mm between pancreas and duodenum was consistent with findings in patients. Carbon ion irradiation revealed homogenous coverage of the virtual tumor at the pancreas in static condition with a 1% deviation from the treatment plan. Instead, the dose delivery during motion with the maximum amplitude yielded an underdosage of 21% at the target and an increased uncertainty by two orders of magnitude.Significance. A dedicated phantom was designed and developed for breathing motion assessment of dose deposition during carbon ion radiotherapy. PPIeT is a unique tool for dose verification in the pancreas and its organs at risk during end-to-end tests.