Design and characterization of a new high-dose-rate brachytherapy Valencia applicator for larger skin lesions.

PURPOSE: The aims of this study were (i) to design a new high-dose-rate (HDR) brachytherapy applicator for treating surface lesions with planning target volumes larger than 3 cm in diameter and up to 5 cm in size, using the microSelectron-HDR or Flexitron afterloader (Elekta Brachytherapy) with a (192)Ir source; (ii) to calculate by means of the Monte Carlo (MC) method the dose distribution for the new applicator when it is placed against a water phantom; and (iii) to validate experimentally the dose distributions in water. METHODS: The penelope2008 MC code was used to optimize dwell positions and dwell times. Next, the dose distribution in a water phantom and the leakage dose distribution around the applicator were calculated. Finally, MC data were validated experimentally for a (192)Ir mHDR-v2 source by measuring (i) dose distributions with radiochromic EBT3 films (ISP); (ii) percentage depth-dose (PDD) curve with the parallel-plate ionization chamber Advanced Markus (PTW); and (iii) absolute dose rate with EBT3 films and the PinPoint T31016 (PTW) ionization chamber. RESULTS: The new applicator is made of tungsten alloy (Densimet) and consists of a set of interchangeable collimators. Three catheters are used to allocate the source at prefixed dwell positions with preset weights to produce a homogenous dose distribution at the typical prescription depth of 3 mm in water. The same plan is used for all available collimators. PDD, absolute dose rate per unit of air kerma strength, and off-axis profiles in a cylindrical water phantom are reported. These data can be used for treatment planning. Leakage around the applicator was also scored. The dose distributions, PDD, and absolute dose rate calculated agree within experimental uncertainties with the doses measured: differences of MC data with chamber measurements are up to 0.8% and with radiochromic films are up to 3.5%. CONCLUSIONS: The new applicator and the dosimetric data provided here will be a valuable tool in clinical practice, making treatment of large skin lesions simpler, faster, and safer. Also the dose to surrounding healthy tissues is minimal.

Design and evaluation of a HDR skin applicator with flattening filter.

The purposes of this study are: (i) to design field flattening filters for the Leipzig applicators of 2 and 3 cm of inner diameter with the source traveling parallel to the applicator contact surface, which are accessories of the microSelectron-HDR afterloader (Nucletron, Veenendaal, The Netherlands). These filters, made of tungsten, aim to flatten the heterogeneous dose distribution obtained with the Leipzig applicators. (ii) To estimate the dose rate distributions for these Leipzig+filter applicators by means of the Monte Carlo (MC) method. (iii) To experimentally verify these distributions for prototypes of these new applicators, and (iv) to obtain the correspondence factors to measure the output of the applicators by the user using an insert into a well chamber. The MC GEANT4 code has been used to design the filters and to obtain the dose rate distributions in liquid water for the two Leipzig+filter applicators. In order to validate this specific application and to guarantee that realistic source-applicator geometry has been considered, an experimental verification procedure was implemented in this study, in accordance with the updated recommendations of the American Association of Physicists in Medicine Task Group No. 43 U1 Report. Thermoluminescent dosimeters, radiochromic film, and a pin-point ionization chamber in a plastic [polymethylmethacrylate (PMMA)] phantom were used to verify the MC results for the two applicators of a microSelectron-HDR afterloader with the mHDR-v2 source. To verify the output of the Leipzig +filter applicators, correspondence factors were deduced for the well chambers HDR100-plus (Standard Imaging, Inc., Middleton, WI) and TM33004 (PTW, Freiburg, Germany) using a specific insert for both applicators. The doses measured in the PMMA phantom agree within experimental uncertainties with the dose obtained by the MC calculations. Percentage depth dose and off-axis profiles were obtained normalized at a depth of 3 mm along the central applicator axis in a cylindrical 20 x 20 cm water phantom. A table of output factors, normalized to 1 U of source air kerma strength at this depth, is presented. Correspondence factors were obtained for the two well chambers considered. The matrix data obtained in the MC simulation with a grid separation of 0.5 mm has been used to build a data set in a convenient format to model these distributions for routine use with a brachytherapy treatment planning system.

