Abstract:We argue that if ultrahigh-energy (E > 10^10 GeV) cosmic rays are heavy nuclei (as indicated by existing data), then the pointing of cosmic rays to their nearest extragalactic sources is expected for 10^10.6 < E/GeV < 10^11. This is because for a nucleus of charge Ze and baryon number A, the bending of the cosmic ray decreases as Z/E with rising energy, so that pointing to nearby sources becomes possible at highest energies. In addition, the maximum energy of acceleration capability of the sources grows linearly in Z, while the energy loss per distance traveled decreases with increasing A. Each of these two points tend to favor heavy nuclei at the highest energies. The traditional bi-dimensional analyses, which simultaneously reproduce Auger data on the spectrum and nuclear composition, may not be capable of incorporating the relative importance of these phenomena. In this paper we propose a multi-dimensional reconstruction of the individual emission spectra (in E, direction, and cross-correlation with nearby putative sources) to study the hypothesis that primaries are heavy nuclei subject to GZK photo-disintegration, and to determine the nature of the extragalactic sources. More specifically, we propose to combine information on nuclear composition and arrival direction to associate a potential clustering of events with a 3-dimensional position in the sky. Actually, both the source distance and maximum emission energy can be obtained through a multi-parameter likelihood analysis to accommodate the observed nuclear composition of each individual event in the cluster. We show that nucleus-emitting-sources exhibit a cepa stratis structure on Earth which could be pealed off by future space-missions, such as POEMMA. We also show that metal-rich starburst galaxies are highly-plausible candidate sources, and we use them as an explicit example of our proposed multi-dimensional analysis.

Generalized statistical mechanics of cosmic rays: Application to positron-electron spectral indices

Possible Primary-Electron-spectrum Hardening at $\sim 240$ GeV: Implications for AMS-02 Observation and the Physical Origin of Cosmic Ray Excesses

Reanalysis of the systematic uncertainties in the cosmic-ray antiproton flux

New estimation of the spectral index of high-energy cosmic rays as determined by the Compton-Getting anisotropy

Reanalysis of the Systematic Uncertainties in Cosmic-Ray Antiproton Flux

Intermediate mass and heavy Galactic cosmic-ray nuclei: the case of new AMS-02 measurements

Simple and less hadronic-model-dependent method for cosmic-ray composition determination

Discriminating Different Scenarios to Account for the Cosmic E^± Excess by Synchrotron and Inverse Compton Radiation

Diffuse Gamma-Rays Produced in Cosmic-Ray Interactions and the TeV-band Spectrum of RX J1713.7-3946

Astrophysical Background and Dark Matter Implication Based on Latest AMS-02 Data

Origins of cosmic positrons and electrons

Statistical Tools for Analyzing the Cosmic Ray Energy Spectrum

Ultrarelativistic generalized Lorentzian thermodynamics and the differential cosmic ray energy flux

Spatial Dependent Diffusion of Cosmic Rays and the Excess of Primary Electrons Derived from High Precision Measurements by AMS-02

Dark Matter with Multiannihilation Channels and the Ams-02 Positron Excess and Antiproton Data

Electron and Positron Fluxes in Primary Cosmic Rays measured with the Alpha Magnetic Spectrometer on the International Space Station

Ams-02 Positron Excess: New Bounds on Dark Matter Models and Hint for Primary Electron Spectrum Hardening

Cosmic rays: the spectrum and chemical composition from $10^{10}$ to $10^{20}$ eV

Elastic scattering at 7 TeV and high energy cross section for cosmic ray studies

Cosmic Mass Spectrometer

QCD parameters and SM-high precisions from $e^+e^-\to$ Hadrons : Summary