Development and characterization of a dedicated dose monitor for ultrahigh-dose-rate scanned carbon-ion beams.

The current monochromatic beam mode (i.e., uHDR irradiation mode) of the scanned carbon-ion beam lacks a dedicated dose monitor, making the beam control challenging. We developed and characterized a dedicated dose monitor for uHDR-scanned carbon-ion beams. Furthermore, a simple measurable dose rate (dose rate per spot (DRspot)) was suggested by using the developed dose monitor and experimentally validating quantities relevant to the uHDR scanned carbon-ion beam. A large plane-parallel ionization chamber (IC) with a smaller electrode spacing was used to reduce uHDR recombination effects, and a dedicated operational amplifier was manufactured for the uHDR-scanned carbon-ion beam. The dose linearity of the IC was within ± 1% in the range of 1.8-12.3 Gy. The spatial inhomogeneity of the dose response of the IC was ± 0.38% inside the ± 40-mm detector area, and a systematic deviation of approximately 2% was measured at the edge of the detector. uHDR irradiation with beam scanning was tested and verified for different doses at the corresponding dose rates (in terms of both the average dose rate and DRspot). We confirmed that the dose monitor can highlight the characteristics (i.e., dose, dose rate, and dose profile) of uHDR-scanned carbon-ion beams at several dose levels in the monochromatic beam mode.

Luminescence imaging of water irradiated by protons under FLASH radiation therapy conditions.

Objective.FLASH radiation therapy with ultrahigh dose rates (UHDR) has the potential to reduce damage to normal tissue while maintaining anti-tumor efficacy. However, rapid and precise dose distribution measurements remain difficult for FLASH radiation therapy with proton beams. To solve this problem, we performed luminescence imaging of water following irradiation by a UHDR proton beam captured using a charge-coupled device camera.Approach. We used 60 MeV proton beams with dose rates of 0.03-837 Gy s-1from a cyclotron. Therapeutic 139.3 MeV proton beams with dose rates of 0.45-4320 Gy s-1delivered by a synchrotron-based proton therapy system were also tested. The luminescent light intensity induced by the UHDR beams was compared with that produced by conventional beams to compare the dose rate dependency of the light intensity and its profile.Main results. Luminescence images of water were clearly visualized under UHDR conditions, with significantly shorter exposure times than those with conventional beams. The light intensity was linearly proportional to the delivered dose, which is similar to that of conventional beams. No significant dose-rate dependency was observed for 0.03-837 Gy s-1. The light-intensity profiles of the UHDR beams agreed with those of conventional beams. The results did not differ between accelerators (synchrotron or cyclotron) and beam energies.Significance. Luminescence imaging of water is achievable with UHDR proton beams as well as with conventional beams. The proposed method should be suitable for rapid and easy quality assurance investigations for proton FLASH therapy, because it facilitates real-time, filmless measurements of dose distributions, and is useful for rapid feedback.

Ultra-high Dose-rate Carbon-ion Scanning Beam With a Compact Medical Synchrotron Contributing to Further Development of FLASH Irradiation.

BACKGROUND/AIM: The focus of this report is establishing an irradiation arrangement to realize an ultra-high dose-rate (uHDR; FLASH) of scanned carbon-ion irradiation possible with a compact commonly available medical synchrotron. MATERIALS AND METHODS: Following adjustments to the operation it became possible to extract ≥1.0×109 carbon ions at 208.3 MeV/u (86 mm in range) per 100 ms. The design takes the utmost care to prevent damage to monitors, particularly in the nozzle, achieved by the uHDR beam not passing through this part of the apparatus. Doses were adjusted by extraction times, using a function generator. After one scan by the carbon-ion beam it became possible to create a field within the extraction time. The Advanced Markus chamber (AMC) and Gafchromic film are then able to measure the absolute dose and field size at a plateau depth, with the operating voltage of the chamber at 400 V at the uHDR for the AMC. RESULTS: The beam scanning utilizing this uHDR irradiation could be confirmed at a dose of 6.5±0.08 Gy (±3% homogeneous) at this volume over at least 16×16 mm2 corresponding to a dose-rate of 92.3 Gy/s (±1.3%). The dose was ca. 0.7, 1.5, 2.9, and 5.4 Gy depending on dose-rate and field size, with the rate of killed cells increasing with the irradiation dose. CONCLUSION: The compact medical synchrotron achieved FLASH dose-rates of >40 Gy/s at different dose levels and in useful field sizes for research with the apparatus and arrangement developed here.

