Estimation of neutron and gamma radiation doses inside the concrete shield wall for 10 and 15 MV medical linear accelerators

Similar documents
Estimate of Photonuclear Reaction in a Medical Linear Accelerator Using a Water-Equivalent Phantom

Progress in Nuclear Science and Technology, Volume 6,

Neutron source strength measurements for Varian, Siemens, Elekta, and General Electric linear accelerators

Calculations of Photoneutrons from Varian Clinac Accelerators and Their Transmissions in Materials*

Neutron Fluence and Energy Spectra Around the Varian Clinac 2lOOC/23OOC Medical Accelerator

Measurement of induced radioactivity in air and water for medical accelerators

Energy response for high-energy neutrons of multi-functional electronic personal dosemeter

In our recent paper [1], applicability of the self-activation of an NaI scintillator has been studied for the photo-neutron monitoring around

Science, Ami, Ami-machi, Inashiki-gun, Ibaraki , Japan. Ibaraki , Japan. Japan. *Corresponding author:

Measurement of Neutron Flux in a Medical Compact cyclotron room with Boron-Containing water Self-Shielding

A Measuring System with Recombination Chamber for Photoneutron Dosimetry at Medical Linear Accelerators

Low energy linear acceleratordriven subcritical assembly

Estimation of Radioactivity and Residual Gamma-ray Dose around a Collimator at 3-GeV Proton Synchrotron Ring of J-PARC Facility

Benchmark Test of JENDL High Energy File with MCNP

Radiation safety of the Danish Center for Proton Therapy (DCPT) Lars Hjorth Præstegaard Dept. of Medical Physics, Aarhus University Hospital

Secondary neutron spectra from modern Varian, Siemens, and Elekta linacs with multileaf collimators

The Influence of Brass Compensator Thickness and Field Size on Neutron Contamination Spectrum in 18MV Elekta SL 75/25 Medical Linear Accelerator

The neutron dose equivalent evaluation and shielding at the maze entrance of a Varian Clinac 23EX treatment room

Comparison with simulations to experimental data for photoneutron reactions using SPring-8 Injector

Convertible source system of thermal neutron and X-ray at Hokkaido University electron linac facility

Monte Carlo Calculations Using MCNP4B for an Optimal Shielding Design. of a 14-MeV Neutron Source * James C. Liu and Tony T. Ng

PERSPECTIVES OF PERSONNEL EXTERNAL DOSIMETRY AT STANFORD LINEAR ACCELERATOR CENTER

Chapter V: Cavity theories

Neutronics Experiments for ITER at JAERI/FNS

Activities of the neutron standardization. at the Korea Research Institute of Standards and Science (KRISS)

Secondary Neutron Dose Measurement for Proton Line Scanning Therapy

Simulation Modeling in Dosimetry

Radiation Protection Dosimetry (2007), Vol. 126, No. 1 4, pp Advance Access publication 11 May 2007

Volume 1 No. 4, October 2011 ISSN International Journal of Science and Technology IJST Journal. All rights reserved

Neutron Sources in the Varian Clinac 2100~2300C Medical Accelerator Calculated by the EGS4 Code

Shielding Design for the Imaging and Medical Beamline at the Australian Synchrotron

ARTICLE. Progress in Nuclear Science and Technology Volume 4 (2014) pp Yoshiko Harima a*, Naohiro Kurosawa b and Yukio Sakamoto c

Secondary Particles Produced by Hadron Therapy

Neutron Shielding Properties Of Concrete With Boron And Boron Containing Mineral

I. INTRODUCTION EXPERIMENTAL

Georgia Institute of Technology. Radiation Detection & Protection (Day 3)

Integral Benchmark Experiments of the Japanese Evaluated Nuclear Data Library (JENDL)-3.3 for the Fusion Reactor Design

Introduction to Radiological Sciences Neutron Detectors. Theory of operation. Types of detectors Source calibration Survey for Dose

THE ACTIVE PERSONNEL DOSIMETER - APFEL ENTERPRISES SUPERHEATED DROP DETECTOR*

Recent Activities on Neutron Calibration Fields at FRS of JAERI

Chapter 4. QUANTIFYING THE HAZARD II: DATA & ANALYSIS. The dose equivalents for spheres in air with 10 cm radius centred at a point in the

Monte Carlo Analyses of X-Ray Absorption, Noise, and Detective Quantum Efficiency Considering Therapeutic X-Ray Spectrum in Portal Imaging Detector

