Simulation and Analytical Solution

Organic scintillators, in both liquid and plastic form, are commonly used in systems for the detection of nuclear materials in nonproliferation and homeland security applications. Neutron detection in this type of detector occurs by multiple scatterings on hydrogen (H) and carbon (C), the main constituents of the scintillator. The analysis of the statistics of neutron collisions is important to understand the mechanism of neutron detection, and to perform subsequent unfolding procedures aimed at determining the incident neutron spectrum. We describe the Monte Carlo simulation of neutron interactions with the detector material for varying detector sizes and varying incident neutron energies. The simulations are performed using the code MCNP-PoliMi, which allows event-resolved predictions of the interactions of neutrons with the detector material. A subsequent post-processing of the simulation results allows us to determine the number of elastic collisions that the neutrons undergo with the nuclei of hydrogen and carbon atoms, together with the amount of energy that is deposited as a function of the number of collisions. The ligh