Conventional versus virtual simulation for radiation treatment planning of prostate cancer: final results.

PURPOSE: Radiotherapy is widely used to treat patients with prostate cancer. Using conventional x-ray simulation is often difficult to accurately localize the extent of the tumor, to cover exactly the lymph nodes at risk and shield the organs at risk. We report on the results of a study comparing target localization with conventional and virtual simulation. METHODS: One hundred prostate cancer patients underwent both conventional and virtual simulation. The conventional simulation films were compared with digitally reconstructed radiographs (DDRs) produced from the computed tomography (CT) data. All patients underwent target localization for radical prostate radiotherapy. The treatment fields were initially marked with a conventional portal film on linear accelerator (LINAC), plain x-ray film and available diagnostic imaging. Each patient then had a CT and these simulated treatment fields were reproduced within the virtual simulation planning system. The treatment fields defined by the clinicians using each modality were compared in terms of field area and implications for target coverage. RESULTS: Virtual simulation showed significantly greater clinical tumor volume coverage and less normal tissue volume irradiated compared with conventional simulation (p <0.001). CONCLUSION: CT localization and virtual simulation allow more accurate definition of the clinical target volume. This could enable a reduction in geographical misses, reducing at the same time treatment-related toxicity.

Comparison of conventional and virtual simulation for radiation treatment planning of prostate cancer.

PURPOSE: radiotherapy is widely used to treat patients with prostate cancer. Using conventional x-ray simulation is often difficult to accurately localize the extent of the tumor, to cover exactly the lymph nodes at risk and shield the organs at risk. We report the initial results of a study conducted to compare target localization with conventional and virtual simulation. METHODS: fifty patients with prostate cancer underwent target localization for radical prostate radiotherapy using conventional and virtual simulation. The treatment fields were initially marked with a conventional portal film on LINAC, plain x-ray film and available diagnostic imaging. Each patient then had a computed tomography (CT) and these simulated treatment fields were reproduced within the virtual simulation planning system. The treatment fields defined by the clinicians using each modality were compared in terms of field area and implications for target coverage. RESULTS: there was significantly greater clinical tumor volume coverage using virtual simulation compared with conventional simulation and less normal tissue volume irradiated (p<0.001). CONCLUSION: CT localization and virtual simulation allow for more accurate definition of the clinical target volume. This could enable a reduction in geographical misses, while also reducing treatment-related toxicity.

Gamma knife output factor measurements using VIP polymer gel dosimetry.

PURPOSE: Water equivalent polymer gel dosimeters and magnetic resonance imaging were employed to measure the output factors of the two smallest treatment fields available in a Gamma Knife model C radiosurgery unit, those formed employing the 4 and 8 mm final collimator helmets. METHODS: Three samples of the VIP normoxic gel formulation were prepared and irradiated so that a single shot of the field whose output factor is to be measured and a single shot of the reference 18 mm field were delivered in each one. Emphasis is given to the development and benchmarking of a refined data processing methodology of reduced uncertainty that fully exploits the 3D dose distributions registered in the dosimeters. RESULTS: Polymer gel results for the output factor of the 8 mm collimator helmet are found to be in close agreement with the corresponding value recommended by the vendor (0.955 +/- 0.007 versus 0.956, respectively). For the 4 mm collimator helmet, however, polymer gel results suggest an output factor 3% lower than the value recommended by the vendor (0.841 +/- 0.009 versus 0.870, respectively). CONCLUSIONS: A comparison with corresponding measurements published in the literature indicates that output factor results of this work are in agreement with those obtained using dosimetric systems which, besides fine spatial resolution and lack of angular and dose rate dependence of the dosimeter's response, share with polymer gels the favorable characteristic of minimal radiation field perturbation.

Radiobiological modelling of dose-gradient effects in low dose rate, high dose rate and pulsed brachytherapy.

