Kyung-Wook Jee

Impacts of gantry angle dependent scanning beam properties on proton PBS treatment

While proton beam models in treatment planning systems are generally assumed invariant with respect to the beam deliveries at different gantry angles. Physical properties of scanning pencil beams can change. The gantry angle dependent properties include the delivered charge to the monitor unit chamber, the spot position and the spot shape. The aim of this study is to investigate the extent of the changes and their dosimetric impacts using historical pencil beam scanning (PBS) treatment data. Online beam delivery records at the time of the patient-specific qualify assurance were retrospectively collected for a total of 34 PBS fields from 28 patients treated at our institution. For each field, proton beam properties at two different gantry angles (the planned and zero gantry angles) were extracted by a newly-developed machine log analysis method and used to reconstruct the delivered dose distributions in the cubic water phantom geometry. The reconstructed doses at the two different angles and a planar dose measurement by a 2D ion-chamber array were compared and the dosimetric impacts of the gantry angle dependency were accessed by a 3D γ-index analysis. In addition, the pencil beam spot size was independently characterized as a function of the gantry angle and the beam energy. The dosimetric effects of the perturbed beam shape were also investigated. Comparisons of spot-by-spot beam positions between both gantry angles show a mean deviation of 0.4 and 0.7 mm and a standard deviation of 0.3 and 0.6 mm for x and y directions, respectively. The delivered giga-protons per spot show a percent mean difference and a standard deviation of 0.01% and 0.3%, respectively, from each planned spot weight. These small deviations lead to an excellent agreement in dose comparisons with an average γ passing rate of 99.1%. When each calculation for both planned and zero gantry angles was compared to the measurement, a high correlation in γ values was also observed, also indicating the dosimetric differences are small when a field is delivered at different gantry angles. Utilizing the online beam delivery records, the gantry angle dependencies of the PBS beam delivery were assessed and quantified. The study confirms the variations of the physical properties to be sufficiently small within the clinical tolerances without taking into account the gantry angle variation.

Investigation of real tissue water equivalent path lengths using an efficient dose extinction method

For proton therapy, an accurate conversion of CT HU to relative stopping power (RSP) is essential. Validation of the conversion based on real tissue samples is more direct than the current practice solely based on tissue substitutes and can potentially address variations over the population. Based on a novel dose extinction method, we measured water equivalent path lengths (WEPL) on animal tissue samples to evaluate the accuracy of CT HU to RSP conversion and potential variations over a population. A broad proton beam delivered a spread out Bragg peak to the samples sandwiched between a water tank and a 2D ion-chamber detector. WEPLs of the samples were determined from the transmission dose profiles measured as a function of the water level in the tank. Tissue substitute inserts and Lucite blocks with known WEPLs were used to validate the accuracy. A large number of real tissue samples were measured. Variations of WEPL over different batches of tissue samples were also investigated. The measured WEPLs were compared with those computed from CT scans with the Stoichiometric calibration method. WEPLs were determined within  ±0.5% percentage deviation (% std/mean) and  ±0.5% error for most of the tissue surrogate inserts and the calibration blocks. For biological tissue samples, percentage deviations were within  ±0.3%. No considerable difference (<1%) in WEPL was observed for the same type of tissue from different sources. The differences between measured WEPLs and those calculated from CT were within 1%, except for some bony tissues. Depending on the sample size, each dose extinction measurement took around 5 min to produce ~1000 WEPL values to be compared with calculations. This dose extinction system measures WEPL efficiently and accurately, which allows the validation of CT HU to RSP conversions based on the WEPL measured for a large number of samples and real tissues.

A minimum-phase Shinnar-Le Roux spectral-spatial excitation RF pulse for simultaneous water and lipid suppression in 1H-MRSI of body extremities

PURPOSE To develop a spectral-spatial (SPSP) excitation RF pulse for simultaneous water and lipid suppression in proton (1H) magnetic resonance spectroscopic imaging (MRSI) of body extremities. METHODS An SPSP excitation pulse is designed to excite Creatine (Cr) and Choline (Cho) metabolite signals while suppressing the overwhelming water and lipid signals. The SPSP pulse is designed using a recently proposed multidimensional Shinnar-Le Roux (SLR) RF pulse design method. A minimum-phase spectral selectivity profile is used to minimize signal loss from T2⁎ decay. RESULTS The performance of the SPSP pulse is evaluated via Bloch equation simulations and phantom experiments. The feasibility of the proposed method is demonstrated using three-dimensional, short repetition-time, free induction decay-based 1H-MRSI in the thigh muscle at 3T. CONCLUSION The proposed SPSP excitation pulse is useful for simultaneous water and lipid suppression. The proposed method enables new applications of high-resolution 1H-MRSI in body extremities.

