Monte Carlo study of Siemens PRIMUS photoneutron production

The dose from photon-induced nuclear particles (neutrons, protons, and alpha particles) generated by high-energy photon beams from medical linacs is investigated. Monte Carlo calculations using the MCNPX code are performed for three different photon beams from two different machines: Siemens 18 MV, Varian 15 MV, and Varian 18 MV. The linac head components are simulated in detail. The dose distributions from photons, neutrons, protons, and alpha particles are calculated in a tissue-equivalent phantom. Neutrons are generated in both the linac head and the phantom. This study includes (a) field size effects, (b) off-axis dose profiles, (c) neutron contribution from the linac head, (d) dose contribution from capture gamma rays, (e) phantom heterogeneity effects, and (f) effects of primary electron energy shift. Results are presented in terms of absolute dose distributions and also in terms of DER (dose equivalent ratio). The DER is the maximum dose from the particle (neutron, proton, or alpha) divided by the maximum photon dose, multiplied by the particle quality factor and the modulation scaling factor. The total DER including neutrons, protons, and alphas is about 0.66 cSv/Gy for the Siemens 18 MV beam (10 cm x 10 cm). The neutron DER decreases with decreasing field size while the proton (or alpha) DER does not vary significantly except for the 1 cm x 1 cm field. Both Varian beams (15 and 18 MV) produce more neutrons, protons, and alphas particles than the Siemens 18 MV beam. This is mainly due to their higher primary electron energies: 15 and 18.3 MeV, respectively, vs 14 MeV for the Siemens 18 MV beam. For all beams, neutrons contribute more than 75% of the total DER, except for the 1 cm x 1 cm field (approximately 50%). The total DER is 1.52 and 2.86 cSv/Gy for the 15 and 18 MV Varian beams (10 cm x 10 cm), respectively. Media with relatively high-Z elements like bone may increase the dose from heavy charged particles by a factor 4. The total DER is sensitive to primary electron energy shift. A Siemens 18 MV beam with 15 MeV (instead of 14 MeV) primary electrons would increase by 40% the neutron DER and by 210% the proton + alpha DER. Comparisons with measurements (neutron yields from different materials and neutron dose equivalent) are also presented. Using the NCRP risk assessment method, we found that the dose equivalent from leakage neutrons (at 50-cm off-axis distance) represent 1.1, 1.1, and 2.0% likelihood of fatal secondary cancer for a 70 Gy treatment delivered by the Siemens 18 MV, Varian 15 MV, and Varian 18 MV beams, respectively.

The shielding calculations for high energy (>10 MV) linear accelerators must include the photoneutron production within the head of the accelerator. Procedures have been described to calculate the treatment room door shielding based on the neutron source strength (Q value) for a specific accelerator and energy combination. Unfortunately, there is currently little data in the literature stating the neutron source strengths for the most widely used linear accelerators. In this study, the neutron fluence for 36 linear accelerators, including models from Varian, Siemens, Elekta/Philips, and General Electric, was measured using gold-foil activation. Several of the models and energy combinations had multiple measurements. The neutron fluence measured in the patient plane was independent of the surface area of the room, suggesting that neutron fluence is more dependent on the direct neutron fluence from the head of the accelerator than from room scatter. Neutron source strength, Q, was determined from the measured neutron fluences. As expected, Q increased with increasing photon energy. The Q values ranged from 0.02 for a 10 MV beam to 1.44(x10(12)) neutrons per photon Gy for a 25 MV beam. The most comprehensive set of neutron source strength values, Q, for the current accelerators in clinical use are presented for use in calculating room shielding.

A combination of three superheated drop detectors with different neutron energy responses was developed to evaluate dose-equivalent and energy distributions of photoneutrons in a phantom irradiated by radiotherapy high-energy x-ray beams. One of the three detectors measures the total neutron dose equivalent and the other two measure the contributions from fast neutrons above 1 and 5.5 MeV, respectively. In order to test the new method, the neutron field produced by the 10 cm X 10 cm x-ray beam of an 18 MV radiotherapy accelerator was studied. Measurements were performed inside a tissue-equivalent liquid phantom, at depths of 1, 5, 10 and 15 cm and at lateral distances of 0, 10, and 20 cm from the central axis. These data were used to calculate the average integral dose to the radiotherapy patient from direct neutrons as well as from neutrons transmitted through the accelerator head. The characteristics of the dosimeters were confirmed by results in excellent agreement with those of prior studies. Track etch detectors were also used and provided an independent verification of the validity of this new technique. Within the primary beam, we measured a neutron entrance dose equivalent of 4.5 mSv per Gy of photons. It was observed that fast neutrons above 1 MeV deliver most of the total neutron dose along the beam axis. Their relative contribution increases with depth, from about 60% at the entrance to over 90% at a depth of 10 cm. Thus, the average energy increases with depth in the phantom as neutron spectra harden.