Showing 1–50 of 107 results for author: Kachelriess, M

Galactic Distribution of Supernovae and OB Associations

Authors: M. Kachelriess, V. Mikalsen

Abstract: We update and extend a previous model by Higdon and Lingenfelter for the longitudinal profile of the N\,II intensity in the Galactic plane. The model is based on four logarithmic spiral arms, to which features like the Local Arm and local sources are added. Connecting then the N\,II to the H\,II emission, we use this model to determine the average spatial distribution of OBassociations in the Milk… ▽ More We update and extend a previous model by Higdon and Lingenfelter for the longitudinal profile of the N\,II intensity in the Galactic plane. The model is based on four logarithmic spiral arms, to which features like the Local Arm and local sources are added. Connecting then the N\,II to the H\,II emission, we use this model to determine the average spatial distribution of OBassociations in the Milky Way. Combined with a stellar mass and cluster distribution function, the model predicts the average spatial and temporal distribution of core-collapse supernovae in the Milky Way. In addition to this average population, we account for supernovae from observed OB associations, providing thereby a more accurate description of the nearby Galaxy. The complete model is made publicly available in the python code SNOB. △ Less

Submitted 11 September, 2024; originally announced September 2024.

Comments: 21 pages, 9 pdf figures

Polarized synchrotron data and the structure of the Galactic magnetic field

Authors: C. T. Becker, M. Kachelriess

Abstract: The polarized synchrotron data of the northern C-BASS survey show a surprisingly low linear polarization fraction, $Πぱい\simeq 3\%$, while the magnetic field polarization features coherent structures over rather large angular scales. The low polarization degree points to a strong dominance of the turbulent magnetic field -- in agreement with the observation that the total synchrotron intensity is muc… ▽ More The polarized synchrotron data of the northern C-BASS survey show a surprisingly low linear polarization fraction, $Πぱい\simeq 3\%$, while the magnetic field polarization features coherent structures over rather large angular scales. The low polarization degree points to a strong dominance of the turbulent magnetic field -- in agreement with the observation that the total synchrotron intensity is much larger than expected from the regular component in current models for the Galactic magnetic field (GMF). In contrast, studies of cosmic ray propagation employing these GMF models suggest that cosmic rays propagate anisotropically, what in turn requires weak turbulent fields. As a solution to this contradicting requirements, we suggest that the GMF consists of three components with different levels of turbulence: a disk dominated by the turbulent field, a halo dominated by the regular field, plus an extended turbulent halo field. △ Less

Submitted 5 August, 2024; originally announced August 2024.

Comments: 9 pages, 7 pdf figures

Energy dependence of the knee in the cosmic ray spectrum across the Milky Way

Authors: C. Prevotat, M. Kachelriess, S. Koldobskiy, A. Neronov, D. Semikoz

Abstract: The all-particle spectrum of cosmic rays measured at Earth has a knee-like feature around 4 PeV. A priori, it is not clear if this is a local feature specific to the Solar neighbourhood in the Milky Way, or if it is a generic property of the galactic cosmic ray spectrum. We argue that combining gamma-ray and cosmic-ray data of LHAASO indicates that the knee is a local feature. In order to demonstr… ▽ More The all-particle spectrum of cosmic rays measured at Earth has a knee-like feature around 4 PeV. A priori, it is not clear if this is a local feature specific to the Solar neighbourhood in the Milky Way, or if it is a generic property of the galactic cosmic ray spectrum. We argue that combining gamma-ray and cosmic-ray data of LHAASO indicates that the knee is a local feature. In order to demonstrate this, we derive a model for the local cosmic-ray spectra and composition, consistent with the recent LHAASO measurements of the all-particle spectrum and the mean logarithmic mass in the knee region. We calculate the spectrum of diffuse gamma-ray emission based on this model and find that the expected spectral shape of the diffuse gamma-ray flux is not consistent with LHAASO measurements of the diffuse gamma-ray emission in the 10-100 TeV energy range. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 8 pages, 3 figures, 3 tables

Lunar response to gravitational waves

Authors: M. Kachelriess, M. P. Nødtvedt

Abstract: It has been suggested to use seismic detectors on the Moon as a tool to search for gravitational waves in an intermediate frequency range between mHz and Hz. Employing three different spherically symmetric models for the lunar interior, we investigate the response of the Moon to gravitational waves in Einstein and Jordan-Brans-Dicke gravity. We find that the first eigenfrequencies of the different… ▽ More It has been suggested to use seismic detectors on the Moon as a tool to search for gravitational waves in an intermediate frequency range between mHz and Hz. Employing three different spherically symmetric models for the lunar interior, we investigate the response of the Moon to gravitational waves in Einstein and Jordan-Brans-Dicke gravity. We find that the first eigenfrequencies of the different models depend only weakly on the model details, with the fundamental frequency $νにゅー_1$ close to 1\,ms both for spheroidal and toroidal oscillations. In contrast, the resulting displacement varies up to a factor two, being in the range $(2.7-5.6)\times 10^{11}/h_0$ cm for spheroidal oscillations with amplitude $h_0$. Toroidal oscillations are suppressed by a factor $2πぱいνにゅーR/c$, both in Einstein gravity and in general scalar-tensor theories. △ Less