Adaptation and dosimetric commissioning of a synchrotron-based proton beamline for FLASH experiments.

Objective. Irradiation with ultra-high dose rates (>40 Gy s-1), also known as FLASH irradiation, has the potential to shift the paradigm of radiation therapy because of its reduced toxicity to normal tissues compared to that of conventional irradiations. The goal of this study was to (1) achieve FLASH irradiation conditions suitable for pre-clinicalin vitroandin vivobiology experiments using our synchrotron-based proton beamline and (2) commission the FLASH irradiation conditions achieved.Approach. To achieve these suitable FLASH conditions, we made a series of adaptations to our proton beamline, including modifying the spill length and size of accelerating cycles, repurposing the reference monitor for dose control, and expanding the field size with a custom double-scattering system. We performed the dosimetric commissioning with measurements using an Advanced Markus chamber and EBT-XD films as well as with Monte Carlo simulations.Main results. Through adaptations, we have successfully achieved FLASH irradiation conditions, with an average dose rate of up to 375 Gy s-1. The Advanced Markus chamber was shown to be appropriate for absolute dose calibration under our FLASH conditions with a recombination factor ranging from 1.002 to 1.006 because of the continuous nature of our synchrotron-based proton delivery within a spill. Additionally, the absolute dose measured using the Advanced Markus chamber and EBT-XD films agreed well, with average and maximum differences of 0.32% and 1.63%, respectively. We also performed a comprehensive temporal analysis for FLASH spills produced by our system, which helped us identify a unique relationship between the average dose rate and the dose in our FLASH irradiation.Significance.We have established a synchrotron-based proton FLASH irradiation platform with accurate and precise dosimetry that is suitable for pre-clinical biology experiments. The unique time structure of the FLASH irradiation produced by our synchrotron-based system may shed new light onto the mechanism behind the FLASH effect.

Dual ring multilayer ionization chamber and theory-based correction technique for scanning proton therapy.

PURPOSE: To develop a multilayer ionization chamber (MLIC) and a correction technique that suppresses differences between the MLIC and water phantom measurements in order to achieve fast and accurate depth dose measurements in pencil beam scanning proton therapy. METHODS: The authors distinguish between a calibration procedure and an additional correction: 1-the calibration for variations in the air gap thickness and the electrometer gains is addressed without involving measurements in water; 2-the correction is addressed to suppress the difference between depth dose profiles in water and in the MLIC materials due to the nuclear interaction cross sections by a semiempirical model tuned by using measurements in water. In the correction technique, raw MLIC data are obtained for each energy layer and integrated after multiplying them by the correction factor because the correction factor depends on incident energy. The MLIC described here has been designed especially for pencil beam scanning proton therapy. This MLIC is called a dual ring multilayer ionization chamber (DRMLIC). The shape of the electrodes allows the DRMLIC to measure both the percentage depth dose (PDD) and integrated depth dose (IDD) because ionization electrons are collected from inner and outer air gaps independently. RESULTS: IDDs for which the beam energies were 71.6, 120.6, 159, 180.6, and 221.4 MeV were measured and compared with water phantom results. Furthermore, the measured PDDs along the central axis of the proton field with a nominal field size of 10 × 10 cm(2) were compared. The spread out Bragg peak was 20 cm for fields with a range of 30.6 and 3 cm for fields with a range of 6.9 cm. The IDDs measured with the DRMLIC using the correction technique were consistent with those that of the water phantom; except for the beam energy of 71.6 MeV, all of the points satisfied the 1% dose/1 mm distance to agreement criterion of the gamma index. The 71.6 MeV depth dose profile showed slight differences in the shallow region, but 94.5% of the points satisfied the 1%/1 mm criterion. The 90% ranges, defined at the 90% dose position in distal fall off, were in good agreement with those in the water phantom, and the range differences from the water phantom were less than ±0.3 mm. The PDDs measured with the DRMLIC were also consistent with those that of the water phantom; 97% of the points passed the 1%/1 mm criterion. CONCLUSIONS: It was demonstrated that the new correction technique suppresses the difference between the depth dose profiles obtained with the MLIC and those obtained from a water phantom, and a DRMLIC enabling fast measurements of both IDD and PDD was developed. The IDDs and PDDs measured with the DRMLIC and using the correction technique were in good agreement with those that of the water phantom, and it was concluded that the correction technique and DRMLIC are useful for depth dose profile measurements in pencil beam scanning proton therapy.

