3D Dosimetry for Electron Flash Radiotherapy: Assessment of Radiochromic Dosimeter Phantoms with Optical CT Scanning as a 3D Dosimetry System

PURPOSE/OBJECTIVE(S)

Many of the dosimeters used in conventional radiation therapy exhibit dose rate dependence which prohibits their use in ultra-high-dose-rate (FLASH) radiation therapy. Radiochromic plastic dosimeter PRESAGE® has been used for 3D dosimetry for many years. We hypothesized that these phantoms would show dose-rate independence throughout both the conventional and FLASH RT regimes, indicating these phantoms exhibit qualities useful for relative 3D dosimetry in FLASH electron beams.

MATERIALS/METHODS

FLASH experiments were performed using a commercially available linear accelerator, converted to deliver an ultra-high-dose-rate 10 MeV electron beam. The LINAC delivered approximately 0.7 Gy/pulse for FLASH irradiations. Dose rate was varied from about 40 Gy/s to 240 Gy/s by changing the repetition rate. PRESAGE phantoms were irradiated en face at six FLASH dose rates: 40 Gy/s, 80 Gy/s, 120 Gy/s, 160 Gy/s, 200 Gy/s, and 240 Gy/s. EBT film and scintillator measurements were used to verify dose delivered. The optical response of the PESAGE phantom versus delivered dose was evaluated with various known doses. A novel parallel-beam optical CT scanner, utilizing fiber optic taper for collimated images, was developed for fast, high resolution, and accurate readout of 3D dosimeters. Percent depth dose curves for various FLASH dose rates and conventional dose rate beams were generated and compared based on the optical response versus dose measurements. Percent depth dose curves from Monte Carlo calculation of the presage phantom were also compared.

RESULTS

As shown in Table 1, the percent depth dose as a function of depth for six FLASH dose rates (240-40 Gy/s) are nearly identical, indicating that optical response of PRESAGE is dose-rate independent. The optical density of PRESAGE phantom was confirmed to be linear with absorbed dose for all FLASH dose rates, consistent with the observation at regular treatment dose rates.

CONCLUSION

PRESAGE phantoms show dose-rate independence in electron beams for a wide range of dose rates from conventional to ultra-high-dose-rates, indicating these phantoms can be useful for relative 3D dose measurements in FLASH electron beams. Future experiments will be undertaken as part of the commissioning of a commercially available FLASH radiotherapy unit.

Validation of a High-Throughput Dicentric Chromosome Assay Using Complex Radiation Exposures

Validation of biodosimetry assays is routinely performed using primarily orthovoltage irradiators at a conventional dose rate of approximately 1 Gy/min. However, incidental/ accidental exposures caused by nuclear weapons can be more complex. The aim of this work was to simulate the DNA damage effects mimicking those caused by the detonation of a several kilotons improvised nuclear device (IND). For this, we modeled complex exposures to: 1. a mixed (photons + IND-neutrons) field and 2. different dose rates that may come from the blast, nuclear fallout, or ground deposition of radionuclides (ground shine). Additionally, we assessed whether myeloid cytokines affect the precision of radiation dose estimation by modulating the frequency of dicentric chromosomes. To mimic different exposure scenarios, several irradiation systems were used. In a mixed field study, human blood samples were exposed to a photon field enriched with neutrons (ranging from 10% to 37%) from a source that mimics Hiroshima's A-bomb's energy spectrum (0.2-9 MeV). Using statistical analysis, we assessed whether photons and neutrons act in an additive or synergistic way to form dicentrics. For the dose rates study, human blood was exposed to photons or electrons at dose rates ranging from low (where the dose was spread over 32 h) to extremely high (where the dose was delivered in a fraction of a microsecond). Potential effects of cytokine treatment on biodosimetry dose predictions were analyzed in irradiated blood subjected to Neupogen or Neulasta for 24 or 48 h at the concentration recommended to forestall manifestation of an acute radiation syndrome in bomb survivors. All measurements were performed using a robotic station, the Rapid Automated Biodosimetry Tool II, programmed to culture lymphocytes and score dicentrics in multiwell plates (the RABiT-II DCA). In agreement with classical concepts of radiation biology, the RABiT-II DCA calibration curves suggested that the frequency of dicentrics depends on the type of radiation and is modulated by changes in the dose rate. The resulting dose-response curves suggested an intermediate dicentric yields and additive effects of photons and IND-neutrons in the mixed field. At ultra-high dose rate (600 Gy/s), affected lymphocytes exhibited significantly fewer dicentrics (P < 0.004, t test). In contrast, we did not find the dose-response modification effects of radiomitigators on the yields of dicentrics (Bonferroni corrected P > 0.006, ANOVA test). This result suggests no bias in the dose predictions should be expected after emergency cytokine treatment initiated up to 48 h prior to blood collection for dicentric analysis.

Fiducial detection and registration for 3D IMRT QA with organ-specific dose information

PURPOSE

Two-dimensional (2D) IMRT QA has been widely performed in Radiation Oncology clinic. However, concerns regarding its sensitivity in detecting delivery errors and its clinical meaning have been raised in publications. In this study, a robust methodology of three-dimensional (3D) IMRT QA using fiducial registration and structure-mapping was proposed to acquire organ-specific dose information.

METHODS

Computed tomography (CT) markers were placed on the PRESAGE dosimeter as fiducials before CT simulation. Subsequently, the images were transferred to the treatment planning system to create a verification plan for the examined treatment plan. Patient's CT images were registered to the CT images of the dosimeter for structure mapping according to the positions of the fiducials. After irradiation, the 3D dose distribution was read-out by an optical-CT (OCT) scanner with fiducials shown on the OCT dose images. An automatic localization algorithm was developed in MATLAB to register the markers in the OCT images to those in the CT images of the dosimeter. SlicerRT was used to show and analyze the results. Fiducial registration error was acquired by measuring the discrepancies in 20 fiducial registrations, and thus the fiducial localization error and target registration error (TRE) was estimated.

RESULTS

Dosimetry comparison between the calculated and measured dose distribution in various forms were presented, including 2D isodose lines comparison, 3D isodose surfaces with patient's anatomical structures, 2D and 3D gamma index, dose volume histogram and 3D view of gamma failing points. From the analysis of 20 fiducial registrations, fiducial registration error was measured to be 0.62 mm and fiducial localization error was calculated to be 0.44 mm. Target registration uncertainty of the proposed methodology was estimated to be within 0.3 mm in the area of dose measurement.

CONCLUSIONS

This study proposed a robust methodology of 3D measurement-based IMRT QA for organ-specific dose comparison and demonstrated its clinical feasibility.

Temperature dependence and temporal stability of stacked radiochromic sheets for three-dimensional dose verification

PURPOSE

Recently a novel radiochromic sheet dosimeter, termed as PRESAGE sheets, consisting of leuco crystal violet dye and radical initiator had been developed and characterized. This study examines the dosimeter's temporal stability and storage temperature dependence postirradiation, and its applicability for dose verification in three dimensions (3D) as a stack dosimeter.

METHODS

PRESAGE sheets were irradiated using 6 MV photons at a dose range of 0-20 Gy with the change in optical density measured using a flatbed scanner. Following their irradiation, PRESAGE sheets were stored in different temperature environments (-18 °C, 4 °C, and 22 °C) and scanned at different time points, ranging from 1 to 168 h postirradiation, to track changes in measured signal and linearity of dose response. Multiple PRESAGE sheets were bound together to create a 12 × 13 × 8.7 cm³ film stack, with EBT3 film inserted between the sheets in the central region of the stack, that was treated using a clinical VMAT plan. Based on the results from the time and storage temperature study, two-dimensional (2D) relative dose distribution measurements in PRESAGE were acquired promptly following irradiation at selected planes in the coronal, sagittal, and axial orientation of the film stack and compared to the treatment planning system calculations in their respective axes. Dose distribution measurements on the coronal axis of the stack dosimeter were also independently verified using EBT3 film.

