Antideuteron and antihelium nuclei have been proposed as promising detection channels for dark matter because of the low astrophysical backgrounds expected. To estimate both potential exotic contributions and their backgrounds, one usually employs the coalescence model in momentum space. Here we use instead a newly developed coalescence model based on the Wigner function representations of the produced nuclei states. This approach includes both the process-dependent size of the formation region of antinuclei, and the momentum correlations of coalescing antinucleons in a semi-classical picture. The model contains a single universal parameter $\sigma$ that we tune to experimental data on antideuteron production in electron-positron, proton-proton and proton-nucleus collisions. The obtained value $\sigma\simeq 1$ fm agrees well with its physical interpretation as the size of the formation region of antinuclei in collisions of point-like particles. This model allows us therefore to calculate in a consistent frame-work the antideuteron and antihelium fluxes both from secondary production and from dark matter annihilations. We find that the antihelium-3 flux falls short by more than an order of magnitude of the detection sensitivity of the AMS-02 experiment, assuming standard cosmic ray propagation parameters, while the antideuteron flux can be comparable to the sensitivities of the AMS-02 and GAPS experiments.

The version 3.01 of ELMAG, a Monte Carlo program for the simulation of electromagnetic cascades initiated by high-energy photons and electrons interacting with extragalactic background light (EBL), is presented. Pair production and inverse Compton scattering on EBL photons as well as synchrotron losses are implemented using weighted sampling of the cascade development. New features include, among others, the implementation of turbulent extragalactic magnetic fields and the calculation of three-dimensional electron and positron trajectories, solving the Lorentz force equation. As final result of the three-dimensional simulations, the program provides two-dimensional source images as function of the energy and the time delay of secondary cascade particles.

Over the past 20 years, KASCADE and its extension KASCADE-Grande were dedicated to measure high-energy cosmic rays with primary energies of 100 TeV to 1 EeV. The data accumulation was fully completed and all experimental components were dismantled, though the analysis of the high-quality data is still continued. E.g., we investigated the validity of the hadronic interaction model of the new SIBYLL version 2.3c. We also published a new result of a search for large-scale anisotropies performed with the KASCADE-Grande data. Investigation of the attenuation length of the muon in the atmosphere is also updated with the predictions of the SIBYLL 2.3 interaction model. We investigated, in addition, the muon content of high-energy air showers and compared them to all post-LHC interaction models. In this contribution, the new and updated results from KASCADE-Grande will be presented. An update of the web-based data center KCDC offering the original scientific data from KASCADE-Grande to the public will be briefly discussed as well.

