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.

A surprisingly large flux of extraterrestrial high-energy neutrinos was discovered by the IceCube experiment. While the flux of muon neutrinos with energies $E>100$ TeV is consistent with the extragalactic gamma-ray background (EBL) determined by Fermi-LAT, the softer component of the cascade neutrino flux at $E<100$ TeV is larger than expected. Moreover, a gamma-ray excess at high Galactic latitudes at energies $E>300$ GeV was found in the data of Fermi-LAT. The gamma-ray excess at TeV energies and the neutrino excess at $E<100$ TeV may have a common Galactic origin. In this work, we study the possibility that both excesses are caused by interactions of cosmic rays (CRs) with energies up to PeV in the wall of the Local Bubble. Source of these CRs may be a recent nearby source like Vela. We show that such a scenario can explain the observed CR flux around the knee, while CR interactions in the bubble wall can generate a substantial fraction of the observed astrophysical high-energy neutrino flux below $\sim {\rm few} \times 100$ TeV.

Charged particles scattering on moving inhomogenities of the magnetised interstellar medium as can gain energy through the process of second-order Fermi acceleration. This energy gain depletes in turn the magnetic wave spectrum around the resonance wave-vector $k\sim 1/R_L$, where $R_L$ is the Larmor radius of the charged particle. This energy transfer can prohibit the cascading of magnetic turbulence to smaller scales, leading to a drop in the diffusion coefficient and allowing the efficient exchange of charged dark matter particles in the disk and the halo. As a result, terrestial limits from direct detection experiments apply to charged dark matter. Together with the no-observation of a drop in the diffusion coefficient, this excludes charged dark matter for $10^2 GeV\lesssim m/q \lesssim 10^{11} GeV$, even if the charged dark matter abundance is only a small part of the total relic abundance.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

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.

We review progress in high-energy cosmic ray physics focusing on recent experimental results and models developed for their interpretation. Emphasis is put on the propagation of charged cosmic rays, covering the whole range from $\sim (20-50)$ GV, i.e. the rigidity when solar modulations can be neglected, up to the highest energies observed. We discuss models aiming to explain the anomalies in Galactic cosmic rays, the knee, and the transition from Galactic to extragalactic cosmic rays.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

We model the contribution of the nearest young supernova remannt Vela to the local cosmic ray flux taking into account both the influence of the Local Superbubble and the effect of anisotropic diffusion. The dominant contribution of this source in the energy region around the cosmic ray knee can naturally explain the observed fluxes of individual groups of nuclei and their total flux. Adding the CR flux from a 2-3 Myr old local CR source suggested earlier, the CR spectra in the whole energy range between 200 GeV and the transition to extragalactic CRs are described well by the combined fluxes from these two local Galactic sources.

We examine the hypothesis of decaying heavy dark matter (HDM) in the context of the IceCube highest energy neutrino events and recent limits on the diffuse flux of high-energy photons. We consider dark matter (DM) particles $X$ of mass $10^{6}\leq M_X\leq~10^{16}$ GeV decaying on tree level into $X \rightarrow \nu \bar{\nu}$, $X \rightarrow e^+e^-$ and $X \rightarrow q \bar{q}$. The full simulation of hadronic and electroweak decay cascades and the subsequent propagation of the decay products through the interstellar medium allows us to determine the permitted values of $M_X$. We show that for leptonic decay channels it is possible to explain the IceCube highest energy neutrino signal without overproducing high-energy photons for $M_X~\lesssim~5.5 \cdot 10^{7}$ GeV and $1.5 \cdot 10^{8}~\lesssim~M_X~\lesssim~1.5 \cdot 10^{9}$ GeV, while hadronic decays contradict the gamma-ray limits for almost the whole range of $M_X$ values considered. The leptonic hypothesis can be probed by operating and planned gamma-ray observatories. For instance, the currently upgrading Carpet experiment will be capable to test a significant part of the remaining parameter window within one year of observation.

We argue that the diffusion of cosmic rays in the Galactic magnetic field has to be strongly anisotropic. As a result, the number of CR sources contributing to the local CR flux is reduced by a factor $\sim 200$. The CR density is therefore less smooth, and the contribution of individual sources to the CR dipole anisotropy becomes more prominent. In the case of anisotropic diffusion, the observed plateau in the CR dipole anisotropy around 2-20 TeV can be explained by a 2-3 Myr old CR source which dominates the local CR flux in this energy range.

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)

We calculate the diffusion coefficients of charged cosmic rays (CR) propagating in regular and turbulent magnetic fields. If the magnetic field is dominated by an isotropic turbulent component, we find that CRs reside too long in the Galactic disc. As a result, CRs overproduce secondary nuclei like boron for any reasonable values of the strength and the coherence length of an isotropic turbulent field. We conclude therefore that the propagation of Galactic CRs has to be strongly anisotropic because of a sufficiently strong regular field and/or of an anisotropy in the turbulent field. As a consequence, the number of sources contributing to the local CR flux is reduced by a factor ${\cal O}(100)$ compared to the case of isotropic CR diffusion.

A signal of high-energy extraterrestrial neutrinos from unknown source(s) was recently discovered by the IceCube experiment. Neutrinos are always produced together with gamma-rays, but the gamma-ray flux from extragalactic sources is suppressed due to attenuation in the intergalactic medium. We report the discovery of a gamma-ray excess at high Galactic latitudes starting at energies 300 GeV in the data of the Fermi telescope. We show that the multi-TeV gamma-ray diffuse emission has spectral characteristics at both low and high Galactic latitudes compatible with those of the IceCube high neutrino signal in the same sky regions. This suggests that these gamma-rays are the counterpart of the IceCube neutrino signal, implying that a sizable part of the IceCube neutrino flux originates from the Milky Way. We argue that the diffuse neutrino and gamma-ray signal at high Galactic latitudes originates either from previously unknown nearby cosmic ray "PeVatron" source(s), an extended Galactic CR halo or from decays of heavy dark matter particles.

