1
|
Nagake Y, Yasui K, Ooe H, Ichihara M, Iwase K, Toshito T, Hayashi N. Investigation of ionization chamber perturbation factors using proton beam and Fano cavity test for the Monte Carlo simulation code PHITS. Radiol Phys Technol 2024; 17:280-287. [PMID: 38261133 DOI: 10.1007/s12194-024-00777-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 12/29/2023] [Accepted: 01/02/2024] [Indexed: 01/24/2024]
Abstract
The reference dose for clinical proton beam therapy is based on ionization chamber dosimetry. However, data on uncertainties in proton dosimetry are lacking, and multifaceted studies are required. Monte Carlo simulations are useful tools for calculating ionization chamber dosimetry in radiation fields and are sensitive to the transport algorithm parameters when particles are transported in a heterogeneous region. We aimed to evaluate the proton transport algorithm of the Particle and Heavy Ion Transport Code System (PHITS) using the Fano test. The response of the ionization chamber f Q and beam quality correction factors k Q were calculated using the same parameters as those in the Fano test and compared with those of other Monte Carlo codes for verification. The geometry of the Fano test consisted of a cylindrical gas-filled cavity sandwiched between two cylindrical walls. f Q was calculated as the ratio of the absorbed dose in water to the dose in the cavity in the chamber. We compared the f Q calculated using PHITS with that of a previous study, which was calculated using other Monte Carlo codes (Geant4, FULKA, and PENH) under similar conditions. The flight mesh, a parameter for charged particle transport, passed the Fano test within 0.15%. This was shown to be sufficiently accurate compared with that observed in previous studies. The f Q calculated using PHITS were 1.116 ± 0.002 and 1.124 ± 0.003 for NACP-02 and PTW-30013, respectively, and the k Q were 0.981 ± 0.008 and 1.027 ± 0.008, respectively, at 150 MeV. Our results indicate that PHITS can calculate the f Q and k Q with high precision.
Collapse
Affiliation(s)
- Yuya Nagake
- Graduate School of Health Sciences, Fujita Health University, Toyoake, Japan
| | - Keisuke Yasui
- School of Medical Sciences, Fujita Health University, Toyoake, Japan.
| | - Hiromu Ooe
- Graduate School of Health Sciences, Fujita Health University, Toyoake, Japan
| | | | - Kaito Iwase
- Graduate School of Health Sciences, Fujita Health University, Toyoake, Japan
| | - Toshiyuki Toshito
- Nagoya Proton Therapy Center, Nagoya City University West Medical Center, Nagoya, Japan
| | - Naoki Hayashi
- School of Medical Sciences, Fujita Health University, Toyoake, Japan
| |
Collapse
|
2
|
Urago Y, Sakama M, Sakata D, Fukuda S, Katayose T, Chang W. Monte Carlo-calculated beam quality and perturbation correction factors validated against experiments for Farmer and Markus type ionization chambers in therapeutic carbon-ion beams. Phys Med Biol 2023; 68:185013. [PMID: 37579752 DOI: 10.1088/1361-6560/acf024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Accepted: 08/14/2023] [Indexed: 08/16/2023]
Abstract
Objective. In current dosimetry protocols, the estimated uncertainty of the measured absorbed dose to waterDwin carbon-ion beams is approximately 3%. This large uncertainty is mainly contributed by the standard uncertainty of the beam quality correction factorkQ. In this study, thekQvalues in four cylindrical chambers and two plane-parallel chambers were calculated using Monte Carlo (MC) simulations in the plateau region. The chamber-specific perturbation correction factorPof each chamber was also determined through MC simulations.Approach.kQfor each chamber was calculated using MC code Geant4. The simulatedkQratios in subjected chambers and reference chambers were validated through comparisons against our measured values. In the measurements in Heavy-Ion Medical Accelerator in Chiba,kQratios were obtained fromDwvalues of60Co, 290- and 400 MeV u-1carbon-ion beams that were measured with the subjected ionization chamber and the reference chamber. In the simulations,fQ(the product of the water-to-air stopping power ratio andP) was acquired fromDwand the absorbed dose to air calculated in the sensitive volume of each chamber.kQvalues were then calculated from the simulatedfQand the literature-extractedWairand compared with previous publications.Main results. The calculatedkQratios in the subjected chambers to the reference chamber agreed well with the measuredkQratios. ThekQuncertainty was reduced from the current recommendation of approximately 3% to 1.7%. ThePvalues were close to unity in the cylindrical chambers and nearly 1% above unity in the plane-parallel chambers.Significance. ThekQvalues of carbon-ion beams were accurately calculated in MC simulations and thekQratios were validated through ionization chamber measurements. The results indicate a need for updating the current recommendations, which assume a constantPof unity in carbon-ion beams, to recommendations that consider chamber-induced differences.
Collapse
Affiliation(s)
- Yuka Urago
- Department of Radiological Science, Graduate School of Human Health Science, Tokyo Metropolitan University, Tokyo, Japan
- QST Hospital, National Institutes for Quantum Science and Technology, Chiba, Japan
| | - Makoto Sakama
- QST Hospital, National Institutes for Quantum Science and Technology, Chiba, Japan
| | | | - Shigekazu Fukuda
- QST Hospital, National Institutes for Quantum Science and Technology, Chiba, Japan
| | | | - Weishan Chang
- Department of Radiological Science, Graduate School of Human Health Science, Tokyo Metropolitan University, Tokyo, Japan
| |
Collapse
|
3
|
Marot M, Jäger F, Greilich S, Karger CP, Jäkel O, Burigo LN. Monte Carlo simulation for proton dosimetry in magnetic fields: Fano test and magnetic field correction factors kBfor Farmer-type ionization chambers. Phys Med Biol 2023; 68:175037. [PMID: 37567226 DOI: 10.1088/1361-6560/acefa1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Accepted: 08/10/2023] [Indexed: 08/13/2023]
Abstract
Objective. In this contribution we present a special Fano test for charged particles in presence of magnetic fields in the MC code TOol for PArticle Simulation (TOPAS), as well as the determination of magnetic field correction factorskBfor Farmer-type ionization chambers using proton beams.Approach. Customized C++ extensions for TOPAS were implemented to model the special Fano tests in presence of magnetic fields for electrons and protons. The Geant4-specific transport parameters,DRoverRandfinalRange,were investigated to optimize passing rate and computation time. ThekBwas determined for the Farmer-type PTW 30013 ionization chamber, and 5 custom built ionization chambers with same geometry but varying inner radius, testing magnetic flux density ranging from 0 to 1.0 T and two proton beam energies of 157.43 and 221.05 MeV.Main results. Using the investigated parameters, TOPAS passed the Fano test within 0.39 ± 0.15% and 0.82 ± 0.42%, respectively for electrons and protons. The chamber response (kB,M,Q) gives a maximum at different magnetic flux densities depending of the chamber size, 1.0043 at 1.0 T for the smallest chamber and 1.0051 at 0.2 T for the largest chamber. The local dose differencecBremained ≤ 0.1% for both tested energies. The magnetic field correction factorkB, for the chamber PTW 30013, varied from 0.9946 to 1.0036 for both tested energies.Significance. The developed extension for the special Fano test in TOPAS MC code with the adjusted transport parameters, can accurately transport electron and proton particles in magnetic field. This makes TOPAS a valuable tool for the determination ofkB. The ionization chambers we tested showed thatkBremains small (≤0.72%). To the best of our knowledge, this is the first calculations ofkBfor proton beams. This work represents a significant step forward in the development of MRgPT and protocols for proton dosimetry in presence of magnetic field.
Collapse
Affiliation(s)
- M Marot
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- University of Heidelberg, Faculty of Medicine, Heidelberg, Germany
- Heidelberg Institute for Radiation Oncology (HIRO) and National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - F Jäger
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- University of Heidelberg, Faculty of Physics and Astronomy, Heidelberg, Germany
| | - S Greilich
- Berthold Technologies GmbH & Co. KG, Business Units Radiation Protection/Bioanalytics, Bad Wildbad, Germany
| | - C P Karger
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- Heidelberg Institute for Radiation Oncology (HIRO) and National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - O Jäkel
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- Heidelberg Institute for Radiation Oncology (HIRO) and National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- Heidelberg Ion-Beam Therapy Center (HIT), University Hospital Heidelberg, Heidelberg, Germany
| | - L N Burigo
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- Heidelberg Institute for Radiation Oncology (HIRO) and National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| |
Collapse
|
4
|
Vedelago J, Karger CP, Jäkel O. A review on reference dosimetry in radiation therapy with proton and light ion beams: status and impact of new developments. RADIAT MEAS 2022. [DOI: 10.1016/j.radmeas.2022.106844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|