This paper presents a generalization of a previously published methodology which quantified the radiobiological consequences of dose-gradient effects in brachytherapy applications. The methodology uses the linear-quadratic (LQ) formulation to identify an equivalent biologically effective dose (BED(eq)) which, if applied uniformly to a specified tissue volume, would produce the same net cell survival as that achieved by a given non-uniform brachytherapy application. Multiplying factors (MFs), which enable the equivalent BED for an enclosed volume to be estimated from the BED calculated at the dose reference surface, have been calculated and tabulated for both spherical and cylindrical geometries. The main types of brachytherapy (high dose rate (HDR), low dose rate (LDR) and pulsed (PB)) have been examined for a range of radiobiological parameters/dimensions. Equivalent BEDs are consistently higher than the BEDs calculated at the reference surface by an amount which depends on the treatment prescription (magnitude of the prescribed dose) at the reference point. MFs are closely related to the numerical BED values, irrespective of how the original BED was attained (e.g., via HDR, LDR or PB). Thus, an average MF can be used for a given prescribed BED as it will be largely independent of the assumed radiobiological parameters (radiosensitivity and alpha/beta) and standardized look-up tables may be applicable to all types of brachytherapy treatment. This analysis opens the way to more systematic approaches for correlating physical and biological effects in several types of brachytherapy and for the improved quantitative assessment and ranking of clinical treatments which involve a brachytherapy component.

Radiation risk assessment in neonatal radiographic examinations of the chest and abdomen: a clinical and Monte Carlo dosimetry study.

Seeking to assess the radiation risk associated with radiological examinations in neonatal intensive care units, thermo-luminescence dosimetry was used for the measurement of entrance surface dose (ESD) in 44 AP chest and 28 AP combined chest-abdominal exposures of a sample of 60 neonates. The mean values of ESD were found to be equal to 44 +/- 16 microGy and 43 +/- 19 microGy, respectively. The MCNP-4C2 code with a mathematical phantom simulating a neonate and appropriate x-ray energy spectra were employed for the simulation of the AP chest and AP combined chest-abdominal exposures. Equivalent organ dose per unit ESD and energy imparted per unit ESD calculations are presented in tabular form. Combined with ESD measurements, these calculations yield an effective dose of 10.2 +/- 3.7 microSv, regardless of sex, and an imparted energy of 18.5 +/- 6.7 microJ for the chest radiograph. The corresponding results for the combined chest-abdominal examination are 14.7 +/- 7.6 microSv (males)/17.2 +/- 7.6 microSv (females) and 29.7 +/- 13.2 microJ. The calculated total risk per radiograph was low, ranging between 1.7 and 2.9 per million neonates, per film, and being slightly higher for females. Results of this study are in good agreement with previous studies, especially in view of the diversity met in the calculation methods.

Technique charts for Kodak EC-L film screen system for portal localization in a 6MV X-ray beam.

Port films are used in radiotherapy for visual evaluation of the radiation fields and subsequent quantitative analysis. Common port films suffer from poor image quality compared to the simulator-diagnostic films and is desirable to determine the appropriate exposure required for the best image contrast. The aim of this work is to generate technique charts for the Kodak EC-L film screen system for use in a 6MV x-ray beam. Three homogeneous water phantoms were used to simulate head-neck, thorax and abdomen dimensions of adult human, correspondingly. The film screen system was calibrated in a 6MV x-ray beam and under various irradiation conditions. The film screen system behavior was studied as a function of phantom thickness, field size and air gap between the phantom and the film screen system. In each case the optimum film exposure which produces the maximum image contrast was determined. The generated technique charts for the EC-L film screen system and for a 6 MV x-ray beam are used in our radiotherapy department for daily quality assurance of the radiotherapy procedure.

Monte Carlo and thermoluminescence dosimetry of the new IsoSeed model I25.S17 125I interstitial brachytherapy seed.

Monte Carlo simulation and experimental thermoluminescence dosimetry were utilized for the dosimetric characterization of the new IsoSeed model I25.S17 125I interstitial brachytherapy seed. The new seed design is similar to that of the selectSeed and 6711 seeds, with the exception of its molybdenum marker. Full dosimetric data are presented following the recommendations in the Update of the AAPM Task Group 43 report (TG-43U1). A difference of 3.3% was found between Monte Carlo dose rate constant results calculated by air kerma strengths from simulations using a point detector and a detector resembling the solid angle subtended to the seed by the Wide Angle Free Air Chamber (WAFAC) in the primary standard calibration geometry. Following the TG-43U1 recommendations, an average value of lambdaMC = (0.929 +/- 0.014) cGy h(-1) U(-1) was adopted for the new seed. This value was then averaged with the measured value of lambdaEXP = (0.951 +/- 0.044) cGy h(-1) U(-1) to yield the proposed dose rate constant for the new seed that is equal to lambda = (0.940 +/- 0.051) cGy h(-1) U(-1). The Monte Carlo calculated radial dose function and two-dimensional (2-D) anisotropy function results for the new seed were found in agreement with experimental results to within statistical uncertainty of repeated measurements. Monte Carlo simulations were also performed for 125I seeds of similar geometry and dimensions for the purpose of comparison. The new seed presents dosimetric characteristics that are very similar to that of the selectSeed. In comparison to the most extensively studied Amersham 6711 seed, the new one presents similar dosimetric characteristics with a slightly reduced dose rate constant (1.5%).

