Multimessenger Constraints on Production Sites of High-Energy Neutrinos from NGC 1068

By analyzing IceCube and Fermi-LAT data, this study establishes multimessenger constraints on NGC 1068's high-energy neutrino production, demonstrating that hadronuclear ($pp$) processes allow for a larger emission region than photohadronic ($p\gamma$) models while challenging certain shock acceleration scenarios and supporting magnetically powered corona models with hard cosmic-ray spectra.

The Gist

Imagine the universe as a giant, chaotic kitchen. In the center of this kitchen sits a massive, hungry chef: a Supermassive Black Hole in a galaxy called NGC 1068. This chef is so powerful that it's been spitting out invisible particles called neutrinos (ghostly particles that can pass through anything) at incredibly high speeds.

For years, scientists have been trying to figure out how this chef is cooking up these neutrinos. Is it a simple recipe, or something much more complex?

This paper is like a team of food critics (the scientists) trying to reverse-engineer the recipe by looking at two things:

1. The Neutrinos: The "ghostly smoke" coming out of the kitchen.
2. The Gamma Rays: The "visible steam and sparks" that should be coming out if the recipe is standard.

Here is the breakdown of their investigation, using simple analogies:

1. The Mystery of the "Missing Sparks"

Usually, when you cook something that creates neutrinos, you also create gamma rays (a type of high-energy light). It's like baking a cake: if you have the batter (neutrinos), you should have the crumbs (gamma rays).

But with NGC 1068, we see a huge amount of neutrinos but very few gamma rays. It's as if the chef is baking a massive cake but somehow hiding all the crumbs. This suggests the "kitchen" (the emission region) is either very small and dense, or the ingredients are interacting in a way we haven't fully understood yet.

2. The Two Main Recipes Tested

The scientists tested two main ways the black hole might be cooking:

• Recipe A: The "Photon Buffet" (Photohadronic)
  Imagine a proton (a particle) running into a wall of light (photons) and smashing into it. This is the "standard" recipe.

  • The Result: For this to work without creating too many visible sparks (gamma rays), the kitchen has to be tiny. The scientists found the cooking area must be very close to the black hole, roughly 3 to 10 times the size of the black hole itself. It's like the chef is cooking in a tiny, cramped closet.

• Recipe B: The "Crowded Room" (Hadronuclear / pp)
  Imagine the proton running into other protons in a crowded room. This is the new focus of this paper.

  • The Result: This recipe is more forgiving. Because the particles are bumping into each other rather than light, the kitchen can be larger (up to 30–70 times the size of the black hole). It's like the chef has a bigger dining room to work in. This suggests that "crowded room" collisions might be a bigger part of the recipe than we thought.

3. The "Energy Budget" Problem

The scientists also looked at the "energy bill."

• The Problem: If the chef is using a "soft" recipe (where the ingredients are spread out over a wide range of energies, like a smoothie), the energy required to make the neutrinos is higher than the total energy the black hole is putting out as light (X-rays).
• The Analogy: It's like trying to power a massive stadium with a single AA battery. It doesn't add up.
• The Conclusion: The recipe must be "hard" (concentrated energy), meaning the chef is using a very specific, high-powered method (like magnetic reconnection) to cook, rather than a slow, shock-wave method.

4. The "Beta Decay" Red Herring

There was a third theory: maybe the neutrinos come from heavy atoms breaking apart (like a radioactive battery leaking).

• The Verdict: The scientists ran the numbers and said, "Nope." Even if they turned the magnetic field up to the absolute maximum allowed by physics (the "Eddington limit"), this recipe still produces too many visible sparks (gamma rays) and requires too much energy. It's a dead end.

The Final Takeaway

So, what is the chef actually doing?

1. It's likely a "Magnetic Corona": The cooking is happening in a super-hot, magnetized cloud right above the black hole.
2. It's crowded: The "crowded room" collisions (protons hitting protons) are likely helping to make the neutrinos, allowing the cooking area to be a bit bigger than previously thought.
3. It's efficient: The energy source is likely magnetic turbulence (like a blender spinning at high speed) rather than simple shockwaves.

