Flattening filter free VMAT for a stereotactic, single-dose of 30 Gy to lung lesion in a 15-min treatment slot.

Cone-beam CT-guided single dose of lung stereotactic body radiotherapy (SBRT) treatment with a flattening filter free (FFF) beam and volumetric modulated arc therapy (VMAT) is a safe and highly effective treatment modality for selective small lung lesions. Four-dimensional (4D) CT-based treatment plans were generated using advanced AcurosXB algorithm for heterogeneity corrections. 6X-FFF beam produced highly conformal radiosurgical dose distribution to the target and reduced lung SBRT fraction duration to less than 10 min for a single dose of 30 Gy, significantly improving patient comfort and clinic workflow. Early follow-up CT imaging results (mean, 8 months) show high local control rates (100%) with no acute lung or rib toxicity. Longer clinical follow-up in a larger patient cohort managed in this fashion is underway to further validate this treatment approach.

Potential reduction of lung dose via VMAT with jaw tracking in the treatment of single-isocenter/two-lesion lung SBRT.

PURPOSE/OBJECTIVES: Due to higher radiosensitivity, non-target normal tissue dose is a major concern in stereotactic body radiation therapy (SBRT) treatment. The aim of this report was to estimate the dosimetric impact, specifically the reduction of normal lung dose in the treatment of single-isocenter/two-lesion lung SBRT via volumetric modulated arc therapy with jaw tracking (JT-VMAT). MATERIALS/METHODS: Twelve patients with two peripherally located early-stage non-small-cell-lung cancer (NSCLC) lung lesions underwent single-isocenter highly conformal non-coplanar JT-VMAT SBRT treatment in our institution. The mean isocenter to tumors distance was 5.6 ± 1.9 (range 4.3-9.5) cm. The mean combined planning target volume (PTV) was 38.7 ± 22.7 (range 5.0-80.9) cc. A single isocenter was placed between the two lesions. Doses were 54 and 50 Gy in three and five fractions, respectively. Plans were optimized in Eclipse with AcurosXB algorithm utilizing jaw tracking options for the Truebeam with a 6 MV-FFF beam and standard 120 leaf millennium multi-leaf collimators. For comparison, the JT-VMAT plans were retrospectively re-computed utilizing identical beam geometry, objectives, and planning parameters, but without jaw tracking (no JT-VMAT). Both plans were normalized to receive the same target coverage. The conformity and heterogeneity indices, intermediate-dose spillage [D2cm , R50, Gradient Index (GI), Gradient Distance (GD)], organs at risks (OAR) doses including normal lung as well as modulation factor (MF) were compared for both plans. RESULTS: For similar target coverage, GI, R50, GD, as well as the normal lung V5, V10, V20, mean lung dose (MLD), and maximum dose received by 1000 cc of lungs were statistically significant. Normal lung doses were reduced by 8%-11% with JT-VMAT. Normal lung dose increased as a function of tumor distance from isocenter. For the other OAR, up to 1%-16% reduction of non-target doses were observed with JT-VMAT. The MF and beam-on time were similar for both plans, however, MF increased as a function of tumors distance, consequently, delivering higher dose to normal lungs. CONCLUSION: Utilizing jaw tracking options during optimization for single-isocenter/two-lesion lung SBRT VMAT plans reduced doses to the normal lung and other OAR, reduced intermediate-dose spillage and provided superior/similar target coverage. Application of jaw tracking did not affect delivery efficiency and provided excellent plan quality with similar MF and beam-on time. Jaw tracking is recommended for future clinical SBRT plan optimization.

Evaluation of plan quality and treatment efficiency for single-isocenter/two-lesion lung stereotactic body radiation therapy.

