Abstract: In this work, we discuss the possibility to test Heavy Neutral Leptons (HNLs) using “Cosmic Ray Beam Dump” experiments. In analogy with terrestrial beam dump experiments, where a beam first hits a target and is then absorbed by a shield, we consider high-energy incident cosmic rays impinging on the Earth's atmosphere and then the Earth's surface. We focus here on HNL production from atmospherically produced kaon, pion and D-meson decays, and discuss the possible explanation of the appearing Cherenkov showers observed by the SHALON Cherenkov telescope and the ultra-high energy events detected by the neutrino experiment ANITA. We show that these observations can not be explained with a long-lived HNL, as the relevant parameter space is excluded by existing constraints. Then we propose two new experimental setups that are inspired by these experiments, namely a Cherenkov telescope pointing at a sub-horizontal angle and shielded by the mountain cliff at Mount Thor, and a geostationary satellite that observes part of the Sahara desert. We show that the Cherenkov telescope at Mount Thor can probe currently untested HNL parameter space for masses below the kaon mass. We also show that the geostationary satellite experiment can significantly increase the HNL parameter space coverage in the whole mass range from 10 MeV up to 2 GeV and test neutrino mixing |U-& alpha;4|(2) down to 10(-11) for masses around 300 MeV.

Abstract: Observations of cosmic rays arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three method can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30, ... , 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.

Abstract: Research in many areas of modern physics such as, e. g., indirect searches for dark matter and particle acceleration in supernova remnant shocks rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma-rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions. The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, which uses astrophysical information, and nuclear and particle data as inputs to self-consistently predict CRs, gamma-rays, synchrotron, and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.

Abstract: The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs “radio-hybrid” measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluorescence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.

Abstract: The existence of diffuse Galactic neutrino production is expected from cosmic-ray interactions with Galactic gas and radiation fields. Thus, neutrinos are a unique messenger offering the opportunity to test the products of Galactic cosmic-ray interactions up to energies of hundreds of TeV. Here we present a search for this production using ten years of Astronomy with a Neutrino Telescope and Abyss environmental RESearch (ANTARES) track and shower data, as well as seven years of IceCube track data. The data are combined into a joint likelihood test for neutrino emission according to the KRA(gamma) model assuming a 5 PeV per nucleon Galactic cosmic-ray cutoff. No significant excess is found. As a consequence, the limits presented in this Letter start constraining the model parameter space for Galactic cosmic-ray production and transport.