Jiankui Yuan

Hypofractionation Regimens for Stereotactic Radiotherapy for Large Brain Tumors

PURPOSE To investigate equivalent regimens for hypofractionated stereotactic radiotherapy (HSRT) for brain tumor treatment and to provide dose-escalation guidance to maximize the tumor control within the normal brain tolerance. METHODS AND MATERIALS The linear-quadratic model, including the effect of nonuniform dose distributions, was used to evaluate the HSRT regimens. The alpha/beta ratio was estimated using the Gammaknife stereotactic radiosurgery (GKSRS) and whole-brain radiotherapy experience for large brain tumors. The HSRT regimens were derived using two methods: (1) an equivalent tumor control approach, which matches the whole-brain radiotherapy experience for many fractions and merges it with the GKSRS data for few fractions; and (2) a normal-tissue tolerance approach, which takes advantages of the dose conformity and fractionation of HSRT to approach the maximal dose tolerance of the normal brain. RESULTS A plausible alpha/beta ratio of 12 Gy for brain tumor and a volume parameter n of 0.23 for normal brain were derived from the GKSRS and whole-brain radiotherapy data. The HSRT prescription regimens for the isoeffect of tumor irradiation were calculated. The normal-brain equivalent uniform dose decreased as the number of fractions increased, because of the advantage of fractionation. The regimens for potential dose escalation of HSRT within the limits of normal-brain tolerance were derived. CONCLUSIONS The designed hypofractionated regimens could be used as a preliminary guide for HSRT dose prescription for large brain tumors to mimic the GKSRS experience and for dose escalation trials. Clinical studies are necessary to further tune the model parameters and validate these regimens.

Comparison of Ray Tracing and Monte Carlo Calculation Algorithms for Thoracic Spine Lesions Treated With CyberKnife-Based Stereotactic Body Radiation Therapy.

Stereotactic body radiation therapy (SBRT) is an emerging technology for the treatment of spinal metastases, although the dosimetric impact of the calculation method on spinal dose distribution is unknown. This study attempts to determine whether CyberKnife (CK)-based SBRT using a Ray Tracing (RyTc) algorithm is comparable dosimetrically to that of Monte Carlo (MC) for thoracic spinal lesions. Our institutional CK-based SBRT database for thoracic spinal lesions was queried and a cohort was generated. Patients were planned using RyTc and MC algorithms using the same beam angles and monitor units. Dose–volume histograms of the planning target volume (PTV), spinal cord, esophagus, and skin were generated, and dosimetric parameters were compared. There were 37 patients in the cohort. The average percentage volume of PTV covered by the prescribed dose with RyTc and MC algorithms was 91.1% and 80.4%, respectively (P < .001). The difference in average maximum spinal cord dose between RyTc and MC plans was significant (1126 vs 1084 cGy, P = .004), with the MC dose ranging from 18.7% below to 13.8% above the corresponding RyTc dose. A small reduction in maximum skin dose was also noted (P = .017), although no difference was seen in maximum esophageal dose (P = .15). Only PTVs smaller than 27 cm3 were found to correlate with large (>10%) changes in dose to 90% of the volume (P = .014), while no correlates with the average percentage volume of PTV covered by the prescribed dose were demonstrated. For thoracic spinal CK-based SBRT, RyTc computation may overestimate the MC calculated average percentage volume of PTV covered by the prescribed dose and have unpredictable effects on doses to organs at risk, particularly the spinal cord. In this setting, use of RyTc optimization should be limited and always verified with MC.

Development of a Monte Carlo model for treatment planning dose verification of the Leksell Gamma Knife Perfexion radiosurgery system