PURPOSE/OBJECTIVES: To evaluate the plan quality and treatment delivery efficiency of single-isocenter/two-lesions volumetric modulated arc therapy (VMAT) lung stereotactic body radiation therapy (SBRT). MATERIALS/METHODS: Eight consecutive patients with two peripherally located early stage nonsmall-cell-lung cancer (NSCLC) lung lesions underwent single-isocenter highly conformal noncoplanar VMAT SBRT treatment in our institution. A single-isocenter was placed between the two lesions. Doses were 54 or 50 Gy in 3 and 5 fractions respectively. Patients were treated every other day. Plans were calculated in Eclipse with AcurosXB algorithm and normalized to at least 95% of the planning target volume (PTV) receiving 100% of the prescribed dose. For comparison, two-isocenter plans (isocenter placed centrally in each target) were retrospectively created. Conformity indices (CIs), heterogeneity index (HI), gradient index (GI), gradient distance (GD), and D2cm were calculated. The normal lung V5, V10, V20, mean lung dose (MLD) and other organs at risk (OARs) doses were evaluated. Total number of monitor units (MUs), beam-on time, and patient-specific quality assurance (QA) results were recorded. RESULTS: The mean isocenter to tumor distance was 6.7 ± 2.3 cm. The mean combined PTV was 44.0 ± 23.4 cc. There was no clinically significant difference in CI, HI, GD, GI, D2cm , and V20 including most of the OARs between single-isocenter and two-isocenter lung SBRT plans, evaluated per RTOG guidelines. However, for single-isocenter plans as the distance between the lesions increased, the V5, V10, and MLD increased, marginally. The total number of MUs and beam-on time was reduced by a factor of 1.5 for a single-isocenter plan compared to a two-isocenter plan. The single-isocenter/two-lesions VMAT lung SBRT QA plans demonstrated an accurate dose delivery of 98.1 ± 3.2% for clinical gamma passing rate of 3%/3 mm. CONCLUSION: The SBRT treatment of two peripherally located lung lesions with a centrally placed single-isocenter was dosimetrically equivalent to two-isocenter plans. Faster treatment delivery for single-isocenter treatment can improve patient compliance and reduce the amount of intrafraction motion errors for well-suited patients.

A feasibility study using TomoDirect for craniospinal irradiation.

The feasibility of delivering craniospinal irradiation (CSI) with TomoDirect is investigated. A method is proposed to generate TomoDirect plans using standard three-dimensional (3D) beam arrangements on Tomotherapy with junctioning of these fields to minimize hot or cold spots at the cranial/spinal junction. These plans are evaluated and compared to a helical Tomotherapy and a three-dimensional conformal therapy (3D CRT) plan delivered on a conventional linear accelerator (linac) for CSI. The comparison shows that a TomoDirect plan with an overlap between the cranial and spinal fields might be preferable over Tomotherapy plans because of decreased low dose to large volumes of normal tissues outside of the planning target volume (PTV). Although the TomoDirect plans were not dosimetrically superior to a 3D CRT linac plan, the patient can be easily treated in the supine position, which is often more comfortable and efficient from an anesthesia standpoint. TomoDirect plans also have only one setup position which obviates the need for matching of fields and feathering of junctions, two issues encountered with conventional 3D CRT plans. TomoDirect plans can be delivered with comparable treatment times to conventional 3D plans and in shorter times than a Tomotherapy plan. In this paper, a method is proposed for creating TomoDirect craniospinal plans, and the dosimetric consequences for choosing different planning parameters are discussed.

Quality assurance of U.S.-guided external beam radiotherapy for prostate cancer: report of AAPM Task Group 154.

Task Group 154 (TG154) of the American Association of Physicists in Medicine (AAPM) was created to produce a guidance document for clinical medical physicists describing recommended quality assurance (QA) procedures for ultrasound (U.S.)-guided external beam radiotherapy localization. This report describes the relevant literature, state of the art, and briefly summarizes U.S. imaging physics. Simulation, treatment planning and treatment delivery considerations are presented in order to improve consistency and accuracy. User training is emphasized in the report and recommendations regarding peer review are included. A set of thorough, yet practical, QA procedures, frequencies, and tolerances are recommended. These encompass recommendations to ensure both spatial accuracy and image quality.

Statistical analysis of dose heterogeneity in circulating blood: implications for sequential methods of total body irradiation.

