The radial dose function of low-energy brachytherapy seeds in different solid phantoms: comparison between calculations with the EGSnrc and MCNP4C Monte Carlo codes and measurements.

The use of low-energy photon emitters for brachytherapy applications, as in the treatment of the prostate or of eye tumours, has drastically increased in the last few years. New seed models for 103Pd and 125I have recently been introduced. The American Association of Physicists in Medicine recommends that measurements are made to obtain the dose rate constant, the radial dose function and the anisotropy function. These results must then be compared with Monte Carlo calculations to finally obtain the dosimetric parameters in liquid water. We have used the results obtained during the characterization of the new InterSource (furnished by IBt, Seneffe, Belgium) palladium and iodine sources to compare two Monte Carlo codes against experiment for these low energies. The measurements have been performed in three different media: two solid water plastics, WT1 and RW1, and polymethylmetacrylate. The Monte Carlo calculations were made using two different codes: MCNP4C and EGSnrc. These codes use photon cross-section data of a different origin. Differences were observed between both sets of input data below 100 keV, especially for the photoelectric effect. We obtained differences in the radial dose functions calculated with each code, which can be explained by the difference between the input data. New cross-section data were then tested for both codes. The agreement between the calculations using these new libraries is excellent. The differences are within the statistical uncertainties of the calculations. These results were compared with the experimental data. A good agreement is reached for both isotopes and in the three phantoms when the measured values are corrected for the presence of the TLDs in the phantom.