We study the propagation of cosmic rays generated by sources residing inside superbubbles. We show that the enhanced magnetic field in the bubble wall leads to an increase of the interior cosmic ray density. Because of the large matter density in the wall, the probability for cosmic ray interactions on gas peaks there. As a result, the walls of superbubbles located near young cosmic ray sources emit efficiently neutrinos. We apply this scenario to the Loop~I and Local Superbubble: These bubbles are sufficiently near such that cosmic rays from a young source as Vela interacting in the bubble wall can generate a substantial fraction of the observed astrophysical high-energy neutrino flux below $\sim$ few $\times 100$ TeV.

The supernova explosion which deposited $^{60}$Fe isotopes on Earth 2-3 million years ago should have also produced cosmic rays which contribute to the locally observed cosmic ray flux. We show that the contribution of this "local source" explains the "anomalies" observed in the positron and antiproton fluxes and explains why their spectral shape agrees with that of the proton flux. At the same time, this local source component accounts for the difference in the slopes of the spectra of cosmic ray nuclei as the result of the slightly varying relative importance of the "local" and the average component for distinct CR nuclei. Such a "local supernova" model for the spectra of nuclei can be tested via a combined measurement of the energy dependence of the boron-to-carbon (primary-to-secondary cosmic rays) ratio and of the antiproton spectrum: While the antiproton spectrum is predicted to extend approximately as a power law into the TeV range without any softening break, the B/C ratio is expected to show a "plateau" at a level fixed by the observed positron excess in the 30-300 GeV range. We discuss the observability of such a plateau with dedicated experiments for the measurement of the cosmic ray composition in 10 TeV energy range (NUCLEON, ISS-CREAM).

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.

Recent 60Fe results have suggested that the estimated distances of supernovae in the last few million years should be reduced from 100 pc to 50 pc. Two events or series of events are suggested, one about 2.7 million years to 1.7 million years ago, and another may at 6.5 to 8.7 million years ago. We ask what effects such supernovae are expected to have on the terrestrial atmosphere and biota. Assuming that the Local Bubble was formed before the event being considered, and that the supernova and the Earth were both inside a weak, disordered magnetic field at that time, TeV-PeV cosmic rays at Earth will increase by a factor of a few hundred. Tropospheric ionization will increase proportionately, and the overall muon radiation load on terrestrial organisms will increase by a factor of 150. All return to pre-burst levels within 10kyr. In the case of an ordered magnetic field, effects depend strongly on the field orientation. The upper bound in this case is with a largely coherent field aligned along the line of sight to the supernova, in which case TeV-PeV cosmic ray flux increases are 10^4; in the case of a transverse field they are below current levels. We suggest a substantial increase in the extended effects of supernovae on Earth and in the lethal distance estimate; more work is needed.This paper is an explicit followup to Thomas et al. (2016). We also here provide more detail on the computational procedures used in both works.

We investigate how the extragalactic proton component derived within the "escape model" can be explained by astrophysical sources. We consider as possible cosmic ray (CR) sources normal/starburst galaxies and radio-loud active galactic nuclei (AGN). We find that the contribution to the total extragalactic proton flux from normal and starburst galaxies is only subdominant and does not fit the spectral shape deduced in the escape model. In the case of radio-loud AGN, we show that the complete extragalactic proton spectrum can be explained by a single source population, BL Lac/FR I, for any of the potential acceleration sites in these sources. We calculate the diffuse neutrino and $\gamma$-ray fluxes produced by these CR protons interacting with gas inside their sources. For a spectral slope of CRs close to $\alpha=2.1-2.2$ as suggested by shock acceleration, we find that these UHECR sources contribute the dominant fraction of both the isotropic $\gamma$-ray background and of the extragalactic part of the astrophysical neutrino signal observed by IceCube.

The locally observed cosmic ray spectrum has several puzzling features, such as the excess of positrons and antiprotons above $\sim 20$ GeV and the discrepancy in the slopes of the spectra of cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features are consistently explained by a nearby source which was active $\sim 2$ Myr ago and has injected $(1-2)\times 10^{50}$ erg in cosmic rays. The transient nature of the source and its overall energy budget point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that the local cosmic ray injection was produced by the same supernova that has deposited $^{60}$Fe isotopes in the deep ocean crust.

We study numerically the anisotropy of the cosmic ray (CR) flux emitted by a single source calculating the trajectories of individual CRs. We show that the contribution of a single source to the observed anisotropy is instead determined solely by the fraction the source contributes to the total CR intensity, its age and its distance,and does not depend on the CR energy at late times. Therefore the observation of a constant dipole anisotropy indicates that a single source dominates the CR flux in the corresponding energy range. A natural explanation for the plateau between 2--20 TeV observed in the CR anisotropy is thus the presence of a single, nearby source. For the source age of 2 Myr, as suggested by the explanation of the antiproton and positron data from PAMELA and AMS-02 through a local source [arXiv:astro-ph/1504.06472], we determine the source distance as $\sim 200$ pc. Combined with the contribution of the global CR sea calculated in the escape model, we can explain qualitatively the data for the dipole anisotropy. Our results suggest that the assumption of a smooth CR source distribution should be abandoned between 200 GeV and 1 PeV.

We show that the cosmic ray (CR) knee can be entirely explained by energy-dependent CR leakage from the Milky Way, with an excellent fit to all existing data. We test this hypothesis calculating the trajectories of individual CRs in the Galactic magnetic field. We find that the CR escape time $\tau_{\rm esc}(E)$ exhibits a knee-like structure around $E/Z={\rm few}\times 10^{15}$ eV for small coherence lengths and strengths of the turbulent magnetic field. The resulting intensities for different groups of nuclei are consistent with the ones determined by KASCADE and KASCADE-Grande, using simple power-laws as injection spectra. The transition from Galactic to extragalactic CRs is terminated at $\approx 2\times 10^{18}$ eV, while extragalactic CRs contribute sizeable to the subdominant proton flux already for $\gtrsim 2\times 10^{16}$ eV. The natural source of extragalactic CRs in the intermediate energy region up to the ankle are in this model normal and starburst galaxies. The escape model provides a good fit to $\ln(A)$ data; it predicts that the phase of the CR dipole varies strongly in the energy range between $1\times 10^{17}$ and $3\times 10^{18}$ eV, while our estimate for the dipole magnitude is consistent with observations.

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.

