Nikos Papanikolaou

Investigation of the convolution method for polyenergetic spectra

The distribution of absolute dose per unit fluence from polyenergetic photon beams impinging upon a water phantom was calculated using two convolution approaches that properly account for beam hardening effects. Dose deposition kernels calculated previously using the EGS4 Monte Carlo code are convolved with the primary terma to give the dose for monoenergetic photon beams of energies ranging from 100 kev to 50 MeV. A polyenergetic dose distribution is composed of separately calculated monoenergetic components, which are appropriately weighted with the fluence spectrum to yield the polyenergetic dose distribution. Alternatively, a single convolution for the polyenergetic beam is considered, where a composite polyenergetic kernel is convolved with the respective polyenergetic terma. The effects of the polyenergetic kernel variance due to beam hardening as well as the effect of tilting the kernels for a diverging beam geometry were also examined. The depth dose data produced using the two proposed methods were compared with measured data and Monte Carlo simulations and showed good agreement.

Dosimetric considerations for patients with HIP prostheses undergoing pelvic irradiation. Report of the AAPM Radiation Therapy Committee Task Group 63.

This document is the report of a task group of the Radiation Therapy Committee of the AAPM and has been prepared primarily to advise hospital physicists involved in external beam treatment of patients with pelvic malignancies who have high atomic number (Z) hip prostheses. The purpose of the report is to make the radiation oncology community aware of the problems arising from the presence of these devices in the radiation beam, to quantify the dose perturbations they cause, and, finally, to provide recommendations for treatment planning and delivery. Some of the data and recommendations are also applicable to patients having implanted high-Z prosthetic devices such as pins, humeral head replacements. The scientific understanding and methodology of clinical dosimetry for these situations is still incomplete. This report is intended to reflect the current state of scientific understanding and technical methodology in clinical dosimetry for radiation oncology patients with high-Z hip prostheses.

Calculation of portal dose using the convolution/superposition method.

The convolution/superposition method was used to predict the dose throughout an extended volume, which includes a phantom and a portal imaging device. From the calculated dose volume, the dose delivered in the portal image plane was extracted and compared to a portal dose image. This comparison aids in verifying the beam configuration or patient setup after delivery of the radiation. The phantoms used to test the accuracy of this method include a solid water cube, a Nuclear Associates CT phantom, and an Alderson Rando thorax phantom. The dose distribution in the image plane was measured with film and an electronic portal imaging device in each case. The calculated portal dose images were within 4% of the measured images for most voxels in the central portion of the field for all of the extended volumes. The convolution/superposition method also enables the determination of the scatter and primary dose contributions using the particular dose deposition kernels for each contribution. The ratio of primary dose to total dose was used to extract the primary dose from the detected portal image, which enhances the megavoltage portal images by removing scatter blurring. By also predicting the primary energy fluence, we can find the ratio of computed primary energy fluence to total dose. Multiplying this ratio by the measured dose image estimates the relative primary energy fluence at the portal imager. The image of primary energy fluence possesses higher contrast and may be used for further quantitative image processing and dose modeling.

Prospective randomized double-blind pilot study of site-specific consensus atlas implementation for rectal cancer target volume delineation in the cooperative group setting

