Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

Probing a cosmogenic origin of astrophysical neutrinos and cosmic rays using gamma-ray observations of TXS 0506+056

A. Acharyya, A. Archer, P. Bangale, J. T. Bartkoske, W. Benbow, J. H. Buckley, Y. Chen, J. L. Christiansen, A. Duerr, M. Errando, M. Escobar Godoy, A. Falcone, S. Feldman, Q. Feng, S. Filbert, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Holder, Z. Hughes, M. Iskakova, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda, T. K. Kleiner, N. Korzoun, M. J. Lang, M. Lundy, G. Maier, C. L. Mooney, R. Mukherjee, W. Ning, R. A. Ong, A. Pandey, M. Pohl, E. Pueschel, J. Quinn, P. L. Rabinowitz, K. Ragan, P. T. Reynolds, D. Ribeiro, E. Roache, I. Sadeh, A. C. Sadun, L. Saha, G. H. Sembroski, R. Shang, M. Splettstoesser, D. Tak, A. K. Talluri, J. V. Tucci, J. Valverde, D. A. Williams, S. L. Wong, T. Yoshikoshi, M. Meyer, J. Müller

TL;DR

This work tests whether cosmogenic proton cascades along the line of sight from TXS 0506+056 can shape the observed gamma-ray spectrum associated with IceCube-170922A. Using CRPropa-based simulations, Fermi-LAT and VERITAS data are fit with a composite model of a primary gamma-ray spectrum plus a cascade component, while exploring proton injection parameters and magnetic-field assumptions. The LP×EBL primary spectrum provides a better overall fit than PLE×EBL, with best-fit proton parameters around $α_{p} \approx 2.0$ and $E_{p,\text{max}} \approx 1.3\times10^{16}$ eV, and 90% CL limits on the proton-escape luminosity $L_{p,\text{esc}}$ in the range $1\times10^{44}$–$3\times10^{45}$ erg s$^{-1}$. Neutrino predictions from the same proton population imply that, for certain luminosities, cosmogenic neutrinos would conflict with IceCube constraints, yielding an upper bound on the escape fraction $f_{esc}$ of order $10^{-6}$–$10^{-7}$. Overall, the results favor a consistent but small role for intergalactic cascades in TXS 0506+056 and highlight the potential of joint gamma-ray–neutrino analyses to probe cosmic-ray propagation and proton acceleration in blazars.

Abstract

In September 2017, a high-energy neutrino event detected by the IceCube Neutrino Observatory (IceCube-170922A) was associated, at the $3σ$ level, with a gamma-ray flare from the blazar TXS 0506+056. Cosmic rays that are accelerated in astrophysical sources can escape from their jets and interact with background radiation fields. Interactions with the extragalactic background light can produce pions and hence neutrinos, while interactions with the cosmic microwave background predominantly drive inverse Compton scattering, contributing to electromagnetic cascades in intergalactic space. The resulting secondary gamma-ray emission can be detected with high-energy gamma-ray telescopes. Here, we report on a new search for such cosmogenic cascade emission from the blazar TXS 0506+056, using a combined data set from the Fermi-Large Area Telescope and VERITAS. We compare the gamma-ray spectrum and neutrino observations with the predictions of cosmic-ray induced cascades in intergalactic space. The observed gamma-ray spectrum is modeled as a combination of the primary spectrum and the cascade spectrum. We apply a Monte Carlo simulation with a $Δχ^2$-based likelihood analysis to jointly determine the best-fit parameters of a proton emission spectrum describing the data and derive constraints on the proton escape luminosity. Assuming a log-parabola primary photon spectrum, we find consistency with a proton injection spectral index of $α_{p} \simeq 2.0$ and a cutoff energy of $E_{p,\text{max}} \simeq 1.3 \times 10^{16}$ eV, and constrain the isotropic proton escape luminosity to $1 \times 10^{44}$ erg s$^{-1}$ $\lesssim L_{p, esc} \lesssim 3 \times 10^{45}$ erg s$^{-1}$ at the 90 % confidence level.

Probing a cosmogenic origin of astrophysical neutrinos and cosmic rays using gamma-ray observations of TXS 0506+056

TL;DR

This work tests whether cosmogenic proton cascades along the line of sight from TXS 0506+056 can shape the observed gamma-ray spectrum associated with IceCube-170922A. Using CRPropa-based simulations, Fermi-LAT and VERITAS data are fit with a composite model of a primary gamma-ray spectrum plus a cascade component, while exploring proton injection parameters and magnetic-field assumptions. The LP×EBL primary spectrum provides a better overall fit than PLE×EBL, with best-fit proton parameters around and eV, and 90% CL limits on the proton-escape luminosity in the range – erg s. Neutrino predictions from the same proton population imply that, for certain luminosities, cosmogenic neutrinos would conflict with IceCube constraints, yielding an upper bound on the escape fraction of order –. Overall, the results favor a consistent but small role for intergalactic cascades in TXS 0506+056 and highlight the potential of joint gamma-ray–neutrino analyses to probe cosmic-ray propagation and proton acceleration in blazars.

Abstract