RESULTS

The dose response was observed to be linear (R² > 0.995) with sheets stored in colder temperatures retaining their signal and dose response sensitivity for extended periods postirradiation. Sheets stored in 22 °C environment should be measured within an hour postirradiation. Sheets stored in a 4 °C and -18 °C environment can be scanned up to 20- and 72 h postirradiation, respectively, while preserving the integrity of their dose response sensitivity and linearity of dose response within a mean absolute percent error of 2.0%. For instance, at 20 h postirradiation the dose response sensitivity for sheets stored in a -18 °C, 4 °C, and 22 °C temperature environment was measured to be 97%, 91%, and 77% of their original values measured within an hour postirradiation, respectively. The 2D gamma pass rate for central slices exceed 95% for PRESAGE film stack compared with treatment planning system on selected planes in the axial, coronal, and sagittal orientation and EBT3 film in the coronal orientation using a 2D gamma index of 2%/2mm. The gamma pass rate in comparing the calculated dose distribution with the measured dose distribution from PRESAGE-LCV was observed to decrease in sheets scanned at later elapsed times postirradiation. In one example, the gamma pass rate for 2%/2mm criteria in the coronal plane was observed to decrease from 97.7% pass rate when scanned within an hour postirradiation to 92.1% pass rate when scanned at 20 h postirradiation under room temperature conditions.

CONCLUSIONS

This is the first study to demonstrate that the temporal stability of PRESAGE sheets can be enhanced through its storage in colder temperature environments postirradiation and that sheets as a film stack dosimeter hold promise for precise relative dose distribution measurements in 3D where advanced optical CT is unavailable.

Focused ultrasound induced-blood-brain barrier opening in mouse brain receiving radiosurgery dose of radiation enhances local delivery of systemic therapy

OBJECTIVE

Investigate the temporal effects of focused ultrasound (FUS)-induced blood-brain barrier (BBB) opening in post-radiotherapy mouse brains.

METHODS AND MATERIALS

C57B6 mice without tumors were used to simulate the scenario after gross total resection (GTR) of brain tumor. Radiation dose of 6 Gy x 5 was delivered to one-hemisphere of the mouse brain. FUS-induced BBB-opening was delivered to the irradiated and non-irradiated brain and was confirmed with MRI. Dynamic MRI was performed to evaluate blood vessel permeability. Two time points were selected: acute (2 days after radiation) and chronic (31 days after radiation).

RESULTS

BBB opening was achieved after FUS in the irradiated field as compared to the contralateral non-irradiated brain without any decrease in permeability. In the acute group, a trend for higher gadolinium concentration was observed in radiated field.

CONCLUSION

Localized BBB-opening can be successfully achieved without loss of efficacy by FUS as early as 2 days after radiotherapy.

ADVANCES IN KNOWLEDGE

Adjuvant radiation after GTR is commonly used for brain tumors. Focused ultrasound facilitated BBB-opening can be achieved without loss of efficacy in the post-irradiated brain as early as 2 days after radiation therapy. This allows for further studies on early application of FUS-mediated BBB-opening.

Dosimetric characterization of a body-conforming radiochromic sheet

Purpose

A novel radiochromic PRESAGE sheet (Heuris Inc.) with 3 mm thickness has been developed as a measurement tool for 2D dosimetry. Its inherent ability to conform to irregular surfaces makes this dosimeter advantageous for patient surface dosimetry. This study is a comprehensive investigation into the PRESAGE sheet’s dosimetric characteristic, accuracy and its potential use as a dosimeter for clinical applications.

Methods

The characterization of the dosimeter included evaluation of the temporal stability of the dose linearity, reproducibility, measurement uncertainties, dose rate, energy, temperature and angular dependence, lateral response artifacts, percent depth dose curve, and 2D dose measurement. Dose distribution measurements were acquired for regular square fields on a flat and irregular surface and an irregular modulated field on the smooth surface. All measurements were performed using an Epson 11000XL high‐resolution scanner.

Results

The examined dosimeters exhibit stable linear response, standard error of repeated measurements within 2%, negligible dose rate, energy, and angular dependence. The same linear dose response was measured while the dosimeter was in contact with a heated water surface. Gamma test and histogram analysis of the dose difference between PRESAGE and EBT3 film, PRESAGE and the treatment planning system (TPS) were used to evaluate the measured dose distributions. The PRESAGE sheet dose distributions showed good agreement with EBT3 film and TPS. A discrepancy smaller than the statistical error of the two dosimeters was reported.

Conclusions

This study established a full dosimetric characterization of the PRESAGE sheets with the purpose of laying the foundation for future clinical uses. The results presented here for the comparison of this novel dosimeter with those currently in use reinforce the possibility of using this dosimeter as an alternative for irregular surface dose measurements.

3D isocentricity analysis for clinical linear accelerators

PURPOSE

To perform a three-dimensional (3D) concurrent isocentricity measurement of a clinical linear accelerator's (linac) using a single 3D dosimeter, PRESAGE.

METHODS

A 3D dosimeter, PRESAGE, set up on the treatment couch of a Varian TrueBeam LINAC using the setup lasers, was irradiated under gantry angles of 0 ∘ , 50 ∘ , 160 ∘ , and 270 ∘ with the couch fixed at 0 ∘ and subsequently, under couch angles of 10 ∘ , 330 ∘ , 300 ∘ , and 265 ∘ with the gantry fixed at 270 ∘ . The 1 cm 2 (at 100 cm SAD) square fields were delivered at 6 MV with 800 MU/field. After irradiation, the dosimeter was scanned using a single-beam optical scanner and images were reconstructed with submillimeter resolution using filtered back-projection. Postprocessing was used to extract views parallel to the star-shot planes from which beam trajectories and the smallest circles enclosing these were drawn and extracted. These circles and information from the view orthogonal to both star-shots were used to represent the rotational centers as spheroids. The linac isocenter was defined by the distribution of midpoints between any, randomly selected, points lying inside the center spheroids defined by the table and gantry rotations; isocenter location and size were defined by the average midpoint and the distribution's semi-axes. Collimator rotations were not included in this study.

RESULTS

Relative to the setup center defined by lasers, the table and gantry rotation center coordinates (lat., long., vert.) were measured in units of millimeters, to be (-0.24, 0.18, -0.52) and (0.10, 0.53, -0.52), respectively. Displacements from the setup center were 0.60 and 0.75 mm for the table and gantry centers, while the distance between them measured 0.49 mm. The linac's radiation isocenter was calculated to be at (-0.07, -0.17, 0.51) relative to the setup lasers and its size was found to be most easily described by a spheroid prolate in vertical direction with semi-axis lengths of 0.13 and 0.23 mm for the lateral-longitudinal and vertical directions, respectively.

CONCLUSIONS

This study demonstrates how to measure the location and sizes of rotational centers in 3D with one setup. The proposed method provides a more comprehensive view on the isocentricity of LINAC than the conventional two-dimensional film measurements. Additionally, a new definition of isocenter and its size was proposed.

Performance of the cone beam computed tomography-based patient positioning system on the Gamma Knife Icon™

PURPOSE

Cone beam computed tomography (CBCT) imaging has been implemented on the Leksell Gamma Knife^® Icon™ for assessing patient positioning in mask-based Gamma Knife radiosurgery. The purpose of this study was to evaluate the performance of the CBCT-based patient positioning system as a tool for frameless Gamma Knife radiosurgery.

METHODS