An evaluation of the TSE MR sequence for time efficient data acquisition in polymer gel dosimetry of applications involving high doses and steep dose gradients.

The use of magnetic resonance imaging as a readout method for polymer gel dosimetry commonly involves long imaging sessions, particularly when high spatial resolution is required in all three dimensions, for the investigation of dose distributions with steep dose gradients and stringent dose delivery specifications. In this work, a volume selective turbo spin echo (TSE) pulse sequence is compared to the established Carr-Purcell-Meiboom-Gill (CPMG) multiecho acquisition with regard to providing accurate dosimetric results in significantly reduced imaging times. Polyethylene glycol diacrylate based (PABIG) gels were irradiated and subsequently scanned to obtain R2 relaxation rate measurements, using a CPMG multiecho sequence and a dual echo TSE utilizing an acceleration (turbo) factor of 64. R2 values, plotted against corresponding Monte Carlo dose calculations, provided calibration data of PABIG gels dose response over a wide dose range. A linear R2 versus dose relationship was demonstrated for both sequences with TSE results presenting reduced dose sensitivity. Although TSE data were found to deviate from linearity at lower doses compared to CPMG data, a relatively wide dynamic dose range of response extending up to approximately 100 Gy was observed for both sequences. The TSE and CPMG sequences were evaluated with a brachytherapy irradiation using a high dose rate 192Ir source and a gamma knife stereotactic radiosurgery irradiation with a single 4 mm collimator helmet shot. Dosimetric results obtained with the TSE and CPMG are shown to compare equally well with the expected dose distributions for these irradiations. The 60-fold scan time reduction achieved with TSE implies that this sequence could prove to be a useful tool for the introduction of polymer gel dosimetry in clinical radiation therapy applications involving high doses and steep dose gradients.

Polymer gel dosimetry close to an 125I interstitial brachytherapy seed.

Despite its advantages, the polymer gel-magnetic resonance imaging (MRI) method has not, as yet, been successfully employed in dosimetry of low energy/low dose rate photon-emitting brachytherapy sources such as 125I or 103Pd interstitial seeds. In the present work, two commercially available 125I seed sources, each of approximately 0.5 U, were positioned at two different locations of a polymer gel filled vial. The gel vial was MR scanned with the sources in place 19 and 36 days after seed implantation. Calibration curves were acquired from the coupling of MRI measurements with accurate Monte Carlo dose calculations obtained simulating the exact experimental setup geometry and materials. The obtained gel response data imply that while linearity of response is sustained, sensitivity (calibration curve slope) is significantly increased (approximately 60%) compared to its typical value for the 192Ir (or 60Co and 6 MV LINAC) photon energies. Water equivalence and relative energy response corrections of the gel cannot account for more than 3-4% of this increase, which, therefore, has to be mainly attributed to physicochemical processes related to the low dose rate of the sources and the associated prolonged irradiation time. The calibration data obtained from one 125I source were used to provide absolute dosimetry results for the other 125I source, which were found to agree with corresponding Monte Carlo calculations within experimental uncertainties. It is therefore suggested that, regardless of the underlying factors accounting for the gel dose response to 125I irradiations, polymer gel dosimetry of new 125I or 103Pd sources should be carried out as originally proposed by Heard and Ibbot (2004 J. Phys.: Conf. Ser. 3 221-3), i.e., by irradiating the same gel sample with the new low dose rate source, as well as with a well-characterized low dose rate source which will provide the dose calibration curve for the same irradiation conditions.

Dose verification of single shot gamma knife applications using VIPAR polymer gel and MRI.

