Source author record

M. Kachelriess

M. Kachelriess appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2022arXiv

Extragalactic cosmic rays

I review the status of ultrahigh-energy cosmic ray (UHECR) physics.After introducing the main experimental results and summarizing possible intepretations, I discuss observational and theoretical constraints on the sources of UHECRs. I comment also briefly on the role of magnetic fields. Combining these constraints, I argue that luminuous and numerous AGN types as FR-I and Seyfert galaxies, or alternatively hypernovae, are the most promising UHECR sources. Finally, I sketch few of the models presented at the conference before concluding.

preprint2022arXiv

On the origin and the detection of characteristic axion wiggles in photon spectra

Photons propagating in an external magnetic field may oscillate into axions or axion-like particles (ALPs). Such oscillations will lead to characteristic features in the energy spectrum of high-energy photons from astrophysical sources that can be used to probe the existence of ALPs. In this work, we revisit the signatures of these oscillations and stress the importance of a proper treatment of turbulent magnetic fields. We implement axions into ELMAG, complementing thereby the usual description of photon-axion oscillations with a Monte Carlo treatment of high-energy photon propagation and interactions. We also propose an alternative method of detecting axions through the discrete power spectrum using as observable the energy dependence of wiggles in the photon spectra.

preprint2020arXiv

Cosmic-ray Antinuclei as Messengers of New Physics: Status and Outlook for the New Decade

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.

preprint2020arXiv

ELMAG 3.01: A three-dimensional Monte Carlo simulation of electromagnetic cascades on the extragalactic background light and in magnetic fields

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.

preprint2020arXiv

High-energy neutrinos from cosmic ray interactions in the Local Bubble

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.

preprint2020arXiv

Revisiting cosmic ray antinuclei fluxes with a new coalescence model

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 $σ$ that we tune to experimental data on antideuteron production in electron-positron, proton-proton and proton-nucleus collisions. The obtained value $σ\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.

preprint2019arXiv

AAfrag: Interpolation routines for Monte Carlo results on secondary production in proton-proton, proton-nucleus and nucleus-nucleus interactions

We provide a compilation of predictions of the QGSJET-II-04m model for the production of secondary species (photons, neutrinos, electrons, positrons, and antinucleons) that are covering a wide range of energies of the beam particles in proton-proton, proton-nucleus, nucleus-proton, and nucleus-nucleus reactions. The current version of QGSJET-II-04m has an improved treatment of the production of secondary particles at low energies: the parameters of the hadronization procedure have been fine-tuned, based on a number of recent benchmark experimental data, notably, from the LHCf, LHCb, and NA61 experiments. Our results for the production spectra are made publicly accessible through the interpolation routines AAfrag which are described below. Besides, we comment on the impact of Feynman scaling violation and isospin symmetry effects on antinucleon production.

preprint2019arXiv

Alternative coalescence model for deuteron, tritium, helium-3 and their antinuclei

Antideuteron and antihelium nuclei have been proposed as a detection channel for dark matter annihilations and decays in the Milky Way, due to the low astrophysical background expected. To estimate both the signal for various dark matter models and the astrophysical background, one employs usually the coalescence model in a Monte Carlo framework. This allows one to treat the production of antinuclei on an event-by-event basis, taking thereby into account momentum correlations between the antinucleons involved in the process. This approach lacks however an underlying microscopic picture, and the numerical value of the coalescence parameter obtained from fits to different reactions varies considerably. Here we propose instead to combine event-by-event Monte Carlo simulations with a microscopic coalescence picture based on the Wigner function representations of the produced antinuclei states. This approach allows us to include in a semi-classical picture both the size of the formation region, which is process dependent, and momentum correlations. The model contains a single, universal parameter which is fixed by fitting the production spectra of antideuterons in proton-proton interactions, measured at the Large Hadron Collider. Using this value, the model describes well the production of various antinuclei both in electron-positron annihilation and in proton-proton collisions.

preprint2019arXiv

Cosmic Ray Models

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.

preprint2015arXiv

