Characterization of a high-resolution 2D transmission ion chamber for independent validation of proton pencil beam scanning of conventional and FLASH dose delivery

Electron beam response corrections for an ultra-high-dose-rate capable diode dosimeter

BACKGROUND

Ultra-high-dose-rate (UHDR) electron beams have been commonly utilized in FLASH studies and the translation of FLASH Radiotherapy (RT) to the clinic. The EDGE diode detector has potential use for UHDR dosimetry albeit with a beam energy dependency observed.

PURPOSE

The purpose is to present the electron beam response for an EDGE detector in dependence on beam energy, to characterize the EDGE detector's response under UHDR conditions, and to validate correction factors derived from the first detailed Monte Carlo model of the EDGE diode against measurements, particularly under UHDR conditions.

METHODS

Percentage depth doses (PDDs) for the UHDR Mobetron were measured with both EDGE detectors and films. A detailed Monte Carlo (MC) model of the EDGE detector has been configured according to the blueprint provided by the manufacturer under an NDA agreement. Water/silicon dose ratios of EDGE detector for a series of mono-energetic electron beams have been calculated. The dependence of the water/silicon dose ratio on depth for a FLASH relevant electron beam was also studied. An analytical approach for the correction of PDD measured with EDGE detectors was established.

RESULTS

Water/silicon dose ratio decreased with decreasing electron beam energy. For the Mobetron 9 MeV UHDR electron beam, the ratio decreased from 1.09 to 1.03 in the build-up region, maintained in range of 0.98-1.02 at the fall-off region and raised to a plateau in value of 1.08 at the tail. By applying the corrections, good agreement between the PDDs measured by the EDGE detector and those measured with film was achieved.

CONCLUSIONS

Electron beam response of an UHDR capable EDGE detector was derived from first principles utilizing a sophisticated MC model. An analytical approach was validated for the PDDs of UHDR electron beams. The results demonstrated the capability of EDGE detector in measuring PDDs of UHDR electron beams.

Beam monitor chamber calibration of a synchro-cyclotron high dose rate per pulse pulsed scanned proton beam

Objective. Ionization chambers, mostly used for beam calibration and for reference dosimetry, can show high recombination effects in pulsed high dose rate proton beams. The aims of this paper are: first, to characterize the linearity response of newly designed asymmetrical beam monitor chambers (ABMC) in a 100-226 MeV pulsed high dose rate per pulse scanned proton beam; and secondly, to calibrate the ABMC with a PPC05 (IBA Dosimetry) plane parallel ionization chamber and compare to calibration with a home-made Faraday cup (FC).Approach. The ABMC response linearity was evaluated with both the FC and a PTW 60019 microDiamond detector. Regarding ionometry-based ABMC calibration, recombination factors were evaluated theoretically, then numerically, and finally experimentally measured in water for a plane parallel ionization chamber PPC05 (IBA Dosimetry) throughkssaturation curves. Finally, ABMC calibration was also achieved with FC and compared to the ionometry method for 7 energies.Main results. Linearity measurements showed that recombination losses in the new ABMC design were well taken into account for the whole range of the machine dose rates. The two-voltage-method was not suitable for recombination correction, but Jaffé's plots analysis was needed, emphasizing the current IAEA TRS-398 reference protocol limitations. Concerning ABMC calibration, FC based absorbed dose estimation and PPC05-based absorbed dose estimation differ by less than 6.3% for the investigated energies.Significance.So far, no update on reference dosimetry protocols is available to estimate the absorbed dose in ionization chambers for clinical high dose rate per pulse pulsed scanned proton beams. This work proposes a validation of the new ABMC design, a method to take into account the recombination effect for ionometry-based ABMC calibration and a comparison with FC dose estimation in this type of proton beams.

Feasibility of Synchrotron-Based Ultra-High Dose Rate (UHDR) Proton Irradiation with Pencil Beam Scanning for FLASH Research

BACKGROUND

