Free automatic software for quality assurance of computed tomography calibration, edges and radiomics metrics reproducibility

PURPOSE

To develop a QA procedure, easy to use, reproducible and based on open-source code, to automatically evaluate the stability of different metrics extracted from CT images: Hounsfield Unit (HU) calibration, edge characterization metrics (contrast and drop range) and radiomic features.

METHODS

The QA protocol was based on electron density phantom imaging. Home-made open-source Python code was developed for the automatic computation of the metrics and their reproducibility analysis. The impact on reproducibility was evaluated for different radiation therapy protocols, and phantom positions within the field of view and systems, in terms of variability (Shapiro-Wilk test for 15 repeated measurements carried out over three days) and comparability (Bland-Altman analysis and Wilcoxon Rank Sum Test or Kendall Rank Correlation Coefficient).

RESULTS

Regarding intrinsic variability, most metrics followed a normal distribution (88% of HU, 63% of edge parameters and 82% of radiomic features). Regarding comparability, HU and contrast were comparable in all conditions, and drop range only in the same CT scanner and phantom position. The percentages of comparable radiomic features independent of protocol, position and system were 59%, 78% and 54%, respectively. The non-significantly differences in HU calibration curves obtained for two different institutions (7%) translated in comparable Gamma Index G (1 mm, 1%, >99%).

CONCLUSIONS

An automated software to assess the reproducibility of different CT metrics was successfully created and validated. A QA routine proposal is suggested.

Evaluation of dosimetric effects of metallic artifact reduction and tissue assignment on Monte Carlo dose calculations for 125 I prostate implants

PURPOSE

Monte Carlo (MC) simulation studies, aimed at evaluating the magnitude of tissue heterogeneity in ¹²⁵ I prostate permanent seed implant brachytherapy (BT), customarily use clinical post-implant CT images to generate a virtual representation of a realistic patient model (virtual patient model). Metallic artifact reduction (MAR) techniques and tissue assignment schemes (TAS) are implemented on the post-implant CT images to mollify metallic artifacts due to BT seeds and to assign tissue types to the voxels corresponding to the bright seed spots and streaking artifacts, respectively. The objective of this study is to assess the combined influence of MAR and TAS on MC absorbed dose calculations in post-implant CT-based phantoms. The virtual patient models used for ¹²⁵ I prostate implant MC absorbed dose calculations in this study are derived from the CT images of an external radiotherapy prostate patient without BT seeds and prostatic calcifications, thus averting the need to implement MAR and TAS.

METHODS

The geometry of the IsoSeed I25.S17plus source is validated by comparing the MC calculated results of the TG-43 parameters for the line source approximation with the TG-43U1S2 consensus data. Four MC absorbed dose calculations are performed in two virtual patient models using the egs_brachy MC code: (1) TG-43-based D_{w,w-TG 43} , (2) D_w,w-MBDC that accounts for interseed scattering and attenuation (ISA), (3) D_m,m that examines ISA and tissue heterogeneity by scoring absorbed dose in tissue, and (4) D_w,m that unlike D_m,m scores absorbed dose in water. The MC absorbed doses (1) and (2) are simulated in a TG-43 patient phantom derived by assigning the densities of every voxel to 1.00 g cm^-3 (water), whereas MC absorbed doses (3) and (4) are scored in the TG-186 patient phantom generated by mapping the mass density of each voxel to tissue according to a CT calibration curve. The MC absorbed doses calculated in this study are compared with VariSeed v8.0 calculated absorbed doses. To evaluate the dosimetric effect of MAR and TAS, the MC absorbed doses of this work (independent of MAR and TAS) are compared to the MC absorbed doses of different ¹²⁵ I source models from previous studies that were calculated with different MC codes using post-implant CT-based phantoms generated by implementing MAR and TAS on post-implant CT images.

RESULTS

The very good agreement of TG-43 parameters of this study and the published consensus data within 3% validates the geometry of the IsoSeed I25.S17plus source. For the clinical studies, the TG-43-based calculations show a D₉₀ overestimation of more than 4% compared to the more realistic MC methods due to ISA and tissue composition. The results of this work generally show few discrepancies with the post-implant CT-based dosimetry studies with respect to the D₉₀ absorbed dose metric parameter. These discrepancies are mainly Type B uncertainties due to the different ¹²⁵ I source models and MC codes.

CONCLUSIONS

