In vivo range verification in particle therapy

RayStation/GATE Monte Carlo simulation framework for verification of proton therapy based on the12N imaging

Objective.12N, having a half-life of 11 ms, is a highly effective positron emitter that can potentially provide near real-time feedback in proton therapy. There is currently no framework for comparing and validating positron emission imaging of12N. This work describes the development and validation of a Monte Carlo (MC) framework to calculate the images of12N, as well as long-lived isotopes, originating from activation by protons.Approach. The available dual-panel Biograph mCT PET scanner was modeled in GATE and validated by comparing the simulated sensitivity map with the measured one. The distributions of12N and long-lived isotopes were calculated by RayStation and used as the input of GATE simulations. The RayStation/GATE combination was verified using proton beam irradiations of homogeneous phantoms. A 120 MeV pulsed pencil beam with 108protons per pulse was used. Two-dimensional images were created from the GATE output and compared with the images based on the measurements and the 1D longitudinal projection of the full 2D image was used to calculate the12N activity range.Main results. The simulated sensitivity in the center of the FoV (5.44%) agrees well with the measured one (5.41%). The simulated and measured 2D sensitivity maps agree in good detail. The relative difference between the measured and simulated positron activity range for both12N and long-lived isotopes is less than 1%. The broadening of the12N images relative to those of the longer-lived isotopes can be understood in terms of the large positron range of12N.Significance. We developed and validated a MC framework based on RayStation/GATE to support the in-beam PET method for quality assurance of proton therapy. The inclusion of the very short-lived isotope12N makes the framework useful for developing near real-time verification. This represents a significant step towards translating12N real-time in vivo verification to the clinic.

Opportunities and challenges of upright patient positioning in radiotherapy

Objective.Upright positioning has seen a surge in interest as a means to reduce radiotherapy (RT) cost, improve patient comfort, and, in selected cases, benefit treatment quality. In particle therapy (PT) in particular, eliminating the need for a gantry can present massive cost and facility footprint reduction. This review discusses the opportunities of upright RT in perspective of the open challenges.Approach.The clinical, technical, and workflow challenges that come with the upright posture have been extracted from an extensive literature review, and the current state of the art was collected in a synergistic perspective from photon and particle therapy. Considerations on future developments and opportunities are provided.Main results.Modern image guidance is paramount to upright RT, but it is not clear which modalities are essential to acquire in upright posture. Using upright MRI or upright CT, anatomical differences between upright/recumbent postures have been observed for nearly all body sites. Patient alignment similar to recumbent positioning was achieved in small patient/volunteer cohorts with prototype upright positioning systems. Possible clinical advantages, such as reduced breathing motion in upright position, have been reported, but limited cohort sizes prevent resilient conclusions on the treatment impact. Redesign of RT equipment for upright positioning, such as immobilization accessories for various body regions, is necessary, where several innovations were recently presented. Few clinical studies in upright PT have already reported promising outcomes for head&neck patients.Significance.With more evidence for benefits of upright RT emerging, several centers worldwide, particularly in PT, are installing upright positioning devices or have commenced upright treatment. Still, many challenges and open questions remain to be addressed to embed upright positioning firmly in the modern RT landscape. Guidelines, professionals trained in upright patient positioning, and large-scale clinical studies are required to bring upright RT to fruition.

Particle Beam Radiobiology Status and Challenges: A PTCOG Radiobiology Subcommittee Report

Particle therapy (PT) represents a significant advancement in cancer treatment, precisely targeting tumor cells while sparing surrounding healthy tissues thanks to the unique depth-dose profiles of the charged particles. Furthermore, their linear energy transfer and relative biological effectiveness enhance their capability to treat radioresistant tumors, including hypoxic ones. Over the years, extensive research has paved the way for PT's clinical application, and current efforts aim to refine its efficacy and precision, minimizing the toxicities. In this regard, radiobiology research is evolving toward integrating biotechnology to advance drug discovery and radiation therapy optimization. This shift from basic radiobiology to understanding the molecular mechanisms of PT aims to expand the therapeutic window through innovative dose delivery regimens and combined therapy approaches. This review, written by over 30 contributors from various countries, provides a comprehensive look at key research areas and new developments in PT radiobiology, emphasizing the innovations and techniques transforming the field, ranging from the radiobiology of new irradiation modalities to multimodal radiation therapy and modeling efforts. We highlight both advancements and knowledge gaps, with the aim of improving the understanding and application of PT in oncology.

In-beam PET treatment monitoring of carbon therapy patients: Results of a clinical trial at CNAO

PURPOSE

Carbon ion therapy treatments can be monitored non-invasively with in-beam Positron Emission Tomography (PET). At CNAO the INSIDE in-beam PET scanner has been used in a clinical trial (NCT03662373) to monitor cancer treatments with proton and carbon therapy. In this work we present the analysis results of carbon therapy data, acquired during the first phase of the clinical trial, analyzing data of nine patients treated at CNAO for various malignant tumors in the head-and-neck region.

MATERIALS AND METHODS

The patient group contained two patients requiring replanning, and seven patients without replanning, based on established protocols. For each patient the PET images acquired along the course of treatment were compared with a reference, applying two analysis methods: the beam-eye-view (BEV) method and the γ-index analysis. Time trends in several parameters were investigated, as well as the agreement with control CTs, if available.

RESULTS

Regarding the BEV-method, the average sigma value σ was 3.7 mm of range difference distributions for patients without changes (sensitivity of the INSIDE detector). The 3D-information obtained from the BEV analysis was partly in agreement with what was observed in the control CT. The data quality and quantity was insufficient for a definite interpretation of the time trends.

CONCLUSION

We analyzed carbon therapy data acquired with the INSIDE in-beam PET detector using two analysis methods. The data allowed to evaluate sensitivity of the INSIDE detector for carbon therapy and to make several recommendations for the future.

The partial adaptation strategy for online-adaptive proton therapy: A proof of concept study in head and neck cancer patients

BACKGROUND

The accuracy of intensity-modulated proton therapy (IMPT) is greatly affected by anatomy variations that might occur during the treatment course. Online plan adaptations have been proposed as a solution to intervene promptly during a treatment session once the anatomy changes are detected. The implementation of online-adaptive proton therapy (OAPT) is still hindered by time-consuming tasks in the workflow.

PURPOSE

The study introduces the novel concept of partial adaptation and aims at investigating its feasibility as a potential solution to parallelize tasks during an OAPT workflow for saving valuable in-room time.

METHODS

The proof-of-principle simulation study includes datasets from six head and neck cancer (HNC) patients, each consisting of one planning CT (pCT) and three contoured control CTs (cCTs). Robust 3-field normo-fractionated initial IMPT plans were generated on the pCTs with a standardized field configuration, delivering 66 Gy and 54 Gy to the high-risk and low-risk clinical target volume (CTVHigh and CTVLow), respectively. For each cCT, a dose-mimicking-based partial adaptation was applied: two fields were adapted on the current anatomy taking into account the background dose of the first non-adapted field supposedly delivered in the meantime. Fraction doses on the cCTs resulting from partially adapted plans with different first (non-adapted) field assignments were compared against those from non-adapted and fully adapted plans regarding target coverage and organs at risk (OARs) sparing. The robustness of partially adapted plans was also evaluated.

RESULTS

Partially adapted plans showed comparable results to fully adapted plans and were superior to non-adapted plans for both target coverage and OAR sparing. Target coverage degradation in the non-adapted plans (median D98%: 95.9% and 97.5% for CTVLow and CTVHigh, respectively) was recovered by both partial (98.0% and 98.5%) and full adaptation (98.2% and 98.7%) in comparison to the initial plans (98.7% and 98.8%). The initial hotspot dose for the CTVHigh (median D2%: 101.8%) increased in the non-adapted plans (102.9%) and was recovered by the adaptive strategies (partial: 102.5%, full: 101.9%). The near-maximum dose (D0.01cc) to brainstem and spinal cord was within clinical constraints for all investigated dose distributions, but clearly increased for no adaptation and improved in the (both partially and fully) adapted plans with respect to the non-adapted ones. The parotids' median doses (D50) were mainly patient-specific depending on the proximity to the target region, but anyway lower for the partially and fully adapted plans compared to the non-adapted ones. OAR sparing was furthermore improved for the partially adapted plans in comparison to full adaptation. Robustness of the target dose metrics was preserved in all evaluated scenarios.

CONCLUSIONS

For OAPT of HNC patients, partial adaptation is able to generate plans of superior conformity to non-adapted plans and of comparable conformity as fully adapted plans, while having the potential to speed up the online-adaptive workflows. Thus, partial adaptation represents an intermediate approach until fast online adaptation workflows become available. Furthermore, it can be applied in workflows where online treatment verification stops the delivery and triggers an online adaptation for the remaining fraction.

Investigation of intra-fractionated range guided adaptive proton therapy (RGAPT): II. Range-shift compensated on-line treatment adaptation and verification