Dosimetric characterization of Ir-192 LDR elongated sources.

Ir-192 wires have been used in low-dose-rate brachytherapy for many years. Commercially available treatment planning systems approximate the dose rate distribution of the straight or curved wires applying the superposition principle using one of the following methods: (i) The wire is modeled as a set of point sources, (ii) the wire is modeled as a set of small straight segment wires, (iii) the values of the parameters and functions of the American Association of Physicists in Medicine (AAPM) Task Group 43 protocol are obtained for wire lengths between 3 and 7 cm assuming some simplifications. The dose rate distributions obtained using these methods for linear wires of different lengths and U-shaped wires present significant deviations compared to those obtained by Monte Carlo. In the present study we propose a new method to model 192Ir wires of any length and shape, named the Two Lengths based Segmented method. This method uses the formalism stated in the AAPM Task Group 43 protocol for two straight wires only, 0.5 and 1 cm, to obtain the dose rate distribution around wires of any length (down to 0.3 cm and up to 10 cm) improving on the results of the aforementioned ones. This method can easily be applied to dose calculations around other wires, such as Pd-103 ones.

Technique for routine output verification of Leipzig applicators with a well chamber.

The H-type Leipzig applicators are accessories of the microSelectron-HDR system (Nucletron, Veenendaal, The Netherlands) for treatment of superficial malignancies. Recently, the dose rate distributions in liquid water for the whole set of applicators using both source models available for the microSelectron-HDR afterloaders have been obtained by means of the experimentally validated Monte Carlo (MC) code GEANT4. Also an output table (cGy/hU) at 3 mm depth on the applicator central axis was provided. The output verification of these applicators by the user, prior to their clinical use, present practical problems: small detectors such as thermoluminescent dosimeters or parallel-plate ionization chambers are not easily used for verification in a clinical environment as they require a rigid setup with the Leipzig applicator and a phantom. In contrast, well-type ionization chambers are readily available in radiotherapy departments. This study presents a technique based on the HDR1000Plus well chamber (Standar Imaging) measurements with a special insert, which allows the output verification of the H-type Leipzig applicators on a routine basis. This technique defines correspondence factors (CF) between the in water dose rate output of the Leipzig applicators (cGy/hU) obtained with MC and the reading on the well chamber with the special insert, normalized to the HDR calibration factor with the HDR insert and to the source strength. To commission the applicators (with the well chamber and the special insert used), the physicist should check if the CF value agrees with its tabulated values presented in this work. If the differences are within 5% the tabulated output values can be used in clinical dosimetry. This technique allows the output validation of the Leipzig applicators with a well chamber widely used for HDR Ir-192 source strength measurements. It can easily be adapted to other types of well chambers for HDR source output verification.

The combined use of the natural and the cumulative dose-volume histograms in planning and evaluation of permanent prostatic seed implants.

BACKGROUND AND PURPOSE: To investigate prostate dose coverage and overdosage in planned and realized permanent iodine seed prostate implants and to explore the use of the natural dose-volume histogram (NDVH) and the cumulative dose-volume histogram (CDVH) as tools to optimize prostate implants. MATERIALS AND METHODS: The optimal prescription dose (PD) or natural prescription dose (NPD) was derived from the NDVH. The mismatch between the NPD and the given PD was called the natural dose ratio (NDR). For an ideal implant the NDR should be 1. The target is overdosed if NDR >1 and underdosed if NDR <1. The NDR and prostate coverage were evaluated in implants of nine patients. Prostate coverage was determined from the CDVH based on pre-implant ultrasound or post-implant MRI for the planned and realized implants, respectively. The use of the NDVH to further optimize the planned prostate implants was also explored. RESULTS: The mean values of the NDRs were 1.30+/-0.34 (range 0.76-1.79), 1.22+/-0.31 (0.76-1.74) and 1.22+/-0.12 (0.98-1.33) for the planned, realized and optimized seed distributions, respectively. The realized prostatic implants showed smaller prostate coverage than the planned implants. The prostate volume fractions receiving 100% of the prescription dose were V(100)=79+/-6% and V(100)=97+/-3% for the realized and the planned implants, respectively. CONCLUSIONS: The NDVH and the CDVH proved to be valuable tools in plan evaluation. The NDVH and its derived parameter NDR quantify the risk of under or overdosage for a given PD. The CDVH is valuable in evaluation of prostate coverage realized prostate. Our strategy to implant just the prostate and not the prostate plus a margin led to NDR values between 1.1 and 1.3 and a prostate coverage of V(100)=79+/-6% in the nine patients. The planned coverage of V(100)=95% was not realized, mainly due to inadequate coverage of the base of the prostate.

