The Origin of Highest Energy Cosmic Rays and Cosmic Magnetic Fields

We summarize the status of the search for the origin of the highest energy cosmic rays. We briefly mention several competing proposals, such as supersymmetric particles, Gamma Ray Bursts also giving rise to energetic protons, interacting high energy neutrinos and cosmological defects, and then concentrate on the possibilities of the propagation of these particles, assuming that they are charged. The distribution of arrival directions of the highest energy particles on the sky ought to reflect the source distribution as well as propagation history. The present status can be summarized as inconclusive. If we were able to confirm any particular theory these particles beyond 10²⁰eV may be turned into tools of high energy physics. This revised version was published online in July 2006 with corrections to the Cover Date.     

Galactic PeV neutrinos

The IceCube experiment has detected two neutrinos with energies between 1 and 10 PeV. They might have originated from Galactic or extragalactic sources of cosmic rays. In the present work we consider hadronic interactions of the diffuse very high energy cosmic rays with the interstellar matter within our Galaxy to explain the PeV neutrino events detected in IceCube. We also expect PeV gamma ray events along with the PeV neutrino events if the observed PeV neutrinos were produced within our Galaxy in hadronic interactions. PeV gamma rays are unlikely to reach us from sources outside our Galaxy due to pair production with cosmic background radiation fields. We suggest that in future with simultaneous detections of PeV gamma rays and neutrinos it would be possible to distinguish between Galactic and extragalactic origins of very high energy neutrinos.     

Dark matter electron anisotropy: A universal upper limit

We study the dipole anisotropy in the arrival directions of high energy CR electrons and positrons (CRE) of dark matter (DM) origin. We show that this quantity is very weakly model dependent and offers a viable criterion to discriminate among CRE from DM or from local discrete sources, like e.g. pulsars. In particular, we find that the maximum anisotropy which DM can provide is to a very good approximation a universal quantity and, as a consequence, if a larger anisotropy is detected, this would constitute a strong evidence for the presence of astrophysical local discrete CRE sources, whose anisotropy, instead, can be naturally larger than the DM upper limit. We further find that the main source of anisotropy from DM is given by the fluctuation in the number density of DM sub-structures in the vicinity of the observer and we thus devote special attention to the study of the variance in the sub-structures realization implementing a dedicated Montecarlo simulation. Such scenarios will be probed in the next years by Fermi-LAT, providing new hints, or constraints, about the nature of DM.     

Potential neutrino signals in a northern hemisphere neutrino telescope from galactic gamma-ray sources

Neutrino energy spectra have been calculated based on the recently measured energy spectra of Galactic very high energy γ-ray sources. Based on these neutrino spectra the expected event rates in the ANTARES neutrino telescope and KM3NeT, a future neutrino telescope in the Mediterranean Sea with an instrumented volume of one km³, have been calculated. For the brightest γ-ray sources we find event rates of the order of one neutrino per year. Although the neutrino event rates are comparable to the background from atmospheric neutrinos the detection of individual sources seems possible.     

The ANTARES telescope neutrino alert system

The ANTARES telescope has the capability to detect neutrinos produced in astrophysical transient sources. Potential sources include gamma-ray bursts, core collapse supernovae, and flaring active galactic nuclei. To enhance the sensitivity of ANTARES to such sources, a new detection method based on coincident observations of neutrinos and optical signals has been developed. A fast online muon track reconstruction is used to trigger a network of small automatic optical telescopes. Such alerts are generated for special events, such as two or more neutrinos, coincident in time and direction, or single neutrinos of very high energy.     

Prospects for radio detection of ultra-high energy cosmic rays and neutrinos