We previously proposed range-guided adaptive proton therapy (RGAPT) that uses mid-range treatment beams as probing beams and intra-fractionated range measurements for online adaptation. In this work, we demonstrated experimental verification and reported the dosimetric accuracy for RGAPT. A STEEV phantom was used for the experiments, and a 3 × 3 × 3 cm3cube inside the phantom was assigned to be the treatment target. We simulated three online range shift scenarios: reference, overshoot, and undershoot, by placing upstream Lucite sheets, 4, 0, and 8 that corresponded to changes of 0, 6.8, and -6.8 mm, respectively, in water-equivalent path length. The reference treatment plan was to deliver single-field uniform target doses in pencil beam scanning mode and generated on the Eclipse treatment planning system. Different numbers of mid-range layers, including single, three, and five layers, were selected as probing beams to evaluate beam range (BR) measurement accuracy in positron emission tomography (PET). Online plans were modified to adapt to BR shifts and compensate for probing beam doses. In contrast, non-adaptive plans were also delivered and compared to adaptive plans by film measurements. The mid-range probing beams of three (5.55MU) and five layers (8.71MU) yielded accurate range shift measurements in 60 s of PET acquisition with uncertainty of 0.5 mm while the single-layer probing (1.65MU) was not sufficient for measurements. The adaptive plans achieved an average gamma (2%/2 mm) passing rate of 95%. In contrast, the non-adaptive plans only had an average passing rate of 69%. RGAPT planning and delivery are feasible and verified by the experiments. The probing beam delivery, range measurements, and adaptive planning and delivery added a small increase in treatment delivery workflow time but resulted in substantial dose improvement. The three-layer mid-range probing was most suitable considering the balance of high range measurement accuracy and the low number of probing beam layers.

Investigation of intra-fractionated range guided adaptive proton therapy: I. On-line PET imaging and range measurement

Objective.Develop a prototype on-line positron emission tomography (PET) scanner and evaluate its capability of on-line imaging and intra-fractionated proton-induced radioactivity range measurement.Approach.Each detector consists of 32 × 32 array of 2 × 2 × 30 mm3Lutetium-Yttrium Oxyorthosilicate scintillators with single-scintillator-end readout through a 20 × 20 array of 3 × 3 mm2Silicon Photomultipliers. The PET can be configurated with a full-ring of 20 detectors for conventional PET imaging or a partial-ring of 18 detectors for on-line imaging and range measurement. All detector-level readout and processing electronics are attached to the backside of the system gantry and their output signals are transferred to a field-programable-gate-array based system electronics and data acquisition that can be placed 2 m away from the gantry. The PET imaging performance and radioactivity range measurement capability were evaluated by both the offline study that placed a radioactive source with known intensity and distribution within a phantom and the online study that irradiated a phantom with proton beams under different radiation and imaging conditions.Main results.The PET has 32 cm diameter and 6.5 cm axial length field-of-view (FOV), ∼2.3-5.0 mm spatial resolution within FOV, 3% sensitivity at the FOV center, 18%-30% energy resolution, and ∼9 ns coincidence time resolution. The offline study shows the PET can determine the shift of distal falloff edge position of a known radioactivity distribution with the accuracy of 0.3 ± 0.3 mm even without attenuation and scatter corrections, and online study shows the PET can measure the shift of proton-induced positron radioactive range with the accuracy of 0.6 ± 0.3 mm from the data acquired with a short-acquisition (60 s) and low-dose (5 MU) proton radiation to a human head phantom.Significance.This study demonstrated the capability of intra-fractionated PET imaging and radioactivity range measurement and will enable the investigation on the feasibility of intra-fractionated, range-shift compensated adaptive proton therapy.

Patient-specific deep learning for 3D protoacoustic image reconstruction and dose verification in proton therapy

BACKGROUND

Protoacoustic (PA) imaging has the potential to provide real-time 3D dose verification of proton therapy. However, PA images are susceptible to severe distortion due to limited angle acquisition. Our previous studies showed the potential of using deep learning to enhance PA images. As the model was trained using a limited number of patients' data, its efficacy was limited when applied to individual patients.

PURPOSE

In this study, we developed a patient-specific deep learning method for protoacoustic imaging to improve the reconstruction quality of protoacoustic imaging and the accuracy of dose verification for individual patients.

METHODS

Our method consists of two stages: in the first stage, a group model is trained from a diverse training set containing all patients, where a novel deep learning network is employed to directly reconstruct the initial pressure maps from the radiofrequency (RF) signals; in the second stage, we apply transfer learning on the pre-trained group model using patient-specific dataset derived from a novel data augmentation method to tune it into a patient-specific model. Raw PA signals were simulated based on computed tomography (CT) images and the pressure map derived from the planned dose. The reconstructed PA images were evaluated against the ground truth by using the root mean squared errors (RMSE), structural similarity index measure (SSIM) and gamma index on 10 specific prostate cancer patients. The significance level was evaluated by t-test with the p-value threshold of 0.05 compared with the results from the group model.

RESULTS

The patient-specific model achieved an average RMSE of 0.014 (p < 0.05 ${{{p}}}<{0.05}$ ), and an average SSIM of 0.981 (p < 0.05 ${{{p}}}<{0.05}$ ), out-performing the group model. Qualitative results also demonstrated that our patient-specific approach acquired better imaging quality with more details reconstructed when comparing with the group model. Dose verification achieved an average RMSE of 0.011 (p < 0.05 ${{{p}}}<{0.05}$ ), and an average SSIM of 0.995 (p < 0.05 ${{{p}}}<{0.05}$ ). Gamma index evaluation demonstrated a high agreement (97.4% [p < 0.05 ${{{p}}}<{0.05}$ ] and 97.9% [p < 0.05 ${{{p}}}<{0.05}$ ] for 1%/3  and 1%/5 mm) between the predicted and the ground truth dose maps. Our approach approximately took 6 s to reconstruct PA images for each patient, demonstrating its feasibility for online 3D dose verification for prostate proton therapy.

CONCLUSIONS

Our method demonstrated the feasibility of achieving 3D high-precision PA-based dose verification using patient-specific deep-learning approaches, which can potentially be used to guide the treatment to mitigate the impact of range uncertainty and improve the precision. Further studies are needed to validate the clinical impact of the technique.

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 of1.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.

A quantitative assessment of Geant4 for predicting the yield and distribution of positron-emitting fragments in ion beam therapy

Objective.To compare the accuracy with which different hadronic inelastic physics models across ten Geant4 Monte Carlo simulation toolkit versions can predict positron-emitting fragments produced along the beam path during carbon and oxygen ion therapy.Approach.Phantoms of polyethylene, gelatin, or poly(methyl methacrylate) were irradiated with monoenergetic carbon and oxygen ion beams. Post-irradiation, 4D PET images were acquired and parent11C,10C and15O radionuclides contributions in each voxel were determined from the extracted time activity curves. Next, the experimental configurations were simulated in Geant4 Monte Carlo versions 10.0 to 11.1, with three different fragmentation models-binary ion cascade (BIC), quantum molecular dynamics (QMD) and the Liege intranuclear cascade (INCL++) - 30 model-version combinations. Total positron annihilation and parent isotope production yields predicted by each simulation were compared between simulations and experiments using normalised mean squared error and Pearson cross-correlation coefficient. Finally, we compared the depth of the maximum positron annihilation yield and the distal point at which the positron yield decreases to 50% of peak between each model and the experimental results.Main results.Performance varied considerably across versions and models, with no one version/model combination providing the best prediction of all positron-emitting fragments in all evaluated target materials and irradiation conditions. BIC in Geant4 10.2 provided the best overall agreement with experimental results in the largest number of test cases. QMD consistently provided the best estimates of both the depth of peak positron yield (10.4 and 10.6) and the distal 50%-of-peak point (10.2), while BIC also performed well and INCL generally performed the worst across most Geant4 versions.Significance.The best predictions of the spatial distribution of positron annihilations and positron-emitting fragment production along the beam path during carbon and oxygen ion therapy was obtained using Geant4 10.2.p03 with BIC or QMD. These version/model combinations are recommended for future heavy ion therapy research.

Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study

GOAL

In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes.

METHODS/MATERIALS

We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans.

RESULTS

Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields.

CONCLUSION

The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

Measurement of the12C(p,n)12N reaction cross section below 150 MeV

Objective. Proton therapy currently faces challenges from clinical complications on organs-at-risk due to range uncertainties. To address this issue, positron emission tomography (PET) of the proton-induced11C and15O activity has been used to provide feedback on the proton range. However, this approach is not instantaneous due to the relatively long half-lives of these nuclides. An alternative nuclide,12N (half-life 11 ms), shows promise for real-timein vivoproton range verification. Development of12N imaging requires better knowledge of its production reaction cross section.Approach. The12C(p,n)12N reaction cross section was measured by detecting positron activity of graphite targets irradiated with 66.5, 120, and 150 MeV protons. A pulsed beam delivery with 0.7-2 × 108protons per pulse was used. The positron activity was measured during the beam-off periods using a dual-head Siemens Biograph mCT PET scanner. The12N production was determined from activity time histograms.Main results. The cross section was calculated for 11 energies, ranging from 23.5 to 147 MeV, using information on the experimental setup and beam delivery. Through a comprehensive uncertainty propagation analysis, a statistical uncertainty of 2.6%-5.8% and a systematic uncertainty of 3.3%-4.6% were achieved. Additionally, a comparison between measured and simulated scanner sensitivity showed a scaling factor of 1.25 (±3%). Despite this, there was an improvement in the precision of the cross section measurement compared to values reported by the only previous study.Significance. Short-lived12N imaging is promising for real-timein vivoverification of the proton range to reduce clinical complications in proton therapy. The verification procedure requires experimental knowledge of the12N production cross section for proton energies of clinical importance, to be incorporated in a Monte Carlo framework for12N imaging prediction. This study is the first to achieve a precise measurement of the12C(p,n)12N nuclear cross section for such proton energies.

Synthetic CT imaging for PET monitoring in proton therapy: a simulation study