Submitted 23 August, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Comments: 18 pages, 17 pdf figures; v2: extended discussion of scalar mode, corrected normalisation of Green function; conclusion unchanged; accepted in PRD

Detecting ALP wiggles at TeV energies

Authors: M. Kachelriess, J. Tjemsland

Abstract: Axions and axion-like-particles (ALPs) are characterised by their two-photon coupling, which entails so-called photon-ALP oscillations as photons propagate through a magnetic field. These oscillations lead to distinctive signatures in the energy spectrum of high-energy photons from astrophysical sources, allowing one to probe the existence of ALPs. In particular, photon-ALP oscillations will induc… ▽ More Axions and axion-like-particles (ALPs) are characterised by their two-photon coupling, which entails so-called photon-ALP oscillations as photons propagate through a magnetic field. These oscillations lead to distinctive signatures in the energy spectrum of high-energy photons from astrophysical sources, allowing one to probe the existence of ALPs. In particular, photon-ALP oscillations will induce energy dependent oscillatory features, or "ALP wiggles", in the photon spectra. We propose to use the discrete power spectrum to search for ALP wiggles and present a model-independent statistical test. By using PKS 2155-304 as an example, we show that the method has the potential to significantly improve the experimental sensitivities for ALP wiggles, and that the ALP wiggles may be detected using the Cherenkov Telescope Array (CTA) for optimistic values of the photon-ALP coupling constant and the magnetic field. Moreover, we discuss how these sensitivities depend on the modelling of the magnetic field. We find that the use of realistic magnetic field models, due to their larger cosmic variance, substantially enhances detection prospects compared to the use of simplified models. △ Less

Submitted 8 December, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 21 pages, 9 figures; v2: clarifications and discussions added, version to appear in JCAP

Impact of the finite life-time of UHECR sources

Authors: Björn Eichmann, Michael Kachelrieß

Abstract: The observational data on ultrahigh energy cosmic rays (UHECR), in particular their mass composition, show strong indications for extremely hard spectra of individual mass groups of CR nuclei at Earth. In this work, we show that such hard spectra can be the result of the finite life-time of UHECR sources, if a few individual sources dominate the UHECR flux at the highest energies. In this case, ti… ▽ More The observational data on ultrahigh energy cosmic rays (UHECR), in particular their mass composition, show strong indications for extremely hard spectra of individual mass groups of CR nuclei at Earth. In this work, we show that such hard spectra can be the result of the finite life-time of UHECR sources, if a few individual sources dominate the UHECR flux at the highest energies. In this case, time delays induced by deflections in the turbulent extragalactic magnetic field as well as from the diffusive or advective escape from the source environment can suppress low-energy CRs, leading to a steepening of the observed spectrum. Considering radio galaxies as the main source of UHECRs, we discuss the necessary conditions that few individual sources dominate over the total contribution from the bulk of sources that have been active in the past. We provide two proof-of-principle scenarios showing that for a turbulent extragalactic magnetic field with a strength $B$ and a coherence length $l_{\rm coh}$, the life-time of a source at a distance $d_{\rm src}$ should satisfy ${t_{\rm act} \sim \left( B/1\,\text{nG} \right)^2\,\left( d_{\rm src}/10\,\text{Mpc} \right)^2\,\left( l_{\rm coh}/1\,\text{Mpc} \right)\,\text{Myr}}$ to obtain the necessary hardening of the CR spectrum at Earth. △ Less

Submitted 28 February, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Journal ref: JCAP02(2023)053

Explaining the UHECR spectrum, composition and large-scale anisotropies with radio galaxies

Authors: Björn Eichmann, Michael Kachelrieß, Foteini Oikonomou