Daily quality assurance (QA) CBCT precision test results from a 12-month period were analyzed for the geometric accuracy and the stability of the imager. The performance of the image acquisition module and the image registration algorithm was evaluated using an anthropomorphic head phantom (CIRS Inc., Norfolk, VA) and a XYZR axis manual positioning stage (TOAUTO Inc., Guangdong, China). The head phantom was fixed on a mask adaptor and manually translated in the X, Y, Z directions or rotated around the X, Y, Z axes in the range of ±10 mm or ±10º. A CBCT scan was performed after each manual position setup followed by an image registration to the reference scan. To assess the overall setup uncertainties in fractionated treatment, two cylindrical Presage phantoms (Heuris Inc., Skillman, NJ) of 15 cm diameter and 10 cm height were irradiated with identical prescription dose and shot placement following standard mask-based treatment workflow according to two different fraction schedules: a single fraction treatment of 7.5 Gy and a 5-fraction treatment with 1.5 Gy per fraction.

RESULTS

The averaged vector deviations of the four marks from their preset values are 0.087, 0.085, 0.095, and 0.079 mm from the 212 daily QA tests. The averaged displacements in the X, Y, Z coordinates and the pitch, yaw, roll angles from the image registration tests are 0.23, 0.27, 0.14, 0.32º, 0.19º, 0.31º from the manual setup. The corresponding maximum differences are 0.41, 0.33, 0.29 mm, 0.45º, 0.31º, and 0.43º, respectively. Compared to the treatment plan using the 2% & 1 mm criteria, the averaged 2D Gamma passing rate is 98.25% for the measured dose distribution from the Presage phantom with 1-fraction irradiation and 95.12% for the 5-fraction irradiation. The averaged Gamma passing rates are 99.53% and 98.16% for the 1-fraction and 5-fraction irradiations using the 2% & 2 mm criteria.

CONCLUSIONS

The CBCT imager and the image registration algorithm can reproduce phantom position with <0.5 mm/0.5º uncertainty. A systematic contribution from the interfraction phantom repositioning procedure was observed in the Gamma analysis over the irradiated volumes of two end-to-end test phantoms.

An investigation of clinical treatment field delivery verification using cherenkov imaging: IMRT positioning shifts and field matching

PURPOSE

Cherenkov light emission has been shown to correlate with ionizing radiation dose delivery in solid tissue. An important clinical application of Cherenkov light is the real-time verification of radiation treatment delivery in vivo. To test the feasibility of treatment field verification, Cherenkov light images were acquired concurrent with radiation beam delivery to standard and anthropomorphic phantoms. Specifically, we tested two clinical treatment scenarios: (a) Observation of field overlaps or gaps in matched 3D fields and (b) Patient positioning shifts during intensity modulated radiation therapy (IMRT) field delivery. Further development of this technique would allow real-time detection of treatment delivery errors on the order of millimeters so that patient safety and treatment quality can be improved.

METHODS

Cherenkov light emission was captured using a PI-MAX4 intensified charge coupled device (ICCD) system (Princeton Instruments). All radiation delivery was performed using a Varian Trilogy linear accelerator (linac) operated at 6 MV or 18 MV for photon and 6 MeV or 16 MeV for electron studies. Field matching studies were conducted with photon and electron beams at gantry angles of 0°, 15°, and 45°. For each modality and gantry angle, a total of three data sets were acquired. Overlap and gap distances of 0, 2, 5, and 10 mm were tested and delivered to solid phantom material of 30 × 30 × 5 cm³ . Phantom materials used were white plastic water and brown solid water. Tests were additionally performed on an anthropomorphic phantom with an irregular surface. Positioning shift studies were performed using IMRT fields delivered to a thoracic anthropomorphic phantom. For thoracic phantom measurements, the camera was placed laterally to observe the entire right side of the phantom. Fields were delivered with known translational patient positioning shifts in four directions. Changes in the Cherenkov fluence were evaluated through the generation of difference maps from unshifted Cherenkov images. All images were evaluated using ImageJ, Python, and MATLAB software packages.

RESULTS

For matched fields, Cherenkov images were able to quantitate matched field separations with discrepancies between 2 and 4 mm, depending on gantry angle and beam energy or modality. For all photon and electron beams delivered at a gantry angle of 0°, image analysis indicated average discrepancies of less than 2 mm for all field gaps and overlaps, with 83% of matched fields exhibiting discrepancies less than 1 mm. Beams delivered obliquely to the phantom surface exhibited average discrepancies as high as 4 mm for electron beams delivered at large oblique angles. Finally, for IMRT field delivery, vertical and lateral patient positioning shifts of 2 mm were detected in some cases, indicating the potential detectability threshold of using this technique alone.

CONCLUSIONS

Our study indicates that Cherenkov imaging can be used to support and bolster current treatment delivery verification techniques, improving our ability to recognize and rectify millimeter-scale delivery and positioning errors.

Velocity-based Adaptive Registration and Fusion for Fractionated Stereotactic Radiosurgery Using the Small Animal Radiation Research Platform

PURPOSE

To implement Velocity-based image fusion and adaptive deformable registration to enable treatment planning for preclinical murine models of fractionated stereotactic radiosurgery (fSRS) using the small animal radiation research platform (SARRP).

METHODS AND MATERIALS

C57BL6 mice underwent 3 unique cone beam computed tomography (CBCT) scans: 2 in the prone position and a third supine. A single T1-weighted post-contrast magnetic resonance imaging (MRI) series of a murine metastatic brain tumor model was selected for MRI-to-CBCT registration and gross tumor volume (GTV) identification. Two arms were compared: Arm 1, where we performed 3 individual MRI-to-CBCT fusions using rigid registration, contouring GTVs on each, and Arm 2, where the authors performed MRI-to-CBCT fusion and contoured GTV on the first CBCT followed by Velocity-based adaptive registration. The first CBCT and associated GTV were exported from MuriPlan (Xstrahl Life Sciences) into Velocity (Varian Medical Systems, Inc, Palo Alto, CA). In Arm 1, the second and third CBCTs were exported similarly along with associated GTVs (Arm 1), while in Arm 2, the first (prone) CBCT was fused separately to the second (prone) and third (supine) CBCTs, performing deformable registrations on initial CBCTs and applying resulting matrices to the contoured GTV. Resulting GTVs were compared between Arms 1 and 2.

RESULTS

Comparing GTV overlays using repeated MRI fusion and GTV delineation (Arm 1) versus those of Velocity-based CBCT and GTV adaptive fusion (Arm 2), mean deviations ± standard deviation in the axial, sagittal, and coronal planes were 0.46 ± 0.16, 0.46 ± 0.22, and 0.37 ± 0.22 mm for prone-to-prone and 0.52 ± 0.27, 0.52 ± 0.36, and 0.68 ± 0.31 mm for prone-to-supine adaptive fusions, respectively.

CONCLUSIONS

Velocity-based adaptive fusion of CBCTs and contoured volumes allows for efficient fSRS planning using a single MRI-to-CBCT fusion. This technique is immediately implementable on current SARRP systems, facilitating advanced preclinical treatment paradigms using existing clinical treatment planning software.

Microdosimetric characteristics of 50 kV X rays at different depths for breast intraoperative radiotherapy

An intraoperative radiation therapy (IORT) device with 50 kV X rays was designed to deliver a single dose to the tumour bed after local excision of breast cancer. The quality of a radiation can be determined by the microscopic distribution of energy transfers along and across the charged particle tracks. The lineal energy, y, serves as an accurate measure of local energy concentration. The dose mean lineal energy, yD, is an indicator of radiation quality. For low linear energy transfer radiation, the ratio of its dose mean lineal energy to that of (60)Co gamma rays can serve as a good indicator of the relative biological effectiveness (RBE) at low doses. In this study, microdosimetric simulations are performed for soft tissue irradiated by 50 kV X rays generated from the IORT device, with a 4-cm breast applicator attached. All energy transfers are recorded with the location coordinates in the tissue. Microdosimetric single events in a sphere of 1 µm in diameter are scored as a function of radial distances from the applicator surface. Single-event spectra are then constructed. From those single-event spectra, dose mean lineal energy is calculated. Compared with dose mean lineal energy of (60)Co gamma rays, the estimated RBEs at low doses are given for the X rays at different depths in the tissue. The RBEs at clinically relevant doses, as a function of depth, are also presented.