Objective.This study addresses a fundamental limitation of in-beam positron emission tomography (IB-PET) in proton therapy: the lack of direct anatomical representation in the images it produces. We aim to overcome this shortcoming by pioneering the application of deep learning techniques to create synthetic control CT images (sCT) from combining IB-PET and planning CT scan data.Approach.We conducted simulations involving six patients who underwent irradiation with proton beams. Leveraging the architecture of a visual transformer (ViT) neural network, we developed a model to generate sCT images of these patients using the planning CT scans and the inter-fractional simulated PET activity maps during irradiation. To evaluate the model's performance, a comparison was conducted between the sCT images produced by the ViT model and the authentic control CT images-serving as the benchmark.Main results.The structural similarity index was computed at a mean value across all patients of 0.91, while the mean absolute error measured 22 Hounsfield Units (HU). Root mean squared error and peak signal-to-noise ratio values were 56 HU and 30 dB, respectively. The Dice similarity coefficient exhibited a value of 0.98. These values are comparable to or exceed those found in the literature. More than 70% of the synthetic morphological changes were found to be geometrically compatible with the ones reported in the real control CT scan.Significance.Our study presents an innovative approach to surface the hidden anatomical information of IB-PET in proton therapy. Our ViT-based model successfully generates sCT images from inter-fractional PET data and planning CT scans. Our model's performance stands on par with existing models relying on input from cone beam CT or magnetic resonance imaging, which contain more anatomical information than activity maps.

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.

Radioactivation effects of titanium caused by clinical proton beam: a simulation study

Objective. In proton beam therapy (PBT), metals in the patient body perturb the dose distribution, and their radioactivation may affect the dose distribution around the metal; however, the radioactivation effect has been not clarified with PBT. In this study, we aimed to evaluate the radioactivation effect of metal depending on proton energies and secondary neutrons with a clinical proton beam using a Monte Carlo (MC) simulation.Approach.The radionuclides produced from a titanium alloy (Ti-6Al-4V) and their radioactivity were calculated using a 210-MeV passive scattering proton beam with a 60-mm Spread-out Bragg Peak, and the deposited doses caused by the radioactivation were computed using the MC simulation. The position of metal was changed according to the proton mean energy in water. To assess neutron effects on the radioactivation, we calculated the radioactivation in following three situations: (i) full MC simulation with neutrons, (ii) simulation without secondary neutrons generated from the beamline components, and (iii) simulation without any secondary neutrons.Main results.Immediately after the irradiation, the radionuclide with the largest activity was Sc-45 m (half-life of 318 ms) regardless of the proton energy and the presence of neutrons. Total radioactivity tended to increase according to the proton energy. The accumulated dose for 24 h caused by the metal activation showed an increasing trend with the proton energy, with a maximum increase rate of 0.045% to the prescribed dose. The accumulated dose at a distance of 10 mm from the metal was lower than 1/10 of that at a distance of 1 mm.Significance.The radioactivation effect of the titanium was comprehensively evaluated in the clinical passive scattering proton beam. We expect that radioactivation effects on the clinical dose distribution would be small. We consider that these results will help the clinical handling of high-Z metals in PBT.

Investigation of contrast mechanisms for MRI phase signal-based proton beam visualization in water phantoms

PURPOSE

The low sensitivity and limitation to water phantoms of convection-dependent MRI magnitude signal-based proton beam visualization hinder its in vivo applicability in MR-integrated proton beam therapy. The purpose of the present study was, therefore, to assess possible contrast mechanisms for MRI phase signal-based proton beam visualization that can potentially be exploited to enhance the sensitivity of the method and extend its applicability to tissue materials.

METHODS

To assess whether proton beam-induced magnetic field perturbations, changes in material susceptibility or convection result in detectable changes in the MRI phase signal, water phantom characteristics, experiment timing, and imaging parameters were varied in combined irradiation and imaging experiments using a time-of-flight angiography pulse sequence on a prototype in-beam MRI scanner. Velocity encoding was used to further probe and quantify beam-induced convection.

RESULTS

MRI phase signal-based proton beam visualization proved feasible. The observed phase difference contrast was evoked by beam-induced buoyant convection with flow velocities in the mm/s range. Proton beam-induced magnetic field perturbations or changes in magnetic susceptibility did not influence the MRI phase signal. Velocity encoding was identified as a means to enhance the detection sensitivity.

CONCLUSION

Because the MRI phase difference contrast observed during proton beam irradiation of water phantoms is caused by beam-induced convection, this method will unlikely be transferable to tightly compartmentalized tissue wherein flow effects are restricted. However, strong velocity encoded pulse sequences were identified as promising candidates for the future development of MRI-based methods for water phantom-based geometric quality assurance in MR-integrated proton beam therapy.

Zinc-Doped Iron Oxide Nanoparticles as a Proton-Activatable Agent for Dose Range Verification in Proton Therapy

Proton therapy allows the treatment of specific areas and avoids the surrounding tissues. However, this technique has uncertainties in terms of the distal dose fall-off. A promising approach to studying the proton range is the use of nanoparticles as proton-activatable agents that produce detectable signals. For this, we developed an iron oxide nanoparticle doped with Zn (IONP@Zn-cit) with a hydrodynamic size of 10 nm and stability in serum. Cytotoxicity, defined as half of the surveillance, was 100 μg Zn/mL in the U251 cell line. The effect on clonogenic cell death was tested after X-ray irradiation, which suggested a radioprotective effect of these nanoparticles at low concentrations (1-10 μg Zn/mL). To evaluate the production of positron emitters and prompt-gamma signals, IONP@Zn-cit was irradiated with protons, obtaining prompt-gamma signals at the lowest measured concentration (10 mg Zn/mL). Finally, 67Ga-IONP@Zn-cit showed accumulation in the liver and spleen and an accumulation in the tumor tissue of 0.95% ID/g in a mouse model of U251 cells. These results suggest the possibility of using Zn nanoparticles as proton-activatable agents to verify the range by prompt gamma detection and face the challenges of prompt gamma detection in a specific biological situation, opening different avenues to go forward in this field.

MC TRIM Algorithm in Mandibula Phantom in Helium Therapy

Helium ion beam therapy, one of the particle therapies developed and studied in the 1950s for cancer treatment, resulted in clinical trials starting at Lawrence Berkeley National Laboratory in 1975. While proton and carbon ion therapies have been implemented in research institutions and hospitals globally after the end of the trials, progress in comprehending the physical, biological, and clinical findings of helium ion beam therapy has been limited, particularly due to its limited accessibility. Ongoing efforts aim to establish programs that evaluate the use of helium ion beams for clinical and research purposes, especially in the treatment of sensitive clinical cases. Additionally, helium ions have superior physical properties to proton beams, such as lower lateral scattering and larger LET. Moreover, they exhibit similar physical characteristics to carbon, oxygen, and neon ions, which are all used in heavy ion therapy. However, they demonstrate a sharper lateral penumbra with a lower radiobiological absence of certainties and lack the degradation of variations in dose distributions caused by excessive fragmenting of heavier-ion beams, especially at greater depths of penetration. In this context, the status and the prospective advancements of helium ion therapy are examined by investigating ionization, recoil, and lateral scattering values using MC TRIM algorithms in mandible plate phantoms designed from both tissue and previously studied biomaterials, providing an overview for dental cancer treatment. An average difference of 1.9% in the Bragg peak positions and 0.211 mm in lateral scattering was observed in both phantoms. Therefore, it is suggested that the 4He ion beam can be used in the treatment of mandibular tumors, and experimental research is recommended using the proposed biomaterial mandible plate phantom.

MRI magnitude signal-based proton beam visualisation in water phantoms reflects composite effects of beam-induced buoyant convection and radiation chemistry

Objective. Local magnetic resonance (MR) signal loss was previously observed during proton beam irradiation of free-floating water phantoms at ambient temperature using a research prototype in-beam magnetic resonance imaging (MRI) scanner. The emergence of this MR signal loss was hypothesised to be dependent on beam-induced convection. The aim of this study was therefore to unravel whether physical conditions allowing the development of convection must prevail for the beam-induced MRI signatures to emerge.Approach. The convection dependence of MRI magnitude signal-based proton beam visualisation was investigated in combined irradiation and imaging experiments using a gradient echo (GE)-based time-of-flight (ToF) angiography pulse sequence, which was first tested for its suitability for proton beam visualisation in free-floating water phantoms at ambient temperature. Subsequently, buoyant convection was selectively suppressed in water phantoms using either mechanical barriers or temperature control of water expansivity. The underlying contrast mechanism was further assessed using sagittal imaging and variation of T1 relaxation time-weighting.Main results. In the absence of convection-driven water flow, weak beam-induced MR signal changes occurred, whereas strong changes did occur when convection was not mechanically or thermally inhibited. Moreover, the degree of signal loss was found to change with the variation of T1-weighting. Consequently, beam-induced MR signal loss in free-floating water phantoms at ambient temperature does not exclusively originate from buoyant convection, but is caused by local composite effects of beam-induced motion and radiation chemistry resulting in a local change in the water T1 relaxation time.Significance. The identification of ToF angiography sequence-based proton beam visualisation in water phantoms to result from composite effects of beam-induced motion and radiation chemistry represents the starting point for the future elucidation of the currently unexplained motion-based MRI contrast mechanism and the identification of the proton beam-induced material change causing T1 relaxation time lengthening.

A Monte-Carlo-based study of a single-2D-detector proton-radiography system

PURPOSE

To assess the feasibility of a proton radiography (pRG) system based on a single thin pixelated detector for water-equivalent path length (WEPL) and relative stopping power (RSP) measurements.

METHODS

A model of a pRG system consisting of a single pixelated detector measuring energy deposition and proton fluence was investigated in a Geant4-based Monte Carlo study. At the position directly after an object traversed by a broad proton beam, spatial 2D distributions are calculated of the energy deposition in, and the number of protons entering the detector. Their ratio relates to the 2D distribution of the average stopping power of protons in the detector. The system response is calibrated against the residual range in water of the protons to provide the 2D distribution of the WEPL of the object. The WEPL distribution is converted into the distribution of the RSP of the object. Simulations have been done, where the system has been tested on 13 samples of homogeneous materials of which the RSPs have been calculated and compared with RSPs determined from simulations of residual-range-in-water, which we refer to as reference RSPs.

