Quantitative, real-time scintillation imaging for experimental comparison of different dose and dose rate estimations in UHDR proton pencil beams.

BACKGROUND: Ultra-high dose rate radiotherapy (UHDR-RT) has demonstrated normal tissue sparing capabilities, termed the FLASH effect; however, available dosimetry tools make it challenging to characterize the UHDR beams with sufficiently high concurrent spatial and temporal resolution. Novel dosimeters are needed for safe clinical implementation and improved understanding of the effect of UHDR-RT. PURPOSE: Ultra-fast scintillation imaging has been shown to provide a unique tool for spatio-temporal dosimetry of conventional cyclotron pencil beam scanning (PBS) deliveries, indicating the potential use for characterization of UHDR PBS proton beams. The goal of this work is to introduce this novel concept and demonstrate its capabilities in recording high-resolution dose rate maps at FLASH-capable proton beam currents, as compared to log-based dose rate calculation, internally developed UHDR beam simulation, and a fast point detector (EDGE diode). METHODS: The light response of a scintillator sheet located at isocenter and irradiated by PBS proton fields (40-210 nA, 250 MeV) was imaged by an ultra-fast iCMOS camera at 4.5-12 kHz sampling frequency. Camera sensor and image intensifier gain were optimized to maximize the dynamic range; the camera acquisition rate was also varied to evaluate the optimal sampling frequency. Large field delivery enabled flat field acquisition for evaluation of system response homogeneity. Image intensity was calibrated to dose with film and the recorded spatio-temporal data was compared to a PPC05 ion chamber, log-based reconstruction, and EDGE diode. Dose and dose rate linearity studies were performed to evaluate agreement under various beam conditions. Calculation of full-field mean and PBS dose rate maps were calculated to highlight the importance of high resolution, full-field information in UHDR studies. RESULTS: Camera response was linear with dose (R2 = 0.997) and current (R22 = 0.98) in the range from 2-22 Gy and 40-210 nA, respectively, when compared to ion chamber readings. The deviation of total irradiation time calculated with the imaging system from the log file recordings decreased from 0.07% to 0.03% when imaging at 12 kfps versus 4.5 kfps. Planned and delivered spot positions agreed within 0.2 ± $\pm$ 0.1 mm and total irradiation time agreed within 0.2 ± $\pm$ 0.2 ms when compared with the log files, indicating the high concurrent spatial and temporal resolution. For all deliveries, the PBS dose rate measured at the diode location agreed between the imaging and the diode within 3% ± $\pm$ 2% and with the simulation within 5% ± $\pm$ 3% CONCLUSIONS: Full-field mapping of dose and dose rate is imperative for complete understanding of UHDR PBS proton dose delivery. The high linearity and various spatiotemporal metric reporting capabilities confirm the continued use of this camera system for UHDR beam characterization, especially for spatially resolved dose rate information.

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.

Characterization of a diode dosimeter for UHDR FLASH radiotherapy.