The `KASCADE Cosmic ray Data Centre' is a web portal (\url{https://kcdc.ikp.kit.edu}), where the data of the astroparticle physics experiment KASCADE-Grande are made available for the interested public. The KASCADE experiment was a large-area detector for the measurement of high-energy cosmic rays via the detection of extensive air showers. The multi-detector installations KASCADE and its extension KASCADE-Grande stopped the active data acquisition in 2013 of all its components end of 2012 after more than 20 years of data taking. In several updates since our first release in 2013 with KCDC we provide the public measured and reconstructed parameters of more than 433 million air showers. In addition, KCDC provides meta data information and documentation to enable a user outside the community of experts to perform their own data analysis. Simulation data from three different high energy interaction models have been made available as well as a compilation of measured and published spectra from various experiments. In addition, detailed educational examples shall encourage high-school students and early stage researchers to learn about astroparticle physics, cosmic radiation as well as the handling of Big Data and about the sustainable and public provision of scientific data.

The evolution of the muon content of very high energy air showers (EAS) in the atmosphere is investigated with data of the KASCADE-Grande observatory. For this purpose, the muon attenuation length in the atmosphere is obtained to $\Lambda_\mu = 1256 \, \pm 85 \, ^{+229}_{-232}(\mbox{syst})\, \mbox{g/cm}^2$ from the experimental data for shower energies between $10^{16.3}$ and $10^{17.0} \, \mbox{eV}$. Comparison of this quantity with predictions of the high-energy hadronic interaction models QGSJET-II-02, SIBYLL 2.1, QGSJET-II-04 and EPOS-LHC reveals that the attenuation of the muon content of measured EAS in the atmosphere is lower than predicted. Deviations are, however, less significant with the post-LHC models. The presence of such deviations seems to be related to a difference between the simulated and the measured zenith angle evolutions of the lateral muon density distributions of EAS, which also causes a discrepancy between the measured absorption lengths of the density of shower muons and the predicted ones at large distances from the EAS core. The studied deficiencies show that all four considered hadronic interaction models fail to describe consistently the zenith angle evolution of the muon content of EAS in the aforesaid energy regime.

Cosmic rays provide an unique approach to study hadronic interactions at high energies in the kinematic forward direction. The KASCADE air shower experiment was the first to conduct quantitative tests of hadronic interactions with air shower data. A brief overview is given on results from KASCADE and its extension KASCADE-Grande with respect to investigations of hadronic interactions and the properties of cosmic rays.

KASCADE and KASCADE-Grande were multi-detector installations to measure individual air showers of cosmic rays at ultra-high energy. Based on data sets measured by KASCADE and KASCADE-Grande, 90% C.L. upper limits to the flux of gamma-rays in the primary cosmic ray flux are determined in an energy range of ${10}^{14} - {10}^{18}$ eV. The analysis is performed by selecting air showers with a low muon content as expected for gamma-ray-induced showers compared to air showers induced by energetic nuclei. The best upper limit of the fraction of gamma-rays to the total cosmic ray flux is obtained at $3.7 \times {10}^{15}$ eV with $1.1 \times {10}^{-5}$. Translated to an absolute gamma-ray flux this sets constraints on some fundamental astrophysical models, such as the distance of sources for at least one of the IceCube neutrino excess models.

We aim to explain in a unified way the experimental data on ultrahigh energy cosmic rays (UHECR) and neutrinos, using a single source class and obeying limits on the extragalactic diffuse gamma-ray background (EGRB). If UHECRs only interact hadronically with gas around their sources, the resulting diffuse CR flux can be matched well to the observed one, providing at the same time large neutrino fluxes. However, air showers in the Earth's atmosphere induced by UHECRs with energies $E>3\times 10^{18}$ eV would reach in such a case their maxima too high. Therefore additional photo-hadronic interactions of UEHCRs close to the accelerator have to be present, in order to modify the nuclear composition of CRs in a relatively narrow energy interval. We include thus both photon and gas backgrounds, and combine the resulting CR spectra with the high-energy part of the Galactic CR fluxes predicted by the escape model. As result, we find a good description of experimental data on the total CR flux, the mean shower depth $X_\max$ and its width $RMS(X_\max)$ in the whole energy range above $E\simeq 10^{17}$ eV. The predicted high-energy neutrino flux matches IceCube measurements, while the contribution to the EGRB is of order 30%.

Local interstellar spectra (LIS) for protons, helium and antiprotons are built using the most recent experimental results combined with the state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod that provides modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameters optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997-2015. The proposed LIS accommodate both the low energy interstellar CR spectra measured by Voyager 1 and the high energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

The present status of high energy cosmic ray interaction models is discussed, concentrating on recent model updates inspired by the data from Run 1 of the LHC. A special attention is devoted to the remaining differences in the model predictions and their relation to the underlying theoretical approaches. Opportunities for the model discrimination by future LHC and cosmic ray experimental studies are analyzed.

We demonstrate that a substantial part of the present uncertainties in model predictions for the average maximum depth of cosmic ray-induced extensive air showers is related to very high energy pion-air collisions. Our analysis shows that the position of the maximum of the muon production profile in air showers is strongly sensitive to the properties of such interactions. Therefore, the measurements of the maximal muon production depth by cosmic ray experiments provide a unique opportunity to constrain the treatment of pion-air interactions at very high energies and to reduce thereby model-related uncertainties for the shower maximum depth.

KCDC, the KASCADE Cosmic-ray Data Centre, is a web portal, where data of astroparticle physics experiments will be made available for the interested public. The KASCADE experiment, financed by public money, was a large-area detector for the measurement of high-energy cosmic rays via the detection of air showers. KASCADE and its extension KASCADE-Grande stopped finally the active data acquisition of all its components including the radio EAS experiment LOPES end of 2012 after more than 20 years of data taking. In a first release, with KCDC we provide to the public the measured and reconstructed parameters of more than 160 million air showers. In addition, KCDC provides the conceptional design, how the data can be treated and processed so that they are also usable outside the community of experts in the research field. Detailed educational examples make a use also possible for high-school students and early stage researchers.

A dramatic increase in the accuracy and statistics of space-borne cosmic ray (CR) measurements has yielded several breakthroughs over the last several years. The most puzzling is the rise in the positron fraction above ~10 GeV over the predictions of the propagation models assuming pure secondary production. The accuracy of the antiproton production cross section is critical for astrophysical applications and searches for new physics since antiprotons in CRs seem to hold the keys to many puzzles including the origin of those excess positrons. However, model calculations of antiproton production in CR interactions with interstellar gas are often employing parameterizations that are out of date or are using outdated physical concepts. That may lead to an incorrect interpretation of antiproton data which could have broad consequences for other areas of astrophysics. In this work, we calculate antiproton production in pp-, pA-, and AA-interactions using EPOS-LHC and QGSJET-II-04, two of the most advanced Monte Carlo (MC) generators tuned to numerous accelerator data including those from the Large Hadron Collider (LHC). We show that the antiproton yields obtained with these MC generators differ by up to an order of magnitude from yields of parameterizations commonly used in astrophysics.

We calculate the neutrino yield from collisions of cosmic ray (CR) nuclei on gas using the event generator QGSJET-II. We present first the general characteristics and numerical results for the neutrino yield assuming power-law fluxes for the primary CR nuclei. Then we use three parameterisations for the Galactic CR flux to derive the neutrino yield for energies around and above the knee. The shape and the normalization of the resulting neutrino flux above $\sim 10^{14}$ eV depends on the composition of the Galactic CR flux employed, but is generally dominated by its proton component. The spectral shape and magnitude of the neutrino flux suggest that the IceCube excess is not connected to interactions of Galactic sea CRs. If a fraction of these events has a Galactic origin, then they may be caused by CR overdensities around recent close-by Galactic sources.

The concept of the nuclear enhancement factor has been used since the beginning of gamma-ray astronomy. It provides a simple and convenient way to account for the contribution of nuclei (A>1) in cosmic rays (CRs) and in the interstellar medium (ISM) to the diffuse gamma-ray emission. An accurate treatment of the dominant emission process, such as hadronic interactions of CRs with the ISM, enables one to study CR acceleration processes, CR propagation in the ISM, and provides a reliable background model for searches of new phenomena. The Fermi Large Area Telescope (Fermi-LAT) launched in 2008 provides excellent quality data in a wide energy range 30 MeV - 1 TeV where the diffuse emission accounts for the majority of photons. Exploiting its data to the fullest requires a new study of the processes of gamma-ray production in hadronic interactions. In this paper we point out that several commonly used studies of the nuclear enhancement factor miss to account for the spectrally averaged energy loss fraction which ensures that the energy fraction transferred to photons is averaged properly with the spectra of CR species. We present a new calculation of the spectrally averaged energy loss fraction and the nuclear enhancement factor using the QGSJET-II-04 and EPOS-LHC interaction models.

KASCADE-Grande is dedicated for investigations of cosmic-ray air showers in the primary energy range from 10 PeV to 1 EeV. The multi-detector system allows us to reconstruct charged particles, electron and muon numbers for individual air showers with high accuracies. Based on the shower size ($N_{ch}$) spectra of the charged particle component, the all-particle energy spectrum of cosmic rays is reconstructed, where attenuation effects in the atmosphere are corrected by applying the constant intensity cut method. The energy calibration is performed by using CORSIKA simulations with high-energy interaction models QGSJET-II-2, QGSJET-II-4, EPOS 1.99 and SIBYLL 2.1, where FLUKA has been used as low-energy interaction model for all cases. In the different hadronic models, different abundances for shower particles are predicted. Such model differences in the observables will be compared and discussed in this contribution. Furthermore, by using data with increasing statistics, the updated energy spectra by means of different interaction models will be presented.

The KASCADE-Grande detector is an air-shower array devoted to the study of primary cosmic rays with very high-energies (E = 10^{16} - 10^{18} eV). The instrument is composed of different particle detector systems suitable for the detailed study of the properties of Extensive Air Showers (EAS) developed by cosmic rays in the atmosphere. Among the EAS observables studied with the detector, the charged number of particles, the muon content (at different energy thresholds), and the number of electrons are found. By comparing the measurements of these air-shower parameters with the expectations from MC simulations, different hadronic interaction models can be tested at the high-energy regime with the KASCADE-Grande experiment. In this work, the results of a study on the evolution of the muon content of EAS with zenith angle, performed with the KASCADE-Grande instrument, is presented. Measurements are compared with predictions from MC simulations based on the QGSJET II, QGSJET II-04, SIBYLL 2.1 and EPOS 1.99 hadronic interaction models. A mismatch between experiment and simulations is observed. A similar problem is found for the evolution of the lateral distribution function of muons in the atmosphere.

The KASCADE-Grande experiment, located at KIT-Karlsruhe, Germany, consists of a large scintillator array for measurements of charged particles, N_ch, and of an array of shielded scintillation counters used for muon counting, N_mu. KASCADE-Grande is optimized for cosmic ray measurements in the energy range 10 PeV to 1000 PeV, thereby enabling the verification of a knee in the iron spectrum expected at approximately 100 PeV. Exploring the composition in this energy range is of fundamental importance for understanding the transition from galactic to extragalactic cosmic rays. Following earlier studies of elemental spectra reconstructed in the knee energy range from KASCADE data, we have now extended these measurements to beyond 100 PeV. By analysing the two-dimensional shower size spectrum N_ch vs. N_mu, we reconstruct the energy spectra of different mass groups by means of unfolding methods. The procedure and its results, giving evidence for a knee-like structure in the spectrum of iron nuclei, will be presented.

The KASCADE-Grande Muon Tracking Detector enables with high accuracy the measurement of directions of EAS muons with energy above 0.8 GeV and up to 700 m distance from the shower centre. Reconstructed muon tracks are used to investigate muon pseudorapidity (eta) distributions. These distributions are nearly identical to the pseudorapidity distributions of their parent mesons produced in hadronic interactions. Comparison of the eta distributions from measured and simulated showers can be used to test the quality of the high energy hadronic interaction models. In this context a comparison of the QGSJet-II-2 and QGSJet-II-4 model will be shown. The pseudorapidity distributions reflect the longitudinal development of EAS and, as such, are sensitive to the mass of the cosmic rays primary particles. With various parameters of the eta distribution, obtained from the MTD data, it is possible to calculate the mean logarithmic mass of CRs. The results of the <lnA> analysis in the primary energy range 10^{16} eV - 10^{17} eV with the 1st quartile (Q1) of eta distribution will be presented.

Data of the Grande extension of the KASCADE experiment allows us to study extensive air showers induced by primary cosmic rays with energies above 10^{16} eV. The energy of an event is estimated in terms of the number of charged particles (Nch ) and the number of muons (N{\mu} ) measured at an altitude of 110 m a.s.l. While a combination of the two numbers is used for the energy, the ratio defines the primary mass (group). The spectrum of the combined light and medium mass components, recently measured with KASCADE-Grande, was found to be compatible with both a single power-law and a broken power-law in the energy range between 10^{16.3} and 10^{18} eV. In this contribution we will present the investigation of possible structures in the spectrum of light primaries with increased statistics both from a larger data set including more recent measurements and by using a larger fiducial area than in the previous study. With the better statistical accuracy and with optimized selection criteria for enhancing light primaries we have found evidence for a hardening (ankle) of the spectrum of the light component at an energy of 10^{17.08+-0.08} eV.

The detection of high-energy cosmic rays above a few hundred TeV is realized by the observation of extensive air-showers. By using the multi-detector setup of KASCADE-Grande, energy spectrum, elemental composition, and anisotropies of high-energy cosmic rays in the energy range from below the knee up to 2 EeV are investigated. In addition, the large high-quality data set permits distinct tests of the validity of hadronic interaction models used in interpreting air-shower measurements. After more than 16 years, the KASCADE-Grande experiment terminated measurements end of 2012. This contribution will give an overview of the main results of the data analysis achieved so far, and will report about the status of KCDC, the KASCADE Cosmic-ray Data Center, where via a web-based interface the data will be made available for the interested public.

The KASCADE-Grande air shower experiment [W. Apel, et al. (KASCADE-Grande collaboration), Nucl. Instrum. Methods A 620 (2010) 202] consists of, among others, a large scintillator array for measurements of charged particles, Nch, and of an array of shielded scintillation counters used for muon counting, Nmu. KASCADE-Grande is optimized for cosmic ray measurements in the energy range 10 PeV to about 2000 PeV, where exploring the composition is of fundamental importance for understanding the transition from galactic to extragalactic origin of cosmic rays. Following earlier studies of the all-particle and the elemental spectra reconstructed in the knee energy range from KASCADE data [T. Antoni, et al. (KASCADE collaboration), Astropart. Phys. 24 (2005) 1], we have now extended these measurements to beyond 200 PeV. By analysing the two-dimensional shower size spectrum Nch vs. Nmu for nearly vertical events, we reconstruct the energy spectra of different mass groups by means of unfolding methods over an energy range where the detector is fully efficient. The procedure and its results, which are derived based on the hadronic interaction model QGSJET-II-02 and which yield a strong indication for a dominance of heavy mass groups in the covered energy range and for a knee-like structure in the iron spectrum at around 80 PeV, are presented. This confirms and further refines the results obtained by other analyses of KASCADE-Grande data, which already gave evidence for a knee-like structure in the heavy component of cosmic rays at about 80 PeV [W. Apel, et al. (KASCADE-Grande collaboration), Phys. Rev. Lett. 107 (2011) 171104].

The aim of this report of the Working Group on Hadronic Interactions and Air Shower Simulation is to give an overview of the status of the field, emphasizing open questions and a comparison of relevant results of the different experiments. It is shown that an approximate overall understanding of extensive air showers and the corresponding hadronic interactions has been reached. The simulations provide a qualitative description of the bulk of the air shower observables. Discrepancies are however found when the correlation between measurements of the longitudinal shower profile are compared to that of the lateral particle distributions at ground. The report concludes with a list of important problems that should be addressed to make progress in understanding hadronic interactions and, hence, improve the reliability of air shower simulations.

Recent results of the KASCADE-Grande experiment provided evidence for a mild knee-like structure in the all-particle spectrum of cosmic rays at $E = 10^{16.92 \pm 0.10} \, \mathrm{eV}$, which was found to be due to a steepening in the flux of heavy primary particles. The spectrum of the combined components of light and intermediate masses was found to be compatible with a single power law in the energy range from $10^{16.3} \, \mathrm{eV}$ to $10^{18} \, \mathrm{eV}$. In this paper, we present an update of this analysis by using data with increased statistics, originating both from a larger data set including more recent measurements and by using a larger fiducial area. In addition, optimized selection criteria for enhancing light primaries are applied. We find a spectral feature for light elements, namely a hardening at $E = 10^{17.08 \pm 0.08} \, \mathrm{eV}$ with a change of the power law index from $-3.25 \pm 0.05$ to $-2.79 \pm 0.08$.

We calculate the diffuse intensity of cosmic ray (CR) nuclei and their secondaries in the Boron-Carbon group produced by supernova remnants (SNR). The trajectories of charged particles in the SNR are modeled as a random walk in the test particle approximation. Secondary production by CRs colliding with gas in the SNR is included as a Monte Carlo process, while we use Galprop to account for the propagation and interactions of CRs in the Galaxy. In the vicinity of a source, we find an approximately constant B/C ratio as a function of energy. As a result, the B/C ratio at Earth does not rise with energy, but flattens instead in the high energy limit. This prediction can be soon tested by the AMS-2 collaboration.

A fundamental problem of cosmic ray (CR) physics is the determination of the average properties of Galactic CRs outside the Solar system. Starting from COS-B data in the 1980's, gamma-ray observations of molecular clouds in the Gould Belt above the Galactic plane have been used to deduce the Galactic CR energy spectrum. We reconsider this problem in view of the improved precision of observational data which in turn require a more precise treatment of photon production in proton-proton scatterings. We show that the spectral shape $dN/dp\propto p^{-2.85}$ of CR protons as determined by the PAMELA collaboration in the energy range 80 GeV<pc<230 GeV is consistent with the photon spectra from molecular clouds observed with Fermi-LAT down to photon energies E\sim 1-2 GeV. Adding a break of the CR flux at 3 GeV, caused by a corresponding change of the diffusion coefficient, improves further the agreement in the energy range 0.2-3 GeV.

The energy spectrum of cosmic rays between 10**16 eV and 10**18 eV, derived from measurements of the shower size (total number of charged particles) and the total muon number of extensive air showers by the KASCADE-Grande experiment, is described. The resulting all-particle energy spectrum exhibits strong hints for a hardening of the spectrum at approximately 2x10**16 eV and a significant steepening at c. 8x10**16 eV. These observations challenge the view that the spectrum is a single power law between knee and ankle. Possible scenarios generating such features are discussed in terms of astrophysical processes that may explain the transition region from galactic to extragalactic origin of cosmic rays.

We show that images of TeV blazars in the GeV energy band should contain, along with point-like sources, degree-scale jet-like extensions. These GeV extensions are the result of electromagnetic cascades initiated by TeV gamma-rays interacting with extragalactic background light and the deflection of the cascade electrons/positrons in extragalactic magnetic fields (EGMF). Using Monte-Carlo simulations, we study the spectral and timing properties of the degree-scale extensions in simulated GeV band images of TeV blazars. We show that the brightness profile of such degree-scale extensions can be used to infer the lightcurve of the primary TeV gamma-ray source over the past 1e7 yr, i.e. over a time scale comparable to the life-time of the parent active galactic nucleus. This implies that the degree-scale jet-like GeV emission could be detected not only near known active TeV blazars, but also from "TeV blazar remnants", whose central engines were switched off up to ten million years ago. Since the brightness profile of the GeV "jets" depends on the strength and the structure of the EGMF, their observation provides additionally information about the EGMF.

A Monte Carlo program for the simulation of electromagnetic cascades initiated by high-energy photons and electrons interacting with extragalactic background light (EBL) is presented. Pair production and inverse Compton scattering on EBL photons as well as synchrotron losses and deflections of the charged component in extragalactic magnetic fields (EGMF) are included in the simulation. Weighted sampling of the cascade development is applied to reduce the number of secondary particles and to speed up computations. As final result, the simulation procedure provides the energy, the observation angle, and the time delay of secondary cascade particles at the present epoch. Possible applications are the study of TeV blazars and the influence of the EGMF on their spectra or the calculation of the contribution from ultrahigh energy cosmic rays or dark matter to the diffuse extragalactic gamma-ray background. As an illustration, we present results for deflections and time-delays relevant for the derivation of limits on the EGMF.

Contributions of the KASCADE-Grande Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

We report the observation of a steepening in the cosmic ray energy spectrum of heavy primary particles at about $8 \times 10^{16}$ eV. This structure is also seen in the all-particle energy spectrum, but is less significant. Whereas the `knee' of the cosmic ray spectrum at $3-5 \times 10^{15}$ eV was assigned to light primary masses by the KASCADE experiment, the new structure found by the KASCADE-Grande experiment is caused by heavy primaries. The result is obtained by independent measurements of the charged particle and muon components of the secondary particles of extensive air showers in the primary energy range of $10^{16}$ to $10^{18}$ eV. The data are analyzed on a single-event basis taking into account also the correlation of the two observables.

