Pith. sign in

REVIEW 3 major objections 5 minor 72 references

The X-ray states of the nearby galaxy NGC 4278 are explained by a radiatively inefficient accretion flow whose rate changes, while its TeV gamma rays likely escape from a jet or wind rather than the inner disk.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-02 18:37 UTC pith:XFSHCA6H

load-bearing objection First NuSTAR hard X-ray spectrum of NGC 4278 is a solid observational result, but the RIAF interpretation and the TeV escape conclusion lean on model-derived inputs that need to be flagged as such. the 3 major comments →

arxiv 2603.07029 v2 pith:XFSHCA6H submitted 2026-03-07 astro-ph.HE

Interpreting Swift and NuSTAR Observations of the Low-Luminosity Active Galactic Nucleus NGC 4278 with Radiatively Inefficient Accretion Flows and Implications for Neutrino Emission

classification astro-ph.HE
keywords NGC 4278low-luminosity AGNradiatively inefficient accretion flowX-ray variabilityTeV gamma raysneutrino emissionmagnetically arrested diskLHAASO
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper reports the first NuSTAR hard X-ray observations of the low-luminosity active galactic nucleus NGC 4278, detecting the source beyond 10 keV with a power-law spectrum and no sign of a high-energy cutoff. The authors argue that the X-ray emission, including a quiescent state, a moderate state, and the brighter 2021 active state, is produced by a radiatively inefficient accretion flow (RIAF) whose accretion rate varies between these states. They further show that within this RIAF, 1–25 TeV gamma rays would be absorbed through pair production on the model's infrared photon field, so the TeV source detected by LHAASO most likely originates in outer regions such as a jet or wind, not the inner disk. They also estimate the hidden PeV neutrino emission from the RIAF and propose that NGC 4278 fits a hard X-ray/neutrino luminosity correlation that extends from Seyfert galaxies down to low-luminosity AGNs. A sympathetic reader would care because this provides a coherent physical picture connecting X-ray variability, the origin of the TeV gamma rays, and the neutrino detectability of one of the closest low-luminosity AGNs.

Core claim

Using a single-zone radiatively inefficient accretion flow model with viscosity parameter α=0.4, plasma β=0.7, and an accretion rate that increases from 6×10⁻⁴ to 10⁻³ of the Eddington rate between the quiescent and moderate states, the paper reproduces the Swift-XRT and NuSTAR X-ray spectra of NGC 4278, including the hard power-law shape and the observed variability. The same model predicts a two-photon annihilation optical depth greater than unity for TeV gamma rays inside roughly 30 Schwarzschild radii, so the LHAASO-detected 1–25 TeV emission cannot emerge from the RIAF disk and must instead come from a jet or wind. The paper also computes RIAF neutrino spectra and argues that NGC 4278 s

What carries the argument

The central object is the radiatively inefficient accretion flow (RIAF) model: a hot, geometrically thick, optically thin inflow in which thermal electrons emit through synchrotron and Comptonization. The paper uses the analytic scalings of this model—density, electron temperature, magnetic field, and Thomson depth as functions of radius, viscosity α, plasma β, and accretion rate ṁ—to fit the X-ray spectral energy distribution. The pair-production optical depth formula for γγ→e⁺e⁻ (Eq. 7) is then used to determine the minimum radius from which TeV gamma rays can escape, providing the key constraint that pushes the TeV emission region outward.

Load-bearing premise

The conclusion that TeV gamma rays cannot escape the inner disk rests on the model-derived infrared luminosity of ∼3×10³⁹ erg/s at 0.1 eV; if the true IR field is more than about a factor of 20 weaker, the pair-production opacity falls below unity and 1–25 TeV photons could escape from within 30 Schwarzschild radii.

What would settle it