In short: NGC 1068 is a cosmic kitchen where a supermassive black hole is using powerful magnetic fields to smash particles together in a crowded, compact space, creating a ghostly neutrino feast while keeping the visible sparks (gamma rays) surprisingly hidden. This paper helps us narrow down exactly how that magic trick is being performed.

Technical Summary

1. Problem Statement

The origin of high-energy astrophysical neutrinos remains a major open question in high-energy astrophysics. While the IceCube collaboration detected a significant neutrino flux from the nearby Seyfert II galaxy NGC 1068 (significance $\sim 4\sigma$), the source is "hidden" or opaque to GeV–TeV gamma rays. This implies that the dominant neutrino production mechanisms must occur in environments where gamma rays are either absorbed or suppressed.

Previous work (e.g., Das et al. 2024) focused primarily on the photohadronic ($p\gamma$) scenario, constraining the emission region to a compact, strongly magnetized corona ($R \sim 3-10 R_S$). However, the contribution of hadronuclear ($pp$) interactions and the beta decay of accelerated nuclei had not been fully explored in the context of NGC 1068's specific multimessenger constraints. Furthermore, the impact of extending the cosmic-ray (CR) injection spectrum down to GeV energies (rather than just TeV) on the required energetics and gamma-ray constraints needed re-evaluation.

2. Methodology

The authors utilized the Astrophysical Multimessenger Emission Simulator (AMES) to model neutrino and gamma-ray production in NGC 1068. Key methodological components include:

• Source Parameters:
  • Black Hole Mass ($M_{BH}$): $10^7 M_\odot$.
  • Distance ($d_L$): 10 Mpc.
  • Eddington Ratio ($\lambda_{Edd}$): $\sim 1$ (near-Eddington).
  • Intrinsic X-ray Luminosity ($L_{2-10 \text{ keV}}$): $\approx 3.4 \times 10^{43}$ erg s$^{-1}$.
• Physical Models:
  • Standard Hadronic Scenario: Neutrinos produced via $pp$ (proton-proton) and $p\gamma$ (proton-photon) interactions.
  • Beta Decay Scenario: Neutrinos produced via photodisintegration of accelerated nuclei (He) followed by neutron beta decay.
  • Target Photon Fields: X-rays from the hot corona, optical/UV from the accretion disk, and infrared from the dust torus.
  • Attenuation: Included gamma-ray absorption via $\gamma\gamma$ pair production (on disk, torus, CMB, and EBL) and Bethe-Heitler pair production.
• Spectral Assumptions:
  • Injected CR spectrum: Power-law with exponential cutoff ($\varepsilon_p^{2-s_{CR}} e^{-\varepsilon_p/\varepsilon_{max}}$).
  • Two Injection Regimes:
    1. Minimal Hadronic: $\varepsilon_{min} = 10$ TeV (mimicking hard spectra).
    2. Broad Spectrum: $\varepsilon_{min} = 1$ GeV (testing softer spectra and shock models).
  • Magnetic Field Parameterization: Defined by $\xi_B = U_B / U_\gamma$ (ratio of magnetic to radiation energy density) and $\beta$ (plasma beta).
• Constraints:
  • Neutrino Data: IceCube 1.5–15 TeV flux (best-fit and 95% contours).
  • Gamma-ray Data: Upper limits from Fermi-LAT and MAGIC.
  • Energetics: Cosmic-ray luminosity ($L_{CR}$) must not exceed the bolometric luminosity ($L_{bol}$) or total X-ray luminosity ($L_{corona}$) in standard models.

3. Key Contributions

• Extension to Hadronuclear ($pp$) Interactions: The paper systematically quantifies how $pp$ interactions relax constraints on the emission radius ($R$) compared to pure $p\gamma$ scenarios, particularly for high Thomson optical depths ($\tau_T$).
• Broad Spectrum Analysis: The authors investigate the consequences of injecting cosmic rays down to 1 GeV. This is critical for testing shock acceleration models which typically predict softer spectra ($s_{CR} > 2$).
• Explicit Beta Decay Constraints: The study explicitly tests the beta decay scenario under the most favorable conditions (maximum magnetic fields allowed by Eddington luminosity), providing a definitive exclusion of this mechanism for NGC 1068.
• Magnetization Constraints: The work derives lower limits on magnetic field strength ($\xi_B$) and upper limits on plasma beta ($\beta$) required to satisfy gamma-ray non-detection.