RESULTS

For both human-tissue- and non-human-tissue-equivalent materials, the RSPs derived with the detector agree with the reference values within 1%.

CONCLUSION

The study shows that a pRG system based on one thin pixelated detection screen has the potential to provide RSP predictions with an accuracy of 1%.

Influence of the background in Compton camera images for proton therapy treatment monitoring

Objective. Background events are one of the most relevant contributions to image degradation in Compton camera imaging for hadron therapy treatment monitoring. A study of the background and its contribution to image degradation is important to define future strategies to reduce the background in the system.Approach. In this simulation study, the percentage of different kinds of events and their contribution to the reconstructed image in a two-layer Compton camera have been evaluated. To this end, GATE v8.2 simulations of a proton beam impinging on a PMMA phantom have been carried out, for different proton beam energies and at different beam intensities.Main results. For a simulated Compton camera made of Lanthanum (III) Bromide monolithic crystals, coincidences caused by neutrons arriving from the phantom are the most common type of background produced by secondary radiations in the Compton camera, causing between 13% and 33% of the detected coincidences, depending on the beam energy. Results also show that random coincidences are a significant cause of image degradation at high beam intensities, and their influence in the reconstructed images is studied for values of the time coincidence windows from 500 ps to 100 ns.Significance. Results indicate the timing capabilities required to retrieve the fall-off position with good precision. Still, the noise observed in the image when no randoms are considered make us consider further background rejection methods.

Direct visualization of proton beam irradiation effects in liquids by MRI

The main advantage proton beams offer over photon beams in radiation therapy of cancer patients is the dose maximum at their finite range, yielding a reduction in the dose deposited in healthy tissues surrounding the tumor. Since no direct method exists to measure the beam's range during dose delivery, safety margins around the tumor are applied, compromising the dose conformality and reducing the targeting accuracy. Here, we demonstrate that online MRI can visualize the proton beam and reveal its range during irradiation of liquid-filled phantoms. A clear dependence on beam energy and current was found. These results stimulate research into novel MRI-detectable beam signatures and already find application in the geometric quality assurance for magnetic resonance-integrated proton therapy systems currently under development.

Image-Guided Proton Therapy: A Comprehensive Review

Image guidance for radiation therapy can improve the accuracy of the delivery of radiation, leading to an improved therapeutic ratio. Proton radiation is able to deliver a highly conformal dose to a target due to its advantageous dosimetric properties, including the Bragg peak. Proton therapy established the standard for daily image guidance as a means of minimizing uncertainties associated with proton treatment. With the increasing adoption of the use of proton therapy over time, image guidance systems for this modality have been changing. The unique properties of proton radiation present a number of differences in image guidance from photon therapy. This paper describes CT and MRI-based simulation and methods of daily image guidance. Developments in dose-guided radiation, upright treatment, and FLASH RT are discussed as well.

A hybrid multi-particle approach to range assessment-based treatment verification in particle therapy

Particle therapy (PT) used for cancer treatment can spare healthy tissue and reduce treatment toxicity. However, full exploitation of the dosimetric advantages of PT is not yet possible due to range uncertainties, warranting development of range-monitoring techniques. This study proposes a novel range-monitoring technique introducing the yet unexplored concept of simultaneous detection and imaging of fast neutrons and prompt-gamma rays produced in beam-tissue interactions. A quasi-monolithic organic detector array is proposed, and its feasibility for detecting range shifts in the context of proton therapy is explored through Monte Carlo simulations of realistic patient models and detector resolution effects. The results indicate that range shifts of [Formula: see text] can be detected at relatively low proton intensities ([Formula: see text] protons/spot) when spatial information obtained through imaging of both particle species are used simultaneously. This study lays the foundation for multi-particle detection and imaging systems in the context of range verification in PT.

Advances in Proton Therapy for the Management of Head and Neck Tumors

Proton therapy (PBRT) is a form of external beam radiotherapy with several dosimetric advantages compared with conventional photon (x-ray) radiotherapy. Unlike x-rays, protons deposit most of their dose over a finite range, with no exit dose, in a pattern known as the Bragg peak. Clinically, this can be exploited to optimize dose to tumors while delivering a lower integral dose to normal tissues. However, the optimal role of PBRT is not as well-defined as advanced x-ray-based techniques such as intensity-modulated radiotherapy.

A feasibility study of enhanced prompt gamma imaging for range verification in proton therapy using deep learning

Background and objective. Range uncertainty is a major concern affecting the delivery precision in proton therapy. The Compton camera (CC)-based prompt-gamma (PG) imaging is a promising technique to provide 3Din vivorange verification. However, the conventional back-projected PG images suffer from severe distortions due to the limited view of the CC, significantly limiting its clinical utility. Deep learning has demonstrated effectiveness in enhancing medical images from limited-view measurements. But different from other medical images with abundant anatomical structures, the PGs emitted along the path of a proton pencil beam take up an extremely low portion of the 3D image space, presenting both the attention and the imbalance challenge for deep learning. To solve these issues, we proposed a two-tier deep learning-based method with a novel weighted axis-projection loss to generate precise 3D PG images to achieve accurate proton range verification.Materials and methods: the proposed method consists of two models: first, a localization model is trained to define a region-of-interest (ROI) in the distorted back-projected PG image that contains the proton pencil beam; second, an enhancement model is trained to restore the true PG emissions with additional attention on the ROI. In this study, we simulated 54 proton pencil beams (energy range: 75-125 MeV, dose level: 1 × 109protons/beam and 3 × 108protons/beam) delivered at clinical dose rates (20 kMU min-1and 180 kMU min-1) in a tissue-equivalent phantom using Monte-Carlo (MC). PG detection with a CC was simulated using the MC-Plus-Detector-Effects model. Images were reconstructed using the kernel-weighted-back-projection algorithm, and were then enhanced by the proposed method.Results. The method effectively restored the 3D shape of the PG images with the proton pencil beam range clearly visible in all testing cases. Range errors were within 2 pixels (4 mm) in all directions in most cases at a higher dose level. The proposed method is fully automatic, and the enhancement takes only ∼0.26 s.Significance. Overall, this preliminary study demonstrated the feasibility of the proposed method to generate accurate 3D PG images using a deep learning framework, providing a powerful tool for high-precisionin vivorange verification of proton therapy.

A high sensitivity Cherenkov detector for prompt gamma timing and time imaging

We recently proposed a new approach for the real-time monitoring of particle therapy treatments with the goal of achieving high sensitivities on the particle range measurement already at limited counting statistics. This method extends the Prompt Gamma (PG) timing technique to obtain the PG vertex distribution from the exclusive measurement of particle Time-Of-Flight (TOF). It was previously shown, through Monte Carlo simulation, that an original data reconstruction algorithm (Prompt Gamma Time Imaging) allows to combine the response of multiple detectors placed around the target. The sensitivity of this technique depends on both the system time resolution and the beam intensity. At reduced intensities (Single Proton Regime-SPR), a millimetric proton range sensitivity can be achieved, provided the overall PG plus proton TOF can be measured with a 235 ps (FWHM) time resolution. At nominal beam intensities, a sensitivity of a few mm can still be obtained by increasing the number of incident protons included in the monitoring procedure. In this work we focus on the experimental feasibility of PGTI in SPR through the development of a multi-channel, Cherenkov-based PG detector with a targeted time resolution of 235 ps (FWHM): the TOF Imaging ARrAy (TIARA). Since PG emission is a rare phenomenon, TIARA design is led by the concomitant optimisation of its detection efficiency and Signal to Noise Ratio (SNR). The PG module that we developed is composed of a small PbF[Formula: see text] crystal coupled to a silicon photoMultiplier to provide the time stamp of the PG. This module is currently read in time coincidence with a diamond-based beam monitor placed upstream the target/patient to measure the proton time of arrival. TIARA will be eventually composed of 30 identical modules uniformly arranged around the target. The absence of a collimation system and the use of Cherenkov radiators are both crucial to increase the detection efficiency and the SNR, respectively. A first prototype of the TIARA block detector was tested with 63 MeV protons delivered from a cyclotron: a time resolution of 276 ps (FWHM) was obtained, resulting in a proton range sensitivity of 4 mm at 2[Formula: see text] with the acquisition of only 600 PGs. A second prototype was also evaluated with 148 MeV protons delivered from a synchro-cyclotron obtaining a time resolution below 167 ps (FWHM) for the gamma detector. Moreover, using two identical PG modules, it was shown that a uniform sensitivity on the PG profiles would be achievable by combining the response of gamma detectors uniformly distributed around the target. This work provides the experimental proof-of-concept for the development of a high sensitivity detector that can be used to monitor particle therapy treatments and potentially act in real-time if the irradiation does not comply to treatment plan.

Technical note: Flat panel proton radiography with a patient specific imaging field for accurate WEPL assessment

BACKGROUND

Proton radiography (PR) uses highly energetic proton beams to create images where energy loss is the main contrast mechanism. Water-equivalent path length (WEPL) measurements using flat panel PR (FP-PR) have potential for in vivo range verification. However, an accurate WEPL measurement via FP-PR requires irradiation with multiple energy layers, imposing high imaging doses.

PURPOSE

A FP-PR method is proposed for accurate WEPL determination based on a patient-specific imaging field with a reduced number of energies (n) to minimize imaging dose.

METHODS