The implementation of MAR and TAS on post-implant CT images have no dosimetric effect on the ¹²⁵ I prostate MC absorbed dose calculation in post-implant CT-based phantoms.

Embedded bleeding detector into a PMMA applicator for electron intraoperative radiotherapy

PURPOSE

The aim of this work is to present a ready to industrialize low-cost and easy-to-install bleeding detector for use in intraoperative electron radiation therapy (IOERT). The detector works in stand-alone mode and is embedded into a translucent polymethylmethacrylate (PMMA) applicator avoiding any contact with the patient, which represent a novelty compared to previous designs. The use of this detector will prevent dose misadministration during irradiation in the event of accumulation of fluids in the applicator.

METHODS

The detector is based on capacitive sensor and wireless power-supply electronics. Both sensor and electronics have been embedded in the applicator, so that any contact with the patient would be avoided. Since access to the tumor can be done through different trajectories, the detector has been calibrated for different tilting angles.

RESULTS

The result of the calibration provides us with a fit curve that allows the interpolation of the results at any angle. Comparison of estimated fluid height vs real height gives an error of 1 mm for tilting angles less than 10° and 2 mm for tilting angles greater than 15°. This accuracy is better than the one required by clinic.

CONCLUSIONS

The performance of the bleeding detector was evaluated in situ. No interference was observed between the detector and the beam. In addition, a user-friendly mobile application has been developed to help the surgical team making decisions before and during irradiation. The measurement provided by the mobile application was stable during the irradiation process.

Phantom development for daily checks in electron intraoperative radiotherapy with a mobile linac

PURPOSE

IORT with mobile linear accelerators is a well-established modality where the dose rate and, therefore, the dose per pulse are very high. The constancy of the dosimetric parameters of the accelerator has to be checked daily. The aim of this work is to develop a phantom with embedded detectors to improve both accuracy and efficiency in the daily test of an IORT linac at the surgery room.

METHODS

The developed phantom is manufactured with transparent polymethyl methacrylate (PMMA), allocating 6 parallel-plate chambers: a central one to evaluate the on-axis beam output, another on-axis one placed at a fixed depth under the previous one to evaluate the energy constancy and four off-axis chambers to evaluate the flatness and symmetry. To analyse the readings a specific application has been developed.

RESULTS

For all chambers and energies, the mean saturation and polarization corrections were smaller than 0.7%. The beam is monitored at different levels of the clinical beam. Output, energy constancy and flatness correlate very well with the correspondent values with the complete applicator. During the first six months of clinical use the beam dosimetric parameters showed excellent stability.

CONCLUSIONS

A phantom has been developed with embedded parallel plate chambers attached to the upper applicator part of an IORT linac. The phantom allows a very efficient setup reducing the time to check the parameters. It provides complete dosimetric information (output, energy and flatness) with just one shot and using ionization chambers with minimum saturation effect, as this highly pulsed beam requires.

Improving bleeding detector features for electron intraoperative radiotherapy

PURPOSE

The aim of this work is to improve the potential bleeding detection during intraoperative radiotherapy with linac polymethyl methacrylate applicators (PMMA), based on one previously developed. The improvements carried out have been focused on: i) minimizing the impact of the detector on the visual through the plastic applicators and ii) avoiding the asymmetry in the detection capability when the applicator is tilted.

METHODS

Simulations have been made to select the geometry that provides a reduced visual impact on the applicator as well as allowing an independent response with the tilting angle of the applicator. A low-noise circuit for signal conditioning has been developed. Measurements have been made on three setups: 10 cm, 7 cm and 4 cm applicator diameters, 0° and 45° tilted.

RESULTS

The detector has a visibility through the applicator greater than 50%. Due to the geometry, optimal detection is ensured regardless of its orientation when the applicator is tilted. It is possible to detect the presence of fluid well below the typical perturbing fluid depth established by the clinic (1-1.5 cm).

CONCLUSIONS

The detector can distinguish the presence of around 0.5 cm of fluid depth while showing a high visual field through the PMMA applicators and providing a measure that does not depend on the detector orientation when the applicator is tilted. The prototype is ready for its industrialization by embedding it into the applicator for clinical use. The detector would have a significant impact on both the quality assurance and the outcome of the treatment.

Review of strategies for MRI based reconstruction of endocavitary and interstitial applicators in brachytherapy of cervical cancer