The workshop on "Hadron-Hadron and Cosmic-Ray Interactions at multi-TeV Energies" held at the ECT* centre (Trento) in Nov.-Dec. 2010 gathered together both theorists and experimentalists to discuss issues of the physics of high-energy hadronic interactions of common interest for the particle, nuclear and cosmic-ray communities. QCD results from collider experiments -- mostly from the LHC but also from the Tevatron, RHIC and HERA -- were discussed and compared to various hadronic Monte Carlo generators, aiming at an improvement of our theoretical understanding of soft, semi-hard and hard parton dynamics. The latest cosmic-ray results from various ground-based observatories were also presented with an emphasis on the phenomenological modeling of the first hadronic interactions of the extended air-showers generated in the Earth atmosphere. These mini-proceedings consist of an introduction and short summaries of the talks presented at the meeting.

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.

Hadronic interaction models at cosmic ray (CR) energies are inherently uncertain due to the lack of a fundamental theoretical description of soft hadronic and nuclear interactions and the large extrapolation required from collider energies to the range of the most energetic cosmic rays observed ($>10^{20}$ eV). Model uncertainties are evaluated within the QGSJET-II model, by varying some of the crucial parameters in the limits allowed by collider data, and between QGSJET-II and other models commonly used in air shower simulations. The crucial parameters relate to hard processes, string fragmentation, diffraction and baryon production. Results on inelastic cross sections, on secondary particle production and on the properties of air showers measured by ground detectors from energies of $10^{12}$ to $10^{19}$ eV are presented.