The origin and nature of the highest energy cosmic ray events is currently the subject of intense investigation by giant air shower arrays and fluorescent detectors. These particles reach energies well beyond what can be achieved in ground-based particle accelerators and hence they are fundamental probes for particle physics as well as astrophysics. One of the main topics today focuses on the high energy end of the spectrum and the potential for the production of high-energy neutrinos. Above about 10²⁰ eV cosmic rays from extragalactic sources are expected to be severely attenuated by pion photoproduction interactions with photons of the cosmic microwave background. Investigating the shape of the cosmic ray spectrum near this predicted cut-off will be very important. In addition, a significant high-energy neutrino background is naturally expected as part of the pion decay chain which also contains much information.Because of the scarcity of these high-energy particles, larger and larger ground-based detectors have been built. The new generation of digital radio telescopes may play an important role in this, if properly designed. Radio detection of cosmic ray showers has a long history but was abandoned in the 1970s. Recent experimental developments together with sophisticated air shower simulations incorporating radio emission give a clearer understanding of the relationship between the air shower parameters and the radio signal, and have led to resurgence in its use. Observations of air showers by the SKA could, because of its large collecting area, contribute significantly to measuring the cosmic ray spectrum at the highest energies. Because of the large surface area of the moon, and the expected excellent angular resolution of the SKA, using the SKA to detect radio Cherenkov emission from neutrino-induced cascades in lunar regolith will be potentially the most important technique for investigating cosmic ray origin at energies above the photoproduction cut-off.     

Could a plasma in quasi-thermal equilibrium be associated to the “orphan” TeV flares?

TeV γ-ray detections in flaring states without activity in X-rays from blazars have attracted much attention due to the irregularity of these “orphan” flares. Although the synchrotron self-Compton model has been very successful in explaining the spectral energy distribution and spectral variability of these sources, it has not been able to describe these atypical flaring events. On the other hand, an electron–positron pair plasma at the base of the AGN jet was proposed as the mechanism of bulk acceleration of relativistic outflows. This plasma in quasi-thermal equilibrium called Wein fireball emits radiation at MeV-peak energies serving as target of accelerated protons. In this work we describe the “orphan” TeV flares presented in blazars 1ES 1959+650 and Mrk 421 assuming geometrical considerations in the jet and evoking the interactions of Fermi-accelerated protons and MeV-peak target photons coming from the Wein fireball. After describing successfully these “orphan” TeV flares, we correlate the TeV γ-ray, neutrino and UHECR fluxes through pγ interactions and calculate the number of high-energy neutrinos and UHECRs expected in IceCube/AMANDA and TA experiment, respectively. In addition, thermal MeV neutrinos produced mainly through electron–positron annihilation at the Wein fireball will be able to propagate through it. By considering two- (solar, atmospheric and accelerator parameters) and three-neutrino mixing, we study the resonant oscillations and estimate the neutrino flavor ratios as well as the number of thermal neutrinos expected on Earth.     

Cherenkov τ shower earth-skimming method for PeV–EeV ντ observation with Ashra

We describe a method of observation for PeV–EeV τ neutrinos using Cherenkov light from the air showers of decayed τs produced by τ neutrino interactions in the Earth. Aiming for the realization of neutrino astronomy utilizing the Earth-skimming τ neutrino detection technique, highly precise determination of arrival direction is key due to the following issues: (1) clear identification of neutrinos by identifying those vertices originating within the Earth’s surface and (2) identification of very high energy neutrino sources. The Ashra detector uses newly developed light collectors which realize both a 42°-diameter field-of-view and arcminute resolution. Therefore, it has superior angular resolution for imaging Cherenkov air showers. In this paper, we estimate the sensitivity of and cosmic-ray background resulting from application of the Ashra-1 Cherenkov τ shower observation method. Both data from a commissioning run and a long-term observation (with fully equipped trigger system and one light collector) are presented. Our estimates are based on a detailed Monte Carlo simulation which describes all relevant shower processes from neutrino interaction to Cherenkov photon detection produced by τ air showers. In addition, the potential to determine the arrival direction of Cherenkov showers is evaluated by using the maximum likelihood method. We conclude that the Ashra-1 detector is a unique probe into detection of very high energy neutrinos and their accelerators.     

