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.Comment: 15 pages, 11 figures, Presented paper at the 13th Meeting of the task-force on Shielding aspects of Accelerators, Targets and Irradiation Facilities (SATIF-13), HZDR, October 10-12, 2016, Dresden, German

Monte Carlo calculations for the ATLAS cavern background

A new application for simulating the ATLAS cavern background was developed. This was done using FLUGG, software that allows Geant4 geometry to be used within the FLUKA simulation framework. A Geant4 description of the ATLAS detector including its cavern was built from scratch for this application. In order to gain computing performance, our geometry is less detailed than that of GeoModel which is used in the full detector simulation, but good enough for the investigation of cavern background. Our geometry can also be used in a standalone Geant4 simulation. Thus it is possible to perform unbiased comparisons between Geant4 and FLUKA using the same complex geometry. We compared neutron and photon fluxes using the FLUKA-FLUGG application with the result of Geant4 simulations based on the QGSP_BERT and QGSP_BERT_HP physics lists. In all cases the same set of initial collision 4-vectors produced by the PHOJET event generator was used. The result from the QGSP_BERT_HP physics list, which uses the High Precision (HP) neutron model, is similar to the result of FLUKA-FLUGG and the differences in the fluxes between them are within 40% in most regions of the ATLAS cavern. The result from the QGSP_BERT physics list, which does not include the HP model, does not agree with either of the previous two results

Diffractive Interaction and Scaling Violation in pp->pi^0 Interaction and GeV Excess in Galactic Diffuse Gamma-Ray Spectrum of EGRET

We present here a new calculation of the gamma-ray spectrum from pp->pi^0 in the Galactic ridge environment. The calculation includes the diffractive pp interaction and incorporates the Feynman scaling violation for the first time. Galactic diffuse gamma-rays come, predominantly, from pi^0->gamma gamma in the sub-GeV to multi-GeV range. Hunter et al. found, however, an excess in the GeV range ("GeV Excess") in the EGRET Galactic diffuse spectrum above the prediction based on experimental pp->pi^0 cross-sections and the Feynman scaling hypothesis. We show, in this work, that the diffractive process makes the gamma-ray spectrum harder than the incident proton spectrum by ~0.05 in power-law index, and, that the scaling violation produces 30-80% more pi^0 than the scaling model for incident proton energies above 100GeV. Combination of the two can explain about a half of the "GeV Excess" with the local cosmic proton (power-law index ~2.7). The excess can be fully explained if the proton spectral index in the Galactic ridge is a little harder (~0.2 in power-law index) than the local spectrum. Given also in the paper is that the diffractive process enhances e^+ over e^- and the scaling violation gives 50-100% higher p-bar yield than without the violation, both in the multi-GeV range.Comment: 35 pages, 11 figures, to appear in Astrophysical Journa

Parameterization of Gamma, e^+/- and Neutrino Spectra Produced by p-p Interaction in Astronomical Environment

We present the yield and spectra of stable secondary particles (gamma, e^+/-, nu_e, nubar_e, nu_mu, and nubar_mu) of p-p interaction in parameterized formulae to facilitate calculations involving them in astronomical environments. The formulae are derived on the up-to-date p-p interaction model by [Kamae05] which incorporates the logarithmically rising inelastic cross-section, the diffraction dissociation process, and the Feynman scaling violation. To improve fidelity to experimental data in lower energies, two baryon resonance contributions have been added: one representing Delta(1232) and the other multiple resonances around 1600MeV/c^2. The parametrized formulae predict that all secondary particle spectra be harder by about 0.05 in power-law indices than that of the incident proton and their inclusive cross-sections be larger than those predicted by p-p interaction models based on the Feynman scaling.Comment: Errors and typos in Eqns 1-4 and Table 1 corrected and editorial changes incorporate

Search for dark matter produced in association with bottom or top quarks in √s = 13 TeV pp collisions with the ATLAS detector

A search for weakly interacting massive particle dark matter produced in association with bottom or top quarks is presented. Final states containing third-generation quarks and miss- ing transverse momentum are considered. The analysis uses 36.1 fb−1 of proton–proton collision data recorded by the ATLAS experiment at √s = 13 TeV in 2015 and 2016. No significant excess of events above the estimated backgrounds is observed. The results are in- terpreted in the framework of simplified models of spin-0 dark-matter mediators. For colour- neutral spin-0 mediators produced in association with top quarks and decaying into a pair of dark-matter particles, mediator masses below 50 GeV are excluded assuming a dark-matter candidate mass of 1 GeV and unitary couplings. For scalar and pseudoscalar mediators produced in association with bottom quarks, the search sets limits on the production cross- section of 300 times the predicted rate for mediators with masses between 10 and 50 GeV and assuming a dark-matter mass of 1 GeV and unitary coupling. Constraints on colour- charged scalar simplified models are also presented. Assuming a dark-matter particle mass of 35 GeV, mediator particles with mass below 1.1 TeV are excluded for couplings yielding a dark-matter relic density consistent with measurements