Dependence of the tray transmission factor on collimator setting and source-surface distance.

When blocks are placed on a tray in megavoltage x-ray beams, generally a single correction factor for the attenuation by the tray is applied for each photon beam quality. In this approach, the tray transmission factor is assumed to be independent of field size and source-surface distance (SSD). Analysis of a set of measurements performed in beams of 13 different linear accelerators demonstrates that there is, however, a slight variation of the tray transmission factor with field size and SSD. The tray factor changes about 1.5% for collimator settings varying between 4x4 cm and 40 x 40 cm for a 1 cm thick PMMA tray and approximately 3% for a 2 cm thick PMMA tray. The variation with field size is smaller if the source-surface distance is increased. The dependence on the collimator setting is not different, within the experimental uncertainty of about 0.5% (1 s.d.), for the nominal accelerating potentials and accelerator types applied in this study. It is shown that the variation of the tray transmission factor with field size and source-surface distance can easily be taken into account in the dose calculation by considering the volume of the irradiated tray material and the position of the tray in the beam. A relation is presented which can be used to calculate the numerical value of the tray transmission factor directly. These calculated values can be checked with only a few measurements using a cylindrical beam coaxial miniphantom.

Effect of electron contamination on scatter correction factors for photon beam dosimetry.

Physical quantities for use in megavoltage photon beam dose calculations which are defined at the depth of maximum absorbed dose are sensitive to electron contamination and are difficult to measure and to calculate. Recently, formalisms have therefore been presented to assess the dose using collimator and phantom scatter correction factors, Sc and Sp, defined at a reference depth of 10 cm. The data can be obtained from measurements at that depth in a miniphantom and in a full scatter phantom. Equations are presented that show the relation between these quantities and corresponding quantities obtained from measurements at the depth of the dose maximum. It is shown that conversion of Sc and Sp determined at a 10 cm depth to quantities defined at the dose maximum such as (normalized) peak scatter factor, (normalized) tissue-air ratio, and vice versa is not possible without quantitative knowledge of the electron contamination. The difference in Sc at dmax resulting from this electron contamination compared with Sc values obtained at a depth of 10 cm in a miniphantom has been determined as a multiplication factor, Scel, for a number of photon beams of different accelerator types. It is shown that Scel may vary up to 5%. Because in the new formalisms output factors are defined at a reference depth of 10 cm, they do not require Scel data. The use of Sc and Sp values, defined at a 10 cm depth, combined with relative depth-dose data or tissue-phantom ratios is therefore recommended. For a transition period the use of the equations provided in this article and Scel data might be required, for instance, if treatment planning systems apply Sc data normalized at d(max).

A consistent formalism for the application of phantom and collimator scatter factors.