Brachytherapy plays an essential role in the curative intent management of locally advanced cervical cancer. The introduction of the magnetic resonance (MR) as a preferred image modality and the development of new type of applicators with interstitial components have further improved its benefits. The aim of this work is to review the current status of one important aspect in the cervix cancer brachytherapy procedure, namely catheter reconstruction. MR compatible intracavitary and interstitial applicators are described. Considerations about the use of MR imaging (MRI) regarding appropriate strategies for applicator reconstruction, technical requirements, MR sequences, patient preparation and applicator commissioning are included. It is recommendable to perform the reconstruction process in the same image study employed by the physician for contouring, that is, T2 weighted (T2W) sequences. Nevertheless, a clear identification of the source path inside the catheters and the applicators is a challenge when using exclusively T2W sequences. For the intracavitary component of the implant, sometimes the catheters may be filled with some substance that produces a high intensity signal on MRI. However, this strategy is not feasible for plastic tubes or titanium needles, which, moreover, induce magnetic susceptibility artifacts. In these situations, the use of applicator libraries available in the treatment planning system (TPS) is useful, since they not only include accurate geometrical models of the intracavitary applicators, but also recent developments have made possible the implementation of the interstitial component. Another strategy to improve the reconstruction process is based on the incorporation of MR markers, such as small pellets, to be used as anchor points. Many institutions employ computed tomography (CT) as a supporting image modality. The registration of CT and MR image sets should be carefully performed, and its uncertainty previously assessed. Besides, an important research work is being carried out regarding the use of ultrasound and electromagnetic tracking technologies.

Collision-kerma conversion between dose-to-tissue and dose-to-water by photon energy-fluence corrections in low-energy brachytherapy

The AAPM TG-43 brachytherapy dosimetry formalism, introduced in 1995, has become a standard for brachytherapy dosimetry worldwide; it implicitly assumes that charged-particle equilibrium (CPE) exists for the determination of absorbed dose to water at different locations, except in the vicinity of the source capsule. Subsequent dosimetry developments, based on Monte Carlo calculations or analytical solutions of transport equations, do not rely on the CPE assumption and determine directly the dose to different tissues. At the time of relating dose to tissue and dose to water, or vice versa, it is usually assumed that the photon fluence in water and in tissues are practically identical, so that the absorbed dose in the two media can be related by their ratio of mass energy-absorption coefficients. In this work, an efficient way to correlate absorbed dose to water and absorbed dose to tissue in brachytherapy calculations at clinically relevant distances for low-energy photon emitting seeds is proposed. A correction is introduced that is based on the ratio of the water-to-tissue photon energy-fluences. State-of-the art Monte Carlo calculations are used to score photon fluence differential in energy in water and in various human tissues (muscle, adipose and bone), which in all cases include a realistic modelling of low-energy brachytherapy sources in order to benchmark the formalism proposed. The energy-fluence based corrections given in this work are able to correlate absorbed dose to tissue and absorbed dose to water with an accuracy better than 0.5% in the most critical cases (e.g. bone tissue).

Consensus on 3D treatment planning in gynaecologic brachytherapy of the Radiation Oncology Spanish Society (SEOR) Brachytherapy Group

Gynaecologic brachytherapy (BT) is a fundamental part of the treatment of gynaecologic carcinomas, especially if it is not possible to perform surgery. Traditionally, low-dose-rate BT and planning by orthogonal X-rays have been used. Nowadays, high-dose-rate BT is available in most centres, enabling it to enhance the dose optimisation, radiation protection and comfort of patients. Technical advantages in image and dosimetry have enabled the switch to CT and MRI as a means to define the tumour and the organs at risk, and to perform realistic anatomical-based planning. However a general consensus about the way the dose should be prescribed has not yet been reached. Some centres still use the Manchester points A as the reference to specify the dose. It was imperative to work towards a consensus and for this reason the VII Consensus Meeting of the Radiation Oncology Spanish Society (SEOR) Brachytherapy Group was organised. We agreed that some imaging system, CT or MRI, rather than the traditional system with orthogonal X-rays, should always be used for treatment planning. The dose should be prescribed by volume instead of points. And the given dose should be reported by adding external beam radiotherapy dose and BT dose as EQD2 (isoeffective dose). Caution is also recommended in changing the mode of planning and prescription, taking into account previous experience with points A.

