⚛️ high-energy theory

Reconstructing the Auger UHECR Dipole: A Hybrid Analysis of 4LAC AGNs and Nearby Starburst Galaxies

This paper demonstrates that neither gamma-ray active galactic nuclei nor nearby starburst galaxies alone can explain the ultra-high-energy cosmic ray dipole observed by the Pierre Auger Observatory, necessitating a hybrid model where starburst galaxies dominate the anisotropic signal while active galactic nuclei provide a subdominant contribution.

Original authors: Swaraj Pratim Sarmah, Umananda Dev Goswami, Sovan Chakraborty

Published 2026-08-25

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Swaraj Pratim Sarmah, Umananda Dev Goswami, Sovan Chakraborty

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The sky above us is constantly bombarded by invisible particles of unimaginable power. These are ultra-high-energy cosmic rays, subatomic messengers that travel across the universe at nearly the speed of light. For decades, scientists have struggled to answer a simple but profound question: where do they come from? Unlike light, which travels in straight lines, these heavy particles are easily knocked off course by magnetic fields that thread through space, making their arrival directions on Earth seem random. However, recent observations have revealed a subtle pattern: a slight excess of these particles arriving from one general direction in the sky, a "dipole" anisotropy. This faint tilt in the cosmic rain offers the first real clue that these particles originate from outside our own galaxy, but pinpointing the exact cosmic factories responsible has remained elusive.

A team of researchers has now taken a fresh look at this mystery by testing two leading theories about the sources of these cosmic rays. One theory points to active galactic nuclei, the supermassive black holes at the centers of distant galaxies that spew out jets of energy. The other theory suggests that starburst galaxies, which are undergoing furious bursts of star formation and supernova explosions, are the primary accelerators. To solve the puzzle, the scientists did not just guess; they built a detailed simulation of the universe. They gathered data on thousands of known active galaxies and combined it with information on the few nearby starburst galaxies that are close enough to influence the cosmic rays reaching Earth. Crucially, they accounted for the fact that as these particles travel across billions of light-years, they lose energy by colliding with the faint afterglow of the Big Bang, a process that effectively filters out the most distant sources.

When the researchers tested each group of galaxies on its own, the results were clear: neither could explain the pattern seen by the Pierre Auger Observatory. If the cosmic rays came only from the distant active galaxies, the predicted tilt in the sky would be far too strong and point in the wrong direction. If they came only from the nearby starburst galaxies, the signal was still too strong and misaligned. The universe, it seemed, was not playing by the rules of a single source. The team then tried a different approach, mixing the contributions from both types of galaxies along with a background of particles coming from everywhere else. This hybrid model worked. By combining the distant active galaxies with the nearby starburst galaxies, the researchers were able to perfectly recreate the observed tilt in the cosmic ray sky.

The most successful version of this model, specifically for the energy range between 8 and 16 exa-electronvolts, revealed a specific recipe for the cosmic rays we detect. The analysis showed that the signal is not dominated by one type of source but is a blend. In this energy range, the nearby starburst galaxies provide the largest share of the directional signal, contributing about 45 percent of the anisotropic component. The distant active galaxies add another 23 percent, while the remaining 32 percent comes from a diffuse background of particles whose directions have been scrambled by magnetic fields. This mixture produced a predicted tilt in the sky that matched the Auger observations with remarkable precision, differing by less than seven degrees in direction and matching the strength of the signal almost exactly.

As the researchers looked at even higher energies, the story shifted slightly but consistently. The influence of the nearby starburst galaxies grew stronger, while the contribution from the distant active galaxies remained a steady, smaller partner. This makes sense physically: the most energetic particles lose energy so quickly during their journey that only the very closest sources can send them to Earth without their signal fading away. The distant active galaxies, while numerous, are too far away to dominate the highest-energy events. The study concludes that the large-scale pattern of cosmic rays is not the signature of a single cosmic monster or a lone galaxy, but the combined voice of the universe's most energetic neighborhoods, speaking in a chorus of nearby starbursts and distant black holes.

Technical Summary: Reconstructing the Auger UHECR Dipole via a Hybrid AGN+SBG Analysis

Problem Statement
The Pierre Auger Observatory has confirmed a large-scale dipole anisotropy in the arrival directions of ultra-high-energy cosmic rays (UHECRs) above 8 EeV, with an amplitude of approximately 6.5% and a direction roughly 125° away from the Galactic Center. This observation establishes an extragalactic origin for these particles. However, the specific astrophysical sources responsible for this dipole remain unidentified. While Active Galactic Nuclei (AGNs) and Starburst Galaxies (SBGs) are physically motivated candidates—satisfying the Hillas criterion and hosting extreme stellar explosions, respectively—previous studies have not definitively determined whether either population alone, or a combination thereof, can reproduce the observed dipole amplitude and direction when accounting for propagation effects.

