Monte Carlo simulations of elemental imaging using the neutron-associated particle technique.

Abstract

The purpose of this study is to develop and employ a Monte Carlo (MC) simulation model of associated particle neutron elemental imaging (APNEI) in order to determine the three-dimensional (3D) imaging resolution of such a system by examining relevant physical and technological parameters and to thereby begin to explore the range of clinical applicability of APNEI to fields such as medical diagnostics, intervention, and etiological research. The presented APNEI model was defined in MCNP by a Gaussian-distributed and isotropic surface source emitting deuterium + deuterium (DD) neutrons, iron as the target element, nine iron-containing voxels (1 cm3 volume each) arranged in a 3-by-3 array as the interrogated volume of interest, and finally, by high-purity germanium (HPGe) gamma-ray detectors anterior and posterior to the 9-voxel array. The MCNP f8 pulse height tally was employed in conjunction with the PTRAC particle tracking function to not only determine the signal acquired from iron inelastic scatter gamma-rays but also to quantitate each of the nine target voxels' contribution to the overall iron signal - each detected iron inelastic scatter gamma-ray being traced to the source neutron which incited its emission. With the spatial, vector, and timing information of the series of events for each relevant neutron history as collected by PTRAC, realistic grayscale images of the distribution of iron concentration in the 9-voxel array were simulated in both the projective and depth dimensions. With an overall 225 ps timing resolution, 6.25 mm2 imaging plate pixels assumed to have well localized scintillation, and a DD neutron, Gaussian-distributed source spot with a diameter of 2 mm, projective and depth resolutions of < 1 cm and <3 cm are achievable, respectively, for iron-containing voxels on the order of 1,000 ppm Fe. The imaging resolution offered by APNEI of target elements such as iron lends itself to potential applications in disease diagnosis and treatment planning (high resolution) as well as to ordnance and contraband detection (low resolution). However, experimental study beyond simulation is required to optimize the layout and electronic configuration of APNEI system components - including realistic shielding and phantom materials - for background signal reduction in order to accurately determine the detection limits and spatial resolution of iron and other elements of interest on a case-by-case basis.