We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a supernova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the trajectories of charged particles as a random walk with a prescribed diffusion coefficient, accelerating the particles at each shock crossing. Secondary production by CRs colliding with gas is included as a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The positron/electron ratio and the antiproton/proton ratio are a few percent and few $\times 10^{-5}$, respectively. Moreover, the obtained positron/electron ratio decreases with energy, while the antiproton/proton ratio rises at most by a factor of two above 10 GeV.

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.

Ultrahigh energy cosmic ray (UHECR) protons interacting with the cosmic microwave background (CMB) produce UHE electrons and gamma-rays that in turn initiate electromagnetic cascades on CMB and infrared photons. As a result, a background of diffuse isotropic gamma radiation is accumulated in the energy range $E\lsim 100$ GeV. The Fermi-LAT collaboration has recently reported a measurement of the extragalactic diffuse background finding it less intense and softer than previously measured by EGRET. We show that this new result constrains UHECR models and the flux of cosmogenic neutrinos. In particular, it excludes models with cosmogenic neutrino fluxes detectable by existing neutrino experiments, while next-generation detectors as e.g. JEM-EUSO can observe neutrinos only for extreme parameters.

High energy photons from blazars can initiate electromagnetic pair cascades interacting with the extragalactic photon background. The charged component of such cascades is deflected and delayed by extragalactic magnetic fields (EGMF), reducing thereby the observed point-like flux and leading potentially to multi degree images in the GeV energy range. We calculate the fluence of 1ES 0229+200 as seen by Fermi-LAT for different EGMF profiles using a Monte Carlo simulation for the cascade development. The non-observation of 1ES 0229+200 by Fermi-LAT suggests that the EGMF fills at least 60% of space with fields stronger than {\cal O}(10^{-16}-10^{-15})G for life times of TeV activity of {\cal O}(10^2-10^4)yr. Thus the (non-) observation of GeV extensions around TeV blazars probes the EGMF in voids and puts strong constraints on the origin of EGMFs: Either EGMFs were generated in a space filling manner (e.g. primordially) or EGMFs produced locally (e.g. by galaxies) have to be efficiently transported to fill a significant volume fraction, as e.g. by galactic outflows.