PURPOSE Variations in target volume delineation represent a significant hurdle in clinical trials involving conformal radiotherapy. We sought to determine the effect of a consensus guideline-based visual atlas on contouring the target volumes. METHODS AND MATERIALS A representative case was contoured (Scan 1) by 14 physician observers and a reference expert with and without target volume delineation instructions derived from a proposed rectal cancer clinical trial involving conformal radiotherapy. The gross tumor volume (GTV), and two clinical target volumes (CTVA, including the internal iliac, presacral, and perirectal nodes, and CTVB, which included the external iliac nodes) were contoured. The observers were randomly assigned to receipt (Group A) or nonreceipt (Group B) of a consensus guideline and atlas for anorectal cancers and then instructed to recontour the same case/images (Scan 2). Observer variation was analyzed volumetrically using the conformation number (CN, where CN = 1 equals total agreement). RESULTS Of 14 evaluable contour sets (1 expert and 7 Group A and 6 Group B observers), greater agreement was found for the GTV (mean CN, 0.75) than for the CTVs (mean CN, 0.46-0.65). Atlas exposure for Group A led to significantly increased interobserver agreement for CTVA (mean initial CN, 0.68, after atlas use, 0.76; p = .03) and increased agreement with the expert reference (initial mean CN, 0.58; after atlas use, 0.69; p = .02). For the GTV and CTVB, neither the interobserver nor the expert agreement was altered after atlas exposure. CONCLUSION Consensus guideline atlas implementation resulted in a detectable difference in interobserver agreement and a greater approximation of expert volumes for the CTVA but not for the GTV or CTVB in the specified case. Visual atlas inclusion should be considered as a feature in future clinical trials incorporating conformal RT.

Treatment planning and delivery of IMRT using 6 and 18 MV photon beams without flattening filter

In light of the increasing use of intensity modulated radiation therapy (IMRT) in modern radiotherapy practice, the use of a flattening filter may no longer be necessary. Commissioning data have been measured for a Varian 23EX linear accelerator with 6 and 18 MV photon energies without a flattening filter. Measurements collected for the commissioning of the linac included percent depth dose curves and profiles for field sizes ranging from 2 x 2 to 40 x 40 cm(2) as defined by the jaws and multileaf collimator. Machine total scatter factors were measured and calculated. Measurements were used to model the unflattened beams with the Pinnacle(3) treatment planning system. IMRT plans for prostate, lung, brain and head and neck cancer cases were generated using the flattening filter and flattening filter-free beams. From our results, no difference in the quality of the treatment plans between the flat and unflattened photon beams was noted. There was however a significant decrease in the number of monitor units required for unflattened beam treatment plans due to the increase in linac output-approximately two times and four times higher for the 6 and 18 MV, respectively.

Flattening filter‐free accelerators: a report from the AAPM Therapy Emerging Technology Assessment Work Group

This report describes the current state of flattening filter‐free (FFF) radiotherapy beams implemented on conventional linear accelerators, and is aimed primarily at practicing medical physicists. The Therapy Emerging Technology Assessment Work Group of the American Association of Physicists in Medicine (AAPM) formed a writing group to assess FFF technology. The published literature on FFF technology was reviewed, along with technical specifications provided by vendors. Based on this information, supplemented by the clinical experience of the group members, consensus guidelines and recommendations for implementation of FFF technology were developed. Areas in need of further investigation were identified. Removing the flattening filter increases beam intensity, especially near the central axis. Increased intensity reduces treatment time, especially for high‐dose stereotactic radiotherapy/radiosurgery (SRT/SRS). Furthermore, removing the flattening filter reduces out‐of‐field dose and improves beam modeling accuracy. FFF beams are advantageous for small field (e.g., SRS) treatments and are appropriate for intensity‐modulated radiotherapy (IMRT). For conventional 3D radiotherapy of large targets, FFF beams may be disadvantageous compared to flattened beams because of the heterogeneity of FFF beam across the target (unless modulation is employed). For any application, the nonflat beam characteristics and substantially higher dose rates require consideration during the commissioning and quality assurance processes relative to flattened beams, and the appropriate clinical use of the technology needs to be identified. Consideration also needs to be given to these unique characteristics when undertaking facility planning. Several areas still warrant further research and development. Recommendations pertinent to FFF technology, including acceptance testing, commissioning, quality assurance, radiation safety, and facility planning, are presented. Examples of clinical applications are provided. Several of the areas in which future research and development are needed are also indicated. PACS number: 87.53.‐j, 87.53.Bn, 87.53.Ly, 87.55.‐x, 87.55.N‐, 87.56.bc

