Cheryl Duzenli

Monte Carlo based, patient-specific RapidArc QA using Linac log files

PURPOSE A Monte Carlo (MC) based QA process to validate the dynamic beam delivery accuracy for Varian RapidArc (Varian Medical Systems, Palo Alto, CA) using Linac delivery log files (DynaLog) is presented. Using DynaLog file analysis and MC simulations, the goal of this article is to (a) confirm that adequate sampling is used in the RapidArc optimization algorithm (177 static gantry angles) and (b) to assess the physical machine performance [gantry angle and monitor unit (MU) delivery accuracy]. METHODS Ten clinically acceptable RapidArc treatment plans were generated for various tumor sites and delivered to a water-equivalent cylindrical phantom on the treatment unit. Three Monte Carlo simulations were performed to calculate dose to the CT phantom image set: (a) One using a series of static gantry angles defined by 177 control points with treatment planning system (TPS) MLC control files (planning files), (b) one using continuous gantry rotation with TPS generated MLC control files, and (c) one using continuous gantry rotation with actual Linac delivery log files. Monte Carlo simulated dose distributions are compared to both ionization chamber point measurements and with RapidArc TPS calculated doses. The 3D dose distributions were compared using a 3D gamma-factor analysis, employing a 3%/3 mm distance-to-agreement criterion. RESULTS The dose difference between MC simulations, TPS, and ionization chamber point measurements was less than 2.1%. For all plans, the MC calculated 3D dose distributions agreed well with the TPS calculated doses (gamma-factor values were less than 1 for more than 95% of the points considered). Machine performance QA was supplemented with an extensive DynaLog file analysis. A DynaLog file analysis showed that leaf position errors were less than 1 mm for 94% of the time and there were no leaf errors greater than 2.5 mm. The mean standard deviation in MU and gantry angle were 0.052 MU and 0.355 degrees, respectively, for the ten cases analyzed. CONCLUSIONS The accuracy and flexibility of the Monte Carlo based RapidArc QA system were demonstrated. Good machine performance and accurate dose distribution delivery of RapidArc plans were observed. The sampling used in the TPS optimization algorithm was found to be adequate.

Image filtering for improved dose resolution in CT polymer gel dosimetry

X-ray computed tomography (CT) has been established as a feasible method of performing dosimetry using polyacrylamide gels (PAGs). A small density change occurs in PAG upon irradiation that provides contrast in PAG CT images. However, low dose resolution limits the clinical usefulness of the technique. This work investigates the potential of using image filtering techniques on PAG CT images in order to reduce image noise and improve dose resolution. CT image noise for the scanner and protocol used for the gel images is analyzed and found to be Gaussian distributed and independent of the contrast level in the images. As a result, several filters for reducing spatially invariant noise are investigated: mean, median, midpoint, adaptive mean, alpha-trimmed mean, sigma mean, and a relatively new filter called SUSAN (smallest univalue segment assimilating nucleus). All filters are applied, using 3x3, 5x5, and 7x7 pixel masks, to a CT image of a PAG irradiated with a stereotactic radiosurgery dose distribution. The dose resolution within 95% confidence (D(delta)95%) is calculated and compared for each filtered image, as well the unfiltered image. In addition, the ability of the filters to maintain the spatial integrity of the dose distribution is evaluated and compared. Results clearly indicate that the filters are not equal in their ability to improve D(delta)95% or in their effect on the spatial integrity of the dose distribution. In general, increasing mask size improves D(delta)95% but simultaneously degrades spatial dose information. The mean filter provides the greatest improvement in D(delta)95%, but also the greatest loss of spatial dose information. The SUSAN, mean adaptive, and alpha-trimmed mean filters all provide comparable, but slightly poorer dose resolution. In addition, the SUSAN and adaptive filters both excel at maintaining the spatial distribution of dose and overall are the best performing filters for this application. The midpoint filter, normally useful for Gaussian noise, is poor all-round, dramatically distorting the dose distribution for masks greater than 3x3. The median filter, a common edge preserving noise reduction filter, performs moderately well, but artificially increases high dose gradients. The sigma filter preserves the spatial distribution of dose very well but is least effective at improving dose resolution. In summary, dose resolution can be significantly improved in CT PAG dosimetry through postprocessing of CT images using spatial noise reduction filters. However, such filters are not equal in their ability to improve dose resolution or to maintain the spatial integrity of the dose distribution and an appropriate filter must be chosen depending on clinical demands of the application.

