Validation of a precision radiochromic film dosimetry system for quantitative two-dimensional imaging of acute exposure dose distributions

We present an evaluation of the precision and accuracy of image-based radiochromic film (RCF) dosimetry performed using a commercial RCF product (Gafchromic MD-55-2, Nuclear Associates, Inc.) and a commercial high-spatial resolution (100 microm pixel size) He-Ne scanning-laser film-digitizer (Personal Densitometer, Molecular Dynamics, Inc.) as an optical density (OD) imaging system. The precision and accuracy of this dosimetry system are evaluated by performing RCF imaging dosimetry in well characterized conformal external beam and brachytherapy high dose-rate (HDR) radiation fields. Benchmarking of image-based RCF dosimetry is necessary due to many potential errors inherent to RCF dosimetry including: a temperature-dependent time evolution of RCF dose response; nonuniform response of RCF; and optical-polarization artifacts. In addition, laser-densitometer imaging artifacts can produce systematic OD measurement errors as large as 35% in the presence of high OD gradients. We present a RCF exposure and readout protocol that was developed for the accurate dosimetry of high dose rate (HDR) radiation sources. This protocol follows and expands upon the guidelines set forth by the American Association of Physicists in Medicine (AAPM) Task Group 55 report. Particular attention is focused on the OD imaging system, a scanning-laser film digitizer, modified to eliminate OD artifacts that were not addressed in the AAPM Task Group 55 report. RCF precision using this technique was evaluated with films given uniform 6 MV x-ray doses between 1 and 200 Gy. RCF absolute dose accuracy using this technique was evaluated by comparing RCF measurements to small volume ionization chamber measurements for conformal external-beam sources and an experimentally validated Monte Carlo photon-transport simulation code for a 192Ir brachytherapy source. Pixel-to-pixel standard deviations of uniformly irradiated films were less than 1% for doses between 10 and 150 Gy; between 1% and 5% for lower doses down to 1 Gy and 1% and 1.5% for higher doses up to 200 Gy. Pixel averaging to form 200-800 microm pixels reduces these standard deviations by a factor of 2 to 5. Comparisons of absolute dose show agreement within 1.5%-4% of dose benchmarks, consistent with a highly accurate dosimeter limited by its observed precision and the precision of the dose standards to which it is compared. These results provide a comprehensive benchmarking of RCF, enabling its use in the commissioning of novel HDR therapy sources.     

Measurements and Monte Carlo calculations to determine the absolute detector response of radiochromic film for brachytherapy dosimetry

GafChromic (MD-55-2) radiochromic film has become increasingly popular for medical applications and has proven to be useful for brachytherapy dosimetry. To measure the absolute dose near a brachytherapy source, the response of the proposed detector in the measurement conditions relative to the response of the detector in calibration conditions must be known. MD-55-2 radiochromic film has been exposed in four different photon beams, a 30 and 40 kVp tungsten anode x-ray beam, a 75 kVp orthovoltage therapy beam, and a 60Co teletherapy beam to measure the relative detector response. These measurements were combined with coupled photon/electron Monte Carlo transport calculations to determine the absolute detector response. The Los Alamos National Laboratory Monte Carlo transport code MCNP4B2 was used. The measured relative response of this batch of MD-55-2 film varies from 8.79 mOD/Gy, measured for the 60Co beam, by as much as 42% for the low-energy x-ray beams. However, the absolute detector response varies from 4.32 mOD/Gy for the 60Co beam by, at most, only 6.3%. In this work we demonstrate that the absolute detector response of MD-55-2 radiochromic film is a constant and independent of beam quality. Further, this work shows that MCNP4B2 accurately simulates the energy response and geometry artifacts of the radiochromic film.     

Experimental validation of dose calculation algorithms for the GliaSite RTS, a novel 125I liquid-filled balloon brachytherapy applicator