What Is the Effect of High-dose Radiation on Bone in Patients With Sacral Chordoma? A CT Study

Background Effects of high-dose radiation using protons and photons on bone are relatively unexplored, but high rates of insufficiency fractures are reported, and the causes of this are incompletely understood. Imaging studies with pre- and postradiation scans can help one understand the effect of radiation on bone. Questions/purposes The purpose of this study was to assess the effects of high-dose radiation on the trabecular density of bone in the sacrum using CT-derived Hounsfield units (HU). Methods Between 2009 and 2015, we treated 57 patients (older then 18 years) with sacral chordoma. Fourteen (25%) of them were treated with radiation only. The general indication for this approach is inoperability resulting from tumor size. Forty-two (74%) patients were treated with transverse sacral resections and high-dose radiotherapy (using either protons or photons or a combination) before surgery and after surgery. During this time period, our indication for this approach generally was symptomatic sacral chordoma in which resection would prevent further growth and reasonable sacrifice of nerve roots was possible. Of those patients, 21 (50%) had CT scans both before and after radiation treatment. We used HU as a surrogate for bone density. CT uses HU to derive information on tissue and bone quantity. A recent study presented reference HU values for normal (mean 133 ± 38 HU), osteoporotic (101 ± 25 HU), and osteopenic bone (79 ± 32 HU). To adjust for scanning protocol-induced changes in HU, we calculated the ratio between bone inside and outside the radiation field rather than using absolute values. To assess the effect of radiation, we tested whether there was a difference in ratio (sacrum/L1) before and after radiation. A control measurement was performed (L2/L1) and also tested for a difference before and after radiation. Statistical analyses were performed using the paired t-test. Results The effects of radiation appeared confined to the intended field, because the bone density outside the treated field was not observed to decrease. The ratio of HU (a surrogate for bone density) in L2 relative to L1 did not change after radiotherapy (preradiation mean: 0.979 ± 0.009, postradiation mean: 0.980 ± 0.009, mean difference outside the radiation field: -0.001, 95% confidence interval [CI], -0.009 to 0.007, p = 0.799). The ratio of HU within the radiation field relative to L1 decreased after radiotherapy (preradiation mean: 0.895 ± 0.050, postradiation mean: 0.658 ± 0.050, mean difference inside the radiation field: 0.237, 95% CI, 0.187-0.287, p < 0.001), suggesting the bone density stayed the same outside the radiation field but decreased inside the radiation field. Conclusions Trabecular bone density decreased after high-dose radiation therapy in a small group of patients with sacral chordoma. High-dose radiation is increasingly gaining acceptance for treating sacral malignancies; further long-term prospective studies using calibrated CT scanners and preferably bone biopsies are needed. Level of Evidence Level IV, therapeutic study.

TU-FG-BRB-02: The Impact of Using Dual-Energy CT for Determining Proton Stopping Powers: Comparison Between Theory and Experiments

PURPOSE To evaluate the clinical performance of dual-energy CT (DECT) in determining proton stopping power ratios (SPR) and demonstrate advantages over conventional single-energy CT (SECT). METHODS SECT and DECT scans of tissue-equivalent plastics as well as animal meat samples are performed with a Siemens SOMATOM Definition Flash. The methods of Schneider et al. (1996) and Bourque et al. (2014) are used to determine proton SPR on SECT and DECT images, respectively. Waterequivalent path length (WEPL) measurements of plastics and tissue samples are performed with a 195 MeV proton beam. WEPL values are determined experimentally using the depth-dose shift and dose extinction methods. RESULTS Comparison between CT-based and experimental WEPL is performed for 12 tissue-equivalent plastic as well as 6 meat boxes containing animal liver, kidney, heart, stomach, muscle and bones. For plastic materials, results show a systematic improvement in determining SPR with DECT, with a mean absolute error of 0.4% compared to 1.7% for SECT. For the meat samples, preliminary results show the ability for DECT to determine WEPL with a mean absolute value of 1.1% over all meat boxes. CONCLUSION This work demonstrates the potential in using DECT for determining proton SPR with plastic materials in a clinical context. Further work is required to show the benefits of DECT for tissue samples. While experimental uncertainties could be a limiting factor to show the benefits of DECT over SECT for the meat samples, further work is required to adapt the DECT formalism in the context of clinical use, where noise and artifacts play an important role.