Since the beam from the JNC linac is a very high current, low energy beam, energy loss induced in the material irradiated by the beam becomes very lar

Radiation protection issues in proton therapy

COMPARISON OF COMPUTER CODES APPLICABILITY IN SHIELDING DESIGN FOR HADRON THERAPY FACILITIES *

Characterization of Low-Energy (6-30 kev) Response of Polish TLDs (MTS-N, MCP-N) with Synchrotron Radiation and Determination of Some TLD Quantities.

5. Gamma and neutron shielding characteristics of concretes containing different colemanite proportions

Study on Radiation Shielding Performance of Reinforced Concrete Wall (2): Shielding Analysis

7. a XV-2 high spatial resolution lm detector (Kodak). Important parameters of these detectors are given in Table1. The ionization chambers and the di

NEUTRON EXPOSURE FROM ELECTROM LINEAR ACCELERATORS AND A PROTON ACCELERATOR: MEASUREMENTS AND SIMULATIONS

Radiation Protection At Synchrotron Radiation Facilities

ACTIVATION ANALYSIS OF DECOMISSIONING OPERATIONS FOR RESEARCH REACTORS

DESIGN OF A TREATMENT ROOM FOR AN 18-MV LINAC

Neutronic Design on a Small Accelerator based 7 Li (p, n) Neutron Source for Neutron Scattering Experiments

Contrabands detection with a low energy electron linac driven photoneutron source

Today, I will present the first of two lectures on neutron interactions.

Introduction. Neutron Effects NSEU. Neutron Testing Basics User Requirements Conclusions

Neutron Interactions Part I. Rebecca M. Howell, Ph.D. Radiation Physics Y2.5321

Available online at ScienceDirect. Physics Procedia 69 (2015 )

Specific Accreditation Criteria Calibration ISO/IEC Annex. Ionising radiation measurements

THE mono-energetic hadron beam such as heavy-ions or

The possibility to use energy plus transmutation set-up for neutron production and transport benchmark studies

APPLIED RADIATION PHYSICS

Design, construction and characterization of a portable irradiator to calibrate installed ambient dose equivalent monitors

Peter J. Biggs Ph.D. Massachusetts General Hospital Harvard Medical School Boston, MA 02114

NEUTRON SPECTROMETRY WITH BUBBLE DETECTORS

Fast-Neutron Production via Break-Up of Deuterons and Fast-Neutron Dosimetry

Induced photonuclear interaction by Rhodotron-TT MeV electron beam

Neutron fluence measurements of the Siemens Oncor linear accelerator utilizing gold foil activation

NEUTRON RADIATION FROM MEDICAL ELECTRON ACCELERATORS*

Nuclear Science Merit Badge Workbook

International Journal of Scientific & Engineering Research, Volume 7, Issue 2, February-2016 ISSN

1.E Neutron Energy (MeV)

New irradiation zones at the CERN-PS

Measurement of 40 MeV Deuteron Induced Reaction on Fe and Ta for Neutron Emission Spectrum and Activation Cross Section

Chapter VIII: Photographic films

Radiation Shielding of a 230 MeV Proton Cyclotron For Cancer Therapy

Fabrication of Radiation Sources for Educational Purposes from Chemical Fertilizers using Compressing and Forming Method

Neutron Metrology Activities at CIAE (2009~2010)

Development of a Dosimetric System using Spectrometric Technique suitable for Operational Radiation Dose Measurements and Evaluation

Linear attenuation coefficient calculation for both pure silicon (Si) and silicone supported with lead

The Radiological Physics Center s standard dataset for small field size output factors

Radiation Shielding of Extraction Absorbers for a Fermilab Photoinjector

Characterization of the 3 MeV Neutron Field for the Monoenergetic Fast Neutron Fluence Standard at the National Metrology Institute of Japan

Performance Test of the Electronic Personal Neutron Dosemeter in Neutron Fields Simulating Workplaces of MOX Fuel Fabrication Facilities

Radiation Therapy Study Guide

ESTIMATION OF 90 SCATTERING COEFFICIENT IN THE SHIELDING CALCULATION OF DIAGNOSTIC X-RAY EQUIPMENT

Evaluation and Measurements of Radioactive Air Emission and Off-Site Doses at SLAC

SLAC-PUB Submitted to Radiation Protection and Dosimetry. Work supported by Department of Energy contract DE-AC02-76SF00515