This paper compares experimentally measured and calculated dose-rate distributions for a novel 125I liquid-filled brachytherapy balloon applicator (the GliaSite RTS), designed for the treatment of malignant brain-tumor resection-cavity margins. This work is intended to comply with the American Association of Physicists in Medicine (AAPM) Radiation Therapy Committee's recommendations [Med. Phys. 25, 2269-2270 (1998)] for dosimetric characterization of low-energy photon interstitial brachytherapy sources. Absolute low dose-rate radiochromic film (RCF) dosimetry measurements were performed in coronal planes about the applicator. The applicator was placed in a solid water phantom, machined to conform to the inflated applicator's surface. The results were used to validate the accuracy of Monte Carlo photon transport (MCPT) simulations and a point-source dose-kernel algorithm in predicting dose to water. The absolute activity of the 125I solution was determined by intercomparing a National Institute of Standards and Technology (NIST) 125I standard with a known mass of radiotherapy solution (Iotrex) in an identical vial and geometry. For the two films not in contact with applicator, the average agreement between RCF and MCPT (specified as the mean absolute deviation in successive 4 mm rings) was found to be within +/-5% at distances 0.2-25 mm from the film centers. For the two films touching the catheter, the mean agreement was +/-14.5% and 7.5% near the balloon surface but improving to 7.5% and 6% by 3.5 mm from the surface. These errors, as large as 20% in isolated pixels, are likely due to trim damage, 125I contamination, and poor conformance with the balloon. At larger distances where the radiation doses were very low, the observed discrepancies were significantly larger than expected. We hypothesize that they are due to a dose-rate dependence of the RCF response. A 1%-10% average difference between a simple one-dimensional path-length semiempirical dose-kernel model and the MCPT calculations was observed over clinically relevant distances.     

Dose rate table for a 32P intravascular brachytherapy source from Monte Carlo calculations

Studies of intravascular brachytherapy to prevent restenosis following angioplasty have shown many promising results. Accurate dose rate tables based on detailed models of the brachytherapy sources are necessary for treatment planning. This work will present an away and along dose rate table for a 27 mm long catheter based 32P beta source. MD-55-2 radiochromic film has been exposed at five different depths (0.5 mm-4 mm) in a polystyrene phantom using a 27 mm long Guidant 32P beta source. The total dose to the active region of the film was determined using the absolute detector response of the MD-55-2 radiochromic film. The Monte Carlo code MCNP4B2 was also used to calculate the dose to the active region of the film using a detailed model of the source, encapsulation, and radiochromic film. The dose to film calculations showed good agreement with the measurements presented in this work with an average difference of 7%. The Monte Carlo calculations were also verified against previously published depth dose in water measurements determined using radiochromic film and plastic scintillator. The depth dose calculations in water showed good agreement with the previously published measurements with the calculations being about 2.5% lower than the film measurements and about 2.5% higher than the scintillator measurements. This work then uses the verified Monte Carlo code to present a dose rate table for the 32P intravascular beta source.     

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 of 192Ir sources used for endovascular brachytherapy

An in-phantom calibration technique for 192Ir sources used for endovascular brachytherapy is presented. Three different source lengths were investigated. The calibration was performed in a solid phantom using a Farmer-type ionization chamber at source to detector distances ranging from 1 cm to 5 cm. The dosimetry protocol for medium-energy x-rays extended with a volume-averaging correction factor was used to convert the chamber reading to dose to water. The air kerma strength of the sources was determined as well. EGS4 Monte Carlo calculations were performed to determine the depth dose distribution at distances ranging from 0.6 mm to 10 cm from the source centre. In this way we were able to convert the absolute dose rate at 1 cm distance to the reference point chosen at 2 mm distance. The Monte Carlo results were confirmed by radiochromic film measurements, performed with a double-exposure technique. The dwell times to deliver a dose of 14 Gy at the reference point were determined and compared with results given by the source supplier (CORDIS). They determined the dwell times from a Sievert integration technique based on the source activity. The results from both methods agreed to within 2% for the 12 sources that were evaluated. A Visual Basic routine that superimposes dose distributions, based on the Monte Carlo calculations and the in-phantom calibration, onto intravascular ultrasound images is presented. This routine can be used as an online treatment planning program.     

