Todd Pawlicki

Clinical implementation of a Monte Carlo treatment planning system

The purpose of this study was to implement the Monte Carlo method for clinical radiotherapy dose calculations. We used the EGS4/BEAM code to obtain the phase-space data for 6-20 MeV electron beams and 4, 6, and 15 MV photon beams for Varian Clinac 1800, 2100C, and 2300CD accelerators. A multiple-source model was used to reconstruct the phase-space data for both electron and photon beams, which retained the accuracy of the Monte Carlo beam data. The multiple-source model reduced the phase-space data storage requirement by a factor of 1000 and the accelerator simulation time by a factor of 10 or more. Agreement within 2% was achieved between the Monte Carlo calculations and measurements of the dose distributions in homogeneous and heterogeneous phantoms for various field sizes, source-surface distances, and beam modulations. The Monte Carlo calculated electron output factors were within 2% of the measured values for various treatment fields while the heterogeneity correction factors for various lung and bone phantoms were within 1% for photon beams and within 2% for electron beams. The EGS4/DOSXYZ Monte Carlo code was used for phantom and patient dose calculations. The results were compared to the dose distributions produced by a conventional treatment planning system and an intensity-modulated radiotherapy inverse-planning system. Significant differences (>5% in dose and >5 mm shift in isodose lines) were found between Monte Carlo calculations and the analytical calculations implemented in the commercial systems. Treatment sites showing the largest dose differences were for head and neck, lung, and breast cases.

A Monte Carlo dose calculation tool for radiotherapy treatment planning

A Monte Carlo user code, MCDOSE, has been developed for radiotherapy treatment planning (RTP) dose calculations. MCDOSE is designed as a dose calculation module suitable for adaptation to host RTP systems. MCDOSE can be used for both conventional photon/electron beam calculation and intensity modulated radiotherapy (IMRT) treatment planning. MCDOSE uses a multiple-source model to reconstruct the treatment beam phase space. Based on Monte Carlo simulated or measured beam data acquired during commissioning, source-model parameters are adjusted through an automated procedure. Beam modifiers such as jaws, physical and dynamic wedges, compensators, blocks, electron cut-outs and bolus are simulated by MCDOSE together with a 3D rectilinear patient geometry model built from CT data. Dose distributions calculated using MCDOSE agreed well with those calculated by the EGS4/DOSXYZ code using different beam set-ups and beam modifiers. Heterogeneity correction factors for layered-lung or layered-bone phantoms as calculated by both codes were consistent with measured data to within 1%. The effect of energy cut-offs for particle transport was investigated. Variance reduction techniques were implemented in MCDOSE to achieve a speedup factor of 10-30 compared to DOSXYZ.

An evidence based review of proton beam therapy: The report of ASTRO’s emerging technology committee

Proton beam therapy (PBT) is a novel method for treating malignant disease with radiotherapy. The purpose of this work was to evaluate the state of the science of PBT and arrive at a recommendation for the use of PBT. The emerging technology committee of the American Society of Radiation Oncology (ASTRO) routinely evaluates new modalities in radiotherapy and assesses the published evidence to determine recommendations for the society as a whole. In 2007, a Proton Task Force was assembled to evaluate the state of the art of PBT. This report reflects evidence collected up to November 2009. Data was reviewed for PBT in central nervous system tumors, gastrointestinal malignancies, lung, head and neck, prostate, and pediatric tumors. Current data do not provide sufficient evidence to recommend PBT in lung cancer, head and neck cancer, GI malignancies, and pediatric non-CNS malignancies. In hepatocellular carcinoma and prostate cancer and there is evidence for the efficacy of PBT but no suggestion that it is superior to photon based approaches. In pediatric CNS malignancies PBT appears superior to photon approaches but more data is needed. In large ocular melanomas and chordomas, we believe that there is evidence for a benefit of PBT over photon approaches. PBT is an important new technology in radiotherapy. Current evidence provides a limited indication for PBT. More robust prospective clinical trials are needed to determine the appropriate clinical setting for PBT.

A quality assurance program for the on-board imager ®