PURPOSE: Improvements in delivery techniques for total body irradiation (TBI) using Tomotherapy and intensity modulated radiation therapy have been proven feasible. Despite the promise of improved dose conformality, the application of these "sequential" techniques has been hampered by concerns over dose heterogeneity to circulating blood. The present study was conducted to provide quantitative evidence regarding the potential clinical impact of this heterogeneity. METHODS: Blood perfusion was modeled analytically as possessing linear, sinusoidal motion in the craniocaudal dimension. The average perfusion period for human circulation was estimated to be approximately 78 s. Sequential treatment delivery was modeled as a Gaussian-shaped dose cloud with a 10 cm length that traversed a 183 cm patient length at a uniform speed. Total dose to circulating blood voxels was calculated via numerical integration and normalized to 2 Gy per fraction. Dose statistics and equivalent uniform dose (EUD) were calculated for relevant treatment times, radiobiological parameters, blood perfusion rates, and fractionation schemes. The model was then refined to account for random dispersion superimposed onto the underlying periodic blood flow. Finally, a fully stochastic model was developed using binomial and trinomial probability distributions. These models allowed for the analysis of nonlinear sequential treatment modalities and treatment designs that incorporate deliberate organ sparing. RESULTS: The dose received by individual blood voxels exhibited asymmetric behavior that depended on the coherence among the blood velocity, circulation phase, and the spatiotemporal characteristics of the irradiation beam. Heterogeneity increased with the perfusion period and decreased with the treatment time. Notwithstanding, heterogeneity was less than +/- 10% for perfusion periods less than 150 s. The EUD was compromised for radiosensitive cells, long perfusion periods, and short treatment times. However, the EUD was unaffected (within 10%) for perfusion periods of less than 150 s or treatment times of 20 min or greater. Treatment over six fractions improved the EUD per fraction such that all parametric combinations resulted in unaffected EUD. The stochastic models confirmed these results. CONCLUSIONS: Dose heterogeneity in circulating blood cells is clinically acceptable for typical treatment times, perfusion rates, and cell types. Development of conformal, sequential TBI treatment techniques should not be withheld based on concerns over circulating blood dose heterogeneity.

Algorithmic determination of source orientations for the CivaSheet directional brachytherapy device.

PURPOSE: The CivaSheet device uses multiple directionally shielded Pd-103 CivaDot sources to produce a directional planar dose distribution. In postplanning, manually digitizing the 3D source orientation is challenging because the 3D vector must be digitized by using 2D displayed images. The aim of this study is to develop an algorithm that will automatically determine the direction of each CivaDot source based on the location of sources adjacent to it. METHODS AND MATERIALS: The algorithm determines the source direction by averaging the normal directions of multiple local planes established by the adjacent sources. The algorithm was tested on a manually constructed CivaSheet-like device that was CT scanned in known flat geometries and two known curved geometries. Algorithmically determined source directions were compared with the known directions. The algorithm was also used on a postplan for a gynecological pelvic sidewall tumor bed implant and compared against manual digitization of the source directions. RESULTS: For the flat and curved test geometries, the average angular difference between the algorithm determined and known orientation was 1.2° ± 0.8° (flat geometry), 1.7° ± 2.1° (curve about vertical axis), and 2.3° ± 2.4° (curve about horizontal axis). For the patient case, results showed on average a 23.1° ± 10.8° difference between the manual digitized orientation and the algorithm's predicted orientation. CONCLUSIONS: The algorithm calculates the source orientation with accuracy better than 2.3° for the controlled experiments. In addition to its accuracy, the algorithm produces consistent results and lessens the difficult challenge of orienting the partially shielded sources.

Impact of backscatter material thickness on the depth dose of orthovoltage irradiators for radiobiology research.

The orthovoltage x-ray energy frequently used in radiation research is prone to dosimetry errors due to insufficient backscatter conditions. In many radiobiology studies, especially for cell irradiations, precise dose calculation algorithms such as Convolution-Superposition or Monte Carlo are impractical and as such, less accurate hand calculation methods are used for dose estimation. These dose estimation methods typically assume full backscatter conditions. The purpose of this study is to demonstrate the magnitude of the dose error that results from insufficient backscatter, and to provide lookup tables to account this issue. The beam spectra of several widely used commercial systems (XRAD-225, XRAD-320, SARRP) were used in Monte Carlo (MC) simulations on a series of phantom setups to investigate the impact of varying backscatter conditions on dosimetry. The depth dose curves for different field sizes, water phantom thicknesses and beam qualities were generated. In addition, depth dependent backscatter factors for different field sizes and different beam qualities were calculated. It is demonstrated that as much as a 50% dose difference exists for different backscatter conditions at the beam qualities studied. The choice of cell dish size as well as other changes in the experiment setup can have more than 10% impact on the dose. The impact of backscatter is reduced with a decrease in field size. Further, the thickness needed to provide full backscatter can be approximated as being equal to the field size. It is imperative to ensure full backscatter conditions during system and dosimeter calibration, or to use the look-up table provided in this study.