Direct aperture optimization for IMRT using Monte Carlo generated beamlets

This work introduces an EGSnrc-based Monte Carlo (MC) beamlet does distribution matrix into a direct aperture optimization (DAO) algorithm for IMRT inverse planning. The technique is referred to as Monte Carlo-direct aperture optimization (MC-DAO). The goal is to assess if the combination of accurate Monte Carlo tissue inhomogeneity modeling and DAO inverse planning will improve the dose accuracy and treatment efficiency for treatment planning. Several authors have shown that the presence of small fields and/or inhomogeneous materials in IMRT treatment fields can cause dose calculation errors for algorithms that are unable to accurately model electronic disequilibrium. This issue may also affect the IMRT optimization process because the dose calculation algorithm may not properly model difficult geometries such as targets close to low-density regions (lung, air etc.). A clinical linear accelerator head is simulated using BEAMnrc (NRC, Canada). A novel in-house algorithm subdivides the resulting phase space into 2.5 X 5.0 mm2 beamlets. Each beamlet is projected onto a patient-specific phantom. The beamlet dose contribution to each voxel in a structure-of-interest is calculated using DOSXYZnrc. The multileaf collimator (MLC) leaf positions are linked to the location of the beamlet does distributions. The MLC shapes are optimized using direct aperture optimization (DAO). A final Monte Carlo calculation with MLC modeling is used to compute the final dose distribution. Monte Carlo simulation can generate accurate beamlet dose distributions for traditionally difficult-to-calculate geometries, particularly for small fields crossing regions of tissue inhomogeneity. The introduction of DAO results in an additional improvement by increasing the treatment delivery efficiency. For the examples presented in this paper the reduction in the total number of monitor units to deliver is approximately 33% compared to fluence-based optimization methods.

A Monte Carlo approach to validation of FFF VMAT treatment plans for the TrueBeam linac

PURPOSE To commission and benchmark a vendor-supplied (Varian Medical Systems) Monte Carlo phase-space data for the 6 MV flattening filter free (FFF) energy mode on a TrueBeam linear accelerator for the purpose of quality assurance of clinical volumetric modulated arc therapy (VMAT) treatment plans. A method for rendering the phase-space data compatible with BEAMnrc/DOSXYZnrc simulation software package is presented. METHODS Monte Carlo (MC) simulations were performed to benchmark the TrueBeam 6 MV FFF phase space data that have been released by the Varian MC Research team. The simulations to benchmark the phase space data were done in three steps. First, the original phase space which was created on a cylindrical surface was converted into a format that was compatible with BEAMnrc. Second, BEAMnrc was used to create field size specific phase spaces located underneath the jaws. Third, doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Calculated percent depth doses (PDD), transverse profiles, and output factors were compared with measurements for all the fields simulated. After completing the benchmarking study, three stereotactic body radiotherapy (SBRT) VMAT plans created with the Eclipse treatment planning system (TPS) were calculated with Monte Carlo. Ion chamber and film measurements were also performed on these plans. 3D gamma analysis was used to compare Monte Carlo calculation with TPS calculations and with film measurement. RESULTS For the benchmarking study, MC calculated and measured values agreed within 1% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2). Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. With the exception of the 1 × 1 cm(2) field, the agreement between measured and calculated output factors was within 1%. It is of note that excellent agreement in output factors for all field sizes including highly asymmetric fields was achieved without accounting for backscatter into the beam monitor chamber. For the SBRT VMAT plans, the agreement between Monte Carlo and ion chamber point dose measurements was within 1%. Excellent agreement between Monte Carlo, treatment planning system and Gafchromic film dose distribution was observed with over 99% of the points in the high dose volume passing the 3%, 3 mm gamma test. CONCLUSIONS The authors have presented a method for making the Varian IAEA compliant 6 MV FFF phase space file of the TrueBeam linac compatible with BEAMnrc/DOSXYZnrc. After benchmarking the modified phase space against measurement, they have demonstrated its potential for use in MC based quality assurance of complex delivery techniques.

Practice Patterns Analysis of Ocular Proton Therapy Centers: The International OPTIC Survey.