Cascades in the Earth’s magnetosphere initiated by photons with the parameters of the highest energy AGASA events

We investigate the cascading effects of extremely high energy (EHE) photons in the Earth’s magnetosphere assuming that these photons arrive with the parameters of the highest energy AGASA events (energies, arrival directions). For the location of the AGASA Observatory, we determine the directions in the sky from which photons can cascade with a high (low) probability. In the case of the primary photons with the parameters of the events with the energies >10²⁰ eV, we compute the average cascade spectra of secondary photons entering the Earth’s atmosphere, and estimate their fluctuations around these average values by selecting the events with the largest and smallest number of secondary cascade photons. It is shown that most photons with the parameters of the highest energy AGASA events should initiate cascades in the Earth’s magnetosphere with a high probability even though they tend to arrive from directions in the sky for which the perpendicular component of the magnetic field is weaker. On the other hand, if these events are caused by the photons with lower energies, then the fluctuations in their shower development in the magnetosphere and the atmosphere should be higher than in the case of photons with the energies estimated by the AGASA experiment.     

New observing modes of a future HTN‐based distributed detection network for high time resolution and quantum optical astrophysics experiments

This contribution aims to introduce the idea that a well‐evolved HTN of the far future, with the anticipated addition of very large apertures, could also be made to incorporate the ability to carry out photonic astronomy observations, particularly Optical VLBI in a revived Hanbury‐Brown and Twiss Intensity Interferometry (HBTII) configuration. Such an HTN could exploit its inherent rapid reconfigurational ability to become a multi‐aperture distributed photon‐counting network able to study higher‐order spatiotemporal photon correlations and provide a unique tool for direct diagnostics of astrophysical emission processes. We very briefly review various considerations associated with the switching of the HTN to a special mode in which single‐photon detection events are continuously captured for a posteriori intercorrelation. In this context, photon arrival times should be determined to the highest time resolution possible and extremely demanding absolute time keeping and absolute time distribution schemes should be devised and implemented in the HTN nodes involved. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Upper limits on neutrino fluxes from point-like sources with AMANDA-II

The AMANDA-II telescope, operated by the IceCube collaboration, is currently the world’s most sensitive telescope to fluxes of neutrinos from individual sources. A data sample of 4282 neutrino induced events collected in 1001 days of detector livetime during the years 2000–2004 have now been analyzed looking for a neutrino signal from point-like sources. A sensitivity to fluxes of of d Φ/dE=1.0×10⁻¹⁰(E/TeV)⁻² TeV⁻¹ cm⁻²s⁻¹ was reached in the energy range between 1.7 TeV and 2.4 PeV. So far no statistically significant localized excess of events over the background of atmospheric neutrinos has been found, which would be ascribed to a neutrino source. However, the flux upper limits derived from the non-observation of a signal are comparable to observed fluxes of high energy gamma rays from blazars and within the range of current models for neutrino emission from selected sources. Possible constraints on these models are discussed.      

LS I +61 303 as a potential neutrino source on the light of magic results

Very high energy γ-rays have recently been detected from the microquasar LS I +61 303 using the MAGIC telescope. A phenomenological study on the concomitant neutrinos that would be radiated if the γ-ray emission is hadronic in origin is herein presented. Neutrino oscillations are considered, and the expected number of events in a km-scale detector such as ICECUBE is computed under different assumptions including orbital periodicity and modulation, as well as different precision in the modeling of the detector. We argue that the upper limits already imposed on the neutrino emission of LS I +61 303 using AMANDA-II and the forthcoming measurements by ICECUBE may significantly constrain – in an independent and unbiased way – the γ-ray to neutrino flux ratio, and thus the possibility of a hadronic origin of the γ-rays. The viability of hadronic models based on wind–jet interactions in the LS +61 303 system after MAGIC measurements is discussed.     