Abstract: Radio galaxies are promising candidates as the sources of ultrahigh energy cosmic rays (UHECRs). In this work, we examine if the stringent constraints imposed by the dipole and quadropole anisotropies as well as the UHECR spectrum and composition allow that radio galaxies are the dominant extragalactic cosmic ray sources. In order to calculate the UHECR flux emitted by individual radio galaxies, w… ▽ More Radio galaxies are promising candidates as the sources of ultrahigh energy cosmic rays (UHECRs). In this work, we examine if the stringent constraints imposed by the dipole and quadropole anisotropies as well as the UHECR spectrum and composition allow that radio galaxies are the dominant extragalactic cosmic ray sources. In order to calculate the UHECR flux emitted by individual radio galaxies, we constrain their properties using information from the radio-CR correlation and a dynamical evolution model. In addition to the UHECR flux from individual, local sources, we include the diffuse flux emitted by the bulk of non-local radio galaxies based on their radio luminosity distribution. Analyzing the source parameters within a range around their expected properties, we finally determine the configurations of local sources describing well the UHECR spectrum, composition and large-scale anisotropies. We obtain a good description of all data even in the case that we include only a small number of local sources. In particular, we find that scenarios where few sources like Fornax A and Virgo A dominate the flux above the ankle, while low-luminosity radio galaxies contribute an isotropic background dominating below the ankle, provide a good fit to the data. △ Less

Submitted 5 July, 2022; v1 submitted 24 February, 2022; originally announced February 2022.

Journal ref: JCAP07(2022)006

Extragalactic cosmic rays

Authors: M. Kachelriess

Abstract: 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 alt… ▽ More 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. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 1+15+4 pages; invited talk at the ICRC 2021

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

Authors: M. Kachelriess, J. Tjemsland

Abstract: 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 tu… ▽ More 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. △ Less

Submitted 13 January, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 1+25 pages, 8 figures; v2: clarifications and references added, improved readability of figures, minor changes, version to appear in JCAP

Energy spectra of secondaries in proton-proton interactions

Authors: S. Koldobskiy, M. Kachelrieß, A. Lskavyan, A. Neronov, S. Ostapchenko, D. V. Semikoz

Abstract: We compare the predictions of AAfrag for the spectra of secondary photons, neutrinos, electrons, and positrons produced in proton-proton collisions to those of the parameterisations of Kamae et al., Kelner et al., and Kafexhiu et al. We find that the differences in the normalisation of the photon energy spectra reach 20$-$50% at intermediate values of the transferred energy fraction $x$, growing u… ▽ More We compare the predictions of AAfrag for the spectra of secondary photons, neutrinos, electrons, and positrons produced in proton-proton collisions to those of the parameterisations of Kamae et al., Kelner et al., and Kafexhiu et al. We find that the differences in the normalisation of the photon energy spectra reach 20$-$50% at intermediate values of the transferred energy fraction $x$, growing up to a factor of two for $x\rightarrow\,$1, while the differences in the neutrino spectra are even larger. We argue that LHCf results on the forward production of photons favor the use of the QGSJET-II-04m model on which AAfrag is based. The differences in the normalisation have important implications in the context of multi-messenger astronomy, in particular, for the prediction of neutrino fluxes based on gamma-ray flux measurements, or regarding the inference of the cosmic ray spectrum, based on gamma-ray data. We note also that the positron-electron ratio from hadronic interactions increases with energy towards the cutoff, an effect which is missed using the average electron-positron spectrum from Kelner et al. Finally, we describe the publicly available python package aafragpy, which provides the secondary spectra of photons, neutrinos, electrons, and positrons. This package complements the AAfrag results for protons with energies above 4 GeV with previous analytical parameterizations of particle spectra for lower energy protons. △ Less

Submitted 17 December, 2021; v1 submitted 1 October, 2021; originally announced October 2021.

Comments: 8 pages, 10 figures

Journal ref: Phys. Rev. D 104, 123027 (2021)

Using Convolutional Neural Networks for the Helicity Classification of Magnetic Fields

Authors: Nicolò Oreste Pinciroli Vago, Ibrahim A. Hameed, Michael Kachelriess

Abstract: The presence of non-zero helicity in intergalactic magnetic fields is a smoking gun for their primordial origin since they have to be generated by processes that break CP invariance. As an experimental signature for the presence of helical magnetic fields, an estimator $Q$ based on the triple scalar product of the wave-vectors of photons generated in electromagnetic cascades from, e.g., TeV blazar… ▽ More The presence of non-zero helicity in intergalactic magnetic fields is a smoking gun for their primordial origin since they have to be generated by processes that break CP invariance. As an experimental signature for the presence of helical magnetic fields, an estimator $Q$ based on the triple scalar product of the wave-vectors of photons generated in electromagnetic cascades from, e.g., TeV blazars, has been suggested previously. We propose to apply deep learning to helicity classification employing Convolutional Neural Networks and show that this method outperforms the $Q$ estimator. △ Less

Submitted 12 June, 2021; originally announced June 2021.

Comments: 14 pages, extended version of a contribution to the proceedings of the 37.th ICRC 2021