Patient-specific FP-PRs were simulated and measured for a head and neck (HN) phantom. An energy selection algorithm estimated spot-wise the lowest energy required to cross the anatomy (Emin) using a water-equivalent thickness map. Starting from Emin, n was restricted to certain values (n = 26, 24, 22, …, 2 for simulations, n = 10 for measurements), resulting in patient-specific FP-PRs. A reference FP-PR with a complete set of energies was compared against patient-specific FP-PRs covering the whole anatomy via mean absolute WEPL differences (MAD), to evaluate the impact of the developed algorithm. WEPL accuracy of patient-specific FP-PRs was assessed using mean relative WEPL errors (MRE) with respect to measured multi-layer ionization chamber PRs (MLIC-PR) in the base of skull, brain, and neck regions.

RESULTS

MADs ranged from 2.1 mm (n = 26) to 21.0 mm (n = 2) for simulated FP-PRs, and 7.2 mm for measured FP-PRs (n = 10). WEPL differences below 1 mm were observed across the whole anatomy, except at the phantom surfaces. Measured patient-specific FP-PRs showed good agreement against MLIC-PRs, with MREs of 1.3 ± 2.0%, -0.1 ± 1.0%, and -0.1 ± 0.4% in the three regions of the phantom.

CONCLUSION

A method to obtain accurate WEPL measurements using FP-PR with a reduced number of energies selected for the individual patient anatomy was established in silico and validated experimentally. Patient-specific FP-PRs could provide means of in vivo range verification.

Analytic prediction of droplet vaporization events to estimate the precision of ultrasound-based proton range verification

BACKGROUND

The safety and efficacy of proton therapy is currently hampered by range uncertainties. The combination of ultrasound imaging with injectable radiation-sensitive superheated nanodroplets was recently proposed for in vivo range verification. The proton range can be estimated from the distribution of nanodroplet vaporization events, which is stochastically related to the stopping distribution of protons, as nanodroplets are vaporized by protons reaching their maximal LET at the end of their range.

PURPOSE

Here, we aim to estimate the range estimation precision of this technique. As for any stochastic measurement, the precision will increase with the sample size, that is, the number of detected vaporizations. Thus, we first develop and validate a model to predict the number of vaporizations, which is then applied to estimate the range verification precision for a set of conditions (droplet size, droplet concentration, and proton beam parameters).

METHODS

Starting from the thermal spike theory, we derived a model that predicts the expected number of droplet vaporizations in an irradiated sample as a function of the droplet size, concentration, and number of protons. The model was validated by irradiating phantoms consisting of size-sorted perfluorobutane droplets dispersed in an aqueous matrix. The number of protons was counted with an ionization chamber, and the droplet vaporizations were recorded and counted individually using high frame rate ultrasound imaging. After validation, the range estimate precision was determined for different conditions using a Monte Carlo algorithm.

RESULTS

A good agreement between theory and experiments was observed for the number of vaporizations, especially for large (5.8 ± 2.2 µm) and medium (3.5 ± 1.1 µm) sized droplets. The number of events was lower than expected in phantoms with small droplets (2.0 ± 0.7 µm), but still within the same order of magnitude. The inter-phantom variability was considerably larger (up to 30x) than predicted by the model. The validated model was then combined with Monte Carlo simulations, which predicted a theoretical range retrieval precision improving with the square-root of the number of vaporizations, and degrading at high beam energies due to range straggling. For single pencil beams with energies between 70 and 240 MeV, a range verification precision below 1% of the range required perfluorocarbon concentrations in the order of 0.3-2.4 µM.

CONCLUSION

We proposed and experimentally validated a model to provide a quick estimate of the number of vaporizations for a given set of conditions (droplet size, droplet concentration, and proton beam parameters). From this model, promising range verification performances were predicted for realistic perfluorocarbon concentrations. These findings are an incentive to move towards preclinical studies, which are critical to assess the achievable droplet distribution in and around the tumor, and hence the in vivo range verification precision.

Design and performance evaluation of a slit-slat camera for 2D prompt gamma imaging in proton therapy monitoring: A Monte Carlo simulation study

PURPOSE

We investigated the design of a prompt gamma camera for real-time dose delivery verification and the partial mitigation of range uncertainties.

METHODS

A slit slat (SS) camera was optimized using the trade-off between the signal-to-noise ratio and spatial resolution. Then, using the GATE Monte Carlo package, the camera performances were estimated by means of target shifts, beam position quantification, changing the camera distance from the beam, and air cavity inserting. A homogeneous PMMA phantom and the air gaps induced PMMA phantom were used. The air gaps ranged from 5 mm to 30 mm by 5 mm increments were positioned in the middle of the beam range. To reduce the simulation time, phase space scoring was used. The batch method with five realizations was used for stochastic error calculations.

RESULTS

The system's detection efficiency was 1.1 × 10 - 4 PGs Emitted PGs ( 1.8 × 10 - 5 $1.1 \times {10}^{-4}\frac{{\rm PGs}}{{\rm Emitted}\ {\rm PGs}}\ (1.8 \times {10}^{-5}$ PGs/proton) for a 10 × 20 cm² detector (source-to-collimator distance = 15.0 cm). Axial and transaxial resolutions were 23 mm and 18 mm, respectively. The SS camera estimated the range as 69.0 ± 3.4 (relative stochastic error 1-sigma is 5%) and 67.6 ± 1.8 mm (2.6%) for the real range of 67.0 mm for 10⁷ and 10⁸ protons of 100 MeV, respectively. Considering 160 MeV, these values are 155.5 ± 3.1 (2%) and 152.2 ± 2.0 mm (1.3%) for the real range of 152.0 mm for 10⁷ and 10⁸ protons, respectively. Considering phantom shift, for a 100 MeV beam, the precision of the quantification (1-sigma) in the axial and lateral phantom shift estimation is 2.6 mm and 1 mm, respectively. Accordingly, the axial and lateral quantification precisions were 1.3 mm and 1 mm for a 160 MeV beam, respectively. Furthermore, the quantification of an air gap formulated as gap d e t = 0.98 × gap real ${{\rm gap}}_{det}=0.98 \times {{\rm gap}}_{{\rm real}}$ , where gap d e t ${{\rm gap}}_{det}$ and gap_real are the estimated and real air gap, respectively. The precision of the air gap quantification is 1.6 mm (1 sigma). Moreover, 2D PG images show the trajectory of the proton beam through the phantom.

CONCLUSION

The proposed slit-slat imaging systems can potentially provide a real-time, in-vivo, and non-invasive treatment monitoring method for proton therapy.

Reconstruction of thermoacoustic emission sources induced by proton irradiation using numerical time reversal

Objective.Mapping of dose delivery in proton beam therapy can potentially be performed by analyzing thermoacoustic emissions measured by ultrasound arrays. Here, a method is derived and demonstrated for spatial mapping of thermoacoustic sources using numerical time reversal, simulating re-transmission of measured emissions into the medium.Approach.Spatial distributions of thermoacoustic emission sources are shown to be approximated by the analytic-signal form of the time-reversed acoustic field, evaluated at the time of the initial proton pulse. Given calibration of the array sensitivity and knowledge of tissue properties, this approach approximately reconstructs the acoustic source amplitude, equal to the product of the time derivative of the radiation dose rate, mass density, and Grüneisen parameter. This approach was implemented using two models for acoustic fields of the array elements, one modeling elements as line sources and the other as rectangular radiators. Thermoacoustic source reconstructions employed previously reported measurements of emissions from proton energy deposition in tissue-mimicking phantoms. For a phantom incorporating a bone layer, reconstructions accounted for the higher sound speed in bone. Dependence of reconstruction quality on array aperture size and signal-to-noise ratio was consistent with previous acoustic simulation studies.Main results.Thermoacoustic source distributions were successfully reconstructed from acoustic emissions measured by a linear ultrasound array. Spatial resolution of reconstructions was significantly improved in the azimuthal (array) direction by incorporation of array element diffraction. Source localization agreed well with Monte Carlo simulations of energy deposition, and was improved by incorporating effects of inhomogeneous sound speed.Significance.The presented numerical time reversal approach reconstructs thermoacoustic sources from proton beam radiation, based on straightforward processing of acoustic emissions measured by ultrasound arrays. This approach may be useful for ranging and dosimetry of clinical proton beams, if acoustic emissions of sufficient amplitude and bandwidth can be generated by therapeutic proton sources.

Ionoacoustic application of an optical hydrophone to detect proton beam range in water

BACKGROUND

Proton range uncertainty has been the main factor limiting the ability of proton therapy to concentrate doses to tumors to their full potential. Ionoacoustic (IA) range verification is an approach to reducing this uncertainty by detecting thermoacoustic waves emitted from an irradiated volume immediately following a pulsed proton beam delivery; however, the signal weakness has been an obstacle to its clinical application. To increase the signal-to-noise ratio (SNR) with the conventional piezoelectric hydrophone (PH), the detector-sensitive volume needs to be large, but it could narrow the range of available beam angles and disturb real-time images obtained during beam delivery.

PURPOSE

To prevent this issue, we investigated a millimeter-sized optical hydrophone (OH) that exploits the laser interferometric principle. For two types of IA waves [γ-wave emitted from the Bragg peak (BP) and a spherical IA wave with resonant frequency (SPIRE) emitted from the gold fiducial marker (GM)], comparisons were made with PH in terms of waveforms, SNR, range detection accuracy, and signal intensity robustness against the small detector misalignment, particularly for SPIRE.

METHODS

A 100-MeV proton beam with a 27 ns pulse width and 4 mm beam size was produced using a fixed-field alternating gradient accelerator and was irradiated to the water phantom. The GM was set on the beam's central axis. Acrylic plates of various thicknesses, up to 12 mm, were set in front of the phantoms to shift the proton range. OH was set distal and lateral to the beam, and the range was estimated using the time-of-flight method for γ-wave and by comparing with the calibration data (SPIRE intensity versus the distance between the GM and BP) derived from an IA wave transport simulation for SPIRE. The BP dose per pulse was 0.5-0.6 Gy. To measure the variation in SPIRE amplitude against the hydrophone misalignment, the hydrophone was shifted by ± 2 mm at a maximum in lateral directions.