Detailed Monte Carlo (MC) modeling of the Leksell Gamma Knife (GK) Perfexion (PFX) collimator system is the only accurate ab initio approach appearing in the literature. As a different approach, in this work, we present a MC model based on film measurement. By adjusting the model parameters and fine-tuning the derived fluence map for each individual source to match the manufacturers ring output factors, we created a reasonable virtual source model for MC simulations to verify treatment planning dose for the GK PFX radiosurgery system. The MC simulation model was commissioned by simple single shots. Dose profiles and both ring and collimator output factors were compared with the treatment planning system (TPS). Good agreement was achieved for dose profiles especially for the region of plateau (<2%), while larger difference (<5%) came from the penumbra region. The maximum difference of the calculated output factor was within 0.7%. The model was further validated by a clinical test case. Good agreement was obtained. The DVHs for brainstem and the skull were almost identical and, for the target, the volume covered by the prescription (12.5 Gy to 50% isodose line) was 95.6% from MC calculation versus 100% from the TPS. PACS number(s): 87.55.dk.Detailed Monte Carlo (MC) modeling of the Leksell Gamma Knife (GK) Perfexion (PFX) collimator system is the only accurate ab initio approach appearing in the literature. As a different approach, in this work, we present a MC model based on film measurement. By adjusting the model parameters and fine‐tuning the derived fluence map for each individual source to match the manufacturers ring output factors, we created a reasonable virtual source model for MC simulations to verify treatment planning dose for the GK PFX radiosurgery system. The MC simulation model was commissioned by simple single shots. Dose profiles and both ring and collimator output factors were compared with the treatment planning system (TPS). Good agreement was achieved for dose profiles especially for the region of plateau (<2%), while larger difference (<5%) came from the penumbra region. The maximum difference of the calculated output factor was within 0.7%. The model was further validated by a clinical test case. Good agreement was obtained. The DVHs for brainstem and the skull were almost identical and, for the target, the volume covered by the prescription (12.5 Gy to 50% isodose line) was 95.6% from MC calculation versus 100% from the TPS. PACS number(s): 87.55.dk

Application of positron emission tomography/computed tomography in radiation treatment planning for head and neck cancers

18-fluorodeoxygluocose positron emission tomography/computed tomography ((18)FDG-PET/CT) provides significant information in multiple settings in the management of head and neck cancers (HNC). This article seeks to define the additional benefit of PET/CT as related to radiation treatment planning for squamous cell carcinomas (SCCs) of the head and neck through a review of relevant literature. By helping further define both primary and nodal volumes, radiation treatment planning can be improved using PET/CT. Special attention is paid to the independent benefit of PET/CT in targeting mucosal primaries as well as in detecting nodal metastases. The utility of PET/CT is also explored for treatment planning in the setting of SCC of unknown primary as PET/CT may help define a mucosal target volume by guiding biopsies for examination under anesthesia thus changing the treatment paradigm and limiting the extent of therapy. Implications of the use of PET/CT for proper target delineation in patients with artifact from dental procedures are discussed and the impact of dental artifact on CT-based PET attenuation correction is assessed. Finally, comment is made upon the role of PET/CT in the high-risk post-operative setting, particularly in the context of radiation dose escalation. Real case examples are used in these settings to elucidate the practical benefits of PET/CT as related to radiation treatment planning in HNCs.

Investigation of Nonuniform Dose Voxel Geometry in Monte Carlo Calculations

The purpose of this work is to investigate the efficacy of using multi-resolution nonuniform dose voxel geometry in Monte Carlo (MC) simulations. An in-house MC code based on the dose planning method MC code was developed in C++ to accommodate the nonuniform dose voxel geometry package since general purpose MC codes use their own coupled geometry packages. We devised the package in a manner that the entire calculation volume was first divided into a coarse mesh and then the coarse mesh was subdivided into nonuniform voxels with variable voxel sizes based on density difference. We name this approach as multi-resolution subdivision (MRS). It generates larger voxels in small density gradient regions and smaller voxels in large density gradient regions. To take into account the large dose gradients due to the beam penumbra, the nonuniform voxels can be further split using ray tracing starting from the beam edges. The accuracy of the implementation of the algorithm was verified by comparing with the data published by Rogers and Mohan. The discrepancy was found to be 1% to 2%, with a maximum of 3% at the interfaces. Two clinical cases were used to investigate the efficacy of nonuniform voxel geometry in the MC code. Applying our MRS approach, we started with the initial voxel size of 5 × 5 × 3 mm3, which was further divided into smaller voxels. The smallest voxel size was 1.25 × 1.25 × 3 mm3. We found that the simulation time per history for the nonuniform voxels is about 30% to 40% faster than the uniform fine voxels (1.25 × 1.25 × 3 mm3) while maintaining similar accuracy.

A Monte Carlo model and its commissioning for the Leksell Gamma Knife Perfexion radiosurgery system