Optimization and evaluation of multiple gating beam delivery in a synchrotron-based proton beam scanning system using a real-time imaging technique.

PURPOSE: To find the optimum parameter of a new beam control function installed in a synchrotron-based proton therapy system. METHODS: A function enabling multiple gated irradiation in the flat top phase has been installed in a real-time-image gated proton beam therapy (RGPT) system. This function is realized by a waiting timer that monitors the elapsed time from the last gate-off signal in the flat top phase. The gated irradiation efficiency depends on the timer value, Tw. To find the optimum Tw value, gated irradiation efficiency was evaluated for each configurable Tw value. 271 gate signal data sets from 58 patients were used for the simulation. RESULTS: The highest mean efficiency 0.52 was obtained in TW=0.2s. The irradiation efficiency was approximately 21% higher than at TW=0s, which corresponds to ordinary synchrotron operation. The irradiation efficiency was improved in 154 (57%) of the 271 cases. The irradiation efficiency was reduced in 117 cases because the TW value was insufficient or the function introduced an unutilized wait time for the next gate-on signal in the flat top phase. In the actual treatment of a patient with a hepatic tumor at Tw=0.2s, 4.48GyE irradiation was completed within 250s. In contrast, the treatment time of ordinary synchrotron operation was estimated to be 420s. CONCLUSIONS: The results suggest that the multiple gated-irradiation function has potential to improve the gated irradiation efficiency and to reduce the treatment time.

Evaluation of the influence of double and triple Gaussian proton kernel models on accuracy of dose calculations for spot scanning technique.

PURPOSE: The main purpose in this study was to present the results of beam modeling and how the authors systematically investigated the influence of double and triple Gaussian proton kernel models on the accuracy of dose calculations for spot scanning technique. METHODS: The accuracy of calculations was important for treatment planning software (TPS) because the energy, spot position, and absolute dose had to be determined by TPS for the spot scanning technique. The dose distribution was calculated by convolving in-air fluence with the dose kernel. The dose kernel was the in-water 3D dose distribution of an infinitesimal pencil beam and consisted of an integral depth dose (IDD) and a lateral distribution. Accurate modeling of the low-dose region was important for spot scanning technique because the dose distribution was formed by cumulating hundreds or thousands of delivered beams. The authors employed a double Gaussian function as the in-air fluence model of an individual beam. Double and triple Gaussian kernel models were also prepared for comparison. The parameters of the kernel lateral model were derived by fitting a simulated in-water lateral dose profile induced by an infinitesimal proton beam, whose emittance was zero, at various depths using Monte Carlo (MC) simulation. The fitted parameters were interpolated as a function of depth in water and stored as a separate look-up table. These stored parameters for each energy and depth in water were acquired from the look-up table when incorporating them into the TPS. The modeling process for the in-air fluence and IDD was based on the method proposed in the literature. These were derived using MC simulation and measured data. The authors compared the measured and calculated absolute doses at the center of the spread-out Bragg peak (SOBP) under various volumetric irradiation conditions to systematically investigate the influence of the two types of kernel models on the dose calculations. RESULTS: The authors investigated the difference between double and triple Gaussian kernel models. The authors found that the difference between the two studied kernel models appeared at mid-depths and the accuracy of predicting the double Gaussian model deteriorated at the low-dose bump that appeared at mid-depths. When the authors employed the double Gaussian kernel model, the accuracy of calculations for the absolute dose at the center of the SOBP varied with irradiation conditions and the maximum difference was 3.4%. In contrast, the results obtained from calculations with the triple Gaussian kernel model indicated good agreement with the measurements within ±1.1%, regardless of the irradiation conditions. CONCLUSIONS: The difference between the results obtained with the two types of studied kernel models was distinct in the high energy region. The accuracy of calculations with the double Gaussian kernel model varied with the field size and SOBP width because the accuracy of prediction with the double Gaussian model was insufficient at the low-dose bump. The evaluation was only qualitative under limited volumetric irradiation conditions. Further accumulation of measured data would be needed to quantitatively comprehend what influence the double and triple Gaussian kernel models had on the accuracy of dose calculations.

