🔭 astrophysics

On the maximum neutrino flux of blazars in the one-zone leptohadronic model

This paper demonstrates through analytical and numerical modeling that the one-zone leptohadronic blazar model cannot produce neutrino fluxes high enough to match IceCube observations due to X-ray constraints, suggesting that multi-zone scenarios or alternative production sites are necessary to explain the origin of high-energy neutrinos.

Original authors: Wei-Jian Li, Rui Xue, Ze-Rui Wang, Dingrong Xiong

Published 2026-05-26

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Wei-Jian Li, Rui Xue, Ze-Rui Wang, Dingrong Xiong

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

Imagine the universe is a vast, dark ocean, and hidden within it are cosmic lighthouses called blazars. These aren't normal lighthouses; they are supermassive black holes at the centers of distant galaxies, shooting out powerful jets of particles at nearly the speed of light. For years, astronomers have been puzzled by high-energy "ghost particles" called neutrinos that rain down on Earth from space. They suspect these ghosts are being created by these cosmic lighthouses, but they can't quite figure out how.

This paper is like a detective trying to solve that mystery using a specific set of rules. Here is the story of what they found, explained simply.

The Detective's Tool: The "One-Room" Theory

To understand how blazars work, scientists often use a model called the "one-zone leptohadronic model."

Think of a blazar's jet as a single, giant, glowing room. Inside this room, there are two types of particles:

Electrons (the "light workers"): They create the light we see (from radio waves to X-rays).
Protons (the "heavy hitters"): They are the ones that smash into things to create neutrinos.

The "one-zone" idea means the scientists assume everything happens in this one room. It's a simple, clean way to do the math.

The Problem: The "X-Ray Alarm"

The detective's job is to see if this "single room" can produce enough neutrinos to match what we detect on Earth. But there's a catch.

When the heavy protons smash into photons (light particles) to make neutrinos, they also create a messy side effect: a cascade of new particles that spews out X-rays.

The Analogy: Imagine you are trying to fill a bucket with water (neutrinos) from a hose. But every time you turn on the hose, it also sprays water everywhere else (X-rays).
The Constraint: We have very sensitive X-ray detectors in space. If the "spray" (X-rays) gets too strong, it would blow past what we actually observe. The universe is saying, "No, the X-rays can't be that bright."

Because of this strict X-ray limit, the "hose" (neutrino production) can't be turned up very high in this "single room" model.

The Investigation: Cracking the Code

The authors of this paper developed a clever mathematical shortcut. Instead of running thousands of complex computer simulations for every blazar, they created an analytical formula.

Think of it like a speed limit calculator. They asked: "Given the amount of X-rays we see (the speed limit), what is the absolute maximum amount of neutrinos this single room could possibly produce?"

They applied this shortcut to a list of known blazar candidates (the suspects) and then double-checked their math with full computer simulations to make sure the shortcut was accurate.

The Verdict: The Single Room Isn't Enough

The results were clear and consistent:

The Ceiling is Low: Even when they cranked the settings to the absolute maximum allowed by the X-ray limits, the "single room" model could only produce a tiny fraction of the neutrinos we actually see.
The Mismatch: The model predicts a "whisper" of neutrinos, but the IceCube telescope on Earth is hearing a "shout."
The Conclusion: The "one-room" theory is too simple. It's like trying to explain a massive concert by saying it's just one guy playing a guitar in a closet. It just doesn't add up.

What Does This Mean?

The paper concludes that if blazars are indeed the source of these high-energy neutrinos, the "single room" model is likely wrong. The action must be happening in a more complex way.

The authors suggest we need to look at multi-zone models.

The New Analogy: Instead of one room, imagine a factory with different departments. Maybe the neutrinos are being made in a hidden basement (the jet base) or a hot attic (the corona), while the X-rays we see are coming from the main floor. This separation would allow the neutrino "hose" to run full blast without triggering the X-ray alarm on the main floor.

Summary

In short, this paper used a new mathematical tool to test a popular theory about how blazars make neutrinos. They found that the theory is too restrictive; the "single room" model hits a wall because of X-ray observations. To explain the high-energy neutrinos we see, we likely need to imagine a more complex, multi-layered cosmic machine.

Technical Summary: On the maximum neutrino flux of blazars in the one-zone leptohadronic model

Problem Statement
The origin of extragalactic high-energy neutrinos remains a significant mystery in astrophysics, with blazars identified as leading candidate sources. While the one-zone leptohadronic jet model is the standard phenomenological framework for modeling blazar radiation, it faces severe challenges when applied to neutrino production. Specifically, the generation of high-energy neutrinos via photohadronic ( $p\gamma$ ) interactions is intrinsically coupled to the production of secondary electron-positron pairs via Bethe-Heitler (BH) processes and internal $\gamma\gamma$ pair production. These secondary cascades generate electromagnetic emission, particularly in the hard X-ray band. Consequently, stringent observational constraints on X-ray fluxes significantly limit the allowed neutrino flux within the one-zone framework, often rendering it insufficient to explain detected neutrino events.