A dosimetric comparison of non-coplanar IMRT versus Helical Tomotherapy for nasal cavity and paranasal sinus cancer.

PURPOSES: To determine if there are clinically significant differences between the dosimetry of sinus tumors delivered by non-coplanar LINAC-based IMRT techniques and Helical Tomotherapy (HT). HT is capable of delivering highly conformal and uniform target dosimetry. However, HT lacks non-coplanar capability, which is commonly used for linear accelerator-based IMRT for nasal cavity and paranasal sinus tumors. METHODS AND MATERIALS: We selected 10 patients with representative early and advanced nasal cavity and paranasal sinus malignancies treated with a preoperative dose of 50 Gy/25 fractions without coverage of the cervical lymphatics for dosimetric comparison. Each plan was independently optimized using either Corvus inverse treatment planning system, commissioned for a Varian 2300 CD linear accelerator with 1cm multileaf collimator (MLC) leaves, or the HT inverse treatment planning system. A non-coplanar seven field technique was used in all Corvus plans with five mid-sagittal fields and two anterior oblique fields as described by Claus et al. [F. Claus, W. De Gersem, C. De Wagter, et al., An implementation strategy for IMRT of ethmoid sinus cancer and bilateral sparing of the optic pathways, Int J Radiat Oncol Biol Phys 51 (2001) 318-331], whereas only coplanar beamlets were used in HT planning. Dose plans were compared using DVHs, the minimum PTV dose to 1cm3 of the PTV, a uniformity index of planned treatment volume (PTV), and a comprehensive quality index (CQI) based on the maximum dose to optical structures, parotids and the brainstem which were deemed as the most critical adjacent structures. RESULTS: Both planning systems showed comparable PTV dose coverage, but HT had significantly higher uniformity (p<0.01) inside the PTV. The CQI for all organs at risk were equivalent except ipsilateral lenses and eyes, which received statistically lower dose from HT plans (p<0.01). CONCLUSIONS: Overall HT provided equivalent or slightly better normal structure avoidance with a more uniform PTV dose for nasal cavity and paranasal sinus cancer treatment than non-coplanar LINAC-based IMRT. The disadvantage of coplanar geometry in HT is apparently counterbalanced by the larger number of fields.

A motion phantom study on helical tomotherapy: the dosimetric impacts of delivery technique and motion.

Helical tomotherapy (HT) can potentially be used for lung cancer treatment including stereotactic radiosurgery because of its advanced image guidance and its ability to deliver highly conformal dose distributions. However, previous theoretical and simulation studies reported that the effect of respiratory motion on statically planned tomotherapy treatments may cause substantial differences between the calculated and actual delivered radiation isodose distribution, particularly when the treatment is hypofractionated. In order to determine the dosimetric effects of motion upon actual HT treatment delivery, phantom film dosimetry measurements were performed under static and moving conditions using a clinical HT treatment unit. The motion phantom system was constructed using a programmable motor, a base, a moving platform and a life size lung heterogeneity phantom with wood inserts representing lung tissue with a 3.0 cm diameter spherical tumour density equivalent insert. In order to determine the effects of different motion and tomotherapy delivery parameters, treatment plans were created using jaw sizes of 1.04 cm and 2.47 cm, with incremental gantry rotation periods between the minimum allowed (10 s) and the maximum allowed (60 s). The couch speed varied from 0.009 cm s(-1) to 0.049 cm s(-1), and delivered to a phantom under static and dynamic conditions with peak-to-peak motion amplitudes of 1.2 cm and 2 cm and periods of 3 and 5 s to simulate human respiratory motion of lung tumours. A cylindrical clinical target volume (CTV) was contoured to tightly enclose the tumour insert. 2.0 Gy was prescribed to 95% of the CTV. Two-dimensional dose was measured by a Kodak EDR2 film. Dynamic phantom doses were then quantitatively compared to static phantom doses in terms of axial dose profiles, cumulative dose volume histograms (DVH), percentage of CTV receiving the prescription dose and the minimum dose received by 95% of the CTV. The larger motion amplitude resulted in more under-dosing at the ends of the CTV in the axis of motion, and this effect was greater for the smaller jaw size plans. Due to the size of the penumbra, the 2.47 cm jaw plans provide adequate coverage for smaller amplitudes of motion, +/-0.6 cm in our experiment, without adding any additional margin in the axis of motion to the treatment volume. The periodic heterogeneous patterns described by previous studies were not observed from the single fraction of the phantom measurement. Besides the jaw sizes, CTV dose coverage is not significantly dependent on machine and phantom motion periods. The lack of adverse synchronization patterns from both results validate that HT is a safe technique for treating moving target and hypofractionation.