RESULTS

Despite its small size, OH could detect γ-wave with a higher SNR than the conventional PH (diameter, 29 mm), and a single measurement was sufficient to detect the beam range with a submillimeter accuracy in water. In the SPIRE measurement, OH was far more robust against the detector misalignment than the focused PH (FPH) used in our previous study [5%/mm (OH) versus 80%/mm (FPH)], and the correlation between the measured SPIRE intensity and the distance between the GM and BP agreed well with the simulation results. However, the OH sensitivity was lower than the FPH sensitivity, and about 5.6-Gy dose was required to decrease the intensity variation among measurements to less than 10%.

CONCLUSION

The miniature OH was found to detect weak IA signals produced by proton beams with a BP dose used in hypofractionated regimens. The OH sensitivity improvement at the MHz regime is worth exploring as the next step.

Proton beam range verification by means of ionoacoustic measurements at clinically relevant doses using a correlation-based evaluation

Purpose

The Bragg peak located at the end of the ion beam range is one of the main advantages of ion beam therapy compared to X-Ray radiotherapy. However, verifying the exact position of the Bragg peak within the patient online is a major challenge. The goal of this work was to achieve submillimeter proton beam range verification for pulsed proton beams of an energy of up to 220 MeV using ionoacoustics for a clinically relevant dose deposition of typically 2 Gy per fraction by i) using optimal proton beam characteristics for ionoacoustic signal generation and ii) improved signal detection by correlating the signal with simulated filter templates.

Methods

A water tank was irradiated with a preclinical 20 MeV proton beam using different pulse durations ranging from 50 ns up to 1 μs in order to maximise the signal-to-noise ratio (SNR) of ionoacoustic signals. The ionoacoustic signals were measured using a piezo-electric ultrasound transducer in the MHz frequency range. The signals were filtered using a cross correlation-based signal processing algorithm utilizing simulated templates, which enhances the SNR of the recorded signals. The range of the protons is evaluated by extracting the time of flight (ToF) of the ionoacoustic signals and compared to simulations from a Monte Carlo dose engine (FLUKA).

Results

Optimised SNR of 28.0 ± 10.6 is obtained at a beam current of 4.5 μA and a pulse duration of 130 ns at a total peak dose deposition of 0.5 Gy. Evaluated ranges coincide with Monte Carlo simulations better than 0.1 mm at an absolute range of 4.21 mm. Higher beam energies require longer proton pulse durations for optimised signal generation. Using the correlation-based post-processing filter a SNR of 17.8 ± 5.5 is obtained for 220 MeV protons at a total peak dose deposition of 1.3 Gy. For this clinically relevant dose deposition and proton beam energy, submillimeter range verification was achieved at an absolute range of 303 mm in water.

Conclusion

Optimal proton pulse durations ensure an ideal trade-off between maximising the ionoacoustic amplitude and minimising dose deposition. In combination with a correlation-based post-processing evaluation algorithm, a reasonable SNR can be achieved at low dose levels putting clinical applications for online proton or ion beam range verification into reach.

Longitudinally Heterogeneous Tumor Dose Optimizes Proton Broadbeam, Interlaced Minibeam, and FLASH Therapy

Simple Summary

The aim of any kind of external radiation therapy is to control a tumor with the highest possible probability of the lowest possible side effects. Here, we study further opportunities of reducing the side effects of proton therapy by applying longitudinally heterogeneous dose distributions in the tumor respecting the delivery of a minimum prescribed dose. In our simulations, the longitudinally heterogeneous dose distributions show a reduced dose in the healthy tissue already in the case of proton broadbeam irradiations, but a much higher (calculated) mean cell survival in the case of proton minibeam irradiation. This demonstrates its potential to substantially reduce side effects at a simultaneously higher tumor control probability, opening new opportunities of easier application when striving for high dose-rate applications of proton beams (>~10 Gy/s), in order to additionally profit from the so-called FLASH effects.

Abstract

The prerequisite of any radiation therapy modality (X-ray, electron, proton, and heavy ion) is meant to meet at least a minimum prescribed dose at any location in the tumor for the best tumor control. In addition, there is also an upper dose limit within the tumor according to the International Commission on Radiation Units (ICRU) recommendations in order to spare healthy tissue as well as possible. However, healthy tissue may profit from the lower side effects when waving this upper dose limit and allowing a larger heterogeneous dose deposition in the tumor, but maintaining the prescribed minimum dose level, particularly in proton minibeam therapy. Methods: Three different longitudinally heterogeneous proton irradiation modes and a standard spread-out Bragg peak (SOBP) irradiation mode are simulated for their depth-dose curves under the constraint of maintaining a minimum prescribed dose anywhere in the tumor region. Symmetric dose distributions of two opposing directions are overlaid in a 25 cm-thick water phantom containing a 5 cm-thick tumor region. Interlaced planar minibeam dose distributions are compared to those of a broadbeam using the same longitudinal dose profiles. Results and Conclusion: All longitudinally heterogeneous proton irradiation modes show a dose reduction in the healthy tissue compared to the common SOBP mode in the case of broad proton beams. The proton minibeam cases show eventually a much larger mean cell survival and thus a further reduced equivalent uniform dose (EUD) in the healthy tissue than any broadbeam case. In fact, the irradiation mode using only one proton energy from each side shows better sparing capabilities in the healthy tissue than the common spread-out Bragg peak irradiation mode with the option of a better dose fall-off at the tumor edges and an easier technical realization, particularly in view of proton minibeam irradiation at ultra-high dose rates larger than ~10 Gy/s (so-called FLASH irradiation modes).

In-vivo range verification analysis with in-beam PET data for patients treated with proton therapy at CNAO

Morphological changes that may arise through a treatment course are probably one of the most significant sources of range uncertainty in proton therapy. Non-invasive in-vivo treatment monitoring is useful to increase treatment quality. The INSIDE in-beam Positron Emission Tomography (PET) scanner performs in-vivo range monitoring in proton and carbon therapy treatments at the National Center of Oncological Hadrontherapy (CNAO). It is currently in a clinical trial (ID: NCT03662373) and has acquired in-beam PET data during the treatment of various patients. In this work we analyze the in-beam PET (IB-PET) data of eight patients treated with proton therapy at CNAO. The goal of the analysis is twofold. First, we assess the level of experimental fluctuations in inter-fractional range differences (sensitivity) of the INSIDE PET system by studying patients without morphological changes. Second, we use the obtained results to see whether we can observe anomalously large range variations in patients where morphological changes have occurred. The sensitivity of the INSIDE IB-PET scanner was quantified as the standard deviation of the range difference distributions observed for six patients that did not show morphological changes. Inter-fractional range variations with respect to a reference distribution were estimated using the Most-Likely-Shift (MLS) method. To establish the efficacy of this method, we made a comparison with the Beam’s Eye View (BEV) method. For patients showing no morphological changes in the control CT the average range variation standard deviation was found to be 2.5 mm with the MLS method and 2.3 mm with the BEV method. On the other hand, for patients where some small anatomical changes occurred, we found larger standard deviation values. In these patients we evaluated where anomalous range differences were found and compared them with the CT. We found that the identified regions were mostly in agreement with the morphological changes seen in the CT scan.

Feasibility of triple gamma ray imaging of 10C for range verification in ion therapy

Objective. In carbon ion therapy, the visualization of the range of incident particles in a patient body is important for treatment verification. In-beam positron emission tomography (PET) imaging is one of the methods to verify the treatment in ion therapy due to the high quality of PET images. We have shown the feasibility of in-beam PET imaging of radioactive 15O and 11C ion beams for range verification using our OpenPET system. Recently, we developed a whole gamma imager (WGI) that can simultaneously work as PET, single gamma ray and triple gamma ray imaging. The WGI has high potential to detect the location of 10C, which emits positrons with a simultaneous gamma ray of 718 keV, within the patient’s body during ion therapy. Approach. In this work, we focus on investigating the performance of WGI for 10C imaging and its feasibility for range verification in carbon ion therapy. First, the performance of the WGI was studied to image a 10C point source using the Geant4 toolkit. Then, the feasibility of WGI was investigated for an irradiated polymethyl methacrylate (PMMA) phantom with a 10C ion beam at the carbon therapy facility of the Heavy Ion Medical Accelerator in Chiba. Main results. The average spatial resolution and sensitivity for the simulated 10C point source at the centre of the field of view were 5.5 mm FWHM and 0.010%, respectively. The depth dose of the 10C ion beam was measured, and the triple gamma image of 10C nuclides for an irradiated PMMA phantom was obtained by applying a simple back projection to the detected triple gammas. Significance. The shift between Bragg peak position and position of the peak of the triple gamma image in an irradiated PMMA phantom was 2.8 ± 0.8 mm, which demonstrates the capability of triple gamma imaging using WGI for range verification of 10C ion beams.

Considerations for Upright Particle Therapy Patient Positioning and Associated Image Guidance