Methodology
The authors develop a dual approach combining analytical derivation and numerical modeling to determine the maximum permissible neutrino flux for a sample of neutrino-candidate blazars within the one-zone leptohadronic framework.

Analytical Approach: The authors derive an analytical method to delineate the permissible physical parameter space (blob radius $R$ , magnetic field $B$ , and Doppler factor $\delta_D$ ) constrained by observed multi-wavelength spectral energy distributions (SEDs).
- They extend previous analytical methods to cover both Thomson (TMS) and Klein-Nishina (KN) regimes for Synchrotron Self-Compton (SSC) and External Compton (EC) dominated scenarios.
- The method utilizes peak frequencies and luminosities of the low- and high-energy SED bumps to constrain $R$ , $B$ , and $\delta_D$ .
- Crucially, the authors impose constraints based on the internal optical depth for $\gamma\gamma$ annihilation and the requirement that the BH-induced cascade emission does not exceed observed X-ray fluxes.
- By maximizing the ratio of $p\gamma$ interaction efficiency ( $f_{p\gamma}$ ) to BH process efficiency ( $f_{BH}$ ) while keeping the BH-induced X-ray emission at the observational upper limit, they derive an analytical upper limit for the neutrino flux.
Numerical Modeling: To validate the analytical results, the authors perform numerical simulations scanning the parameter space identified by the analytical approach.
- The model assumes primary relativistic electrons and protons injected with power-law distributions.
- It calculates steady-state distributions considering radiative cooling (synchrotron, IC, $p\gamma$ , BH, proton synchrotron) and particle escape.
- The code computes differential spectra for decayed $\gamma$ -rays, secondary pairs, and neutrinos, including contributions from up to three generations of secondary pairs from internal $\gamma\gamma$ annihilation.
- The proton injection luminosity is constrained primarily by keV-band observations rather than the Eddington limit.

Key Contributions

Analytical Upper Limits: The paper provides a rigorous analytical framework to derive the maximum neutrino flux as a function of observed X-ray flux and SED parameters, applicable to both SSC- and EC-dominated blazars in TMS and KN regimes.
Parameter Space Constraints: The study systematically maps the physical parameter space for a sample of eight neutrino-candidate blazars (including TXS 0506+056, PKS 0735+178, and FSRQs like GB6 J2113+1121), identifying the specific combinations of $R$ , $B$ , and $\delta_D$ that yield the maximum neutrino flux without violating X-ray constraints.
Validation: The authors demonstrate consistency between their analytical estimates and numerical modeling, confirming that the analytical method provides a reliable upper bound for the one-zone scenario.

Results

Flux Discrepancy: Both analytical and numerical results consistently show that the maximum neutrino flux predicted by the one-zone leptohadronic model does not significantly exceed values found in previous one-zone studies and remains below the flux levels inferred from IceCube observations for most sources.
Detection Rates: The predicted annual neutrino detection rates for the sample sources are generally low (e.g., $< 0.1$ for most cases), suggesting that the one-zone scenario alone is unlikely to fully account for the observed neutrino-blazar associations.
Physical Constraints: To achieve the maximum flux, the models often require proton injection luminosities that are highly super-Eddington (up to $10^5$ times the Eddington luminosity) for SSC-dominated cases, particularly for compact emission regions. EC-dominated solutions are physically more plausible, requiring lower luminosities (sub-Eddington to $\sim 27 \times$ Eddington), yet still fail to reach observed flux levels.
Spectral Shifts: The study notes that while optimized parameters can shift the peak of the secondary pair cascade away from the X-ray band (potentially correlating with $\gamma$ -ray flares), this does not resolve the fundamental issue of insufficient neutrino flux.

Significance and Claims
The paper concludes that the one-zone leptohadronic framework, while useful as a first-order description, has limited explanatory power for current neutrino-blazar associations due to the intrinsic coupling between neutrino production and electromagnetic cascade constraints. The authors claim that their analytical approach provides a practical tool to assess the maximum neutrino yield under the one-zone assumption, explicitly tied to X-ray flux measurements.

The primary significance of the work lies in its demonstration that the tension between neutrino production and X-ray constraints is robust within the one-zone paradigm. Consequently, the authors argue that explaining high-energy neutrino origins in blazars likely requires moving beyond this simplified picture to consider more complex scenarios, such as multi-zone jet models (e.g., inner-outer blob models) or alternative production sites like the jet base, hot corona, or accretion flow. The paper does not claim to solve the neutrino origin problem but rather defines the boundaries of the one-zone solution space, highlighting the necessity for more sophisticated models.