Fetal radiation monitoring and dose minimization during intensity modulated radiation therapy for glioblastoma in pregnancy

We examined the fetal dose from irradiation of glioblastoma during pregnancy using intensity modulated radiation therapy (IMRT), and describe fetal dose minimization using mobile shielding devices. A case report is described of a pregnant woman with glioblastoma who was treated during the third trimester of gestation with 60 Gy of radiation delivered via a 6 MV photon IMRT plan. Fetal dose without shielding was estimated using an anthropomorphic phantom with ion chamber and diode measurements. Clinical fetal dose with shielding was determined with optically stimulated luminescent dosimeters and ion chamber. Clinical target volume (CTV) and planning target volume (PTV) coverage was 100 and 98 % receiving 95 % of the prescription dose, respectively. Normal tissue tolerances were kept below quantitative analysis of normal tissue effects in the clinic (QUANTEC) recommendations. Without shielding, anthropomorphic phantom measurements showed a cumulative fetal dose of 0.024 Gy. In vivo measurements with shielding in place demonstrated a cumulative fetal dose of 0.016 Gy. The fetal dose estimated without shielding was 0.04 % and with shielding was 0.026 % of the target dose. In vivo estimation of dose equivalent received by the fetus was 24.21 mSv. Using modern techniques, brain irradiation can be delivered to pregnant patients in the third trimester with very low measured doses to the fetus, without compromising target coverage or normal tissue dose constraints. Fetal dose can further be reduced with the use of shielding devices, in keeping with the principle of as low as reasonably achievable.

Performance of an improved first generation optical CT scanner for 3D dosimetry

Performance analysis of a modified 3D dosimetry optical scanner based on the first generation optical CT scanner OCTOPUS is presented. The system consists of PRESAGE dosimeters, the modified 3D scanner, and a new developed in-house user control panel written in Labview program which provides more flexibility to optimize mechanical control and data acquisition technique. The total scanning time has been significantly reduced from initial 8 h to ∼2 h by using the modified scanner. The functional performance of the modified scanner has been evaluated in terms of the mechanical integrity uncertainty of the data acquisition process. Optical density distribution comparison between the modified scanner, OCTOPUS and the treatment plan system has been studied. It has been demonstrated that the agreement between the modified scanner and treatment plans is comparable with that between the OCTOPUS and treatment plans.

Recommended ethics curriculum for medical physics graduate and residency programs: report of Task Group 159

The AAPM Professional Council approved the formation of a task group in 2007, whose purpose is to develop recommendations for an ethics curriculum for medical physics graduate and residency programs. Existing program's ethics curricula range in scope and content considerably. It is desirable to have a more uniform baseline curriculum for all programs. Recommended subjects areas, suggested ethics references, and a sample curriculum are included. This report recommends a reasonable ethics course time to be 15-30 h while allowing each program the flexibility to design their course.

Sensitivity calibration procedures in optical-CT scanning of BANG 3 polymer gel dosimeters

The dose response of the BANG 3 polymer gel dosimeter (MGS Research Inc., Madison, CT) was studied using the OCTOPUS laser CT scanner (MGS Research Inc., Madison, CT). Six 17 cm diameter and 12 cm high Barex cylinders, and 18 small glass vials were used to house the gel. The gel phantoms were irradiated with 6 and 10 MV photons, as well as 12 and 16 MeV electrons using a Varian Clinac 2100EX. Three calibration methods were used to obtain the dose response curves: (a) Optical density measurements on the 18 glass vials irradiated with graded doses from 0 to 4 Gy using 6 or 10 MV large field irradiations; (b) optical-CT scanning of Barex cylinders irradiated with graded doses (0.5, 1, 1.5, and 2 Gy) from four adjacent 4 x 4 cm2 photon fields or 6 x 6 cm2 electron fields; and (c) percent depth dose (PDD) comparison of optical-CT scans with ion chamber measurements for 6 x 6 cm2, 12 and 16 MeV electron fields. The dose response of the BANG3 gel was found to be linear and energy independent within the uncertainties of the experimental methods (about 3%). The slopes of the linearly fitted dose response curves (dose sensitivities) from the four field irradiations (0.0752 +/- 3%, 0.0756 +/- 3%, 0.0767 +/- 3%, and 0.0759 +/- 3% cm(-1) Gy(-1)) and the PDD matching methods (0.0768 +/- 3% and 0.0761 +/- 3% cm(-1) Gy(-1)) agree within 2.2%, indicating a good reproducibility of the gel dose response within phantoms of the same geometry. The dose sensitivities from the glass vial approach are different from those of the cylindrical Barex phantoms by more than 30%, owing probably to the difference in temperature inside the two types of phantoms during gel formation and irradiation, and possible oxygen contamination of the glass vial walls. The dose response curve obtained from the PDD matching approach with 16 MeV electron field was used to calibrate the gel phantom irradiated with the 12 MeV, 6 x 6 cm2 electron field. Three-dimensional dose distributions from the gel measurement and the Eclipse planning system (Varian Corporation, Palo Alto, CA) were compared and evaluated using 3% dose difference and 2 mm distance-to-agreement criteria.

Improved Biochemical Control and Clinical Disease-Free Survival with Intraoperative Versus Preoperative Preplanning for Transperineal Interstitial Permanent Prostate Brachytherapy

PURPOSE

We hypothesized that intraoperative preplanning for transperineal interstitial permanent prostate brachytherapy may yield better prostate cancer control than preoperative preplanning. We tested this hypothesis by comparing treatment outcomes of patients who underwent implantation using these two preplanning methods.

PATIENTS AND METHODS

We analyzed the data of 135 consecutive patients with localized prostate cancer treated from 1996 to 2001 with transperineal interstitial permanent prostate brachytherapy+/-preimplantation hormonal therapy: 42 received preoperative preplanning (group 1), and 93 underwent intraoperative preplanning (group 2). Biochemical status was assessed using two failure definitions: American Society for Therapeutic Radiology and Oncology (ASTRO) (three consecutive rises in prostate-specific antigen level) and Houston (prostate-specific antigen level>or=current nadir+2 ng/mL). Clinical disease-free survival and postimplantation dosimetry were also examined.

RESULTS

All disease control outcomes were superior for group 2. The 4-year ASTRO biochemical no evidence of disease rate was 80% for group 1 versus 94% for group 2. The 4-year Houston biochemical no evidence of disease rate was 82% for group 1 versus 96% for group 2. The 4-year clinical disease-free survival rate was 87% for group 1 versus 99% for group 2. Preplanning method (preoperative versus intraoperative) remained predictive of disease control outcomes in multivariate analyses with the covariates of pretreatment prostate-specific antigen level, Gleason score, clinical stage, and case sequence number (proxy for brachytherapist experience and "stage migration"). Dosimetric prostate coverage was superior for group 2. The mean percentage of the prescription dose delivered to 90% of the prostate volume (%D90) was 75% for group 1 versus 90% for group 2. A %D90>or=70% predicted for improved disease control; fewer group 1 than 2 patients met this dosimetric criterion (55% versus 87%).

DISCUSSION