To develop a quality assurance (QA) program for the On-Board Imager (OBI) system and to summarize the results of these QA tests over extended periods from multiple institutions. Both the radiographic and cone-beam computed tomography (CBCT) mode of operation have been evaluated. The QA programs from four institutions have been combined to generate a series of tests for evaluating the performance of the On-Board Imager. The combined QA program consists of three parts: (1) safety and functionality, (2) geometry, and (3) image quality. Safety and functionality tests evaluate the functionality of safety features and the clinical operation of the entire system during the tube warm-up. Geometry QA verifies the geometric accuracy and stability of the OBI/CBCT hardware/software. Image quality QA monitors spatial resolution and contrast sensitivity of the radiographic images. Image quality QA for CBCT includes tests for Hounsfield Unit (HU) linearity, HU uniformity, spatial linearity, and scan slice geometry, in addition. All safety and functionality tests passed on a daily basis. The average accuracy of the OBI isocenter was better than 1.5mm with a range of variation of less than 1mm over 8 months. The average accuracy of arm positions in the mechanical geometry QA was better than 1mm, with a range of variation of less than 1mm over 8 months. Measurements of other geometry QA tests showed stable results within tolerance throughout the test periods. Radiographic contrast sensitivity ranged between 2.2% and 3.2% and spatial resolution ranged between 1.25 and 1.6lp∕mm. Over four months the CBCT images showed stable spatial linearity, scan slice geometry, contrast resolution (1%; <7mm disk) and spatial resolution (>6lp∕cm). The HU linearity was within ±40HU for all measurements. By combining test methods from multiple institutions, we have developed a comprehensive, yet practical, set of QA tests for the OBI system. Use of the tests over extended periods show that the OBI system has reliable mechanical accuracy and stable image quality. Nevertheless, the tests have been useful in detecting performance deficits in the OBI system that needed recalibration. It is important that all tests are performed on a regular basis.

Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system

The purpose of this work was to use Monte Carlo simulations to verify the accuracy of the dose distributions from a commercial treatment planning optimization system (Corvus, Nomos Corp., Sewickley, PA) for intensity-modulated radiotherapy (IMRT). A Monte Carlo treatment planning system has been implemented clinically to improve and verify the accuracy of radiotherapy dose calculations. Further modifications to the system were made to compute the dose in a patient for multiple fixed-gantry IMRT fields. The dose distributions in the experimental phantoms and in the patients were calculated and used to verify the optimized treatment plans generated by the Corvus system. The Monte Carlo calculated IMRT dose distributions agreed with the measurements to within 2% of the maximum dose for all the beam energies and field sizes for both the homogeneous and heterogeneous phantoms. The dose distributions predicted by the Corvus system, which employs a finite-size pencil beam (FSPB) algorithm, agreed with the Monte Carlo simulations and measurements to within 4% in a cylindrical water phantom with various hypothetical target shapes. Discrepancies of more than 5% (relative to the prescribed target dose) in the target region and over 20% in the critical structures were found in some IMRT patient calculations. The FSPB algorithm as implemented in the Corvus system is adequate for homogeneous phantoms (such as prostate) but may result in significant under or over-estimation of the dose in some cases involving heterogeneities such as the air-tissue, lung-tissue and tissue-bone interfaces.

Photon beam characterization and modelling for Monte Carlo treatment planning.

Photon beams of 4, 6 and 15 MV from Varian Clinac 2100C and 2300C/D accelerators were simulated using the EGS4/BEAM code system. The accelerators were modelled as a combination of component modules (CMs) consisting of a target, primary collimator, exit window, flattening filter, monitor chamber, secondary collimator, ring collimator, photon jaws and protection window. A full phase space file was scored directly above the upper photon jaws and analysed using beam data processing software, BEAMDP, to derive the beam characteristics, such as planar fluence, angular distribution, energy spectrum and the fractional contributions of each individual CM. A multiple-source model has been further developed to reconstruct the original phase space. Separate sources were created with accurate source intensity, energy, fluence and angular distributions for the target, primary collimator and flattening filter. Good agreement (within 2%) between the Monte Carlo calculations with the source model and those with the original phase space was achieved in the dose distributions for field sizes of 4 cm x 4 cm to 40 cm x 40 cm at source surface distances (SSDs) of 80-120 cm. The dose distributions in lung and bone heterogeneous phantoms have also been found to be in good agreement (within 2%) for 4, 6 and 15 MV photon beams for various field sizes between the Monte Carlo calculations with the source model and those with the original phase space.

Quality and safety considerations in stereotactic radiosurgery and stereotactic body radiation therapy: Executive summary

In summary, SRS and SBRT require a team-based approach, staffed by appropriately trained and credentialed specialists. SRS and SBRT training should become a required part of radiation oncology residency training and of Accreditation of Medical Physics Educational Programs accredited clinical medical physics training. SRS and SBRT require significant resources in personnel, specialized technology, and implementation time. A thorough feasibility analysis of resources required to achieve the clinical and technical goals must be performed and discussed with all personnel, including medical center administration. Because various disease sites may have different clinical and technical requirements, feasibility and planning discussions are needed prior to undertaking new disease sites. Treatment of SRS/SBRT patients should adhere to established national guidelines. Acceptance and commissioning protocols and tests must be developed to explore in detail every aspect of the individual and integrated systems with the goal of ensuring safe and effective operation. A comprehensive quality assurance program, encompassing all clinical, technical, and patient-specific treatment aspects, must be developed to ensure SRS and SBRT are performed in a safe and effective manner. Patient safety in radiation therapy is everyones responsibility. Professional organizations, regulators, vendors, and end-users must demonstrate a clear commitment to working closely together to ensure the highest levels of safety and efficacy in stereotactic radiosurgery and stereotactic body radiation therapy.