Development of an adjoint sensitivity field-based treatment-planning technique for the use of newly designed directional LDR sources in brachytherapy

The development and application of an automated 3D greedy heuristic (GH) optimization algorithm utilizing the adjoint sensitivity fields for treatment planning to assess the advantage of directional interstitial prostate brachytherapy is presented. Directional and isotropic dose kernels generated using Monte Carlo simulations based on Best Industries model 2301 I-125 source are utilized for treatment planning. The newly developed GH algorithm is employed for optimization of the treatment plans for seven interstitial prostate brachytherapy cases using mixed sources (directional brachytherapy) and using only isotropic sources (conventional brachytherapy). All treatment plans resulted in V100 > 98% and D90 > 45 Gy for the target prostate region. For the urethra region, the D10(Ur), D90(Ur) and V150(Ur) and for the rectum region the V100cc, D2cc, D90(Re) and V90(Re) all are reduced significantly when mixed sources brachytherapy is used employing directional sources. The simulations demonstrated that the use of directional sources in the low dose-rate (LDR) brachytherapy of the prostate clearly benefits in sparing the urethra and the rectum sensitive structures from overdose. The time taken for a conventional treatment plan is less than three seconds, while the time taken for a mixed source treatment plan is less than nine seconds, as tested on an Intel Core2 Duo 2.2 GHz processor with 1GB RAM. The new 3D GH algorithm is successful in generating a feasible LDR brachytherapy treatment planning solution with an extra degree of freedom, i.e. directionality in very little time.     

Photon spectrometry for the determination of the dose-rate constant of low-energy photon-emitting brachytherapy sources

Accurate determination of dose-rate constant (lambda) for interstitial brachytherapy sources emitting low-energy photons (< 50 keV) has remained a challenge in radiation dosimetry because of the lack of a suitable absolute dosimeter for accurate measurement of the dose rates near these sources. Indeed, a consensus value of lambda taken as the arithmetic mean of the dose-rate constants determined by different research groups and dosimetry techniques has to be used at present for each source model in order to minimize the uncertainties associated with individual determinations of lambda. Because the dosimetric properties of a source are fundamentally determined by the characteristics of the photons emitted by the source, a new technique based on photon spectrometry was developed in this work for the determination of dose-rate constant. The photon spectrometry technique utilized a high-resolution gamma-ray spectrometer to measure source-specific photon characteristics emitted by the low-energy sources and determine their dose-rate constants based on the measured photon-energy spectra and known dose-deposition properties of mono-energetic photons in water. This technique eliminates many of the difficulties arising from detector size, the energy dependence of detector sensitivity, and the use of non-water-equivalent solid phantoms in absolute dose rate measurements. It also circumvents the uncertainties that might be associated with the source modeling in Monte Carlo simulation techniques. It was shown that the estimated overall uncertainty of the photon spectrometry technique was less than 4%, which is significantly smaller than the reported 8-10% uncertainty associated with the current thermo-luminescent dosimetry technique. In addition, the photon spectrometry technique was found to be stable and quick in lambda determination after initial setup and calibration. A dose-rate constant can be determined in less than two hours for each source. These features make it ideal to determine the dose-rate constant of each source model from a larger and more representative sample of actual sources and to use it as a quality assurance resource for periodic monitoring of the constancy of lambda for brachytherapy sources used in patient treatments.     

Comparison of measured and calculated dose rates in water near I-125 and Ir-192 seeds