This work describes an experimental procedure with potential to assess the overall accuracy associated with gamma knife clinical applications, from patient imaging and dosimetry planning to patient positioning and dose delivery using the automated positioning system of a Leksell Gamma Knife model C. The VIPAR polymer gel-MRI dosimetry method is employed due to its inherent three-dimensional feature and linear dose response over the range of gamma knife applications. Different polymer gel vials were irradiated with single shot gamma knife treatment plans using each of the four available collimator helmets to deliver a maximum dose of 30 Gy. Percentage relative dose results are presented not only in the form of one-dimensional profiles but also planar isocontours and isosurfaces in three dimensions. Experimental results are compared with corresponding Gammaplan treatment planning system calculations as well as acceptance test radiochromic film measurements. A good agreement, within the experimental uncertainty, is observed between measured and expected dose distributions. This experimental uncertainty is of the order of one imaging pixel in the MRI gel readout session (<1 mm) and allows for the verification of single shot gamma knife applications in terms of acceptance specifications for precision in beam alignment and accuracy. Averaging net R(2) results in the dose plateau of the 4 mm and 18 mm collimator irradiated gel vials, which were MR scanned in the same session, provides a crude estimate of the 4 mm output factor which agrees within errors with the default value of 0.870.

Three-dimensional dose verification of the clinical application of gamma knife stereotactic radiosurgery using polymer gel and MRI.

This work seeks to verify multi-shot clinical applications of stereotactic radiosurgery with a Leksell Gamma Knife model C unit employing a polymer gel-MRI based experimental procedure, which has already been shown to be capable of verifying the precision and accuracy of dose delivery in single-shot gamma knife applications. The treatment plan studied in the present work resembles a clinical treatment case of pituitary adenoma using four 8 mm and one 14 mm collimator helmet shots to deliver a prescription dose of 15 Gy to the 50% isodose line (30 Gy maximum dose). For the experimental dose verification of the treatment plan, the same criteria as those used in the clinical treatment planning evaluation were employed. These included comparison of measured and GammaPlan calculated data, in terms of percentage isodose contours on axial, coronal and sagittal planes, as well as 3D plan evaluation criteria such as dose-volume histograms for the target volume, target coverage and conformity indices. Measured percentage isodose contours compared favourably with calculated ones despite individual point fluctuations at low dose contours (e.g., 20%) mainly due to the effect of T2 measurement uncertainty on dose resolution. Dose-volume histogram data were also found in a good agreement while the experimental results for the percentage target coverage and conformity index were 94% and 1.17 relative to corresponding GammaPlan calculations of 96% and 1.12, respectively. Overall, polymer gel results verified the planned dose distribution within experimental uncertainties and uncertainty related to the digitization process of selected GammaPlan output data.

Technical note: evaluation of dosimetric performance in a commercial 3D treatment planning system.

The aim of this work was to evaluate the dosimetric performance of a commercial treatment planning system (TPS) which employs a three-dimensional calculation algorithm (Nucletron Plato version 2.2.3), following the guidelines of the AAPM Task Group 23 (TG23). Seven test cases were used to test the TPS dosimetric performance in homogeneous water. These cases involved absolute dose measurements on central as well as off-axis points situated at various depths, using simple field arrangements, and comparison with corresponding TPS calculations. This comparison yielded differences within +/-2% at all points, for all test cases. To test the ability of the TPS to account for tissue inhomogeneities, corresponding comparisons were performed with the presence of a low-density material in the beam to resemble an air inhomogeneity. Absolute dose measurements and corresponding TPS calculations showed a mean deviation of the order of +/-3.5%, reaching a maximum of 11.5% for small field sizes (5 cm x 5 cm). In summary, observed deviations are well within the set tolerance levels while comparison with previous TPS versions showed that Plato version 2.2.3 is significantly improved, especially in dose calculations in the presence of low density inhomogeneities.

Entrance radiation doses during paediatric cardiac catheterisations performed for diagnosis or the treatment of congenital heart disease.