Analysis of the sources of uncertainty for EDR2 film-based IMRT quality assurance

In our institution, patient‐specific quality assurance (QA) for intensity‐modulated radiation therapy (IMRT) is usually performed by measuring the dose to a point using an ion chamber and by measuring the dose to a plane using film. In order to perform absolute dose comparison measurements using film, an accurate calibration curve should be used. In this paper, we investigate the film response curve uncertainty factors, including film batch differences, film processor temperature effect, film digitization, and treatment unit. In addition, we reviewed 50 patient‐specific IMRT QA procedures performed in our institution in order to quantify the sources of error in film‐based dosimetry. Our study showed that the EDR2 film dosimetry can be done with less than 3% uncertainty. The EDR2 film response was not affected by the choice of treatment unit provided the nominal energy was the same. This investigation of the different sources of uncertainties in the film calibration procedure can provide a better understanding of the film‐based dosimetry and can improve quality control for IMRT QA. PACS numbers: 87.86.Cd, 87.53.Xd, 87.57.Nk

Treatment plan comparison between helical tomotherapy and MLC-based IMRT using radiobiological measures

The rapid implementation of advanced treatment planning and delivery technologies for radiation therapy has brought new challenges in evaluating the most effective treatment modality. Intensity-modulated radiotherapy (IMRT) using multi-leaf collimators (MLC) and helical tomotherapy (HT) are becoming popular modes of treatment delivery and their application and effectiveness continues to be investigated. Presently, there are several treatment planning systems (TPS) that can generate and optimize IMRT plans based on user-defined objective functions for the internal target volume (ITV) and organs at risk (OAR). However, the radiobiological parameters of the different tumours and normal tissues are typically not taken into account during dose prescription and optimization of a treatment plan or during plan evaluation. The suitability of a treatment plan is typically decided based on dosimetric criteria such as dose-volume histograms (DVH), maximum, minimum, mean and standard deviation of the dose distribution. For a more comprehensive treatment plan evaluation, the biologically effective uniform dose (D) is applied together with the complication-free tumour control probability (P(+)). Its utilization is demonstrated using three clinical cases that were planned with two different forms of IMRT. In this study, three different cancer types at different anatomical sites were investigated: head and neck, lung and prostate cancers. For each cancer type, a linac MLC-based step-and-shoot IMRT plan and a HT plan were developed. The MLC-based IMRT treatment plans were developed on the Philips treatment-planning platform, using the Pinnacle 7.6 software release. For the tomotherapy HiArt plans, the dedicated tomotherapy treatment planning station was used, running version 2.1.2. By using D as the common prescription point of the treatment plans and plotting the tissue response probabilities versus D for a range of prescription doses, a number of plan trials can be compared based on radiobiological measures. The applied plan evaluation method shows that in the head and neck cancer case the HT treatment gives better results than MLC-based IMRT in terms of expected clinical outcome P(+) of 62.2% and 46.0%, D to the ITV of 72.3 Gy and 70.7 Gy, respectively). In the lung cancer and prostate cancer cases, the MLC-based IMRT plans are better over the clinically useful dose prescription range. For the lung cancer case, the HT and MLC-based IMRT plans give a P(+) of 66.9% and 72.9%, D to the ITV of 64.0 Gy and 66.9 Gy, respectively. Similarly, for the prostate cancer case, the two radiation modalities give a P(+) of 68.7% and 72.2%, D to the ITV of 86.0 Gy and 85.9 Gy, respectively. If a higher risk of complications (higher than 5%) could be allowed, the complication-free tumour control could increase by over 40%, 2% and 30% compared to the initial dose prescription for the three cancer cases, respectively. Both MLC-based IMRT and HT can encompass the often-large ITV required while they minimize the volume of the organs at risk receiving high doses. Radiobiological evaluation of treatment plans may provide an improved correlation of the delivered treatment with the clinical outcome by taking into account the dose-response characteristics of the irradiated targets and normal tissues. There may exist clinical cases, which may look dosimetrically similar but in radiobiological terms may be quite different. In such situations, traditional dose-based evaluation tools can be complemented by the use of P(+)--D diagrams to effectively evaluate and compare treatment plans.