Intraoperative preplanning yielded superior disease control outcomes in this analysis, likely due at least in part to improved dosimetric prostate coverage with this method. Although not mandatory for obtaining high prostate brachytherapy efficacy, intraoperative preplanning nevertheless may offer an excellent means of improving dosimetric prostate coverage and therefore disease control outcomes.

Three-dimensional dose verification for intensity modulated radiation therapy using optical CT based polymer gel dosimetry

Dose distributions generated from intensity-modulated-radiation-therapy (IMRT) treatment planning present high dose gradient regions in the boundaries between the target and the surrounding critical organs. Dose accuracy in these areas can be critical, and may affect the treatment. With the increasing use of IMRT in radiotherapy, there is an increased need for a dosimeter that allows for accurate determination of three-dimensional (3D) dose distributions with high spatial resolution. In this study, polymer gel dosimetry and an optical CT scanner have been employed to implement 3D dose verification for IMRT. A plastic cylinder of 17 cm diameter and 12 cm height, filled with BANG3 polymer gels (MGS Research, Inc., Madison, CT) and modified to optimal dose-response characteristics, was used for IMRT dose verification. The cylindrical gel phantom was immersed in a 24 x 24 x 20 cm water tank for an IMRT irradiation. The irradiated gel sample was then scanned with an optical CT scanner (MGS Research Inc., Madison, CT) utilizing a single He-Ne laser beam and a single photodiode detector. Similar to the x-ray CT process, filtered back-projection was used to reconstruct the 3D dose distribution. The dose distributions measured from the gel were compared with those from the IMRT treatment planning system. For comparative dosimetry, a solid water phantom of 24 x 24 x 20 cm, having the same geometry as the water tank for the gel phantom, was used for EDR2 film and ion chamber measurements. Root mean square (rms) deviations for both dose difference and distance-to-agreement (DTA) were used in three-dimensional analysis of the dose distribution comparison between treatment planning calculations and the gel measurement. Comparison of planar dose distributions among gel dosimeter, film, and the treatment planning system showed that the isodose lines were in good agreement on selected planes in axial, coronal, and sagittal orientations. Absolute point-dose verification was performed with ion chamber measurements at four different points, varying from 48% to 110% of the prescribed dose. The measured and calculated doses were found to agree to within 4.2% at all measurement points. For the comparison between the gel measurement and treatment planning calculations, rms deviations were 2%-6% for dose difference and 1-3 mm for DTA, at 60%-110% doses levels. The results from this study show that optical CT based polymer gel dosimetry has the potential to provide a high resolution, accurate, three-dimensional tool for IMRT dose distribution verification.

Calculated microdosimetric characteristics of 125I and 103Pd brachytherapy seeds at different depths in water

Both (125)I and (103)Pd sources have been widely used in the permanent prostate implant. An important consideration for the choice of brachytherapy sources is the relative biological effectiveness (RBE) for the source/seed used in the implantation. As RBE is closely related to the microdosimetric parameter, it is desirable to calculate the dose mean lineal energies for both (125)I and (103)Pd at various radial distances to the seed surface. Monte Carlo simulation was performed for photons emitted from (125)I and (103)Pd. Energy depositions from photons and all their secondary electrons were tracked. Dose distributions of lineal energy, d(y), were calculated for spheres of 1 microm in diameter and at various radial distances to the seed surface. From the dose distribution of lineal energy, the dose mean lineal energy, y(D), was derived. The results showed that the radiation qualities are constant in the distance range from 0.5 to 5 cm. In this distance range, the quality factor, relative to gamma rays from (60)Co, is 2.2 for (125)I and 2.5 for (103)Pd.

Performance of a commercial optical CT scanner and polymer gel dosimeters for 3-D dose verification

Performance analysis of a commercial three-dimensional (3-D) dose mapping system based on optical CT scanning of polymer gels is presented. The system consists of BANG 3 polymer gels (MGS Research, Inc., Madison, CT), OCTOPUS laser CT scanner (MGS Research, Inc., Madison, CT), and an in-house developed software for optical CT image reconstruction and 3-D dose distribution comparison between the gel, film measurements and the radiation therapy treatment plans. Various sources of image noise (digitization, electronic, optical, and mechanical) generated by the scanner as well as optical uniformity of the polymer gel are analyzed. The performance of the scanner is further evaluated in terms of the reproducibility of the data acquisition process, the uncertainties at different levels of reconstructed optical density per unit length and the effects of scanning parameters. It is demonstrated that for BANG 3 gel phantoms held in cylindrical plastic containers, the relative dose distribution can be reproduced by the scanner with an overall uncertainty of about 3% within approximately 75% of the radius of the container. In regions located closer to the container wall, however, the scanner generates erroneous optical density values that arise from the reflection and refraction of the laser rays at the interface between the gel and the container. The analysis of the accuracy of the polymer gel dosimeter is exemplified by the comparison of the gel/OCT-derived dose distributions with those from film measurements and a commercial treatment planning system (Cadplan, Varian Corporation, Palo Alto, CA) for a 6 cm x 6 cm single field of 6 MV x rays and a 3-D conformal radiotherapy (3DCRT) plan. The gel measurements agree with the treatment plans and the film measurements within the "3%-or-2 mm" criterion throughout the usable, artifact-free central region of the gel volume. Discrepancies among the three data sets are analyzed.

Determining optimal gel sensitivity in optical CT scanning of gel dosimeters

A method for determining the gel sensitivity that is necessary for obtaining optimal image contrast in optical CT scanning of gel dosimeters is presented. The effective dynamic range of the OCTOPUS-ONE research scanner (MGS Research, Inc., Madison, CT) is analyzed. Optical density increments for selected straight-line paths across a gel cylinder to be scanned are calculated based on the optical properties of the polymer gel and the dose distribution from a commercial treatment planning system (Cadplan, Varian Corporation, Palo Alto, CA). Maximum optical density increment across the entire gel is obtained by searching the gel cylinder over a set of transverse planes at different rotational angles. The application of this quantity as a criterion for optimizing the quality of the optical CT scanning is demonstrated through dose verification of two representative treatment plans. When the MU dependence of the dose distribution for a treatment plan is linear, as is the case for static field irradiation, it is possible to scale the treatment plan such that the intensity variation of the signals received by the photodetector spans its entire dynamic range. For treatment plans that are possibly nonlinear, IMRT plans, for example, modification of the sensitivity of the gel material is necessary for the high-dose signals to be collected at a certain signal-to-noise ratio. Results obtained using the optimized CT scanning approach are compared with those from the treatment planning system and the film measurement.

Radiation-induced second cancers: the impact of 3D-CRT and IMRT

Information concerning radiation-induced malignancies comes from the A-bomb survivors and from medically exposed individuals, including second cancers in radiation therapy patients. The A-bomb survivors show an excess incidence of carcinomas in tissues such as the gastrointestinal tract, breast, thyroid, and bladder, which is linear with dose up to about 2.5 Sv. There is great uncertainty concerning the dose-response relationship for radiation-induced carcinogenesis at higher doses. Some animal and human data suggest a decrease at higher doses, usually attributed to cell killing; other data suggest a plateau in dose. Radiotherapy patients also show an excess incidence of carcinomas, often in sites remote from the treatment fields; in addition there is an excess incidence of sarcomas in the heavily irradiated in-field tissues. The transition from conventional radiotherapy to three-dimensional conformal radiation therapy (3D-CRT) involves a reduction in the volume of normal tissues receiving a high dose, with an increase in dose to the target volume that includes the tumor and a limited amount of normal tissue. One might expect a decrease in the number of sarcomas induced and also (less certain) a small decrease in the number of carcinomas. All around, a good thing. By contrast, the move from 3D-CRT to intensity-modulated radiation therapy (IMRT) involves more fields, and the dose-volume histograms show that, as a consequence, a larger volume of normal tissue is exposed to lower doses. In addition, the number of monitor units is increased by a factor of 2 to 3, increasing the total body exposure, due to leakage radiation. Both factors will tend to increase the risk of second cancers. Altogether, IMRT is likely to almost double the incidence of second malignancies compared with conventional radiotherapy from about 1% to 1.75% for patients surviving 10 years. The numbers may be larger for longer survival (or for younger patients), but the ratio should remain the same.