The purpose of this study was to estimate the radiation exposure of children, during cardiac catheterisations for the diagnosis or treatment of congenital heart disease. Radiation doses were estimated for 45 children aged from 1 d to 13 y old. Thermoluminescent dosemeters (TLDs) were used to estimate the posterior entrance dose (DP), the lateral entrance dose (DLAT), the thyroid dose and the gonads dose. A dose-area product (DAP) meter was also attached externally to the tube of the angiographic system and gave a direct value in mGy cm2 for each procedure. Posterior and lateral entrance dose values during cardiac catheterisations ranged from 1 to 197 mGy and from 1.1 to 250.3 mGy, respectively. Radiation exposure to the thyroid and the gonads ranged from 0.3 to 8.4 mGy to 0.1 and 0.7 mGy, respectively. Finally, the DAP meter values ranged between 360 and 33,200 mGy cm2. Radiation doses measured in this study are comparable with those reported to previous studies. Moreover, strong correlation was found between the DAP values and the entrance radiation dose measured with TLDs.

Fetal dose evaluation during breast cancer radiotherapy.

PURPOSE: The aim of the work was to estimate the radiation dose delivered to the fetus in a pregnant patient irradiated for breast cancer. METHODS AND MATERIALS: A 45-year woman was treated for left breast cancer using a 6 MV photon beam with two isocentric opposing tangential unwedged fields. Daily dose was 2.3 Gy at 95% isodose line given by two fields/day, 5 days/week. A total dose of 46 Gy was given in 20 fractions over a 4-week period. Pregnancy confirmed during the second therapeutic week. Treatment lasted between the second and sixth gestation week. Radiation dose to fetus was estimated from in vivo and phantom measurements using thermoluminescence dosimeters and an ionization chamber. In vivo measurements were performed by inserting either a catheter with TL dosimeters or ionization chamber into the patient's rectum. Phantom measurements were performed by simulating the treatment conditions on an anthropomorphic phantom. RESULTS: TLD measurements (in vivo and phantom) revealed fetal dose to be 0.085% of the tumor dose, corresponding to a cumulative fetal dose of 3.9 cGy for the entire treatment of 46 Gy. Chamber measurements (in vivo and phantom) revealed a fetal dose less than the TLD result: 0.079 and 0.083% of the tumor dose corresponding to cumulative fetal dose of 3.6 cGy and 3.8 cGy for in vivo and phantom measurement, respectively. CONCLUSIONS: It was concluded that the cumulative dose delivered to the unshielded fetus was 3.9 cGy for a 46 Gy total tumor dose. The estimated fetal dose is low compared to the total tumor dose given due to the early stage of pregnancy, the large distance between fundus-radiation field, and the fact that no wedges and/or lead blocks were used. No deterministic biological effects of radiation on the live-born embryo are expected. The lifetime risk for radiation-induced fatal cancer is higher than the normal incidence, but is considered as inconsequential.

Combined teletherapy and intracavitary brachytherapy boost for the treatment of nasopharyngeal carcinoma.

For the non-metastatic nasopharyngeal carcinoma, external beam radiation therapy (median dose 64 Gy) and a boost of intracavitary irradiation (ICRT) has been given. Caesium-137 pellets of 40 mCi were used at a dose rate of 3-3.5 Gy/h, 1 cm from the sources. The median dose was 8.5 Gy. Overall 5-year actuarial survival for the 48 studied patients was 60.4% and LRFS was 64%. The procedure was well tolerated by our patients.

Limitations of the point and line source approximations for the determination of geometry factors around brachytherapy sources.

Geometry factors were calculated around commercially available pellets and elongated brachytherapy sources taking into account their actual active core geometries. These calculations were compared with corresponding ones derived using the point and line source approximations commonly used for the determination of geometry factors, as proposed by AAPM Task Group 43. The point source approximation was found to be efficient for the determination of geometry factors around single active pellets, even at radial distances very close to the source. It is also valid for the determination of geometry factors at radial distances r>2L from elongated brachytherapy sources of length L. For smaller radial distances, however, this approximation introduces significant errors, >50%, around elongated source designs, thus being unacceptable for the determination of geometry factors in this case. The line source approximation was found to accurately reproduce geometry factors around elongated brachytherapy source designs. Errors greater than 3%, due to the fact that the line source approximation ignores the radial dimension d of the source, are observed only at radial distances very close to the source (r< or =L/2), at polar angles far away from their transverse bisectors. These errors depend on the ratio d/L and increase as this ratio increases.

Monte Carlo dosimetry of a new 192Ir high dose rate brachytherapy source.