We investigate the possibility that the cosmic ray (CR) knee is entirely explained by the energy-dependent CR leakage from the Milky Way. We test this hypothesis calculating the trajectories of individual CRs with energies between $E/Z=10^{14}$ eV and $10^{17}$ eV propagating them in the regular and turbulent Galactic magnetic field. We find a knee-like structure of the CR escape time $\tau_{\rm esc}(E)$ around $E/Z={\rm few}\times 10^{15}$ eV for a coherence length $l_{\rm c} \simeq 2$ pc of the turbulent field, while the decrease of $\tau_{\rm esc}(E)$ slows down around $E/Z\simeq 10^{16}$ eV in models with a weak turbulent magnetic field. Assuming that the injection spectra of CR nuclei are power-laws, the resulting CR intensities in such a turbulence are consistent with the energy spectra of CR nuclei determined by KASCADE and KASCADE-Grande. We calculate the resulting CR dipole anisotropy as well as the source rate in this model.

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.

Analyses of TeV-PeV cosmic ray (CR) diffusion around their sources usually assume either isotropic diffusion or anisotropic diffusion due to the regular Galactic magnetic field. We show that none of them are adequate on distances smaller than the maximal scale Lmax ~ 100 pc of fluctuations in the turbulent interstellar magnetic field. As a result, we predict anisotropic gamma-ray emissions around CR proton and electron sources, even for uniform densities of target gas. The centers of extended emission regions may have non-negligible offsets from their sources, leading to risks of misidentification. Gamma-rays from CR filaments have steeper energy spectra than those from surrounding regions. We point out that gamma-ray telescopes can be used in the future as a new way to probe and deduce the parameters of the interstellar magnetic field.

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.

We constrain the energy at which the transition from Galactic to extragalactic cosmic rays occurs by computing the anisotropy at Earth of cosmic rays emitted by Galactic sources. Since the diffusion approximation starts to loose its validity for E/Z >~ 10^(16-17) eV, we propagate individual cosmic rays (CRs) using Galactic magnetic field models and taking into account both their regular and turbulent components. The turbulent field is generated on a nested grid which allows spatial resolution down to fractions of a parsec. Assuming sufficiently frequent Galactic CR sources, the dipole amplitude computed for a mostly light or intermediate primary composition exceeds the dipole bounds measured by the Auger collaboration around E ~ 10^18 eV. Therefore, a transition at the ankle or above would require a heavy composition or a rather extreme Galactic magnetic field with strength >~ 10 muG. Moreover, the fast rising proton contribution suggested by KASCADE-Grande data between 10^17 eV and 10^18 eV should be of extragalactic origin. In case heavy nuclei dominate the flux at E >~ 10^18 eV, the transition energy can be close to the ankle, if Galactic CRs are produced by sufficiently frequent transients as e.g. magnetars.

We investigate the diffusion of cosmic rays (CR) close to their sources. Propagating individual CRs in purely isotropic turbulent magnetic fields with maximal scale of spatial variations Lmax, we find that CRs diffuse anisotropically at distances r <~ Lmax from their sources. As a result, the CR densities around the sources are strongly irregular and show filamentary structures. We determine the transition time t* to standard diffusion as t* ~ 10^4 yr (Lmax/150 pc)^b (E/PeV)^(-g) (Brms/4 muG)^g, with b ~ 2 and g = 0.25-0.5 for a turbulent field with Kolmogorov power spectrum. We calculate the photon emission due to CR interactions with gas and the resulting irregular source images.

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.

In this work we study how the turbulent component of the Galactic magnetic field (GMF) affects the propagation of ultrahigh energy heavy nuclei. We investigate first how the images of individual sources and of the supergalactic plane depend on the properties of the turbulent GMF. Then we present a quantitative study of the impact of the turbulent field on (de-) magnification of source fluxes, due to magnetic lensing effects. We also show that it is impossible to explain the Pierre Auger data assuming that all ultrahigh energy nuclei are coming from Cen A, even in the most favorable case of a strong, extended turbulent field in the Galactic halo.

In this work, we present detailed simulations for propagation of ultra-high energy (UHE) heavy nuclei, with E > 60 EeV, within recent Galactic Magnetic Field (GMF) models. We investigate the impacts of the regular and turbulent components of the GMF. We show that with UHE heavy nuclei, there is no one-to-one correspondence between the arrival directions of cosmic rays (CR) measured at Earth and the direction of their extragalactic sources. Sources can have several distorted images on the sky. We compute images of galaxy clusters and of the supergalactic plane in recent GMF models and show the challenges, and possibilities, of "UHECR astronomy" with heavy nuclei. Finally, we present a quantitative study of the impact of the GMF on the (de-)magnification of source fluxes, due to magnetic lensing effects. We find that for 60 EeV iron nuclei, sources located in up to about one fifth of the sky would have their fluxes so strongly demagnified that they would not be detectable at Earth, even by the next generation of UHECR experiments.

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.

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.

Observations are consistent with a significant fraction of heavy nuclei in the cosmic ray flux above a few times 10^19 eV. Such nuclei can be deflected considerably in the Galactic magnetic field, with important implications for the search of their sources. We perform detailed simulations of heavy nuclei propagation within recent Galactic magnetic field models. While such models are not yet sufficiently constrained to predict deflection maps in detail, we find general features of the distribution of (de-) magnified flux from sources. Since in most theoretical models sources of heavy nuclei are located in the local large scale structure of galaxies, we show examples of images of several nearby galaxy clusters and of the supergalactic plane. Such general features may be useful to develop efficient methods for source reconstruction from observed ultrahigh energy cosmic ray arrival directions.

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.

Superheavy dark matter (SHDM) exchanges energy with its environment much slower than particles with masses close to the electroweak (EW) scale and has therefore different small-scale clustering properties. Using the neutralino as candidate for the SHDM, we find that free-streaming allows the formation of DM clumps of all masses down to $\sim 260 m_\chi$ in the case of bino. If small-scale clumps evolve from a non-standard, spiky spectrum of perturbations, DM clumps may form during the radiation dominated era. These clumps are not destroyed by tidal interactions and can be extremely dense. In the case of a bino, a "gravithermal catastrophe" can develop in the central part of the most dense clumps, increasing further the central density and thus the annihilation signal. In the case of a higgsino, the annihilation signal is enhanced by the Sommerfeld effect. As a result annihilations of superheavy neutralinos in dense clumps may lead to observable fluxes of annihilation products in the form of ultrahigh energy particles, for both cases, higgsinos and binos, as lightest supersymmetric particles.

