Reevaluation of Photon Activation Yields of 11C, 13N, and 15O for the Estimation of Activity in Gas and Water Induced by the Operation of Electron Accelerators for Medical Use
Article information
Abstract
Background
Activation of air and water in the electron linear accelerator for medical use has not been considered severely. By the new Japanese regulation for protection of radiation hazard, it became indispensable to evaluate of activation of air and water in the accelerator room. The measurement of induced activity in air and water components in the electron energy region of 10 to 20 MeV is very difficult, because this energy region is close to the threshold energy region of photonuclear reactions. Then, we measured the photonuclear reaction yields of 13N, 15O, and 11C by using the electron linear accelerator. Obtained data were compared with the data calculated by the Monte Carlo method.
Materials and Methods
An activation experiment was performed at the Research Center for Electron Photon Science, Tohoku University. Highly purified SiO2, Si3N4, and carbon disks were irradiated for 10 minutes by bremsstrahlung converted by a tungsten plate. Induced activity from C, N, and O was obtained. Monte Carlo calculation was performed using MCNP5 and AERY (DCHAIN-SP) to simulate the experimental condition. Cross section data were adopted the KAERI dataset.
Results and Discussion
In our experiment in hospital, calculated values were not agreed with experimental values. It might be three possible reasons as the cause of this deference, such as irradiation energy, calculation procedure and cross section data. Obtained data of this work, calculated and experimental values were good agreement with each other within one order. In this work, we used KAERI dataset of photonuclear reaction instead of JENDL. Therefore, it was found that the photonuclear cross section data of light elements are most important for yield calculation in these reactions.
Conclusion
Further improvement for calculation using a new dataset JENDL/PD-2015 and considering electron energy spreading will be needed.
Introduction
About 900 electron linear accelerators are currently installed in hospitals in Japan, and electron acceleration energies from 6 to 15 MeV are popularly used in cancer therapy. Bremsstrahlung irradiation produced by hitting the target with accelerated electrons is the primary form of irradiation used in cancer therapies. Recently, it has become necessary to evaluate the activation of air and water to confirm that the induced activity is lower than the limit set by Japanese regulations in the new radiation protection law of radiation safety. The activation of air was induced by bremsstrahlung irradiation emitted from the target. The major activity in air is 13N (T1/2=9.96 min) from the 14N(γ,n)13N reaction above electron energies of 10 MeV. If the bremsstrahlung irradiation hits the cooling water or water in a phantom, 15O (T1/2=122 sec) will be produced by the 16O(γ,n)15O reaction above electron energies of 15 MeV [1].1)
For medical use, the irradiation energy is near the threshold energy of the photonuclear reaction, and the irradiation time is very short. Therefore, the induced activity of air and water does not generally cause an internal dose effect for patients and therapists because of the low induced activity and short-lived radioisotopes. We studied the activation of air and water in various accelerator facilities, which used energies between 6 and 20 MeV, in order to confirm the safety levels of activity in air in the treatment room. The concentration of radioisotopes in air and water was confirmed to be negligible below 15 MeV [2].
In this experiment, the measured activities of 13N and 15O in air and 15O in water are not consistent with the results of Monte Carlo calculations obtained using the Monte Carlo N-particle (MCNP) transport code and DCHAIN-SP code for high energy particle induced radioactivity calculation code. The measured activity was higher than calculated value. Therefore, it is important to confirm the major source of this discrepancy; for example, the accuracy of the electron energy, the Monte Carlo calculation procedure, the nuclear reaction cross-section for activation calculation, and the radioactivity measurement could be the major sources of this discrepancy. We plan to study the yields of the photonuclear reactions of light elements at the electron energies of 10 to 20 MeV at the facility for nuclear physics and to compare the results of measurement and calculation.
Materials and Methods
1. Materials
The samples were highly purified SiO2, Si3N4, and carbon disks (10 mm in diameter and 4 mm in thickness, Tokyo-Denshi Kogyo Co. Ltd., Kawasaki, Japan), which were used for monitoring the reaction yields from O, N, and C, respectively. The photon flux on each sample was monitored using Au foils (20 μm in thickness, Nilaco Co. Tokyo, Japan) by monitoring the yield of the 197Au(γ,n)196Au reaction. The samples were sandwiched between Au monitors. The stacked samples were encapsulated in a quartz ampoule (1 mm in thickness).
2. Irradiation
An activation experiment was performed at the Research Center for Electron Photon Science, Tohoku University. The electron energy varied from 10 to 20 MeV. The electron beam current ranged from 20 to 70 μA, and the irradiation time was 10 minutes. The beam was extracted through the window of Ti foil (50 μm in thickness) and the BeO beam profile monitor (0.5 mm in thickness). The distance from the Ti window to the irradiation chamber was 19.6 cm in air. Samples and a tungsten converter were set in an irradiation chamber of stainless steel, which had cooling water circulating in the interior. Electron was converted to bremsstrahlung by hitting a tungsten plate converter (2 mm in thickness) inside the irradiation chamber, and a quartz ampoule was also set at a position 1 cm behind the converter. Irradiation setup is shown in Figure 1.
Layout of irradiation setup.
3. Measurement of radioactivity
Irradiated ampoules were transported to the radioisotope laboratory and unsealed to separate each sample and the monitors. Induced activity was measured by a portable Ge detector (GR2018, Canberra Industry Inc., Meriden, CT) coupled with MCA (Inspector 2000, Canberra Industry Inc., Meriden, CT). The ISOCS (LabSOCS) software (Canberra Industry Inc., Meriden, CT) was used for efficiency calculation. Counting was performed repeatedly. After the decay curve analysis of the full energy peak of the positron annihilation gamma ray (511 keV), the activity of 13N, 15O, and 11C was obtained at the end of irradiation.
4. Monte Carlo calculation
Monte Carlo calculation was performed by MCNP5 and Induced radioactivity calculation code, AERY (renovated DCAIN-SP2001). The Korea Atomic Energy Research Institute (KAERI) database of photonuclear cross-section data was primarily used [3]. Calculations were repeated for 1012 to 1019 electron events, which had a beam size of 3 mm in diameter. The irradiation geometry was depicted in detail. To adjust the calculated results, the measured radioactivity of the sample was corrected to the calculated position.
Results and Discussion
1. Induced activity
In case of SiO2 and Si3N4, the radioactivity induced from Si could not be observed. However, 18F activity was observed from Si3N4 irradiated from 10 to 20 MeV, and 13N activity was observed from SiO2 in the case of irradiation lower than 15 MeV. The activity of 18F seems to be caused by the impurity of fluorine, and 13N seems to be caused by contamination. Therefore, the decay analysis is important. After least squares fitting, each activity was corrected as the activity at the end of irradiation.
2. Yield of pertinent nuclear reaction
Figure 2 shows the photonuclear reaction yield curves of 13N, 15O, 11C, and 196Au from 10 to 20 MeV. All data are normalized the yield value of 196Au at 20 MeV. The measured yield curve of 196Au was in good agreement with the calculated yield. The Q-values of the 14N(γ,n)13N, 16O(γ,n)15O, and 12C(γ,n)11C reactions are −10.6, −15.7, and −18.7 MeV, respectively. However, very low activities of 13N, 15O, and 11C were observed below the threshold energies. In this experiment, the width of the electron energy was about 5%. Therefore, the small contribution of higher energy electrons should be considered.
Photonuclear reaction yield of 13N, 15O, 11C, and 196Au.
Figure 3 shows the relative yield to the 196Au yield. The measured yield of 13N was three to four times higher than the calculated value, and the measured yield of 15O was one order higher than the calculated value. The measured and calculated relative yield curves showed similar trends.
Relative yields of 3N, 15O, 11C to 196Au.
3. Activity in air components
In our previous work measured at the hospital [2], the activity of 13N induced in the air of the 250-L glove box under the beam exit window was 0.87 Bq·cm−3 at 15 MeV (50 Gy, 10 minutes) irradiation. On the other hand, the calculated activity was 0.009 Bq·cm−3. The calculation of the activity was performed with GENDLE/PD-2004 for the 14N(γ,n)13N reaction and LA-150 for the 16O(γ,n)15O reaction. Saeed et al. measured the activity of 13N induced in the Marinelli box filled with N2 gas at 15 MeV (45 Gy, 15 minutes) irradiation [4]. The obtained result was 0.765 Bq·cm−3. This data is in good agreement with our data.
Figures 4 and 5 show the production cross-section of the 14N(γ,n)13N and 16O(γ,n)15O reaction, respectively. To fit the calculated data to the experimental data, the selection of cross-section data is important. In this work, we used the data set of KAERI. In this case, the discrepancy between the experiment and calculation results became small.
Production cross section of 14N(γ,n)13N reaction from the photonuclear evaluation files.
Production cross section of 16O(γ,n)15O reaction from the photonuclear evaluation files.
Conclusion
The disagreement between the experiment and calculation results might be caused by the cross-section data near the threshold energy. The experiments that obtained the photonuclear reaction cross-section for light elements were performed between 1950 and 1970. In this work using the electron accelerator for nuclear physics, it was found that the cross-section data set provided by KAERI were suited for reproducing our experimental results. Therefore, reevaluation of these results is important, especially in the threshold energy regions. Now we have a plan to recalculate by using the JENDL/PD-2015, which is newly evaluated in Japan.2) Also a spread of electron acceleration energy is unknown in case of electron accelerators for medical use. The activation of air in the accelerator room might be considered the energy spreading, because the threshold energy of photonuclear reactions of light element is near the electron acceleration energy.
Acknowledgements
The authors would like to thank the staff of the Research Center for Electron Photon Science, Tohoku University, for their kind help in the irradiation experiment.
Notes
Swanson WP. An effective biological safety program for electron linear accelerators. SLAC-PUB-2643. 1980;126.
Iwamoto N, Kosako K, Murata T. Photonuclear data file. Proceedings of the 2015 Symposium of Nuclear Data. Tokai Japan. November 19–20, 2015.