Dosimetry study of Re-188 liquid balloon for intravascular brachytherapy using polymer gel dosimeters and laser-beam optical CT scanner

Angioplasty balloons inflated with a solution of the beta-emitter Re-188 have been used for intravascular brachytherapy to prevent restenosis. Coronary stents are in extensive clinical use for the treatment of de novo atherosclerotic stenoses. In this study, the effect of an interposed stent on the dose distribution has been measured for Re-188 balloon sources using the proprietary BANG polymer gel dosimeters and He-Ne laser-beam optical CT scanner. In polymer gels, after ionizing radiation is absorbed, free-radical chain-polymerization of soluble acrylic monomers occurs to form an insoluble polymer. The BANG polymer gel dosimeters used in these measurements allow high resolution, precise, and accurate three-dimensional determination of dosimetry from a given source. Re-188 liquid balloons, with or without an interposed metallic stent, were positioned inside thin walled tubes placed in such a polymer dosimeter to deliver a prescribed dose (e.g., 15 Gy at 0.5 mm). After removing the balloon source, each irradiated sample was mounted in the optical scanner for scanning, utilizing a single compressed He-Ne laser beam and a single photodiode. In the absence of a stent, doses at points along the balloon axis, at radial distance 0.5 mm from the balloon surface and at least 2.5 mm from the balloon ends, are within 90% of the maximum dose. This uniformity of axial dose is independent of the balloon diameter and length. Dose rate and dose uniformity for intravascular brachytherapy with Re-188 balloon are altered by the presence of stent. The dose reduction by the stent is rather constant (13%-15%) at different radial distances. However, dose inhomogeneity caused by the stent decreases rapidly with radial distance.

3D dosimetry study of 188Re liquid balloon for intravascular brachytherapy using bang polymer gel dosemeters

It has been suggested that the combination of intravascular brachytherapy and coronary stent implantation may result in further reduction of restenosis after percutaneous balloon angioplasty. The use of an angioplasty balloon filled with a 188Re liquid beta source for intravascular brachytherapy provides the advantages of accurate source positioning and uniform dose distribution to the coronary vessel wall. The effect of source edge and stent on the dose distribution of the target tissue may be clinically important. In BANG gels, the absorbed radiation produces free-radical chain polymerisation of acrylic monomers that are initially dissolved in the gel. The number of polymer particles is proportional to the absorbed dose. In this study, 3D dose distributions are presented for 188Re balloons, with and without stents, using a prototype He-Ne laser CT scanner and the proprietary BANG polymer gel dosemeters.

Ultrasonic spectrum-analysis and neural-network classification as a basis for ultrasonic imaging to target brachytherapy of prostate cancer

Conventional B-mode ultrasound is the standard means of imaging the prostate for guiding prostate biopsies and planning brachytherapy of prostate cancer. Yet B-mode images do not allow adequate visualization of cancerous lesions of the prostate. Ultrasonic tissue-typing imaging based on spectrum analysis of radiofrequency echo signals has shown promise for overcoming the limitations of B-mode imaging for visualizing prostate tumors. Tissue typing based on radiofrequency spectrum analysis uses nonlinear methods, such as neural networks, to classify tissue by using spectral-parameter and clinical-variable values. Two- and three-dimensional images based on these methods show potential for improving the guidance of prostate biopsies and the targeting of radiotherapy of prostate cancer. Two-dimensional images have been imported into instrumentation for real-time biopsy guidance and into commercial dose-planning software for brachytherapy planning. Three-dimensional renderings seem to be capable of depicting locations and volumes of cancer foci.

Spectrum-analysis and neural networks for imaging to detect and treat prostate cancer

Conventional B-mode ultrasound currently is the standard means of imaging the prostate for guiding prostate biopsies and planning brachytherapy to treat prostate cancer. Yet B-mode images do not adequately display cancerous lesions of the prostate. Ultrasonic tissue-type imaging based on spectrum analysis of radiofrequency (rf) echo signals has shown promise for overcoming the limitations of B-mode imaging for visualizing prostate tumors. This method of tissue-type imaging utilizes nonlinear classifiers, such as neural networks, to classify tissue based on values of spectral parameter and clinical variables. Two- and three-dimensional images based on these methods demonstrate potential for guiding prostate biopsies and targeting radiotherapy of prostate cancer. Two-dimensional images are being generated in real time in ultrasound scanners used for real-time biopsy guidance and have been incorporated into commercial dosimetry software used for brachytherapy planning. Three-dimensional renderings show promise for depicting locations and volumes of cancer foci for disease evaluation to assist staging and treatment planning, and potentially for registration or fusion with CT images for targeting external-beam radiotherapy.

Factors predicting for postimplantation urinary retention after permanent prostate brachytherapy

PURPOSE

Urinary retention requiring catheterization is a known complication among prostate cancer patients treated with permanent interstitial radioactive seed implantation. However, the factors associated with this complication are not well known. This study was conducted to determine these factors.

METHODS AND MATERIALS

Ninety-one consecutive prostate cancer patients treated with permanent interstitial implantation at our institution from 1996 to 1999 were evaluated. All patients underwent pre-implant ultrasound and postimplant CT volume studies. Isotopes used were (125)I (54 patients) or (103)Pd (37 patients). Twenty-three patients were treated with a combination of 45 Gy of external beam radiation therapy as well as seed implantation, of which only 3 patients were treated with (125)I. Mean pretreatment prostate ultrasound volume was 35.4 cc (range, 10.0-70.2 cc). The mean planning ultrasound target volume (PUTV) was 39.6 cc (range, 16.1-74.5 cc), whereas the mean posttreatment CT target volume was 55.0 cc (range, 20.2-116 cc). Patient records were reviewed to determine which patients required urinary catheterization for relief of urinary obstruction. The following factors were analyzed as predictors for urinary retention: clinical stage; Gleason score; prostate-specific antigen; external beam radiation therapy; hormone therapy; pre-implant urinary symptoms (asymptomatic/nocturia x 1 vs. more significant urinary symptoms); pretreatment ultrasound prostate volume; PUTV; PUTV within the 125%, 150%, 200%, 250%, 300% isodose lines; postimplant CT volume within the 125%, 150%, 200%, 250%, 300% isodose lines; D90; D80; D50; ratio of post-CT volume to the PUTV; the absolute change in volume between the CT volume and PUTV; number of needles used; activity per seed; and the total activity of the implant. Statistical analyses using logistic regression and chi2 were performed.

RESULTS

Eleven of 91 (12%) became obstructed. Significant factors predicting for urinary retention were the total number of needles used (p < 0.038); the pretreatment ultrasound prostate volume (p < 0.048); the PUTV (p < 0.02); and the posttreatment CT volume (p < 0.021). Two of 51 patients (3.9%) requiring 33 or fewer needles (median) experienced obstruction vs. 9 of 40 (22.5%) requiring more than 33 (p < 0.007). If the pretreatment ultrasound prostate volume was 35 cc or less (median), 3 of 43 (7%) vs. 8 of 36 (22%) with a volume greater than 35 cc experienced obstruction (p < 0.051).

CONCLUSION

The number of needles required (perhaps related to trauma to the prostate) and the prostate volumes were significant factors predicting for urinary retention after permanent prostate seed implantation.

Dosimetric and volumetric criteria for selecting a source activity and a source type ((125)I or (103)Pd) in the presence of irregular seed placement in permanent prostate implants

