Monte Carlo calculation of beam quality correction factors in proton beams using FLUKA.

Purpose.To provide Monte Carlo calculated beam quality correction factors (kQ) for monoenergetic proton beams using the Monte Carlo codefluka.Materials and methods.The Monte Carlo codeflukawas used to calculate the dose absorbed in a water-filled reference volume and the air-filled cavities of six plane-parallel and four cylindrical ionization chambers. The chambers were positioned at the entrance region of monoenergetic proton beams with energies between 60 and 250 MeV. Based on these dose values,fQas well askQfactors were calculated whilefQ0factors were taken from Andreoet al(2020Phys. Med. Biol.65095011).Results. kQfactors calculated in this work were found to agree with experimentally determinedkQfactors on the 1%-level, with only two exceptions with deviations of 1.4% and 1.9%. The comparison offQfactors calculated usingflukawithfQfactors calculated using the Monte Carlo codesgeant4 andpenhshowed a general good agreement for low energies, while differences for higher energies were pronounced. For high energies, in most cases the Monte Carlo codesflukaandgeant4 lead to comparable results while thefQfactors calculated withpenhare larger.Conclusion.flukacan be used to calculatekQfactors in clinical proton beams. The divergence of Monte Carlo calculatedkQfactors for high energies suggests that the role of nuclear interaction models implemented in the different Monte Carlo codes needs to be investigated in more detail.

Baumann KS ，Derksen L ，Witt M ，Michael Burg J ，Engenhart-Cabillic R ，Zink K ... -  《-》

Monte Carlo calculated ionization chamber correction factors in clinical proton beams - deriving uncertainties from published data.

Baumann KS ，Gomà C ，Wulff J ，Kretschmer J ，Zink K ... -  《-》

Monte Carlo calculation of beam quality correction factors in proton beams using TOPAS/GEANT4.

Baumann KS ，Kaupa S ，Bach C ，Engenhart-Cabillic R ，Zink K ... -  《-》

Monte Carlo simulated beam quality and perturbation correction factors for ionization chambers in monoenergetic proton beams.

Beam quality correction factors provided in current codes of practice for proton beams are approximated using the water-to-air mass stopping power ratio and by assuming the proton beam quality related perturbation correction factors to be unity. The aim of this work is to use Monte Carlo simulations to calculate energy dependent beam quality and perturbation correction factors for a set of nine ionization chambers in proton beams. The Monte Carlo code EGSnrc was used to determine the ratio of the absorbed dose to water and the absorbed dose to the sensitive air volume of ionization chambers f Q 0 related to the reference photon beam quality (60 Co). For proton beams, the quantity f Q was simulated with GATE/Geant4 for five monoenergetic beam energies between 70 MeV and 250 MeV. The perturbation correction factors for the air cavity, chamber wall, chamber stem, central electrode, and displacement effect in proton radiation were investigated separately. Additionally, the correction factors of cylindrical chambers were investigated with and without consideration of the effective point of measurement. The perturbation factors p Q were shown to deviate from unity for the investigated chambers, contradicting the assumptions made in dosimetry protocols. The beam quality correction factors for both plane-parallel and cylindrical chambers positioned with the effective point of measurement at the measurement depth were constant within 0.8%. An increase of the beam quality correction factors determined for cylindrical ionization chambers placed with their reference point at the measurement depth with decreasing energy is attributed to the displacement perturbation correction factors p dis , which were up to 1.045 ± 0.1% for the lowest energy and 1.005 ± 0.1% for the highest energy investigated. Besides p dis , the largest perturbation was found for the chamber wall where the smallest p wall determined was 0.981 ± 0.3%. Beam quality correction factors applied in dosimetry with cylindrical chambers in monoenergetic proton beams strongly depend on the positioning method used. We found perturbation correction factors different from unity. Consequently, the approximation of ionization chamber perturbations in proton beams by the respective water-to-air mass stopping power ratio shall be revised.

Kretschmer J ，Dulkys A ，Brodbek L ，Stelljes TS ，Looe HK ，Poppe B ... -  《-》

Comparison of penh, fluka, and Geant4/topas for absorbed dose calculations in air cavities representing ionization chambers in high-energy photon and proton beams.

The purpose of this work is to analyze whether the Monte Carlo codes penh, fluka, and geant4/topas are suitable to calculate absorbed doses and f Q / f Q 0 ratios in therapeutic high-energy photon and proton beams. We used penh, fluka, geant4/topas, and egsnrc to calculate the absorbed dose to water in a reference water cavity and the absorbed dose to air in two air cavities representative of a plane-parallel and a cylindrical ionization chamber in a 1.25 MeV photon beam and a 150 MeV proton beam - egsnrc was only used for the photon beam calculations. The physics and transport settings in each code were adjusted to simulate the particle transport as detailed as reasonably possible. From these absorbed doses, f Q 0 factors, f Q factors, and f Q / f Q 0 ratios (which are the basis of Monte Carlo calculated beam quality correction factors k Q , Q 0 ) were calculated and compared between the codes. Additionally, we calculated the spectra of primary particles and secondary electrons in the reference water cavity, as well as the integrated depth-dose curve of 150 MeV protons in water. The absorbed doses agreed within 1.4% or better between the individual codes for both the photon and proton simulations. The f Q 0 and f Q factors agreed within 0.5% or better for the individual codes for both beam qualities. The resulting f Q / f Q 0 ratios for 150 MeV protons agreed within 0.7% or better. For the 1.25 MeV photon beam, the spectra of photons and secondary electrons agreed almost perfectly. For the 150 MeV proton simulation, we observed differences in the spectra of secondary protons whereas the spectra of primary protons and low-energy delta electrons also agreed almost perfectly. The first 2 mm of the entrance channel of the 150 MeV proton Bragg curve agreed almost perfectly while for greater depths, the differences in the integrated dose were up to 1.5%. penh, fluka, and geant4/topas are capable of calculating beam quality correction factors in proton beams. The differences in the f Q 0 and f Q factors between the codes are 0.5% at maximum. The differences in the f Q / f Q 0 ratios are 0.7% at maximum.

Baumann KS ，Horst F ，Zink K ，Gomà C ... -  《-》