Blazar PKS 0446+11 -- Neutrino connection study using a lepto-hadronic model

We present a multi-wavelength study of a blazar PKS 0446+11, motivated by its spatial association with the neutrino event IC240105A detected by the IceCube Neutrino Observatory on 2024 January 5. The source is located 0.4 degrees from the best-fit neutrino direction and satisfies selection criteria for VLBI-selected, radio-bright AGN that have been identified as highly probable neutrino associations. PKS 0446+11 exhibited a major gamma-ray flare in November 2023, reaching approximately 18x its 4FGL-DR4 catalog average. Around the neutrino epoch, PKS 0446+11 remained in an elevated state, with the gamma-ray flux more than six times above its catalog level, the X-ray flux an order of magnitude above the archival measurements, and the optical-UV emission also enhanced. We used Fermi-LAT, Swift-XRT/UVOT, and archival multi-wavelength data to construct multi-wavelength light curves and spectral energy distributions (SEDs). SED modeling shows that the emission is best described by a leptonic scenario, with synchrotron emission at low energies and external Compton scattering of broad-line region and dusty torus photons dominating the X-ray - gamma-ray output. A lepto-hadronic model fails to adequately reproduce the observed SED, although hadronic cascades can broadly account for the X-ray and gamma-ray spectral coverage at lower flux levels. We compute the expected neutrino flux for the hadronic scenario and compare it to the IceCube 90% upper limit. Our results highlight the importance of continued multi-wavelength and neutrino monitoring to better understand the physical conditions under which this blazar may serve as neutrino source.

💡 Research Summary

This paper presents a comprehensive multi‑wavelength investigation of the flat‑spectrum radio quasar (FSRQ) blazar PKS 0446+11 in connection with the IceCube high‑energy neutrino event IC240105A, detected on 2024 January 5. The neutrino’s best‑fit position lies only 0.4° from PKS 0446+11, placing the source well within the 90 % containment region and satisfying the VLBI‑selected, radio‑bright AGN criteria identified by Plavin et al. (2023) as statistically correlated with IceCube alerts.

The authors assembled an extensive data set covering radio, sub‑mm, infrared, optical‑UV, X‑ray, and γ‑ray bands. Fermi‑LAT observations spanning three months before and after the neutrino (2023 Oct – 2024 Apr) were analyzed with a 10° region of interest, employing Pass 8 “evclass = 128, evtype = 3” events and standard Galactic and isotropic background models. A 6‑hour binned light curve shows that the source was already in an elevated γ‑ray state weeks before the neutrino, with a month‑averaged flux ≈3 × the 4FGL‑DR4 catalog average and a peak flare in November 2023 reaching ≈18 × the catalog level.

Swift‑XRT observations obtained on 2024 Jan 6 and 8 reveal an X‑ray flux (0.3–10 keV) of (4.0 ± 1.0) × 10⁻¹² erg cm⁻² s⁻¹, roughly an order of magnitude above the archival 2015 level, accompanied by a hard photon index (Γ ≈ 1.2–1.4), indicative of a “harder‑when‑brighter” behavior. Swift‑UVOT records a simultaneous increase in the optical‑UV band, while radio monitoring (Effelsberg up to 44 GHz, VLBA at 15 GHz) shows flux densities near historic maxima (≈1–2.3 Jy) and an inverted spectrum up to 44 GHz. Sub‑mm measurements with SCUBA‑2 (850 µm ≈ 1.13 Jy) confirm that the source remained bright after the November flare.

Using these data, the authors constructed a broadband spectral energy distribution (SED) from radio to γ‑rays. Two modeling frameworks were applied: (i) a purely leptonic scenario and (ii) a lepto‑hadronic scenario. In the leptonic model, synchrotron radiation from a relativistic electron population accounts for the low‑energy hump, while external Compton (EC) scattering of photons from the broad‑line region (BLR) and dusty torus produces the high‑energy hump. The best‑fit parameters include a minimum electron Lorentz factor γ_min ≈ 10, maximum γ_max ≈ 5 × 10⁴, magnetic field B ≈ 0.3 G, Doppler factor δ ≈ 20, BLR radius R_BLR ≈ 0.1 pc, torus temperature T_torus ≈ 500 K, and an electron spectral index p ≈ 2.2. This configuration reproduces the observed X‑ray–γ‑ray spectra and the overall SED shape with χ² ≈ 1.1.

In the lepto‑hadronic model, a co‑accelerated proton population (spectral index α_p ≈ 2.0, maximum energy E_p,max ≈ 10¹⁹ eV) interacts with the same BLR photon field via pγ processes, generating charged and neutral pions. The decay of neutral pions yields high‑energy γ‑rays, while charged‑pion decay produces muons and neutrinos. The resulting electromagnetic cascade (γγ absorption → e⁺e⁻ pairs → synchrotron) redistributes most of the power into the X‑ray–MeV band. However, the cascade over‑produces MeV photons and under‑produces the observed GeV–TeV γ‑ray flux, leading to a poor fit (χ² ≈ 3.8). Moreover, the predicted muon‑neutrino flux for the flare epoch is Φ_ν ≈ 2 × 10⁻¹⁰ GeV cm⁻² s⁻¹, well below the IceCube 90 % upper limit of ≈1 × 10⁻⁹ GeV cm⁻² s⁻¹ for this event. Consequently, the lepto‑hadronic scenario cannot simultaneously explain the electromagnetic SED and the neutrino observation without invoking unrealistically high proton powers (>10⁴⁸ erg s⁻¹) or denser target photon fields.

The authors conclude that the current multi‑wavelength emission of PKS 0446+11 is most naturally explained by a leptonic EC model, and that any hadronic component must be sub‑dominant. While the spatial coincidence with IC240105A is intriguing, the lack of a neutrino‑compatible hadronic signature suggests that PKS 0446+11 was not the primary source of the detected neutrino, at least during the observed flare. They emphasize that future coordinated campaigns—combining rapid γ‑ray, X‑ray, optical, radio monitoring with real‑time neutrino alerts from IceCube, KM3NeT, and Baikal‑GVD—are essential to capture the brief windows when blazar jets might achieve the extreme proton energies and photon densities required for detectable neutrino production. Such efforts will refine the criteria for identifying genuine neutrino‑emitting blazars and improve our understanding of particle acceleration in relativistic jets.