Hot spots in the neutrino flux created by cosmic rays from Cygnus and Vela?

Authors: M. Bouyahiaoui, M. Kachelriess, D. V. Semikoz

Abstract: An analysis of 7.5 years of data in the high-energy starting event sample has been recently published by the IceCube collaboration. The hottest spot in a search for neutrino sources was found far above the Galactic plane and is thus, at first sight, difficult to reconcile with a Galactic origin. In this work, we calculate the cosmic ray (CR) density around nearby, young supernova remnants assuming… ▽ More An analysis of 7.5 years of data in the high-energy starting event sample has been recently published by the IceCube collaboration. The hottest spot in a search for neutrino sources was found far above the Galactic plane and is thus, at first sight, difficult to reconcile with a Galactic origin. In this work, we calculate the cosmic ray (CR) density around nearby, young supernova remnants assuming anisotropic diffusion. Combining the obtained CR densities with the matter distribution deduced from extinction maps, we find two prominent hot spots: The one close to the most significant point in the IceCube search for point sources is created by CRs from the Cygnus loop and has an intensity corresponding to two to four neutrino events. Another, more extended one may be caused by CRs from Vela if CR trajectories are sufficiently disturbed by the magnetic field in the shell around the superbubble Loop I. △ Less

Submitted 30 August, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Comments: v2: 9 pages, 6 figures, extended version, added discussion of a possible extended hot region at low b connected to Vela; to appear in PRD

Journal ref: Phys.Rev.D 106 (2022) 6, 063004

Meson production in air showers and the search for light exotic particles

Authors: M. Kachelriess, J. Tjemsland

Abstract: Decays of mesons produced in cosmic ray induced air showers in Earth's atmosphere can lead to a flux of light exotic particles which can be detected in underground experiments. We evaluate the energy spectra of the light neutral mesons $πぱい^0$, $ηいーた$, $ρろー^0$, $ωおめが$, $φふぁい$ and $J/ψぷさい$ produced in interactions of cosmic ray protons and helium nuclei with air using QCD inspired event generators. Summing up the… ▽ More Decays of mesons produced in cosmic ray induced air showers in Earth's atmosphere can lead to a flux of light exotic particles which can be detected in underground experiments. We evaluate the energy spectra of the light neutral mesons $πぱい^0$, $ηいーた$, $ρろー^0$, $ωおめが$, $φふぁい$ and $J/ψぷさい$ produced in interactions of cosmic ray protons and helium nuclei with air using QCD inspired event generators. Summing up the mesons produced in the individual hadronic interactions of air showers, we obtain the resulting fluxes of undecayed mesons. As an application, we re-consider the case of millicharged particles created in the electromagnetic decay channels of neutral mesons. △ Less

Submitted 22 June, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 13 pages, 6 pdf figures, ancillary files contain tabulated fluxes of undecayed atmospheric mesons; v2: added clarifications, improved readablility of figures, accepted version

Cosmic ray positrons from compact binary millisecond pulsars

Authors: Manuel Linares, Michael Kachelriess

Abstract: A new population of neutron stars has emerged during the last decade: compact binary millisecond pulsars (CBMSPs). Because these pulsars and their companion stars are in tight orbits with typical separations of $10^{11}$ cm, their winds interact strongly forming an intrabinary shock. Electron-positron pairs reaccelerated at the shock can reach energies of about 10 TeV, which makes this new populat… ▽ More A new population of neutron stars has emerged during the last decade: compact binary millisecond pulsars (CBMSPs). Because these pulsars and their companion stars are in tight orbits with typical separations of $10^{11}$ cm, their winds interact strongly forming an intrabinary shock. Electron-positron pairs reaccelerated at the shock can reach energies of about 10 TeV, which makes this new population a potential source of GeV-TeV cosmic ray positrons. We present an analytical model for the fluxes and spectra of positrons from intrabinary shocks of CBMSPs. We find that the minimum energy $E_{\min}$ of the pairs that enter the shock is critical to quantify the energy spectrum with which positrons are injected into the interstellar medium. We measure for the first time the Galactic scale height of CBMSPs, $z_e=0.4\pm0.1$ kpc, after correcting for an observational bias against finding them close to the Galactic plane. From this, we estimate a local density of 5-9 kpc$^{-3}$ and an extrapolated total of 2-7 thousand CBMSPs in the Galaxy. We then propagate the pairs in the isotropic diffusion approximation and find that the positron flux from the total population is about two times higher than that from the 52 currently known systems. For $E_{\min}$ between 1 and 50 GeV, our model predicts only a minor contribution from CBMSPs to the diffuse positron flux at 100 GeV observed at Earth. We also quantify the effects of anisotropic transport due to the ordered Galactic magnetic field, which can change the diffuse flux from nearby sources drastically. Finally, we find that a single "hidden" CBMSP close to the Galactic plane can yield a positron flux comparable to the AMS-02 measurements at 600 GeV if its line-of-sight to Earth is along the ordered Galactic field lines, while its combined electron and positron flux at higher energies would be close to the measurements of CALET, DAMPE and Fermi-LAT. △ Less