Recent theoretical and experimental work indicates that currently accepted 125I dosimetry data may overestimate dose in water at 1 cm by 10%-24%. Among the most comprehensive measurements are those of the NCI-sponsored brachytherapy contract participants. Absolute dose rates in water calculated by the Monte Carlo method have been compared with the NCI dose measurements about 125I and 192Ir seeds embedded in solid-water phantoms. The photon transport code allows realistic geometric simulation of the complex internal seed structure, the National Institute of Standards and Technology air-kerma strength standardization geometry, and the dose measurement setup. When the appropriate measurement medium and geometry are assumed, agreement between theory and measurement is excellent, within 3% at 1 cm and averaging 3% at larger distances. However, the data do not support the water equivalence of solid water at 125I energies indicating that solid-water measurements underestimate 125I specific dose-rate constants in water by 4.3%. Because of its higher ratio of absorption to scatter, 125I dose distributions measured in solid water are less penetrating (by 35% at 10 cm) than those measured in liquid water. For model 6711, model 6702, and steel-clad 192Ir seeds, Monte Carlo calculations yielded specific dose-rate constants (assuming liquid water medium) of 0.877, 0.932, and 1.122 cGy cm2 h-1 per unit air-kerma strength, respectively. For 125I, currently accepted values are 18% and 11% larger for the two seed models.     

A standard dosimetry procedure for 192Ir sources used for endovascular brachytherapy

The experimental dosimetry of a high dose rate (HDR) 192Ir source used for the brachytherapy of peripheral vessels is reported. The direct determination of the reference air kerma rate Kr agrees, within the experimental uncertainty, with the results obtained by a well ionization chamber calibrated at the NIST and the manufacturer's certification. A highly sensitive (HS) radiochromic film (RCF), that presents only one active layer, was used for the source dosimetry in a water phantom. The adopted experimental set-up, with the source in its catheter positioned on the RCF plane, seems to have given better accuracy of the RCF optical density measurements. The agreement between the measurement of the dose rate constant DKr (10 mm, pi/2) and the literature data confirmed the coherence of the HS RCF calibration obtained by the kerma in air measurements. The RCF measurements supplied dosimetric information about the dose to water per reference air kerma rate D(r, theta)/Kr along the source transverse bisector axis, the radial dose function g(r) and the anisotropy function F(r, theta). The value D(2 mm, pi/2)/Kr = 22.4 +/- 1.2 cGy h(-1)/(microGy h(-1)) is supplied with a dose uncertainty that is essentially due to the indeterminacy of the source position in the catheter. The data of the radial and anisotropy functions have been compared with Monte Carlo determinations reported in the literature.     

Dependence of radiochromic film optical density post-exposure kinetics on dose and dose fractionation

Radiochromic film (RCF) has been shown to be a precise and accurate secondary planar dosimeter for acute exposure radiation fields. However, its application to low dose-rate brachytherapy has been questioned because of possible dose-rate effects. To address this concern, we have measured the optical density (OD) of Model 55-2 RCF as a function of time (interval between the completion of irradiation and densitometry using a 633 nm laser scanner) following exposure (from less than 1 hour to 90 days) for single and split doses from 1 Gy to 100 Gy. Our work demonstrates that film darkening as a function of post-irradiation time depends significantly on total dose, with films exposed to lower doses developing faster than films given higher doses. At 1 Gy, the OD 90 days after exposure is 200% larger than that measured 1 h after exposure compared to a 20% increase over 90 days for doses larger than 20 Gy. An empirical model with time-independent, fast and slow growth terms was used to fit single exposure data. The dependence of the resulting best-fit parameters on dose was investigated. Splitting the dose into two fractions (20 Gy followed by doses of 1-80 Gy 24 h later) results in modest post-irradiation time-dependent changes in the total optical density (at most 15% at small doses), which dissipates within 20 hours following the second exposure. This experimental finding is consistent with the predictions of a simple cumulative dose superposition model. Overall, both experimental and empirical modeling suggest that dose-rate effects may be relatively small despite the strong dependence of film darkening kinetics on total dose. However, more experimental evaluation of radiochromic film response dependence on dose rate and dose-time-fractionation patterns is needed.     