Particle therapy is a rapidly growing field in cancer therapy. Worldwide, over 100 centers are in operation, and more are currently in construction phase. The interest in particle therapy is founded in the superior target dose conformity and healthy tissue sparing achievable through the particles’ inverse depth dose profile. This physical advantage is, however, opposed by increased complexity and cost of particle therapy facilities. Particle therapy, especially with heavier ions, requires large and costly equipment to accelerate the particles to the desired treatment energy and steer the beam to the patient. A significant portion of the cost for a treatment facility is attributed to the gantry, used to enable different beam angles around the patient for optimal healthy tissue sparing. Instead of a gantry, a rotating chair positioning system paired with a fixed horizontal beam line presents a suitable cost-efficient alternative. Chair systems have been used already at the advent of particle therapy, but were soon dismissed due to increased setup uncertainty associated with the upright position stemming from the lack of dedicated image guidance systems. Recently, treatment chairs gained renewed interest due to the improvement in beam delivery, commercial availability of vertical patient CT imaging and improved image guidance systems to mitigate the problem of anatomical motion in seated treatments. In this review, economical and clinical reasons for an upright patient positioning system are discussed. Existing designs targeted for particle therapy are reviewed, and conclusions are drawn on the design and construction of chair systems and associated image guidance. Finally, the different aspects from literature are channeled into recommendations for potential upright treatment layouts, both for retrofitting and new facilities.

Roadmap: helium ion therapy

Helium ion beam therapy for the treatment of cancer was one of several developed and studied particle treatments in the 1950s, leading to clinical trials beginning in 1975 at the Lawrence Berkeley National Laboratory. The trial shutdown was followed by decades of research and clinical silence on the topic while proton and carbon ion therapy made debuts at research facilities and academic hospitals worldwide. The lack of progression in understanding the principle facets of helium ion beam therapy in terms of physics, biological and clinical findings persists today, mainly attributable to its highly limited availability. Despite this major setback, there is an increasing focus on evaluating and establishing clinical and research programs using helium ion beams, with both therapy and imaging initiatives to supplement the clinical palette of radiotherapy in the treatment of aggressive disease and sensitive clinical cases. Moreover, due its intermediate physical and radio-biological properties between proton and carbon ion beams, helium ions may provide a streamlined economic steppingstone towards an era of widespread use of different particle species in light and heavy ion therapy. With respect to the clinical proton beams, helium ions exhibit superior physical properties such as reduced lateral scattering and range straggling with higher relative biological effectiveness (RBE) and dose-weighted linear energy transfer (LET_d) ranging from ∼4 keVμm^-1to ∼40 keVμm^-1. In the frame of heavy ion therapy using carbon, oxygen or neon ions, where LET_dincreases beyond 100 keVμm^-1, helium ions exhibit similar physical attributes such as a sharp lateral penumbra, however, with reduced radio-biological uncertainties and without potentially spoiling dose distributions due to excess fragmentation of heavier ion beams, particularly for higher penetration depths. This roadmap presents an overview of the current state-of-the-art and future directions of helium ion therapy: understanding physics and improving modeling, understanding biology and improving modeling, imaging techniques using helium ions and refining and establishing clinical approaches and aims from learned experience with protons. These topics are organized and presented into three main sections, outlining current and future tasks in establishing clinical and research programs using helium ion beams-A. Physics B. Biological and C. Clinical Perspectives.

Comparison of reconstructed prompt gamma emissions using maximum likelihood estimation and origin ensemble algorithms for a Compton camera system tailored to proton range monitoring

Compton-based prompt gamma (PG) imaging is being investigated by several groups as a potential solution for in vivo range monitoring in proton therapy. The performance of this technique depends on the detector system as well as the ability of the reconstruction method to obtain good spatial resolution to establish a quantitative correlation between the PG emission and the proton beam range in the patient. To evaluate the feasibility of PG imaging for range monitoring, we quantitatively evaluated the emission distributions reconstructed by a Maximum Likelihood Expectation Maximization (MLEM) and a Stochastic Origin Ensemble (SOE) algorithm. To this end, we exploit experimental and Monte Carlo (MC) simulation data acquired with the Polaris-J Compton Camera (CC) prototype. The differences between the proton beam range (R_D) defined as the 80% distal dose fall-off and the PG range (R_PG), obtained by fitting the distal end of the reconstructed profile with a sigmoid function, were quantified. A comparable performance of both reconstruction algorithms was found. For both experimental and simulated irradiation scenarios, the correlation between R_D and R_PG was within 5 mm. These values were consistent with the ground truth distance (R_D-R_PGg≈ 3 mm) calculated by using the expected PG emission available from MC simulation. Furthermore, shifts of 3 mm in the proton beam range were resolved with the MLEM algorithm by calculating the relative difference between the R_PG for each reconstructed profile. In non-homogeneous targets, the spatial changes in the PG emission due to the different materials could not be fully resolved from the reconstructed profiles; however, the fall-off region still resembled the ground truth emission. For this scenario, the PG correlation (R_D-R_PG) varied from 0.1 mm to 4 mm, which is close to the ground truth correlation (3 mm). This work provides a framework for the evaluation of the range monitoring capabilities of a CC device for PG imaging. The two investigated image reconstruction algorithms showed a comparable and consistent performance for homogeneous and heterogeneous targets.

Image guidance for FLASH radiotherapy

FLASH radiotherapy (FLASH-RT) is an emerging ultra-high dose (>40 Gy/s) delivery that promises to improve the therapeutic potential by limiting toxicities compared to conventional RT while maintaining similar tumor eradication efficacy. Image guidance is an essential component of modern RT that should be harnessed to meet the special emerging needs of FLASH-RT and its associated high risks in planning and delivering of such ultra-high doses in short period of times. Hence, this contribution will elaborate on the imaging requirements and possible solutions in the entire chain of FLASH-RT treatment, from the planning, through the setup and delivery with online in vivo imaging and dosimetry, up to the assessment of biological mechanisms and treatment response. In patient setup and delivery, higher temporal sampling than in conventional RT should ensure that the short treatment is delivered precisely to the targeted region. Additionally, conventional imaging tools such as cone-beam computed tomography will continue to play an important role in improving patient setup prior to delivery, while techniques based on magnetic resonance imaging or positron emission tomography may be extremely valuable for either linear accelerator (Linac) or particle FLASH therapy, to monitor and track anatomical changes during delivery. In either planning or assessing outcomes, quantitative functional imaging could supplement conventional imaging for more accurate utilization of the biological window of the FLASH effect, selecting for or verifying things such as tissue oxygen and existing or transient hypoxia on the relevant timescales of FLASH-RT delivery. Perhaps most importantly at this time, these tools might help improve the understanding of the biological mechanisms of FLASH-RT response in tumor and normal tissues. The high dose deposition of FLASH provides an opportunity to utilize pulse-to-pulse imaging tools such as Cherenkov or radiation acoustic emission imaging. These could provide individual pulse mapping or assessing the 3D dose delivery superficially or at tissue depth, respectively. In summary, the most promising components of modern RT should be used for safer application of FLASH-RT, and new promising developments could be advanced to cope with its novel demands but also exploit new opportunities in connection with the unique nature of pulsed delivery at unprecedented dose rates, opening a new era of biological image guidance and ultrafast, pulse-based in vivo dosimetry.

Acoustic Modulation Enables Proton Detection With Nanodroplets at Body Temperature

Superheated nanodroplet (ND) vaporization by proton radiation was recently demonstrated, opening the door to ultrasound-based in vivo proton range verification. However, at body temperature and physiological pressures, perfluorobutane nanodroplets (PFB-NDs), which offer a good compromise between stability and radiation sensitivity, are not directly sensitive to primary protons. Instead, they are vaporized by infrequent secondary particles, which limits the precision for range verification. The radiation-induced vaporization threshold (i.e., sensitization threshold) can be reduced by lowering the pressure in the droplet such that ND vaporization by primary protons can occur. Here, we propose to use an acoustic field to modulate the pressure, intermittently lowering the proton sensitization threshold of PFB-NDs during the rarefactional phase of the ultrasound wave. Simultaneous proton irradiation and sonication with a 1.1 MHz focused transducer, using increasing peak negative pressures (PNPs), were applied on a dilution of PFB-NDs flowing in a tube, while vaporization was acoustically monitored with a linear array. Sensitization to primary protons was achieved at temperatures between [Formula: see text] and 40 °C using acoustic PNPs of relatively low amplitude (from 800 to 200 kPa, respectively), while sonication alone did not lead to ND vaporization at those PNPs. Sensitization was also measured at the clinically relevant body temperature (i.e., 37 °C) using a PNP of 400 kPa. These findings confirm that acoustic modulation lowers the sensitization threshold of superheated NDs, enabling a direct proton response at body temperature.

Detection and discrimination of neutron capture events for NCEPT dose quantification

Neutron Capture Enhanced Particle Therapy (NCEPT) boosts the effectiveness of particle therapy by capturing thermal neutrons produced by beam-target nuclear interactions in and around the treatment site, using tumour-specific [Formula: see text]B or [Formula: see text]Gd-based neutron capture agents. Neutron captures release high-LET secondary particles together with gamma photons with energies of 478 keV or one of several energies up to 7.94 MeV, for [Formula: see text]B and [Formula: see text]Gd, respectively. A key requirement for NCEPT's translation is the development of in vivo dosimetry techniques which can measure both the direct ion dose and the dose due to neutron capture. In this work, we report signatures which can be used to discriminate between photons resulting from neutron capture and those originating from other processes. A Geant4 Monte Carlo simulation study into timing and energy thresholds for discrimination of prompt gamma photons resulting from thermal neutron capture during NCEPT was conducted. Three simulated [Formula: see text] mm[Formula: see text] cubic PMMA targets were irradiated by [Formula: see text]He or [Formula: see text]C ion beams with a spread out Bragg peak (SOBP) depth range of 60 mm; one target is homogeneous while the others include [Formula: see text] mm[Formula: see text] neutron capture inserts (NCIs) of pure [Formula: see text]B or [Formula: see text]Gd located at the distal edge of the SOBP. The arrival times of photons and neutrons entering a simulated [Formula: see text] mm[Formula: see text] ideal detector were recorded. A temporal mask of 50-60 ns was found to be optimal for maximising the discrimination of the photons resulting from the neutron capture by boron and gadolinium. A range of candidate detector and thermal neutron shielding materials were simulated, and detections meeting the proposed acceptance criteria (i.e. falling within the target energy window and arriving 60 ns post beam-off) were classified as true or false positives, depending on their origin. The ratio of true/false positives ([Formula: see text]) was calculated; for targets with [Formula: see text]B and [Formula: see text]Gd NCIs, the detector materials which resulted in the highest [Formula: see text] were cadmium-shielded CdTe and boron-shielded LSO, respectively. The optimal irradiation period for both carbon and helium ions was 1 µs for the [Formula: see text]B NCI and 1 ms for the [Formula: see text]Gd NCI.