A large area (128 m^2) Muon Tracking Detector (MTD), located within the KASCADE experiment, has been built with the aim to identify muons (E_mu > 0.8 GeV) and their directions in extensive air showers by track measurements under more than 18 r.l. shielding. The orientation of the muon track with respect to the shower axis is expressed in terms of the radial- and tangential angles. By means of triangulation the muon production height H_mu is determined. By means of H_mu, a transition from light to heavy cosmic ray primary particle with increasing shower energy Eo from 1-10 PeV is observed. Muon pseudorapidity distributions for the first interactions above 15 km are studied and compared to Monte Carlo simulations.

The present status of extensive air shower (EAS) simulation procedures is reviewed. The advantages of combining numerical and Monte Carlo methods for the description of EAS development are discussed. Physics content of cosmic ray interaction models is briefly described and their predictions are compared to the first LHC data. Finally, some outstanding puzzles related to cosmic ray composition at the "ankle" energies are analyzed.

KASCADE-Grande is a large detector array dedicated for studies of high-energy cosmic rays in the primary energy range from 100 TeV to 1 EeV. The multi-detector concept of the experimental set-up offers the possibility to measure simultaneously various observables related to the electromagnetic, muonic, and hadronic air shower components. The experimental data are compared to predictions of CORSIKA simulations using high-energy hadronic interaction models (e.g. QGSJET or EPOS), as well as low-energy interaction models (e.g. FLUKA or GHEISHA). This contribution will summarize the results of such investigations. In particular, the validity of the new EPOS version 1.99 for EAS with energy around 100 PeV will be discussed.