Submitted 16 February, 2021; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: Matches version published in JCAP, after minor corrections. Added acknowledgements to the PHAROS COST Action (CA16214) and to the Spanish Ministry of Science, Innovation and Universities (grant CAS19/00228). 31 pages, 10 figures, 5 tables

Journal ref: JCAP02(2021)030

Searching for primordial helical magnetic fields

Authors: M. Kachelriess, B. C. Martinez

Abstract: The presence of non-zero helicity in intergalactic magnetic fields (IGMF) has been suggested as a clear signature for their primordial origin. We extend a previous analysis of diffuse Fermi-LAT gamma-ray data from 2.5 to more than 11 years and show that a hint for helical magnetic fields in the 2.5 year data was a statistical fluctuation. Then we examine the detection prospects of helical magnetic… ▽ More The presence of non-zero helicity in intergalactic magnetic fields (IGMF) has been suggested as a clear signature for their primordial origin. We extend a previous analysis of diffuse Fermi-LAT gamma-ray data from 2.5 to more than 11 years and show that a hint for helical magnetic fields in the 2.5 year data was a statistical fluctuation. Then we examine the detection prospects of helical magnetic fields using individual sources as, e.g., TeV gamma-ray blazars. We find that a detection is challenging employing realistic models for the cascade evolution, the IGMF and the detector resolution in our simulations. △ Less

Submitted 5 November, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: 10 pages, 14 png figures; v2: improved quality of the figures, matches published version

Journal ref: Phys. Rev. D 102, 083001 (2020)

Reacceleration of charged dark matter

Authors: M. Kachelriess, J. Tjemsland

Abstract: Charged particles scattering on moving inhomogenities of the magnetised interstellar medium 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 turbulenc… ▽ More Charged particles scattering on moving inhomogenities of the magnetised interstellar medium 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^3 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. △ Less

Submitted 18 September, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 1+11 pages, 8 pdf figures; v2: references added, version to appear in JCAP

Revisiting cosmic ray antinuclei fluxes with a new coalescence model

Authors: M. Kachelriess, S. Ostapchenko, J. Tjemsland

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 8 July, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 1+21 pages, 14 pdf figures; v2: added additional calculations of antihelium production and few clarifications, accepted version

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

Authors: P. von Doetinchem, K. Perez, T. Aramaki, S. Baker, S. Barwick, R. Bird, M. Boezio, S. E. Boggs, M. Cui, A. Datta, F. Donato, C. Evoli, L. Fabris, L. Fabbietti, E. Ferronato Bueno, N. Fornengo, H. Fuke, C. Gerrity, D. Gomez Coral, C. Hailey, D. Hooper, M. Kachelriess, M. Korsmeier, M. Kozai, R. Lea , et al. (37 additional authors not shown)

Abstract: 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, essential… ▽ More 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. △ Less

Submitted 18 August, 2020; v1 submitted 10 February, 2020; originally announced February 2020.

Comments: 45 pages, 14 figures

Journal ref: JCAP08(2020)035

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

Authors: M. Bouyahiaoui, M. Kachelriess, D. V. Semikoz

Abstract: 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 latitu… ▽ More 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. △ Less

Submitted 24 June, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: 9 pages, 10 figures; v2: added discussion of the importance of other sources

Journal ref: Phys. Rev. D 101, 123023 (2020)

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

Authors: M. Blytt, M. Kachelriess, S. Ostapchenko

Abstract: 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… ▽ More 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. △ Less

Submitted 14 January, 2020; v1 submitted 19 September, 2019; originally announced September 2019.

Comments: 9 pages, 5 figures. v2: minor changes, to appear in CPC; arxiv version contains an addendum describing version 3.02 which includes plasma energy losses; program available at http://elmag.sourceforge.net

Cosmic Ray Models

Authors: M. Kachelriess, D. V. Semikoz

Abstract: 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 Gal… ▽ More 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. △ Less

Submitted 29 August, 2019; v1 submitted 17 April, 2019; originally announced April 2019.

Comments: 89 pages, 56 figures; v2: minor additions, version to appear

Journal ref: Prog. Part. Nucl. Phys. 109, 103710 (2019)

Near-Earth Supernova Explosions: Evidence, Implications, and Opportunities