Measure the nuclear infrared luminosity of NGC 4278 at ∼0.1 eV (e.g., with ALMA or JWST); if it is below ∼1.5×10³⁸ erg/s, the optical depth in Eq. 7 becomes less than one in the LHAASO energy band, so the inner RIAF would no longer be ruled out as the TeV source.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The first NuSTAR detection above 10 keV establishes the hard X-ray spectrum of NGC 4278 and rules out an exponential cutoff below roughly 30 keV, constraining the electron energy distribution in the inner flow.
  • A single variable parameter—the accretion rate—can explain both the factor-of-two month-scale X-ray variability and the roughly fivefold higher X-ray flux during the 2021 active state, without invoking an additional soft X-ray component.
  • If the TeV gamma rays originate in a jet or wind at radii beyond ∼30 Schwarzschild radii, the RIAF's infrared photons naturally serve as target photons for external inverse-Compton scattering, making the jet interpretation self-consistent with the X-ray modeling.
  • The predicted RIAF neutrino flux is comparable to the observed TeV gamma-ray flux, providing a concrete target for neutrino telescopes and strengthening the case for a hard X-ray/neutrino luminosity correlation connecting LLAGNs and Seyferts.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the true sub-mm/far-IR luminosity of NGC 4278 is more than roughly a factor of 20 lower than the model's ∼3×10³⁹ erg/s at 0.1 eV, the pair-production opacity would drop below unity and 1–25 TeV photons could escape from the inner RIAF, removing the need for a jet/wind origin. This could be tested with ALMA or JWST measurements of the nuclear IR–sub-mm SED.
  • A direct measurement of the IR field would also sharpen the neutrino flux prediction, since the photohadronic and hadronuclear neutrino yields depend on the same target-photon population; a lower IR luminosity would shift the relative importance of pp versus pγ channels.
  • If future simultaneous X-ray and TeV observations catch a flare that tracks the accretion-rate timescale, the inner-disk versus jet/wind distinction could be settled decisively; a lag or lack of correlation would support the two-zone picture, while rapid correlated variability would favor a compact inner origin.
  • The same single-zone RIAF fitting strategy could be applied to other LLAGNs with sparse hard-X-ray coverage to identify which are hidden TeV sources and the most promising targets for IceCube and future neutrino telescopes; extending the L_X–L_ν relation to lower luminosities is a testable prediction for stacked LLAGN searches.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper presents the first NuSTAR hard X-ray observations of the LLAGN NGC 4278, clearly detecting the source beyond 10 keV with a power-law spectrum (photon index ~2.2–2.5) and no evidence for a high-energy cutoff. Two NuSTAR epochs separated by about a month show a factor-of-two flux change, and quasi-simultaneous Swift-XRT data define a moderate state; comparison with a 2021 Swift-XRT exposure during the LHAASO-reported active period shows a higher, harder X-ray state. The authors interpret the X-ray and archival broadband SED with a single-zone RIAF model (Kimura et al. 2019a, 2021), varying the accretion rate mdot from 6e-4 (quiescent) to 1e-3 (moderate) with fixed alpha=0.4, beta=0.7, and R=10 R_S. They argue that TeV gamma rays observed by LHAASO cannot escape from the inner RIAF because of two-photon pair-production opacity against IR disk photons, and therefore likely originate in an outer jet or wind. They further estimate RIAF neutrino fluxes and place NGC 4278 on a hypothetical hard-X-ray/neutrino correlation extending from Seyferts to LLAGNs.

Significance. The NuSTAR data are new and valuable: they provide the first hard X-ray spectrum of NGC 4278 above 10 keV, and the variability analysis is carefully done with standard HEASOFT/NuSTARDAS reduction and a fractional variability measurement. If the RIAF interpretation is robust, it would support hot accretion-flow models as the dominant hard X-ray source in LLAGNs and tie X-ray state changes to accretion-rate variations. The pair-production opacity argument, if secure, would locate the LHAASO TeV emission outside the inner disk, favoring jet/wind scenarios, and the neutrino discussion is a useful, clearly speculative extension. However, the central conclusions are currently conditional on model-dependent quantities and lack quantitative fit diagnostics, so the paper's strongest contribution is the data themselves, while the interpretive claims need further support or explicit hedging.