This study aims to present the feasibility of developing a synchrotron-based proton ultra-high dose rate (UHDR) pencil beam scanning (PBS) system.

METHODS

The RF extraction power in the synchrotron system was increased to generate 142.4 MeV pulsed proton beams for UHDR irradiation at ~100 nA beam current. The charge per spill was measured using a Faraday cup. The spill length and microscopic time structure of each spill was measured with a 2D strip transmission ion chamber. The measured UHDR beam fluence was used to derive the spot dwell time for pencil beam scanning. Absolute dose distributions at various depths and spot spacings were measured using Gafchromic films in a solid-water phantom.

RESULTS

For proton UHDR beams at 142.4 MeV, the maximum charge per spill is 4.96 ± 0.10 nC with a maximum spill length of 50 ms. This translates to an average beam current of approximately 100 nA during each spill. Using a 2 × 2 spot delivery pattern, the delivered dose per spill at 5 cm and 13.5 cm depth is 36.3 Gy (726.3 Gy/s) and 56.2 Gy (1124.0 Gy/s), respectively.

CONCLUSIONS

The synchrotron-based proton therapy system has the capability to deliver pulsed proton UHDR PBS beams. The maximum deliverable dose and field size per pulse are limited by the spill length and extraction charge.

FLASH dose rate calculation based on log files in proton pencil beam scanning therapy

BACKGROUND

In radiation therapy, irradiating healthy normal tissues in the beam trajectories is inevitable. This unnecessary dose means that patients undergoing treatment risk developing side effects. Recently, FLASH radiotherapy delivering ultra-high-dose-rate beams has been re-examined because of its normal-tissue-sparing effect. To confirm the mean and instantaneous dose rates of the FLASH beam, stable and accurate dosimetry is required.

PURPOSE

Detailed verification of the FLASH effect requires dosimeters and a method to measure the average and instantaneous dose rate stably for 2- or 3-dimensional dose distributions. To verify the delivered FLASH beam, we utilized machine log files from the built-in monitor chamber to develop a dosimetry method to calculate the dose and average/instantaneous dose rate distributions in two or three dimensions in a phantom.

METHODS

To create a spread-out Bragg peak (SOBP) and deliver a uniform dose in a target, a mini-ridge filter was created with a 3D printer. Proton pencil beam line scanning plans of 2 × 2 cm2 , 3 × 3 cm2 , 4 × 4 cm2 , and round shapes with 2.3 cm diameter patterns delivering 230 MeV energy protons were created. The absorbed dose in the solid water phantom of each plan was measured using a PPC05 ionization chamber (IBA Dosimetry, Virginia, USA) in the SOBP region, and the log files for each plan were exported from the treatment control system console. Using these log files, the delivered dose and average dose rate were calculated using two methods: a direct method and a Monte Carlo (MC) simulation method that uses log file information. The computed and average dose rates were compared with the ionization chamber measurements. Additionally, instantaneous dose rates in user-defined volumes were calculated using the MC simulation method with a temporal resolution of 5 ms.

RESULTS

Compared to ionization chamber dosimetry, 10 of 12 cases using the direct calculation method and 9 of 11 cases using the MC method had a dose difference below ±3%. Nine of 12 cases using the direct calculation method and 8 of 11 cases using the MC method had dose rate differences below ±3%. The average and maximum dose differences for the direct calculation and MC method were-0.17, +0.72%, and -3.15, +3.32%, respectively. For the dose rate difference, the average and maximum for the direct calculation and MC method were +1.26, +1.12%, and +3.75, +3.15%, respectively. In the instantaneous dose rate calculation with the MC simulation, a large fluctuation with a maximum of 163 Gy/s and a minimum of 4.29 Gy/s instantaneous dose rate was observed in a specific position, whereas the mean dose rate was 62 Gy/s.

CONCLUSIONS

We successfully developed methods in which machine log files are used to calculate the dose and the average and instantaneous dose rates for FLASH radiotherapy and demonstrated the feasibility of verifying the delivered FLASH beams.