Evaluation of the reproducibility of lung motion probability distribution function (PDF) using dynamic MRI.

Treatment planning based on probability distribution function (PDF) of patient geometries has been shown a potential off-line strategy to incorporate organ motion, but the application of such approach highly depends upon the reproducibility of the PDF. In this paper, we investigated the dependences of the PDF reproducibility on the imaging acquisition parameters, specifically the scan time and the frame rate. Three healthy subjects underwent a continuous 5 min magnetic resonance (MR) scan in the sagittal plane with a frame rate of approximately 10 f s-1, and the experiments were repeated with an interval of 2 to 3 weeks. A total of nine pulmonary vessels from different lung regions (upper, middle and lower) were tracked and the dependences of their displacement PDF reproducibility were evaluated as a function of scan time and frame rate. As results, the PDF reproducibility error decreased with prolonged scans and appeared to approach equilibrium state in subjects 2 and 3 within the 5 min scan. The PDF accuracy increased in the power function with the increase of frame rate; however, the PDF reproducibility showed less sensitivity to frame rate presumably due to the randomness of breathing which dominates the effects. As the key component of the PDF-based treatment planning, the reproducibility of the PDF affects the dosimetric accuracy substantially. This study provides a reference for acquiring MR-based PDF of structures in the lung.

Intensity-modulated radiation therapy (IMRT) dosimetry of the head and neck: a comparison of treatment plans using linear accelerator-based IMRT and helical tomotherapy.

PURPOSE: To date, most intensity-modulated radiation therapy (IMRT) delivery has occurred using linear accelerators (linacs), although helical tomotherapy has become commercially available. To quantify the dosimetric difference, we compared linac-based and helical tomotherapy-based treatment plans for IMRT of the oropharynx. METHODS AND MATERIALS: We compared the dosimetry findings of 10 patients who had oropharyngeal carcinoma. Five patients each had cancers in the base of the tongue and tonsil. Each plan was independently optimized using either the CORVUS planning system (Nomos Corporation, Sewickly, PA), commissioned for a Varian 2300 CD linear accelerator (Varian Medical Systems, Palo Alto, CA) with 1-cm multileaf collimator leaves, or helical tomotherapy. The resulting treatment plans were evaluated by comparing the dose-volume histograms, equivalent uniform dose (EUD), dose uniformity, and normal tissue complication probabilities. RESULTS: Helical tomotherapy plans showed improvement of critical structure avoidance and target dose uniformity for all patients. The average equivalent uniform dose reduction for organs at risk (OARs) surrounding the base of tongue and the tonsil were 17.4% and 27.14% respectively. An 80% reduction in normal tissue complication probabilities for the parotid glands was observed in the tomotherapy plans relative to the linac-based plans. The standard deviation of the planning target volume dose was reduced by 71%. In our clinic, we use the combined dose-volume histograms for each class of plans as a reference goal for helical tomotherapy treatment planning optimization. CONCLUSIONS: Helical tomotherapy provides improved dose homogeneity and normal structure dose compared with linac-based IMRT in the treatment of oropharyngeal carcinoma resulting in a reduced risk for complications from focal hotspots within the planning target volume and for the adjacent parotid glands.

A computer simulated phantom study of tomotherapy dose optimization based on probability density functions (PDF) and potential errors caused by low reproducibility of PDF.

