Characteristic analysis of a photon converter in a medical LINAC using GEANT4

1. Introduction

Medical LINAC (linear accelerator), which is an equipment used in radiotherapy, utilizes a high energy photon beam to treat deeply located tumors. Generally, neutrons are produced by photonuclear reactions in the head of the medical LINAC using photon beams over 8 MV [1]. Neutrons with high relative biological effectiveness (RBE) and permeability can cause many problems in terms of radiation protection [2]. The threshold energies of photons producing neutrons in the head of a medical LINAC are 13.06 MV (Aluminum), 13.38 MV (Iron), 8.06 MV (Gold), and 7.41 MV (Tungsten) [3].

The medical LINAC used in clinics cannot fundamentally block the production of photoneutrons because of the electron acceleration principle and their physical characteristics. Therefore, a study is required to analyze the characteristics and effects of photoneutrons. However, it is difficult to measure neutrons in clinics because their energy is distributed continuously and it is hard to separate them from other secondary particles [4]. The GEANT4 simulation can calculate the particle transport by considering the detailed characteristics of the material and can predict the dose distribution with higher accuracy than the conventional dose calculation algorithm. Furthermore, the scope of its application to the medical physics field is expanding [5] and the reliability of the simulation has been verified [6].

In this paper, we aim to analyze the characteristics of a photon converter in a medical LINAC using GEANT4 (Geometry and Tracking). For this, the photon and neutron spectrums of the optimized photon converter [7] with tungsten were simulated in 2~24 MeV electron beams, and the ratio of the photons with energies over 7.41 MV was analyzed. Furthermore, we analyzed the ratio of photons with energies over 7.41 MV that can produce neutrons.

2. Design of GEANT4 simulation environment converter and photon converter

An electron incident on a target loses its kinetic energy, which is moved until it comes to a stop in the material. Photon transport induces a bremsstrahlung spectrum with the maximum energy of the electron beam. The factors influencing the bremsstrahlung cross-sectional area are the energy of the electron beam, the atomic number of the target material, and the thickness of the target [8]. The optimized photon converter [7] is made of tungsten, and is cylindrical in shape with a radius of 1.5 cm and thickness of 1 mm as shown in Fig. 1.

Fig. 1The concept and the design of the photon converter

a) The designed photon converter

b) The concept of the photon converter

GEANT4 9.10 p02 [9] was used for the simulation to analyze the characteristics of the optimized photon converter. The cross-sectional data and the particle transport method were sourced from GEANT4 physics reference manual [9], and materials are defined by NIST [10]. The bremsstrahlung model utilizes the Begert and Seltzer model [11]. During the simulation, the electron beam was vertically incident on the photon converter at a position 2 mm away from the center point of the incidence plane. The simulation used 107 electron beams to minimize the statistical uncertainty. Fig. 2 shows a GEANT4 simulation result screenshot for electron beams incident on the photon converter.

Fig. 2GEANT4 simulation result screenshot

3. Analysis of GEANT4 simulation results and photon converter characteristics

Fig. 3 shows the photon yield spectra of the photon converters for 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, and 24 MeV electron beams. The characteristic radiation at all beam energies was produced at the 520 kV. The results are due to the characteristics of the photon, which is less affected by the characteristic radiation in photon beam over 300 kV. Furthermore, it was observed that the characteristic radiation appears at relatively high energy only 520 kV because the characteristic radiation with low energy is absorbed in the photon converter.

Table 1 shows the results of the analysis of the photon beam energy spectrum, which serves as a reference for investigating the photoneutrons produced by photonuclear reactions during the interactions between photons and matter. The number and content ratios of photons obtained over 7.41 MV are provided against the cumulative number of measured photons at 1 cm below the photon converter.

Fig. 3Photon yield spectrum of photon converter for different electron beam energies

Fig. 4 shows the neutron yield spectra of the photon converters at 2, 4, 6, 8, 10, 12, 14, 16, 18, 20 22, and 24 MeV electron beams. The maximum neutron yield was measured as 5.0-7 and 6.7-7 neutron yield/e at 10 MeV and 22 MeV electron beams respectively.

Fig. 4The neutron yield of the photon converter for different electron beam energies

4. Conclusions

The medical LINAC cannot block the production of photoneutrons because of the electron acceleration principle and their physical characteristics. Therefore, this paper proposes basic data regarding radiation protection by analyzing the characteristics of photoneutrons. To perform this analysis, we modeled the photon converter in the medical LINAC head using GEANT4, which is a Monte Carlo simulation. Simulations were performed to analyze the photoneutron characteristics at 2-24 MeV electron beams, which are typically used in clinics.

In the simulation results, we confirmed that the characteristic radiation is observed at 520 kV and the ratio of photons obtained with energy over 7.41 MV increased as the electron energy increased. Furthermore, the neutron. generation from the photon was maximum at 10 MeV electron beam, and the production of neutrons increased as the electron energy increased. As the neutrons have a high RBE, even low doses can be harmful to the patient. There is also the problem of causing exposure to the operator by radiating to the surrounding environment. Therefore, the neutrons generated from the photon converters must be considered when designing the shields. The results of this study are expected to be used to improve medical LINACs by solving these problems.