Authors: Brian D. Fields, John R. Ellis, Walter R. Binns, Dieter Breitschwerdt, Georgia A. de Nolfo, Roland Diehl, Vikram V. Dwarkadas, Adrienne Ertel, Thomas Faestermann, Jenny Feige, Caroline Fitoussi, Priscilla Frisch, David Graham, Brian Haley, Alexander Heger, Wolfgang Hillebrandt, Martin H. Israel, Thomas Janka, Michael Kachelriess, Gunther Korschinek, Marco Limongi, Maria Lugaro, Franciole Marinho, Adrian Melott, Richard A. Mewaldt , et al. (14 additional authors not shown)

Abstract: 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 o… ▽ More 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. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures. Astro2020 Science White Paper submitted to the 2020 Decadal Survey on Astronomy and Astrophysics

Vela as the Source of Galactic Cosmic Rays above 100 TeV

Authors: M. Bouyahiaoui, M. Kachelriess, D. V. Semikoz

Abstract: 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 C… ▽ More 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. △ Less

Submitted 9 December, 2018; originally announced December 2018.

Comments: 11 pages, 6 figures

Anisotropic diffusion and the cosmic ray anisotropy

Authors: M. Kachelriess

Abstract: 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… ▽ More 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. △ Less

Submitted 6 November, 2018; originally announced November 2018.

Comments: 6 pages, 5 eps figures; invited talk at the 26th Extended European Cosmic Ray Symposium, 2018

Heavy decaying dark matter and IceCube high energy neutrinos

Authors: M. Kachelriess, O. E. Kalashev, M. Yu. Kuznetsov

Abstract: 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 νにゅー\barνにゅー$, $X \rightarrow e^+e^-$ and $X \rightarrow q \bar{q}$. The full simulation of ha… ▽ More 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 νにゅー\barνにゅー$, $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. △ Less

Submitted 30 November, 2018; v1 submitted 11 May, 2018; originally announced May 2018.

Comments: 20 pages, 6 figures, journal version

Report number: INR-TH-2018-012

Journal ref: Phys. Rev. D 98, 083016 (2018)

IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays

Authors: Francois Boulanger, Torsten Ensslin, Andrew Fletcher, Philipp Girichides, Stefan Hackstein, Marijke Haverkorn, Joerg R. Hoerandel, Tess R. Jaffe, Jens Jasche, Michael Kachelriess, Kumiko Kotera, Christoph Pfrommer, Jorg P. Rachen, Luiz F. S. Rodrigues, Beatriz Ruiz-Granados, Amit Seta, Anvar Shukurov, Gunter Sigl, Theo Steininger, Valentina Vacca, Ellert van der Velden, Arjen van Vliet, Jiaxin Wang

Abstract: 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 fie… ▽ More 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.) △ Less

Submitted 6 September, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 56 pages, 15 figures, minor updates to match published JCAP version

Journal ref: JCAP08(2018)049

Multi-messenger gamma-ray counterpart of the IceCube neutrino signal

Authors: A. Neronov, M. Kachelriess, D. V. Semikoz

Abstract: 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 th… ▽ More 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. △ Less

Submitted 6 July, 2018; v1 submitted 27 February, 2018; originally announced February 2018.

Journal ref: Phys. Rev. D 98, 023004 (2018)

High-energy neutrinos from Galactic superbubbles

Authors: K. J. Andersen, M. Kachelriess, D. V. Semikoz

Abstract: 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 em… ▽ More 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. △ Less

Submitted 25 June, 2018; v1 submitted 8 December, 2017; originally announced December 2017.

Comments: 5 pages, 3 eps figures; v2: explanations and 2 figures added, to appear in ApJL

Reconciling cosmic ray diffusion with Galactic magnetic field models

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz

Abstract: 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 f… ▽ More 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. △ Less

Submitted 10 July, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: 13 pages, 4 figures; v2: extended discussion of (toy) models, 16 pages, 7 figures, to appear in JCAP

Journal ref: JCAP 07, 051 (2018)

Cosmic ray signatures of a 2-3 Myr old local supernova

Authors: M. Kachelriess, A. Neronov, D. V. Semikoz

Abstract: 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 ti… ▽ More 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). △ Less

Submitted 26 February, 2018; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: v2: added section on electrons, references and extended discussion; to appear in PRD; 10 pages, 7 pdf figures

Journal ref: Phys. Rev. D 97, 063011 (2018)

A minimal model for extragalactic cosmic rays and neutrinos

Authors: M. Kachelriess, O. Kalashev, S. Ostapchenko, D. V. Semikoz

Abstract: 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 fl… ▽ More 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%. △ Less

Submitted 13 September, 2017; v1 submitted 23 April, 2017; originally announced April 2017.