Lung tumor motion trajectories measured by four-dimensional CT or dynamic MRI can be converted to a probability density function (PDF), which describes the probability of the tumor at a certain position, for PDF based treatment planning. Using this method in simulated sequential tomotherapy, we study the dose reduction of normal tissues and more important, the effect of PDF reproducibility on the accuracy of dosimetry. For these purposes, realistic PDFs were obtained from two dynamic MRI scans of a healthy volunteer within a 2 week interval. The first PDF was accumulated from a 300 s scan and the second PDF was calculated from variable scan times from 5 s (one breathing cycle) to 300 s. Optimized beam fluences based on the second PDF were delivered to the hypothetical gross target volume (GTV) of a lung phantom that moved following the first PDF The reproducibility between two PDFs varied from low (78%) to high (94.8%) when the second scan time increased from 5 s to 300 s. When a highly reproducible PDF was used in optimization, the dose coverage of GTV was maintained; phantom lung receiving 10%-20% prescription dose was reduced by 40%-50% and the mean phantom lung dose was reduced by 9.6%. However, optimization based on PDF with low reproducibility resulted in a 50% underdosed GTV. The dosimetric error increased nearly exponentially as the PDF error increased. Therefore, although the dose of the tumor surrounding tissue can be theoretically reduced by PDF based treatment planning, the reliability and applicability of this method highly depend on if a reproducible PDF exists and is measurable. By correlating the dosimetric error and PDF error together, a useful guideline for PDF data acquisition and patient qualification for PDF based planning can be derived.

A phantom with reduced complexity for spatial 3-D ultrasound calibration.

The design of a new phantom for 3-D ultrasound calibration is presented. The phantom provides a viable alternative to existing phantoms that are significantly more complex and require high precision fabrication. The phantom, referred to as a "plane-of-wires" phantom, consists of two wires mounted at the same fixed height above the bottom of a water tank. Data collection for calibration involved rotating and translating the phantom so that the wires remained in a single plane parallel to the tank bottom. The mean reconstruction accuracy of the plane-of-wires calibration is 0.66 mm at a mean depth of 12.3 mm, with a precision of 1.23 mm at the same mean depth. The calibration was used to determine the volume of a cube with known volume with an error of 2.51%. The calibration performance achieved is comparable with that of existing approaches.

A method to compare supra-pubic ultrasound and CT images of the prostate: technique and early clinical results.

We describe a unique method that allows the comparison of spatially registered ultrasound (SRUS) images and computed tomography-derived contours (CTDCs) that were acquired with a minimal time lapse. As such, we have a tool that will provide validation of the spatial accuracy of the US system and that will allow comparison of anatomical boundaries derived via the two different imaging modalities. We describe the method by which the commercial US system is mechanically registered to a CT simulator and a unique data processing procedure. This data processing procedure circumvents the standard data acquisition and manual contouring sequence, thus reducing the time lapse from CT to US image acquisition to 10 minutes on average. Verification using a phantom demonstrated the method to be spatially accurate to within +/- 1 mm in the anterior-posterior (AP) and lateral directions and +/- 3 mm in the inferior-superior (IS) direction. Early clinical results gathered on 8 patients demonstrated alignment between the US and the CTDCs to be 0 mm in the AP and lateral directions and 2 mm in the IS direction, on average. The technique was used to compare the appearance of the prostate using US and CT imaging. The lateral dimension of the prostate indicated by the CTDCs was larger than that indicated by US imaging in all cases and on average by 0.9 cm. The height of the prostate in the AP direction was larger on average by 0.3 cm using CTDCs than US, and was larger by 5 mm or more in 3 out of 7 cases. The role of uncertainties in the determination of the CTDCs is examined as a possible cause and implications for treatment planning are described.

Segmentation of the prostate from suprapubic ultrasound images.

We present a technique for semiautomated segmentation of human prostates using suprapubic ultrasound (US) images. In this approach, a speckle reducing anisotropic diffusion (SRAD) is applied to enhance the images and the instantaneous coefficient of variation (ICOV) is utilized for edge detection. Segmentation is accomplished via a parametric active contour model in a polar coordinate system that is tailored to the application. The algorithm initially approximates the prostate boundary in two stages. First a primary contour is detected using an elliptical model, followed by a primary contour optimization using an area-weighted mean-difference binary flow geometric snake model. The algorithm was assessed by comparing the computer-derived contours with contours produced manually by three sonographers. The proposed method has application in radiation therapy planning and delivery, as well as in automated volume measurements for ultrasonic diagnosis. The average root mean square discrepancy between computed and manual outlines is less than the inter-observer variability. Furthermore, 76% of the computer-outlined contour is less than 1 sigma manual outline variance away from "true" boundary of prostate. We conclude that the methods developed herein possess acceptable agreement with manually contoured prostate boundaries and that they are potentially valuable tools for radiotherapy treatment planning and verification.

