Hillas meets Eddington: The case for blazars as ultra-high-energy neutrino sources

This paper presents a new physically motivated, multi-zone leptohadronic model demonstrating that sub-Eddington blazar jets can efficiently accelerate protons to EeV energies to produce the ultra-high-energy neutrino fluxes observed by IceCube and KM3NeT, such as those from TXS 0506+056 and PKS 0605-085, while remaining consistent with multi-wavelength emission data and energy budget constraints.

The Gist

Imagine the universe as a vast, dark ocean. For decades, scientists have been trying to figure out where the most energetic particles in existence come from. These particles, called cosmic rays, are like tiny bullets moving at nearly the speed of light. When they hit the Earth's atmosphere, they create a shower of other particles, including neutrinos—ghostly particles that can pass through planets without stopping.

For a long time, we knew these "ghosts" existed, but we didn't know who was throwing them. This paper is like a detective story where the authors finally catch the suspect: Blazars.

Here is the story of their discovery, explained simply.

The Suspect: The Blazar

Think of a Blazar as a cosmic lighthouse. It's a supermassive black hole at the center of a galaxy that is shooting out a massive, high-speed jet of particles. Usually, these jets point in random directions, but a Blazar is special because its jet is pointed almost directly at Earth. It's like someone shining a laser pointer straight into your eyes.

For years, scientists thought these jets were mostly made of electrons (tiny, light particles). But electrons are like ping-pong balls; they can't carry enough energy to create the ultra-powerful neutrinos we see. The authors of this paper asked: "What if the jet is actually carrying heavy cannonballs (protons) instead?"

The New Theory: The "Hillas meets Eddington" Model

The authors built a new computer model to test this idea. They combined two famous physics concepts:

1. Hillas: How big and strong a magnetic field needs to be to accelerate a particle.
2. Eddington: The maximum amount of energy a black hole can safely output before it blows itself apart.

The Analogy of the Jet:
Imagine the jet is a giant, expanding highway.

• Near the Black Hole (The Start): The highway is narrow and crowded with magnetic fields. It's like a high-pressure water hose. Here, the magnetic fields act like a giant slingshot, grabbing protons and flinging them to incredible speeds.
• Further Out (The Expansion): As the jet moves away, it widens and slows down. The magnetic "slingshot" weakens, and the jet becomes more like a river of moving water (kinetic energy).

The authors' model shows that in this "magnetic slingshot" zone, protons can be accelerated to energies so high they are quadrillions of times more energetic than the particles in the Large Hadron Collider.

The Smoking Gun: Connecting the Dots

The real breakthrough of this paper is how it explains three different things happening at once, which previous models couldn't do:

1. The Neutrino: The super-fast protons smash into light particles (photons) floating around the black hole. This collision creates a neutrino. The authors show that the energy of these neutrinos matches what the IceCube detector in Antarctica has seen.
2. The Gamma Rays (The Flash): When those protons get hit, they also glow with high-energy light (gamma rays). The authors show that this "proton glow" explains the bright flashes of light we see from these galaxies.
3. The Optical Light (The Color): The electrons in the jet are also moving fast, but not as fast as the protons. They create the visible light (optical) we see.

The "Aha!" Moment:
Previous models struggled because they tried to explain the neutrinos and the light using only electrons, or they required the black hole to use more energy than it physically has (like trying to fill a bathtub with a firehose).
This new model says: "The protons do the heavy lifting for the neutrinos and the gamma rays, while the electrons handle the visible light." It fits perfectly within the energy budget of the black hole.

The Case Studies: Two Culprits Caught

The authors tested their model on two specific "suspects":

1. TXS 0506+056: This is the famous blazar that sent a neutrino to Earth in 2017, coinciding with a bright flash of light. The authors' model perfectly recreates this event. They found that a temporary "burst" of particles near the black hole (like a sudden traffic jam on our cosmic highway) explains why the neutrino and the light arrived together.
2. PKS 0605-085: This is a newer suspect. A different detector in the Mediterranean (KM3NeT) recently spotted a neutrino with even higher energy. The authors applied their model to this blazar and found it fits perfectly, predicting that the neutrino came from protons hitting dust clouds far out in the jet.

Why This Matters

This paper is a game-changer because it solves a puzzle that has been stuck for years. It tells us that:

• Blazars are the factories: They are the cosmic accelerators creating the most energetic particles in the universe.
• Protons are the stars: It's the heavy protons, not the light electrons, that are doing the heavy lifting.
• The future is bright: As our detectors get better (like the upcoming IceCube-Gen2), we will be able to see these "ghosts" more clearly, confirming that our cosmic lighthouses are indeed the engines of the universe's most violent fireworks.

In short: The authors built a new map of the cosmic highway. They proved that the black holes are powerful enough to launch proton "bullets" that create the ghostly neutrinos we detect, solving the mystery of where the universe's most energetic particles come from.

Technical Summary

Here is a detailed technical summary of the paper "Hillas meets Eddington: The case for blazars as ultra-high-energy neutrino sources" by Rodrigues et al.

1. Problem Statement

The origin of Ultra-High-Energy (UHE) cosmic rays and neutrinos (energies $>100$ PeV) remains a major open question in astroparticle physics. While blazars (active galactic nuclei with jets aligned to Earth) are prime candidates, existing leptohadronic models face significant challenges:

• Energetic Constraints: Models where protons only reach sub-PeV energies often require proton powers exceeding the Eddington limit to produce the observed neutrino fluxes, which is physically unrealistic.
• Predictive Limitations: Standard models often rely on a "single-zone" approximation and phenomenological parameters, failing to self-consistently link particle acceleration to jet dynamics and energetics.
• Multi-messenger Tension: In sub-PeV models, protons contribute negligibly to $\gamma$-ray emission, weakening the causal link between observed $\gamma$-ray flares and neutrino events. Conversely, models where protons reach EeV energies often struggle to satisfy X-ray constraints due to electromagnetic cascades.