Comparison of Monte Carlo calculations around a Fletcher Suit Delclos ovoid with radiochromic film and normoxic polymer gel dosimetry

The Fletcher Suit Delclos (FSD) ovoids employed in intracavitary brachytherapy (ICB) for cervical cancer contain shields to reduce dose to the bladder and rectum. Many treatment planning systems (TPS) do not include the shields and other ovoid structures in the dose calculation. Instead, TPSs calculate dose by summing the dose contributions from the individual sources and ignoring ovoid structures such as the shields. The goal of this work was to calculate the dose distribution with Monte Carlo around a Selectron FSD ovoid and compare these calculations with radiochromic film (RCF) and normoxic polymer gel dosimetry. Monte Carlo calculations were performed with MCNPX 2.5.c for a single Selectron FSD ovoid with and without shields. RCF measurements were performed in a plane parallel to and displaced laterally 1.25 cm from the long axis of the ovoid. MAGIC gel measurements were performed in a polymethylmethacrylate phantom. RCF and MAGIC gel were irradiated with four 33 microGy m2 h(-1) Cs-137 pellets for a period of 24 h. Results indicated that MCNPX calculated dose to within +/- 2% or 2 mm for 98% of points compared with RCF measurements and to within +/- 3% or 3 mm for 98% of points compared with MAGIC gel measurements. It is concluded that MCNPX 2.5.c can calculate dose accurately in the presence of the ovoid shields, that RCF and MAGIC gel can demonstrate the effect of ovoid shields on the dose distribution and the ovoid shields reduce the dose by as much as 50%.     

Dosimetric properties of the new 125I BrachySeed model LS-1 source

The BrachySeed model LS-1 is one of the latest in a series of new brachytherapy 125I seeds that have recently become available commercially for interstitial implants. The dosimetric properties of the seed were investigated analytically, experimentally, and by Monte Carlo simulation. Following the AAPM Task Group 43 formalism, the radial dose function, dose rate constant, and anisotropy parameters were determined. Experimental measurements were made in solid water-equivalent phantoms using GafChromic MD-55-2 films, with correction for the low energy film response. Analyses were carried out from absolute measurements, as well as relative measurements against the Nycomed Amersham OncoSeed Model 6711, which also served to validate our experimental methodology. A small, but systematic difference in the absolute measurements was observed depending on the duration of the irradiation. Monte Carlo simulation was performed using the Integrated Tiger Series CYLTRAN code. We benchmarked the code by comparing the dose parameters of Model 6702 with published values. The radial dose function, g(r), of the Model LS-1 seed was computed at distances from 0.25 to 10 cm by analytical and Monte Carlo calculations with reasonably good agreement. The suggested dose rate constant, A, based on the Monte Carlo simulation is 0.90+/-0.03 cGy h(-1) U(-1). This value is smaller than, but overlaps the experimental determination of 0.98+/-0.06 cGy h(-1) U(-1). The anisotropy function, F(r, theta), and anisotropy factor, phi(an)(r), compared favorably with those of the Model 6711.     

Monte Carlo dose calculations using MCNP4C and EGSnrc/BEAMnrc codes to study the energy dependence of the radiochromic film response to beta-emitting sources

The energy dependence of the radiochromic film (RCF) response to beta-emitting sources was studied by dose theoretical calculations, employing the MCNP4C and EGSnrc/BEAMnrc Monte Carlo codes. Irradiations with virtual monochromatic electron sources, electron and photon clinical beams, a (32)P intravascular brachytherapy (IVB) source and other beta-emitting radioisotopes ((188)Re, (90)Y, (90)Sr/(90)Y,(32)P) were simulated. The MD-55-2 and HS radiochromic films (RCFs) were considered, in a planar or cylindrical irradiation geometry, with water or polystyrene as the surrounding medium. For virtual monochromatic sources, a monotonic decrease with energy of the dose absorbed to the film, with respect to that absorbed to the surrounding medium, was evidenced. Considering the IVB (32)P source and the MD-55-2 in a cylindrical geometry, the calibration with a 6 MeV electron beam would yield dose underestimations from 14 to 23%, increasing the source-to-film radial distance from 1 to 6 mm. For the planar beta-emitting sources in water, calibrations with photon or electron clinical beams would yield dose underestimations between 5 and 12%. Calibrating the RCF with (90)Sr/(90)Y, the MD-55-2 would yield dose underestimations between 3 and 5% for (32)P and discrepancies within +/-2% for (188)Re and (90)Y, whereas for the HS the dose underestimation would reach 4% with (188)Re and 6% with (32)P.     