BACKGROUND: Ultra-high dose rate (UHDR) FLASH beams typically deliver dose at rates of  >40 Gy/sec. Characterization of these beams with respect to dose, mean dose rate, and dose per pulse requires dosimeters which exhibit high temporal resolution and fast readout capabilities. PURPOSE: A diode EDGE Detector with a newly designed electrometer has been characterized for use in an UHDR electron beam and demonstrated appropriateness for UHDR FLASH radiotherapy dosimetry. METHODS: Dose linearity, mean dose rate, and dose per pulse dependencies of the EDGE Detector were quantified and compared with dosimeters including a W1 scintillator detector, radiochromic film, and ionization chamber that were irradiated with a 10 MeV UHDR beam. The dose, dose rate, and dose per pulse were controlled via an in-house developed scintillation-based feedback mechanism, repetition rate of the linear accelerator, and source-to-surface distance, respectively. Depth-dose profiles and temporal profiles at individual pulse resolution were compared to the film and scintillation measurements, respectively. The radiation-induced change in response sensitivity was quantified via irradiation of ∼5kGy. RESULTS: The EDGE Detector agreed with film measurements in the measured range with varying dose (up to 70 Gy), dose rate (nearly 200 Gy/s), and dose per pulse (up to 0.63 Gy/pulse) on average to within 2%, 5%, and 1%, respectively. The detector also agreed with W1 scintillation detector on average to within 2% for dose per pulse (up to 0.78 Gy/pulse). The EDGE Detector signal was proportional to ion chamber (IC) measured dose, and mean dose rate in the bremsstrahlung tail to within 0.4% and 0.2% respectively. The EDGE Detector measured percent depth dose (PDD) agreed with film to within 3% and per pulse output agreed with W1 scintillator to within -6% to +5%. The radiation-induced response decrease was 0.4% per kGy. CONCLUSIONS: The EDGE Detector demonstrated dose linearity, mean dose rate independence, and dose per pulse independence for UHDR electron beams. It can quantify the beam spatially, and temporally at sub millisecond resolution. It's robustness and individual pulse detectability of treatment deliveries can potentially lead to its implementation for in vivo FLASH dosimetry, and dose monitoring.

The probe-format graphite calorimeter, Aerrow, for absolute dosimetry in ultrahigh pulse dose rate electron beams.

