Self-Consistent Modelling of Neutrino Production in Turbulent Black Hole Coronae

1. Problem Statement

The IceCube Collaboration has identified a significant excess of high-energy neutrinos (1–10 TeV) spatially associated with the nearby Seyfert 2 galaxy NGC 1068. This discovery suggests that non-jetted Active Galactic Nuclei (AGN) are significant contributors to the extragalactic neutrino flux. However, a major theoretical challenge exists:

• The Multi-Messenger Mismatch: NGC 1068 emits a strong neutrino flux but very weak gamma-ray flux (orders of magnitude below the neutrino flux). This implies the neutrino production region is optically thick to high-energy photons (due to $\gamma\gamma$ absorption), likely a compact, dense corona around the supermassive black hole (SMBH).
• Modeling Limitations: Previous models often relied on steady-state assumptions, semi-analytical spectra, or prescribed radiative losses. They frequently failed to self-consistently couple stochastic particle acceleration with the dynamical evolution of turbulence and radiative/hadronic cascades. Specifically, the feedback of accelerated particles on the turbulence (damping) and the time-dependent nature of the acceleration process were often neglected, despite being crucial in dense, magnetized environments.
• Scale Discrepancy: Direct kinetic simulations (Particle-in-Cell or PIC) cannot bridge the gap between microscopic acceleration scales and the macroscopic source size (requiring $\sim$5 orders of magnitude in energy), while standard transport codes often lack the necessary microphysical details.

2. Methodology: The Turb-AM3 Framework

The authors introduce Turb-AM3, a novel hybrid numerical code designed to solve these issues self-consistently.

• Core Architecture: The code couples the state-of-the-art time-dependent lepto-hadronic radiative solver AM3 with a custom stochastic acceleration module.
• Physical Environment:
  • Models a spherical, magnetized, turbulent corona around an SMBH (mass $M_{BH} \sim 10^{7.2} M_\odot$ for NGC 1068).
  • Assumes a hot, sub-relativistic electron population (Maxwellian) and a turbulent magnetic field.
  • Includes a soft photon field from the accretion disk (multi-color blackbody) which serves as the target for Inverse Compton scattering and photohadronic ($p\gamma$) interactions.
• Particle Acceleration:
  • Protons are accelerated via stochastic interactions with turbulent electric fields.
  • The code implements two transport operators: a standard Fokker-Planck diffusive scheme and a Generalized Fermi operator (based on Lemoine 2022) to handle large-amplitude turbulence and non-local momentum jumps.
  • Injection: Protons are injected from the thermal pool via magnetic reconnection in current sheets, with an injection fraction $\xi_p$ treated as a free parameter.
• Self-Consistent Feedback (Key Innovation):
  • The code solves the turbulent cascade equation simultaneously with the particle transport equation.
  • It accounts for turbulent damping: As protons extract energy from the turbulence to accelerate, the turbulent power spectrum is damped. This reduces the acceleration rate for subsequent particles, creating a self-regulating mechanism.
  • This ensures energy conservation between the non-thermal proton population and the turbulent magnetic energy.
• Radiative & Hadronic Processes:
  • Tracks protons, electrons/positrons, photons, neutrinos, pions, and muons.
  • Includes all relevant losses: Synchrotron, Inverse Compton, $pp$ collisions, $p\gamma$ interactions, and Bethe-Heitler pair production.
  • Models the electromagnetic cascade: High-energy gamma rays produced by hadronic interactions are absorbed via $\gamma\gamma \to e^+e^-$, generating secondary pairs that cool and re-emit, effectively reprocessing the energy to lower frequencies.
• Spatial Transport:
  • Uses a "one-zone" approximation where spatial transport is mapped to time evolution via an advection velocity ($v_{adv}$), representing the inflow of plasma toward the black hole or outflow in a jet.

3. Key Contributions