Applications of Machine Learning to Improve the Clinical Viability of Compton Camera Based in vivo Range Verification in Proton Radiotherapy

We studied the application of a deep, fully connected Neural Network (NN) to process prompt gamma (PG) data measured by a Compton camera (CC) during the delivery of clinical proton radiotherapy beams. The network identifies 1) recorded "bad" PG events arising from background noise during the measurement, and 2) the correct ordering of PG interactions in the CC to help improve the fidelity of "good" data used for image reconstruction. PG emission from a tissue-equivalent target during irradiation with a 150 MeV proton beam delivered at clinical dose rates was measured with a prototype CC. Images were reconstructed from both the raw measured data and the measured data that was further processed with a neural network (NN) trained to identify "good" and "bad" PG events and predict the ordering of individual interactions within the good PG events. We determine if NN processing of the CC data could improve the reconstructed PG images to a level in which they could provide clinically useful information about the in vivo range and range shifts of the proton beams delivered at full clinical dose rates. Results showed that a deep, fully connected NN improved the achievable contrast to noise ratio (CNR) in our images by more than a factor of 8x. This allowed the path, range, and lateral width of the clinical proton beam within a tissue equivalent target to easily be identified from the PG images, even at the highest dose rates of a 150 MeV proton beam used for clinical treatments. On average, shifts in the beam range as small as 3 mm could be identified. However, when limited by the amount of PG data measured with our prototype CC during the delivery of a single proton pencil beam (~1 × 109 protons), the uncertainty in the reconstructed PG images limited the identification of range shift to ~5 mm. Substantial improvements in CC images were obtained during clinical beam delivery through NN pre-processing of the measured PG data. We believe this shows the potential of NNs to help improve and push CC-based PG imaging toward eventual clinical application for proton RT treatment delivery verification.

A novel proton-integrating radiography system design using a monolithic scintillator detector: experimental studies

Research on proton-based imaging systems aims to improve treatment planning, internal anatomy visualization, and patient alignment for proton radiotherapy. The purpose of this study was to demonstrate a new proton radiography system design consisting of a monolithic plastic scintillator volume and two optical cameras for use with scanning proton pencil beams. Unlike the thin scintillating plates currently used for proton radiography, the plastic scintillator volume (20 × 20 × 20 cm3) captures a wider distribution of proton beam energy depositions and avoids proton-beam modulation. The proton imaging system's characteristics were tested using image uniformity (2.6% over a 5 × 5 cm2 area), stability (0.37%), and linearity (R2 = 1) studies. We used the light distribution produced within the plastic scintillator to generate proton radiographs via two different approaches: (a) integrating light by using a camera placed along the beam axis, and (b) capturing changes to the proton Bragg peak positions with a camera placed perpendicularly to the beam axis. The latter method was used to plot and evaluate relative shifts in percentage depth light (PDL) profiles of proton beams with and without a phantom in the beam path. A curvelet minimization algorithm used differences in PDL profiles to reconstruct and refine the phantom water-equivalent thickness (WET) map. Gammex phantoms were used to compare the proton radiographs generated by these two methods. The relative accuracies in calculating WET of the phantoms using the calibration-based beam-integration (and the PDL) methods were -0.18 ± 0.35% (-0.29 ± 3.11%), -0.11 ± 0.51% (-0.15 ± 2.64%), -2.94 ± 1.20% (-0.75 ± 6.11%), and -1.65 ± 0.35% (0.36 ± 3.93%) for solid water, adipose, cortical bone, and PMMA, respectively. Further exploration of this unique multicamera-based imaging system is warranted and could lead to clinical applications that improve treatment planning and patient alignment for proton radiotherapy.

The effects of Compton camera data acquisition and readout timing on PG imaging for proton range verification

The purpose of this study was to determine how the characteristics of the data acquisition (DAQ) electronics of a Compton camera (CC) affect the quality of the recorded prompt gamma (PG) interaction data and the reconstructed images, during clinical proton beam delivery. We used the Monte-Carlo-plus-Detector-Effect (MCDE) model to simulate the delivery of a 150 MeV clinical proton pencil beam to a tissue-equivalent plastic phantom. With the MCDE model we analyzed how the recorded PG interaction data changed as two characteristics of the DAQ electronics of a CC were changed: (1) the number of data readout channels; and (2) the active charge collection, readout, and reset time. As the proton beam dose rate increased, the number of recorded PG single-, double-, and triple-scatter events decreased by a factor of 60× for the current DAQ configuration of the CC. However, as the DAQ readout channels were increased and the readout/reset timing decreased, the number of recorded events decreased by <5× at the highest clinical dose rate. The increased number of readout channels and reduced readout/reset timing also resulted in higher quality recorded data. That is, a higher percentage of the recorded double- and triple-scatters were "true" events (caused by a single incident gamma) and not "false" events (caused by multiple incident gammas). The increase in the number and the quality of recorded data allowed higher quality PG images to be reconstructed even at the highest clinical dose rates.

Dictionary-based software for proton dose reconstruction and submilimetric range verification

Objective. This paper presents a new method for fast reconstruction (compatible with in-beam use) of deposited dose during proton therapy using data acquired from a PET scanner. The most innovative feature of this novel method is the production of noiseless reconstructed dose distributions from which proton range can be derived with high precision. Approach. A new MLEM & simulated annealing (MSA) algorithm, developed especially in this work, reconstructs the deposited dose distribution from a realistic pre-calculated activity-dose dictionary. This dictionary contains the contribution of each beam in the plan to the 3D activity and dose maps, as calculated by a Monte Carlo simulation. The MSA algorithm, using a priori information of the treatment plan, seeks for the linear combination of activities of the precomputed beams that best fits the observed PET data, obtaining at the same time the deposited dose. Main results. the method has been tested using simulated data to determine its performance under 4 different test cases: (1) dependency of range detection accuracy with delivered dose, (2) in-beam versus offline verification, (3) ability to detect anatomical changes and (4) reconstruction of a realistic spread-out Bragg peak. The results show the ability of the method to accurately reconstruct doses from PET data corresponding to 1 Gy irradiations, both in intra-fraction and inter-fraction verification scenarios. For this dose level (1 Gy) the method was able to spot range variations as small as 0.6 mm. Significance. out method is able to reconstruct dose maps with remarkable accuracy from clinically relevant dose levels down to 1 Gy. Furthermore, due to the noiseless nature of reconstructed dose maps, an accuracy better than one millimeter was obtained in proton range estimates. These features make of this method a realistic option for range verification in proton therapy.

Proton Therapy for Prostate Cancer: Challenges and Opportunities

Simple Summary

Reported clinical outcomes of proton therapy (PT) for localized prostate cancer are similar to photon-based external beam radiotherapy. Apparently, the dosimetric advantages of PT have yet to be translated to clinical benefits. The suboptimal clinical outcomes of PT might be attributable to inadequate dose prescription, as indicated by the ASCENDE-RT trial. Moreover, uncertainties involved in the treatment planning and delivery processes, as well as technological limitations in PT treatment systems, may lead to discrepancies between planned doses and actual doses delivered to patients. In this article, we reviewed the current status of PT for prostate cancer and discussed different clinical implementations that could potentially improve the clinical outcome of PT for prostate cancer. Various technological advancements under which uncertainties in dose calculations can be minimized, including MRI-guided PT, dual-energy photon-counting CT and high-resolution Monte Carlo-based treatment planning systems, are highlighted.

Abstract

The dosimetric advantages of proton therapy (PT) treatment plans are demonstrably superior to photon-based external beam radiotherapy (EBRT) for localized prostate cancer, but the reported clinical outcomes are similar. This may be due to inadequate dose prescription, especially in high-risk disease, as indicated by the ASCENDE-RT trial. Alternatively, the lack of clinical benefits with PT may be attributable to improper dose delivery, mainly due to geometric and dosimetric uncertainties during treatment planning, as well as delivery procedures that compromise the dose conformity of treatments. Advanced high-precision PT technologies, and treatment planning and beam delivery techniques are being developed to address these uncertainties. For instance, external magnetic resonance imaging (MRI)-guided patient setup rooms are being developed to improve the accuracy of patient positioning for treatment. In-room MRI-guided patient positioning systems are also being investigated to improve the geometric accuracy of PT. Soon, high-dose rate beam delivery systems will shorten beam delivery time to within one breath hold, minimizing the effects of organ motion and patient movements. Dual-energy photon-counting computed tomography and high-resolution Monte Carlo-based treatment planning systems are available to minimize uncertainties in dose planning calculations. Advanced in-room treatment verification tools such as prompt gamma detector systems will be used to verify the depth of PT. Clinical implementation of these new technologies is expected to improve the accuracy and dose conformity of PT in the treatment of localized prostate cancers, and lead to better clinical outcomes. Improvement in dose conformity may also facilitate dose escalation, improving local control and implementation of hypofractionation treatment schemes to improve patient throughput and make PT more cost effective.