• ### The hadronic interaction model SIBYLL 2.3c and Feynman scaling(1709.07227)

Sept. 21, 2017 hep-ph
The Monte Carlo model Sibyll has been designed for efficient simulation of hadronic multiparticle production up to the highest energies as needed for interpreting cosmic ray measurements. For more than 15 years, version 2.1 of Sibyll has been one of the standard models for air shower simulation. Motivated by data of LHC and fixed-target experiments and a better understanding of the phenomenology of hadronic interactions, we have developed an improved version of this model, version 2.3, which has been released in 2016. In this contribution we present a revised version of this model, called Sibyll 2.3c, that is further improved by adjusting particle production spectra to match the expectation of Feynman scaling in the fragmentation region. After a brief introduction to the changes implemented in Sibyll 2.3 and 2.3c with respect to Sibyll 2.1, the current predictions of the model for the depth of shower maximum, the number of muons at ground, and the energy spectrum of muons in extensive air showers are presented.
• We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.
• Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
• A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.
• ### The FUNK search for Hidden Photon Dark Matter in the eV range(1510.05869)

Oct. 20, 2015 hep-ph, hep-ex, physics.ins-det
We give a brief update on the search for Hidden Photon Dark Matter with FUNK. The experiment uses a large spherical mirror, which, if Hidden Photon Dark Matter exists in the accessible mass and coupling parameter range, would yield an optical signal in the mirror's center in an otherwise dark environment. After a test run with a CCD, preparations for a run with a low-noise PMT are under way and described in this proceedings.
• ### A new version of the event generator Sibyll(1510.00568)

Oct. 2, 2015 hep-ph
The event generator Sibyll can be used for the simulation of hadronic multiparticle production up to the highest cosmic ray energies. It is optimized for providing an economic description of those aspects of the expected hadronic final states that are needed for the calculation of air showers and atmospheric lepton fluxes. New measurements from fixed target and collider experiments, in particular those at LHC, allow us to test the predictive power of the model version 2.1, which was released more than 10 years ago, and also to identify shortcomings. Based on a detailed comparison of the model predictions with the new data we revisit model assumptions and approximations to obtain an improved version of the interaction model. In addition a phenomenological model for the production of charm particles is implemented as needed for the calculation of prompt lepton fluxes in the energy range of the astrophysical neutrinos recently discovered by IceCube. After giving an overview of the new ideas implemented in Sibyll and discussing how they lead to an improved description of accelerator data, predictions for air showers and atmospheric lepton fluxes are presented.
• ### The ultra-high energy cosmic rays image of Virgo A(1509.09033)

Sept. 30, 2015 astro-ph.HE
Arrival directions of ultra-high energy cosmic rays from the direction of ten brightest radio sources lying within 50 Mpc from our Galaxy were studied by using recent models of the largescale Galactic magnetic field. A detailed study, where also small-scale turbulent magnetic field component was implemented, is presented for the radiogalaxy Virgo A. This radiogalaxy is located far from the Galactic plane which leads to a unique image of this UHECR source candidate, if the flux is composed from a mixture of intermediate mass nuclei. We present a method suitable for identifying cosmic rays arriving from this close-by radiogalaxy.
• ### Calculation of conventional and prompt lepton fluxes at very high energy(1503.00544)

March 3, 2015 hep-ph, astro-ph.HE
An efficient method for calculating inclusive conventional and prompt atmospheric leptons fluxes is presented. The coupled cascade equations are solved numerically by formulating them as matrix equation. The presented approach is very flexible and allows the use of different hadronic interaction models, realistic parametrizations of the primary cosmic-ray flux and the Earth's atmosphere, and a detailed treatment of particle interactions and decays. The power of the developed method is illustrated by calculating lepton flux predictions for a number of different scenarios.
• ### Charm production in SIBYLL(1502.06353)