SU-C-207A-05: Feature Based Water Equivalent Path Length (WEPL) Determination for Proton Radiography by the Technique of Time Resolved Dose Measurement

PURPOSE Studies show that WEPL can be determined from modulated dose rate functions (DRF). However, the previous calibration method based on statistics of the DRF is sensitive to energy mixing of protons due to scattering through different materials (termed as range mixing here), causing inaccuracies in the determination of WEPL. This study intends to explore time-domain features of the DRF to reduce the effect of range mixing in proton radiography (pRG) by this technique. METHODS An amorphous silicon flat panel (PaxScan™ 4030CB, Varian Medical Systems, Inc., Palo Alto, CA) was placed behind phantoms to measure DRFs from a proton beam modulated by a specially designed modulator wheel. The performance of two methods, the previously used method based on the root mean square (RMS) and the new approach based on time-domain features of the DRF, are compared for retrieving WEPL and RSP from pRG of a Gammex phantom. RESULTS Calibration by T80 (the time point for 80% of the major peak) was more robust to range mixing and produced WEPL with improved accuracy. The error of RSP was reduced from 8.2% to 1.7% for lung equivalent material, with the mean error for all other materials reduced from 1.2% to 0.7%. The mean error of the full width at half maximum (FWHM) of retrieved inserts was decreased from 25.85% to 5.89% for the RMS and T80 method respectively. Monte Carlo simulations in simplified cases also demonstrated that the T80 method is less sensitive to range mixing than the RMS method. CONCLUSION WEPL images have been retrieved based on single flat panel measured DRFs, with inaccuracies reduced by exploiting time-domain features as the calibration parameter. The T80 method is validated to be less sensitive to range mixing and can thus retrieve the WEPL values in proximity of interfaces with improved numerical and spatial accuracy for proton radiography.

TU‐FG‐BRB‐10: A New Approach to Proton Radiography Using the Beamline X‐Ray Flat Panel

PURPOSE Proton radiography, which images the patients with the same type of particles that they are to be treated with, is a promising approach for image guidance and range uncertainties reduction. This study aimed to realize quality proton radiography by measuring dose rate functions (DRF) in time domain using a single flat panel and retrieve water equivalent path length (WEPL) from them. METHODS An amorphous silicon flat panel (PaxScan™ 4030CB, Varian Medical Systems, Inc., Palo Alto, CA) was placed behind phantoms to measure DRFs from a proton beam modulated by the modulator wheel. To retrieve WEPL and RSP, calibration models based on the intensity of DRFs only, root mean square (RMS) of DRFs only and the intensity weighted RMS were tested. The quality of obtained WEPL images (in terms of spatial resolution and level of details) and the accuracy of WEPL were compared. RESULTS RSPs for most of the Gammex phantom inserts were retrieved within ± 1% errors by calibration models based on the RMS and intensity weighted RMS. The mean percentage error for all inserts was reduced from 1.08% to 0.75% by matching intensity in the calibration model. In specific cases such as the insert with a titanium rod, the calibration model based on RMS only fails while the that based on intensity weighted RMS is still valid. The quality of retrieved WEPL images were significantly improved for calibration models including intensity matching. CONCLUSION For the first time, a flat panel, which is readily available in the beamline for image guidance, was tested to acquire quality proton radiography with WEPL accurately retrieved from it. This technique is promising to be applied for image-guided proton therapy as well as patient specific RSP determination to reduce uncertainties of beam ranges.

SU-G-TeP4-04: An Automated Monte Carlo Based QA Framework for Pencil Beam Scanning Treatments