The all-particle energy spectrum as measured by the KASCADE-Grande experiment for E = 10^{16} - 10^{18} eV is presented within the framework of the QGSJET II/FLUKA hadronic interaction models. Three different methods were applied based on the muon size and the total number of charged particles individually and in combination. From the study it is found that the spectrum cannot be completely described by a smooth power law due to the presence of characteristic features.

We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a supernova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the trajectories of charged particles as a random walk with a prescribed diffusioncoefficient, accelerating the particles at each shock crossing. Secondary production by CRs colliding with gas is included as a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The positron/electron ratio and the antiproton/proton ratio are a few percent and few $\times 10^{-5}$, respectively. Both ratios do not rise with energy.

TeV gamma-rays have been observed from blazars as well as from radio galaxies like M87 and Cen A. In leptonic models, gamma-rays above the pair production threshold can escape from the ultra-relativistic jet, since large Lorentz factors reduce the background photon densities compared to those required for isotropic emission. Here we discuss an alternative scenario, where VHE photons are generated as secondaries from UHECR interaction in the AGN core. We show that TeV gamma-rays can escape from the core despite large IR and UV backgrounds. For the special case of Cen A, we study if the various existing observations from the far infra-red to the UHE range can be reconciled within this picture.

High energy photons from blazars interact within tens of kpc with the extragalactic photon background, initiating electromagnetic pair cascades. The charged component of such cascades is deflected by extragalactic magnetic fields (EGMF), leading to halos even around initially point-like sources. We calculate the intensity profile of the resulting secondary high-energy photons for different assumptions on the initial source spectrum and the strength of the EGMF, employing also fields found earlier in a constrained simulation of structure formation including MHD processes. We find that the observation of halos around blazars like Mrk~180 probes an interesting range of EGMF strengths and acceleration models: In particular, blazar halos test if the photon energy spectrum at the source extends beyond \sim 100 TeV and how anisotropic this high energy component is emitted.