PURPOSE: The purpose of this investigation is to evaluate the use of a probe-format graphite calorimeter, Aerrow, as an absolute and relative dosimeter of high-energy pulse dose rate (UHPDR) electron beams for in-water reference and depth-dose-type measurements, respectively. METHODS: In this paper, the calorimeter system is used to investigate the potential influence of dose per pulses delivered up to 5.6 Gy, the number of pulses delivered per measurement, and its potential for relative measurement (depth-dose curve measurement). The calorimeter system is directly compared against an Advanced Markus ion chamber. The finite element method was used to calculate heat transfer corrections along the percentage depth dose of a 20-MeV electron beam. Monte Carlo-calculated dose conversion factors necessary to calculate absorbed dose-to-water at a point from the measured dose-to-graphite are also presented. RESULTS: The comparison of Aerrow against a fully calibrated Advanced Markus chamber, corrected for the saturation effect, has shown consistent results in terms of dose-to-water determination. The measured reference depth is within 0.5 mm from the expected value from Monte Carlo simulation. The relative standard uncertainty estimated for Aerrow was 1.06%, which is larger compared to alanine dosimetry (McEwen et al. https://doi.org/10.1088/0026-1394/52/2/272) but has the advantage of being a real-time detector. CONCLUSION: In this investigation, it was demonstrated that the Aerrow probe-type graphite calorimeter can be used for relative and absolute dosimetries in water in an UHPDR electron beam. To the author's knowledge, this is the first reported use of an absorbed dose calorimeter for an in-water percentage depth-dose curve measurement. The use of the Aerrow in quasi-adiabatic mode has greatly simplified the signal readout, compared to isothermal mode, as the resistance was directly measured with a high-stability digital multimeter.

Validation of an improved helical diode array and dose reconstruction software using TG-244 datasets and stringent dose comparison criteria.

The original helical ArcCHECK (AC) diode array and associated software for 3D measurement-guided dose reconstruction were characterized and validated; however, recent design changes to the AC required that the subject be revisited. The most important AC change starting in 2014 was a significant reduction in the overresponse of diodes to scattered radiation outside of the direct beam, accom-plished by reducing the amount of high-Z materials adjacent to the diodes. This change improved the diode measurement accuracy, but in the process invalidated the dose reconstruction models that were assembled based on measured data acquired with the older version of the AC. A correction mechanism was intro-duced in the reconstruction software (3DVH) to accommodate this and potential future design changes without requiring updating model parameters. For each permutation of AC serial number and beam model, the user can define in 3DVH a single correction factor which will be used to compensate for the difference in the out-of-field response between the new and original AC designs. The exact value can be determined by minimizing the dose-difference with an ionization chamber or another independent dosimeter. A single value of 1.17, corresponding to the maximum measured out-of-field response difference between the new and old AC, provided satisfactory results for all studied energies (6X, 15X, and flatten-ing filter-free 10XFFF). A library of standard cases recommended by the AAPM TG-244 Report was used for reconstructed dose verification. The overall difference between reconstructed dose and an ion chamber in a water-equivalent phantom in the targets was 0.0% ± 1.4% (1 SD). The reconstructed dose on a homogeneous phantom was also compared to a biplanar diode dosimeter (Delta4) using gamma analysis with 2% (local dose-error normalization) / 2 mm / 10% cutoff criteria. The mean agreement rate was 96.7% ± 3.7%. For the plans common with the previous comparison, the mean agreement rate was 98.3% ± 0.8%, essentially unchanged. We conclude that the proposed software modification adequately addresses the change in the dosimeter response.

Evaluating dosimetric accuracy of flattening filter free compensator-based IMRT: measurements with diode arrays.

PURPOSE: Compensator-based IMRT coupled with the high dose rate flattening filter free (FFF) beams offers an intriguing possibility of delivering an intensity modulated radiation field in just a few seconds. As a first step, the authors evaluate the dosimetric accuracy of the treatment planning system (TPS) FFF beam model with compensators. METHODS: A 6 MV FFF beam from a TrueBeam accelerator (Varian Medical Systems, Palo Alto CA) was modeled in PINNACLE TPS (v. 9.0, Philips Radiation Oncology, Fitchburg WI). Flat brass slabs from 0.3 to 7 cm thick and an 18° brass wedge were used to adjust the beam model. A 2D (MAPCHECK) and 3D (ARCCHECK) diode arrays (Sun Nuclear Corp, Melbourne FL), were investigated for use with the compensator FFF beams. Corrections for diode sensitivity caused by the spectral changes in the beam were introduced. Four compensator plans based on the AAPM TG-119 report were developed. A composite ion chamber measurement, beam by beam MAPCHECK measurements, and a composite ARCCHECK measurement were performed. The array results were analyzed with the same thresholds as in TG-119 report-3%/3 mm with global dose normalization-as well as with the more stringent combinations of the gamma analysis criteria. RESULTS: The FFF beam shows a greater variation of the effective attenuation coefficient with brass thickness due to the prevalence of the low energy photons compared to the conventional 6X beam. As a result, a compromise had to be made while trying to achieve dose agreement for a combination of field sizes, brass thicknesses, and measurement depths (≥5 cm in water). An agreement of measured and calculated dose to within 1% was observed for brass thicknesses up to 2 cm. For the 3 cm slab, an error of up to 2.8% was noted for the field sizes above 10 × 10 cm(2), and up to 3.8% for the 5 × 5 cm(2) field. Both diode arrays exhibit a substantial sensitivity drop as the compensator thickness increases, reaching 10% for a 7 cm brass slab. A simple correction based on the brass thickness along the ray was introduced to counteract this effect. Pooled for five profiles, the average ratio of uncorrected and corrected MAPCHECK to ion chamber readings are 0.966 and 1.008, respectively. With the proper correction, all MAPCHECK measurement to calculation comparisons exhibit 100% γ(3%/3 mm) passing rates with global dose-error normalization. For the TG-119-type plans, the average γ(2%/2 mm) passing rate with local normalization is 94% (range 87.8%-98.3%). The lower ARCCHECK γ-analysis passing rates (corrected for diode sensitivity) are predictable based on the observed PDD discrepancies. However, with the 3%/3 mm thresholds and global normalization, the average γ-analysis passing rate is 96.4% (range 89.9%-100%). CONCLUSIONS: MAPCHECK analysis demonstrates high passing rates with the stringent γ(2%/2 mm) and local normalization criteria combination. The geometry of the ARCCHECK array creates a stress test for the FFF TPS model because of the shallow depth of the entrance diodes and large air cavity. Hence, the ARCCHECK γ-analysis passing rates are lower than with the MAPCHECK, while still on par with TG-119.

Optimizing the accuracy of a helical diode array dosimeter: a comprehensive calibration methodology coupled with a novel virtual inclinometer.

PURPOSE: The goal of any dosimeter is to be as accurate as possible when measuring absolute dose to compare with calculated dose. This limits the uncertainties associated with the dosimeter itself and allows the task of dose QA to focus on detecting errors in the treatment planning (TPS) and/or delivery systems. This work introduces enhancements to the measurement accuracy of a 3D dosimeter comprised of a helical plane of diodes in a volumetric phantom. METHODS: We describe the methods and derivations of new corrections that account for repetition rate dependence, intrinsic relative sensitivity per diode, field size dependence based on the dynamic field size determination, and positional correction. Required and described is an accurate "virtual inclinometer" algorithm. The system allows for calibrating the array directly against an ion chamber signal collected with high angular resolution. These enhancements are quantitatively validated using several strategies including ion chamber measurements taken using a "blank" plastic shell mimicking the actual phantom, and comparison to high resolution dose calculations for a variety of fields: static, simple arcs, and VMAT. A number of sophisticated treatment planning algorithms were benchmarked against ion chamber measurements for their ability to handle a large air cavity in the phantom. RESULTS: Each calibration correction is quantified and presented vs its independent variable(s). The virtual inclinometer is validated by direct comparison to the gantry angle vs time data from machine log files. The effects of the calibration are quantified and improvements are seen in the dose agreement with the ion chamber reference measurements and with the TPS calculations. These improved agreements are a result of removing prior limitations and assumptions in the calibration methodology. Average gamma analysis passing rates for VMAT plans based on the AAPM TG-119 report are 98.4 and 93.3% for the 3%/3 mm and 2%/2 mm dose-error/distance to agreement threshold criteria, respectively, with the global dose-error normalization. With the local dose-error normalization, the average passing rates are reduced to 94.6 and 85.7% for the 3%/3 mm and 2%/2 mm criteria, respectively. Some algorithms in the convolution/superposition family are not sufficiently accurate in predicting the exit dose in the presence of a 15 cm diameter air cavity. CONCLUSIONS: Introduction of the improved calibration methodology, enabled by a robust virtual inclinometer algorithm, improves the accuracy of the dosimeter's absolute dose measurements. With our treatment planning and delivery chain, gamma analysis passing rates for the VMAT plans based on the AAPM TG-119 report are expected to be above 91% and average at about 95% level for γ(3%/3 mm) with the local dose-error normalization. This stringent comparison methodology is more indicative of the true VMAT system commissioning accuracy compared to the often quoted dose-error normalization to a single high value.

Characterization of a multi-axis ion chamber array.

PURPOSE: The aim of this work was to characterize a multi-axis ion chamber array (IC PROFILER; Sun Nuclear Corporation, Melbourne, FL, USA) that has the potential to simplify the acquisition of LINAC beam data. METHODS: The IC PROFILER (or panel) measurement response was characterized with respect to radiation beam properties, including dose, dose per pulse, pulse rate frequency (PRF), and energy. Panel properties were also studied, including detector-calibration stability, power-on time, backscatter dependence, and the panel's agreement with water tank measurements [profiles, fractional depth dose (FDD), and output factors]. RESULTS: The panel's relative deviation was typically within (+/-) 1% of an independent (or nominal) response for all properties that were tested. Notable results were (a) a detectable relative field shape change of approximately 1% with linear accelerator PRF changes; (b) a large range in backscatter thickness had a minimal effect on the measured dose distribution (typically less than 1%); (c) the error spread in profile comparison between the panel and scanning water tank (Blue Phantom, CC13; IBA Schwarzenbruck, DE) was approximately (+/-) 0.75%. CONCLUSIONS: The ability of the panel to accurately reproduce water tank profiles, FDDs, and output factors is an indication of its abilities as a dosimetry system. The benefits of using the panel versus a scanning water tank are less setup time and less error susceptibility. The same measurements (including device setup and breakdown) for both systems took 180 min with the water tank versus 30 min with the panel. The time-savings increase as the measurement load is increased.

Calibration of a novel four-dimensional diode array.

PURPOSE: The aim of this work is to develop effective calibration methods for a novel fourdimensional (4D) diode array for pretreatment verification of intensity-modulated radiation therapy (IMRT) and rotational therapy. METHODS: A novel 4D diode array (ArcCHECK, Sun Nuclear, Melbourne, FL) was developed to meet the needs of appropriate and efficient quality assurance for IMRT and especially rotational radiotherapy. The diode array presents a consistent detector image in beam's eye view at arbitrary gantry angles due to isotropic arrangement of diodes in a three-dimensional (3D) cylindrical phantom. The 50 ms simultaneous update of all diodes on the detector array (fourth dimension) makes it capable of time-resolved beam delivery analysis with any rotational delivery techniques. The calibration procedure consisted of delivering and measuring a series of calibration beams with 5.8 degrees angular spacing surrounding the cylindrical diode array. Correction factors for diode intrinsic sensitivity and directional response dependence were derived from these measurements. A real-time algorithm to derive gantry angles based on the detector signal was developed to interpolate and apply the corresponding angular correction factors. RESULTS: The calibration was validated with ion chamber scanned beam profiles in a 3D water tank. Excellent agreement was observed between diode array measurement and treatment planning system calculation. The accuracy of the gantry angle derivation algorithm was within 1 degree which caused a less than 0.2% dosimetric uncertainty. CONCLUSIONS: With the proposed calibration method and the automatic gantry angle derivation algorithm, the 4D diode array achieved isotropic detector response and is suitable for both IMRT and rotational therapy pretreatment verification.

Novel dosimetric phantom for quality assurance of volumetric modulated arc therapy.

The objective of this work is to assess the suitability and performance of a new dosimeter system with a novel geometry for the quality assurance (QA) of volumetric modulated arc therapy (VMAT). The new dosimeter system consists of a hollow cylinder (15 and 25 cm inner and outer diameters) with 124 diodes embedded in the phantom's cylindrical wall forming four rings of detectors. For coplanar beams, the cylindrical geometry and the ring diode pattern offer the advantage of invariant perpendicular incidence on the beam central axis for any gantry angle and also have the benefit of increasing the detector density as both walls of the cylinder sample the beam. Other advantages include real-time readout and reduced weight with the hollow phantom shape. A calibration method taking into account the variation in radiation sensitivity of the diodes as a function of gantry angle was developed and implemented. In this work, the new dosimeter system was used in integrating mode to perform composite dose measurements along the cylindrical surface supporting the diodes. The reproducibility of the dosimeter response and the angular dependence of the diodes were assessed using simple 6 MV photon static beams. The performance of the new dosimeter system for VMAT QA was then evaluated using VMAT plans designed for a head and neck, an abdominal sarcoma, and a prostate patient. These plans were optimized with 90 control points (CPs) and additional versions of each plan were generated by increasing the number of CPs to 180 and 360 using linear interpolation. The relative dose measured with the dosimeter system for the VMAT plans was compared to the corresponding TPS dose map in terms of relative dose difference (% deltaD) and distance to agreement (DTA). The dosimeter system's sensitivity to gantry rotation offset and scaling errors as well as setup errors was also evaluated. For static beams, the dosimeter system offered good reproducibility and demonstrated small residual diode angular dependence after calibration. For VMAT deliveries, the agreement between measured and calculated doses was good with > or = 86.4% of the diodes satisfying 3% of % deltaD or 2 mm DTA for the 180 CP plans. The phantom offered sufficient sensitivity for the detection of small gantry rotation offset (3 degrees) and scaling errors (1 degree) as well as phantom setup errors of 1 mm, although the results were plan dependent. With its novel geometry, the dosimeter system was also able to experimentally demonstrate the discretization effect of the number of CPs used in the TPS to simulate a continuous arc. These results demonstrate the suitability of the new dosimeter system for the patient-specific QA of VMAT plans and suggest that the dosimeter system can be an effective tool in the routine QA and commissioning of treatment machines capable of VMAT delivery and cone-beam CT image guidance.

A prototype quantitative film scanner for radiochromic film dosimetry.

We have developed a high resolution, quantitative, two-dimensional optical film scanner for use with a commercial high sensitivity radiochromic film (RCF) for measuring single fraction external-beam radiotherapy dose distributions. The film scanner was designed to eliminate artifacts commonly observed in RCF dosimetry. The scanner employed a stationary light source and detector with a moving antireflective glass film platen attached to a high precision computerized X-Y translation stage. An ultrabright red light emitting diode (LED) with a peak output at 633 nm and full width at half maximum (FWHM) of 16 nm was selected as the scanner light source to match the RCF absorption peak. A dual detector system was created using two silicon photodiode detectors to simultaneously measure incident and transmitted light. The LED light output was focused to a submillimeter (FWHM 0.67 mm) spot size, which was determined from a scanning knife-edge technique for measuring Gaussian optical beams. Data acquisition was performed with a 16-bit A/D card in conjunction with commercial software. The linearity of the measured densities on the scanner was tested using a calibrated neutral-density step filter. Sensitometric curves and three IMRT field scans were acquired with a spatial resolution of 1 mm for both radiographic film and RCF. The results were compared with measurements taken with a commercial diode array under identical delivery conditions. The RCF was rotated by 90 deg and rescanned to study orientation effects. Comparison between the RCF and the diode array measurements using percent dose difference and distance-to-agreement criteria produced average passing rates of 99.0% using 3%/3 mm criteria and 96.7% using 2%/2 mm criteria. The same comparison between the radiographic film and diode array measurements resulted in average passing rates 96.6% and 91.6% for the above two criteria, respectively. No measurable light-scatter or interference scanner artifacts were observed. The RCF rotated by 90 deg showed no measurable orientation effect. A scan of a 15 x 15 cm2 area with 1 mm resolution required 22 min to acquire. The LED densitometer provides an accurate film dosimetry system with 1 mm or better resolution. The scanner eliminates the orientation dependence of RCF dosimetry that was previously reported with commercial flatbed scanners.

Important considerations for radiochromic film dosimetry with flatbed CCD scanners and EBT GAFCHROMIC film.

In this study, we present three significant artifacts that have the potential to negatively impact the accuracy and precision of film dosimetry measurements made using GAFCHROMIC EBT radiochromic film when read out with CCD flatbed scanners. Films were scanned using three commonly employed instruments: a Macbeth TD932 spot densitometer, an Epson Expression 1680 CCD array scanner, and a Microtek ScanMaker i900 CCD array scanner. For the two scanners we assessed the variation in optical density (OD) of GAFCHROMIC EBT film with scanning bed position, angular rotation of the film with respect to the scan line direction, and temperature inside the scanner due to repeated scanning. Scanning uniform radiochromic films demonstrated a distinct bowing effect in profiles in the direction of the CCD array with a nonuniformity of up to 17%. Profiles along a direction orthogonal to the CCD array demonstrated a 7% variation. A strong angular dependence was found in measurements made with the flatbed scanners; the effect could not be reproduced with the spot densitometer. An IMRT quality assurance film was scanned twice rotating the film 90' between the scans. For films scanned on the Epson scanner, up to 12% variation was observed in unirradiated EBT films rotated between 0 degrees and 90 degrees, which decreased to approximately 8% for EBT films irradiated to 300 cGy. Variations of up to 80% were observed for films scanned with the Microtek scanner. The scanners were found to significantly increase the film temperature with repeated scanning. Film temperature between 18 and 33 degrees C caused OD changes of approximately 7%. Considering these effects, we recommend adherence to a strict scanning protocol that includes: maintaining the orientation of films scanned on flatbed scanners, limiting scanning to the central portion of the scanner bed, and limiting the number of consecutive scans to minimize changes in OD caused by film heating.