Units S H I E L D I N H = D Q. H: Dose equivalent (Sv) D: Dose (Gy) Q: Quality Factor. 1Sv = 1 J/Kg. 1Gy = 1 J/Kg

Shielding Design Considerations for Proton Therapy Facilities

Investigation Of The Effects Of Variation Of Neutron Source-Detector Distance On The Emitted Neutron Dose Equivalent

City University of Hong Kong

The scanning microbeam PIXE analysis facility at NIRS

Radiation Quantities and Units

TITLE: Air Kerma Primary Standard: Experimental and Simulation Studies on Cs-137

Altitude Variation of cosmic-ray neutron energy spectrum and ambient dose equivalent at Mt.Fuji in Japan

PHITS calculation of the radiation field in HIMAC BIO

neutrons in the few kev to several MeV Neutrons are generated over a wide range of energies by a variety of different processes.

Implementing the ISO Reference Radiations at a Brazilian Calibration Laboratory

Transcription:

DOI: 10.15669/pnst.4.280 Progress in Nuclear Science and Technology Volume 4 (2014) pp. 280-284 ARTICLE Estimation of neutron and gamma radiation doses inside the concrete shield wall for 10 and medical linear accelerators Toshioh Fujibuchi a*, Satoshi Obara b, Masaru Nakajima c, Nozomi Kitaura c and Tomoharu Sato c a Ibaraki Prefectural University of Health Sciences, 4669-2, Ami, Ami-machi, Inashiki-gun, Ibaraki-ken, 300-0394, Japan; b National Institute of Radiological Sciences, 4-9-1, Anagawa, Inage-ku, Chiba,Chiba-ken, 263-8555, Japan; c Center Institute Hospital, Japanese Foundation for Cancer Research, 3-8-31,Ariake, Koto-ku, Tokyo, 135-8550, Japan High-energy electron linear accelerators (linacs) are used in radiotherapy for tumor treatment. By operating at electron energies above the thresholds of the (γ, n) photonuclear reaction, undesirable photoneutrons and radioisotopes are generated by high-energy x-rays. These photonuclear reactions occur on the linac accelerator head components, which are made from various materials. In these reactions, one can assume that secondary neutrons are dispersed all over the treatment rooms. The present work observed ambient dose equivalent of photon and neutron, and number of etch pits of wide-energy range Neupit produced by thermal and fast neutrons inside a model of iron shielding wall. The neutron dose around the surface of the shielding wall was detected less than three orders of magnitude of photons at linac. The total neutron dose in a linac was about one order higher than that in the linac. In addition, the linac incident photons clearly generated secondary neutrons inside the iron shielding wall. Keywords: medical linear accelerator; iron shield wall; photonuclear reaction; neutron 1. Introduction 1 External beam radiotherapy is the most common form of radiotherapy for tumor treatment. High-energy electron linear accelerators (linacs) are used in external beam radiotherapy. Since the threshold for photoneutron production on lead is about 8 MeV, contamination due to neutrons is not a problem for therapy beams below 10 MV. At this energy, the production of neutrons is still quite low, but should be taken into account. By, neutron production increases by a factor of ten for 10 MV, and by 18 MV a further factor of two for [1]. Few measurements and calculations for the neutron dose rate distribution around linac shielding walls have been conducted previously [2-4]. Considering the large discrepancy reported between the radiation dose results of calculations and measurements [5], the generation of photoneutrons and several radioisotopes inside the shielding wall should be taken into account when the therapeutic beam energy is higher than the nuclear reaction thresholds [6]. By further increase of the incident photon energy, photoneutrons are produced in the structure materials such as iron of shielding wall. As far as radiation dose from the target is concerned, the dose of photoneutrons is smaller than the dose of photon, and can be ignored. However, shielding calculations considering only photons from the target *Corresponding author. Email: fujibuch@hs.med.kyushu-u.ac.jp have been performed and found out that the measured dose was greater than the calculated ones by about two orders of magnitude [6]. The literature evaluated behavior of radiation inside of wall by Monte-Carlo calculation and measured leakage dose using ion chamber and rem meter. Oishi et al. measured and dose rate distribution in the laminated shielding assembly irradiated by 18, 25 and 34 MeV photons [7]. But they did not measured photon and neutron dose inside shielding wall from 10 and linac. The purpose of this study was to increase the understanding of neutron behavior inside the shielding wall of a radiotherapy treatment room. In this study, we measured photon and neutron exposure dose inside a model of iron shielding wall. 2. Material and methods 2.1. Experiments in 10 and linacs Photoneutrons are mainly induced in the high-z materials on the linac head such as the target, the flattening filter, the jaws, the multi leaf collimators, and the head shielding [8]. Experiments were performed for two types of linacs to evaluate the energy dependence of the secondary neutrons generated inside the model of shielding wall. The accelerating electron energy of the linac at Ibaraki Prefectural University of Health Sciences was (EXL-15SP: Mitsubishi Electric 2014 Atomic Energy Society of Japan. All rights reserved.

Progress in Nuclear Science and Technology, Volume 4, 2014 281 Figure 1. Experimental shield array of the concrete and iron shielding (cross-sectional plane view). Co.) with copper target and tungsten collimators, and that of the linac at the Japanese Foundation for Cancer Research, Center Institute Hospital of Radiation Oncology, Tokyo was (CLINACix: Varian Medical systems) with tungsten and copper target, and tungsten collimators. Shielding for neutrons requires hydrogen-containing material and that for photons requires high-z material. Experiments were performed for a pair of models to allow a comparison between concrete-block shielding and shielding composed of iron sandwiched between concrete blocks. The experimental arrangement is shown in Figure 1. The layout II is configured to simulate a typical concrete shielding wall with built-in iron gamma-ray shielding. The photon beam and other scattered or photonuclear production particles were irradiated on the surface of the shielding wall perpendicularly at its center. The irradiated area was 10 cm 10 cm, and the surface cross section of the shielding was 30 cm 30 cm. The total thickness of the concrete iron shielding wall was 56 cm, where two 18-cm-thick concrete layers sandwiched a 20-cm-thick iron layer. Each concrete layer consisted of three slabs cut to 6 cm thickness with an area of 30 cm 30 cm and density is 2.1 g/cm 3. The concrete is a composition of four components blended in a certain proportion: cement, crushed stones, sand, and water. The main weight compositions of concrete are calcium (41.6%), silicon (11.0%), aluminum (2.3%), iron (2.2%), sulfur (1.4%), potassium (0.56%), phosphorus (0.6%) and sodium (0.08%)[9]. The weight compositions of hydrogen, the carbon, and oxygen are not measured. The array of iron blocks (5 cm 10 cm 20 cm, 7.9 g/cm 3 ) was 20 cm in thickness. 2.2. Measurements of photon and neutron dose rate distributions inside the shielding wall In this study, we used Glass Badge using silver activated phosphate glass dosimeter (GD-450: Chiyoda Technol Co.) for photon dose measurement, and wide-energy range Neupit (WNP: Chiyoda Technol Co.) that is allyl diglycol carbonate plastic solid state nuclear track detector for neutron dose measurement. They were specially developed for large scale personal monitoring. GD-450 has two plastic filters and three metal filters (aluminum, copper and tin), and this composition enables to measure 10 kev to 10 MeV photon and 300 kev to 3 MeV beta radiation [10]. GD-450 size is 45 mm 13 mm 5 mm and 5 g weight. WNP has two kind of charged particle radiator boron nitride (BN) and polyethylene to cover wide energy range of neutron measurement such as 0.025 ev to 15 MeV continuously[10, 11]. WNP size is 27 mm 9 mm 2 mm and 0.5 g weight. GD-450 and WNP were placed between the concrete slabs and in the iron array. The measurement points were shown in Figure 1. Use of the WNP that is commonly employed to detect wide energy ranges of neutrons and various particles is suitable for simple and easy measurements. The analysis of measured GD-450 and WNP were executed by Chiyoda technol corporation. 3. Results and discussion 3.1. Photon and neutron ambient dose equivalent inside a model of shielding wall The experimental values for both photon and neutron dose equivalents inside a model of shielding wall were normalized to a monitor unit (MU). A MU is the amount of charge recorded in the ionization chamber mounted in the head of linac. Linac is calibrated to give a particular absorbed dose under particular conditions. 100 MU gives an absorbed dose of 1 Gy at isocenter, for a 10 cm 10 cm field. The isocenter is the point in space where radiation beams intersect when the gantry is rotated during beam-on. GD-450 and WNP were calibrated on a condition of free air, and they gave the ambient dose equivalents in this measurement. We observed that the attenuation of photon and neutron doses depended on the observation distance from the surface of the shielding wall. Figure 2 shows the photon ambient dose equivalent inside the shielding walls of two configurations and two energies. The photon doses of 10 and were in agreement. The photon attenuation profiles of the dose equivalent inside the concrete-only shielding wall decreased exponentially. Figure 3 shows the neutron ambient dose equivalent inside the shielding walls. The neutron dose equivalent increased inside the iron shielding wall for both electron energies. The neutron dose equivalent was increased by the inserted iron slab about ten times in and 3 times in. The reason the neutron dose increases includes photonuclear reaction of iron. We did not measure energy spectrum of incident photons; however, the dose increase suggested that even the linac emitted photons of energies higher than the threshold