A comparative dosimetric study on tangential photon beams, intensity-modulated radiation therapy (IMRT) and modulated electron radiotherapy (MERT) for breast cancer treatment

Recently, energy- and intensity-modulated electron radiotherapy (MERT) has garnered a growing interest for the treatment of superficial targets. In this work. we carried out a comparative dosimetry study to evaluate MERT, photon beam intensity-modulated radiation therapy (IMRT) and conventional tangential photon beams for the treatment of breast cancer. A Monte Carlo based treatment planning system has been investigated, which consists of a set of software tools to perform accurate dose calculation, treatment optimization, leaf sequencing and plan analysis. We have compared breast treatment plans generated using this home-grown treatment optimization and dose calculation software forthese treatment techniques. The MERT plans were planned with up to two gantry angles and four nominal energies (6, 9, 12 and 16 MeV). The tangential photon treatment plans were planned with 6 MV wedged photon beams. The IMRT plans were planned using both multiple-gantry 6 MV photon beams or two 6 MV tangential beams. Our results show that tangential IMRT can reduce the dose to the lung, heart and contralateral breast compared to conventional tangential wedged beams (up to 50% reduction in high dose volume or 5 Gy in the maximum dose). MERT can reduce the maximum dose to the lung by up to 20 Gy and to the heart by up to 35 Gy compared to conventional tangential wedged beams. Multiple beam angle IMRT can significantly reduce the maximum dose to the lung and heart (up to 20 Gy) but it induces low and medium doses to a large volume of normal tissues including lung, heart and contralateral breast. It is concluded that MERT has superior capabilities to achieve dose conformity both laterally and in the depth direction, which will be well suited for treating superficial targets such as breast cancer.

Validation of a Monte Carlo dose calculation tool for radiotherapy treatment planning

A new EGS4/PRESTA Monte Carlo user code, MCDOSE, has been developed as a routine dose calculation tool for radiotherapy treatment planning. It is suitable for both conventional and intensity modulated radiation therapy. Two important features of MCDOSE are the inclusion of beam modifiers in the patient simulation and the implementation of several variance reduction techniques. Before this tool can be used reliably for clinical dose calculation, it must be properly validated. The validation for beam modifiers has been performed by comparing the dose distributions calculated by MCDOSE and the well-benchmarked EGS4 user codes BEAM and DOSXYZ. Various beam modifiers were simulated. Good agreement in the dose distributions was observed. The differences in electron cutout factors between the results of MCDOSE and measurements were within 2%. The accuracy of MCDOSE with various variance reduction techniques was tested by comparing the dose distributions in different inhomogeneous phantoms with those calculated by DOSXYZ without variance reduction. The agreement was within 1.0%. Our results demonstrate that MCDOSE is accurate and efficient for routine dose calculation in radiotherapy treatment planning, with or without beam modifiers.

Energy- and intensity-modulated electron beams for radiotherapy

This work investigates the feasibility of optimizing energy- and intensity-modulated electron beams for radiation therapy. A multileaf collimator (MLC) specially designed for modulated electron radiotherapy (MERT) was investigated both experimentally and by Monte Carlo simulations. An inverse-planning system based on Monte Carlo dose calculations was developed to optimize electron beam energy and intensity to achieve dose conformity for target volumes near the surface. The results showed that an MLC with 5 mm leaf widths could produce complex field shapes for MERT. Electron intra- and inter-leaf leakage had negligible effects on the dose distributions delivered with the MLC, even at shallow depths. Focused leaf ends reduced the electron scattering contributions to the dose compared with straight leaf ends. As anticipated, moving the MLC position toward the patient surface reduced the penumbra significantly. There were significant differences in the beamlet distributions calculated by an analytic 3-D pencil beam algorithm and the Monte Carlo method. The Monte Carlo calculated beamlet distributions were essential to the accuracy of the MERT dose distribution in cases involving large air gaps, oblique incidence and heterogeneous treatment targets (at the tissue-bone and bone-lung interfaces). To demonstrate the potential of MERT for target dose coverage and normal tissue sparing for treatment of superficial targets, treatment plans for a hypothetical treatment were compared using photon beams and MERT.