Framework for Quality Assurance of Ultrahigh Dose Rate Clinical Trials Investigating FLASH Effects and Current Technology Gaps

FLASH radiation therapy (FLASH-RT), delivered with ultrahigh dose rate (UHDR), may allow patients to be treated with less normal tissue toxicity for a given tumor dose compared with currently used conventional dose rate. Clinical trials are being carried out and are needed to test whether this improved therapeutic ratio can be achieved clinically. During the clinical trials, quality assurance and credentialing of equipment and participating sites, particularly pertaining to UHDR-specific aspects, will be crucial for the validity of the outcomes of such trials. This report represents an initial framework proposed by the NRG Oncology Center for Innovation in Radiation Oncology FLASH working group on quality assurance of potential UHDR clinical trials and reviews current technology gaps to overcome. An important but separate consideration is the appropriate design of trials to most effectively answer clinical and scientific questions about FLASH. This paper begins with an overview of UHDR RT delivery methods. UHDR beam delivery parameters are then covered, with a focus on electron and proton modalities. The definition and control of safe UHDR beam delivery and current and needed dosimetry technologies are reviewed and discussed. System and site credentialing for large, multi-institution trials are reviewed. Quality assurance is then discussed, and new requirements are presented for treatment system standard analysis, patient positioning, and treatment planning. The tables and figures in this paper are meant to serve as reference points as we move toward FLASH-RT clinical trial performance. Some major questions regarding FLASH-RT are discussed, and next steps in this field are proposed. FLASH-RT has potential but is associated with significant risks and complexities. We need to redefine optimization to focus not only on the dose but also on the dose rate in a manner that is robust and understandable and that can be prescribed, validated, and confirmed in real time. Robust patient safety systems and access to treatment data will be critical as FLASH-RT moves into the clinical trials.

Development of a Real-Time Pixel Array-Type Detector for Ultrahigh Dose-Rate Beams

Although research into ultrahigh dose-rate (UHDR) radiation therapy is ongoing, there is a significant lack of experimental measurements for two-dimensional (2D) dose-rate distributions. Additionally, conventional pixel-type detectors result in significant beam loss. In this study, we developed a pixel array-type detector with adjustable gaps and a data acquisition system to evaluate its effectiveness in measuring UHDR proton beams in real time. We measured a UHDR beam at the Korea Institute of Radiological and Medical Sciences using an MC-50 cyclotron, which produced a 45-MeV energy beam with a current range of 10-70 nA, to confirm the UHDR beam conditions. To minimize beam loss during measurement, we adjusted the gap and high voltage on the detector and determined the collection efficiency of the developed detector through Monte Carlo simulation and experimental measurements of the 2D dose-rate distribution. We also verified the accuracy of the real-time position measurement using the developed detector with a 226.29-MeV PBS beam at the National Cancer Center of the Republic of Korea. Our results indicate that, for a current of 70 nA with an energy beam of 45 MeV generated using the MC-50 cyclotron, the dose rate exceeded 300 Gy/s at the center of the beam, indicating UHDR conditions. Simulation and experimental measurements show that fixing the gap at 2 mm and the high voltage at 1000 V resulted in a less than 1% loss of collection efficiency when measuring UHDR beams. Furthermore, we achieved real-time measurements of the beam position with an accuracy of within 2% at five reference points. In conclusion, our study developed a beam monitoring system that can measure UHDR proton beams and confirmed the accuracy of the beam position and profile through real-time data transmission.

Failure Mode and Effects Analysis for Experimental Use of FLASH on a Clinical Accelerator

PURPOSE

The use of a linear accelerator (LINAC) in ultrahigh-dose-rate (UHDR) mode can provide a conduit for wider access to UHDR FLASH effects, sparing normal tissue, but care needs to be taken in the use of such systems to ensure errors are minimized. The failure mode and effects analysis was carried out in a team that has been involved in converting a LINAC between clinical use and UHDR experimental mode for more than 1 year after the proposed methods of TG100.

