skip to main content

Title: Soft gamma rays from low accreting supermassive black holes and connection to energetic neutrinos
Abstract

The Universe is filled with a diffuse background of MeV gamma-rays and PeV neutrinos, whose origins are unknown. Here, we propose a scenario that can account for both backgrounds simultaneously. Low-luminosity active galactic nuclei have hot accretion flows where thermal electrons naturally emit soft gamma rays via Comptonization of their synchrotron photons. Protons there can be accelerated via turbulence or reconnection, producing high-energy neutrinos via hadronic interactions. We demonstrate that our model can reproduce the gamma-ray and neutrino data. Combined with a contribution by hot coronae in luminous active galactic nuclei, these accretion flows can explain the keV – MeV photon and TeV – PeV neutrino backgrounds. This scenario can account for the MeV background without non-thermal electrons, suggesting a higher transition energy from the thermal to nonthermal Universe than expected. Our model is consistent with X-ray data of nearby objects, and testable by future MeV gamma-ray and high-energy neutrino detectors.

Authors:
; ;
Award ID(s):
1908689
Publication Date:
NSF-PAR ID:
10306129
Journal Name:
Nature Communications
Volume:
12
Issue:
1
ISSN:
2041-1723
Publisher:
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Recent multimessenger studies have provided evidence for high-energy neutrino sources that are opaque to GeV–TeV gamma rays. We present model-independent studies on the connection between neutrinos and gamma rays in the active galaxy NGC 1068, and find that the neutrinos most likely come from regions within ∼30–100 Schwarzschild radii. This is especially the case if neutrinos are produced via the photomeson production process, although the constraints could be alleviated if hadronuclear interactions are dominant. We consider the most favorable neutrino production regions, and discuss coronae, jets, winds, and their interactions with dense material. The results strengthen the importance of understanding dissipation mechanisms near the coronal region and the outflow base. There could be a connection between active galactic nuclei with near-Eddington accretion and tidal disruptions events, in that neutrinos are produced in the obscured vicinity of supermassive black holes.

  2. Abstract

    We investigate acceleration and propagation processes of high-energy particles inside hot accretion flows. The magnetorotational instability (MRI) creates turbulence inside accretion flows, which triggers magnetic reconnection and may produce non-thermal particles. They can be further accelerated stochastically by the turbulence. To probe the properties of such relativistic particles, we perform magnetohydrodynamic simulations to obtain the turbulent fields generated by the MRI, and calculate orbits of the high-energy particles using snapshot data of the MRI turbulence. We find that the particle acceleration is described by a diffusion phenomenon in energy space with a diffusion coefficient of the hard-sphere type: Dε ∝ ε2, where ε is the particle energy. Eddies in the largest scale of the turbulence play a dominant role in the acceleration process. On the other hand, the stochastic behaviour in configuration space is not usual diffusion but superdiffusion: the radial displacement increases with time faster than that in the normal diffusion. Also, the magnetic field configuration in the hot accretion flow creates outward bulk motion of high-energy particles. This bulk motion is more effective than the diffusive motion for higher energy particles. Our results imply that typical active galactic nuclei that host hot accretion flows can accelerate CRs up tomore »ε ∼ 0.1−10 PeV.

    « less
  3. Abstract

    We present a bottom-up calculation of the flux of ultrahigh-energy cosmic rays (UHECRs) and high-energy neutrinos produced by powerful jets of active galactic nuclei (AGNs). By propagating test particles in 3D relativistic magnetohydrodynamic jet simulations, including a Monte Carlo treatment of sub-grid pitch-angle scattering and attenuation losses due to realistic photon fields, we study the spectrum and composition of the accelerated UHECRs and estimate the amount of neutrinos produced in such sources. We find that UHECRs may not be significantly affected by photodisintegration in AGN jets where theespressomechanism efficiently accelerates particles, consistent with Auger’s results that favor a heavy composition at the highest energies. Moreover, we present estimates andupper boundsfor the flux of high-energy neutrinos expected from AGN jets. In particular, we find that (i) source neutrinos may account for a sizable fraction, or even dominate, the expected flux of cosmogenic neutrinos; (ii) neutrinos from theβ-decay of secondary neutrons produced in nucleus photodisintegration end up in the teraelectronvolt to petaelectronvolt band observed by IceCube, but can hardly account for the observed flux; (iii) UHECRs accelerated via theespressomechanism lead to nearly isotropic neutrino emission, which suggests that nearby radio galaxies may be more promising as potential sources. We discuss ourmore »results in light of multimessenger astronomy and current/future neutrino experiments.

    « less
  4. Abstract Galaxy clusters are considered to be gigantic reservoirs of cosmic rays (CRs). Some of the clusters are found with extended radio emission, which provides evidence for the existence of magnetic fields and CR electrons in the intra-cluster medium. The mechanism of radio halo (RH) emission is still under debate, and it has been believed that turbulent reacceleration plays an important role. In this paper, we study the reacceleration of CR protons and electrons in detail by numerically solving the Fokker–Planck equation, and show how radio and gamma-ray observations can be used to constrain CR distributions and resulting high-energy emission for the Coma cluster. We take into account the radial diffusion of CRs and follow the time evolution of their one-dimensional distribution, by which we investigate the radial profile of the CR injection that is consistent with the observed RH surface brightness. We find that the required injection profile is nontrivial, depending on whether CR electrons have a primary or secondary origin. Although the secondary CR electron scenario predicts larger gamma-ray and neutrino fluxes, it is in tension with the observed RH spectrum for hard injection indexes, α < 2.45. This tension is relaxed if the turbulent diffusion of CRsmore »is much less efficient than the fiducial model, or the reacceleration is more efficient for lower-energy CRs. In both the secondary and primary scenario, we find that galaxy clusters can make a sizable contribution to the all-sky neutrino intensity if the CR energy spectrum is nearly flat.« less
  5. Abstract

    Flat-spectrum radio quasars (FSRQs) are the most luminous blazars at GeV energies but only rarely emit detectable fluxes of TeV gamma rays, typically during bright GeV flares. We explore the gamma-ray variability and spectral characteristics of three FSRQs that have been observed at GeV and TeV energies by Fermi-LAT and VERITAS, making use of almost 100 hr of VERITAS observations spread over 10 yr: 3C 279, PKS 1222+216, and Ton 599. We explain the GeV flux distributions of the sources in terms of a model derived from a stochastic differential equation describing fluctuations in the magnetic field in the accretion disk and estimate the timescales of magnetic flux accumulation and stochastic instabilities in their accretion disks. We identify distinct flares using a procedure based on Bayesian blocks and analyze their daily and subdaily variability and gamma-ray energy spectra. Using observations from VERITAS, as well as Fermi, Swift, and the Steward Observatory, we model the broadband spectral energy distributions of PKS 1222+216 and Ton 599 during very high energy (VHE)–detected flares in 2014 and 2017, respectively, strongly constraining the jet Doppler factors and gamma-ray emission region locations during these events. Finally, we place theoretical constraints on the potential production ofmore »PeV-scale neutrinos during these VHE flares.

    « less