major comments (3)
  1. [Sec. 5.1, Eq. (7)] The TeV-opacity conclusion that the RIAF disk is transparent only for R >~ 30 R_S rests on the IR target luminosity L_disk ~ 3e39 erg/s at ~0.1 eV. As the paper itself notes in Sec. 4.2, the archival 0.1 eV IR excess is dominated by heated dust, so the nuclear RIAF IR component is a model output, not a measured quantity; it is degenerate with alpha, beta, and mdot and was rescaled by up to a factor of 4.1 to accommodate the lower NuSTAR flux. Since tau_gamma_gamma is proportional to L_IR, a nuclear IR luminosity lower by only ~20 would drop tau below unity at R=10 R_S, allowing 1-25 TeV photons to escape from R < 30 R_S. The claim that LHAASO gamma rays 'cannot escape from the RIAF disk' is therefore not empirically anchored. The authors should either provide an observational upper limit on the nuclear IR luminosity or explicitly reframe the conclusion as conditional with the scaling sho
  2. [Sec. 4.3 and Table 2] The best-fit RIAF parameters are quoted without uncertainties, and no goodness-of-fit statistic is reported. The model is explicitly acknowledged to have large degeneracy (Sec. 4.2), and the claimed fit to the broadband SED involves rescaling the optical flux by up to a factor of 4.1. Thus the statement 'We find the best-fit values' is not quantitatively supported. A grid-search with confidence regions, or at least a chi-square/delta-chi-square analysis over the parameter ranges stated in Sec. 4.3, is needed to justify the RIAF interpretation and the inferred mdot values. As written, the fit cannot be distinguished from a visual overlay.
  3. [Sec. 4.3] The inference that the flow is magnetically arrested (MAD) is based mainly on the chosen beta=0.7, which is described as typical of MAD-like disks, while the paper also notes that beta alone is not enough to determine the magnetic topology. This claim is used later to support consistency with jet models, but the supporting evidence is thin. The authors should soften this or provide additional diagnostics (e.g., variability properties or jet-power comparison) before asserting that the model 'favors' a MAD.
minor comments (5)
  1. [Abstract/Introduction] Typographical issues: 'lager viewing angle' should be 'larger viewing angle' (Sec. 1), and 'comsic-ray' should be 'cosmic-ray' (Sec. 5.2).
  2. [Sec. 3.1] The phrase 'almost flat X-ray spectrum' is slightly misleading for a photon index of 2.2-2.5 in E F_E space; consider saying 'hard power-law spectrum' or explicitly defining 'flat' in nu F_nu space.
  3. [Eq. (7)] The notation with tildes (e.g., L_tilde_disk, epsilon_tilde) is not defined in the text. Define these quantities explicitly before use.
  4. [Fig. 2] The top panel's shaded 'LHAASO-quiet period' and 'LHAASO-active period' regions are helpful but the boundary definitions should be stated in the caption or text for reproducibility.
  5. [Sec. 5.2] The discussion of the L_nu-L_X correlation is clearly labeled as speculative, but the figure caption for Fig. 9 should also explicitly repeat the caveat that no neutrino detection from NGC 4278 exists; currently this is only in the text.

Circularity Check

0 steps flagged

No significant circularity: the RIAF model is an external framework, and model-dependent neutrino/opacity statements are clearly conditional, not identity-level reductions.

full rationale

The paper's central X-ray interpretation is a fit of a previously published multi-purpose RIAF model (Kimura et al. 2019a, 2021) to new Swift/NuSTAR spectra; the model is not derived from the NGC 4278 data, and the fit parameters (α=0.4, β=0.7, ṁ=6×10⁻⁴/10⁻³) are presented as best-fit values, not as predictions forced by a prior step. The TeV-opacity argument (Sec. 5.1, Eq. 7) uses a model-derived IR luminosity (L̃_disk ~3×10³⁹ erg/s at 0.1 eV) rather than a direct measurement; this makes the gamma-ray escape conclusion model-dependent, but it is not circular because L̃_disk is an output of the SED model, not a quantity tuned to the LHAASO detection, and the paper explicitly notes the attenuation could be compensated by increasing dissipation power. The neutrino flux estimate and the Lν–LX point (Sec. 5.2, Fig. 9) are computed from the same ṁ fitted to the X-ray data, so they are not an independent empirical confirmation; however, the authors label this a 'rough estimation' and explicitly caveat that no neutrino emission from NGC 4278 has ever been reported. This is model extrapolation, not an equation-level identity or a fitted parameter renamed as a prediction. Self-citations to the authors' RIAF/neutrino methodology are present, but they cite a published framework with stated assumptions; no uniqueness theorem or unverified ansatz is imported under disguise. No step in the derived chain reduces, by the paper's own equations, to its own input, so no circular step is established.

Axiom & Free-Parameter Ledger

8 free parameters · 6 axioms · 0 invented entities

The ledger shows that the paper's conclusions rest on a previously published RIAF framework (largely by the same authors) plus several hand-set parameters. The X-ray data are externally anchored, so the central claim is not definitionally circular, but the neutrino flux and the Lnu-LX placement are derived from the same mdot fitted to X-ray data, and the gamma-ray escape argument depends on a model-derived IR luminosity.