A proton beam therapy system dedicated to spot-scanning increases accuracy with moving tumors by real-time imaging and gating and reduces equipment size.

PURPOSE: A proton beam therapy (PBT) system has been designed which dedicates to spot-scanning and has a gating function employing the fluoroscopy-based real-time-imaging of internal fiducial markers near tumors. The dose distribution and treatment time of the newly designed real-time-image gated, spot-scanning proton beam therapy (RGPT) were compared with free-breathing spot-scanning proton beam therapy (FBPT) in a simulation. MATERIALS AND METHODS: In-house simulation tools and treatment planning system VQA (Hitachi, Ltd., Japan) were used for estimating the dose distribution and treatment time. Simulations were performed for 48 motion parameters (including 8 respiratory patterns and 6 initial breathing timings) on CT data from two patients, A and B, with hepatocellular carcinoma and with clinical target volumes 14.6 cc and 63.1 cc. The respiratory patterns were derived from the actual trajectory of internal fiducial markers taken in X-ray real-time tumor-tracking radiotherapy (RTRT). RESULTS: With FBPT, 9/48 motion parameters achieved the criteria of successful delivery for patient A and 0/48 for B. With RGPT 48/48 and 42/48 achieved the criteria. Compared with FBPT, the mean liver dose was smaller with RGPT with statistical significance (p<0.001); it decreased from 27% to 13% and 28% to 23% of the prescribed doses for patients A and B, respectively. The relative lengthening of treatment time to administer 3 Gy (RBE) was estimated to be 1.22 (RGPT/FBPT: 138 s/113 s) and 1.72 (207 s/120 s) for patients A and B, respectively. CONCLUSIONS: This simulation study demonstrated that the RGPT was able to improve the dose distribution markedly for moving tumors without very large treatment time extension. The proton beam therapy system dedicated to spot-scanning with a gating function for real-time imaging increases accuracy with moving tumors and reduces the physical size, and subsequently the cost of the equipment as well as of the building housing the equipment.

Integration of a real-time tumor monitoring system into gated proton spot-scanning beam therapy: an initial phantom study using patient tumor trajectory data.