Cosmic ray anisotropy analysis with a full-sky observatory

A cosmic ray observatory with full-sky coverage can exploit standard anisotropy analysis methods that do not work if part of the celestial sphere is never seen. In particular, the distribution of arrival directions can be fully characterized by a list of spherical harmonic coefficients. The dipole vector and quadrupole tensor are of special interest, but the full set of harmonic coefficients constitutes the anisotropy fingerprint that may be needed to reveal the identity of the cosmic ray sources. The angular power spectrum is a coordinate-independent synopsis of that fingerprint. The true cosmic ray anisotropy can be measured despite non-uniformity in celestial exposure, provided the observatory is not blind to any region of the sky. This paper quantitatively examines how the accuracy of anisotropy measurement depends on the number of arrival directions in a data set.     

Source structure and dynamics of Pc 1 pulsations at low latitudes

Structured Pc 1 signals propagate in the ionospheric F2 region duct from their secondary sources at high latitudes to lower latitudes. Propagation directions to low latitude stations can be inferred from measurements of polarization parameters. The analysis of five events recorded at two low latitude stations (L = 1.9) are presented. Direction of arrival measurements are used to investigate the spatial and temporal structure of Pc 1 sources. Results show a close relationship between the structure of events identified in the frequency-time representation and direction of arrival measurement patterns. Multiple sources are sometimes indicated.     

3XMM J185246.6+003317 as transient neutron star

We performed an X‐ray timing and spectral analysis of the variable source 3XMM J185246.6+003317 to investigate its physical nature. The data from all observations of 3XMMJ185246.6+003317 conducted by XMM‐Newton EPIC MOS1 and MOS2 with the same instrumental setup in 2004–2009 were reprocessed to form a homogenous data set of solar barycenter corrected photon arrival times and high S/N spectra of 3XMM J185246.6+003317. A Bayesian method for the search, detection, and estimation of the parameters of a periodic signal of unknown shape was employed, as developed by Gregory & Loredo (1992, 1993). The results show that 3XMM J185246.6+003317 is a transient neutron star with the genuine spin‐period of 23.11722 (23.11711–23.11727) s and its derivative of 5.3(0.3–5.5)×10^–11 s s^–1, implying a characteristic age of 7 (6–104) kyr, if the period derivative can be ascribed to the genuine spin‐down rate of the neutron star. The rotational‐phase averaged X‐ray spectra at the different brightness periods can be fitted with a highly absorbed blackbody model with different temperatures. The phase‐folded light curves in different energy bands with high S/N ratio show a double‐peaked profile; the variations depend on time and energy, indicating that radiation emerges from at least two emitting areas. The spectra at the phases corresponding to the maxima in the phase‐folded light curve show different spectral parameters of absorbed blackbody radiation, i.e. the hotter one has a smaller size. The source is detected only from September 2008 to April 2009 with persistently decreasing brightness, but not before, even though it was observed by XMM. Hence, it is a transient neutron star or a binary system hosting it. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Detecting neutrino transients with optical follow-up observations

A novel method is presented which will enhance the sensitivity of neutrino telescopes to identify transient sources such as Gamma-Ray Bursts (GRBs) and core-collapse Supernovae (SNe). Triggered by the detection of high energy neutrino events from IceCube or other large scale neutrino telescopes, an optical follow-up program will allow the identification of the transient neutrino source. We show that once the follow-up program is implemented, the achievable sensitivity of IceCube to neutrinos from SNe and GRBs would increase by a factor of 2–3. The program can be realized with a small network of automated 1–2 m telescopes and has rather modest observing time requirements.     

GRBs on probation: Testing the UHE CR paradigm with IceCube