PURPOSE To assess the planning, treatment, and follow-up strategies worldwide in dedicated proton therapy ocular programs. METHODS AND MATERIALS Ten centers from 7 countries completed a questionnaire survey with 109 queries on the eye treatment planning system (TPS), hardware/software equipment, image acquisition/registration, patient positioning, eye surveillance, beam delivery, quality assurance (QA), clinical management, and workflow. RESULTS Worldwide, 28,891 eye patients were treated with protons at the 10 centers as of the end of 2014. Most centers treated a vast number of ocular patients (1729 to 6369). Three centers treated fewer than 200 ocular patients. Most commonly, the centers treated uveal melanoma (UM) and other primary ocular malignancies, benign ocular tumors, conjunctival lesions, choroidal metastases, and retinoblastomas. The UM dose fractionation was generally within a standard range, whereas dosing for other ocular conditions was not standardized. The majority (80%) of centers used in common a specific ocular TPS. Variability existed in imaging registration, with magnetic resonance imaging (MRI) rarely being used in routine planning (20%). Increased patient to full-time equivalent ratios were observed by higher accruing centers (P=.0161). Generally, ophthalmologists followed up the post-radiation therapy patients, though in 40% of centers radiation oncologists also followed up the patients. Seven centers had a prospective outcomes database. All centers used a cyclotron to accelerate protons with dedicated horizontal beam lines only. QA checks (range, modulation) varied substantially across centers. CONCLUSIONS The first worldwide multi-institutional ophthalmic proton therapy survey of the clinical and technical approach shows areas of substantial overlap and areas of progress needed to achieve sustainable and systematic management. Future international efforts include research and development for imaging and planning software upgrades, increased use of MRI, development of clinical protocols, systematic patient-centered data acquisition, and publishing guidelines on QA, staffing, treatment, and follow-up parameters by dedicated ocular programs to ensure the highest level of care for ocular patients.

Monte Carlo modeling of HD120 multileaf collimator on Varian TrueBeam linear accelerator for verification of 6X and 6X FFF VMAT SABR treatment plans

A Monte Carlo (MC) validation of the vendor‐supplied Varian TrueBeam 6 MV flattened (6X) phase‐space file and the first implementation of the Siebers‐Keall MC MLC model as applied to the HD120 MLC (for 6X flat and 6X flattening filterfree (6X FFF) beams) are described. The MC model is validated in the context of VMAT patient‐specific quality assurance. The Monte Carlo commissioning process involves: 1) validating the calculated open‐field percentage depth doses (PDDs), profiles, and output factors (OF), 2) adapting the Siebers‐Keall MLC model to match the new HD120‐MLC geometry and material composition, 3) determining the absolute dose conversion factor for the MC calculation, and 4) validating this entire linac/MLC in the context of dose calculation verification for clinical VMAT plans. MC PDDs for the 6X beams agree with the measured data to within 2.0% for field sizes ranging from 2 × 2 to 40 × 40 cm2. Measured and MC profiles show agreement in the 50% field width and the 80%‐20% penumbra region to within 1.3 mm for all square field sizes. MC OFs for the 2 to 40 cm2 square fields agree with measurement to within 1.6%. Verification of VMAT SABR lung, liver, and vertebra plans demonstrate that measured and MC ion chamber doses agree within 0.6% for the 6X beam and within 2.0% for the 6X FFF beam. A 3D gamma factor analysis demonstrates that for the 6X beam, > 99% of voxels meet the pass criteria (3%/3 mm). For the 6X FFF beam, > 94% of voxels meet this criteria. The TrueBeam accelerator delivering 6X and 6X FFF beams with the HD120 MLC can be modeled in Monte Carlo to provide an independent 3D dose calculation for clinical VMAT plans. This quality assurance tool has been used clinically to verify over 140 6X and 16 6X FFF TrueBeam treatment plans. PACS number: 87.55.K‐

The use of modified single pencil beam dose kernels to improve IMRT dose calculation accuracy

Intensity modulated radiation therapy (IMRT) is used to deliver highly conformal radiation doses to tumors while sparing nearby sensitive tissues. Discrepancies between calculated and measured dose distributions have been reported for regions of high dose gradients corresponding to complex radiation fluence patterns. For the single pencil beam convolution dose calculation algorithm, the ability to resolve areas of high dose structure is partly related to the shape of the pencil beam dose kernel (similar to how a photon detectors point spread function relates to imaging resolution). Improvements in dose calculation accuracy have been reported when the treatment planning system (TPS) is recommissioned using high-resolution measurement data as input. This study proposes to improve the dose calculation accuracy for IMRT planning by modifying clinical dose kernel shapes already present in the TPS, thus avoiding the need to reacquire higher resolution commissioning data. The in-house optimization program minimizes a cost-function based on a two-dimensional composite dose subtraction/distance-to-agreement (gamma) analysis. The final modified kernel shapes are reintroduced into the treatment planning system and improvements to the dose calcula tion accuracy for complex IMRT dose distributions evaluated. The central kernel value (radius =0 cm) has the largest effect on the dose calculation resolution and is the focus of this study.