Feb. 23, 2015 hep-ph
SIBYLL 2.1 is an event generator for hadron interactions at the highest energies. It is commonly used to analyze and interpret extensive air shower measurements. In light of the first detection of PeV neutrinos by the IceCube collaboration the inclusive fluxes of muons and neutrinos in the atmosphere have become very important. Predicting these fluxes requires understanding of the hadronic production of charmed particles since these contribute significantly to the fluxes at high energy through their prompt decay. We will present an updated version of SIBYLL that has been tuned to describe LHC data and extended to include the production of charmed hadrons.
• ### Hidden Photon Dark Matter Search with a Large Metallic Mirror(1410.0200)

Oct. 1, 2014 hep-ph, hep-ex, physics.ins-det
If Dark Matter is composed of hidden-sector photons that kinetically mix with photons of the visible sector, then Dark Matter has a tiny oscillating electric field component. Its presence would lead to a small amount of visible radiation being emitted from a conducting surface, with the photon frequency given approximately by the mass of the hidden photon. Here, we report on experimental efforts that have started recently to search for such hidden photon Dark Matter in the (sub-)eV regime with a prototype mirror for the Auger fluorescence detector at the Karlsruhe Institute for Technology.
• ### Asymmetry of the angular distribution of Cherenkov photons of extensive air showers induced by the geomagnetic field(1405.4671)

May 19, 2014 astro-ph.IM, astro-ph.HE
The angular distribution of Cherenkov light in an air shower is closely linked to that of the shower electrons and positrons. As charged particles in extensive air showers are deflected by the magnetic field of the Earth, a deformation of the angular distribution of the Cherenkov light, that would be approximately symmetric about the shower axis if no magnetic field were present, is expected. In this work we study the variation of the Cherenkov light distribution as a function of the azimuth angle in the plane perpendicular to shower axis. It is found that the asymmetry induced by the geomagnetic field is most significant for early stages of shower evolution and for showers arriving almost perpendicular to the vector of the local geomagnetic field. Furthermore, it is shown that ignoring the azimuthal asymmetry of Cherenkov light might lead to a significant under- or overestimation of the Cherenkov light signal especially at sites where the local geomagnetic field is strong. Based on CORSIKA simulations, the azimuthal distribution of Cherenkov light is parametrized in dependence on the magnetic field component perpendicular to the shower axis and the local air density. This parametrization provides an efficient approximation for estimating the asymmetry of the Cherenkov light distribution for shower simulation and reconstruction in cosmic ray and gamma-ray experiments in which the Cherenkov signal of showers with energies above 10^14 eV is observed.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### Simulation of ultra-high energy photon propagation with PRESHOWER 2.0(1302.6939)

Feb. 27, 2013 astro-ph.IM, astro-ph.HE
In this paper we describe a new release of the PRESHOWER program, a tool for Monte Carlo simulation of propagation of ultra-high energy photons in the magnetic field of the Earth. The PRESHOWER program is designed to calculate magnetic pair production and bremsstrahlung and should be used together with other programs to simulate extensive air showers induced by photons. The main new features of the PRESHOWER code include a much faster algorithm applied in the procedures of simulating the processes of gamma conversion and bremsstrahlung, update of the geomagnetic field model, and a minor correction. The new simulation procedure increases the flexibility of the code so that it can also be applied to other magnetic field configurations such as, for example, encountered in the vicinity of the sun or neutron stars.
• ### Disentangling the Air Shower Components Using Scintillation and Water Cherenkov Detectors(1208.2167)

Aug. 10, 2012 astro-ph.HE
We consider a ground array of scintillation and water Cherenkov detectors with the purpose of determining the muon content of air showers. The different response characteristics of these two types of detectors to the components of the air shower provide a way to infer their relative contributions. We use a detailed simulation to estimate the impact of parameters, such as scintillation detector size, in the determination of the size of the muon component.
• ### Anomalous Longitudinal Shower Profiles and Hadronic Interactions(1111.0504)