Comments: v2: 7 pages, 6 eps figures; slightly extended version to appear in PRD

Journal ref: Phys. Rev. D 96, 083006 (2017)

Solution of heliospheric propagation: unveiling the local interstellar spectra of cosmic ray species

Authors: M. J. Boschini, S. Della Torre, M. Gervasi, D. Grandi, G. Johannesson, M. Kachelriess, G. La Vacca, N. Masi, I. V. Moskalenko, E. Orlando, S. S. Ostapchenko, S. Pensotti, T. A. Porter, L. Quadrani, P. G. Rancoita, D. Rozza, M. Tacconi

Abstract: 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 leve… ▽ More 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. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 22 pages, 23 figures, 8 tables, and 3 machine-readable tables; ApJ, in press

Journal ref: The Astrophysical Journal, 840:115 (21pp), 2017 May 10

A supernova at 50 pc: Effects on the Earth's atmosphere and biota

Authors: A. L Melott, B. C. Thomas, M. Kachelriess, D. V. Semikoz, A. C. Overholt

Abstract: 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. Assum… ▽ More 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. △ Less

Submitted 7 April, 2017; v1 submitted 14 February, 2017; originally announced February 2017.

Comments: 5 figures. Revised version to be published in ApJ

Journal ref: ApJ 840 105 (2017)

Terrestrial Effects Of Nearby Supernovae In The Early Pleistocene

Authors: Brian C. Thomas, E. E. Engler, M. Kachelrieß, A. L. Melott, A. C. Overholt, D. V. Semikoz

Abstract: Recent results have strongly confirmed that multiple supernovae happened at distances ~100 pc consisting of two main events: one at 1.7 to 3.2 million years ago, and the other at 6.5 to 8.7 million years ago. These events are said to be responsible for excavating the Local Bubble in the interstellar medium and depositing 60Fe on Earth and the Moon. Other events are indicated by effects in the loca… ▽ More Recent results have strongly confirmed that multiple supernovae happened at distances ~100 pc consisting of two main events: one at 1.7 to 3.2 million years ago, and the other at 6.5 to 8.7 million years ago. These events are said to be responsible for excavating the Local Bubble in the interstellar medium and depositing 60Fe on Earth and the Moon. Other events are indicated by effects in the local cosmic ray (CR) spectrum. Given this updated and refined picture, we ask whether such supernovae are expected to have had substantial effects on the terrestrial atmosphere and biota. In a first cut at the most probable cases, combining photon and cosmic ray effects, we find that a supernova at 100 pc can have only a small effect on terrestrial organisms from visible light and that chemical changes such as ozone depletion are weak. However, tropospheric ionization right down to the ground due to the penetration of $\geq$TeV cosmic rays will increase by nearly an order of magnitude for thousands of years, and irradiation by muons on the ground and in the upper ocean will increase 20-fold, which will approximately triple the overall radiation load on terrestrial organisms. Such irradiation has been linked to possible changes in climate and increased cancer and mutation rates. This may be related to a minor mass extinction around the Pliocene-Pleistocene boundary, and further research on the effects is needed. △ Less

Submitted 28 June, 2016; v1 submitted 16 May, 2016; originally announced May 2016.

Comments: Revised version accepted at ApJL

Journal ref: The Astrophysical Journal Letters, Volume 826, Number 1, 2016

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

Authors: G. Giacinti, M. Kachelriess, O. Kalashev, A. Neronov, D. V. Semikoz

Abstract: 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 spect… ▽ More 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. △ Less

Submitted 15 October, 2015; v1 submitted 27 July, 2015; originally announced July 2015.

Comments: 10 pages, 9 eps figures; v2: minor additions, matches version to appear in PRD

Journal ref: Phys. Rev. D 92, 083016 (2015)

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

Authors: V. Savchenko, M. Kachelriess, D. V. Semikoz

Abstract: 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… ▽ More 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. △ Less

Submitted 21 July, 2015; v1 submitted 11 May, 2015; originally announced May 2015.

Comments: 4 pages, 4 eps figures; v2: minor changes, to appear in ApJL

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

Authors: M. Kachelriess, A. Neronov, D. V. Semikoz

Abstract: 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… ▽ More 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. △ Less

Submitted 22 September, 2015; v1 submitted 24 April, 2015; originally announced April 2015.

Comments: 5+2 pages, 3+1 eps figures; v2: reshuffled figures, minor rewording and added supplementary material; to appear in PRL

Journal ref: Phys. Rev. Lett. 115, 181103 (2015)

New calculation of antiproton production by cosmic ray protons and nuclei

Authors: Michael Kachelriess, Igor V. Moskalenko, Sergey S. Ostapchenko