A coherent system for the use of scatter correction factors, determined at 10 cm depth, is described for dose calculations on the central axis of arbitrarily shaped photon beams. The system is suitable for application in both the fixed source-surface distance (SSD) and in the isocentric treatment set-up. This is in contrast to some other proposals where only one of these approaches forms the basis of the calculation system or where distinct quantities and data sets are needed. In order to derive the relations in the formalism, we introduced a separation of the phenomena related to the energy fluence in air and to the phantom scatter contribution to the dose. Both are used relative to quantities defined for the reference irradiation set-up. It is shown that dose calculations can be performed with only one set of basic beam data, obtained at a reference depth of 10 cm. These data consist for each photon beam quality of measured collimator and phantom scatter correction factors, in combination with a set of (percentage/relative) depth-dose or tissue-phantom ratio values measured along the central axis of the beam. Problems related to measurements performed at the depth of maximum absorbed dose, due to the electron contamination of the beam, are avoided in this way. Collimator scatter correction factors are obtained by using a mini-phantom, while phantom scatter correction factors are derived from measurements in a full scatter phantom in combination with the results of the mini-phantom measurements. For practical reasons the fixed SSD system was chosen to determine the data. Then, dose calculations in a fixed SSD treatment set-up itself are straightforward. Application in the isocentric treatment set-up needs simple conversion steps, while the inverse approach, from isocentric to fixed SSD, is described as well. Differences between the two approaches are discussed and the equations for the conversions are given.

A three-Gaussian fit of phantom scatter correction data.

The phantom scatter correction factor Sp of megavoltage photon beams can be accurately described using a three-Gaussian fit. The model leads to six parameters, with which Sp(r) is described as a smooth function of the field radius r for beam qualities in the range from 60Co up to 25 MV. The parameters allow Sp values to be calculated at intermediate beam energies and for any field shape. Calculated Sp(X, Y) values for rectangular fields (X, Y) can be subsequently used as reference values to compare with measured Sp(X, Y) values, for example when appraising a new beam.

CT-based software for 3-D localization and reconstruction in stepping source brachytherapy.

This paper describes innovative software for catheter localization and three-dimensional (3-D) reconstruction in stepping source brachytherapy applications. Patient information is a set of computed tomography (CT) slices scanned during the implantation of brachytherapy catheters. Catheter geometry and patient anatomy are exported for use with dose calculation software modules. The errors produced by the system are also encouragingly low. Time saving was achieved, in terms of other traditional reconstruction techniques. Various automated procedures, 3-D graphics and a user-friendly GUI, have contributed to providing a powerful, comprehensive software tool, directly useable in the clinical practice.

Comparison of parametrization methods of the collimator scatter correction factor for open rectangular fields of 6-25 MV photon beams.

PURPOSE: To facilitate the use of the collimator scatter correction factor, Sc, parametrization methods that relate Sc to the field size by fitting were investigated. MATERIALS AND METHODS: Sc was measured with a mini-phantom for five types of dual photon energy accelerators with energies varying between 6 and 25 MV. Using these Sc-data six methods of parametrizing Sc for square fields were compared, including a third-order polynomial of the natural logarithm of the field size normalized to the field size of 10 cm2. Also five methods of determining Sc for rectangular fields were considered, including one which determines the equivalent field size by extending Sterling's method. RESULTS: The deviations between measured and calculated Sc-values were determined for all photon beams and methods investigated in this study. The resulting deviations of the most accurate method varied between 0.07 and 0.42% for square fields and between 0.26 and 0.79% for rectangular fields. A recommendation is given as to how to limit the number of fields for which Sc should be measured in order to be able to accurately predict it for an arbitrary field size.

Is there a need for a revised table of equivalent square fields for the determination of phantom scatter correction factors?

The use of the British Journal of Radiology (BJR) (supplement 17) tables of equivalent square fields for dose calculations is widespread. A revised version of the supplement was published recently, with a more elaborate discussion, but without changes in data given in these tables (Br. J. Radiol. suppl 25). The tables were generated for use in dose calculations, with relative beam data such as PDD, BSF, PSF, all with d(max) as the reference depth. However, the current philosophy in dose calculational methods is based on quantities defined at a reference depth, d(ref) = 10 cm, on a separation of phantom and head scatter, and on the use of the relative depth-dose or tissue-phantom ratios normalized at d(ref). By using these quantities as a starting point, problems at shallow depths related to the influence of contaminating electrons in the beam can be eliminated. Recently, a comprehensive set of phantom scatter factor data with d(ref) = 10 cm has been published for a set of square field sizes and a wide range of photon beam energies, showing that phantom scatter is a smoothly varying function of field size and quality index. It is not a priori evident that the conventional concept of equivalent squares for rectangular fields is also fully applicable for phantom scatter factors and phantom scatter related quantities at a depth of 10 cm. It was questioned whether or not new tables of equivalent square fields are needed for this purpose. In this paper, new tables have been constructed for four photon beam energies in the range of Co-60 to 25 MV (quality index from 0.572 to 0.783). The small differences between the outcome of these new tables allowed the construction of one averaged table of equivalent square fields. Phantom scatter factors were calculated for rectangular fields based on the use of the BJR table and on the use of the newly constructed tables and the differences were quantified. For Co-60 no improvements could be shown when using the new averaged table, but for beam energies of 6 to 10 MV small improvements of the order of 0.5 to 1.0% were found. For a higher beam energy of 25 MV the improvement is smaller. Deviations resulting from the BJR table are within the limits of accuracy as stated by the authors. Therefore, for clinical use, the continued use of the BJR table of equivalent squares for phantom scatter factors and phantom scatter related quantities of rectangular fields is justified, irrespective of photon beam energy.