Validation of a dose deposited by low-energy photons using GATE/GEANT4

The GATE Monte Carlo simulation platform based on the Geant4 toolkit has now become a diffused tool for simulating PET and SPECT imaging devices. In this paper, we explore its relevance for dosimetry of low-energy 125I photon brachytherapy sources used to treat prostate cancers. To that end, three 125-iodine sources widely used in prostate cancer brachytherapy treatment have been modelled. GATE simulations reproducing dosimetric reference observables such as radial dose function g(r), anisotropy function F(r, theta) and dose-rate constant (Lambda) were performed in liquid water. The calculations were splitted on the EGEE grid infrastructure to reduce the computing time of the simulations. The results were compared to other relevant Monte Carlo results and to measurements published and fixed as recommended values by the AAPM Task Group 43. GATE results agree with consensus values published by AAPM Task Group 43 with an accuracy better than 2%, demonstrating that GATE is a relevant tool for the study of the dose induced by low-energy photons.     

Monte Carlo calculations and measurements of absorbed dose per monitor unit for the treatment of uveal melanoma with proton therapy

The treatment of uveal melanoma with proton radiotherapy has provided excellent clinical outcomes. However, contemporary treatment planning systems use simplistic dose algorithms that limit the accuracy of relative dose distributions. Further, absolute predictions of absorbed dose per monitor unit are not yet available in these systems. The purpose of this study was to determine if Monte Carlo methods could predict dose per monitor unit (D/MU) value at the center of a proton spread-out Bragg peak (SOBP) to within 1% on measured values for a variety of treatment fields relevant to ocular proton therapy. The MCNPX Monte Carlo transport code, in combination with realistic models for the ocular beam delivery apparatus and a water phantom, was used to calculate dose distributions and D/MU values, which were verified by the measurements. Measured proton beam data included central-axis depth dose profiles, relative cross-field profiles and absolute D/MU measurements under several combinations of beam penetration ranges and range-modulation widths. The Monte Carlo method predicted D/MU values that agreed with measurement to within 1% and dose profiles that agreed with measurement to within 3% of peak dose or within 0.5 mm distance-to-agreement. Lastly, a demonstration of the clinical utility of this technique included calculations of dose distributions and D/MU values in a realistic model of the human eye. It is possible to predict D/MU values accurately for clinical relevant range-modulated proton beams for ocular therapy using the Monte Carlo method. It is thus feasible to use the Monte Carlo method as a routine absolute dose algorithm for ocular proton therapy.     

The radial depth-dose distribution of a 188W/188Re beta line source measured with novel,ultra-thin TLDs in a PMMA phantom: comparison with Monte Carlo simulations

The radial depth-dose distribution of a prototype 188W/188Re beta particle line source of known activity has been measured in a PMMA phantom, using a novel, ultra-thin type of LiF:Mg,Cu,P thermoluminescent detector (TLD). The measured radial dose function of this intravascular brachytherapy source agrees well with MCNP4C Monte Carlo simulations, which indicate that 188Re accounts for > or = 99% of the dose between 1 mm and 5 mm radial distance from the source axis. The TLDs were calibrated using a 90Sr/90Y beta secondary standard. Several correction factors are calculated using analytical and Monte Carlo methods. An analysis of the measurement uncertainty is made. Since it is partly determined by components of uncertainty arising from random effects, repeated measurements yield a lower uncertainty. The expanded uncertainty in the absolute dose at 2 mm radial distance equals 11%, 10%, 9% and 8% for 1, 2, 3 and 5 measurements, respectively. After a correction for source non-uniformity, the measured dose rate per unit source activity at 2 mm radial distance equals (1.53 +/- 0.16) Gy min(-1) GBq(-1) (2sigma), in agreement with the value of (1.45 +/- 0.01) Gy min(-1) GBq(-1) (2sigma) predicted by the MCNP4C simulations.     