PURPOSE

The dosimetric merit of a permanent prostate implant relies on two factors: the quality of the plan itself, and the fidelity of its implementation. The former factor depends on source type and on source strength, while the latter is a combination of skill and experience. The purpose of this study is to offer criteria by which to select a source type ((125)I or (103)Pd) and activity.

METHODS AND MATERIALS

Given a prescription dose and potential seed positions along needles, treatment plans were designed for a number of seed types and activities, specifically for (125)I with activities ranging from 0.3 to 0.7 mCi, and for (103)Pd with activities in the range of 0.8 to 1.6 mCi. To avoid human planner bias, an automated computerized planning system based on integer programming was used to obtain optimal seed configurations for each seed type and activity. To simulate the effect of seed-placement inaccuracies, random seed-displacement "errors" were generated for all plans. The displacement errors were assumed to be uniformly distributed within a cube with side equal to 2sigma. The resulting treatment plans were assessed using two volumetric and two dosimetric indices.

RESULTS

For (125)I implants a coverage index (CI) of 98.5% or higher can be achieved for all activities (CI is the fraction of the target volume receiving the prescribed or larger dose). The external volume index (EI) (i.e., the amount of healthy tissue, as percentage of the target volume, receiving the prescribed or larger dose) increases from 13.9% to 20% as the activity increases from 0.3 to 0.7 mCi. For implants using (103)Pd, the external volume index increases from 10. 2% to 13.9% whenever CI exceeds 98.5%. Volumetric and dosimetric indices (coverage index, external volume index, D90, and D80) are all sensitive to seed displacement, although the activity dependence of these indices is more pronounced for (125)I than for (103)Pd implants.

CONCLUSIONS

For both isotopes, the lower activities studied systematically result in lower EIs. If seeds can be placed within approximately 0.5 cm of their intended position (103)Pd should be preferred because its EI is lower than that of (125)I. For all activities the coverage indices and D90 are within the required range. If seed placement uncertainties are larger than 0.5 cm, (125)I provides slightly better target coverage; however, in terms of external volume (healthy tissue) covered, (103)Pd is superior to (125)I.

Potential decreased morbidity of interstitial brachytherapy for gynecologic malignancies using laparoscopy: A pilot study

OBJECTIVES

This pilot study was designed to prospectively assess whether the addition of laparoscopy at the time of interstitial brachytherapy is safe, provides verification and/or guidance of needle placement, and results in a reduction of treatment-related morbidity.

METHODS

Between 7/93 and 2/97 15 consecutive eligible patients were entered into this study. All patients received external pelvic radiation to a dose range between 45 and 61.20 Gy using 1.8-Gy fractions. In each patient the minimum prescribed dose for the brachytherapy portion was 20 Gy. Minimum cumulative doses to sites of gross disease ranged from 71.8 to 115.3 Gy. A Syed-Neblett afterloading perineal template was used in all the procedures. Laparoscopy using established guidelines was performed during placement of interstitial needles. During template placement, verification of interstitial needles on laparoscopy and any subsequent changes or needle rearrangement were noted.

RESULTS

No acute radiation toxicity greater than Grade 2 was noted during the external beam portion of treatment, and no perioperative complications were encountered. These needles were withdrawn under laparoscopic guidance to just below the peritoneal reflection, avoiding proximity to the bowel and improving tumor coverage. Median follow-up time was 26 months. No late radiation morbidity greater than Grade 2 nor any laparoscopic-related complications were noted. To date, one patient has died of disease; six are alive with disease; and eight are alive free of disease with a mean disease-free survival of 17.3 months.

CONCLUSION

Laparoscopy at the time of interstitial brachytherapy appears to be safe. No radiation toxicity greater than Grade 2 has developed. No perioperative complications were seen with the addition of laparoscopy. The addition of laparoscopy to the placement of transperineal interstitial implants impacted needle arrangement and/or loading of sources in 50% of patients.

Treatment planning for brachytherapy: an integer programming model, two computational approaches and experiments with permanent prostate implant planning

An integer linear programming model is proposed as a framework for optimizing seed placement and dose distribution in brachytherapy treatment planning. The basic model involves using 0/1 indicator variables to describe the placement or non-placement of seeds in a prespecified three-dimensional grid of potential locations. The dose delivered to each point in a discretized representation of the diseased organ and neighbouring healthy tissue can then be modelled as a linear combination of the indicator variables. A system of linear constraints is imposed to attempt to keep the dose level at each point to within specified target bounds. Since it is physically impossible to satisfy all constraints simultaneously, each constraint uses a variable to either record when the target dose level is achieved, or to record the deviation from the desired level. These additional variables are embedded into an objective function to be optimized. Variations on this model are discussed and two computational approaches--a branch-and-bound algorithm and a genetic algorithm--for finding 'optimal' seed placements are described. Results of computational experiments on a collection of prostate cancer cases are reported. The results indicate that both optimization algorithms are capable of producing good solutions within 5 to 15 min, and that small variations in model parameters can have a measurable effect on the dose distribution of the resulting plans.

A calculation of the relative biological effectiveness of 125I and 103Pd brachytherapy sources using the concept of proximity function

The clinical application of encapsulated radioactive sources in brachytherapy plays an important role in the treatment of malignancy. 125I and 103Pd sources have been widely used in the permanent implant of prostate cancer. An important consideration for the choice of brachytherapy sources is their relative biological effectiveness (RBE). Previous calculations of this quantity have used the dose-averaged lineal energy, yD, as a measure of biological effectiveness. In this approach, however, the selection of a relevant site size remains an open question. Here we avoid this problem by using the generalized theory of dual radiation action to calculate the initial slope, alpha, of the dose-effect curves using the proximity function, t(x), and the biological response function, gamma(x). At low doses and/or low dose rates (e.g., prostate implants) the parameter alpha determines the RBE. Proximity function, t(x), is the probability distribution function of distances between pairs of sublesions; and the biological function, gamma(x), is the probability that two sublesions at a distance x apart results in a lesion. Functions t(x) have been calculated for each source using the Monte Carlo transport codes PHOEL and PROTON5. The function gamma(x) has been taken from a published analysis. The RBE values thus obtained are: 1.5 for 125I and 1.6 for 103Pd. The question of whether an "effective" site size exists where yD approximates best the variation of alpha with radiation quality is also addressed.

Microdosimetric evaluation of relative biological effectiveness for 103Pd, 125I, 241Am, and 192Ir brachytherapy sources

PURPOSE

To determine the microdosimetric-derived relative biological effectiveness (RBE) of 103Pd, 125I, 241Am, and 192Ir brachytherapy sources at low doses and/or low dose rates.

METHODS AND MATERIALS

The Theory of Dual Radiation Action can be used to predict expected RBE values based on the spatial distribution of energy deposition at microscopic levels from these sources. Single-event lineal energy spectra for these isotopes have been obtained both experimentally and theoretically. A grid-defined wall-less proportional counter was used to measure the lineal energy distributions. Unlike conventional Rossi proportional counters, the counter used in these measurements has a conducting nylon fiber as the central collecting anode and has no metal parts. Thus, the Z-dependence of the photoelectric effect is eliminated as a source of measurement error. Single-event spectra for these brachytherapy sources have been also calculated by: (a) the Monte Carlo code MCNP to generate the electron slowing down spectrum, (b) transport of monoenergetic electron tracks, event by event, with our Monte Carlo code DELTA, (c) using the concept of associated volume to obtain the lineal energy distribution f(y) for each monoenergetic electron, and (d) obtaining the composite lineal energy spectrum for a given brachytherapy source based on the electron spectrum calculated at step (a).

RESULTS