The determination of phantom and collimator scatter components of the output of megavoltage photon beams: measurement of the collimator scatter part with a beam-coaxial narrow cylindrical phantom.

The separation of the total scatter correction factor Sc,p in a collimator scatter component, Sc, and a phantom scatter component, Sp, has proven to be an useful concept in megavoltage photon beam dose calculations in situations which differ from the standard treatment geometry. A clinically applicable method to determine Sc is described. Measurements are carried out with an ionization chamber, placed at a depth beyond the range of contaminant electrons, in a narrow cylindrical polystyrene phantom with a diameter of 4 cm of which the axis coincides with the beam axis. Sc,p is measured in a full-scatter phantom and Sp can be derived from Sc,p and Sc. In order to obtain a reliable separation, i.e. excluding the influence of contaminant electrons, measurements of Sc,p have been carried out at depths of 5 cm for photon beams with a quality index (QI) up to and including 0.75 and a depth of 10 cm with QI larger than 0.75. These depths are in accordance with recommendations given in recent dosimetry protocols. The consistency of the method was checked by comparing calculated and measured values of Sc,p for a set of blocked fields for a range of photon beam energies from 60Co up to 25 MV showing a maximum deviation of 2%. The method can easily be implemented in existing procedures for the calculation of the number of monitor units to deliver a specified dose to a target volume.

Discrepancies in dose and dose specification in interstitial implants.

The first step in the execution of an interstitial implant is the decision on size and location of the target volume. Several implant systems, e.g. the Paterson-Parker system and the Paris system, give instructions for the optimal arrangement of sources to assure that the planned target volume is adequately covered. They also give guidelines to calculate the reference dose rate encompassing the planned target volume. These systems provide different solutions for the source arrangement for the same planned target volume, and vice versa, resulting in different reference dose rates. The problem of dose specification is discussed. For a number of theoretical implants predicted reference dose rates for the planned target volume were compared with the computer calculated dose rates for that volume. Discrepancies increase when moderate digressions from the adopted implant system rules are allowed, such as could commonly occur clinically. For a number of examples the degree of change in dose rate, if over 10%, and the position where this deviation is likely to occur are described. For optimal results the clinician should be well aware of these variations.

[Standardized low-dose-rate afterloading therapy of cervical cancer using a ring applicator]. / Zur Frage einer standardisierten Low-dose-rate Afterloadingtherapie des Kollumkarzinoms unter Verwendung eines Ringapplikators.

At our hospital, the Munich method has already proved its value in gynecologic brachytherapy. In order to use this method for the afterloading technique (LDR), too, we developed a computer program for a standardized therapy with the circular applicator. A certain arrangement of sources is indicated for all applicator sizes. So the dose prescribed (point A) as well as the treatment times and the maximum rectum doses can be given in tabular form. Furthermore, another reduction of the exposure of the rectum could be achieved by means of a certain symmetric arrangement of the sources. To our opinion, this standardized treatment scheme allows in most cases a sufficiently exact definition of the dose distribution in the pelvis minor. If necessary, a more individualized irradiation planning can be established later using CT scans.