Radiation transmission data for radionuclides and materials relevant to brachytherapy facility shielding

To address the limited availability of radiation shielding data for brachytherapy as well as some disparity in existing data, Monte Carlo simulation was used to generate radiation transmission data for 60Co, 137CS, 198Au, 192Ir 169Yb, 170Tm, 131Cs, 125I, and 103pd photons through concrete, stainless steel, lead, as well as lead glass and baryte concrete. Results accounting for the oblique incidence of radiation to the barrier, spectral variation with barrier thickness, and broad beam conditions in a realistic geometry are compared to corresponding data in the literature in terms of the half value layer (HVL) and tenth value layer (TVL) indices. It is also shown that radiation shielding calculations using HVL or TVL values could overestimate or underestimate the barrier thickness required to achieve a certain reduction in radiation transmission. This questions the use of HVL or TVL indices instead of the actual transmission data. Therefore, a three-parameter model is fitted to results of this work to facilitate accurate and simple radiation shielding calculations.

Monte Carlo study of the dose rate distributions for the Ir2.A85-2 and Ir2.A85-1 Ir-192 afterloading sources

The two most commonly used modalities of cancer treatment in clinical brachytherapy practice today are high dose rate (HDR) and pulsed dose rate (PDR) brachytherapy. In a clinical treatment, quality dose rate distribution data sets of the brachytherapy sources are required for each source model. The purpose of this study is to obtain detailed dose rate distributions around the new BEBIG HDR and PDR Ir-192 brachytherapy sources. These distributions will then be used as input data in the treatment planning systems dedicated to brachytherapy and its calculations can be verified. The Monte Carlo method was used to obtain the dose rate distributions around the sources studied, taking into account the AAPM-ESTRO recent recommendations. A complete dosimetric data set for the BEBIG Ir-192 HDR and PDR sources, types Ir2.A85-2 and Ir2.A85-1, were obtained. This dosimetric data set is composed of the TG-43 dosimetric functions and parameters and along-away dose rate table to facilitate quality control of treatment planning systems.

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.

Technical note: Dosimetric study of a new Co-60 source used in brachytherapy

The purpose of this study is to obtain the dosimetric parameters of a new Co-60 source used in high dose rate brachytherapy and manufactured by BEBIG (Eckert & Ziegler BEBIG GmbH, Germany). The Monte Carlo method has been used to obtain the dose rate distribution in the updated TG-43U1 formalism of the American Association of Physicists in Medicine. In addition, to aid the quality control process on treatment planning systems (TPS), a two-dimensional rectangular dose rate table, coherent with the TG-43U1 dose calculation formalism, is given. These dosimetric data sets can be used as input data of the TPS calculations and to validate them.

Dosimetric prerequisites for routine clinical use of photon emitting brachytherapy sources with average energy higher than 50 keva)

This paper presents the recommendations of the American Association of Physicists in Medicine (AAPM) and the European Society for Therapeutic Radiology and Oncology (ESTRO) on the dosimetric parameters to be characterized, and dosimetric studies to be performed to obtain them, for brachytherapy sources with average energy higher than 50 keV that are intended for routine clinical use. In addition, this document makes recommendations on procedures to be used to maintain vendor source strength calibration accuracy. These recommendations reflect the guidance of the AAPM and the ESTRO for its members, and may also be used as guidance to vendors and regulatory agencies in developing good manufacturing practices for sources used in routine clinical treatments.

A dosimetric study on the Ir-192 high dose rate Flexisource

In this work, the dose rate distribution of a new Ir-192 high dose rate source (Flexisource used in the afterloading Flexitron system, Isodose Control, Veenendaal, The Netherlands) is studied by means of Monte Carlo techniques using the GEANT4 code. The dosimetric parameters of the Task Group No. 43 Report (TG43) formalism and two-dimensional rectangular look-up tables have been obtained.

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.

Broad-beam transmission data for new brachytherapy sources, Tm-170 and Yb-169

The characteristics of the radionuclides (170)Tm and (169)Yb are highly interesting for their use as high dose-rate brachytherapy sources. The introduction of brachytherapy equipment containing these sources will lead to smaller required thicknesses of the materials used in radiation protection barriers compared with the use of conventional sources such as (192)Ir and (137)Cs. The purpose of this study is to determine the required thicknesses of protection material for the design of the protecting walls. Using the Monte Carlo method, transmission data were derived for broad-beam geometries through lead and concrete barriers, from which the first half value layer and tenth value layer are obtained. In addition, the dose reduction in a simulated patient was studied to determine whether transmission in the patient is a relevant factor in radiation protection calculations.