free parameters (8)
  • Viscosity parameter alpha = 0.4
    Chosen as best-fit to the multi-wavelength SED; no uncertainty quoted.
  • Plasma beta (gas-to-magnetic pressure ratio) = 0.7
    Best-fit value; used to argue for a MAD-like disk, although the text admits beta alone cannot decide MAD vs SANE.
  • Normalized accretion rate mdot (quiescent) = 6e-4
    Set so the single-zone RIAF SED matches the December 2024 NuSTAR flux.
  • Normalized accretion rate mdot (moderate) = 1e-3
    Increased to fit the January 2025 NuSTAR/Swift flux; the 'variability' explanation is this fitted change, not an independent prediction.
  • Multi-zone radial index s = 0.5
    Hand-chosen from literature (Yuan et al. 2012; Guo et al. 2024) in Sec. 4.1; affects multi-zone emission profile, though outer zones are subdominant.
  • Cosmic-ray injection efficiency eta_CR = 0.01
    Assumed in Sec. 5.2 for the neutrino spectra; neutrino luminosity scales directly with this number.
  • Acceleration efficiency eta_acc = 10 (Model A), 1e4 (Model B)
    Two ad hoc scenarios in Sec. 5.2 bracket the neutrino predictions; not constrained by data in this paper.
  • Emission radius R = 10 R_S (single-zone)
    Chosen as the typical emitting region; the gamma-ray escape conclusion depends on where RIAF photons are produced.
axioms (6)
  • domain assumption The RIAF model of Kimura et al. (2019a, 2021) correctly describes the accretion flow structure and radiative processes in NGC 4278.
    Invoked in Sec. 4.1; the entire SED interpretation rests on this prior model, which is not re-derived here.
  • domain assumption Electron heating fraction f_e follows the Chael et al. (2018) fitting formula; half of released gravitational energy heats the plasma.
    Used in Sec. 4.1 to set electron temperature; if f_e is different, the SED changes.
  • domain assumption Black hole mass M_BH = 3e8 M_sun (Wang & Zhang 2003).
    Used throughout for R_S, Eddington luminosity, and opacity scaling.
  • ad hoc to paper The quiescent and moderate/active states differ only by the accretion rate mdot, with alpha and beta fixed.
    This is the paper's variability interpretation (Sec. 4.2, Table 2), not derived from independent data.
  • ad hoc to paper Cosmic rays are injected with a power-law spectral index 1 and do not affect the thermal electron spectrum except through a subdominant cascade.
    Assumed in Sec. 5.2 to compute neutrino spectra; not constrained by X-ray data.
  • domain assumption The IR photon field for pair-production opacity is dominated by RIAF disk emission with L_disk ~ 3e39 erg/s at ~0.1 eV and Gamma_disk ~ 1.
    Used in Eq. (7) to conclude tau>1 and R > 30 R_S; if the IR luminosity were lower, TeV gamma rays could escape from the inner disk.

pith-pipeline@v1.3.0-alltime-deepseek · 16303 in / 15683 out tokens · 129765 ms · 2026-08-02T18:37:27.449855+00:00 · methodology

0 comments
read the original abstract

We report the first NuSTAR hard X-ray observations of the low-luminosity active galactic nucleus NGC 4278. The source is clearly detected beyond 10 keV with a hard X-ray spectrum consistent with a power law of photon index between 2.2 and 2.5 without evidence for a high-energy cutoff. The X-ray flux is low compared to the active state in 2021, but exhibits variability by a factor of ~2 on a timescale of a month. We discuss the origin of the hard X-ray emission and explore its connection to gamma rays and high-energy neutrinos. We explain the X-ray data, including both quiescent and active states, using a radiatively inefficient accretion flow (RIAF) model with a variable accretion rate. We also show that TeV gamma rays cannot escape from the RIAF disk, and very high-energy gamma rays observed in LHAASO are likely to originate from outer regions such as jets and winds, which is consistent with our results favoring a magnetically arrested disk. We also discuss hidden neutrino emission from RIAFs together with possible connections to coronae of active galactic nuclei with standard, radiatively efficient disks.

Figures

Figures reproduced from arXiv: 2603.07029 by Abhishek Das, Ashwani Pandey, Eleanor Young, Kohta Murase, Qi Feng, Shigeo S. Kimura.

