High-energy neutrinos are expected to be produced by the interaction of accelerated particles near the acceleration sites. For this reason, it is intresting to search for correlation in the arrival ...directions of ultra–high energy cosmic rays (UHECRs) and HE neutrinos. We present here the results of a search for correlations between UHECR events measured by the Pierre Auger Observatory and Telescope Array and high-energy neutrino candidate events from IceCube and ANTARES. We perform a cross-correlation analysis, where the angular separation between the arrival directions of UHECRs and neutrinos is scanned. When comparing the results with the expectations from a null hypothesis contemplating an isotropic distribution of neutrinos or of UHECR we obtain post-trial p-values of the order of ~ 10
2
.
Context.
In recent years, evidence for an anisotropic distribution of ultra-high-energy cosmic rays (UHECRs) has been claimed, notably a dipole modulation in right ascension has been reported by the ...Auger collaboration above the 5
σ
significance threshold.
Aims.
We investigate the implications of the current data regarding large-scale anisotropies, including higher order multipoles, and we examine to what extent they can be used to shed some light on the origin of UHECRs and constrain the astrophysical and/or physical parameters of the source scenarios. We investigate the possibility of observing an associated anisotropy of the UHECR composition and discuss the potential benefit of a good determination of the composition and of the separation of the different nuclear components. We also discuss the interest and relevance of observing the UHECR sky with larger exposure future observatories.
Methods.
We simulated realistic UHECR sky maps for a wide range of astrophysical scenarios satisfying the current observational constraints, taking into account the energy losses and the photo-dissociation of the UHE protons and nuclei, as well as their deflexions by intervening magnetic fields. We investigated scenarios in which the UHECR source distribution follows that of the galaxies in the Universe (with possible biases), varying the UHECR source composition and spectrum, as well as the source density and the magnetic field models. For each of them, we simulated 300 realizations of independent datasets corresponding to various assumptions for the statistics and sky coverage, and we applied similar analyses as those used by the Auger collaboration for the search of large-scale anisotropies.
Results.
We find the following. First, reproducing the amplitude of the first-order (dipole) anisotropy observed in the Auger data, as well as its evolution as a function of energy, is relatively easy within our general assumptions. Second, this general agreement can be obtained with different sets of assumptions on the astrophysical and physical parameters, and thus it cannot be used, at the present stage, to derive strong constraints on the UHECR source scenarios or draw model-independent constraints on the various parameters individually. Third, the actual direction of the dipole modulation reconstructed from the Auger data, in the energy bin where the signal is most significant, appears highly unnatural in essentially all scenarios investigated, and this calls for their main assumptions to be reconsidered, either regarding the source distribution itself or the assumed magnetic field configuration, especially in the Galaxy. Fourth, the energy evolution of the reconstructed dipole direction contains potentially important information, which may become constraining for specific source models when larger statistics is collected. Fifth, for such high-statistics datasets, most of our investigated scenarios predict a significant quadrupolar modulation, especially if the light component of UHECRs can be extracted from the all-particle dataset. Sixth, except for protons, the energy range in which the GZK horizon strongly reduces is a key target for anisotropy searches for each given nuclear species. Seventh, although a difference in the average composition of the UHECRs in regions having a different count rate is naturally expected in our models, it is unlikely that the composition anisotropy recently reported by Auger can be explained by this effect, unless the reported amplitude is a strong positive statistical fluctuation of an intrinsically weaker signal.
Context. Various signals of anisotropy of the ultra-high-energy cosmic rays (UHECRs) have recently been reported, whether at large angular scales, with a dipole modulation in right ascension observed ...in the data of the Pierre Auger observatory ( Auger ), as discussed in the first paper accompanying the present one, or at intermediate angular scales, with flux excesses identified in specific directions by Auger and the Telescope Array (TA) Collaborations. Aims. We investigated the implications of the current data regarding these intermediate scale anisotropies, and examined to what extent they can be used to shed light on the origin of UHECRs, and constrain the astrophysical and/or physical parameters of the viable source scenarios. We also investigated what could be learnt from the study of the evolution of the various UHECR anisotropy signals, and discussed the expected benefit of an increased exposure of the UHECR sky using future observatories. Methods. We simulated realistic UHECR sky maps for a wide range of astrophysical scenarios satisfying the current observational constraints, with the assumption that the UHECR source distribution follows that of the galaxies in the Universe, also implementing possible biases towards specific classes of sources. In each case, several scenarios were explored with different UHECR source compositions and spectra, a range of source densities and different models of the Galactic magnetic field. We also implemented the Auger sky coverage, and explored various levels of statistics. For each scenario, we produced 300 independent datasets on which we applied similar analyses as those recently used by the Auger Collaboration, searching for flux excesses through either blind or targeted searches and quantifying correlations with predefined source catalogues through a likelihood analysis. Results. We find the following. First, with reasonable choices of the parameters, the investigated astrophysical scenarios can easily account for the significance of the anisotropies reported by Auger , even with large source densities. Second, the direction in which the maximum flux excess is found in the Auger data differs from the region where it is found in most of our simulated datasets, although an angular distance as large as that between the Auger direction and the direction expected from the simulated models at infinite statistics, of the order of ∼20°, occurs in ∼25% of the cases. Third, for datasets simulated with the same underlying astrophysical scenario, and thus the same actual UHECR sources, the significance with which the isotropy hypothesis is rejected through the Auger likelihood analysis can be largest either when ‘all galaxies’ or when only ‘starburst’ galaxies are used to model the signal, depending on which model is used to model the Galactic magnetic field and the resulting deflections. Fourth, the study of the energy evolution of the anisotropy patterns can be very instructive and provide new astrophysical insight about the origin of the UHECRs. Fifth, the direction in which the most significant flux excess is found in the Auger dataset above 8 EeV appears to essentially disappear in the dataset above 32 EeV, and, conversely, the maximum excess at high energy has a much reduced significance in the lower energy dataset. Sixth, both of these appear to be very uncommon in the simulated datasets, which could point to a failure of some generic assumption in the investigated astrophysical scenarios, such as the dominance of one type of source with essentially the same composition and spectrum in the observed UHECR flux above the ankle. Seventh, given the currently observed level of anisotropy signals, a meaningful measurement of their energy evolution, say from 10 EeV to the highest energies, will require a significant increase in statistics and a new generation of UHECR observatories.
We develop a method for studying the anisotropy of a cosmic-ray angular distribution, using both the right ascension and the declination of the arrival directions. It generalises the full-sky ...coverage method of Sommers (2001, Astropart. Phys., 14, 271) to partial-sky coverage experiments. When the angular distribution consists of a dipolar modulation of an otherwise isotropic flux, the method allows one to reconstruct the dipole amplitude and the dipole orientation in 3D space. We analyse the statistical properties of the method in detail, introducing the concept of reconstruction power, and show that it is generally more powerful than the standard Rayleigh analysis in right ascension. We clarify the link between the traditionally-used first harmonic amplitude and the true physical dipole amplitude, and we investigate the variation of the reconstruction powers as a function of the dipole orientation. We illustrate the method by computing the amplitude and angular reconstruction powers of the Pierre Auger Observatory, with the Southern site alone and with both Southern and Northern sites. In this particular case, we find that with an additional similar site in the Northern hemisphere the time needed for the method to reveal a significant departure from an isotropic cosmic-ray distribution would be reduced by a factor of about eight.
We present the results of three searches for correlations between ultra-high energy cosmic ray events (UHECRs) measured by Telescope Array and the Pierre Auger Observatory and high-energy neutrino ...candidate events from IceCube. Two cross-correlation analyses of UHECRs are done: one with 28 "cascades" from the IceCube 'high-energy starting events' sample and the other one with 12 high-energy "tracks". The angular separation between the arrival directions of neutrinos and UHECRs is scanned. The same events are also used in a separate search stacking the neutrino arrival directions and using a maximum likelihood approach. We assume that UHECR magnetic deflections are inversely proportional to the energy with values 3°, 6° and 9° at 100 EeV to account for the various scenarios of the magnetic field strength and UHECR charges. A similar analysis is performed on stacked UHECR arrival directions and the IceCube 4-year sample of through-going muon-track events that was optimized for neutrino point source searches.
We present the results of three searches for correlations between ultra-high energy cosmic ray events measured by Telescope Array and the Pierre Auger Observatory and high-energy neutrino candidate ...events from IceCube. Two cross-correlation analyses of ultra-high energy cosmic rays are done: one with 39 “cascades” from the IceCube “high-energy starting events” sample and the other one with 16 high-energy “tracks”. The angular separation between the arrival directions of neutrinos and UHECRs is scanned. The same events are also used in a separate search stacking the neutrino arrival directions and using a maximum likelihood approach. We assume that UHECR magnetic deflections are inversely proportional to the energy with values 3∘, 6∘ and 9∘ at 100 EeV to account for the uncertainties in the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube 4-year sample of through-going muon-track events that was optimized for neutrino point source searches.
We consider the energy spectrum of cosmic-rays (CRs) from a purely phenomenological point of view and investigate the possibility that they all be produced by the same type of sources with a single ...power-law spectrum, in $E^{-x}$, from thermal to ultra-high energies. We show that the relative fluxes of the Galactic (GCR) and extra-galactic (EGCR) components are compatible with such a holistic model, provided that the index of the source spectrum be $x \simeq 2.23\pm 0.07$. This is compatible with the best-fit indices for both GCRs and EGCRs, assuming that their source composition is the same, which is indeed the case in a holistic model. It is also compatible with theoretical expectations for particle acceleration at relativistic shocks.
ABSTRACT
Addressing the origin of the astrophysical neutrino flux observed by IceCube is of paramount importance. Gamma-Ray Bursts (GRBs) are among the few astrophysical sources capable of achieving ...the required energy to contribute to such neutrino flux through pγ interactions. In this work, ANTARES data have been used to search for upward going muon neutrinos in spatial and temporal coincidence with 784 GRBs occurred from 2007 to 2017. For each GRB, the expected neutrino flux has been calculated in the framework of the internal shock model and the impact of the lack of knowledge on the majority of source redshifts and on other intrinsic parameters of the emission mechanism has been quantified. It is found that the model parameters that set the radial distance where shock collisions occur have the largest impact on neutrino flux expectations. In particular, the bulk Lorentz factor of the source ejecta and the minimum variability time-scale are found to contribute significantly to the GRB-neutrino flux uncertainty. For the selected sources, ANTARES data have been analysed by maximizing the discovery probability of the stacking sample through an extended maximum-likelihood strategy. Since no neutrino event passed the quality cuts set by the optimization procedure, 90 per cent confidence level upper limits (with their uncertainty) on the total expected diffuse neutrino flux have been derived, according to the model. The GRB contribution to the observed diffuse astrophysical neutrino flux around 100 TeV is constrained to be less than 10 per cent.
A
bstract
The ANTARES neutrino telescope has an energy threshold of a few tens of GeV. This allows to study the phenomenon of atmospheric muon neutrino disappearance due to neutrino oscillations. In ...a similar way, constraints on the 3+1 neutrino model, which foresees the existence of one sterile neutrino, can be inferred. Using data collected by the ANTARES neutrino telescope from 2007 to 2016, a new measurement of Δ
m
32
2
and
θ
23
has been performed — which is consistent with world best-fit values — and constraints on the 3+1 neutrino model have been derived.