14 KASCADE-Grande reports submitted to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009

We calculated for the nearest active galactic nucleus (AGN), Centaurus A (Cen A), the flux of high energy cosmic rays and of accompanying secondary photons and neutrinos expected from hadronic interactions in the source. We used as two basic models for the generation of ultrahigh energy cosmic rays (UHECR) shock acceleration in the radio jet and acceleration in the regular electromagnetic field close to the core of the AGN, normalizing the UHECR flux to the observations of the Auger experiment. Here we compare the previously obtained photon fluxes with the recent data reported by the Fermi LAT and H.E.S.S. collaborations. In the case of the core model, we find good agreement both for the predicted spectral shape and the overall normalization between our earlier results and the H.E.S.S. observations for a primary proton energy $dN/dE\propto E^{-\alpha}$ with $\alpha\sim 2$ or smaller. A broken-power law with high-energy part $\alpha=-2.7$ leads to photon fluxes in excess of the Fermi measurements. The energy spectrum of the photon fluxes obtained by us for the jet scenario is in all cases at variance with the H.E.S.S. and Fermi observations.

We calculate for the nearest active galactic nucleus (AGN), Centaurus A, the flux of high energy cosmic rays and of accompanying secondary photons and neutrinos expected from hadronic interactions in the source. We use as two basic models for the generation of ultrahigh energy cosmic rays (UHECR) shock acceleration in the radio jet and acceleration in the regular electromagnetic field close to the core of the AGN. While scattering on photons dominates in scenarios with acceleration close to the core, scattering on gas becomes more important if acceleration takes place along the jet. Normalizing the UHECR flux from Centaurus A to the observations of the Auger experiment, the neutrino flux may be marginally observable in a 1 km$^3$ neutrino telescope, if a steep UHECR flux $\d N/\d E\propto E^{-\alpha}$ with $\alpha=2.7$ extends down to $10^{17}$ eV. The associated photon flux is close to or exceeds the observational data of atmospheric Cherenkov and $\gamma$-ray telescopes for $\alpha\gsim 2$. In particular, we find that already present data favour either a softer UHECR injection spectrum than $\alpha=2.7$ for Centaurus A or a lower UHECR flux than expected from the normalization to the Auger observations.

On basis of the theorem of a universal shower development stating that a hadronically generated extensive air shower is completely described by the primary energy, the position of the shower maximum and a parameter related to the total muon number, the so-called correlation curve method is developed and applied to KASCADE data. Correlation information of the muon and electron content of showers measured by the KASCADE experiment are used for the reconstruction of energy and mass of primary cosmic rays. Systematic uncertainties of the method and the results are discussed in detail. It is shown that by this method general tendencies in spectrum and composition indeed can be revealed, but the absolute normalization in energy and mass scale requires much more detailed simulations.

Interpretation of EAS measurements strongly depends on detailed air shower simulations. The uncertainty in the prediction of shower observables for different primary particles and energies is currently dominated by differences between hadronic interaction models. The new models QGSJETII-3 and EPOS 1.6, which reproduce all major results of existing accelerator data (including detailed data of RHIC experiments for EPOS), have been implemented in the air shower simulation programs CORSIKA and CONEX. We show predictions of these new models and compare them with those from older models such as QGSJET01 or SIBYLL. Results for important air shower observables are discussed in detail.

We discuss the signatures of the transition from galactic to extragalactic cosmic rays in different scenarios, giving most attention to the dip scenario. The dip is a feature in the diffuse spectrum of ultra-high energy (UHE) protons in the energy range $1\times 10^{18} - 4\times 10^{19}$ eV, which is caused by electron-positron pair production on the cosmic microwave background (CMB) radiation. The dip scenario provides a simple physical description of the transition from galactic to extragalactic cosmic rays. Here we summarize the signatures of the pair production dip model for the transition, most notably the spectrum, the anisotropy and the chemical composition. The main focus of our work is however on the description of the features that arise in the elongation rate and in the distribution of the depths of shower maximum $X_{\rm max}$ in the dip scenario. We find that the curve for $X_{\max}(E)$ shows a sharp increase with energy, which reflects a sharp transition from an iron dominated flux at low energies to a proton dominated flux at $E\sim 10^{18}$ eV. We also discuss in detail the shape of the $X_{\max}$ distributions for cosmic rays of given energy and demonstrate that this represents a powerful tool to discriminate between the dip scenario and other possible models of the transition.

We calculate the yield of high energy neutrinos produced in astrophysical sources for arbitrary interaction depths $\tau_0$ and magnetic field strengths $B$. We take into account energy loss processes like synchrotron radiation and diffusion of charged particles in turbulent magnetic fields as well as the scattering of secondaries on background photons and the direct production of charm neutrinos. Meson-photon interactions are simulated with an extended version of the SOPHIA model. Diffusion leads to an increased path-length before protons leave the source of size R_s and therefore magnetized sources lose their transparency below the energy $E\sim 10^{18}{\rm eV} (R_s/{\rm pc}) (B/{\rm mG}) \tau_0^{1/\alpha}$, with $\alpha=1/3$ and 1 for Kolmogorov and Bohm diffusion, respectively. Moreover, the neutrino flux is suppressed above the energy where synchrotron energy losses become important for charged particles. As a consequence, the energy spectrum and the flavor composition of neutrinos are strongly modified both at low and high energies even for sources with $\tau_0\lsim 1$.