Relative to 60Co, the RBE values obtained from this study are: 2.3 for 103Pd, 2.1 for 125I, 2.1 for 241Am, and 1.3 for 192Ir.

CONCLUSIONS

These values are consistent with available data from in vitro cell survival experiments. We suggest that, at least for these brachytherapy sources, microdosimetry may be used as a credible alternative to time-consuming (and often uncertain) radiobiological experiments to obtain information on radiation quality and make reliable predictions of RBE in low dose rate brachytherapy.

The combined effects of sublethal damage repair, cellular repopulation and redistribution in the mitotic cycle. I. Survival probabilities after exposure to radiation

An analytical model is presented that describes radiation-induced cellular inactivation in the presence of sublethal damage repair, cellular repopulation and redistribution in the mitotic cycle (the 3 Rs). The parameters of the model are measurable experimentally. Also taken into account are the initial age distribution of the cell population, the fact that subgroups of cells progress through the cycle at different speeds, the effects of a dose of radiation on the duration of the four phases of the cycle (G1, S, G2, M), the possibility that a certain fraction of the cells are quiescent, and cell loss and/or cell removal from the proliferating population. Survival probabilities are expressed as linear-quadratic functions of dose where the coefficient alpha and beta as well as the recovery constant (t0) are taken to depend on the position of the cell in the mitotic cycle. Explicit analytical expressions for inactivation probability are given for clonogenic cells exposed to continuous or fractionated radiation. Two model calculations are used to illustrate the formalism: in one, the redistribution of cells during fractionated therapy is examined. In the other calculation, it is shown that it is sufficient to take into account differences in proliferation rates and the change in the ratio alpha/beta within the generation cycle for cells that may have otherwise equal response to acute exposures to explain that in a fractionated treatment protocol late-responding cells are more sensitive to the dose per fraction than early-responding cells. It is not necessary to invoke differences in radiosensitivity between these two classes of cells.

Application of the HSEF to assessing radiation risks in the practice of radiation protection

The primary risk coefficients upon which exposure limits for radiation protection purposes are currently based are derived almost exclusively from cancer-induction data obtained from human populations exposed to radiations of low linear energy transfer. The question of higher linear energy transfer radiations is handled by means of quality factors derived from values for relative biological effectiveness obtained from animal data. However, the advent of microdosimetry has made it possible to establish hit size effectiveness functions from single-cell systems, both in vitro and in vivo. This type of function can substitute completely for the concept of relative biological effectiveness, Q and equivalent dose. A common basis for risk coefficients and the hit size effectiveness function lies in the fact that human cancers are monoclonal and thus single cell in origin. The present communication utilizes this common base as a means of extending the present low-linear energy transfer based risk coefficients to include carcinogenic responses from exposure in radiation fields of any one or mixed qualities, extending from the smallest to the largest linear energy transfers of practical consequence. In doing so, risks from ionizing radiations of any linear energy transfer may be predicted more accurately than at present.

The effects of sublethal damage recovery and cell cycle progression on the survival probability of cells exposed to radioactive sources

Cell progression through the mitotic cycle during low dose rate irradiation may alter notably the survival probability, particularly when a fraction of the dose is delivered during a sensitive phase of the cycle. In this paper we indicate that the consequences of this phenomenon, commonly believed to lead to an "inverse dose rate effect", may be significantly modulated (and even cancelled) as a result of (a) interactions among sublethal lesions produced in different phases of the mitotic cycle, and (b) variations in these lesions' production rates and repair ability from one phase of the cycle to another. The mathematical model presented (and accompanying numerical examples) takes into account the possibility of changes (e.g. radioactive decay) in the dose rate during exposure.

Rapid dose calculations for stereotactic radiosurgery

A three-dimensional dose calculation algorithm is described for stereotactic radiosurgery using multiple noncoplanar beam arcs. Precalculated dose libraries of 20-deg arc segments, or mini arcs, are stored in computer memory which permits rapid calculation of complete, high resolution, three-dimensional isodose distributions and dose volume histograms. Three-dimensional patient contours and target volumes are obtained from CT scans and angiographic x rays. Rapid dose calculations are made possible by the use of arc libraries and an improved algorithm for mapping beam doses to the dose calculation grid. This permits more flexibility in designing optimum treatment plans, as five-six complete plans can be generated in less than 1 h. Thus many possible treatment options can be tested in the 3-4-h time period typically available in stereotactic procedures between CT scanning and treatment.

Microdosimetry for boron neutron capture therapy

Preclinical studies for boron neutron capture therapy (BNCT) using epithermal neutrons are ongoing at several laboratories. The absorbed dose in tumor cells is a function of the thermal neutron flux at depth, the microscopic boron concentration, and the size of the cell. Dosimetry is therefore complicated by the admixture of thermal, epithermal, and fast neutrons, plus gamma rays, and the array of secondary high-linear-energy-transfer particles produced within the patient from neutron interactions. Microdosimetry can be a viable technique for determining absorbed dose and radiation quality. A 2.5-cm-diameter tissue-equivalent gas proportional counter has been built with 50 parts per million (ppm) 10B incorporated into the walls and counting gas to simulate the boron uptake anticipated in tumors. Measurements of lineal energy (y) spectra for BNCT in simulated volumes of 1-10 microns diameter show a dose enhancement factor of 4.3 for 30 ppm boron, and a "y" of 250 keV/microns for the boron capture process. Chamber design plus details of experimental and calculated linear energy spectra will be presented.

Determination of the kerma factors in tissue-equivalent plastic, C, Mg, and Fe for 14.7-MeV neutrons

Microdosimetric measurements were made with tissue-equivalent plastic (TEP), C-, Mg-, and Fe-walled proportional counters filled with propane-based tissue equivalent (TE) gas and Ar gas and irradiated with 14.7-MeV neutrons. A theoretical model was used for the analysis of energy deposition in spherical detectors. An effective average mass stopping-power ratio and a W correction were calculated to convert the gas ionization to the kerma in the wall material. The neutron fluence at the position of microdosimetric measurements was determined with an associated particle chamber mounted with surface barrier detectors. The experimental measurements along with the calculated correction factors yielded kerma factors of 0.660 X 10(-8) cGy cm2 for TEP, 0.219 X 10(-8) cGy cm2 for C, 0.122 X 10(-8) cGy cm2 for Mg, and 0.479 X 10(-9) cGy cm2 for Fe. The estimated uncertainties are 8.0% for TEP, 10.5% for C, and 9.3% for Mg and Fe.

The direct determination of dose-to-water using a water calorimeter

A flexible, temperature-regulated, water calorimeter has been constructed which consists of three nested cylinders. The innermost "core" is a 10 X 10 cm right cylinder made of glass, the contents of which are isolated from the environment. It has two Teflon-washered glass valves for filling, and two thermistors are supported at the center by glass capillary tubes. Surrounding the core is a "jacket" that provides approximately 2 cm of air insulation between the core and the "shield." The shield surrounds the jacket with a 2.5-cm layer of temperature-regulated water flowing at 51/min. The core is filled with highly purified water the gas content of which is established prior to filling. Convection currents, which may be induced by dose gradients or thermistor power dissipation, are eliminated by operating the calorimeter at 4 degrees C. Depending upon the power level of the thermistors, 15-200 microW, and the insulation provided by the glass capillary tubing, the temperature of the thermistors is higher than that of the surrounding water. To minimize potential errors caused by differences between calibration curves obtained at finite power levels, the zero-power-level calibration curve obtained by extrapolation is employed. Also the calorimeter response is corrected for the change in power level, and therefore thermistor temperature, that follows the resistance change caused by irradiation. The response of the calorimeter to 4-MV x rays has been compared to that of an ionization chamber irradiated in an identical geometry.(ABSTRACT TRUNCATED AT 250 WORDS)