This work provides full dosimetric data for a new high-strength 192Ir source currently launched by Varian Oncology Systems for use in their high dose rate remote afterloading systems. The active core length of the new source is reduced to 5 mm compared to a value of 10 mm for the existing VariSource source design, with all other geometric source and encapsulation details being similar. Dose-rate constant, radial dose functions, geometry factors, and anisotropy functions, utilized in the AAPM Task Group 43 dose calculation formalism, were calculated using Monte Carlo simulation. Results are compared with corresponding data published for the existing VariSource and microSelectron high dose rate sources. The dose-rate constant for the new Varian source was found to be equal to 1.101 +/- 0.006cGyh(-1) U(-1), compared to values of 1.043 +/- 0.005 and 1.116 +/- 0.006 cGyh(-1) U(-1) calculated for the existing VariSource and microSelectron sources, respectively. The radial dose functions between the three sources are similar with the exception of their values at radial distances very close to the source (r approximately 2 mm) where differences of approximately 3% are observed. The new Varian source demonstrates a smaller anisotropy relative to the existing VariSource source design for polar angles close to the source longitudinal axis, due to its smaller active core length.

Dosimetry of 192Ir wires for LDR interstitial brachytherapy following the AAPM TG-43 dosimetric formalism.

Implementation of the AAPM Task Group 43 dosimetric formalism for 192Ir wires used as interstitial sources in low dose-rate (LDR) brachytherapy applications is investigated. Geometry factors, dose-rate constant values, radial dose functions, and anisotropy functions, utilized in this formalism, were calculated for various lengths of all commercially available wire source designs by means of a well-established Monte Carlo simulation code and an improved modification of the Sievert integral method. Results are presented in the form of look up tables that allow interpolation for dose-rate calculations around all practically used wire lengths, with accuracy acceptable for clinical applications.

Dosimetry comparison of 192Ir sources.

192Ir sources besides being widely utilized in the field of conventional brachytherapy also find use in contemporary peripheral and coronal intravascular applications. In this study, the same Monte Carlo simulation code and input data were used to investigate differences between the dose rate distributions of the most commonly used 192Ir sources in the cm and mm distance range. Findings are discussed in view of differences in source and encapsulation dimensions as well as structural details. Results are presented in the AAPM TG-43 formalism, as generalized by AAPM TG-60, for five 192Ir HDR source designs as well as an LDR seed and an LDR wire source. Dose rate constants of the sources at r0 = 1 cm and r0 = 2 mm were found proportional to the corresponding geometry factors along the transverse source bisectors and an equation of the form lambda r0(cGyh(-1) U(-1)) = 1.12 x G(r0,90 degrees) provides results within clinical accuracy (less than 2%) for any 192Ir source. Radial dose functions do not depend significantly on source and encapsulation geometry and agree within 2% with that of a point 192Ir source. Anisotropy is of importance for accurate dosimetry at the cm distance range but it does not affect dose rate in the mm distance range significantly. At such short radial distances the source geometry factor defines the shape of isodose lines. Dose uniformity at given distances from the sources is strongly dependent on source dimensions as indicated by dose rate profiles in polar and Cartesian coordinates.

Thermoluminescent dosimetry of the selectseed 125I interstitial brachytherapy seed.

This work presents experimental dosimetry results for the new selectSeed 125I prostate seed design for use with the seedSelectron afterloading device, in accordance with the AAPM advisory that all new low energy interstitial brachytherapy seeds should undergo one Monte Carlo (MC) and at least one experimental dosimetry characterization. TLD dosimetry was performed using 120 cylindrical LiF TLD type-100 rods calibrated using a 6 MV photon beam. They were irradiated in solid water phantoms for the experimental determination of the seed dose rate constant, radial dose functions and anisotropy functions. MC simulations were performed for the determination of the TLDs relative energy response that was found position independent and equal to 1.40+/-0.03, and for the calculation of the ratio of dose in liquid water to dose in solid water that was found to be well described by Dliquidwater/Dsolidwater= 1.013*r+0.030 presenting only a minor dependence on polar angle. The selectSeed dose rate constant in liquid water was found equal to 0.938+/-0.065 cGy h(-1) U(-1), which agrees within experimental uncertainties with corresponding MC results of lambdaselect Seed=0.954+/-0.005 cGy h(-1) U(-1). The experimental radial dose and anisotropy function results were also found in good agreement with corresponding MC calculations.