METHODS AND MATERIALS

A team of 9 professionals with extensive experience were polled to outline the process map and workflow for analysis, and developed fault trees for potential errors, as well as failure modes that would result. The team scored the categories of severity magnitude, occurrence likelihood, and detectability potential in a scale of 1 to 10, so that a risk priority number (RPN = severity×occurrence×detectability) could be assessed for each.

RESULTS

A total of 46 potential failure modes were identified, including 5 with an RPN >100. These failure modes involved (1) patient set up, (2) gating mechanisms in delivery, and (3) detector in the beam stop mechanism. The identified methods to mitigate errors included the (1) use of a checklist post conversion, (2) use of robust radiation detectors, (3) automation of quality assurance and beam consistency checks, and (4) implementation of surface guidance during beam delivery.

CONCLUSIONS

The failure mode and effects analysis process was considered critically important in this setting of a new use of a LINAC, and the expert team developed a higher level of confidence in the ability to safely move UHDR LINAC use toward expanded research access.

An Integrated Physical Optimization Framework for Proton Stereotactic Body Radiation Therapy FLASH Treatment Planning Allows Dose, Dose Rate, and Linear Energy Transfer Optimization Using Patient-Specific Ridge Filters

PURPOSE

Patient-specific ridge filters provide a passive means to modulate proton energy to obtain a conformal dose. Here we describe a new framework for optimization of filter design and spot maps to meet the unique demands of ultrahigh-dose-rate (FLASH) radiation therapy. We demonstrate an integrated physical optimization Intensity-modulated proton therapy (IMPT) (IPO-IMPT) approach for optimization of dose, dose-averaged dose rate (DADR), and dose-averaged linear energy transfer (LET_d).

METHODS AND MATERIALS

We developed an inverse planning software to design patient-specific ridge filters that spread the Bragg peak from a fixed-energy, 250-MeV beam to a proximal beam-specific planning target volume. The software defines patient-specific ridge filter pin shapes and uses a Monte Carlo calculation engine, based on Geant4, to provide dose and LET influence matrices. Plan optimization, using matRAD, accommodates the IPO-IMPT objective function considering dose, dose rate, and LET simultaneously with minimum monitor unit constraints. The framework enables design of both regularly spaced and sparse-optimized ridge filters, from which some pins are omitted to allow faster delivery and selective LET optimization. To demonstrate the framework, we designed ridge filters for 3 example patients with lung cancer and optimized the plans using IPO-IMPT.

RESULTS

The IPO-IMPT framework selectively spared the organs at risk by reducing LET and increasing dose rate, relative to IMPT planning. Sparse-optimized ridge filters were superior to regularly spaced ridge filters in dose rate. Depending on which parameter is prioritized, volume distributions and histograms for dose, DADR, and LET_d, using evaluation structures specific to heart, lung, and esophagus, show high levels of FLASH dose-rate coverage and/or reduced LET_d, while maintaining dose coverage within the beam specific planning target volume.

CONCLUSIONS

This proof-of-concept study demonstrates the feasibility of using an IPO-IMPT framework to accomplish proton FLASH stereotactic body proton therapy, accounting for dose, DADR, and LET_d simultaneously.

A phenomenological model of proton FLASH oxygen depletion effects depending on tissue vasculature and oxygen supply

Introduction

Radiation-induced oxygen depletion in tissue is assumed as a contributor to the FLASH sparing effects. In this study, we simulated the heterogeneous oxygen depletion in the tissue surrounding the vessels and calculated the proton FLASH effective-dose-modifying factor (FEDMF), which could be used for biology-based treatment planning.

Methods