A Monte Carlo study of intersource effects in dome-type applicators loaded with LDR Cs-137 sources

In this study, the dose rate distributions produced by low dose rate Cs-137 sources loaded in afterloadable dome applicators are studied using the Monte Carlo method. Dose differences between Monte Carlo results and calculations done using the superposition principle are within 1-3% in front of the applicator and between 3 and 10% near and along the longitudinal source axis. Consequently, the real doses to lateral vaginal wall, rectum and bladder are very close to the doses estimated applying the superposition principle, while the dose to the vaginal cuff has been overestimated by up to 10%.

Monte Carlo dosimetric study of the BEBIG Co-60 HDR source

Although not as widespread as Ir-192, Co-60 is also available on afterloading equipment devoted to high dose rate brachytherapy, mainly addressed to the treatment of gynaecological lesions. The purpose of this study is to obtain the dosimetric parameters of the Co-60 source used by the BEBIG MultiSource remote afterloader (BEBIG GmbH, Germany) for which there are no dosimetric data available in the literature. The Monte Carlo code GEANT4 has been used to obtain the TG43 parameters and the 2D dose rate table in Cartesian coordinates of the BEBIG Co-60 HDR source. The dose rate constant, radial dose function and anisotropy function have been calculated and are presented in a tabular form as well as a detailed 2D dose rate table in Cartesian coordinates. These dosimetric datasets can be used as input data and to validate the treatment planning system calculations.

Technical note: Monte Carlo derivation of TG-43 dosimetric parameters for radiation therapy resources and 3M Cs-137 sources

In clinical brachytherapy dosimetry, a detailed dose rate distribution of the radioactive source in water is needed in order to plan for quality treatment. Two Cs-137 sources are considered in this study; the Radiation Therapy Resources 67-800 source (Radiation Therapy Resources Inc., Valencia, CA) and the 3M model 6500/6D6C source. A complete dosimetric dataset for both sources has been obtained by means of the Monte Carlo GEANT4 code. Dose rate distributions are presented in two different ways; following the TG43 formalism and in a 2D rectangular dose rate table. This 2D dose rate table is helpful for the TPS quality control and is fully consistent with the TG43 dose calculation formalism. In this work, several improvements to the previously published data for these sources have been included: the source asymmetries were taken explicitly into account in the MC calculations, TG43 data were derived directly from MC calculations, the data radial range was increased, the angular grid in the anisotropy function was increased, and TG43 data is now consistent with the along and away dose rate table as recommended by the TG43 update.

Monte Carlo calculation of the TG-43 dosimetric parameters of a new BEBIG Ir-192 HDR source

BACKGROUND AND PURPOSE

High dose rate (HDR) brachytherapy is a highly extended practice in clinical brachytherapy today. Quality dose rate distribution datasets of the HDR sources used in a clinical treatment are required. Because of the different source designs, a specific dosimetry dataset is required for each source model. In the recently published BRAPHYQS-ESTRO Report, an overview of available dosimetric data for all HDR Ir-192 sources is given, pointing out the lack of data for one of the sources that is used by the BEBIG MultiSource afterloading system (BEBIG GmbH, Germany). The purpose of this study is to obtain detailed dose rate distributions in liquid water media around this source.

MATERIAL AND METHODS

The Monte Carlo code GEANT4 was used to estimate dose rate in water and air-kerma strength around the Ir-192 source. All the details of the stainless steel encapsulated BEBIG HDR 1.1mm in external diameter has been included in the simulation.

RESULTS

A complete dosimetric dataset for the BEBIG Ir-192 HDR source is presented. TG43 dosimetric functions and parameters have been obtained as well as a 2D rectangular dose rate table, consistent with the TG43 dose calculation formalism. The dosimetric parameters and functions obtained for the BEBIG HDR source have been compared with that obtained in the literature for others HDR sources, showing that the use of specific datasets for this new source is justified.

CONCLUSIONS

This dataset can be used as input in the TPS and to validate its calculations. As policy of BRAPHYQS-ESTRO task group, this dataset will be incorporated to the website: available to users in excel format.

A dosimetric study of Leipzig applicators

