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

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 Ca and 60 Co is studied as an example. In addition, shifts in the stable isotope ratios 157 Gd/ 160 Gd, 158 Gd/ 160 Gd, 149 Sm/ 152 Sm, and 150 Sm/ 152 Sm (and combinations thereof) are modeled and compared to experimental data. The model describes experimental 41 Ca activity concentrations in different types of meteorites and the lunar surface within the uncertainties. In contrast, it fails to describe 60 Co activity concentrations. In addition, it is difficult to consistently model the isotope shifts 157 Gd/ 160 Gd and 150 Sm/ 152 Sm in Apollo 15 drill core samples. The observed trends depend on the temperature of the irradiated object and are more pronounced for colder temperatures. Since the observed discrepancies are likely related to the shape of the neutron spectra, self-shielding effects by, e.g., 56 Fe, might be of importance and some of the consequences are discussed.