ALGEBRA: ALgorithm for the heterogeneous dosimetry based on GEANT4 for BRAchytherapy

Task group 43 (TG43)-based dosimetry algorithms are efficient for brachytherapy dose calculation in water. However, human tissues have chemical compositions and densities different than water. Moreover, the mutual shielding effect of seeds on each other (interseed attenuation) is neglected in the TG43-based dosimetry platforms. The scientific community has expressed the need for an accurate dosimetry platform in brachytherapy. The purpose of this paper is to present ALGEBRA, a Monte Carlo platform for dosimetry in brachytherapy which is sufficiently fast and accurate for clinical and research purposes. ALGEBRA is based on the GEANT4 Monte Carlo code and is capable of handling the DICOM RT standard to recreate a virtual model of the treated site. Here, the performance of ALGEBRA is presented for the special case of LDR brachytherapy in permanent prostate and breast seed implants. However, the algorithm is also capable of handling other treatments such as HDR brachytherapy.     

Dosimetric evaluation of the Fletcher-Williamson ovoid for pulsed-dose-rate brachytherapy: a Monte Carlo study

We used radiochromic film dosimetry to validate a Monte Carlo (MC) model of a 192Ir pulsed-dose-rate (PDR) source inside a Fletcher-Williamson ovoid. MD-55-2 radiochromic film was placed in a high-impact polystyrene phantom in a plane parallel to and displaced 2.0 cm medially from the long axis of the ovoid. MC N-particle transport code (MCNPX) version 2.4 was used to model the ovoid and the 192Ir source. Energy deposition was calculated using a track-length estimator modified by an energy-dependent heating function, which is a good approximation of the collision kerma. To convert the estimates of the MC dose per simulated particle to clinically relevant absolute dosimetry, additional MC models of an actual and a virtual 192Ir source in dry air were simulated to determine air kerma strength for the penetrating part of the photon spectrum (>11.3 keV). The absolute dose distributions predicted by MCNPX agreed with the film results and were within +/-9.4% (k = 2) and within +/-2% or within a distance to agreement of 2 mm for 94% of the dose grid. Additional MC models characterized the uncertainty resulting from source positioning inside the ovoid. For a worst-case scenario of 1 mm off centre from the nominal source position in the 3 mm diameter ovoid shaft, the average dose deviation over the film plane was +/-5% (1sigma = +/-4%), with maximum deviation near the sharp dose-gradient provided by the shields of -20% to + 26%. A validated MC model is the first requirement to simulate common LDR clinical loadings (5-20 mgRaEq) and, thus, will aid in the transition from the current 137Cs Selectron LDR ICBT to PDR for treatment of gynecologic cancers.     

On the development of consensus values of reference dosimetry parameters for interstitial brachytherapy sources

The American Association of Physicists in Medicine recommends that the reference dose-rate distribution, used for treatment planning for low-energy photon brachytherapy sources in routine clinical use, must be based on at least two independent determinations: one using experimentally measured dose rates and one using Monte Carlo simulation dosimetry techniques. In this work, we present an approach for developing consensus dosimetry parameters from various independent reference dosimetry studies for interstitial brachytherapy sources. This approach is applied to four recently published papers on the dosimetric properties of the BrachySeed Model LS-1 125I seed. Consensus values for the dose-rate constant, radial dose function, and anisotropy parameters are presented for the LS-1 Model 125I seed.