PURPOSE Prior to treating new PBS field, multiple (three) patient-field-specific QA measurements are performed: two 2D dose distributions at shallow depth (M1) and at the tumor depth (M2) with treatment hardware at zero gantry angle; one 2D dose distribution at iso-center (M3) without patient specific devices at the planned gantry angle. This patient-specific QA could be simplified by the use of MC model. The results of MC model commissioning for a spot-scanning system and the fully automated TOPAS/MC-based QA framework will be presented. METHODS We have developed in-house MC interface to access a TPS (Astroid) database from a computer cluster remotely. Once a plan is identified, the interface downloads information for the MC simulations, such as patient images, apertures points, and fluence maps and initiates calculations in both the patient and QA geometries. The resulting calculations are further analyzed to evaluate the TPS dose accuracy and the PBS delivery. RESULTS The Monte Carlo model of our system was validated within 2.0 % accuracy over the whole range of the dose distribution (proximal/shallow part, as well as target dose part) due to the location of the measurements. The averaged range difference after commissioning was 0.25 mm over entire treatment ranges, e.g., 6.5 cm to 31.6 cm. CONCLUSION As M1 depths range typically from 1 cm to 4 cm from the phantom surface, The Monte Carlo model of our system was validated within +- 2.0 % in absolute dose level over a whole treatment range. The averaged range difference after commissioning was 0.25 mm over entire treatment ranges, e.g., 6.5 cm to 31.6 cm. This work was supported by NIH/NCI under CA U19 21239.

SU-F-T-177: Impacts of Gantry Angle Dependent Scanning Beam Properties for Proton Treatment

PURPOSE In pencil beam scanning (PBS), the delivered spot MU, position and size are slightly different at different gantry angles. We investigated the level of delivery uncertainty at different gantry angles through a log file analysis. METHODS 34 PBS fields covering full 360 degrees gantry angle spread were collected retrospectively from 28 patients treated at our institution. All fields were delivered at zero gantry angle and the prescribed gantry angle, and measured at isocenter with the MatriXX 2D array detector at the prescribed gantry angle. The machine log files were analyzed to extract the delivered MU per spot and the beam position from the strip ionization chambers in the treatment nozzle. The beam size was separately measured as a function of gantry angle and beam energy. Using this information, the dose was calculated in a water phantom at both gantry angles and compared to the measurement using the 3D γ-index at 2mm/2%. RESULTS The spot-by-spot difference between the beam position in the log files from the delivery at the two gantry angles has a mean of 0.3 and 0.4 mm and a standard deviation of 0.6 and 0.7 mm for × and y directions, respectively. Similarly, the spot-by-spot difference between the MU in the log files from the delivery at the two gantry angles has a mean 0.01% and a standard deviation of 0.7%. These small deviations lead to an excellent agreement in dose calculations with an average γ pass rate for all fields being approximately 99.7%. When each calculation is compared to the measurement, a high correlation in γ was also found. CONCLUSION Using machine logs files, we verified that PBS beam delivery at different gantry angles are sufficiently small and the planned spot position and MU. This study brings us one step closer to simplifying our patient-specific QA.

SU-F-J-193: Efficient Dose Extinction Method for Water Equivalent Path Length (WEPL) of Real Tissue Samples for Validation of CT HU to Stopping Power Conversion

PURPOSE For proton therapy, an accurate model of CT HU to relative stopping power (RSP) conversion is essential. In current practice, validation of these models relies solely on measurements of tissue substitutes with standard compositions. Validation based on real tissue samples would be much more direct and can address variations between patients. This study intends to develop an efficient and accurate system based on the concept of dose extinction to measure WEPL and retrieve RSP in biological tissue in large number of types. METHODS A broad AP proton beam delivering a spread out Bragg peak (SOBP) is used to irradiate the samples with a Matrixx detector positioned immediately below. A water tank was placed on top of the samples, with the water level controllable in sub-millimeter by a remotely controlled dosing pump. While gradually lowering the water level with beam on, the transmission dose was recorded at 1 frame/sec. The WEPL were determined as the difference between the known beam range of the delivered SOBP (80%) and the water level corresponding to 80% of measured dose profiles in time. A Gammex 467 phantom was used to test the system and various types of biological tissue was measured. RESULTS RSP for all Gammex inserts, expect the one made with lung-450 material (<2% error), were determined within ±0.5% error. Depends on the WEPL of investigated phantom, a measurement takes around 10 min, which can be accelerated by a faster pump. CONCLUSION Based on the concept of dose extinction, a system was explored to measure WEPL efficiently and accurately for a large number of samples. This allows the validation of CT HU to stopping power conversions based on large number of samples and real tissues. It also allows the assessment of beam uncertainties due to variations over patients, which issue has never been sufficiently studied before.