Origin of ultra-high-energy galactic cosmic rays: The isotropy problem

We study the propagation of ultra-high-energy cosmic rays (UHECRs) in the Galaxy, concentrating on the energy range below the ankle in the spectrum at 4 EeV. A Monte Carlo method, based on analytical solutions to the time-dependent diffusion problem, is used to account for intermittency by placing sources at random locations. Assuming a source population that scales with baryon mass density or star formation (e.g., long GRB), we derive constraints arising from intermittency and the observational limits on the composition and anisotropy. It is shown that the composition and anisotropy at 10¹⁸ eV are difficult to reproduce and require that either (1) the particle mean free path is much smaller than a gyroradius, implying the escape time is very long, (2) the composition is heavier than suggested by recent Auger data, (3) the ultra-high-energy sub-ankle component is mostly extragalactic, or (4) we are living in a rare lull in the UHECR production, and the current UHECR intensity is far below the Galactic time average. We therefore recommend a strong observational focus on determining the UHECR composition around 10¹⁸ eV.