In September 2017, a high-energy neutrino event detected by the IceCube Neutrino Observatory (IceCube-170922A) was associated, at the level, with a gamma-ray flare from the blazar TXS 0506+056. Cosmic rays that are accelerated in astrophysical sources can escape from their jets and interact with background radiation fields. Interactions with the extragalactic background light can produce pions and hence neutrinos, while interactions with the cosmic microwave background predominantly drive inverse Compton scattering, contributing to electromagnetic cascades in intergalactic space. The resulting secondary gamma-ray emission can be detected with high-energy gamma-ray telescopes. Here, we report on a new search for such cosmogenic cascade emission from the blazar TXS 0506+056, using a combined data set from the Fermi-Large Area Telescope and VERITAS. We compare the gamma-ray spectrum and neutrino observations with the predictions of cosmic-ray induced cascades in intergalactic space. The observed gamma-ray spectrum is modeled as a combination of the primary spectrum and the cascade spectrum. We apply a Monte Carlo simulation with a -based likelihood analysis to jointly determine the best-fit parameters of a proton emission spectrum describing the data and derive constraints on the proton escape luminosity. Assuming a log-parabola primary photon spectrum, we find consistency with a proton injection spectral index of and a cutoff energy of eV, and constrain the isotropic proton escape luminosity to erg s erg s at the 90 % confidence level.

Paper Structure

This paper contains 13 sections, 4 equations, 5 figures, 1 table.

Table of Contents

  1. Introduction
  2. Observations and Data Analysis
  3. VERITAS
  4. Fermi-LAT
  5. Swift
  6. The Role of proton Cascade Emission
  7. Results
  8. Constraints on the spectrum of protons escaping the source
  9. Log parabola primary photon spectrum
  10. Power law with exponential cut-off primary photon spectrum
  11. Simulated neutrino spectrum
  12. Constraints on the fraction of escaped protons
  13. Summary

Figures (5)

  • Figure 1: Multiwavelength lightcurves of TXS 0506+056 between MJD 54683 and MJD 60218, which corresponds to midnight on August 4, 2008 until midnight on October 1, 2023. The blue narrow band corresponds to the IceCube alert, IceCube-170922A, observed on September 22, 2017 (MJD 58018). The orange shaded region corresponds to the time interval between September 23, 2017 and February 6, 2018, investigated in VERITAS_TXS0506_original. (A) Monthly-binned VHE lightcurve for VERITAS observations above an energy threshold of 190 GeV. Upper limits at 95% confidence level are quoted for observations with a significance of lower than $2\sigma$. (B)Fermi-LAT lightcurve in 30 day bins in the 0.1--300 GeV band. (C) Swift-XRT energy flux in the 0.3–10 keV band. (D and E) Optical/UV flux observations with the Swift-UVOT in the UW1, UW2, UM2, U, B, and V band filters, split into the UV and optical components for enhanced readability.
  • Figure 2: The simulated photon cascade flux per solid angle for the best-fit spectrum i.e the subdominant cascade, assuming a LP $\times$ EBL primary photon spectrum (top) and PLE $\times$ EBL primary photon spectrum (bottom). The values assumed for the maximum blazar activity time and magnetic field are $t_{\text{max}} = 6$ months and $B = 10^{-16} \text{G}$ respectively. The pink dotted line indicates the 90$\%$ confidence-level region for the neutrino event IceCube-170922A (from 2018_Icecube). The green, cyan, and purple dashed lines represent the PSFs of the VERITAS and the Fermi-LAT observations and the VERITAS OFF region, respectively.
  • Figure 3: The observed Fermi-LAT and VERITAS spectra of TXS 0506+056 (red and magenta markers, respectively), together with the total best-fit spectrum (dashed blue line, including uncertainties). This is the sum of the primary photon spectra LP $\times$ EBL (top) and PLE $\times$ EBL (bottom) (dashed green lines, including uncertainties) and the cascade components (dashed orange lines, including uncertainties) respectively. Spectral points are shown as upper limits if the detection significance is below 2$\sigma$. The SEDs of the predicted cosmogenic neutrinos is shown in black. The spectral curves for the proton cascades and the neutrino spectrum correspond to the best-fit proton injection spectral parameters when considering a LP $\times$ EBL primary photon spectrum and PLE $\times$ EBL primary photon spectrum respectively (see Table. \ref{['tab:table']}). The values assumed for the maximum blazar activity time and magnetic field are $t_{\text{max}} = 6$ months and $B = 10^{-16} \text{G}$ respectively. Also shown are the representative neutrino flux upper limits that produce on average one detection like IceCube-170922A over a period of 0.5 and 7.5 years assuming an injection spectrum of $dN/dE \propto E^{-2}$ at the most probable neutrino energy, $E_{\nu} = 311$ TeV 2018_Icecube.
  • Figure 4: The stats-surfaces obtained over the investigated parameter space of the proton spectrum, $E_{p, \text{max}}$ and $\alpha_{p}$, assuming a LP $\times$ EBL primary photon spectrum (top) and PLE $\times$ EBL primary photon spectrum (bottom). The values assumed for the maximum blazar activity time and magnetic field are $t_{\text{max}} = 6$ months and $B = 10^{-16} \text{G}$ respectively. The color scale denotes the $\Delta \chi^{2}$, the difference between the $\chi^{2}$ values at each point and the minimum $\chi^{2}$ value obtained over the entire parameter space. The blue stars mark the best-fit parameter combination corresponding to the global minimum $\chi^{2}$, and listed in Table \ref{['tab:table']}. The dashed purple curves show 2$\sigma$ uncertainty contours. The dashed green curves represent points in the parameter space which correspond to the same proton escape luminosity, as labeled. The dark green curves mark the lower and upper bounds on the proton escape luminosity, denoted as $L_{\min}$ and $L_{\max}$, respectively and shown in the legend.
  • Figure 5: The $\Delta \chi^{2}$ values obtained as a function of the proton escape luminosity $L_{p, esc}$. The high-energy cut-off is fixed at $E_{p, \text{max}} = 3.6 \times 10^{17}$ eV and the proton spectral index are fixed in the interval $\alpha_{p} \in [1.8, 2.6]$ at intervals of 0.1.