Purpose To develop and commission a Monte Carlo (MC) simulation model for the Leksell Gamma Knife (LGK) Perfexion (PFX) radiosurgery system. Method We previously established a source model for MC simulations of the LGK PFX for the purpose of the treatment planning system (TPS) dose verification and plan evaluation. To make practical and effective use of the model in clinic, several issues need to be addressed. First, thorough commissioning procedures are needed to ensure the validity of the model parameters, such as the source‐to‐focus (STF) distance, the source solid angle. Second, an efficient source particle sampling method is required to facilitate dose calculations for multitarget and multishot configurations in patient treatment plans. Third, inseparably, it is interesting to know the dose difference between the two GK TPS algorithms (TMR and convolution) and the MC method in extreme heterogeneous cases resulting from the inhomogeneous effect. We report our recent development in addressing these issues. Phantoms with the frame fiducials were manually created in the format of DICOM CT image to eliminate the uncertainties associated with scanner artifacts and image registration. The created homogeneous phantom was used to calibrate the model parameters to match the output factors with the manufacturer provided data, and the heterogeneous phantom with multilayer materials was used to study the inhomogeneous effect. Results The agreement between the MC calculation and TPS was very good for the homogeneous spherical phantom. The difference of the full width at half maximum (FWHM) of the profiles was less than 1 mm except for the profile for 16 mm collimator along z‐axis (less than 2 mm). For the extreme heterogeneous test case, it was shown that the TMR algorithm can overestimate the target dose by up to 22% using the measure of dose volume parameter D95. The agreement between the MC method and the TPS convolution method was better (within 3.6%) for the target near the center of phantom, however, discrepancy (up to 10.7%) existed for the target close to the skull. The difference between the two TPS dose algorithms was about 11%. Conclusions Considerable dose difference may result from the effect of heterogeneity, such as in the regions of the air cavities and bones. As the MC method has been extensively used in conventional external beams, it is worthwhile for further investigation in applying the MC method to accurate dose planning in the new GK PFX radiosurgery platform.

Technical Note: An approach to building a Monte Carlo simulation model for a double scattering proton beam system

PURPOSE The purpose of this study was to demonstrate and develop a Monte Carlo (MC) simulation model for a passive double scattering compact proton therapy system based on limited information of the mechanical components. METHOD We built a virtual machine source model (VMSM) which included a detailed definition of each beam-modifying component in the nozzle. Conceptually, it is similar to the conventional virtual analytical source model (VASM), except that the numerical machine nozzle or beamline is constructed in the VMSM, whereas in the VASM analytical parameters characterizing the energy spectrum and source fluence distribution are sought. All major beam shaping components were included in the VMSM and the model simulates interactions of the beam with a rotating range modulation wheel (RMW) combined with the beam current modulation. The RMWs, the first and second scatterer in the system were generated and tuned to reproduce measurement data as closely as possible. To validate the model, we compared the percent depth dose curves, spread out Bragg peaks (SOBPs) and lateral profiles against measured commissioning beam data. RESULTS The agreement of beam range between the MC calculation and measurement was within 1 mm for all beam options. The distal-falloff length was in good agreement as well (<1 mm for the large and deep groups, <1.5 mm for the small group). Agreement to within 2.5 mm of measured SOBP widths was obtained for all MC calculations. For lateral profiles, differences were found to be less than 2 mm. CONCLUSIONS We demonstrated that with limited geometrical information it is possible to build an acceptable source model for MC simulations of a passive double scattering compact proton therapy system. The agreement between the measurements and the MC model provides validation for use of the model for further studies of the dosimetric effects in patient treatments.

Development and Validation of a Small Animal Immobilizer and Positioning System for the Study of Delivery of Intracranial and Extracranial Radiotherapy Using the Gamma Knife System