Abstract: 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 appl… ▽ More 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. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: 16 pages (emulateapj style), 18 figures, 2 tables; The Astrophysical Journal, in press

Journal ref: ApJ 2015, 803, 54

Escape model for Galactic cosmic rays and an early extragalactic transition

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz

Abstract: 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 coherenc… ▽ More 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. △ Less

Submitted 6 May, 2015; v1 submitted 5 February, 2015; originally announced February 2015.

Comments: 14 pages, 16 eps figures; v2: minor updates, matches published version

Journal ref: Phys. Rev. D 91, 083009 (2015)

Nuclear enhancement of the photon yield in cosmic ray interactions

Authors: Michael Kachelriess, Igor V. Moskalenko, Sergey S. Ostapchenko

Abstract: 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… ▽ More 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. △ Less

Submitted 30 May, 2014; originally announced June 2014.

Comments: 7 pages (emulateapj style), 1 figure, 5 tables; The Astrophysical Journal, July 2014 volume, in press

Journal ref: ApJ 2014, 789, 136

Neutrino yield from Galactic cosmic rays

Authors: M. Kachelriess, S. Ostapchenko

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2014; v1 submitted 15 May, 2014; originally announced May 2014.

Comments: 9 pages; 7 eps figures; v2: Hillas model corrected, contains bonus material compared to the published version

Report number: NORDITA-2014-47

Journal ref: Phys. Rev. D 90, 083002 (2014)

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

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz

Abstract: 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… ▽ More 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. △ Less

Submitted 5 August, 2014; v1 submitted 13 March, 2014; originally announced March 2014.

Comments: 5 pages, 4 figures; v2: added comparisons with energy spectra of KASCADE (Grande), to be published

Journal ref: Phys. Rev. D 90, 041302(R) (2014)

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

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz

Abstract: 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… ▽ More 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. △ Less

Submitted 26 July, 2013; v1 submitted 13 June, 2013; originally announced June 2013.

Comments: 13 pages (2 columns), 9 figures. Published in Physical Review D

Journal ref: Phys. Rev. D 88, 023010 (2013)

B/C ratio and the PAMELA positron excess

Authors: M. Kachelriess, S. Ostapchenko

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 1 March, 2013; v1 submitted 5 November, 2012; originally announced November 2012.

Comments: 4 pages, 4 eps figures, v2: one reference added

Journal ref: Phys. Rev. D 87, 047301 (2013)

Deriving the cosmic ray spectrum from gamma-ray observations

Authors: M. Kachelriess, S. Ostapchenko

Abstract: 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… ▽ More 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. △ Less

Submitted 24 July, 2012; v1 submitted 20 June, 2012; originally announced June 2012.

Comments: 10 pages; v2: matches version to be published, plus 2 appendices. Photon and antiproton fragmentation functions available from http://sourceforge.net/projects/ppfrag

Filamentary Diffusion of Cosmic Rays on Small Scales

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz

Abstract: 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 t… ▽ More 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. △ Less

Submitted 28 June, 2012; v1 submitted 5 April, 2012; originally announced April 2012.

Comments: 5 pages (2 columns), 4 figures. Published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 108, 261101 (2012)

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

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz, G. Sigl

Abstract: 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 turbul… ▽ More 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. △ Less

Submitted 16 July, 2012; v1 submitted 23 December, 2011; originally announced December 2011.

Comments: 22 pages, 10 figures. Published in JCAP

Journal ref: JCAP 1207, 031 (2012)

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

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz, G. Sigl

Abstract: 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 th… ▽ More 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. △ Less

Submitted 5 September, 2011; v1 submitted 1 August, 2011; originally announced August 2011.

Comments: 4 pages (2 columns), 5 figures. To appear in Proceedings of the 32nd International Cosmic Ray Conference, Beijing, China, 11-18 August 2011

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

Authors: M. Kachelriess, S. Ostapchenko, R. Tomas

Abstract: 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… ▽ More 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. △ Less

Submitted 11 January, 2012; v1 submitted 27 June, 2011; originally announced June 2011.

Comments: 10 pages, 5 eps figures; program available at http://elmag.sourceforge.net; v2 matches published version

Ultrahigh Energy Nuclei in the Turbulent Galactic Magnetic Field

Authors: G. Giacinti, M. Kachelriess, D. V. Semikoz, G. Sigl

Abstract: 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… ▽ More 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. △ Less

Submitted 5 September, 2011; v1 submitted 6 April, 2011; originally announced April 2011.

Comments: 10 pages (2 columns), 8 figures. Published in Astroparticle Physics

Journal ref: Astropart. Phys. 35, 192-200 (2011)