The formation and evolution of superdense clumps (or subhalos) is studied. Such clumps of dark matter (DM) can be produced by many mechanisms, most notably by spiky features in the spectrum of inflationary perturbations and by cosmological phase transitions. Being produced very early during the radiation dominated epoch, superdense clumps evolve as isolated objects. They do not belong to hierarchical structures for a long time after production, and therefore they are not destroyed by tidal interactions during the formation of larger structures. For DM particles with masses close to the electroweak (EW) mass scale, superdense clumps evolve towards a power-law density profile $\rho(r) \propto r^{-1.8}$ with a central core. Superdense clumps cannot be composed of standard neutralinos, since their annihilations would overproduce the diffuse gamma radiation. If the clumps are constituted of superheavy DM particles and develop a sufficiently large central density, the evolution of their central part can lead to a 'gravithermal catastrophe.' In such a case, the initial density profile turns into an isothermal profile with $\rho \propto r^{-2}$ and a new, much smaller core in the center. Superdense clumps can be bserved by gamma radiation from DM annihilations and by gravitational wave detectors, while the production of primordial black holes and cascade nucleosynthesis constrain this scenario.

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.

Interpretations of indirect searches for dark matter (DM) require theoretical predictions for the annihilation or decay rates of DM into stable particles of the standard model. These predictions include usually only final states accessible as lowest order tree-level processes, with electromagnetic bremsstrahlung and the loop-suppressed two gamma-ray line as exceptions. We show that this restriction may lead to severely biased results for DM tailored to produce only leptons in final states and with mass in the TeV range. For such models, unavoidable electroweak bremsstrahlung of Z and W-bosons has a significant influence both on the branching ratio and the spectral shape of the final state particles. We work out the consequences for two situations: Firstly, the idealized case where DM annihilates at tree level with 100% branching ratio into neutrinos. For a given cross section, this leads eventually to "minimal yields" of photons, electrons, positrons and antiprotons. Secondly, the case where the only allowed two-body final states are electrons. The latter case is typical of models aimed at fitting cosmic ray e^- and e^+ data. We find that the multimessenger signatures of such models can be significantly modified with respect to results presented in the literature.

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.

We perform an autocorrelation study of the Auger data with the aim to constrain the number density n_s of ultrahigh energy cosmic ray (UHECR) sources, estimating at the same time the effect on $n_s$ of the systematic energy scale uncertainty and of the distribution of UHECR. The use of global analysis has the advantage that no biases are introduced, either in n_s or in the related error bar, by the a priori choice of a single angular scale. The case of continuous, uniformly distributed sources is nominally disfavored at 99% C.L. and the fit improves if the sources follow the large-scale structure of matter in the universe. The best fit values obtained for the number density of proton sources are within a factor ~2 around n_s ~ 10^{-4}/Mpc^3 and depend mainly on the Auger energy calibration scale, with lower densities being preferred if the current scale is correct. The data show no significant small-scale clustering on scales smaller than a few degrees. This might be interpreted as a signature of magnetic smearing of comparable size, comparable with the indication of a ~3 degrees magnetic deflection coming from cross-correlation results. The effects on the above results of some approximations done is also discussed.

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 discuss how lensing by magnetic fields in galaxy clusters affects ultrahigh energy cosmic ray (UHECR) observations. As specific example, we use Virgo together with the cluster magnetic fields obtained earlier in a constrained simulation of structure formation including MHD processes. We find that, if M87 is the single source of UHECRs from Virgo, the emitted flux is strongly anisotropic in the most interesting energy range, (50-100)EeV, and differs from the average value by a factor five or more for a significant fraction of observers. Since magnetic lensing is energy dependent, the external energy spectrum as seen by different observers varies strongly too. These anisotropies are averaged out in the case that all active galactic nuclei in Virgo emit UHECRs. In both cases, the anisotropies of the emitted UHECR flux may introduce an important bias in the interpretation of UHECR data like, e.g., the determination of the source density n_s and the source energy spectrum of UHECRs.

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.

The motivation and the current status of top-down models as sources of ultrahigh energy cosmic rays (UHECR) are reviewed. Stimulated by the AGASA excess, they were proposed as the main source of UHECRs beyond the GZK cutoff. Meanwhile searches for their signatures have limited their contribution to the UHECR flux to be subdominant, while the theoretical motivation for these searches remained strong: Topological defects are a generic consequence of Grand Unified Theories and superheavy particles are a creditable dark matter candidate. While Fermi/GLAST results should help to improve soon bounds on topological defects from the diffuse gamma-ray background, the most promising detection method are UHE neutrino searches. Superheavy dark matter can be restricted or detected by its characteristic galactic anisotropy combined with searches for UHE photons.

We discuss the GZK horizon of protons and present a method to constrain the injection spectrum of ultrahigh energy cosmic rays (UHECRs) from supposedly identified extragalactic sources. This method can be applied even when only one or two events per source are observed and is based on the analysis of the probability for a given source to populate different energy bins, depending on the actual CR injection spectral index. In particular, we show that for a typical source density of $4\times 10^{-5} Mpc^{-3}$, a data set of 100 events above $6\times 10^{19}$ eV allows one in 97% of all cases to distinguish a source spectrum $dN/dE\propto E^{-1.1}$ from one with $E^{-2.7}$ at 95% confidence level.

The arrival directions of ultrahigh energy cosmic rays (UHECRs) may show anisotropies on all scales, from just above the experimental angular resolution up to medium scales and dipole anisotropies. We find that a global comparison of the two-point auto-correlation function of the data with the one of catalogues of potential sources is a powerful diagnostic tool. In particular, this method is far less sensitive to unknown deflections in magnetic fields than cross-correlation studies while keeping a strong discrimination power among source candidates. We illustrate these advantages by considering ordinary galaxies, gamma ray bursts and active galactic nuclei as possible sources. Already the sparse publicly available data suggest that the sources of UHECRs may be a strongly clustered sub-sample of galaxies or of active galactic nuclei. We present forecasts for various cases of source distributions which can be checked soon by the Pierre Auger Observatory.

