Secondary neutrons inside a proton therapy facility: MCNPX simulations compared to measurements performed with a Bonner Sphere Spectrometer and neutron H*(10) monitors

Neutron spectrometry measurements with an extended-range Bonner Sphere Spectrometer (BSS), as well as neutron H*(10) measurements using an extended-range rem meter WENDI-2, a conventional rem meter LB 6411 and a tissue-equivalent proportional counter, were performed inside and around the Fixed-Beam Treatment Room at the proton therapy facility of Essen, in Germany. The WENDI-2 stood out as the easiest detector for making accurate neutron H*(10) measurements, since its direct measurements were equivalent to the H*(10) rates obtained with the BSS. The measurements were also compared to simulation results obtained with MCNPX 2.7.0 using two different selections of physics models for the hadron interactions above 150 MeV: the Bertini & Dresner models and the CEM03 model. For neutron H*(10) rates outside the treatment room, factors of 1.6–1.8 were obtained between the results of the two simulations, the Bertini & Dresner models yielding the largest values in all positions. The comparison of the simulation results with the WENDI-2 and BSS measurements for positions inside the treatment room showed that the Bertini & Dresner models reproduce the global neutron production in the water phantom relatively well, whereas the CEM03 model underestimates it by a factor of ∼1.3. At the most-forward angle, however, the Bertini model (unlike the CEM03 model) seemed to overestimate the production of neutrons with energies above 100 MeV. Outside the shielding, the simulated H*(10) overestimated the WENDI-2 measurements by factors of 2–3 with the Bertini & Dresner models, and 1.1–1.7 with the CEM03 model. Both simulations were thus conservative with respect to the neutron fluxes transmitted through the concrete walls. This conservative behaviour is probably caused by a combination of several uncertainties, including for instance uncertainties on the proton and neutron interaction cross-sections and uncertainties on the concrete composition and density.