A modern Monte Carlo investigation of the TG-43 dosimetry parameters for an 125I seed already having AAPM consensus data.

PURPOSE: To investigate potential causes for differences in TG-43 brachytherapy dosimetry parameters in the existent literature for the model IAI-125A(125)I seed and to propose new standard dosimetry parameters. METHODS: The MCNP5 code was used for Monte Carlo (MC) simulations. Sensitivity of dose distributions, and subsequently TG-43 dosimetry parameters, was explored to reproduce historical methods upon which American Association of Physicists in Medicine (AAPM) consensus data are based. Twelve simulation conditions varying(125)I coating thickness, coating mass density, photon interaction cross-section library, and photon emission spectrum were examined. RESULTS: Varying(125)I coating thickness, coating mass density, photon cross-section library, and photon emission spectrum for the model IAI-125A seed changed the dose-rate constant by up to 0.9%, about 1%, about 3%, and 3%, respectively, in comparison to the proposed standard value of 0.922 cGy h(-1) U(-1). The dose-rate constant values by Solberg et al. ["Dosimetric parameters of three new solid core (125)I brachytherapy sources," J. Appl. Clin. Med. Phys. 3, 119-134 (2002)], Meigooni et al. ["Experimental and theoretical determination of dosimetric characteristics of IsoAid ADVANTAGE™ (125)I brachytherapy source," Med. Phys. 29, 2152-2158 (2002)], and Taylor and Rogers ["An EGSnrc Monte Carlo-calculated database of TG-43 parameters," Med. Phys. 35, 4228-4241 (2008)] for the model IAI-125A seed and Kennedy et al. ["Experimental and Monte Carlo determination of the TG-43 dosimetric parameters for the model 9011 THINSeed™ brachytherapy source," Med. Phys. 37, 1681-1688 (2010)] for the model 6711 seed were +4.3% (0.962 cGy h(-1) U(-1)), +6.2% (0.98 cGy h(-1) U(-1)), +0.3% (0.925 cGy h(-1) U(-1)), and -0.2% (0.921 cGy h(-1) U(-1)), respectively, in comparison to the proposed standard value. Differences in the radial dose functions between the current study and both Solberg et al. and Meigooni et al. were <10% for r ≤ 5 cm, and increased for r > 5 cm with a maximum difference of 29% at r = 9 cm. In comparison to Taylor and Rogers, these differences were lower (maximum of 2% at r = 9 cm). For the similarly designed model 6711 (125)I seed, differences did not exceed 0.5% for 0.5 ≤ r ≤ 10 cm. Radial dose function values varied by 1% as coating thickness and coating density were changed. Varying the cross-section library and source spectrum altered the radial dose function by 25% and 12%, respectively, but these differences occurred at r = 10 cm where the dose rates were very low. The 2D anisotropy function results were most similar to those of Solberg et al. and most different to those of Meigooni et al. The observed order of simulation condition variables from most to least important for influencing the 2D anisotropy function was spectrum, coating thickness, coating density, and cross-section library. CONCLUSIONS: Several MC radiation transport codes are available for calculation of the TG-43 dosimetry parameters for brachytherapy seeds. The physics models in these codes and their related cross-section libraries have been updated and improved since publication of the 2007 AAPM TG-43U1S1 report. Results using modern data indicated statistically significant differences in these dosimetry parameters in comparison to data recommended in the TG-43U1S1 report. Therefore, it seems that professional societies such as the AAPM should consider reevaluating the consensus data for this and others seeds and establishing a process of regular evaluations in which consensus data are based upon methods that remain state-of-the-art.

Independent dosimetric assessment of the model EP917 episcleral brachytherapy plaque.