I give a concise introduction into high energy cosmic ray physics, including also few related aspects of high energy gamma-ray and neutrino astrophysics. The main emphasis is placed on astrophysical questions, and the level of the presentation is kept basic.

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$.

Deviations from isotropy have been a key tool to identify the sources and the primary type of cosmic rays (CRs) at low energies. We argue that anisotropies due to blind regions induced by the Galactic magnetic field, the cosmological Compton-Getting effect, medium-scale anisotropies reflecting the large-scale distribution of CR sources and the small-scale clustering of the CR arrival directions at the highest energies may play the same role for extragalactic CRs.

We calculate the yield of high energy neutrinos produced in astrophysical sources with negligible magnetic fields varying their interaction depth from nearly transparent to opaque. We take into account the scattering of secondaries on background photons as well as the direct production of neutrinos in decays of charm mesons. If multiple scattering of nucleons becomes important, the neutrino spectra from meson and muon decays are strongly modified with respect to transparent sources. Characteristic for neutrino sources containing photons as scattering targets is a strong energy-dependence of the ratio $R^0$ of $\nu_\mu$ and $\nu_e$ fluxes at the sources, ranging from $R^0=\phi_\mu/\phi_e\sim 0$ below threshold to $R^0\sim 4$ close to the energy where the decay length of charged pions and kaons equals their interaction length on target photons. Above this energy, the neutrino flux is strongly suppressed and depends mainly on charm production.

Deviations from isotropy have been a key tool to identify the origin and the primary type of cosmic rays at low energies. We suggest that the Compton-Getting effect can play a similar role at ultra-high energies: If at these energies the cosmic ray flux is dominated by sources at cosmological distances, then the movement of the Sun relative to the cosmic microwave background frame induces a dipole anisotropy at the 0.6% level. The energy dependence and the orientation of this anisotropy provide important information about the transition between galactic and extragalactic cosmic rays, the charge of the cosmic ray primaries, the galactic magnetic field and, at the highest energies, the energy-loss horizon of cosmic rays. A 3-sigma detection of this effect requires around 10^6 events in the considered energy range and is thus challenging but not impossible with present detectors. As a corollary we note that the Compton-Getting effect allows one also to constrain the fraction of the diffuse gamma-ray background emitted by sources at cosmological distance, with promising detection possibilities for the GLAST satellite.

We study the influence of the regular component of the Galactic magnetic field (GMF) on the arrival directions of ultra-high energy cosmic rays (UHECRs). We find that, if the angular resolution of current experiments has to be fully exploited, deflections in the GMF cannot be neglected even for E=10^20 eV protons, especially for trajectories along the Galactic plane or crossing the Galactic center region. On the other hand, the GMF could be used as a spectrograph to discriminate among different source models and/or primaries of UHECRs, if its structure would be known with sufficient precision. We compare several GMF models introduced in the literature and discuss for the example of the AGASA data set how the significance of small-scale clustering or correlations with given astrophysical sources are affected by the GMF. We point out that the non-uniform exposure to the extragalactic sky induced by the GMF should be taken into account estimating the significance of potential (auto-)correlation signals.

Superheavy particles are a natural candidate for the dark matter in the universe and our galaxy, because they are produced generically during inflation in cosmologically interesting amounts. The most attractive model for the origin of superheavy dark matter (SHDM) is gravitational production at the end of inflation. The observed cosmological density of dark matter determines the mass of the SHDM particle as $m_X=$(a few)$\:\times 10^{13}$ GeV, promoting it to a natural candidate for the source of the observed ultra-high energy cosmic rays (UHECR). After a review of the theoretical aspects of SHDM, we up-date its predictions for UHECR observations: no GZK cutoff, flat energy spectrum with $dN/dE\approx 1/E^{1.9}$, photon dominance and galactic anisotropy. We analyze the existing data and conclude that SDHM as explanation for the observed UHECRs is at present disfavored but not yet excluded. We calculate the anisotropy relevant for future Auger observations that should be the conclusive test for this model. Finally, we emphasize that negative results of searches for SHDM in UHECR do not disfavor SHDM as a dark matter candidate. Therefore, UHECRs produced by SHDM decays and with the signatures as described should be searched for in the future as subdominant effect.

The two-point autocorrelation function of ultra-high energy cosmic ray (UHECR) arrival directions has a broad maximum around 25 degrees, combining the data with energies above $4\times 10^{19}$ eV (in the HiRes energy scale) of the HiRes stereo, AGASA, Yakutsk and SUGAR experiments. This signal is not or only marginally present analyzing events of a single experiment, but becomes significant when data from several experiments are added. Both the energy dependence of the signal and its angular scale might be interpreted as first signatures of the large-scale structure of UHECR sources and of intervening magnetic fields.

The energy spectrum of ultra-high energy cosmic rays (UHECR) is usually calculated for sources with identical properties. Assuming that all sources can accelerate UHECR protons to the same extremely high maximal energy E_max > 10^{20} eV and have the steeply falling injection spectrum 1/E^{2.7}, one can reproduce the measured cosmic ray flux above E > 10^{18} eV. We show that relaxing the assumption of identical sources and using a power-law distribution of their maximal energy allows one to explain the observed UHECR spectrum with the injection 1/E^2 predicted by Fermi shock acceleration.

Brane inflation predicts the production of cosmic superstrings with tension 10^{-12}<G\mu<10^{-7}. Superstring theory predicts also the existence of a dilaton with a mass that is at most of the order of the gravitino mass. We show that the emission of dilatons imposes severe constraints on the allowed evolution of a cosmic superstring network. In particular, the detection of gravitational wave burst from cosmic superstrings by LIGO is only possible if the typical length of string loops is much smaller than usually assumed.