2. Methodology

The authors present a new extended leptohadronic model that moves beyond the single-zone approximation by simulating the continuous evolution of a blazar jet from the black hole base to the parsec scale.

• Jet Dynamics & Energetics:

  • The jet is modeled as sub-Eddington ($L_j < L_{\text{Edd}}$) and magnetically launched.
  • It evolves from a Poynting-flux dominated regime (near the black hole) to a kinetically dominated regime further out.
  • Energy conservation is enforced: $L_B(r) + L_k(r) = \text{const} < L_{\text{Edd}}$.
  • The jet accelerates as magnetic energy converts to kinetic energy, reaching bulk Lorentz factors $\Gamma_j \sim 30-40$ at the parsec scale.

• Particle Acceleration & Injection:

  • A constant fraction ($f_{NT} \sim 10^{-6} - 10^{-8}$) of thermal protons and electrons picked up from the ambient medium are continuously accelerated.
  • Acceleration Mechanism: Stochastic Fermi acceleration (second-order) in turbulent magnetic fields is the dominant mechanism in the inner jet.
  • Diffusion: The model tests diffusion scaling $D' \propto E^\delta$. The authors find that Kolmogorov-like diffusion ($\delta \approx 0.3$) is required to simultaneously explain optical synchrotron peaks and UHE proton acceleration, ruling out Bohm diffusion ($\delta=1$).
  • Spectral Index: A power-law index of $p \approx 1.8$ is adopted, consistent with feedback between particles and turbulence.

• Radiative Framework:

  • The model numerically solves time- and energy-dependent continuity equations using the AM3 code across 30 zones per decade in radius.
  • It accounts for interactions with external photon fields: the Broad Line Region (BLR) and the dusty torus.
  • It calculates full electromagnetic cascades (pair production, synchrotron, inverse Compton) and neutrino production via photohadronic interactions ($p\gamma$).

3. Key Contributions

• Unified Extended Jet Model: The paper bridges the gap between single-zone models and full jet hydrodynamics, deriving particle properties directly from jet physics (magnetic field profiles, turbulence, density) rather than fitting free parameters arbitrarily.
• Resolution of the "Hillas vs. Eddington" Tension: The model demonstrates that protons can reach EeV energies within a sub-Eddington jet power budget, provided acceleration occurs in the highly magnetized inner jet ($r \sim 0.1$ pc) where the magnetization $\sigma \gtrsim 1$.
• Diffusion Constraint: By jointly fitting optical, $\gamma$-ray, and neutrino data, the model constrains the diffusion coefficient scaling to $\delta \approx 0.3$, providing a physical constraint on turbulence in relativistic jets.
• Spatial Emission Mapping: The model predicts the specific spatial origin of different emission components along the jet (e.g., X-rays from cascades at 0.2 pc, optical from electrons at parsec scales).

4. Results

The model was applied to two specific sources:

A. TXS 0506+056 (IceCube Candidate)

• Neutrino Flux: Protons accelerated in the inner jet ($r < 0.1$ pc) interact with infrared photons from the dusty torus, producing a neutrino flux peaking at $\sim 100$ PeV. This is consistent with the 10-year IceCube point-source limits.
• Multiwavelength Fit:
  • $\gamma$-rays: Dominated by proton synchrotron emission (peaking at GeV energies), establishing a direct causal link between the 2017 $\gamma$-ray flare and the neutrino event.
  • Optical: Dominated by electron synchrotron emission from the parsec-scale jet.
  • X-rays: A combination of electron synchrotron and hadronic cascades.
• Flare Modeling: The 2017 flare is reproduced by a temporary increase in particle injection near the BLR ($r \approx 0.05$ pc), boosting the proton luminosity to $\sim 35\%$ of $L_{\text{Edd}}$ locally, while maintaining the global jet power below the Eddington limit.

B. PKS 0605-085 (KM3NeT Candidate)

• Applied to a blazar associated with a recent KM3NeT event ($>100$ PeV).
• The model predicts a neutrino peak at 200 PeV, consistent with the detected muon energy ($120^{+110}_{-60}$ PeV).
• Similar to TXS 0506+056, the $\gamma$-ray spectrum is dominated by proton synchrotron, and the neutrino production is driven by interactions with the dusty torus.

C. 3HSP J1528+2004 (Northern Sky Candidate)

• Applied to a source well below the IceCube horizon.
• The model successfully fits the data by selecting parameters that limit proton acceleration to the PeV range (due to faster magnetic field dissipation), demonstrating the model's flexibility to accommodate different source geometries and neutrino energies.

5. Significance

• Blazars as UHE Sources: The study provides strong evidence that blazars are efficient emitters of UHE neutrinos ($>100$ PeV), a regime where current detectors like IceCube have limited sensitivity but where the flux is theoretically viable.
• Future Detection: The predicted neutrino spectra are well within the reach of next-generation detectors like IceCube-Gen2 and KM3NeT, which are designed to probe the UHE range.
• Physical Constraints: The work imposes strict constraints on jet physics, specifically ruling out Bohm diffusion in favor of Kolmogorov-like turbulence and confirming that efficient UHE acceleration requires a transition from magnetically to kinetically dominated flow.
• Multi-messenger Consistency: The model resolves the tension between neutrino and $\gamma$-ray observations by showing that protons can dominate the high-energy $\gamma$-ray flux without violating X-ray constraints, provided the acceleration occurs in the specific inner-jet environment described.

In conclusion, the paper establishes a physically motivated, extended jet framework that successfully explains multimessenger data from blazars, identifying them as viable sources of the highest-energy cosmic particles.