Gamma ray burst (GRB) fireballs provide one of very few astrophysical environments where one can contemplate the acceleration of cosmic rays to energies that exceed 10²⁰ eV. The assumption that GRBs are the sources of the observed cosmic rays generates a calculable flux of neutrinos produced when the protons interact with fireball photons. With data taken during construction IceCube has already reached a sensitivity to observe neutrinos produced in temporal coincidence with individual GRBs provided that they are the sources of the observed extra-galactic cosmic rays. We here point out that the GRB origin of cosmic rays is also challenged by the IceCube upper limit on a possible diffuse flux of cosmic neutrinos which should not be exceeded by the flux produced by all GRB over Hubble time. Our alternative approach has the advantage of directly relating the diffuse flux produced by all GRBs to measurements of the cosmic ray flux. It also generates both the neutrino flux produced by the sources and the associated cosmogenic neutrino flux in a synergetic way.     

A possible universal origin of hadronic cosmic rays from ultrarelativistic ejecta of bipolar supernovae

Based on the “cannonball model” for gamma-ray bursts of Dar and De Rújula, it is proposed that masses of baryonic plasma (“cannonballs”), ejected in bipolar supernova explosions in our Galaxy are the sources of hadronic Galactic cosmic rays (CRs) at all energies. The propagation of the cannonballs in the Galactic disk and halo is studied. Two mechanisms for the acceleration of the observed CRs are proposed. The first is based on ultrarelativistic shocks in the interstellar medium and could accelerate the bulk of CRs up to the “knee” energy of 4×10¹⁵ eV. The second operates with second-order Fermi acceleration within the cannonball. If the total initial energy of the ejected plasmoids in a SN explosion is 10⁵³ erg or higher, this second mechanism may explain the CR spectrum above the knee up to the highest observed energies. It is shown that together with plausible assumptions about CR propagation in the Galactic confinement volume, the observed spectral indices of the CR spectrum can be theoretically understood to first order. The model allows a natural understanding of various basic CR observations like the absence of the Greisen–Zatsepin cutoff, the anisotropy of arrival directions as function of energy and the small Galactocentric gradient of the CR density.     

TeV γ-ray fluxes from the long campaigns on Mrk 421 as constraints on the emission of TeV–PeV neutrinos and UHECRs

Long TeV γ-ray campaigns have been carried out to study the spectrum, variability and duty cycle of the BL Lac object Markarian 421. These campaigns have given some evidence of the presence of protons in the jet: (i) Its spectral energy distribution which shows two main peaks; one at low energies (∼1 keV) and the other at high energies (hundreds of GeV), has been described by using synchrotron proton blazar model. (ii) The study of the variability at GeV γ-rays and X-rays has indicated no significant correlation. (iii) TeV γ-ray detections without activity in X-rays, called “orphan flares” have been observed in this object.Recently, The Telescope Array Collaboration reported the arrival of 72 ultra-high-energy cosmic rays with some of them possibly related to the direction of Markarian 421. The IceCube Collaboration reported the detection of 37 extraterrestrial neutrinos in the TeV–PeV energy range collected during three consecutive years. In particular, no neutrino track events were associated with this source. In this paper, we consider the proton–photon interactions to correlate the TeV γ-ray fluxes reported by long campaigns with the neutrino and ultra-high-energy cosmic ray observations around this blazar. Considering the results reported by The IceCube and Telescope Array Collaborations, we found that only from ∼25% to 70% of TeV fluxes described with a power law function with exponential cutoff can come from the proton–photon interactions.     

Search for point-like sources of cosmic rays with energies above 1018.5 eV in the HiRes-I monocular data set

We report the results of a search for point-like deviations from isotropy in the arrival directions of ultra-high energy cosmic rays in the northern hemisphere. In the monocular data set collected by the High-Resolution Fly’s Eye, consisting of 1525 events with energy exceeding 10^18.5 eV, we find no evidence for point-like excesses. We place a 90% c.l. upper limit of 0.8 hadronic cosmic rays/km² yr on the flux from such sources for the northern hemisphere and place tighter limits as a function of position in the sky.