The dose and dose-weighted linear energy transfer (LET) of a small animal proton irradiator was simulated with Monte Carlo simulation. We deployed a parabolic partial differential equation to account for the generalized radiation oxygen depletion, tissue oxygen diffusion, and metabolic processes to investigate oxygen distribution in 1D, 2D, and 3D solution space. Dose and dose rates, particle LET, vasculature spacing, and blood oxygen supplies were considered. Using a similar framework for the hypoxic reduction factor (HRF) developed previously, the FEDMF was derived as the ratio of the cumulative normoxic-equivalent dose (CNED) between CONV and UHDR deliveries.

Results

Dynamic equilibrium between oxygen diffusion and tissue metabolism can result in tissue hypoxia. The hypoxic region displayed enhanced radio-resistance and resulted in lower CNED under UHDR deliveries. In 1D solution, comparing 15 Gy proton dose delivered at CONV 0.5 and UHDR 125 Gy/s, 61.5% of the tissue exhibited ≥20% FEDMF at 175 μm vasculature spacing and 18.9 μM boundary condition. This percentage reduced to 34.5% and 0% for 8 and 2 Gy deliveries, respectively. Similar trends were observed in the 3D solution space. The FLASH versus CONV differential effect remained at larger vasculature spacings. A higher FLASH dose rate showed an increased region with ≥20% FEDMF. A higher LET near the proton Bragg peak region did not appear to alter the FLASH effect.

Conclusion

We developed 1D, 2D, and 3D oxygen depletion simulation process to obtain the dynamic HRF and derive the proton FEDMF related to the dose delivery parameters and the local tissue vasculature information. The phenomenological model can be used to simulate or predict FLASH effects based on tissue vasculature and oxygen concentration data obtained from other experiments.

A potential revolution in cancer treatment: A topical review of FLASH radiotherapy

FLASH radiotherapy (RT) is a novel technique in which the ultrahigh dose rate (UHDR) (≥40 Gy/s) is delivered to the entire treatment volume. Recent outcomes of in vivo studies show that the UHDR RT has the potential to spare normal tissue without sacrificing tumor control. There is a growing interest in the application of FLASH RT, and the ultrahigh dose irradiation delivery has been achieved by a few experimental and modified linear accelerators. The underlying mechanism of FLASH effect is yet to be fully understood, but the oxygen depletion in normal tissue providing extra protection during FLASH irradiation is a hypothesis that attracts most attention currently. Monte Carlo simulation is playing an important role in FLASH, enabling the understanding of its dosimetry calculations and hardware design. More advanced Monte Carlo simulation tools are under development to fulfill the challenge of reproducing the radiolysis and radiobiology processes in FLASH irradiation. FLASH RT may become one of standard treatment modalities for tumor treatment in the future. This paper presents the history and status of FLASH RT studies with a focus on FLASH irradiation delivery modalities, underlying mechanism of FLASH effect, in vivo and vitro experiments, and simulation studies. Existing challenges and prospects of this novel technique are discussed in this manuscript.

Advanced pencil beam scanning Bragg peak FLASH-RT delivery technique can enhance lung cancer planning treatment outcomes compared to conventional multiple-energy proton PBS techniques

PURPOSE

To investigate the dosimetric characteristics between an advanced proton pencil beam scanning (PBS) Bragg peak FLASH technique and conventional PBS planning technique in lung tumors. To evaluate the "FLASHness" of single-field in a multiple-field delivery scheme for a hypofractionation regimen and move a step forward to clinical application.

METHODS

Single-energy PBS Bragg peak FLASH treatment plans were optimized based on a novel Bragg peak tracking technique to enable Bragg peaks to stop at the distal edge of the target. Inverse treatment planning using multiple-field optimization (MFO) can achieve sufficient FLASH dose rate and intensity-modulated proton therapy (IMPT)-equivalent dosimetric quality. The dose rate of organs-at-risk (OARs) and the target were calculated under FLASH machine parameters. A group of 10 consecutive lung SBRT patients was optimized to 34 Gy/fraction using a standard treatment of PBS technique with multiple energy layers as references to the Bragg peak plans. The dosimetric quality was compared between Bragg peak FLASH and conventional plans based on RTOG0915 dose metrics. FLASH dose rate ratios (V_40Gy/s) were calculated as a metric of the FLASH-sparing effect.