PURPOSE: In spot-scanning proton therapy, the interplay effect between tumor motion and beam delivery leads to deterioration of the dose distribution. To mitigate the impact of tumor motion, gating in combination with repainting is one of the most promising methods that have been proposed. This study focused on a synchrotron-based spot-scanning proton therapy system integrated with real-time tumor monitoring. The authors investigated the effectiveness of gating in terms of both the delivered dose distribution and irradiation time by conducting simulations with patients' motion data. The clinically acceptable range of adjustable irradiation control parameters was explored. Also, the relation between the dose error and the characteristics of tumor motion was investigated. METHODS: A simulation study was performed using a water phantom. A gated proton beam was irradiated to a clinical target volume (CTV) of 5 × 5 × 5 cm(3), in synchronization with lung cancer patients' tumor trajectory data. With varying parameters of gate width, spot spacing, and delivered dose per spot at one time, both dose uniformity and irradiation time were calculated for 397 tumor trajectory data from 78 patients. In addition, the authors placed an energy absorber upstream of the phantom and varied the thickness to examine the effect of changing the size of the Bragg peak and the number of required energy layers. The parameters with which 95% of the tumor trajectory data fulfill our defined criteria were accepted. Next, correlation coefficients were calculated between the maximum dose error and the tumor motion characteristics that were extracted from the tumor trajectory data. RESULTS: With the assumed CTV, the largest percentage of the data fulfilled the criteria when the gate width was ± 2 mm. Larger spot spacing was preferred because it increased the number of paintings. With a prescribed dose of 2 Gy, it was difficult to fulfill the criteria for the target with a very small effective depth (the sum of an assumed energy absorber's thickness and the target depth in the phantom) because of the sharpness of the Bragg peak. However, even shallow targets could be successfully irradiated by employing an adequate number of paintings and by placing an energy absorber of sufficient thickness to make the effective target depth more than 12 cm. The authors also observed that motion in the beam direction was the main cause of dose distortion, followed by motion in the lateral plane perpendicular to the scan direction. CONCLUSIONS: The results suggested that by properly adjusting irradiation control parameters, gated proton spot-scanning beam therapy can be robust to target motion. This is an important first step toward establishing treatment plans in real patient geometry.

The precision of respiratory-gated delivery of synchrotron-based pulsed beam proton therapy.

A synchrotron-based proton therapy system operates in a low repetition rate pulsed beam delivery mode. Unlike cyclotron-based beam delivery, there is no guarantee that a synchrotron beam can be delivered effectively or precisely under the respiratory-gated mode. To evaluate the performance of gated synchrotron treatment, we simulated proton beam delivery in the synchrotron-based respiratory-gated mode using realistic patient breathing signals. Parameters used in the simulation were respiratory motion traces (70 traces from 24 patients), respiratory gate levels (10%, 20% and 30% duty cycles at the exhalation phase) and synchrotron magnet excitation cycles (T(cyc)) (fixed T(cyc) mode: 2.7, 3.0-6.0 s and each patient breathing cycle, and variable T(cyc) mode). The simulations were computed according to the breathing trace in which the proton beams were delivered. In the shorter fixed T(cyc) (<4 s), most of the proton beams were delivered uniformly to the target during the entire expiration phase of the respiratory cycle. In the longer fixed T(cyc) (>4 s) and the variable T(cyc) mode, the proton beams were not consistently delivered during the end-expiration phase of the respiratory cycle. However we found that the longer and variable T(cyc) operation modes delivered proton beams more precisely during irregular breathing.

Efficiency of respiratory-gated delivery of synchrotron-based pulsed proton irradiation.

Significant differences exist in respiratory-gated proton beam delivery with a synchrotron-based accelerator system when compared to photon therapy with a conventional linear accelerator. Delivery of protons with a synchrotron accelerator is governed by a magnet excitation cycle pattern. Optimal synchronization of the magnet excitation cycle pattern with the respiratory motion pattern is critical to the efficiency of respiratory-gated proton delivery. There has been little systematic analysis to optimize the accelerator's operational parameters to improve gated treatment efficiency. The goal of this study was to estimate the overall efficiency of respiratory-gated synchrotron-based proton irradiation through realistic simulation. Using 62 respiratory motion traces from 38 patients, we simulated respiratory gating for duty cycles of 30%, 20% and 10% around peak exhalation for various fixed and variable magnet excitation patterns. In each case, the time required to deliver 100 monitor units in both non-gated and gated irradiation scenarios was determined. Based on results from this study, the minimum time required to deliver 100 MU was 1.1 min for non-gated irradiation. For respiratory-gated delivery at a 30% duty cycle around peak exhalation, corresponding average delivery times were typically three times longer with a fixed magnet excitation cycle pattern. However, when a variable excitation cycle was allowed in synchrony with the patient's respiratory cycle, the treatment time only doubled. Thus, respiratory-gated delivery of synchrotron-based pulsed proton irradiation is feasible and more efficient when a variable magnet excitation cycle pattern is used.