Escape model for Galactic cosmic rays and an early extragalactic transition

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 $τ_{\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.

preprint2015arXiv

Imprint of a 2 Myr old source on the cosmic ray anisotropy

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.

preprint2015arXiv

Signatures of a two million year old supernova in the spectra of cosmic ray protons, antiprotons and positrons

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.

preprint2015arXiv

Unified model for cosmic rays above $10^{17}$ eV and the diffuse gamma-ray and neutrino backgrounds

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 $γ$-ray fluxes produced by these CR protons interacting with gas inside their sources. For a spectral slope of CRs close to $α=2.1-2.2$ as suggested by shock acceleration, we find that these UHECR sources contribute the dominant fraction of both the isotropic $γ$-ray background and of the extragalactic part of the astrophysical neutrino signal observed by IceCube.

preprint2014arXiv

Explaining the Spectra of Cosmic Ray Groups above the Knee by Escape from the Galaxy

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 $τ_{\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 $τ_{\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.

preprint2014arXiv

Neutrino yield from Galactic cosmic rays

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.

preprint2013arXiv

Anisotropic Cosmic Ray Diffusion and its Implications for Gamma-Ray Astronomy

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.

preprint2013arXiv

B/C ratio and the PAMELA positron excess

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.

preprint2012arXiv

Antideuterons from dark matter annihilations and hadronization model dependence

We calculate the antideuteron yield in dark matter annihilations on an event-by-event basis using the HERWIG++ Monte Carlo generator. We present the resulting antideuteron fluxes for quark and gauge boson final states. As deuteron production in the coalescence model depends on momentum differences between nucleons that are small compared to Λ_{QCD}, it is potentially very sensitive to the hadronization model employed. We therefore compare our antideuteron yields to earlier results based on PYTHIA, thereby estimating their uncertainties. We also briefly discuss the importance of n>2 final states for annihilations of heavy DM particles.

preprint2012arXiv

Cosmic Ray Anisotropy as Signature for the Transition from Galactic to Extragalactic Cosmic Rays

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.

preprint2012arXiv

Degree-scale GeV "jets" from active and dead TeV blazars

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.

preprint2012arXiv

Deriving the cosmic ray spectrum from gamma-ray observations

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.

preprint2012arXiv

ELMAG: A Monte Carlo simulation of electromagnetic cascades on the extragalactic background light and in magnetic fields

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.

preprint2012arXiv

Filamentary Diffusion of Cosmic Rays on Small Scales

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.

preprint2011arXiv

Antimatter production in supernova remnants

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.

preprint2011arXiv

Propagation of Ultrahigh Energy Nuclei in the Magnetic Field of our Galaxy

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.

preprint2011arXiv

Ultrahigh Energy Nuclei in the Turbulent Galactic Magnetic Field

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.

preprint2010arXiv

Annihilations of superheavy dark matter in superdense clumps

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

preprint2010arXiv

Antimatter spectra from a time-dependent modeling of supernova remnants

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.

preprint2010arXiv

Lower limit on the strength and filling factor of extragalactic magnetic fields

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.

preprint2010arXiv

Restricting UHECRs and cosmogenic neutrinos with Fermi-LAT

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.

preprint2010arXiv

Superdense cosmological dark matter clumps

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 $ρ(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 $ρ\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.

preprint2010arXiv

TeV gamma-rays from UHECR interactions in AGN cores: Lessons from Centaurus A

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.

preprint2010arXiv

Ultrahigh Energy Nuclei in the Galactic Magnetic Field

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.

preprint2009arXiv

On the role of electroweak bremsstrahlung for indirect dark matter signatures

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.

preprint2008arXiv

Supersymmetric superheavy dark matter

We propose the lightest supersymmetric particle (LSP) as a well-suited candidate for superheavy dark matter (SHDM). Various production mechanisms at the end of inflation can produce SHDM with the correct abundance, $Ω_{LSP} h^2 \sim 0.1$, if its mass is sufficiently high. In particular, gravitational production requires that the mass $m_{LSP}$ of the LSP is above $3\times 10^{11} GeV$. Weak interactions remain perturbative despite the large mass hierarchy, $m_{LSP}\gg m_Z$, because of the special decoupling properties of supersymmetry. As a result the model is predictive and we discuss the relevant cosmological processes for the case of a superheavy neutralino within this scheme.

preprint2008arXiv

UHECR observations and lensing in the magnetic field of the Virgo cluster

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.

preprint2006arXiv

On the status of superheavy dark matter

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.

preprint2005arXiv

Exploiting the neutronization burst of a galactic supernova

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 $ν_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 $θ_{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.

preprint2005arXiv

Gravitational radiation from rotating monopole-string systems

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μ\sim 10^{-13}$.

preprint2005arXiv

Measuring the 13-mixing angle and the CP phase with neutrino telescopes

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, ϕ(\barν_e):ϕ(\barν_μ):ϕ(\barν_τ)=1:0:0, offers the opportunity to perform a combined \barν_μ/\barν_τappearance and \barν_e disappearance experiment. The observable ratio ϕ(\barν_μ)/ϕ(\barν_e+\barν_τ) of fluxes arriving on Earth depends appreciably on the 13-mixing angle θ_13 and the leptonic CP phase δ_CP, opening thus a new experimental avenue to measure these two quantities.

preprint2004arXiv

Fragmentation functions in SUSY QCD and UHECR spectra produced in top-down models

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.

preprint2003arXiv

New hadrons as ultra-high energy cosmic rays

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.

preprint2003arXiv

Superheavy dark matter as UHECR source versus the SUGAR data

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.

preprint1998arXiv

Limiting SUSY-QCD spectrum and its application for decays of superheavy particles

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.

preprint1997arXiv

Ultra-high energy cosmic rays without GZK cutoff

We study the decays of ultraheavy (m_X \geq 10^{13}GeV) and quasistable (lifetime τ_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 ξ_X of CDM in the Universe, with ξ_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.

preprint1996arXiv

Radiative neutrino decays in very strong magnetic fields

The radiative decay $ν_H \to ν_L + γ$ of massive neutrinos is analyzed in the framework of the standard model with lepton mixing for very strong magnetic fields $B\gg B_{cr} = m_e^2/e \sim 4.14 \times 10^{13}$G. The analysis is based on the approximate decay amplitude obtained by Gvozdev et al. Numerical results as well as analytical approximations for the decay rate are obtained for energies of the initial neutrino below and {\em above} the electron-positron pair creation threshold $2m_e$.

preprint1995arXiv

Curvature Radiation by Ultrarelativistic Protons

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.