PURPOSE

To obtain the absolute dose-rate distribution in liquid water for all six cup-shaped Leipzig applicators by means of an experimentally validated Monte Carlo (MC) code. These six applicators were used in high-dose-rate (HDR) afterloaders with the "classic" and v2 (192)Ir sources. The applicators have an inner diameter of 1, 2, and 3 cm, with the source traveling parallel or perpendicular to the contact surface.

METHODS AND MATERIALS

The MC GEANT4 code was used to obtain the dose-rate distribution in liquid water for the six applicators and the two HDR source models. To normalize the applicator output factors, a MC simulation for the "classic" and v2 sources in air was performed to estimate the air-kerma strength. To validate this specific application and to guarantee that realistic source-applicator geometry was considered, an experimental verification procedure was implemented in this study, in accordance with the TG43U1 recommendations. Thermolumniscent dosimeter chips and a parallel plate ionization chamber in a polymethyl methacrylate (PMMA) phantom were used to verify the MC results for the six applicators in a microSelectronHDR afterloader with the "classic" source. Dose-rate distributions dependence on phantom size has been evaluated using two different phantom sizes.

RESULTS

Percentage depth dose and off-axis profiles were obtained normalized at a depth of 3 mm along the central axis for both phantom sizes. A table of output factors, normalized to 1 U of source kerma strength at this depth, is presented. The dose measured in the PMMA phantom agrees within experimental uncertainties with the dose obtained by the MC GEANT4 code calculations. The phantom size influence on dose-rate distributions becomes significant at depths greater than 5 cm.

CONCLUSIONS

MC-detailed simulation was performed for the Nucletron Leipzig HDR applicators. The matrix data obtained, with a grid separation of 0.5 mm, can be used to build a dataset in a convenient format to model these distributions for routine use with a brachytherapy treatment planning system.

Monte Carlo and experimental derivation of TG43 dosimetric parameters for CSM-type Cs-137 sources

In this study, complete dosimetric datasets for the CSM2 and CSM3 Cs-137 sources were obtained using the Monte Carlo GEANT4 code. The application of this calculation method was experimentally validated with thermoluminescent dosimetry (TLD). Functions and parameters following the TG43 formalism are presented: the dose rate constant, the radial dose functional, and the anisotropy function. In addition, to aid the quality control process on treatment planning systems, a two-dimensional (2D) rectangular dose rate table (the traditional along-away table), coherent with the TG43 dose calculation formalism, is given. The data given in this study complement existing information for both sources on the following aspects: (i) the source asymmetries were considered explicitly in the Monte Carlo calculations, (ii) TG43 data were derived directly from Monte Carlo calculations, (iii) the radial range of the different tables was increased as well as the angular resolution in the anisotropy function, including angles close to the longitudinal source axis. The CSM2 source TG-43 data of Liu et al. [Med. Phys. 31, 477-483 (2004)] are not consistent with the Williamson 2D along-away data [Int. J. Radiat. Oncol., Biol., Phys. 15, 227-237 (1988)] at distances closer than approximately 2 cm from the source. The data presented here for this source are consistent with this 2D along-away table, and are suitable for use in clinical practice.

Monte Carlo evaluation of kerma in an HDR brachytherapy bunker

In recent years, the use of high dose rate (HDR) after-loader machines has greatly increased due to the shift from traditional Cs-137/Ir-192 low dose rate (LDR) to HDR brachytherapy. The method used to calculate the required concrete and, where appropriate, lead shielding in the door is based on analytical methods provided by documents published by the ICRP, the IAEA and the NCRP. The purpose of this study is to perform a more realistic kerma evaluation at the entrance maze door of an HDR bunker using the Monte Carlo code GEANT4. The Monte Carlo results were validated experimentally. The spectrum at the maze entrance door, obtained with Monte Carlo, has an average energy of about 110 keV, maintaining a similar value along the length of the maze. The comparison of results from the aforementioned values with the Monte Carlo ones shows that results obtained using the albedo coefficient from the ICRP document more closely match those given by the Monte Carlo method, although the maximum value given by MC calculations is 30% greater.

Dosimetric study of the 15mm ROPES eye plaque