The purpose of this research is to establish a process of irradiating mice using the Gamma Knife as a versatile system for small animal irradiation and to validate accurate intracranial and extracranial dose delivery using this system. A stereotactic immobilization device was developed for small animals for the Gamma Knife head frame allowing for isocentric dose delivery. Intercranial positional reproducibility of a reference point from a primary reference animal was verified on an additional mouse. Extracranial positional reproducibility of the mouse aorta was verified using 3 mice. Accurate dose delivery was validated using film and thermoluminescent dosimeter measurements with a solid water phantom. Gamma Knife plans were developed to irradiate intracranial and extracranial targets. Mice were irradiated validating successful targeted radiation dose delivery. Intramouse positional variability of the right mandible reference point across 10 micro-computed tomography scans was 0.65 ± 0.48 mm. Intermouse positional reproducibility across 2 mice at the same reference point was 0.76 ± 0.46 mm. The accuracy of dose delivery was 0.67 ± 0.29 mm and 1.01 ± 0.43 mm in the coronal and sagittal planes, respectively. The planned dose delivered to a mouse phantom was 2 Gy at the 50% isodose with a measured thermoluminescent dosimeter dose of 2.9 ± 0.3 Gy. The phosphorylated form of member X of histone family H2A (γH2AX) staining of irradiated mouse brain and mouse aorta demonstrated adjacent tissue sparing. In conclusion, our system for preclinical studies of small animal irradiation using the Gamma Knife is able to accurately deliver intracranial and extracranial targeted focal radiation allowing for preclinical experiments studying focal radiation.

Experimental Validation of Monte Carlo Simulations Based on a Virtual Source Model for TomoTherapy in a RANDO Phantom

A virtual source model for Monte Carlo simulations of helical TomoTherapy has been developed previously by the authors. The purpose of this work is to perform experiments in an anthropomorphic (RANDO) phantom with the same order of complexity as in clinical treatments to validate the virtual source model to be used for quality assurance secondary check on TomoTherapy patient planning dose. Helical TomoTherapy involves complex delivery pattern with irregular beam apertures and couch movement during irradiation. Monte Carlo simulation, as the most accurate dose algorithm, is desirable in radiation dosimetry. Current Monte Carlo simulations for helical TomoTherapy adopt the full Monte Carlo model, which includes detailed modeling of individual machine component, and thus, large phase space files are required at different scoring planes. As an alternative approach, we developed a virtual source model without using the large phase space files for the patient dose calculations previously. In this work, we apply the simulation system to recompute the patient doses, which were generated by the treatment planning system in an anthropomorphic phantom to mimic the real patient treatments. We performed thermoluminescence dosimeter point dose and film measurements to compare with Monte Carlo results. Thermoluminescence dosimeter measurements show that the relative difference in both Monte Carlo and treatment planning system is within 3%, with the largest difference less than 5% for both the test plans. The film measurements demonstrated 85.7% and 98.4% passing rate using the 3 mm/3% acceptance criterion for the head and neck and lung cases, respectively. Over 95% passing rate is achieved if 4 mm/4% criterion is applied. For the dose–volume histograms, very good agreement is obtained between the Monte Carlo and treatment planning system method for both cases. The experimental results demonstrate that the virtual source model Monte Carlo system can be a viable option for the accurate dose calculation of helical TomoTherapy.

SU-F-T-467: A Cross-Checking Approach for Dosimetric Verification of Beam- Matched Elekta Linear Accelerators

PURPOSE To verify the similarity of the dosimetric characteristics between two Elekta linear accelerators (linacs) in order to treat patients interchangeably on these two machines without re-planning. METHODS To investigate the viability of matching the 6 MV flattened beam on an existing linac (Elekta Synergy with Agility head) with a recently installed new linca (Elekta Versa HD), percent depth doses (PDD), flatness and symmetry output factors were compared for both machines. To validate the beam matching among machines, we carried out two approaches to cross-check the dosimetrical equivalence: 1) the prior treatment plans were re-computed based on the newly built Versa HD treatment planning system (TPS) model without changing the beam control points; 2) The same plans were delivered on both machines and the radiation dose measurements on a MapCheck2 were compared with TPS calculations. Three VMAT plans (Head and neck, lung, and prostate) were used in the study. RESULTS The difference between the PDDs for 10×10 cm2 field at all depths was less than 0.8%. The difference of flatness and symmetry for 30×30 cm2 field was less than 0.8%, and the measured output factors varies by less than 1% for each field size ranging from 2×2 cm2 to 40×40 cm2 . For the same plans, the maximum difference of the two calculated dose distributions is 2% of prescription. For the QA measurements, the gamma index passing rates were above 99% for 3%/3mm criteria with 10% threshold for all three clinical plans. CONCLUSION A beam modality matching between two Elekta linacs is demonstrated with a cross-checking approach.