Po-production in lead: A benchmark between Geant4, FLUKA and MCNPX

On the last SATIF a comparison between the measured activities of the polonium isotopes Po-208, Po-209 and Po-210 and the simulated results using MCNPX2.7.0 was presented. The lead samples were cut from the SINQ spallation target at the Paul Scherrer Institut (PSI) and irradiated in 2000/2001 by 575 MeV protons. The Po-isotopes were separated using radiochemical methods by the group of D. Schumann at PSI and measured. Choosing the default model in MCNPX, Bertini-Dresner, the prediction underestimated the measured activities by up to several orders of magnitude. Therefore the Li`ege intranuclear-cascade model (INCL4.6) coupled to the de-excitation model ABLA07 were implemented into MCNPX2.7.0 and very good agreement was found to the measurement. The reason for the disagreement was traced back to a suppression of alpha reactions on the lead isotopes leading to Po and neglecting the triton capture on Pb-208, which leads to Pb-210 and decays into Po-210 with a much longer life time (22.3 years) than the decay of Po-210 itself (138 days). The prediction of the Po-isotope activities turns out to be a sensitive test for models and codes as it requires the accurate treatment of reaction channels not only with neutrons, protons and pions but also with alphas and tritons, which are not considered in intra-nuclear cascade models of the first generation. Therefore it was decided to perform a benchmark by comparing the results obtained with MCNPX2.7.0 using INCL4.6/ABLA07 to the predictions of FLUKA and Geant4. Since the model of the SINQ spallation source requires an elaborate geometry a toy model was setup. The toy model has a simplified geometry preserving the main features of the original geometry. The results for the activities of the three Po-isotopes and Pb-210 as well as the energy spectra for alphas and tritons obtained with the three particle transport Monte Carlo codes are presented.