PURPOSE: To investigate the influence of slot design on dose distributions and dose-volume histograms (DVHs) for the model EP917 plaque for episcleral brachytherapy. METHODS: Dimensions and orientations of the slots were measured for three model EP917 plaques and compared to data in the Plaque Simulator (PS) treatment planning software (version 5.7.6). These independently determined coordinates were incorporated into the MCNP Monte Carlo simulation environment to obtain dose from the plaques in a water environment and in a clinical environment with ocular structures. A tumor volume was simulated as 5 mm in apical height and 11 mm in basal diameter. Variations in plaque mass density and composition; slot length, width, and depth; seed positioning; and Ag-marker rod positioning were simulated to examine their influence on plaque central axis (CAX) and planar dose distributions, and DVHs. RESULTS: Seed shifts in a single slot toward the eye and shifts of the(125)I-coated Ag rod within the capsule had the greatest impact on CAX dose distribution. A shift of 0.0994 mm toward the eye increased dose by 14%, 9%, 4.3%, and 2.7% at 1, 2, 5, and 10 mm, respectively, from the inner sclera. When examining the fully-modeled plaque in the ocular geometry, the largest dose variations were caused by shifting the Ag rods toward the sclera and shifting the seeds away from the globe when the slots were made 0.51 mm deeper, causing +34.3% and -69.4% dose changes to the outer sclera, respectively. At points along the CAX, dose from the full plaque geometry using the measured slot design was 2.4%±1.1% higher than the manufacturer-provided slot design and 2.2%±2.3% higher than the homogeneous calculation of PS treatment planning results. The ratio of D10 values for the measured slot design to the D10 values for the manufacturer-provided slot design was higher by 9%, 10%, and 19% for the tumor, inner sclera, and outer sclera, respectively. In comparison to the measured slot design, a theoretical plaque having narrower and deeper slots delivered 30%, 37%, and 62% lower D10 doses to the tumor, inner sclera, and outer sclera, respectively. CONCLUSIONS: While the measured positions of the slots on the model EP917 plaque were in close agreement (<0.7 mm) with the PS values, small differences in the slot shape caused substantial differences in dose distributions and DVH metrics. Increasing slot depth by 0.1 mm decreased outer scleral dose by 20%, yet shifting the Ag rods in the seeds toward the globe by 0.1 mm increased outer scleral dose by 35%. The clinical medical physicist is advised to measure these types of plaques upon acceptance testing before clinical use to inspect slot shape and position for comparison with data used for treatment planning purposes.

Determination of prescription dose for Cs-131 permanent implants using the BED formalism including resensitization correction.

PURPOSE: The current widely used biological equivalent dose (BED) formalism for permanent implants is based on the linear-quadratic model that includes cell repair and repopulation but not resensitization (redistribution and reoxygenation). The authors propose a BED formalism that includes all the four biological effects (4Rs), and the authors propose how it can be used to calculate appropriate prescription doses for permanent implants with Cs-131. METHODS: A resensitization correction was added to the BED calculation for permanent implants to account for 4Rs. Using the same BED, the prescription doses with Au-198, I-125, and Pd-103 were converted to the isoeffective Cs-131 prescription doses. The conversion factor F, ratio of the Cs-131 dose to the equivalent dose with the other reference isotope (Fr: with resensitization, Fn: without resensitization), was thus derived and used for actual prescription. Different values of biological parameters such as α, ß, and relative biological effectiveness for different types of tumors were used for the calculation. RESULTS: Prescription doses with I-125, Pd-103, and Au-198 ranging from 10 to 160 Gy were converted into prescription doses with Cs-131. The difference in dose conversion factors with (Fr) and without (Fn) resensitization was significant but varied with different isotopes and different types of tumors. The conversion factors also varied with different doses. For I-125, the average values of Fr/Fn were 0.51/0.46, for fast growing tumors, and 0.88/0.77 for slow growing tumors. For Pd-103, the average values of Fr/Fn were 1.25/1.15 for fast growing tumors, and 1.28/1.22 for slow growing tumors. For Au-198, the average values of Fr/Fn were 1.08/1.25 for fast growing tumors, and 1.00/1.06 for slow growing tumors. Using the biological parameters for the HeLa/C4-I cells, the averaged value of Fr was 1.07/1.11 (rounded to 1.1), and the averaged value of Fn was 1.75/1.18. Fr of 1.1 has been applied to gynecological cancer implants with expected acute reactions and outcomes as expected based on extensive experience with permanent implants. The calculation also gave the average Cs-131 dose of 126 Gy converted from the I-125 dose of 144 Gy for prostate implants. CONCLUSIONS: Inclusion of an allowance for resensitization led to significant dose corrections for Cs-131 permanent implants, and should be applied to prescription dose calculation. The adjustment of the Cs-131 prescription doses with resensitization correction for gynecological permanent implants was consistent with clinical experience and observations. However, the Cs-131 prescription doses converted from other implant doses can be further adjusted based on new experimental results, clinical observations, and clinical outcomes.