Characterization of emergent leakage neutrons from multiple layers of hydrogen/water in the lunar regolith by Monte Carlo simulation

Introduction: The leakage flux of lunar neutrons produced by precipitation of galactic cosmic ray (GCR) particles in the upper layer of the lunar regolith and measured by orbital instruments such as the Lunar Exploration Neutron Detector (LEND) is investigated by Monte Carlo simulation. Previous Monte Carlo (MC) simulations have been used to investigate neutron production and leakage from the lunar surface to assess the elemental composition of lunar soil [1-6] and its effect on the leakage neutron flux. We investigate effects on the emergent flux that depend on the physical distribution of hydrogen within the regolith. We use the software package GEANT4 [7] to calculate neutron production from spallation by GCR particles [8,9] in the lunar soil. Multiple layers of differing hydrogen/water at different depths in the lunar regolith model are introduced to examine enhancement or suppression of leakage neutron flux. We find that the majority of leakage thermal and epithermal neutrons are produced in 25 cm to 75 cm deep from the lunar surface. Neutrons produced in the shallow top layer retain more of their original energy due to fewer scattering interactions and escape from the lunar surface mostly as fast neutrons. This provides a diagnostic tool in interpreting leakage neutron flux enhancement or suppression due to hydrogen concentration distribution in lunar regolith. We also find that the emitting angular distribution of thermal and epithermal leakage neutrons can be described by cos3/2(theta) where the fast neutrons emitting angular distribution is cos(theta). The energy sensitivity and angular response of the LEND detectors SETN and CSETN are investigated using the leakage neutron spectrum from GEANT4 simulations. A simplified LRO model is used to benchmark MCNPX[10] and GEANT4 on CSETN absolute count rate corresponding to neutron flux from bombardment of 120MV solar potential GCR particles on FAN lunar soil. We are able to interpret the count rates of SETN and CSETN from the leakage neutron spectrum, emission angle, detector energy sensitivity and angular response. Reference: [1] W. C. Feldman, et al., Science 4 September 1998: Vol. 281 no. 5382 pp. 1496-1500. [2] Gasnault, O., et al., (2000) J. Geophys. Res., 105(E2), 4263-4271. [3] Little, R. C., et al. (2003), J. Geophys. Res., 108(E5), 5046. [4] McKinney et al., (2006), J. Geophys. Res., 111, E06004. [5] Lawrence et al., (2006), J. Geophys. Res., 111, E08001. [6] Looper et al, (2013), Space Weather, VOL. 11, 142-152. [7] J. Allison, et al, (2006) IEEE TRANS. ON NUCL SCI, VOL. 53, NO. 1. [8] J. Masarik and R. Reedy (1996), J. Geophys. Res., 101, 18,891-18,912. [9] P. O'Neil (2010) IEEE Trans. Nucl. Sci., 57(6), 3148-3153. [10] D. Pelowitz, (2005), Rep. LA-CP-05-0369, LANL, Los Alamos, NM.