Monte Carlo calculation of perturbation correction factors for air-filled ionization chambers in clinical proton beams using TOPAS/GEANT

Abstract

IntroductionCurrent dosimetry protocols for clinical protons using air-filled ionization chambers assume that the perturbation correction factor is equal to unity for all ionization chambers and proton energies. Since previous Monte Carlo based studies suggest that perturbation correction factors might be significantly different from unity this study aims to determine perturbation correction factors for six plane-parallel and four cylindrical ionization chambers in proton beams at clinical energies. Materials and methodsThe dose deposited in the air cavity of the ionization chambers was calculated with the help of the Monte Carlo code TOPAS/Geant4 while specific constructive details of the chambers were removed step by step. By comparing these dose values the individual perturbation correction factors pcel, pstem, psleeve, pwall, pcav⋅pdis as well as the total perturbation correction factor pQ were derived for typical clinical proton energies between 80 and 250MeV. ResultsThe total perturbation correction factor pQ was smaller than unity for almost every ionization chamber and proton energy and in some cases significantly different from unity (deviation larger than 1%). The maximum deviation from unity was 2.0% for cylindrical and 1.5% for plane-parallel ionization chambers. Especially the factor pwall was found to differ significantly from unity. It was shown that this is due to the fact that secondary particles, especially alpha particles and fragments, are scattered from the chamber wall into the air cavity resulting in an overresponse of the chamber. ConclusionPerturbation correction factors for ionization chambers in proton beams were calculated using Monte Carlo simulations. In contrast to the assumption of current dosimetry protocols the total perturbation correction factor pQ can be significantly different from unity. Hence, beam quality correction factors kQ,Q0 that are calculated with the help of perturbation correction factors that are assumed to be unity come with a corresponding additional uncertainty.