Monte Carlo validation of the TrueBeam 10XFFF phase–space files for applications in lung SABR

PURPOSE To establish the clinical acceptability of universal Monte Carlo phase-space data for the 10XFFF (flattening filter free) photon beam on the Varian TrueBeam Linac, including previously unreported data for small fields, output factors, and inhomogeneous media. The study was particularly aimed at confirming the suitability for use in simulations of lung stereotactic ablative radiotherapy treatment plans. METHODS Monte Carlo calculated percent depth doses (PDDs), transverse profiles, and output factors for the TrueBeam 10 MV FFF beam using generic phase-space data that have been released by the Varian MC research team were compared with in-house measurements and published data from multiple institutions (ten Linacs from eight different institutions). BEAMnrc was used to create field size specific phase-spaces located underneath the jaws. Doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Particular attention was paid to small fields (down to 1 × 1 cm(2)) and dose per pulse effects on dosimeter response for high dose rate 10XFFF beams. Ion chamber measurements were corrected for changes in ion collection efficiency (P(ion)) with increasing dose per pulse. MC and ECLIPSE ANISOTROPIC ANALYTICAL ALGORITHM (AAA) calculated PDDs were compared to Gafchromic film measurement in inhomogeneous media (water, bone, lung). RESULTS Measured data from all machines agreed with Monte Carlo simulations within 1.0% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2) in a homogeneous water phantom. Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. For all the field sizes considered, the agreement between their measured and calculated output factors was within 1.1%. Monte Carlo results for dose to water at water/bone, bone/lung, and lung/water interfaces as well as within lung agree with film measurements to within 2.8% for 10 × 10 and 3 × 3 cm(2) field sizes. This represents a significant improvement over the performance of the ECLIPSE AAA. CONCLUSIONS The 10XFFF phase-space data offered by the Varian Monte Carlo research team have been validated for clinical use using measured, interinstitutional beam data in water and with film dosimetry in inhomogeneous media.

Characterization of a terbium activated gadolinium oxysulfide plastic optical fibre sensor in photons and protons

A characterization study was carried out to determine if a novel, millimetre sized Terbium-activated Gadolinium Oxysulfide optical fibre detector has potential for future use in dosimetry. 1.25 MeV photons from Co-60 decay and 74 MeV protons from the TRIUMF facility studies were used. Dose response, field size response and Cerenkov contributions, as well as raw and spread-out Bragg peak depth doses were investigated.

Poster - 37: Pre-clinical geometric, dosimetric and timing assessment of head and neck OARs using an in-house atlas-based auto-segmentation (ABAS) tool

Purpose: This study aims to validate the geometric and dosimetric performance of our in-house ABAS tool using manual inter-observer variation as benchmark data. Materials and methods: An in-house ABAS constructed, using MIM MaestroTM version 6.5, from 36 previously treated head and neck cases. 15 OARs of eight nasopharynx patients were segmented via three observers from the same institution and using the in-house ABAS. Percentage of volume differences (ΔV%,), degree of overlap (DICE), distance-to-agreement and the standard deviation of absolute dose difference ΔD SD(Gy) among the eight cases [ΔDmax SD(Gy) for serial organs and ΔDmean SD(Gy) for parallel organs] were compared between the manual segmentation and original ABAS contours for each OAR. Results: The geometric results indicated that ABAS ΔV% was within 1SD from the manual segmentation. DICE showed that manual segmentation marginally outperformed ABAS in the delineation of parotid glands, submandibular glands, and laryngopharynx, but ABSA performed as well as the manual segmentation or better for all other structures. The distance-to-agreement was <1.5 cm for 87% of ABAS structures. From a dosimetric perspective, only ABAS salivary glands demonstrated higher ΔD SD(Gy) compared to manual segmentation. For all other structures ABAS results were similar or better than the manual segmentation. The average time to segment a complete H&N OAR set was <3minutes versus 30minutes for manual segmentation. Conclusion: ABAS is a practical time-saving tool. This study indicated up 90% time-saving of the operator time per case. ABAS geometric and dosimetric results were within ±1SD inter-observer variation for most of the structures.