Nov. 2, 2011 astro-ph.HE
The bulk of air showers initiated by very high energy cosmic rays exhibits a longitudinal development in depth with a single well-defined shower maximum. However, a small fraction of showers has a profile that differs considerably from this average behaviour. In extreme cases, such anomalous longitudinal profiles can even have two distinct shower maxima. We discuss the properties of the primary interactions that lead to such profiles. Simulations are used to estimate the rate of anomalous profiles in dependence of primary energy, mass, and characteristic features of hadronic multiparticle production at very high energies.
• ### Mass Composition Sensitivity of an Array of Water Cherenkov and Scintillation Detectors(1109.1692)

Sept. 8, 2011 astro-ph.IM, astro-ph.HE
We consider a hybrid array composed of scintillation and water Cherenkov detectors designed to measure the cosmic ray primary mass composition at energies of about 1 EeV. We have developed a simulation and reconstruction chain to study the theoretical performance of such an array. In this work we investigate the sensitivity of mass composition observables in relation to the geometry of the array. The detectors are arranged in a triangular grid with fixed 750 m spacing and the configuration of the scintillator detectors is optimized for mass composition sensitivity. We show that the performance for composition determination can be compared favorably to that of Xmax measurements after the difference in duty cycles is considered.
• ### The strong interaction at the collider and cosmic-rays frontiers(1106.2453)

June 13, 2011 hep-ph, hep-ex, astro-ph.HE
First data on inclusive particle production measured in proton-proton collisions at the Large Hadron Collider (LHC) are compared to predictions of various hadron-interaction Monte Carlos (QGSJET, EPOS and SIBYLL) used commonly in high-energy cosmic-ray physics. While reasonable overall agreement is found for some of the models, none of them reproduces consistently the sqrt(s) evolution of all the measured observables. We discuss the implications of the new LHC data for the modeling of the non-perturbative and semihard parton dynamics in hadron-hadron and cosmic-rays interactions at the highest energies studied today.
• ### Constraints from the first LHC data on hadronic event generators for ultra-high energy cosmic-ray physics(1101.5596)

May 4, 2011 hep-ph, hep-ex, astro-ph.HE
The determination of the primary energy and mass of ultra-high-energy cosmic-rays (UHECR) generating extensive air-showers in the Earth's atmosphere, relies on the detailed modeling of hadronic multiparticle production at center-of-mass (c.m.) collision energies up to two orders of magnitude higher than those studied at particle colliders. The first Large Hadron Collider (LHC) data have extended by more than a factor of three the c.m. energies in which we have direct proton-proton measurements available to compare to hadronic models. In this work we compare LHC results on inclusive particle production at energies sqrt(s) = 0.9, 2.36, and 7 TeV to predictions of various hadronic Monte Carlo (MC) models used commonly in cosmic-ray (CR) physics (QGSJET, EPOS and SIBYLL). As a benchmark with a standard collider physics model we also show PYTHIA (and PHOJET) predictions with various parameter settings. While reasonable overall agreement is found for some of the MC, none of them reproduces consistently the sqrt(s) evolution of all the observables. We discuss implications of the new LHC data for the description of cosmic-ray interactions at the highest energies.
• ### Sibyll with charm(1102.5705)

Feb. 28, 2011 hep-ph, hep-ex, astro-ph.HE
The cosmic ray interaction event generator Sibyll is widely used in extensive air shower simulations for cosmic ray and neutrino experiments. Charmed particle production has been added to the Monte Carlo with a phenomenological, non-perturbative model that properly accounts for charm production in the forward direction. As prompt decays of charm can become a significant background for neutrino detection, proper simulation of charmed particles is very important. We compare charmed meson and baryon production to accelerator data.
• ### Hadronic Multiparticle Production at Ultra-High Energies and Extensive Air Showers(1010.4310)

