Sources Of Ultra-high Energy Cosmic Rays Research Articles

How isotropic can the UHECR flux be?

Modern observatories of ultra-high energy cosmic rays (UHECR) have collected over 10^4 events with energies above 10 EeV, whose arrival directions appear to be nearly isotropically distributed. On the other hand, the distribution of matter in the nearby Universe -- and, therefore, presumably also that of UHECR sources -- is not homogeneous. This is expected to leave an imprint on the angular distribution of UHECR arrival directions, though deflections by cosmic magnetic fields can confound the picture. In this work, we investigate quantitatively this apparent inconsistency. To this end we study observables sensitive to UHECR source inhomogeneities but robust to uncertainties on magnetic fields and the UHECR mass composition. We show, in a rather model-independent way, that if the source distribution tracks the overall matter distribution, the arrival directions at energies above 30 EeV should exhibit a sizeable dipole and quadrupole anisotropy, detectable by UHECR observatories in the very near future. Were it not the case, one would have to seriously reconsider the present understanding of cosmic magnetic fields and/or the UHECR composition. Also, we show that the lack of a strong quadrupole moment above 10 EeV in the current data already disfavours a pure proton composition, and that in the very near future measurements of the dipole and quadrupole moment above 60 EeV will be able to provide evidence about the UHECR mass composition at those energies.

