No Arabic abstract
We present a novel method to search for structures of coherently aligned patterns in ultra-high energy cosmic-ray arrival directions simultaneously across the entire sky. This method can be used to obtain information on the Galactic magnetic field, in particular the integrated component perpendicular to the line of sight, from cosmic-ray data only. Using a likelihood-ratio approach, neighboring cosmic rays are related by rotatable, elliptically shaped density distributions and the significance of their alignment with respect to circular distributions is evaluated. In this way, a vector field tangential to the celestial sphere is fitted which approximates the local deflections in cosmic magnetic fields if significant deflection structures are detected. The sensitivity of the method is evaluated on the basis of astrophysical simulations of the ultra-high energy cosmic-ray sky, where a discriminative power between isotropic and signal-induced scenarios is found.
The sources of ultrahigh-energy cosmic rays (UHECRs) have been difficult to catch. It was recently pointed out that while sources of UHECR protons exhibit anisotropy patterns that become denser and compressed with rising energy, nucleus-emitting-sources give rise to a cepa stratis (onion-like) structure with layers that become more distant from the source position with rising energy. The peculiar shape of the hot spots from nucleus-accelerators is steered by the competition between energy loss during propagation and deflection on the Galactic magnetic field (GMF). Here, we run a full-blown simulation study to accurately characterize the deflections of UHECR nuclei in the GMF. We show that while the cepa stratis structure provides a global description of anisotropy patterns produced by UHECR nuclei en route to Earth, the hot spots are elongated depending on their location in the sky due to the regular structure of the GMF. We demonstrate that with a high-statistics sample at the high-energy-end of the spectrum, like the one to be collected by NASAs POEMMA mission, the energy dependence of the hot-spot contours could become a useful observable to identify the nuclear composition of UHECRs. This new method to determine the nature of the particle species is complementary to those using observables of extensive air showers, and therefore is unaffected by the large systematic uncertainties of hadronic interaction models.
Anisotropy in the arrival directions of cosmic rays with energies above 10$^{17}$eV is studied using data from the Akeno 20 km$^2$ array and the Akeno Giant Air Shower Array (AGASA), using a total of about 117,000 showers observed during 11 years. In the first harmonic analysis, we have found strong anisotropy of $sim$ 4% around 10$^{18}$eV, corresponding to a chance probability of 0.2%. With two dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region.
We present a new approach for the identification of ultra-high energy cosmic rays from sources using dynamic graph convolutional neural networks. These networks are designed to handle sparsely arranged objects and to exploit their short- and long-range correlations. Our method searches for patterns in the arrival directions of cosmic rays, which are expected to result from coherent deflections in cosmic magnetic fields. The network discriminates astrophysical scenarios with source signatures from those with only isotropically distributed cosmic rays and allows for the identification of cosmic rays that belong to a deflection pattern. We use simulated astrophysical scenarios where the source density is the only free parameter to show how density limits can be derived. We apply this method to a public data set from the AGASA Observatory.
The IceCube Neutrino Observatory has accumulated a total of 318 billion cosmic-ray induced muon events between May 2009 and May 2015. This data set was used for a detailed analysis of the cosmic-ray arrival direction anisotropy in the TeV to PeV energy range. The observed global anisotropy features large regions of relative excess and deficit, with amplitudes on the order of $10^{-3}$ up to about 100 TeV. A decomposition of the arrival direction distribution into spherical harmonics shows that most of the power is contained in the low-multipole ($ellleq 4$) moments. However, higher multipole components are found to be statistically significant down to an angular scale of less than $10^{circ}$, approaching the angular resolution of the detector. Above 100 TeV, a change in the morphology of the arrival direction distribution is observed, and the anisotropy is characterized by a wide relative deficit whose amplitude increases with primary energy up to at least 5,PeV, the highest energies currently accessible to IceCube. No time dependence of the large- and small-scale structures is observed in the six-year period covered by this analysis. The high-statistics data set reveals more details on the properties of the anisotropy and is potentially able to shed light on the various physical processes that are responsible for the complex angular structure and energy evolution.
Energy-dependent patterns in the arrival directions of cosmic rays are searched for using data of the Pierre Auger Observatory. We investigate local regions around the highest-energy cosmic rays with $E geq 6 cdot 10^{19}$ eV by analyzing cosmic rays with energies above $E = 5 cdot 10^{18}$ eV arriving within an angular separation of approximately $15{deg}$. We characterize the energy distributions inside these regions by two independent methods, one searching for angular dependence of energy-energy correlations and one searching for collimation of energy along the local system of principal axes of the energy distribution. No significant patterns are found with this analysis. The comparison of these measurements with astrophysical scenarios can therefore be used to obtain constraints on related model parameters such as strength of cosmic-ray deflection and density of point sources.