Measurements of the electron dose distribution near inhomogeneities using a plastic scintillation detector

PURPOSE Accurate measurement of the electron dose distribution near an inhomogeneity is difficult with traditional dosimeters which themselves perturb the electron field. We tested the performance of a new high resolution, water-equivalent plastic scintillation detector which has ideal properties for this application. METHODS AND MATERIALS A plastic scintillation detector with a 1 mm diameter, 3 mm long cylindrical sensitive volume was used to measure the dose distributions behind standard benchmark inhomogeneities in water phantoms. The plastic scintillator material is more water equivalent than polystyrene in terms of its mass collision stopping power and mass scattering power. Measurements were performed for beams of electrons having initial energies of 6 and 18 MeV at depths from 0.2-4.2 cm behind the inhomogeneities. RESULTS The detector reveals hot and cold spots behind heterogeneities at resolutions equivalent to typical film digitizer spot sizes. Plots of the dose distributions behind air, aluminum, lead, and formulations for cortical and inner bone-equivalent materials are presented. CONCLUSION The plastic scintillation detector is suited for measuring the electron dose distribution near an inhomogeneity.

Characterization of a novel 2D array dosimeter for patient‐specific quality assurance with volumetric arc therapy

PURPOSE In this study, the authors are evaluating a new, commercially available 2D array that offers 3D dose reconstruction for patient specific intensity modulated radiation therapy quality assurance (IMRT QA). METHODS The OCTAVIUS 4D system and its accompanying software (VERISOFT) by PTW were evaluated for the accuracy of the dose reconstruction for patient specific pretreatment IMRT QA. OCTAVIUS 4D measures the dose plane at the linac isocenter as the phantom rotates synchronously with the gantry, maintaining perpendicularity with the beam, by means of an inclinometer and a motor. The measurements collected during a volumetric modulated arc therapy delivery (VMAT) are reconstructed into a 3D dose volume. The VERISOFT application is used to perform the analysis, by comparing the reconstructed dose against the 3D dose matrix from the treatment planning system (TPS) that is computed for the same geometry and beam arrangement as that of the measurement. In this study, the authors evaluated the 3D dose reconstruction algorithm of this new system using a series of tests. Using the Octavius 4D phantom as the patient, dose distributions for various field sizes, beam orientations, shapes, and combination of fields were calculated using the Pinnacle3, TPS, and the respective DICOMRT dose was exported to the VERISOFT analysis software. Measurements were obtained by delivering the test treatment plans and comparisons were made based on gamma index, dose profiles, and isodose distribution analysis. In addition, output factors were measured and the dose linearity of the array was assessed. Those measurements were compared against measurements in water using a single, calibrated ionization chamber as well as calculations from Pinnacle for the same delivery geometries. RESULTS The number of voxels that met the 3%/3 mm criteria for the volumetric 3D gamma index analysis ranged from 92.3% to 98.9% for all the patient plans that the authors evaluated. 2D gamma analysis in the axial, sagittal, and coronal planes produced similar results to those in the 3D gamma analysis. The new detector system does not require an angular dependence correction because it rotates in synchrony with the gantry and the detector array maintains a constant SAD while always perpendicular to the beam axis. Output factors were within 2% when compared to ionization chamber measurements and Pinnacle calculations. Similar agreement was observed when testing the MU linearity (for MU values above 2) as well as dose rate effect. CONCLUSIONS The OCTAVIUS 4D system has some unique characteristics that can potentially improve the patient specific pretreatment IMRT QA data collection and analysis. The ability of the software to reconstruct from the measurements the true 3D dose distribution in the phantom, provides a unique perspective for the medical physicist that evaluates a patients QA plan.