Premium
Monte Carlo calculations of thermal neutron capture in gadolinium: A comparison of GEANT4 and MCNP with measurements
Author(s) -
Enger Shirin A.,
Munck af Rosenschöld Per,
Rezaei Arash,
Lundqvist Hans
Publication year - 2006
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1118/1.2150787
Subject(s) - neutron temperature , neutron , monte carlo method , gadolinium , neutron scattering , nuclear physics , scattering , fluence , physics , materials science , neutron transport , neutron capture , radiochemistry , irradiation , chemistry , optics , statistics , mathematics , metallurgy
GEANT4 is a Monte Carlo code originally implemented for high‐energy physics applications and is well known for particle transport at high energies. The capacity of GEANT4 to simulate neutron transport in the thermal energy region is not equally well known. The aim of this article is to compare MCNP , a code commonly used in low energy neutron transport calculations and GEANT4 with experimental results and select the suitable code for gadolinium neutron capture applications. To account for the thermal neutron scattering from chemically bound atoms [ S ( α , β ) ] in biological materials a comparison of thermal neutron fluence in tissue‐like poly(methylmethacrylate) phantom is made with MCNP4B , GEANT4 6.0 patch1, and measurements from the neutron capture therapy (NCT) facility at the Studsvik, Sweden. The fluence measurements agreed with MCNP calculated results considering S ( α , β ) . The location of the thermal neutron peak calculated with MCNP without S ( α , β )and GEANT4 is shifted by about 0.5 cm towards a shallower depth and is 25%–30% lower in amplitude. Dose distribution from the gadolinium neutron capture reaction is then simulated by MCNP and compared with measured data. The simulations made by MCNP agree well with experimental results. As long as thermal neutron scattering from chemically bound atoms are not included in GEANT4 it is not suitable for NCT applications.