282 T. Fujibuchi et al. Ambient dose equivalent (msv/mu) 10MV Etch pits (/cm 2 /MU) Figure 2. Photon ambient dose equivalent inside a model of shielding wall. Figure 4. Number of etch pits of WNP produced by fast neutrons inside a model of shielding wall. Ambient dose equivalent (msv/mu) 1.E-03 1.E-04 1.E-05 (Concrete & Iron) (Concrete & Iron) Etch pits (/cm 2 /MU) 10MV Figure 3. Neutron ambient dose equivalent inside a model of shielding wall. Figure 5. Number of etch pits of WNP produced by thermal neutrons inside a model of shielding wall. energy of iron (γ, n) reaction (7.6 MeV for 57 Fe [abundance:2.1%], 10 MeV for 58 Fe [0.28%], 11.2 MeV for 56 Fe [91.6%], 13.4 MeV for 54 Fe [5.6%]). Moreover, attenuation of the neutron in concrete is larger than iron [7]. It is thought that the difference of attenuation with the medium also influences. In this study, the neutron dose equivalent rate decreased to about 1/100 in the front and the rear shielding block at 10 and X ray. The block was composed of the concrete of 36 cm and the iron of 20 cm. According to Oishi et al. thesis, the neutron dose equivalent rate decreased to about 1/100 in the front and