RESULTS

For higher dose thresholds, the Bragg peak plans achieved greater V_40Gy/s FLASH coverage for all major OARs. The V40Gy/s was close to 100% for all OARs when the dose thresholds were > 5 Gy for full plan and single beam evaluations. The less "FLASHness" region was associated with a low dose distribution, mainly occurring in the PBS field penumbra region. The conventional IMPT treatment plans yielded slightly superior target dose uniformity with a D_2%(%) of 108.02% versus that of Bragg peak 300 MU plans of 111.81% (p < 0.01) and that of Bragg peak 1200 MU plans of 115.95% (p < 0.01). No significant difference in dose metrics was found between Bragg peak and IMPT treatment plans for the spinal cord, esophagus, heart, or lung-GTV (all p > 0.05).

CONCLUSION

Hypofractionated lung Bragg peak plans can maintain comparable plan quality to conventional PBS while achieving sufficient FLASH dose rate coverage for major OARs for each field under the multiple-field delivery scheme. The novel Bragg peak FLASH technique has the potential to enhance lung cancer planning treatment outcomes compared to standard PBS treatment techniques.

A high spatiotemporal resolution 2D strip ionization chamber array for proton pencil beam scanning FLASH radiotherapy

PURPOSE

Experimental measurements of 2D dose rate distributions in proton pencil beam scanning (PBS) FLASH radiation therapy (RT) are currently lacking. In this study, we characterize a newly designed 2D strip-segmented ionization chamber array (SICA) with high spatial and temporal resolution and demonstrate its applications in a modern proton PBS delivery system at both conventional and ultra-high dose rates.

METHODS

A dedicated research beamline of the Varian ProBeam system was employed to deliver a 250 MeV proton PBS beam with nozzle currents up to 215 nA. In the research and clinical beamlines, the spatial, temporal, and dosimetric performance of the SICA was characterized and compared with measurements using parallel-plate ion chambers (IBA PPC05 and PTW Advanced Markus chamber), a 2D scintillator camera (IBA Lynx), Gafchromic films (EBT-XD), and a Faraday Cup. A novel reconstruction approach was proposed to enable the measurement of 2D dose and dose rate distributions using such a strip-type detector.

RESULTS

The SICA demonstrated a position accuracy of 0.12 ± 0.02 mm at a 20 kHz sampling rate (50 μs per event) and a linearity of R² > 0.99 for both dose and dose rate with nozzle beam currents ranging from 1 nA to 215 nA. The 2D dose comparison to the film measurement resulted in a gamma passing rate of 99.8% (2 mm/2%). A measurement-based proton PBS 2D FLASH dose rate distribution was compared to simulation results and showed a gamma passing rate of 97.3% (2 mm/2%).

CONCLUSIONS

The newly designed SICA demonstrated excellent spatial, temporal, and dosimetric performance and is well suited for commissioning, quality assurance (QA), and a wide range of clinical applications in proton PBS clinical and FLASH radiotherapy. This article is protected by copyright. All rights reserved.

The current status of preclinical proton FLASH radiation and future directions

We review the current status of proton FLASH experimental systems, including preclinical physical and biological results. Technological limitations on preclinical investigation of FLASH biological mechanisms and determination of clinically relevant parameters are discussed. A review of the biological data reveals no reproduced proton FLASH effect in vitro and a significant in vivo FLASH sparing effect of normal tissue toxicity observed with multiple proton FLASH irradiation systems. Importantly, multiple studies suggest little or no difference in tumor growth delay for proton FLASH when compared to conventional dose rate proton radiation. A discussion follows on future areas of development with a focus on the determination of the optimal parameters for maximizing the therapeutic ratio between tumor and normal tissue response and ultimately clinical translation of proton FLASH radiation.