Erratum: Constraints on Ultrahigh-Energy Cosmic-Ray Sources from a Search for Neutrinos Above 10PeV with IceCube [Phys. Rev. Lett. 117, 241101 (2016)

This corrects the article DOI: 10.1103/PhysRevLett.117.241101.

Measurement of anisotropy and the search for ultra high energy cosmic ray sources

Ultra high energy cosmic rays (UHECRs) are particles, likely protons and/or nuclei, with energies up to |$10^{20}$| eV that are observed through the giant air showers they produce in the atmosphere. These particles carry information on the most extreme phenomena in the Universe. At these energies, even charged particles could be magnetically rigid enough to keep track of, or even point directly to, the original positions of their sources on the sky. The discovery of anisotropy of UHECRs would thus signify the opening of an entirely new window onto the Universe. With the construction and operation of the new generation of cosmic ray experiments—the Pierre Auger Observatory in the Southern hemisphere and the Telescope Array in the Northern one—the study of these particles, the most energetic ever detected, has experienced a jump in statistics as well as in data quality, allowing for much better sensitivity in searching for their sources. In this review, we summarize the searches for anisotropies and the efforts to identify the sources of UHECRs which have been carried out using these new data.

The Extragalactic Ultra-high-energy Cosmic-Ray Dipole

Abstract We explore the possibility that the recently detected dipole anisotropy in the arrival directions of >8 EeV ultra-high-energy cosmic-rays (UHECRs) arises due to the large-scale structure. We assume that the cosmic-ray sources follow the matter distribution and calculate the flux-weighted UHECRs’ rms dipole amplitude taking into account the diffusive transport in the intergalactic magnetic field (IGMF). We find that the flux-weighted rms dipole amplitude is ∼8% before entering the Galaxy. The amplitude in the [4–8] EeV is only slightly lower ∼5%. The required IGMF is of the order of 5–30 nG, and the UHECR sources must be relatively nearby, within ∼300 Mpc. The absence of a statistically significant signal in the lower-energy bin can be explained if the same nuclei specie dominates the composition in both energy bins and diffusion in the Galactic magnetic field reduces the dipole of these lower-rigidity particles. Photodisintegration of higher-energy UHECRs could also reduce somewhat the lower-energy dipole.

SimProp v2r4: Monte Carlo simulation code for UHECR propagation

We introduce the new version of SimProp, a Monte Carlo code for simulating the propagation of ultra-high energy cosmic rays in intergalactic space. This version, SimProp v2r4, together with an overall improvement of the code capabilities with a substantial reduction in the computation time, also computes secondary cosmogenic particles such as electron-positron pairs and gamma rays produced during the propagation of ultra-high energy cosmic rays. As recently pointed out by several authors, the flux of this secondary radiation and its products, within reach of the current observatories, provides useful information about models of ultra-high energy cosmic ray sources which would be hard to discriminate otherwise.

What have we learned about the sources of ultrahigh-energy cosmic rays via neutrino astronomy?

Observations of TeV–PeV-energy cosmic neutrinos by the IceCube observatory have suggested that extragalactic cosmic-ray sources should have an optical depth greater than ∼0.01 and contribute to more than 10% of the observed bulk of cosmic rays at 10 PeV. If the spectrum of cosmic rays from these extragalactic sources extends well beyond 1 EeV, the neutrino flux indicates that extragalactic cosmic-ray protons are dominant in the observed total cosmic-ray flux at 1 EeV. Among known powerful astronomical objects, including gamma-ray bursters (GRBs), only flat-spectrum radio quasars could (barely) satisfy these conditions. On the other hand, the null detection of neutrinos with energies well beyond PeV has excluded the possibility that radio-loud active galactic nuclei (AGNs) and/or GRBs, the popular source candidates discussed in the literature, are the origins of the highest-energy cosmic rays (∼ 100EeV) if they are composed mainly of protons. Their origins must be objects that have evolved on time scales comparable to or slower than the star formation rate. These considerations indicate that none of the known extragalactic astronomical objects can be simultaneously a source of both PeV- and trans-EeV-energy cosmic rays. As a result of the stringent limits on EeV-energy neutrino fluxes, a significant part of the parameter space for the AGN and new-born pulsar models is starting to seem unfavorable, even for scenarios of mixed and heavy cosmic-ray compositions at the highest energies.

High-energy cosmic ray nuclei from tidal disruption events: Origin, survival, and implications

Tidal disruption events (TDEs) by supermassive or intermediate mass black holes have been suggested as candidate sources of ultrahigh-energy cosmic rays (UHECRs) and high-energy neutrinos. Motivated by the recent measurements from the Pierre Auger Observatory, which indicates a metal-rich cosmic-ray composition at ultrahigh energies, we investigate the fate of UHECR nuclei loaded in TDE jets. First, we consider the production and survival of UHECR nuclei at internal shocks, external forward and reverse shocks, and nonrelativistic winds. Based on the observations of Swift J1644+57, we show that the UHECRs can survive for external reverse and forward shocks, and disk winds. On the other hand, UHECR nuclei are significantly disintegrated in internal shocks, although they could survive for low-luminosity TDE jets. Assuming that UHECR nuclei can survive, we consider implications of different composition models of TDEs. We find that the tidal disruption of main sequence stars or carbon-oxygen white dwarfs does not successfully reproduce UHECR observations, namely the observed composition or spectrum. The observed mean depth of the shower maximum and its deviation could be explained by oxygen-neon-magnesium white dwarfs, but they may be too rare to be the sources of UHECRs.

Probing the Extragalactic Cosmic-Ray Origin with Gamma-Ray and Neutrino Backgrounds

Abstract GeV–TeV gamma-rays and PeV–EeV neutrino backgrounds provide a unique window on the nature of the ultra-high-energy cosmic rays (UHECRs). We discuss the implications of the recent Fermi-LAT data regarding the extragalactic gamma-ray background and related estimates of the contribution of point sources as well as IceCube neutrino data on the origin of the UHECRs. We calculate the diffuse flux of cosmogenic γ-rays and neutrinos produced by the UHECRs and derive constraints on the possible cosmological evolution of UHECR sources. In particular, we show that the mixed-composition scenario considered in Globus et al., which is in agreement with both (i) Auger measurements of the energy spectrum and composition up to the highest energies and (ii) the ankle-like feature in the light component detected by KASCADE-Grande, is compatible with both the Fermi-LAT measurements and with current IceCube limits. We also discuss the possibility for future experiments to detect associated cosmogenic neutrinos and further constrain the UHECR models, including possible subdominant UHECR proton sources.

Cosmogenic photons strongly constrain UHECR source models

With the newest version of our Monte Carlo code for ultra-high-energy cosmic ray (UHECR) propagation, CRPropa 3, the flux of neutrinos and photons due to interactions of UHECRs with extragalactic background light can be predicted. Together with the recently updated data for the isotropic diffuse gamma-ray background (IGRB) by Fermi LAT, it is now possible to severely constrain UHECR source models. The evolution of the UHECR sources especially plays an important role in the determination of the expected secondary photon spectrum. Pure proton UHECR models are already strongly constrained, primarily by the highest energy bins of Fermi LAT's IGRB, as long as their number density is not strongly peaked at recent times.

Constraints on Ultrahigh-Energy Cosmic-Ray Sources from a Search for Neutrinos above 10PeV with IceCube.

We report constraints on the sources of ultrahigh-energy cosmic rays (UHECRs) above 10^{9} GeV, based on an analysis of seven years of IceCube data. This analysis efficiently selects very high- energy neutrino-induced events which have deposited energies from 5×10^{5} GeV to above 10^{11} GeV. Two neutrino-induced events with an estimated deposited energy of (2.6±0.3)×10^{6} GeV, the highest neutrino energy observed so far, and (7.7±2.0)×10^{5} GeV were detected. The atmospheric background-only hypothesis of detecting these events is rejected at 3.6σ. The hypothesis that the observed events are of cosmogenic origin is also rejected at >99% CL because of the limited deposited energy and the nonobservation of events at higher energy, while their observation is consistent with an astrophysical origin. Our limits on cosmogenic neutrino fluxes disfavor the UHECR sources having a cosmological evolution stronger than the star formation rate, e.g., active galactic nuclei and γ-ray bursts, assuming proton-dominated UHECRs. Constraints on UHECR sources including mixed and heavy UHECR compositions are obtained for models of neutrino production within UHECR sources. Our limit disfavors a significant part of parameter space for active galactic nuclei and new-born pulsar models. These limits on the ultrahigh-energy neutrino flux models are the most stringent to date.

Searches for correlation between UHECR events and high-energy gamma-ray Fermi-LAT data

The astrophysical sources responsible for ultra high-energy cosmic rays (UHECRs) continue to be one of the most intriguing mysteries in astrophysics. We present a comprehensive search for correlations between high-energy (≳ 1 GeV) gamma-ray events from the Fermi Large Area Telescope (LAT) and UHECRs (≳ 60 EeV) detected by the Telescope Array and the Pierre Auger Observatory. We perform two separate searches. First, we conduct a standard cross-correlation analysis between the arrival directions of 148 UHECRs and 360 gamma-ray sources in the Second Catalog of Hard Fermi-LAT sources (2FHL). Second, we search for a possible correlation between UHECR directions and unresolved Fermi-LAT gamma-ray emission. For the latter, we use three different methods: a stacking technique with both a model-dependent and model-independent background estimate, and a cross-correlation function analysis. We also test for statistically significant excesses in gamma rays from signal regions centered on Cen A and the Telescope Array hotspot. No significant correlation is found in any of the analyses performed, except a weak (≲ 2σ) hint of signal with the correlation function method on scales ∼ 1°. Upper limits on the flux of possible power-law gamma-ray sources of UHECRs are derived.

Infinite efficiency of the collisional Penrose process: Can a overspinning Kerr geometry be the source of ultrahigh-energy cosmic rays and neutrinos?

The origin of the ultra-high-energy particles we receive on the Earth from the outer space such as EeV cosmic rays and PeV neutrinos remains an enigma. All mechanisms known to us currently make use of electromagnetic interaction to accelerate charged particles. In this paper we propose a mechanism exclusively based on gravity rather than electromagnetic interaction. We show that it is possible to generate ultra-high-energy particles starting from particles with moderate energies using the collisional Penrose process in an overspinning Kerr spacetime transcending the Kerr bound only by an infinitesimal amount, i.e., with the Kerr parameter $a=M(1+\epsilon)$, where we take the limit $\epsilon \rightarrow 0^+$. We consider two massive particles starting from rest at infinity that collide at $r=M$ with divergent center-of-mass energy and produce two massless particles. We show that massless particles produced in the collision can escape to infinity with the ultra-high energies exploiting the collisional Penrose process with the divergent efficiency $\eta \sim {1}/{\sqrt{\epsilon}} \rightarrow \infty$. Assuming the isotropic emission of massless particles in the center-of-mass frame of the colliding particles, we show that half of the particles created in the collisions escape to infinity with the divergent energies. To a distant observer, ultra-high-energy particles appear to originate from a bright spot which is at the angular location $\xi \sim {2M}/{r_{obs}}$ with respect to the singularity on the side which is rotating towards the observer. We show that the anisotropy in emission in the center-of-mass frame, which is dictated by the differential cross-section of underlying particle physics process, leaves a district signature on the spectrum of ultra-high-energy massless particles. Thus, it provides a unique probe into fundamental particle physics.

ON ULTRA-HIGH-ENERGY COSMIC RAYS AND THEIR RESULTANT GAMMA-RAYS

ABSTRACT The Fermi Large Area Telescope collaboration has recently reported on 50 months of measurements of the isotropic extragalactic gamma-ray background (EGRB) spectrum between 100 MeV and 820 GeV. Ultra-high-energy cosmic ray (UHECR) protons interact with the cosmic microwave background photons and produce cascade photons of energies 10 MeV–1 TeV that contribute to the EGRB flux. We examine seven possible evolution models for UHECRs and find that UHECR sources that evolve as the star formation rate (SFR), medium low luminosity active galactic nuclei type-1 (L = 1043.5 erg s−1 in the [0.5–2] KeV band), and BL Lacertae objects (BL Lacs) are the most acceptable given the constraints imposed by the observed EGRB. Other possibilities produce too much secondary γ-radiation. In all cases, the decaying dark matter (DM) contribution improves the fit at high energy, but the contribution of still unresolved blazars, which would leave the smallest role for decaying DM, may yet provide an alternative improvement. The possibility that the entire EGRB can be fitted with resolvable but not-yet-resolved blazars, as recently claimed by Ajello et al., would leave little room in the EGRB to accommodate γ-rays from extragalactic UHECR production, even for many source evolution rates that would otherwise be acceptable. We find that under the assumption of UHECRs being mostly protons, there is not enough room for producing extragalactic UHECRs with active galactic nucleus, gamma-ray burst, or even SFR source evolution. Sources that evolve as BL Lacs, on the other hand, would produce much less secondary γ-radiation and would remain a viable source of UHECRs, provided that they dominate.

ISOTROPY CONSTRAINTS ON POWERFUL SOURCES OF ULTRAHIGH-ENERGY COSMIC RAYS AT 1019 eV

ABSTRACT Anisotropy in the arrival direction distribution of ultrahigh-energy cosmic rays (UHECRs) produced by powerful sources is numerically evaluated. We show that nondetection of significant anisotropy at ≈ 10 19 ?> eV at present and in future experiments imposes general upper limits on UHECR proton luminosity of steady sources as a function of source redshifts. The upper limits constrain the existence of typical steady 10 19 ?> eV UHECR sources in the local universe and limit their local density to ≳ 10 − 3 ?> Mpc − 3 , ?> assuming average intergalactic magnetic fields less than 10 − 9 ?> G. This isotropy, being stronger than that measured at the highest energies, may indicate the transient generation of UHECRs. Our calculations are applied for extreme high-frequency-peaked BL Lacertae objects 1ES 0229+200, 1ES 1101–232, and 1ES 0347–121, to test the UHECR-induced cascade model, in which beamed UHECR protons generate TeV radiation in transit from sources. While the magnetic-field structure surrounding the sources affects the required absolute cosmic-ray luminosity of the blazars, the magnetic-field structure surrounding the Milky Way directly affects the observed anisotropy. If these magnetic fields are weak enough, significant UHECR anisotropy from these blazars is detectable by the Pierre Auger Observatory unless the maximum energy of UHECR protons is below 1019 eV. Furthermore, if these are the sources of UHECRs above 1019 eV, a local magnetic structure surrounding the Milky Way is needed to explain the observed isotropy at ∼ 10 19 ?> eV, which may be incompatible with large magnetic structures around all galaxies for the UHECR-induced cascade model to work with reasonable jet powers.

Multi-messenger aspects of cosmic neutrinos*

The recent observation of TeV-PeV neutrinos by IceCube has opened a new window to the high-energy Universe. I will discuss this signal in the context of multi-messenger astronomy. For extragalactic source scenarios the corresponding gamma-rays are not directly observable due to interactions with the cosmic radiation backgrounds. Nevertheless, the isotropic sub-TeV gamma ray background observed by Fermi -LAT contains indirect information from secondary emission produced in electromagnetic cascades. On the other hand, observation of PeV gamma rays would provide a smoking-gun signal for Galactic emission. Interestingly, the overall energy density of the observed neutrino flux is close to a theoretical limit for neutrino production in ultra-high energy cosmic ray sources and might indicate a common origin of these phenomena. I will highlight various multi-messenger relations and their implications for neutrino source scenarios.

Relativistic Shocks: Particle Acceleration and Magnetization

We review the physics of relativistic shocks, which are often invoked as the sources of non-thermal particles in pulsar wind nebulae (PWNe), gamma-ray bursts (GRBs), and active galactic nuclei (AGN) jets, and as possible sources of ultra-high energy cosmic-rays. We focus on particle acceleration and magnetic field generation, and describe the recent progress in the field driven by theory advances and by the rapid development of particle-in-cell (PIC) simulations. In weakly magnetized or quasi parallel-shocks (where the magnetic field is nearly aligned with the flow), particle acceleration is efficient. The accelerated particles stream ahead of the shock, where they generate strong magnetic waves which in turn scatter the particles back and forth across the shock, mediating their acceleration. In contrast, in strongly magnetized quasi-perpendicular shocks, the efficiencies of both particle acceleration and magnetic field generation are suppressed. Particle acceleration, when efficient, modifies the turbulence around the shock on a long time scale, and the accelerated particles have a characteristic energy spectral index of ~ 2.2 in the ultra-relativistic limit. We discuss how this novel understanding of particle acceleration and magnetic field generation in relativistic shocks can be applied to high-energy astrophysical phenomena, with an emphasis on PWNe and GRB afterglows.

NEAR-EQUIPARTITION JETS WITH LOG-PARABOLA ELECTRON ENERGY DISTRIBUTION AND THE BLAZAR SPECTRAL-INDEX DIAGRAMS

Fermi-LAT analyses show that the gamma-ray photon spectral indices Gamma_gamma of a large sample of blazars correlate with the vFv peak synchrotron frequency v_s according to the relation Gamma_gamma = d - k log v_s. The same function, with different constants d and k, also describes the relationship between Gamma_gamma and peak Compton frequency v_C. This behavior is derived analytically using an equipartition blazar model with a log-parabola description of the electron energy distribution (EED). In the Thomson regime, k = k_EC = 3b/4 for external Compton processes and k = k_SSC = 9b/16 for synchrotron self-Compton (SSC) processes, where b is the log-parabola width parameter of the EED. The BL Lac object Mrk 501 is fit with a synchrotron/SSC model given by the log-parabola EED, and is best fit away from equipartition. Corrections are made to the spectral-index diagrams for a low-energy power-law EED and departures from equipartition, as constrained by absolute jet power. Analytic expressions are compared with numerical values derived from self-Compton and external Compton scattered gamma-ray spectra from Ly alpha broad-line region and IR target photons. The Gamma_gamma vs. v_s behavior in the model depends strongly on b, with progressively and predictably weaker dependences on gamma-ray detection range, variability time, and isotropic gamma-ray luminosity. Implications for blazar unification and blazars as ultra-high energy cosmic-ray sources are presented. Arguments by Ghisellini et al. (2014) that the jet power exceeds the accretion luminosity depend on the doubtful assumption that we are viewing at the Doppler angle.

Simulations for a next-generation UHECR observatory

We explore the potential of a future, ultra-high energy cosmic ray (UHECR) experiment, that is able to overcome the limitation of low statistics, to detect anisotropy in the arrival directions of UHECRs. We concentrate on the lower energy range of future instruments (E≳ 50 EeV), where, if the UHECR source number density is not too low, the sources should be numerous enough to imprint a clustering pattern in the sky, and thus possibly in the UHECR arrival directions. Under these limits, the anisotropy signal should be dominated by the clustering of astrophysical sources per se in the large-scale structures, and not the clustering of events around individual sources. We study the potential for a statistical discrimination between different astrophysical models which we parametrise by the number density of UHECR sources, the possible bias of the UHECR accelerators with respect to the galaxy distribution, and the unknown fraction of UHECRs that have been deflected by large angles. We demonstrate that an order-of-magnitude increase in statistics would allow to discriminate between a variety of astrophysical models, provided that a sub-sample of light elements can be extracted, and that it represents a fraction ≳ 70% of the overall flux, sensitive to the UHECR source number density. Discrimination should be possible even without knowledge of the composition of the UHECRs, as long as the data are proton-dominated. We find that an anisotropy at the 99.7% level should be detectable when the number of detected events exceeds 2000 beyond 50 EeV, as long as the composition is proton dominated, and the number density of UHECR sources is relatively high, n̄ ⩾ 10−3 Mpc−3. If the UHECR sources are strongly biased relative to the galaxy distribution, as are for example galaxy clusters, in which the sources might be embedded, an anisotropy at the 99.7% level should be detectable once the number of detected events exceeds 1000, if the fraction of protons at the highest energies is ≳ 60%.

Testing the newborn pulsar origin of ultrahigh energy cosmic rays with EeV neutrinos

Fast-spinning newborn pulsars are intriguing candidate sources of ultrahigh-energy cosmic rays (UHECRs). The acceleration of particles with a given composition in a fraction of the extragalactic pulsar population can give a consistent explanation for the measurements of the Auger Observatory. We calculate the associated diffuse neutrino flux produced while particles cross the supernova ejecta surrounding the stars. We show that in the minimal pulsar scenarios that are compatible with the UHECR data, the effective optical depth to hadronuclear interactions is larger than unity at ultrahigh energies. Thus, even in the most pessimistic case, one expects energy fluxes of $\ensuremath{\sim}0.1--1\text{ }\text{ }\mathrm{EeV}$ neutrinos that should be detectable with IceCube or the Askaryan Radio Array within a decade.

The JEM-EUSO instrument

In this paper we describe the main characteristics of the JEM-EUSO instrument. The Extreme Universe Space Observatory on the Japanese Experiment Module (JEM-EUSO) of the International Space Station (ISS) will observe Ultra High-Energy Cosmic Rays (UHECR) from space. It will detect UV-light of Extensive Air Showers (EAS) produced by UHECRs traversing the Earth’s atmosphere. For each event, the detector will determine the energy, arrival direction and the type of the primary particle. The advantage of a space-borne detector resides in the large field of view, using a target volume of about 1012 tons of atmosphere, far greater than what is achievable from ground. Another advantage is a nearly uniform sampling of the whole celestial sphere. The corresponding increase in statistics will help to clarify the origin and sources of UHECRs and characterize the environment traversed during their production and propagation. JEM-EUSO is a 1.1 ton refractor telescope using an optics of 2.5 m diameter Fresnel lenses to focus the UV-light from EAS on a focal surface composed of about 5,000 multi-anode photomultipliers, for a total of ≃3⋅105 channels. A multi-layer parallel architecture handles front-end acquisition, selecting and storing valid triggers. Each processing level filters the events with increasingly complex algorithms using FPGAs and DSPs to reject spurious events and reduce the data rate to a value compatible with downlink constraints.