Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images

F. Botta, A. Mairani, R. F. Hobbs, A. Vergara Gil, M. Pacilio, K. Parodi, M. Cremonesi, M. A.Coca Pérez, A. Di Dia, M. Ferrari, F. Guerriero, G. Battistoni, G. Pedroli, G. Paganelli, L. A.Torres Aroche, G. Sgouros

Research output: Contribution to journalArticle

Abstract

Patient-specific absorbed dose calculation for nuclear medicine therapy is a topic of increasing interest. 3D dosimetry at the voxel level is one of the major improvements for the development of more accurate calculation techniques, as compared to the standard dosimetry at the organ level. This study aims to use the FLUKA Monte Carlo code to perform patient-specific 3D dosimetry through direct Monte Carlo simulation on PET-CT and SPECT-CT images. To this aim, dedicated routines were developed in the FLUKA environment. Two sets of simulations were performed on model and phantom images. Firstly, the correct handling of PET and SPECT images was tested under the assumption of homogeneous water medium by comparing FLUKA results with those obtained with the voxel kernel convolution method and with other Monte Carlo-based tools developed to the same purpose (the EGS-based 3D-RD software and the MCNP5-based MCID). Afterwards, the correct integration of the PET/SPECT and CT information was tested, performing direct simulations on PET/CT images for both homogeneous (water) and non-homogeneous (water with air, lung and bone inserts) phantoms. Comparison was performed with the other Monte Carlo tools performing direct simulation as well. The absorbed dose maps were compared at the voxel level. In the case of homogeneous water, by simulating 108 primary particles a 2% average difference with respect to the kernel convolution method was achieved; such difference was lower than the statistical uncertainty affecting the FLUKA results. The agreement with the other tools was within 3-4%, partially ascribable to the differences among the simulation algorithms. Including the CT-based density map, the average difference was always within 4% irrespective of the medium (water, air, bone), except for a maximum 6% value when comparing FLUKA and 3D-RD in air. The results confirmed that the routines were properly developed, opening the way for the use of FLUKA for patient-specific, image-based dosimetry in nuclear medicine.

Original languageEnglish (US)
Pages (from-to)8099-8120
Number of pages22
JournalPhysics in medicine and biology
Volume58
Issue number22
DOIs
StatePublished - Nov 21 2013

    Fingerprint

ASJC Scopus subject areas

  • Radiological and Ultrasound Technology
  • Radiology Nuclear Medicine and imaging

Cite this

Botta, F., Mairani, A., Hobbs, R. F., Vergara Gil, A., Pacilio, M., Parodi, K., Cremonesi, M., Pérez, M. A. C., Di Dia, A., Ferrari, M., Guerriero, F., Battistoni, G., Pedroli, G., Paganelli, G., Aroche, L. A. T., & Sgouros, G. (2013). Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images. Physics in medicine and biology, 58(22), 8099-8120. https://doi.org/10.1088/0031-9155/58/22/8099