1. Development of Turb-AM3: A fully time-dependent, self-consistent numerical framework that bridges the gap between microphysical PIC simulations and macroscopic source modeling.
2. Inclusion of Turbulent Damping: A rigorous treatment of the back-reaction of accelerated protons on the turbulence cascade, demonstrating that this feedback naturally limits the non-thermal proton energy density to a level comparable to the turbulent magnetic pressure.
3. Unified Multi-Messenger Modeling: The ability to simultaneously predict proton, photon, and neutrino spectra while satisfying observational constraints (e.g., Fermi-LAT and MAGIC gamma-ray upper limits) and physical energy budgets.
4. Template Spectra: The derivation of generic "template" neutrino spectral shapes for turbulent AGN coronae, which can be used to interpret current and future neutrino data.

4. Key Results (NGC 1068 Case Study)

Using a fiducial set of parameters ($R_{cor} \approx 15 r_g$, $v_A \approx 0.25c$, $\xi_p = 0.1$, $\tau_T \approx 0.5$), the model yields the following:

• Neutrino Flux Reproduction: The model successfully reproduces the IceCube neutrino flux in the 1–10 TeV range. The spectral shape is characterized by a broad peak around a few TeV.
• Spectral Shape Mechanism:
  • Low Energy ($<1$ TeV): Dominated by $pp$ interactions, resulting in a relatively flat spectral extension ($E^2 \Phi \propto E^0$).
  • High Energy ($>1$ TeV): Dominated by $p\gamma$ interactions with coronal X-ray photons.
  • Cutoff: A sharp cutoff occurs at $\sim$30 TeV due to the dominance of photohadronic losses over acceleration.
• Gamma-Ray Consistency: The model predicts that high-energy gamma rays are efficiently absorbed and reprocessed within the corona. The resulting reprocessed emission (peaking in the MeV–GeV range) remains well below the observed gamma-ray upper limits from Fermi-LAT and MAGIC, resolving the "neutrino-gamma mismatch."
• Self-Regulation: The non-thermal proton energy density self-regulates to be comparable to the turbulent magnetic energy density ($\sim 10^{43}$ erg s$^{-1}$), consistent with energy budget constraints and Eddington limits.
• Parameter Sensitivity:
  • Injection Fraction ($\xi_p$): For $\xi_p \gtrsim 10^{-2}$, the spectral shape becomes universal and insensitive to the exact value of $\xi_p$ due to strong turbulent damping. Low $\xi_p$ leads to peaked spectra (linear regime).
  • Timescale Ratios: The ratio of escape to acceleration time ($t_{esc}/t_{acc}$) must be $>1$ for efficient high-energy acceleration. The ratio of advection to acceleration ($t_{adv}/t_{acc}$) sets the maximum energy cutoff.
• Extension to Other AGN: The model successfully reproduces neutrino data for other Seyfert galaxies (NGC 4151, NGC 7469) by simply scaling the disk luminosity (and thus the target photon density), explaining the observed trend where higher luminosity sources emit neutrinos at lower energies due to stronger cooling.

5. Significance

• Validation of Turbulent Acceleration: The results strongly support the hypothesis that stochastic acceleration in magnetized, turbulent coronae is the dominant mechanism for high-energy neutrino production in non-jetted AGN.
• New Diagnostic Tool: The specific spectral shape (flat low-energy extension vs. peaked spectrum) serves as a diagnostic for the pair-loading factor ($n_p/n_e$) and the strength of the turbulent feedback. Future high-sensitivity neutrino observations could distinguish between different coronal microphysics scenarios.
• Multi-Messenger Constraints: The framework demonstrates that AGN coronae can be "hidden" neutrino sources where the gamma-ray signal is suppressed by opacity, yet the neutrino signal remains detectable.
• Future Applications: The Turb-AM3 framework is applicable to a broad range of compact astrophysical environments, including X-ray binaries, tidal disruption events, and hidden neutrino sources, providing a robust tool for interpreting the multi-messenger sky.

In conclusion, this paper provides a rigorous, self-consistent physical explanation for the NGC 1068 neutrino signal, linking the observed flux and spectral shape directly to the microphysics of turbulence and the macroscopic conditions of the black hole corona.