Radiation dose around a PET scanner installation: comparison of Monte Carlo simulations, analytical calculations and experimental results.

PURPOSE: Monte Carlo study of radiation transmission around areas surrounding a PET room. METHODS: An extended population of patients administered with (18)F-FDG for PET-CT investigations was studied, collecting air kerma rate and gamma ray spectra measurements at a reference distance. An MC model of the diagnostic room was developed, including the scanner and walls with variable material and thickness. MC simulations were carried out with the widely used code GEANT4. RESULTS: The model was validated by comparing simulated radiation dose values and gamma ray spectra produced by a volumetric source with experimental measurements; ambient doses in the surrounding areas were assessed for different combinations of wall materials and shielding and compared with analytical calculations, based on the AAPM Report 108. In the range 1.5-3.0 times of the product between the linear attenuation coefficient and thickness of an absorber (µ x), it was observed that the effectiveness of different combinations of shielding is roughly equivalent. An extensive tabulation of results is given in the text. CONCLUSIONS: The validation tests performed showed a satisfactory agreement between the simulated and expected results. The simulated dose rates incident on, and transmitted by the walls in our model of PET scanner room, are generally in good agreement with analytical estimates performed using the AAPM Publication No. 108 method. This provides an independent confirmation of AAPM's approach. Even in this specific field of application, GEANT4 proved to be a relevant and accurate tool for dosimetry estimates, shielding evaluation and for general radiation protection use.

A free database of radionuclide voxel S values for the dosimetry of nonuniform activity distributions.

The increasing availability of SPECT/CT devices with advanced technology offers the opportunity for the accurate assessment of the radiation dose to the biological target volume during radionuclide therapy. Voxel dosimetry can be performed employing direct Monte Carlo radiation transport simulations, based on both morphological and functional images of the patient. On the other hand, for voxel dosimetry calculations the voxel S value method can be considered an easier approach than patient-specific Monte Carlo simulations, ensuring a good dosimetric accuracy at least for anatomic regions which are characterized by uniform density tissue. However, this approach has been limited because of the lack of tabulated S values for different voxel dimensions and radionuclides. The aim of this work is to provide a free dataset of values which can be used for voxel dosimetry in targeted radionuclide studies. Seven different radionuclides (89Sr, 90Y, 131I, 153Sm, 177Lu, 186Re, 188Re), and 13 different voxel sizes (2.21, 2.33, 2.4, 3, 3.59, 3.9, 4, 4.42, 4.8, 5, 6, 6.8 and 9.28 mm) are considered. Voxel S values are calculated performing simulations of monochromatic photon and electron sources in two different homogeneous tissues (soft tissue and bone) with DOSXYZnrc code, and weighting the contributions on the basis of the radionuclide emission spectra. The outcomes are validated by comparison with Monte Carlo simulations obtained with other codes (PENELOPE and MCNP4c) performing direct simulation of the radionuclide emission spectra. The differences among the different Monte Carlo codes are of the order of a few per cent when considering the source voxel and the bremsstrahlung tail, whereas the highest differences are observed at a distance close to the maximum continuous slowing down approximation range of electrons. These discrepancies would negligibly affect dosimetric assessments. The dataset of voxel S values can be freely downloaded from the website www.medphys.it.

Simulation of radioactive decay in GEANT Monte Carlo codes: comparison between spectra and efficiencies computed with sch2for and G4RadioactiveDecay.

Two CERN Monte Carlo codes, i.e. GEANT3.21 and GEANT4, were compared. The specific routine (sch2for), implemented in GEANT3.21 to simulate a disintegration process, and the G4RadioactiveDecay class, provided by GEANT4, were used for the computation of the full-energy-peak and total efficiencies of several radionuclides. No reference to experimental data was involved. A level of agreement better than 1% for the total efficiency and a deviation lower than 3.5% for the full-energy-peak efficiencies were found.

Differences among Monte Carlo codes in the calculations of voxel S values for radionuclide targeted therapy and analysis of their impact on absorbed dose evaluations.

Several updated Monte Carlo (MC) codes are available to perform calculations of voxel S values for radionuclide targeted therapy. The aim of this work is to analyze the differences in the calculations obtained by different MC codes and their impact on absorbed dose evaluations performed by voxel dosimetry. Voxel S values for monoenergetic sources (electrons and photons) and different radionuclides (90Y, 131I, and 188Re) were calculated. Simulations were performed in soft tissue. Three general-purpose MC codes were employed for simulating radiation transport: MCNP4C, EGSnrc, and GEANT4. The data published by the MIRD Committee in Pamphlet No. 17, obtained with the EGS4 MC code, were also included in the comparisons. The impact of the differences (in terms of voxel S values) among the MC codes was also studied by convolution calculations of the absorbed dose in a volume of interest. For uniform activity distribution of a given radionuclide, dose calculations were performed on spherical and elliptical volumes, varying the mass from 1 to 500 g. For simulations with monochromatic sources, differences for self-irradiation voxel S values were mostly confined within 10% for both photons and electrons, but with electron energy less than 500 keV, the voxel S values referred to the first neighbor voxels showed large differences (up to 130%, with respect to EGSnrc) among the updated MC codes. For radionuclide simulations, noticeable differences arose in voxel S values, especially in the bremsstrahlung tails, or when a high contribution from electrons with energy of less than 500 keV is involved. In particular, for 90Y the updated codes showed a remarkable divergence in the bremsstrahlung region (up to about 90% in terms of voxel S values) with respect to the EGS4 code. Further, variations were observed up to about 30%, for small source-target voxel distances, when low-energy electrons cover an important part of the emission spectrum of the radionuclide (in our case, for 131I). For 90Y and 188Re, the differences among the various codes have a negligible impact (within few percents) on convolution calculations of the absorbed dose; thus either one of the MC programs is suitable to produce voxel S values for radionuclide targeted therapy dosimetry. However, if a low-energy beta-emitting radionuclide is considered, these differences can affect also dose depositions at small source-target voxel distances, leading to more conspicuous variations (about 9% for 1311) when calculating the absorbed dose in the volume of interest.