The observed excess of high-energy cosmic rays from the Galactic plane in the energy range \sim 10^18 eV may be naturally explained by neutron primaries generated in the photo-dissociation of heavy nuclei. In this scenario, neutrons with lower energy decay before reaching the Earth and produce a detectable flux in a 1 km^3 neutrino telescope. The initial flavor composition of these neutrinos, \phi(\bar\nu_e):\phi(\bar\nu_\mu):\phi(\bar\nu_\tau)=1:0:0, offers the opportunity to perform a combined \bar\nu_\mu/\bar\nu_\tau appearance and \bar\nu_e disappearance experiment. The observable ratio \phi(\bar\nu_\mu)/\phi(\bar\nu_e+\bar\nu_\tau) of fluxes arriving on Earth depends appreciably on the 13-mixing angle \theta_13 and the leptonic CP phase \delta_CP, opening thus a new experimental avenue to measure these two quantities.

We have calculated the probability that the clustering of arrival directions of ultra-high energy cosmic rays (UHECRs) is consistent with a finite number of uniformly distributed proton sources. The case of a continuous source distribution is reached only for an unrealisticly high source density, $n_s\gg 10^{-2}/Mpc^3$. Even for densities as large as $n_s=10^{-3}/Mpc^3$, less than half of the observed cluster are on average by chance. For the best-fit value $n_s=(1-3)\times 10^{-5}/Mpc^3$ derived from the AGASA data, the probability that at least one observed cluster is from a true point source is larger than 99.97%, while on average almost all observed clusters are true. The best-fit value found is comparable to the density of AGNs and consistent with the recent HiRes stereo data. In this scenario, the Pierre Auger Observatory will not only establish the clustering of UHECRs but also determine the density of UHECR sources within a factor of a few after one year of data taking.

We study the gravitational radiation from a rotating monopole-antimonopole pair connected by a string. While at not too high frequencies the emitted gravitational spectrum is described asymptotically by $P_n\propto n^{-1}$, the spectrum is exponentially suppressed in the high-frequency limit, $P_n\propto \exp(-n/n_{\rm cr})$. Below $n_{\rm cr}$, the emitted spectrum of gravitational waves is very similar to the case of an oscillating monopole pair connected by a string, and we argue therefore that the spectrum found holds approximately for any moving monopole-string system. As application, we discuss the stochastic gravitational wave background generated by monopole-antimonopole pairs connected by strings in the early Universe and gravitational wave bursts emitted at present by monopole-string networks. We confirm that advanced gravitational wave detectors have the potential to detect a signal for string tensions as small as $G\mu\sim 10^{-13}$.

One of the robust features found in simulations of core-collapse supernovae (SNe) is the prompt neutronization burst, i.e. the first $\sim 25$ milliseconds after bounce when the SN emits with very high luminosity mainly $\nu_e$ neutrinos. We examine the dependence of this burst on variations in the input of current SN models and find that recent improvements of the electron capture rates as well as uncertainties in the nuclear equation of state or a variation of the progenitor mass have only little effect on the signature of the neutronization peak in a megaton water Cherenkov detector for different neutrino mixing schemes. We show that exploiting the time-structure of the neutronization peak allows one to identify the case of a normal mass hierarchy and large 13-mixing angle $\theta_{13}$, where the peak is absent. The robustness of the predicted total event number in the neutronization burst makes a measurement of the distance to the SN feasible with a precision of about 5%, even in the likely case that the SN is optically obscured.

A few seconds after bounce in a core-collapse supernova, the shock wave passes the density region corresponding to resonant neutrino oscillations with the ``atmospheric'' neutrino mass difference. The transient violation of the adiabaticity condition manifests itself in an observable modulation of the neutrino signal from a future galactic supernova. In addition to the shock wave propagation effects that were previously studied, a reverse shock forms when the supersonically expanding neutrino-driven wind collides with the slower earlier supernova ejecta. This implies that for some period the neutrinos pass two subsequent density discontinuities, giving rise to a ``double dip'' feature in the average neutrino energy as a function of time. We study this effect both analytically and numerically and find that it allows one to trace the positions of the forward and reverse shocks. We show that the energy dependent neutrino conversion probabilities allow one to detect oscillations even if the energy spectra of different neutrino flavors are the same as long as the fluxes differ. These features are observable in the \bar\nu_e signal for an inverted and in the \nu_e signal for a normal neutrino mass hierarchy, provided the 13-mixing angle is ``large'' (sin^2\theta_{13}\gg 10^{-5}).

In this work we present a detailed computation of the spectra of UHECR in the top-down scenario. We compare the spectra of hadrons obtained by two different methods in QCD and supersymmetric (SUSY) QCD with large primary energies $\sqrt{s}$ up to $10^{16}$ GeV. The two methods discussed are a Monte Carlo (MC) simulation and the evolution of the hadron fragmentation functions as described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. The hadron spectra obtained by the two methods agree fairly well in the interesting energy range $10^{-5}M_{X}<E<0.3M_{X}$ ($M_X$ is the energy scale of the process $M_{X}\ge 10^{12}$ GeV). We have also computed the spectra of photons, neutrinos and nucleons obtaining a good agreement with other published results. The consistency of the spectra computed by different methods allows us to consider the spectral shape as a signature of the production model for UHECR, such as the decay of super heavy relic particles or topological defects.

The status of solutions to the ultra-high energy cosmic ray puzzle that involve particle physics beyond the standard model is reviewed. Signatures and experimental constraints are discussed for most proposals like the Z burst model and topological defects (both allowed only as subdominant contribution), supermassive dark matter (no positive evidence from its key signatures galactic anisotropy and photon dominance), strongly interacting neutrinos or new primaries (no viable models known), and violation of Lorentz invariance (possible).

We present results from two different methods for the calculation of hadron spectra in QCD and SUSY QCD with large primary energies $\sqrt{s}$ up to $10^{16}$ GeV. The two methods considered are a Monte Carlo (MC) simulation and the evolution of fragmentation functions described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. We find that the pion, nucleon and all-hadron spectra calculated with the two methods agree well. The MC simulation is performed with new hadronization functions (in comparison with our previous work), motivated by low-energy ($\sqrt s < M_Z$) data and DGLAP. The hadron spectra calculated with both sets of hadronization functions agree well, which indicates that our method for calculating the hadronization function works successfully. The small difference in the calculated hadron spectra characterizes the uncertainties of this method. We calculate also the spectra of photons, neutrinos and nucleons and compare them with other published results. The agreement is good for all $x$ from $\sim 10^{-5}$ up to $x\leq 0.3$. The consistency of the spectra calculated by different methods allows to consider the spectral shape as a signature of models with decays or annihilations of superheavy particles, such as topological defects or superheavy DM. The UHECR spectra from these sources are calculated.