Progress in Nuclear Science and Technology, Volume 4, 2014 283 the rear concrete wall at 18 MV X ray [6]. The wall was composed of the concrete of 54 cm and the iron of 66 cm. The decrease in the dose rate was almost the same though our case was about a half in the thickness of the shielding. Oishi et al. calculates the attenuation of the neutron and photon dose rate in the shielding wall to the incidence of the photon of 18 MV. It agreed in the point that the rise of the neutron dose in front side of iron, and the sudden attenuation in a rear concrete block. The neutron dose equivalent rate in the front and the rear shielding is thought that the difference of the cross section of an photonuclear reaction by the energy of entering X rays influences. The neutron dose measured in present work was 2.0 10-5 Sv/MU in the linac, and 0.4 10-5 Sv/MU in the linac at the position of the isocenter. Radiation Shielding Information Center Code Package [12] shows 3.0 10-5 Sv/MU in the and 0.3 10-5 Sv/MU in the 10MV, it agreed almost. This kernel integration code system for multi- group gamma-ray scattering calculations using the GP buildup factor was contributed by Japan Atomic Energy Research Institute which based on CCC-75/G3-6ED, contributed by Los Alamos National Laboratory. The Q values ranged from 0.06 for a beam to 0.76 ( 10 12 ) neutrons per photon Gy for a beam of Varian. They measured number of fast neutron fluence using the gold-foil activation technique [13]. 3.2. Fast and thermal neutron behaviors inside a model of shielding wall Figures 4 shows the number of etch pits of WNP produced by fast neutrons inside a model of shielding wall. The etch pits are decreased exponentially at the model of concrete-only shielding wall except deepest position of. The reason of this phenomenon is thought that the photoneutrons generated at the target and the linac head entered through the back surface of the shield due to the room scattering. It is thought that the scattering radiation can be prevented from contributing by arranging the fast neutron moderator with polyethylene and the thermal neutron absorption material with boron behind the concrete brocks. The attenuation of the fast neutron was different in 10 MV and. The reason is thought that the neutron energy by the model is different. According to Howell et. al [14], the neutron average energies of linacs are 0.27 MeV at Varian, 0.20 MeV at Siemens, and 0.23 MeV at Elekta. Energy were almost the same though the manufactures were different. According to Kim et. al [15], they found the neutron average energy at isocenter for a 0 0 cm 2 closed field is 0.38 MeV for and is 0.45 MeV for. The results for the concrete iron shielding wall showed the energy dependence of the incident photons. On the basis of this information, we conclude that a photonuclear reaction occurred inside the iron shielding wall, in which fast neutrons were generated by incident photons mainly from the target. Figure 5 shows the number of etch pits of WNP produced by the thermal neutrons are decreased exponentially in the concrete-only shield. The decrease of thermal neutrons in iron slabs was remarkable, however, they increased again in back concrete part. The reason can be thought that the fast neutrons permeated the rear concrete slab were moderated to thermal neutrons, and the thermal neutrons permeated backward the iron slab again. And thermal neutron scatters in the room and incident to the rear slab, because the space exists during the blocks where the CR-39s were arranged in this experiment system. 4. Conclusion In the linac, the dose contribution of secondary neutrons to the total dose was negligibly small. The number of WNP etch pits produced by fast neutrons increased in and near the iron slab with the photoneutron production, however, that produced by thermal neutrons decreased. References [1] M.M. Elsalim, Characterization of the neutron environment around Varian medical electron accelerators, M.S. thesis, San Jose State University, (1994). [2] Y. Uwamino, T. Nakamura, T. Ohkubo and A. Hara, Measurement and calculation of neutron leakage from a medical electron accelerator, Medical Physics 13 (1986), pp.374-384. [3] T. Nakamura and Y. Uwamino, Shielding design and detection of neutrons from medical and industrial electron accelerators--simple method of design calculation for neutron shielding, Radioisotopes 35 (1986), pp.51-56. [4] P. H. McGinley, Shielding Techniques for Radiation Oncology Facilities, Medical Physics Pub Corp, 2 edition, (2002). [5] K. Masumoto, A. Toyoda, K. Eda and T. Ishihara, Measurement of the Spatial Distribution of Neutrons in an Accelerator Room by the Combination of Activation Detectors and an Imaging Plate. Radiation Safety Management 1 (2002), pp.12-16. [6] K. Oishi, K. Kosako, Y.Kobayashi, I. Sonoki and T. Nakamura, Consideration to Shielding Design Calculation of Radiation Oncology Electron Linac Facility whose Energy is Greater than 10MeV. Radioisotopes 57 (2008), pp.277-286. [7] K. Oishi, K. Kosako, T. Nakamura, Y. Kobayashi, I. Sonoki, K. Sato, T. Kamiyama and Y. Kiyanagi: Measurement and Analysis of Dose Rate Distribution in the Laminated Shielding Assembly Irradiated by High Energy Photons, J. Nucl. Sci. and Technol., Supplement 5 (2008), pp.32-35. [8] H. R. Vega-Carrillo, S. A. Martı nez-ovalle, A. M. Lallena, G. A. Mercado and J. L. Benites-Rengifo,