The status of present theoretical description of very high energy hadronic interactions is reviewed. The impact of new results of accelerator and cosmic ray experiments on hadronic interaction model constructions is discussed in detail. Special attention is payed to remaining uncertainties in model extrapolations into the ultra-high energy domain, in particular, concerning model predictions for the muon component of extensive air showers. New promising theoretical approaches are outlined and future experimental prospects are discussed.

An efficient scheme for one-dimensional extensive air shower simulation and its implementation in the program CONEX are presented. Explicit Monte Carlo simulation of the high-energy part of hadronic and electromagnetic cascades in the atmosphere is combined with a numeric solution of cascade equations for smaller energy sub-showers to obtain accurate shower predictions. The developed scheme allows us to calculate not only observables related to the number of particles (shower size) but also ionization energy deposit profiles which are needed for the interpretation of data of experiments employing the fluorescence light technique. We discuss in detail the basic algorithms developed and illustrate the power of the method. It is shown that Monte Carlo, numerical, and hybrid air shower calculations give consistent results which agree very well with those obtained within the CORSIKA program.

The new hadronic Monte Carlo model QGSJET-II is applied for extensive air shower (EAS) calculations. The obtained results are compared to the predictions of the original QGSJET and of the SIBYLL 2.1 interaction models. It is shown that non-linear effects change substantially model predictions for hadron-nucleus interactions and produce observable effects for calculated EAS characteristics. Finally the impact of the new model on the interpretation of air shower array data is discussed.

A hybrid simulation code is developed that is suited for fast one-dimensional simulations of shower profiles, including fluctuations. It combines the Monte Carlo simulation of high energy interactions with a fast numerical solution of cascade equations for the resulting distributions of secondary particles. Results obtained with this new code, called CONEX, are presented and compared to CORSIKA predictions.

We present an introductory lecture on general features of cosmic rays, for non-experts, and some recent developments concerning cascade equations for air shower developments.

We study decays of superheavy particles $X$ into leptons. We show that they initiate cascades similar to QCD parton jets, if $m_X\gsim 10^6$ GeV. Electroweak cascading is studied and the energy spectra of the produced leptons are calculated in the framework of a broken SU(2) model of weak interactions. As application, important for the Z-burst model for ultrahigh energy cosmic rays, we consider decays of superheavy particles coupled on tree-level only to neutrinos and derive stringent limit for these decays from the observed diffuse extragalactic $\gamma$-ray flux.

We study the proposal that the cosmic ray primaries above the Greisen-Zatsepin-Kuzmin (GZK) cutoff are gluino-containing hadrons ($\tilde g$-hadrons). We describe the interaction of $\tilde g$-hadrons with nucleons in the framework of the Gribov-Regge approach using a modified version of the hadronic interaction model QGSJET for the generations of Extensive Air Showers (EAS). There are two mass windows marginally allowed for gluinos: $m_{\tilde g}\lsim 3$ GeV and $25\lsim m_{\tilde g}\lsim 35$ GeV. Gluino-containing hadrons corresponding to the second window produce EAS very different from the observed ones. Light $\tilde g$-hadrons corresponding to the first gluino window produce EAS similar to those initiated by protons, and only future detectors can marginally distinguish them. We propose a beam-dump accelerator experiment to search for $\tilde g$-hadrons in this mass window. We emphasize the importance of this experiment: it can discover (or exclude) the light gluino and its role as a cosmic ray primary at ultra high energies.

Practically all serious calculations of exclusive particle production in ultra-relativistic nuclear or hadronic interactions are performed in the framework of Gribov-Regge theory or the eikonalized parton model scheme. It is the purpose of this paper to point out serious inconsistencies in the above-mentioned approaches. We will demonstrate that requiring theoretical self-consistency reduces the freedom in modeling high energy nuclear scattering enormously. We will introduce a fully self-consistent formulation of the multiple-scattering scheme in the framework of a Gribov-Regge type effective theory. In addition, we develop new computational techniques which allow for the first time a satisfactory solution of the problem in the sense that calculation s of observable quantities can be done strictly within a self-consistent formalism.

The purpose of this paper is twofold: first we want to introduce a new type of hadronic interaction model (NEXUS), which has a much more solid theoretical basis as, for example, presently used models like QGSJET and VENUS, and ensures therefore a much more reliable extrapolation towards high energies. Secondly, we want to promote an extensive air shower (EAS) calculation scheme, based on cascade equations rather than explicit Monte Carlo simulations, which is very accurate in calculations of main EAS characteristics and extremely fast concerning computing time. We employ the NEXUS model to provide the necessary data on particle production in hadron-air collisions and present the average EAS characteristics for energies 10^14 - 10^17 eV. The experimental data of the casa-blanka group are analyzed in the framework of the new model.