Three Monte Carlo photon/electron transport codes (GEPTS, EGSnrc, and MCNP) are bench-marked against dose measurements in homogeneous (both low- and high-Z) media as well as at interfaces. A brief overview on physical models used by each code for photon and electron (positron) transport is given. Absolute calorimetric dose measurements for 0.5 and 1 MeV electron beams incident on homogeneous and multilayer media are compared with the predictions of the three codes. Comparison with dose measurements in two-layer media exposed to a 60Co gamma source is also performed. In addition, comparisons between the codes (including the EGS4 code) are done for (a) 0.05 to 10 MeV electron beams and positron point sources in lead, (b) high-energy photons (10 and 20 MeV) irradiating a multilayer phantom (water/steel/air), and (c) simulation of a 90Sr/90Y brachytherapy source. A good agreement is observed between the calorimetric electron dose measurements and predictions of GEPTS and EGSnrc in both homogeneous and multilayer media. MCNP outputs are found to be dependent on the energy-indexing method (Default/ITS style). This dependence is significant in homogeneous media as well as at interfaces. MCNP(ITS) fits more closely the experimental data than MCNP(DEF), except for the case of Be. At low energy (0.05 and 0.1 MeV), MCNP(ITS) dose distributions in lead show higher maximums in comparison with GEPTS and EGSnrc. EGS4 produces too penetrating electron-dose distributions in high-Z media, especially at low energy (<0.1 MeV). For positrons, differences between GEPTS and EGSnrc are observed in lead because GEPTS distinguishes positrons from electrons for both elastic multiple scattering and bremsstrahlung emission models. For the 60Co source, a quite good agreement between calculations and measurements is observed with regards to the experimental uncertainty. For the other cases (10 and 20 MeV photon sources and the 90Sr/90Y beta source), a good agreement is found between the three codes. In conclusion, differences between GEPTS and EGSnrc results are found to be very small for almost all media and energies studied. MCNP results depend significantly on the electron energy-indexing method.

The EGS4 Monte Carlo radiation transport code was used to systematically study the dose perturbation near planar and cylindrical air cavities in a water medium irradiated by megavoltage x-ray beams. The variables of the problem included x-ray energy, cavity shape and dimension, and depth of the cavity in water. The Monte Carlo code was initially validated against published measurements and its results were found to agree within 2% with the published measurements. The study results indicate that the dose perturbation is strongly dependent on x-ray energy, field size, depth, and size of cavity in water. For example, the Monte Carlo calculations show dose reductions of 42% and 18% at 0.05 and 2 mm, respectively, beyond the air-water interface distal to the radiation source for a 3 cm thick air slab irradiated by a single 5x5 cm2 15 MV beam. The dose reductions are smaller for a parallel-opposed pair of 5x5 cm2 15 MV x-ray beams, being 21% and 11% for the same depths. The combined set of Monte Carlo calculations showed that the dose reduction near an air cavity is greater for: (a) Smaller x-ray field size, (b) higher x-ray energy, (c) larger air-cavity size, and (d) smaller depth in water where the air cavity is situated. A potential clinical application of these results to the treatment of prostate cancer is discussed.

The radiation treatment with catheter-based beta-emitter sources is under clinical trials to prevent restenosis following interventional coronary procedures. There are still large uncertainties in the dose calculation due to the complicated treatment geometry. We present the Monte Carlo simulations to account for the dosimetric perturbations due to neighboring trained seeds, proximal/distal gold markers, and a stainless steel stent. A catheter-based beta-emitter system is modeled using the Monte Carlo code, MCNP4B. Dose distributions and dose rates are calculated in voxels (0.64x0.64x0.5 mm3) around the long cylindrical trains of 90Sr/Y source with and without the stent (at 1.92 mm from the source axis). For the total activity of 70 mCi (2.59x10(9) Bq), the dose around most of the source length (except for edge seeds and gold markers) varies from 40 to 0.23 cGy/s as the radial distance from the source axis (r) increases from 0.64 to 6.4 mm. At the prescription range of r = 1.5-4.0 mm, the dose gradient is very steep and the contribution of neighboring seeds to the dose is significant. The dose enhancement due to neighboring seeds (the so-called "train effect") varies from 9% to 64% as r increases from 0.64 to 5.2 mm. The doses at r = 2 mm from the last edge seed and the gold marker are about 80% and 40% of that of the nonedge seed (8.7 cGy/s), respectively. The dose enhancement due to the secondary electrons and the primary electrons scattered with the stent is shown to be about 9.3% in the voxel including the stent. However, as r increases beyond the stent (r = 2.0-6.4 mm), the dose is slightly reduced by 4%-12%, compared to that without the stent.

The dose parameters for the beta-particle emitting 90Sr/90Y source for intravascular brachytherapy (IVBT) have been calculated by different investigators. At a distant distance from the source, noticeable differences are seen in these parameters calculated using different Monte Carlo codes. The purpose of this work is to quantify as well as to understand these differences. We have compared a series of calculations using an EGS4, an EGSnrc, and the MCNP Monte Carlo codes. Data calculated and compared include the depth dose curve for a broad parallel beam of electrons, and radial dose distributions for point electron sources (monoenergetic or polyenergetic) and for a real 90Sr/90Y source. For the 90Sr/90Y source, the doses at the reference position (2 mm radial distance) calculated by the three code agree within 2%. However, the differences between the dose calculated by the three codes can be over 20% in the radial distance range interested in IVBT. The difference increases with radial distance from source, and reaches 30% at the tail of dose curve. These differences may be partially attributed to the different multiple scattering theories and Monte Carlo models for electron transport adopted in these three codes. Doses calculated by the EGSnrc code are more accurate than those by the EGS4. The two calculations agree within 5% for radial distance <6 mm.