284 T. Fujibuchi et al. Neutron and photon spectrain LINACs, Applied Radiation and Isotopes, (2012). [9] J. C. Evans, E. L. Lepel, R. W. Sanders, C. L. Wilkerson, W. Silker, C. W. Thomas, K. H. Abel and D. R Robertson, Long-lived activation products in reactor materials, NUREG/CR-3474, (1984). [10] N. Juto, The large scale personal monitoring service using the latest personal monitor glass badge, paper for proceedings of AOCRP-1, (2002). [11] H. Ohguchi, and T. Nakamura, Development of Wide-Energy Range Personal Neutron Dosimeter Using CR-39 Track Detector, Appl. Radiat. Isot. 46 (1995), pp.509-510. [12] Radiation Shielding Information Center Code Package, CCC-494/G33-GP, (1986). [13] D.S. Followill, M.S. Stovall and S.F. Kry, Neutron source strength measurements for Varian, Siemens, Elekta, and General Electric linear accelerators, Journal of Applied Clinical Medical Physics 4(3) (2003), pp.189-194. [14] R. M. Howell, S. F. Kry, E. Burgett, N. E. Hertel and D. S. Followill, Secondary neutron spectra from modern Varian, Siemens, and Elekta linacs with multileaf collimators, Medical Physics 36 (2009), pp.4027-4038. [15] H. S. Kim, Y.H. Park, B. C. Koo, J. W. Kwon, J. S. Lee and H. S. Choi, Evaluation of the photoneutron field produced in a medical linear accelerator, Radiat Prot Dosimetry 123(3), (2007), pp.323-328.