The main aim of this paper is to make a study of dose-rate distributions obtained around the 15 mm, radiation oncology physics and engineering services, Australia (ROPES) eye plaque loaded with 125I model 6711 radioactive seeds. In this study, we have carried out a comparison of the dose-rate distributions obtained by the algorithm used by the Plaque Simulator (PS) (BEBIG GmbH, Berlin, Germany) treatment planning system with those obtained by means of the Monte Carlo method for the ROPES eye plaque. A simple method to obtain the dose-rate distributions in a treatment planning system via the superposition of the dose-rate distributions of a seed placed in the eye plaque has been developed. The method uses eye plaque located in a simplified geometry of the head anatomy and distributions obtained by means of the Monte Carlo code GEANT4. The favorable results obtained in the development of this method suggest that it could be implemented on a treatment planning system to improve dose-rate calculations. We have also found that the dose-rate falls sharply along the eye and that outside the eye the dose-rate is very low. Furthermore, the lack of backscatter photons from the air located outside the eye-head phantom produces a dose reduction negligible for distances from the eye-plaque r<1 cm but reaches up to 20% near the air-eye interface. Results showed that the treatment planning system lacks accuracy around the border of the eye (in the sclera and the surrounding area) due to the simplicity of the algorithm used. The BEBIG treatment planning system uses a global attenuation factor that takes into account the effect of the eye plaque seed carrier and the lack of backscatter photons caused by the metallic cover, which in the case of a ROPES eye plaque has a default value of T= 1 (no correction). In the present study, a global attenuation factor T=0.96 and an air-interface correction factor which improve on treatment planning system calculations were obtained.

Monte Carlo dosimetric study of Best Industries and Alpha Omega Ir-192 brachytherapy seeds

Ir-192 seeds are widely used in the USA for low dose rate interstitial brachytherapy. There are two commercially available models: those manufactured by Best Industries filtered with stainless steel, and those manufactured by Alpha-Omega seeds filtered with Pt. Newly developed 3D correction algorithms for brachytherapy are based on dosimetry data obtained on unbounded phantom size, allowing corrections for heterogeneities and actual tissue boundaries. Published dosimetric datasets for both seeds have been obtained under bounded conditions. The aim of the present study is to obtain dosimetric datasets for these seeds under full scatter conditions. The Monte Carlo GEANT4 code has been used to estimate air-kerma strength and dose rate in water around the Ir-192 seeds. Functions and parameters following the TG43 formalism are obtained and presented in tabular forms: the dose rate constant, the radial dose function, and the anisotropy function. Tables for the anisotropy factor have been obtained in order to apply punctual approximation. Differences between dose rate distributions for both seeds show that specific dataset must be used for each type of seed in clinical dosimetry. The data in the present study improve on published data in the following aspects: (i) dosimetric data were obtained under full scatter conditions, which affect dose values at distances greater than 4-5 cm from the source; (ii) the dose rate tables are given at greater distances from the source; and (iii) the spatial resolution in high dose gradient areas, such as those near the longitudinal source axis, has been improved.

Phantom size in brachytherapy source dosimetric studies

An important point to consider in a brachytherapy dosimetry study is the phantom size involved in calculations or experimental measurements. As pointed out by Williamson [Med. Phys. 18, 776-786 (1991)] this topic has a relevant influence on final dosimetric results. Presently, one-dimensional (1-D) algorithms and newly-developed 3-D correction algorithms are based on physics data that are obtained under full scatter conditions, i.e., assumed infinite phantom size. One can then assume that reference dose distributions in source dosimetry for photon brachytherapy should use an unbounded phantom size rather than phantom-like dimensions. Our aim in this paper is to study the effect of phantom size on brachytherapy for radionuclide 137Cs, 192Ir, 125I and 103Pd, mainly used for clinical purposes. Using the GEANT4 Monte Carlo code, we can ascertain effects on derived dosimetry parameters and functions to establish a distance dependent difference due to the absence of full scatter conditions. We have found that for 137Cs and 192Ir, a spherical phantom with a 40 cm radius is the equivalent of an unbounded phantom up to a distance of 20 cm from the source, as this size ensures full scatter conditions at this distance. For 125I and 103Pd, the required radius for the spherical phantom in order to ensure full scatter conditions at 10 cm from the source is R = 15 cm. A simple expression based on fits of the dose distributions for various phantom sizes has been developed for 137Cs and 192Ir in order to compare the dose rate distributions published for different phantom sizes. Using these relations it is possible to obtain radial dose functions for unbounded medium from bounded phantom ones.