Abstract:For many decades, ultrahigh energy charged particles of unknown origin that can be observed from the ground have been a puzzle for particle physicists and astrophysicists. As an attempt to discriminate among several possible production scenarios, astrophysicists try to test the statistical isotropy of the directions of arrival of these cosmic rays. At the highest energies, they are supposed to point toward their sources with good accuracy. However, the observations are so rare that testing the distribution of such samples of directional data on the sphere is nontrivial. In this paper, we choose a non-parametric framework that makes weak hypotheses on the alternative distributions and allows in turn to detect various and possibly unexpected forms of anisotropy. We explore two particular procedures. Both are derived from fitting the empirical distribution with wavelet expansions of densities. We use the wavelet frame introduced by [SIAM J. Math. Anal. 38 (2006b) 574-594 (electronic)], the so-called needlets. The expansions are truncated at scale indices no larger than some J(star), and the L-p distances between those estimates and the null density are computed. One family of tests (called MULTIPLE) is based on the idea of testing the distance from the null for each choice of J = 1, ... , J(star), whereas the so-called PLUGIN approach is based on the single full J(star) expansion, but with thresholded wavelet coefficients. We describe the practical implementation of these two procedures and compare them to other methods in the literature. As alternatives to isotropy, we consider both very simple toy models and more realistic nonisotropic models based on Physics-inspired simulations. The Monte Carlo study shows good performance of the MULTIPLE test, even at moderate sample size, for a wide sample of alternative hypotheses and for different choices of the parameter J(star). On the 69 most energetic events published by the Pierre Auger Collaboration, the needlet-based procedures suggest statistical evidence for anisotropy. Using several values for the parameters of the methods, our procedures yield p-values below 1%, but with uncontrolled multiplicity issues. The flexibility of this method and the possibility to modify it to take into account a large variety of extensions of the problem make it an interesting option for future investigation of the origin of ultrahigh energy cosmic rays.

Techniques for Searching for Anisotropy in the Arrival Directions of Ultra-High Energy Cosmic Rays

Simulations for a Next-Generation UHECR Observatory

What can be learnt from UHECR anisotropies observations Paper II: intermediate-scale anisotropies

Search for correlations of high-energy neutrinos and ultra-high energy cosmic rays

A Search for Correlation of Ultra-High Energy Cosmic Rays with IRAS-PSCz and 2MASS-6dF Galaxies

Constraints on UHECR Sources and Extragalactic Magnetic Fields from Directional Anisotropies

Analysis of the Arrival Directions of Ultrahigh Energy Cosmic Rays

Data-driven analysis for understanding ultrahigh energy cosmic ray source spectra

Testing the Correlation of Ultra-High Energy Cosmic Rays with High Redshift Sources

Testing the Isotropy of High Energy Cosmic Rays Using Spherical Needlets

Estimation Prospects of the Source Number Density of Ultra-high-energy Cosmic Rays

Anisotropies of Diffusive Ultra-high Energy Cosmic Rays in $f(R)$ Gravity Theory

Modification of the Dipole in Arrival Directions of Ultra-high-energy Cosmic Rays due to the Galactic Magnetic Field

Reconstructing Galactic magnetic fields from local measurements for backtracking ultra-high-energy cosmic rays

Identifying the nature of high energy Astroparticles

Honing cross-correlation tools for inference on ultra-high-energy cosmic-ray composition

Particle and Astrophysics Aspects of Ultrahigh Energy Cosmic Rays

Images of the Ultra-High Energy Cosmic Rays from Point Sources

Measurement of Energy Spectrum and Elemental Composition of PeV Cosmic Rays: Open Problems and Prospects

High energy cosmic ray and neutrino astronomy