4. Key Results

A. Standard Hadronic Scenario ($pp + p\gamma$)

• Emission Radius ($R$):
  • Photohadronic ($p\gamma$) only: Favors compact regions, $R \lesssim 15 R_S$ (low $\beta$) and $R \lesssim 3 R_S$ (high $\beta$).
  • Hadronuclear ($pp$) included: The constraints are significantly relaxed. For low $\beta$ plasma ($\xi_B=1$), $R \lesssim 30-70 R_S$. For high $\beta$ plasma ($\xi_B=0.01$), $R \lesssim 5-50 R_S$.
  • Transition: The transition from $p\gamma$-dominated to $pp$-dominated neutrino production occurs at $R_{pp=p\gamma} \approx 33 R_S$ (for typical parameters).
• Cosmic-Ray Luminosity ($L_{CR}$) and Spectral Index ($s_{CR}$):
  • When the injection spectrum extends to 1 GeV, the required $L_{CR}$ increases drastically for softer spectra ($s_{CR} > 2$).
  • For $s_{CR} \gtrsim 2$, $L_{CR}$ exceeds the total X-ray luminosity ($L_{corona}$) and approaches the Eddington limit.
  • Implication: This strongly disfavors shock acceleration models that attribute coronal X-rays to shock dissipation (which predict $L_{CR} < L_{corona}$ and soft spectra).
• Magnetic Field Constraints:
  • To satisfy gamma-ray limits with $s_{CR} \gtrsim 2$, a strong magnetic field is required ($\xi_B \gtrsim 0.1$, implying $\beta \lesssim 10$).
  • For $s_{CR} \lesssim 2$ (hard spectra), $L_{CR} < L_{corona}$ is achievable, supporting magnetically powered corona models (turbulent/reconnection-driven) with $\beta \sim 1-3$.

B. Beta Decay Scenario

• Exclusion: The beta decay scenario is excluded as the dominant origin of NGC 1068 neutrinos.
• Reasoning: Even assuming the maximum possible magnetic field ($B_{max}$, limited by Eddington luminosity), the required cosmic-ray luminosity ($L_{CR}$) exceeds the Eddington limit ($L_{CR} \gg L_{Edd}$) by orders of magnitude.
• Gamma-Ray Constraints: The inverse Compton cascade from secondary pairs produces a gamma-ray flux that violates Fermi-LAT and MAGIC upper limits. Strong synchrotron cooling (at $B \approx B_{max}$) suppresses GeV flux but is insufficient to bring the total flux below observational limits.

5. Significance and Conclusion

• Support for Standard Hadronic Models: The results provide strong evidence for the standard hadronic scenario in NGC 1068, specifically favoring magnetically powered coronae with hard cosmic-ray spectra ($s_{CR} \lesssim 2$) and strong magnetic fields ($\beta \lesssim 10$).
• Refutation of Beta Decay: The study definitively rules out beta decay as the primary mechanism for NGC 1068, even under extreme magnetic field assumptions, resolving recent debates (e.g., Yasuda et al. 2025).
• Constraints on Shock Models: The analysis challenges simple shock acceleration models (failed winds, accretion shocks) that predict soft spectra and low magnetic fields, as these would require super-Eddington cosmic-ray powers or violate gamma-ray limits.
• Future Directions: The paper highlights the need for further investigation into $pp$ contributions and suggests that future multimessenger constraints should bridge phenomenological models with first-principle kinetic simulations (PIC/MHD) to better understand particle injection and acceleration in AGN coronae.

In summary, the paper establishes that NGC 1068 is a hidden neutrino source best explained by a compact, magnetically dominated corona where cosmic rays are accelerated to hard spectra, with hadronuclear interactions playing a crucial role in relaxing spatial constraints on the emission region.