Oct. 20, 2010 hep-ph, astro-ph.HE
Studies of the nature of cosmic ray particles at the highest energies are based on the measurement of extensive air showers. Most cosmic ray properties can therefore only be obtained from the interpretation of air shower data and are thus depending on predictions of hadronic interaction models at ultra-high energies. We discuss different scenarios of model extrapolations from accelerator data to air shower energies and investigate their impact on the corresponding air shower predictions. To explore the effect of different extrapolations by hadronic interaction models we developed an ad hoc model. This ad hoc model is based on the modification of the output of standard hadronic interaction event generators within the air shower simulation process and allows us to study the impact of changing interaction features on the air shower development. In a systematic study we demonstrate the resulting changes of important air shower observables and also discuss them in terms of the predictions of the Heitler model of air shower cascades. It is found that the results of our ad hoc modifications are, to a large extend, independent of the choice of the underlying hadronic interaction model.
• ### Cosmic ray event generator Sibyll 2.1(0906.4113)

Nov. 30, 2009 hep-ph, hep-ex, astro-ph.HE
The cosmic ray interaction event generator Sibyll is widely used in extensive air shower simulations. We describe in detail the properties of Sibyll 2.1 and the differences with the original version 1.7. The major structural improvements are the possibility to have multiple soft interactions, introduction of new parton density functions, and an improved treatment of diffraction. Sibyll 2.1 gives better agreement with fixed target and collider data, especially for the inelastic cross sections and multiplicities of secondary particles. Shortcomings and suggestions for future improvements are also discussed.
• ### Proton-Air Cross Section and Extensive Air Showers(0906.3075)

June 17, 2009 hep-ph, astro-ph.HE
Hadronic cross sections at ultra-high energy have a significant impact on the development of extensive air shower cascades. Therefore the interpretation of air shower data depends critically on hadronic interaction models that extrapolate the cross section from accelerator measurements to the highest cosmic ray energies. We discuss how extreme scenarios of cross section extrapolations can affect the interpretation of air shower data. We find that the theoretical uncertainty of the extrapolated proton-air cross section at ultra-high energies is much larger than suggested by the existing spread of available Monte Carlo model predictions. The impact on the depth of the shower maximum is demonstrated.
• ### Sensitivity of Extensive Air Showers to Features of Hadronic Interactions at Ultra-High Energies(0906.0418)

June 2, 2009 astro-ph.HE
We study the dependence of extensive air shower development on the first hadronic interactions at ultra-high energies occurring in the startup phase of the air shower cascade. The interpretation of standard air shower observables depends on the characteristics of these interactions. Thus, it is currently difficult to draw firm conclusions for example on the primary cosmic ray mass composition from the analysis of air shower data. On the other hand, a known primary mass composition would allow us to study hadronic interactions at center of mass energies well above the range that is accessible to accelerators measurements.
• ### On the measurement of the proton-air cross section using air shower data(0903.0404)

May 14, 2009 hep-ph, astro-ph.HE
The analysis of high-energy air shower data allows one to study the proton-air cross section at energies beyond the reach of fixed target and collider experiments. The mean depth of the first interaction point and its fluctuations are a measure of the proton-air particle production cross section. Since the first interaction point in air cannot be measured directly, various methods have been developed in the past to estimate the depth of the first interaction from air shower observables in combination with simulations. As the simulations depend on assumptions made for hadronic particle production at energies and phase space regions not accessible in accelerator experiments, the derived cross sections are subject to significant systematic uncertainties. The focus of this work is the development of an improved analysis technique that allows a significant reduction of the model dependence of the derived cross section at very high energy. Performing a detailed Monte Carlo study of the potential and the limitations of different measurement methods, we quantify the dependence of the measured cross section on the used hadronic interaction model. Based on these results, a general improvement to the analysis methods is proposed by introducing the actually derived cross section already in the simulation of reference showers. The reduction of the model dependence is demonstrated for one of the measurement methods.