Model calculations for neutron-induced reactions in meteorites and planetary surfaces

Abstract This study investigates the capabilities of the GEANT4 Monte Carlo toolkit to quantitatively predict neutron production, neutron transport, and nuclide production by neutron capture reactions in cosmochemical relevant objects. The model reproduces neutron densities measured in the lunar surface within the experimental uncertainties, which is a major improvement compared to earlier studies. Since, for many applications in meteorites and planetary surfaces, nuclide production by neutron capture is of importance, the production of $$^{41}$$ 41 Ca and $$^{60}$$ 60 Co is studied as an example. In addition, shifts in the stable isotope ratios $$^{157}$$ 157 Gd/ $$^{160}$$ 160 Gd, $$^{158}$$ 158 Gd/ $$^{160}$$ 160 Gd, $$^{149}$$ 149 Sm/ $$^{152}$$ 152 Sm, and $$^{150}$$ 150 Sm/ $$^{152}$$ 152 Sm (and combinations thereof) are modeled and compared to experimental data. The model describes experimental $$^{41}$$ 41 Ca activity concentrations in different types of meteorites and the lunar surface within the uncertainties. In contrast, it fails to describe $$^{60}$$ 60 Co activity concentrations. In addition, it is difficult to consistently model the isotope shifts $$^{157}$$ 157 Gd/ $$^{160}$$ 160 Gd and $$^{150}$$ 150