Figure 1
Figure 1. Figure 1: The X-ray spectra of NGC 4278. Quasi-simul￾taneous spectra measured by Swift-XRT (blue squares) and NuSTAR (dark blue down triangles for FPMA and light blue plus signs for FPMB) on January 13, 2025, show a moder￾ate flux level. The Swift-XRT spectrum measured during the LHAASO-reported active state (pink circles) shows an X-ray active state that is more prominent above 2 keV. The NuSTAR spectrum (orange up… view at source ↗
Figure 2
Figure 2. Figure 2: (Top) The X-ray light curve of NGC 4278 be￾tween 2020 and 2025, showing an active X-ray flux state during the LHAASO-reported active periods. (Bottom) A zoomed-in view of the variability during the X-ray campaign in 2024 and 2025, showing a low to moderate flux state with variability by a factor of ∼ 2 over timescales of weeks to months. The 10 recent Swift-XRT observations between May 2024 and January 202… view at source ↗
Figure 3
Figure 3. Figure 3: shows the broadband SED of NGC 4278 during the quiescent and moderate X-ray flux states in 2024 and 2025, together with the soft X-ray SED during the active state in 2021, as well as the archival SEDs across many wavelengths. The TeV gamma-ray spectra mea￾sured by LHAASO (Cao et al. 2024b) during the ac￾tive state (April 10, 2021 – August 28, 2021) and the quiet state (March 5, 2021 – October 31, 2023, exc… view at source ↗
Figure 4
Figure 4. Figure 4: Left: SEDs with varying α for constant values of β and ˙m. Middle: Same with β for constant values of α and ˙m. Right: Same with ˙m for constant values of α and β. The archival and X-ray data points have the same legends as fig. 3. The purple and black data points are from Nemmen et al. (2014) and Bandyopadhyay et al. (2019) respectively. The bolometric luminosity of the SED is computed with given physical… view at source ↗
Figure 5
Figure 5. Figure 5: Multi-zone SED from different radii with R = (10 − 40) RS, compared to a contribution from each radius. The archival and X-ray data points have the same legends as [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Best fit results for a single-zone model using a typical emission radius of 10 RS. The best fit parameters for the SED are listed in table 2. The archival and X-ray data points have the same legends as [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Two-photon annihilation optical depths for the single-zone moderate state SED at different emission radii. The red shaded region is the LHAASO energy band. 5. DISCUSSIONS 5.1. Incapability of the escape of very high-energy gamma rays from the RIAF disk It is interesting to ask whether the RIAF model ex￾plains very high-energy gamma rays from the RIAF. For the RIAF model to explain the LHAASO data, the sour… view at source ↗
Figure 8
Figure 8. Figure 8: Neutrino spectra from our RIAF model. The solid line is the photon spectrum emitted by thermal elec￾trons. The dotted and dashed lines are neutrino spectra by pp and pγ interactions, respectively. The blue and red lines are for model A (ηacc = 10) and model B (ηacc = 1.0 × 104 ), respectively (see text for the definition of ηacc). gamma rays from the Comptonization process of ther￾mal electrons (Kimura et … view at source ↗
Figure 9
Figure 9. Figure 9: (a) shows a possible scaling relation between the intrinsic 2-10 keV hard X-ray luminosity and the neu￾trino luminosity from Seyfert galaxies. The estimation of neutrino luminosities is sensitive to the templates of neutrino spectra, and the power-law assumption is not valid in general, or at least can be misleading for the model comparison. Thus, for the neutrino luminosi￾ties of NGC 1068 and NGC 4151, we… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

72 extracted references · 7 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    =6kls|"=46 - 9U|ׁt<y0g

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

  4. [4]

    P., Tollerud , E

    Astropy Collaboration , Robitaille , T. P., Tollerud , E. J., et al. 2013, , 558, A33, 10.1051/0004-6361/201322068

  5. [5]

    M., Sip o cz , B

    Astropy Collaboration , Price-Whelan , A. M., Sip o cz , B. M., et al. 2018, , 156, 123, 10.3847/1538-3881/aabc4f

  6. [6]

    A., & Hawley, J

    Balbus, S. A., & Hawley, J. F. 1998, Rev. Mod. Phys., 70, 1, 10.1103/RevModPhys.70.1

  7. [7]

    2021, The Astrophysical Journal Supplement Series, 256, 13

    Baldini, L., Ballet, J., Bastieri, D., et al. 2021, The Astrophysical Journal Supplement Series, 256, 13

  8. [8]

    M., et al

    Bandyopadhyay, B., Xie, F.-G., Nagar, N. M., et al. 2019, Monthly Notices of the Royal Astronomical Society, 490, 4606, 10.1093/mnras/stz2874

  9. [9]

    D., & Begelman, M

    Blandford, R. D., & Begelman, M. C. 1999, Monthly Notices of the Royal Astronomical Society, 303, L1, 10.1046/j.1365-8711.1999.02358.x

  10. [10]

    J., Tr ""u mper , J., et al

    Boller , T., Freyberg , M. J., Tr ""u mper , J., et al. 2016, aap, 588, A103, 10.1051/0004-6361/201525648

  11. [11]

    2024, , 977, L16, 10.3847/2041-8213/ad93cf

    Bronzini , E., Grandi , P., Torresi , E., & Buson , S. 2024, , 977, L16, 10.3847/2041-8213/ad93cf

  12. [12]

    N., Hill , J

    Burrows , D. N., Hill , J. E., Nousek , J. A., et al. 2005, , 120, 165, 10.1007/s11214-005-5097-2

  13. [13]

    2024 a , , 271, 25, 10.3847/1538-4365/acfd29

    Cao , Z., Aharonian , F., An , Q., et al. 2024 a , , 271, 25, 10.3847/1538-4365/acfd29

  14. [14]

    2024 b , , 971, L45, 10.3847/2041-8213/ad5e6d

    Cao , Z., Aharonian , F., Axikegu , et al. 2024 b , , 971, L45, 10.3847/2041-8213/ad5e6d

  15. [15]

    A., Kheirandish, A., & Murase, K

    Carpio, J. A., Kheirandish, A., & Murase, K. 2026, to be submitted

  16. [16]

    E., Narayan, R., Johnson, M

    Chael, A., Rowan, M. E., Narayan, R., Johnson, M. D., & Sironi, L. 2018, Mon. Not. Roy. Astron. Soc., 478, 5209, 10.1093/mnras/sty1261

  17. [17]

    T., et al

    Chen, S., Das, A., Zhang, B. T., et al. 2026, Physical origin of very-high-energy gamma rays from the low-luminosity active galactic nucleus NGC 4278 and implications for neutrino observations. 2601.23242

  18. [18]

    J., Cotton , W

    Condon , J. J., Cotton , W. D., Greisen , E. W., et al. 1998, , 115, 1693, 10.1086/300337

  19. [19]

    T., & Murase, K

    Das, A., Zhang, B. T., & Murase, K. 2024, The Astrophysical Journal, 972, 44, 10.3847/1538-4357/ad5a04

  20. [20]

    D., Murase, K., & Inoue, Y

    Dermer, C. D., Murase, K., & Inoue, Y. 2014, JHEAp, 3-4, 29, 10.1016/j.jheap.2014.09.001

  21. [21]

    Dixon , R. S. 1970, , 20, 1, 10.1086/190216

  22. [22]

    2004, , 611, 1005, 10.1086/422091

    Gehrels , N., Chincarini , G., Giommi , P., et al. 2004, , 611, 1005, 10.1086/422091

  23. [23]

    B., & Giovannini , G

    Giroletti , M., Taylor , G. B., & Giovannini , G. 2005, , 622, 178, 10.1086/427898

  24. [24]

    C., Scott , W

    Gregory , P. C., Scott , W. K., Douglas , K., & Condon , J. J. 1996, , 103, 427, 10.1086/192282

  25. [25]

    M., Quataert, E., & Kim, C.-G

    Guo, M., Stone, J. M., Quataert, E., & Kim, C.-G. 2024, Astrophys. J., 973, 141, 10.3847/1538-4357/ad5fe7

  26. [26]

    H. E. S. S. Collaboration , Aharonian , F., Ait Benkhali , F., et al. 2022, Science, 376, 77, 10.1126/science.abn0567

  27. [27]

    A., Craig, W

    Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, , 770, 103, 10.1088/0004-637x/770/2/103

  28. [28]

    E., Romani , R

    Healey , S. E., Romani , R. W., Taylor , G. B., et al. 2007, , 171, 61, 10.1086/513742

  29. [29]

    Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90, 10.1109/MCSE.2007.55

  30. [30]

    A., Browne , I

    Jackson , N., Battye , R. A., Browne , I. W. A., et al. 2007, , 376, 371, 10.1111/j.1365-2966.2007.11442.x

  31. [31]

    Joint Iras Science , W. G. 1994, VizieR Online Data Catalog, 2125, 0

  32. [32]

    2001, SciPy : Open source scientific tools for Python

    Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy : Open source scientific tools for Python . http://www.scipy.org/

  33. [33]

    Kheirandish, A., Murase, K., & Kimura, S. S. 2021, Astrophys. J., 922, 45, 10.3847/1538-4357/ac1c77

  34. [34]

    S., Kashiyama, K., & Hotokezaka, K

    Kimura, S. S., Kashiyama, K., & Hotokezaka, K. 2021, Astrophys. J. Lett., 922, L15, 10.3847/2041-8213/ac35dc

  35. [35]

    S., Murase, K., & M\'esz\'aros, P

    Kimura, S. S., Murase, K., & M\'esz\'aros, P. 2019 a , Phys. Rev. D, 100, 083014, 10.1103/PhysRevD.100.083014

  36. [36]

    S., Murase , K., & M \'e sz \'a ros , P

    Kimura , S. S., Murase , K., & M \'e sz \'a ros , P. 2021, Nature Communications, 12, 5615, 10.1038/s41467-021-25111-7

  37. [37]

    S., Murase, K., & Toma, K

    Kimura, S. S., Murase, K., & Toma, K. 2015, Astrophys.J., 806, 159, 10.1088/0004-637X/806/2/159

  38. [38]

    S., & Toma, K

    Kimura, S. S., & Toma, K. 2020, Astrophys. J., 905, 178, 10.3847/1538-4357/abc343

  39. [39]

    S., Tomida, K., & Murase, K

    Kimura, S. S., Tomida, K., & Murase, K. 2019 b , Mon. Not. Roy. Astron. Soc., 485, 163, 10.1093/mnras/stz329

  40. [40]

    2024, arXiv e-prints, arXiv:2404.06867, 10.48550/arXiv.2404.06867

    Kun , E., Bartos , I., Becker Tjus , J., et al. 2024, arXiv e-prints, arXiv:2404.06867, 10.48550/arXiv.2404.06867

  41. [41]

    1997, Astrophys

    Mahadevan, R. 1997, Astrophys. J., 477, 585, 10.1086/303727

  42. [42]

    1990, in IRAS Faint Source Catalogue, version 2.0 (1990), 0

    Moshir , M., & et al. 1990, in IRAS Faint Source Catalogue, version 2.0 (1990), 0

  43. [43]

    2022, , 941, L17, 10.3847/2041-8213/aca53c

    Murase , K. 2022, , 941, L17, 10.3847/2041-8213/aca53c

  44. [44]

    Murase, K., Inoue, Y., & Dermer, C. D. 2014, Phys. Rev. D, 90, 023007, 10.1103/PhysRevD.90.023007

  45. [45]

    M., Kimura, S

    Murase, K., Karwin, C. M., Kimura, S. S., Ajello, M., & Buson, S. 2024, Astrophys. J. Lett., 961, L34, 10.3847/2041-8213/ad19c5

  46. [46]

    S., & Meszaros, P

    Murase, K., Kimura, S. S., & Meszaros, P. 2020, Phys. Rev. Lett., 125, 011101, 10.1103/PhysRevLett.125.011101

  47. [47]

    S., Mukhopadhyay, M., & Bhattacharya, M

    Murase, K., Kimura, S. S., Mukhopadhyay, M., & Bhattacharya, M. 2026. 2602.20145

  48. [49]

    V., & Abramowicz, M

    Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, Publications of the Astronomical Society of Japan, 55, L69, 10.1093/pasj/55.6.L69

  49. [50]

    F., & Kulkarni, A

    Narayan, R., Sadowski, A., Penna, R. F., & Kulkarni, A. K. 2012, Monthly Notices of the Royal Astronomical Society, 426, 3241, 10.1111/j.1365-2966.2012.22002.x

  50. [51]

    1994, , 428, L13, 10.1086/187381

    Narayan , R., & Yi , I. 1994, , 428, L13, 10.1086/187381

  51. [52]

    2014, HEAsoft: Unified Release of FTOOLS and XANADU , Astrophysics Source Code Library, record ascl:1408.004

    Nasa High Energy Astrophysics Science Archive Research Center (Heasarc) . 2014, HEAsoft: Unified Release of FTOOLS and XANADU , Astrophysics Source Code Library, record ascl:1408.004. 1408.004

  52. [53]

    S., Storchi-Bergmann, T., & Eracleous, M

    Nemmen, R. S., Storchi-Bergmann, T., & Eracleous, M. 2014, Monthly Notices of the Royal Astronomical Society, 438, 2804, 10.1093/mnras/stt2388

  53. [54]

    S., Storchi-Bergmann, T., Eracleous, M., & Yuan, F

    Nemmen, R. S., Storchi-Bergmann, T., Eracleous, M., & Yuan, F. 2009, Proceedings of the International Astronomical Union, 5, 313–318, 10.1017/S1743921310006538

  54. [55]

    2012, , 758, 94, 10.1088/0004-637X/758/2/94

    Pellegrini , S., Wang , J., Fabbiano , G., et al. 2012, , 758, 94, 10.1088/0004-637X/758/2/94

  55. [56]

    W., Palumbo, D

    Pesce, D. W., Palumbo, D. C. M., Narayan, R., et al. 2021, Astrophys. J., 923, 260, 10.3847/1538-4357/ac2eb5

  56. [57]

    2011, arXiv e-prints, arXiv:1103.0749, 10.48550/arXiv.1103.0749

    Stratta , G., Capalbi , M., Giommi , P., et al. 2011, arXiv e-prints, arXiv:1103.0749, 10.48550/arXiv.1103.0749

  57. [58]

    1987, Mon

    Svensson, R. 1987, Mon. Not. Roy. Astron. Soc., 227, 403

  58. [59]

    2011, Monthly Notices of the Royal Astronomical Society, 414, 1827, 10.1111/j.1365-2966.2011.18267.x

    Tang, Y., Gu, Q., Zhang, S., & Tang, B. 2011, Monthly Notices of the Royal Astronomical Society, 414, 1827, 10.1111/j.1365-2966.2011.18267.x

  59. [60]

    Tchekhovskoy , A., Narayan , R., & McKinney , J. C. 2011, , 418, L79, 10.1111/j.1745-3933.2011.01147.x

  60. [61]

    E., Taylor , G

    Tremblay , S. E., Taylor , G. B., Ortiz , A. A., et al. 2016, , 459, 820, 10.1093/mnras/stw592

  61. [62]

    C., & Varoquaux, G

    van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Computing in Science Engineering, 13, 22, 10.1109/MCSE.2011.37

  62. [63]

    S., & Uttley , P

    Vaughan , S., Edelson , R., Warwick , R. S., & Uttley , P. 2003, , 345, 1271, 10.1046/j.1365-2966.2003.07042.x

  63. [64]

    1999, , 349, 389

    Voges , W., Aschenbach , B., Boller , T., et al. 1999, , 349, 389

  64. [65]

    2003, , 340, 793, 10.1046/j.1365-8711.2003.06336.x

    Wang , T.-G., & Zhang , X.-G. 2003, , 340, 793, 10.1046/j.1365-8711.2003.06336.x

  65. [66]

    L., & Becker , R

    White , R. L., & Becker , R. H. 1992, , 79, 331, 10.1086/191656

  66. [67]

    L., Becker , R

    White , R. L., Becker , R. H., Helfand , D. J., & Gregg , M. D. 1997, , 475, 479

  67. [68]

    Willingale , R., Starling , R. L. C., Beardmore , A. P., Tanvir , N. R., & O'Brien , P. T. 2013, , 431, 394, 10.1093/mnras/stt175

  68. [69]

    L., Eisenhardt , P

    Wright , E. L., Eisenhardt , P. R. M., Mainzer , A. K., et al. 2010, , 140, 1868, 10.1088/0004-6256/140/6/1868

  69. [70]

    , Porquet, D

    Younes, G. , Porquet, D. , Sabra, B. , et al. 2010, , 517, A33, 10.1051/0004-6361/201014371

  70. [71]

    2026, TeV Gamma-Rays from the Low-Luminosity Active Galactic Nucleus NGC 4278: Implications for the Diffuse Neutrino Background

    Yuan, C., & Liu, R.-Y. 2026, TeV Gamma-Rays from the Low-Luminosity Active Galactic Nucleus NGC 4278: Implications for the Diffuse Neutrino Background. 2601.21411

  71. [72]

    2003, Astrophys

    Yuan, F., Quataert, E., & Narayan, R. 2003, Astrophys. J., 598, 301, 10.1086/378716

  72. [73]

    2012, The Astrophysical Journal, 761, 129, 10.1088/0004-637X/761/2/129

    Yuan, F., Wu, M., & Bu, D. 2012, The Astrophysical Journal, 761, 129, 10.1088/0004-637X/761/2/129