Methodology
The authors construct a hybrid analysis framework to reconstruct the UHECR dipole using two primary source populations:

AGN Sample: The Fermi-LAT Fourth AGN Catalog (4LAC DR3) is utilized, providing 1,806 valid sources with measured sky positions, redshifts, and gamma-ray energy fluxes (100 MeV–100 GeV). These fluxes serve as proxies for UHECR acceleration capability.
SBG Sample: A restricted sample of four nearby starburst galaxies (M82, NGC 253, NGC 4945, and M83) within 5 Mpc is selected, as their proximity minimizes propagation attenuation.

The analysis incorporates a mixed-composition propagation model based on Auger data and the EPOS-LHC hadronic interaction model. Key physical effects included are:

GZK Attenuation: Energy-dependent losses due to interactions with the Cosmic Microwave Background (CMB), including photopion production, pair production, and photodisintegration for nuclei.
Magnetic Deflections: Rigidity-dependent smearing of arrival directions due to the Galactic Magnetic Field (GMF), modeled with a sky-averaged uniform smearing angle (baseline $\sigma_0 = 5^\circ$ at 60 EeV).
Hybrid Vector Construction: The total dipole vector is modeled as a linear combination of the AGN dipole vector, the SBG dipole vector, and an isotropic background component ( $\vec{d}_{total} = f_{AGN}\vec{d}_{AGN} + f_{SBG}\vec{d}_{SBG}$ ). The isotropic fraction is defined as $f_{iso} = 1 - f_{AGN} - f_{SBG}$ .

The model is tested across four energy intervals (4–8, 8–16, 16–32, and >32 EeV) by minimizing a joint $\chi^2$ metric that accounts for discrepancies in both dipole amplitude and angular separation from the observed Auger direction.

Key Results

Failure of Single Populations: Analysis of single-source populations reveals that neither AGNs nor SBGs alone can simultaneously match the observed Auger dipole amplitude and direction.
- The AGN-only model in the 8–16 EeV bin predicts a dipole amplitude of 51.87% (far exceeding the observed ~5.8%) with a direction of $(167.3^\circ, +73.6^\circ)$ .
- The SBG-only model predicts an amplitude of 30.49% with a direction of $(48.0^\circ, -65.9^\circ)$ .
Hybrid Model Success: A hybrid model combining both populations and an isotropic background successfully reproduces the Auger observations.
- Best-Fit (8–16 EeV): The optimal solution requires a 23% AGN contribution and a 45% SBG contribution, with the remaining 33% being isotropic. This yields a reconstructed dipole amplitude of 5.03% and a direction of $(104.1^\circ, -25.1^\circ)$ , resulting in an angular separation of only 7.10° from the observed Auger dipole.
- Energy Evolution: As energy increases, the SBG contribution remains dominant (peaking at 47% in the 16–32 EeV bin), while the isotropic fraction decreases. The reconstructed dipole direction systematically shifts toward southern declinations, consistent with the increasing influence of nearby southern SBGs (specifically NGC 253 and NGC 4945) as GZK attenuation suppresses distant sources.
- Low Energy Behavior: At the lowest energy bin (4–8 EeV), the model is dominated by the isotropic component (97%), with negligible contributions from AGNs (0%) and SBGs (3%), reflecting stronger magnetic deflections and contributions from more distant sources.

Significance and Claims
The paper claims that a single source class is insufficient to explain the large-scale UHECR anisotropy observed by the Auger Observatory. The primary significance of this work is the demonstration that a hybrid framework, incorporating gamma-ray AGNs, nearby starburst galaxies, and an isotropic background, can simultaneously reproduce the observed dipole amplitude and directional evolution across multiple energy intervals.

The authors emphasize that their analysis differs from intermediate-scale clustering searches (which focus on localized flux excesses) by addressing the global large-scale dipole. They conclude that nearby southern starburst galaxies provide the dominant anisotropic contribution, particularly at higher energies, while AGNs provide a subdominant contribution alongside a quasi-isotropic background. The results suggest that the observed dipole is a composite signal driven by the interplay of local source dominance and propagation effects, rather than a single population of extragalactic accelerators.

More like this