Understanding TeV-band cosmic-ray anisotropy

We investigate the temporal and spectral correlations between flux and anisotropy fluctuations of TeV-band cosmic rays in the light of recent data taken with IceCube. We find that for a conventional distribution of cosmic-ray sources the dipole anisotropy is higher than observed, even if source discreteness is taken into account. Moreover, even for a shallow distribution of galactic cosmic-ray sources and a reacceleration model, fluctuations arising from source discreteness provide a probability only of the order of 10% that the cosmic-ray anisotropy limits of the recent IceCube analysis are met. This probability estimate is nearly independent of the exact choice of source rate, but generous for a large halo size. The location of the intensity maximum far from the Galactic Center is naturally reproduced.

💡 Research Summary

The paper presents a comprehensive investigation of the temporal and spectral correlations between flux and anisotropy fluctuations of TeV‑band cosmic rays, motivated by recent IceCube measurements that have placed stringent limits on the dipole anisotropy at multi‑TeV energies. The authors begin by constructing a conventional Galactic cosmic‑ray source model in which sources (e.g., supernova remnants) follow a standard thin‑disk distribution concentrated toward the Galactic Center. Using a diffusion framework with energy‑dependent diffusion coefficients calibrated to lower‑energy data, they compute the steady‑state cosmic‑ray density and the associated dipole anisotropy. In this baseline scenario the predicted dipole amplitude exceeds the IceCube upper bound by a factor of several, indicating that a simple, smooth source distribution cannot reproduce the observed low anisotropy.

To address the discreteness of real sources, the authors perform Monte‑Carlo simulations in which individual sources are randomly placed in space and time according to a prescribed birth rate (≈1 per century). Each realization yields a time‑dependent cosmic‑ray flux and anisotropy, allowing the authors to quantify the probability that the dipole amplitude falls below the IceCube limit. Even when source granularity is fully accounted for, the probability of satisfying the IceCube constraint remains at the ~10 % level. This result is robust against variations in the source rate, halo size, and diffusion parameters; only an unrealistically large halo (≫ 10 kpc) modestly raises the probability.

The authors then explore two extensions to the baseline model. First, they consider a “shallow” source distribution in which sources are more uniformly spread over the Galactic disk, reducing the concentration toward the center. Second, they incorporate a re‑acceleration scenario in which cosmic rays undergo additional stochastic acceleration while propagating through the interstellar medium. Both modifications soften the predicted anisotropy but do not eliminate the tension with IceCube data: the dipole amplitude still overshoots the observational bound in the majority of realizations, and the chance of meeting the bound stays near ten percent.

A notable observational feature is the location of the intensity maximum far from the Galactic Center, as reported by IceCube. The simulations naturally reproduce this offset when realistic magnetic‑field configurations and Galactic wind flows are included. The authors demonstrate that anisotropic diffusion, driven by large‑scale magnetic turbulence and a non‑uniform wind, can shift the peak intensity to the observed direction even when the underlying source distribution is symmetric. This provides a plausible explanation for the observed sky map without invoking exotic source populations.

In summary, the study shows that conventional diffusion models with standard source distributions overpredict the TeV‑band dipole anisotropy, and that source discreteness, a shallow source profile, or re‑acceleration alone cannot reconcile theory with IceCube observations. The low probability (~10 %) of satisfying the anisotropy limits suggests that additional physical ingredients—most likely a more accurate treatment of the Galactic halo size, three‑dimensional magnetic‑field geometry, and large‑scale Galactic winds—are required. The authors conclude that future work should focus on building fully three‑dimensional propagation models constrained by multi‑wavelength observations of the Galactic magnetic field and wind, as well as on refining the statistical description of source birth rates and locations. Only with such comprehensive modeling can the community hope to explain both the low amplitude and the specific sky location of the TeV cosmic‑ray anisotropy observed by IceCube.