The Earth matter effects on supernova (SN) neutrinos can be identified at a single detector through peaks in the Fourier transform of their ``inverse energy'' spectrum. The positions of these peaks are independent of the SN models and therefore the peaks can be used as a robust signature of the Earth matter effects, which in turn can distinguish between different neutrino mixing scenarios. Whereas only one genuine peak is observable when the neutrinos traverse only the Earth mantle, traversing also the core gives rise to multiple peaks. We calculate the strengths and positions of these peaks analytically and explore their features at a large scintillation detector as well as at a megaton water Cherenkov detector through Monte Carlo simulations. We propose a simple algorithm to identify the peaks in the actual data and quantify the chances of a peak identification as a function of the location of the SN in the sky.

A future galactic SN can be located several hours before the optical explosion through the MeV-neutrino burst, exploiting the directionality of $\nu$-$e$-scattering in a water Cherenkov detector such as Super-Kamiokande. We study the statistical efficiency of different methods for extracting the SN direction and identify a simple approach that is nearly optimal, yet independent of the exact SN neutrino spectra. We use this method to quantify the increase in the pointing accuracy by the addition of gadolinium to water, which tags neutrons from the inverse beta decay background. We also study the dependence of the pointing accuracy on neutrino mixing scenarios and initial spectra. We find that in the ``worst case'' scenario the pointing accuracy is $8^\circ$ at 95% C.L. in the absence of tagging, which improves to $3^\circ$ with a tagging efficiency of 95%. At a megaton detector, this accuracy can be as good as $0.6^\circ$. A TeV-neutrino burst is also expected to be emitted contemporaneously with the SN optical explosion, which may locate the SN to within a few tenths of a degree at a future km$^2$ high-energy neutrino telescope. If the SN is not seen in the electromagnetic spectrum, locating it in the sky through neutrinos is crucial for identifying the Earth matter effects on SN neutrino oscillations.

Decay or annihilation products of superheavy dark matter (SHDM) could be responsible for the end of the Ultra-High Energy Cosmic Ray (UHECR) spectrum. In this case, the south array of the Pierre Auger Observatory should observe in the future a significant anisotropy of UHECR arrival directions towards the galactic center. Here we use the already existing data of the SUGAR array to test this possibility. If decaying SHDM is distributed according a Navarro-Frenk-White (NFW) dark matter profile with core radius R_c=15 kpc and is responsible only for UHECRs above 6*10^{19} eV, i.e. the AGASA excess, then the arrival directions measured by the SUGAR array have a probability of ~10% to be consistent with this model. By contrast, the model of annihilating SHDM is disfavoured at least at 99% CL by the SUGAR data, if the smooth component of the DM dominates the signal.

Ultra-high energy cosmic ray (UHECR) protons produced by uniformly distributed astrophysical sources contradict the energy spectrum measured by both the AGASA and HiRes experiments, assuming the small scale clustering of UHECR observed by AGASA is caused by point-like sources. In that case, the small number of sources leads to a sharp exponential cutoff at the energy E<10^{20} eV in the UHECR spectrum. New hadrons with mass 1.5-3 GeV can solve this cutoff problem. For the first time we discuss the production of such hadrons in proton collisions with infrared/optical photons in astrophysical sources. This production mechanism, in contrast to proton-proton collisions, requires the acceleration of protons only to energies E<10^{21} eV. The diffuse gamma-ray and neutrino fluxes in this model obey all existing experimental limits. We predict large UHE neutrino fluxes well above the sensitivity of the next generation of high-energy neutrino experiments. As an example we study hadrons containing a light bottom squark. These models can be tested by accelerator experiments, UHECR observatories and neutrino telescopes.

We study neutrino oscillations and the level-crossing probability PLZ in power-law potential profiles A(r)\propto r^n. We give local and global adiabaticity conditions valid for all mixing angles theta and discuss different representations for PLZ. For the 1/r^3 profile typical of supernova envelopes we compare our analytical to numerical results and to earlier approximations used in the literature. We then perform a combined likelihood analysis of the observed SN1987A neutrino signal and of the latest solar neutrino data, including the recent SNO CC measurement. We find that, unless all relevant supernova parameters (released binding energy, \bar\nu_e and \bar\nu_{\mu,\tau} temperatures) are near their lowest values found in simulations, the status of large mixing type solutions deteriorates considerably compared to fits using only solar data. This is sufficient to rule out the vacuum-type solutions for most reasonable choices of astrophysics parameters. The LOW solution may still be acceptable, but becomes worse than the SMA-MSW solution which may, in some cases, be the best combined solution. On the other hand the LMA-MSW solution can easily survive as the best overall solution, although its size is generally reduced when compared to fits to the solar data only.

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.

The current status of the ultrahigh energy cosmic ray (UHE CR) enigma and several proposed solutions involving particle physics beyond the standard model are discussed. Emphasis is given to top--down models, and as a main example, supermassive dark matter as galactic source for UHE CR and the status of its experimental signatures (galactic anisotropy, chemical composition and clustering) is reviewed. Then different approaches to calculate fragmentation spectra of supermassive particles are discussed. Finally, it is argued that UHE neutrinos cannot be - neither directly or indirectly - responsible for the observed vertical air showers.

It has been suggested that ultrahigh energy neutrinos can acquire cross-sections approaching hadronic size if the string scale is as low as 1-10 TeV. In this case, the vertical air showers observed with energies above the Greisen-Zatsepin-Kuzmin cutoff at E approximately 6x10^{19} eV could be initiated by neutrinos which are the only known primaries able to travel long distances unimpeded. We have calculated the neutrino-nucleon cross-section due to the exchange of Kaluza-Klein excitations of the graviton in a field theoretical framework. We have found that the neutrino-nucleon cross section and the transferred energy per interaction are too small to explain vertical showers even in the most optimistic scenario.

We present results from a new Monte Carlo simulation for jet fragmentation in QCD and SUSY QCD for large primary energies $\sqrt s$ up to $10^{16}$ GeV. In the case of SUSY QCD the simulation takes into account not only gluons and quarks as cascading particles, but also their supersymmetric partners. A new model-independent hadronization scheme is developed, in which the hadronization functions are found from LEP data. An interesting feature of SUSY QCD is the prediction of a sizeable flux of the lightest supersymmetric particles (LSPs), if R-parity is conserved. About 10% of the jet energy is transferred to LSPs which, owing to their harder spectra, constitute an important part of the spectra for large $x=E/E_{jet}$. Spectra of protons and of secondary particles, photons and neutrinos, are also calculated. These results have implications for the decay of superheavy particles with masses up to the GUT scale, which have been suggested as a source of ultrahigh energy cosmic rays.

withdrawn - wrong selection rules used.

The supersymmetric generalization of the limiting and Gaussian QCD spectra is obtained. These spectra are valid for $x \ll 1$, when the main contribution to the parton cascade is given by gluons and gluinos. The derived spectra are applied to decaying superheavy particles with masses up to the GUT scale. These particles can be relics from the Big Bang or produced by topological defects and could give rise to the observed ultrahigh energy cosmic rays. General formulae for the fluxes of protons, photons and neutrinos due to decays of superheavy particles are obtained.

Recently, models proposing superheavy particles $X$ as source of ultrahigh energy cosmic rays have attracted some interest. The $X$-particles are either metastable relic particles from the early Universe or are released by topological defects. In these models, the detected air-showers are produced by primaries originating from the fragmentation of the $X$-particles. We present the fragmentation spectrum of superheavy particles calculated in SUSY-QCD. Then we discuss the status of the lightest supersymmetric particle as possible ultrahigh energy primary.

We argue that the lightest supersymmetric particles (LSP) can be produced with extremely high energies $E\gsim 10^{10}$ GeV in the Universe at the present epoch. Their most probable sources are decaying superheavy particles produced by topological defects or as relic Big Bang particle. We discuss the mechanisms of production of LSP at ultra-high energies (UHE) and the interaction of the UHE LSP with matter. The most attention is given to the neutralino as LSP, although the gluino is also considered as a phenomenological possibility.

We considered recently as a new axion production mechanism the process $e^- \to e^- +a$ in a strong magnetic field $B$. Requiring that for a strongly magnetized neutron star the axion luminosity is smaller than the neutrino luminosity we obtained the bound $g_{ae} \lsim 10^{-10}$ for the axion electron coupling constant. This limit is considerably weaker than the bound derived earlier by Borisov and Grishina using the same method. Applying a similar argument to magnetic white dwarf stars we obtained $g_{ae}\lsim 9 \cdot 10^{-13} (T/10^7)^{5/4} (B/10^{10})^{-2}$, where $T$ is the internal temperature of the white dwarf. Here we note that the observed lack of magnetized white dwarfs with low-temperature in the galactic disc could also be interpreted as a signature of axion emission. Moreover, we speculate that axion emission could explain why the putative galactic halo population of white dwarfs is so dim.

We study the decays of ultraheavy (m_X \geq 10^{13}GeV) and quasistable (lifetime \tau_X much larger than the age of the Universe t_0) particles as the source of Ultra High Energy Cosmic Rays (UHE CR). These particles are assumed to constitute a tiny fraction \xi_X of CDM in the Universe, with \xi_X being the same in the halo of our Galaxy and in the intergalactic space. The elementary-particle and cosmological scenarios for these particles are briefly outlined. The UHE CR fluxes produced at the decays of X- particles are calculated. The dominant contribution is given by fluxes of photons and nucleons from the halo of our Galaxy and thus they do not exibit the GZK cutoff. The extragalactic components of UHE CR are suppressed by the smaller extragalactic density of X-particles and hence the cascade limit is relaxed. We discuss the spectrum of produced Extensive Air Showers (EAS) and a signal from Virgo cluster as signatures of this model.

The energy loss rate of a magnetized electron gas emitting axions a due to the process $e^- \to e^- +a$ is derived for arbitrary magnetic field strength B. Requiring that for a strongly magnetized neutron star the axion luminosity is smaller than the neutrino luminosity we obtain the bound $g_{ae}\lsim 10^{-10}$ for the axion electron coupling constant. This limit is considerably weaker than the bound derived earlier by Borisov and Grishina using the same method. Applying a similar argument to magnetic white dwarf stars results in the more stringent bound $g_{ae}\lsim 9x10^{-13} (T/10^7 K)^{5/4} (B/10^{10} G)^{-2}$ where T is the internal temperature of the white dwarf.

We question the use of stable asymptotic scattering states in QED of strong magnetic fields. To correctly describe excited Landau states and photons above the pair creation threshold the asymptotic fields are chosen as generalized Licht fields. In this way the off-shell behavior of unstable particles is automatically taken into account, and the resonant divergences that occur in scattering cross sections in the presence of a strong external magnetic field are avoided. While in a limiting case the conventional electron propagator with Breit-Wigner form is obtained, in this formalism it is also possible to calculate $S$-matrix elements with external unstable particles.

We study pion curvature radiation by a proton, i.e. pion emission by a proton moving along a curved trajectory. We suggest an approximate semiclassical solution and the exact solution for which we assume that a proton moves in a fictitious magnetic field with the Larmor radius equal to the curvature radius of the real trajectory. As possible application we consider the pion radiation by ultrahigh energy protons moving along curved magnetic field lines. Such situation can occur in the magnetosphere of a young pulsar, in the magnetosphere of the accretion disk around a black hole, and in the vicinity of a superconducting cosmic string. The decay products of these pions, such as high energy photons or neutrinos, can give the observable consequences of the considered mechanism.

The infrared behavior of QED changes drastically in the presence of a strong magnetic field: the electron self-energy and the vertex function are infrared {\em finite}, in contrast with field-free QED, while new infrared divergences appear that are absent in free space. One famous example of the latter is the infrared catastrophe of magnetic Compton scattering, where the cross section for scattering of photons from electrons which undergo a transition to the Landau ground state {\em diverges} as the frequency of the incoming photon goes to zero. We examine this divergence in more detail and prove that the singularity of the cross section is {\em removed} as soon as proper account is taken of all quantum electrodynamical processes that become indistinguishable from Compton scattering in the limit of vanishing frequency of the incident photon.