pith. sign in

arxiv: 2401.05534 · v2 · submitted 2024-01-10 · 🌌 astro-ph.HE

A closer look at the electromagnetic signatures of Bethe-Heitler pair production process in blazars

Pith reviewed 2026-05-24 04:30 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords Bethe-Heitler pair productionblazarsgamma-ray emissionsynchrotron radiationjet energeticspair production spectrum
0
0 comments X

The pith

Gamma-rays in low- and intermediate-peaked blazars may arise from synchrotron emission by Bethe-Heitler pairs in jets of size ~10^15 cm with magnetic fields 5-500 G.

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

The Bethe-Heitler process creates electron-positron pairs from photon-nucleus interactions and transfers energy from relativistic protons to pairs whose maximum energy follows from process kinematics rather than radiative losses. The paper supplies an empirical function that approximates the resulting pair spectrum for monoenergetic protons and photons, replacing a double numerical integration over the differential cross section. The authors then examine blazar jet conditions and conclude that synchrotron radiation from these pairs can produce the observed gamma rays in low- and intermediate-peaked sources once broadband spectra, jet energetics, and ambient radiation fields are taken into account.

Core claim

After deriving an empirical kernel for the Bethe-Heitler pair spectrum and considering broadband spectra, jet energetics, and radiation fields, the paper concludes that gamma-rays in low- and intermediate-peaked blazars may arise from Bethe-Heitler pairs in jet regions with typical transverse size ∼10^{15} cm and co-moving magnetic field 5-500 G.

What carries the argument

The empirical function describing the spectrum of Bethe-Heitler pairs from single-energy proton-photon interactions, which acts as a kernel replacing double numerical integration over the differential cross section.

Load-bearing premise

That the broadband spectral characteristics, jet energetics, and properties of radiation fields in the blazar environment allow the synchrotron emission from Bethe-Heitler pairs to dominate or be distinguishable from inverse Compton emission.

What would settle it

A measured gamma-ray spectrum in a low-peaked blazar that cannot be reproduced by synchrotron emission from the empirical Bethe-Heitler pair distribution but requires inverse Compton scattering, or a direct measurement showing the emission-region magnetic field lies outside 5-500 G.

read the original abstract

The ''twin birth'' of a positron and an electron by a photon in the presence of a nucleus, known as Bethe-Heitler pair production, is a key process in astroparticle physics. The Bethe-Heitler process offers a way of channeling energy stored in a population of relativistic protons (or nuclei) to relativistic pairs with extended distributions. Contrary to accelerated leptons, whose maximum energy is limited by radiative losses, the maximal energy of pairs is determined by the kinematics of the process and can be as high as the parent proton energy. We take a closer look at the features of the injected pair distribution, and provide a novel empirical function that describes the spectrum of pairs produced by interactions of single-energy protons with single-energy photons. The function is the kernel of the Bethe-Heitler pair production spectrum that replaces a double numerical integration involving the complex differential cross section of the process, and can be easily implemented in numerical codes. We further examine under which conditions Bethe-Heitler pairs produced in blazar jets can emit $\gamma$-ray photons via synchrotron radiation, thus providing an alternative to the inverse Compton scattering process for high-energy emission in jetted active galactic nuclei. After taking into consideration the broadband spectral characteristics of the source, the jet energetics, and the properties of radiation fields present in the blazar environment, we conclude that $\gamma$-rays in low- and intermediate-peaked blazars may arise from Bethe-Heitler pairs in regions of the jet with typical transverse size $\sim 10^{15}$ cm and co-moving magnetic field $5-500$ G.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript introduces a novel empirical function approximating the injected electron-positron pair spectrum from the Bethe-Heitler process for monoenergetic protons interacting with monoenergetic photons, intended to replace double numerical integration over the differential cross section for use in numerical codes. It then examines conditions under which synchrotron radiation from these pairs can produce the observed γ-ray emission in low- and intermediate-peaked blazars, concluding that this is possible in jet regions with transverse size ∼10^{15} cm and co-moving magnetic field 5–500 G after considering broadband spectra, jet energetics, and ambient radiation fields.

Significance. If the empirical kernel proves accurate and the stated jet parameters permit the Bethe-Heitler-pair synchrotron component to be both luminous enough and spectrally separable from inverse-Compton emission, the work supplies a practical computational tool for astroparticle modeling and offers a concrete alternative channel for high-energy emission in a subclass of blazars, with potential consequences for jet-power budgets and proton-acceleration constraints.

major comments (2)
  1. [Abstract] Abstract (final paragraph) and the section applying the kernel to blazars: the conclusion that synchrotron from Bethe-Heitler pairs can account for γ-rays rests on the assertion that broadband spectral characteristics, jet energetics, and radiation-field properties permit the component to dominate or be distinguishable from inverse Compton; however, no quantitative luminosity ratios, spectral-shape comparisons, or viability thresholds are supplied to demonstrate this is realized rather than merely possible in principle under the quoted R ∼ 10^{15} cm and B = 5–500 G window.
  2. [Section on empirical function] Section presenting the empirical function: the claim that the function is a ready-to-use kernel replacing the double integration is load-bearing for the entire analysis, yet the manuscript provides neither the fitting procedure, direct validation against the exact Bethe-Heitler cross section, nor an error characterization or uncertainty envelope on the approximation.
minor comments (1)
  1. Notation for the empirical kernel should be introduced with an explicit functional form and range of validity (proton and photon energies) to allow immediate implementation by readers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback. We address each major comment below and will incorporate revisions to strengthen the manuscript.

read point-by-point responses
  1. Referee: [Abstract] Abstract (final paragraph) and the section applying the kernel to blazars: the conclusion that synchrotron from Bethe-Heitler pairs can account for γ-rays rests on the assertion that broadband spectral characteristics, jet energetics, and radiation-field properties permit the component to dominate or be distinguishable from inverse Compton; however, no quantitative luminosity ratios, spectral-shape comparisons, or viability thresholds are supplied to demonstrate this is realized rather than merely possible in principle under the quoted R ∼ 10^{15} cm and B = 5–500 G window.

    Authors: We agree that explicit quantitative demonstrations would better support the conclusion. The revised manuscript will include calculations of luminosity ratios between the Bethe-Heitler synchrotron and inverse-Compton components for representative parameter values, along with direct spectral-shape comparisons and viability thresholds within the R ∼ 10^{15} cm and B = 5–500 G range. revision: yes

  2. Referee: [Section on empirical function] Section presenting the empirical function: the claim that the function is a ready-to-use kernel replacing the double integration is load-bearing for the entire analysis, yet the manuscript provides neither the fitting procedure, direct validation against the exact Bethe-Heitler cross section, nor an error characterization or uncertainty envelope on the approximation.

    Authors: We acknowledge that the fitting procedure, validation, and error analysis were not sufficiently detailed. The revised version will describe the fitting method used to derive the empirical kernel, present direct comparisons to numerical evaluations of the Bethe-Heitler differential cross section, and include quantitative error characterization with uncertainty estimates. revision: yes

Circularity Check

0 steps flagged

No circularity: empirical kernel is an independent approximation to the cross-section integral; blazar conclusion rests on external jet parameters

full rationale

The paper derives an empirical fitting function for the Bethe-Heitler pair spectrum as a replacement for the double integral over the known differential cross section; this is a numerical convenience whose validity can be checked against the underlying QED formula without reference to blazar data. The subsequent claim that synchrotron from these pairs can produce gamma-rays in low/intermediate-peaked blazars is conditioned on independently stated jet size, magnetic field, and radiation-field properties rather than on any quantity fitted or defined from the same dataset. No self-citation is invoked as a uniqueness theorem, and no prediction reduces by construction to a fitted input. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The empirical function is constructed by fitting to the exact Bethe-Heitler kernel, introducing at least one free functional form and associated parameters whose values are not reported in the abstract. The viability conclusion further depends on unstated assumptions about the ambient photon fields and jet geometry that are treated as given inputs from the broader blazar literature.

pith-pipeline@v0.9.0 · 5836 in / 1211 out tokens · 18068 ms · 2026-05-24T04:30:33.796044+00:00 · methodology

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

64 extracted references · 64 canonical work pages · 23 internal anchors

  1. [1]

    Bethe and W

    H. Bethe and W. Heitler, On the stopping of fast particles and on the creation of positive electrons, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 146 (1934) 83

  2. [2]

    Penzias and R.W

    A.A. Penzias and R.W. Wilson, A Measurement of Excess Antenna Temperature at 4080 Mc/s., The Astrophysical Journal 142 (1965) 419

  3. [3]

    Linsley, Evidence for a primary cosmic-ray particle with energy 1020 ev, Phys

    J. Linsley, Evidence for a primary cosmic-ray particle with energy 1020 ev, Phys. Rev. Lett. 10 (1963) 146

  4. [4]

    Greisen, End to the Cosmic-Ray Spectrum?, PhRvL 16 (1966) 748

    K. Greisen, End to the Cosmic-Ray Spectrum?, PhRvL 16 (1966) 748

  5. [5]

    Zatsepin and V .A

    G.T. Zatsepin and V .A. Kuz’min,Upper Limit of the Spectrum of Cosmic Rays, Soviet Journal of Experimental and Theoretical Physics Letters 4 (1966) 78

  6. [6]

    Blumenthal, Energy Loss of High-Energy Cosmic Rays in Pair-Producing Collisions with Ambient Photons, PhRvD 1 (1970) 1596

    G.R. Blumenthal, Energy Loss of High-Energy Cosmic Rays in Pair-Producing Collisions with Ambient Photons, PhRvD 1 (1970) 1596

  7. [7]

    R. Svensson, Non-thermal pair production in compact X-ray sources : first-order Compton cascades in soft radiation fields., Monthly Notices of the Royal Astronomical Society 227 (1987) 403

  8. [8]

    Zdziarski, G

    A.A. Zdziarski, G. Ghisellini, I.M. George, R. Svensson, A.C. Fabian and C. Done, Electron-Positron Pairs, Compton Reflection, and the X-Ray Spectra of Active Galactic Nuclei, The Astrophysical Journal Letters 363 (1990) L1

  9. [9]

    Sikora, J.G

    M. Sikora, J.G. Kirk, M.C. Begelman and P. Schneider, Electron Injection by Relativistic Protons in Active Galactic Nuclei, The Astrophysical Journal Letters 320 (1987) L81

  10. [10]

    Begelman, B

    M.C. Begelman, B. Rudak and M. Sikora, Consequences of Relativistic Proton Injection in Active Galactic Nuclei, The Astrophysical Journal 362 (1990) 38

  11. [11]

    Kirk and A

    J.G. Kirk and A. Mastichiadis, X-ray flares from runaway pair production in active galactic nuclei, Natur 360 (1992) 135

  12. [12]

    Energy spectra of gamma-rays, electrons and neutrinos produced at interactions of relativistic protons with low energy radiation

    S.R. Kelner and F.A. Aharonian, Energy spectra of gamma-rays, electrons and neutrinos produced at interactions of relativistic protons with low energy radiation, 0803.0688 [astro-ph] . – 28 –

  13. [13]

    Petropoulou and A

    M. Petropoulou and A. Mastichiadis, Bethe-Heitler emission in BL Lacs: filling the gap between X-rays and gamma-rays, Monthly Notices of the Royal Astronomical Society 447 (2015) 36

  14. [14]

    Blazar Flares as an Origin of High-Energy Cosmic Neutrinos?

    K. Murase, F. Oikonomou and M. Petropoulou, Blazar Flares as an Origin of High-energy Cosmic Neutrinos?, The Astrophysical Journal 865 (2018) 124 [1807.04748]

  15. [15]

    Leptohadronic blazar models applied to the 2014-15 flare of TXS 0506+056

    X. Rodrigues, S. Gao, A. Fedynitch, A. Palladino and W. Winter, Leptohadronic Blazar Models Applied to the 2014-2015 Flare of TXS 0506+056, The Astrophysical Journal Letters 874 (2019) L29 [1812.05939]

  16. [16]

    Zhang, M

    B.T. Zhang, M. Petropoulou, K. Murase and F. Oikonomou, A Neutral Beam Model for High-energy Neutrino Emission from the Blazar TXS 0506+056, The Astrophysical Journal 889 (2020) 118 [1910.11464]

  17. [17]

    Zhang and K

    B.T. Zhang and K. Murase, Nuclear and electromagnetic cascades induced by ultra-high-energy cosmic rays in radio galaxies: implications for Centaurus A, Monthly Notices of the Royal Astronomical Society 524 (2023) 76 [2302.14048]

  18. [18]

    Dimitrakoudis, S., Mastichiadis, A., Protheroe, R. J. and Reimer, A., The time-dependent one-zone hadronic model - first principles, A&A 546 (2012) A120

  19. [19]

    A hadronic origin for ultra-high-frequency-peaked BL Lac objects

    M. Cerruti, A. Zech, C. Boisson and S. Inoue, A hadronic origin for ultra-high-frequency-peaked BL Lac objects, Monthly Notices of the Royal Astronomical Society 448 (2015) 910 [1411.5968]

  20. [20]

    Zacharias, A

    M. Zacharias, A. Reimer, C. Boisson and A. Zech, ExHaLe-jet: An extended hadro-leptonic jet model for blazars. I. Code description and initial results, Monthly Notices of the Royal Astronomical Society 512 (2022) 3948

  21. [21]

    Stathopoulos, M

    S.I. Stathopoulos, M. Petropoulou, G. Vasilopoulos and A. Mastichiadis, LeHaMoC: a versatile time-dependent lepto-hadronic modeling code for high-energy astrophysical sources, arXiv e-prints (2023) arXiv:2308.06174 [2308.06174]

  22. [22]

    Klinger, A

    M. Klinger, A. Rudolph, X. Rodrigues, C. Yuan, G. Fichet de Clairfontaine, A. Fedynitch et al., AM3: An Open-Source Tool for Time-Dependent Lepto-Hadronic Modeling of Astrophysical Sources, arXiv e-prints (2023) arXiv:2312.13371 [2312.13371]

  23. [23]

    Propagation of Ultra High Energy Protons over Cosmological Distances and Implications for Topological Defect Models

    R.J. Protheroe and P.A. Johnson, Propagation of ultra high energy protons and gamma rays over cosmological distances and implications for topological defect models, Astroparticle Physics 4 (1996) 253 [astro-ph/9506119]

  24. [24]

    Mastichiadis and J

    A. Mastichiadis and J. Kirk, Synchrotron radiation of high-energy electrons in blazars, Astronomy & Astrophysics 295 (1995) 613

  25. [25]

    S. Gao, M. Pohl and W. Winter, On the Direct Correlation between Gamma-Rays and PeV Neutrinos from Blazars, The Astrophysical Journal 843 (2017) 109 [1610.05306]

  26. [26]

    Spectral and temporal signatures of ultrarelativistic protons in compact sources

    A. Mastichiadis, R.J. Protheroe and J.G. Kirk, Spectral and temporal signatures of ultrarelativistic protons in compact sources, astro-ph/0501156

  27. [27]

    Inoue and F

    S. Inoue and F. Takahara, Electron Acceleration and Gamma-Ray Emission from Blazars, ApJ 463 (1996) 555

  28. [28]

    Rybicki and A.P

    G.B. Rybicki and A.P. Lightman, Radiative processes in astrophysics, Haruard-Smithsonian Center for Astrophysics (1979)

  29. [29]

    Maraschi, G

    L. Maraschi, G. Ghisellini and A. Celotti, A Jet Model for the Gamma-Ray–emitting Blazar 3C 279, The Astrophysical Journal Letters 397 (1992) L5

  30. [30]

    Bloom and A.P

    S.D. Bloom and A.P. Marscher, An Analysis of the Synchrotron Self-Compton Model for the Multi–Wave Band Spectra of Blazars, The Astrophysical Journal 461 (1996) 657

  31. [31]

    Variability in the synchrotron self-Compton model of blazar emission

    A. Mastichiadis and J.G. Kirk, Variability in the synchrotron self-Compton model of blazar emission., Astronomy & Astrophysics 320 (1997) 19 [astro-ph/9610058]. – 29 –

  32. [32]

    Sikora, M.C

    M. Sikora, M.C. Begelman and M.J. Rees, Comptonization of Diffuse Ambient Radiation by a Relativistic Jet: The Source of Gamma Rays from Blazars?, The Astrophysical Journal 421 (1994) 153

  33. [33]

    Dermer, R

    C.D. Dermer, R. Schlickeiser and A. Mastichiadis, High-energy gamma radiation from extragalactic radio sources., Astronomy & Astrophysics 256 (1992) L27

  34. [34]

    TeV Gamma Rays from BL Lac Objects due to Synchrotron Radiation of Extremely High Energy Protons

    F.A. Aharonian, TeV gamma rays from BL Lac objects due to synchrotron radiation of extremely high energy protons, New Astronomy 5 (2000) 377 [astro-ph/0003159]

  35. [35]

    A Proton Synchrotron Blazar Model for Flaring in Markarian~501

    A. Mücke and R.J. Protheroe, A proton synchrotron blazar model for flaring in Markarian 501, Astroparticle Physics 15 (2001) 121 [astro-ph/0004052]

  36. [36]

    Mannheim and P.L

    K. Mannheim and P.L. Biermann, Gamma-ray flaring of 3C 279 : a proton-initiated cascade in the jet ?, Astronomy & Astrophysics 253 (1992) L21

  37. [37]

    Photohadronic origin of $\gamma$-ray BL Lac emission: implications for IceCube neutrinos

    M. Petropoulou, S. Dimitrakoudis, P. Padovani, A. Mastichiadis and E. Resconi, Photohadronic origin of γ -ray BL Lac emission: implications for IceCube neutrinos, Monthly Notices of the Royal Astronomical Society 448 (2015) 2412 [1501.07115]

  38. [38]

    A.A. Abdo, M. Ackermann, I. Agudo, M. Ajello, H.D. Aller, M.F. Aller et al., The Spectral Energy Distribution of Fermi Bright Blazars, The Astrophysical Journal 716 (2010) 30 [0912.2040]

  39. [39]

    Petropoulou and C.D

    M. Petropoulou and C.D. Dermer, Properties of Blazar Jets Defined by an Economy of Power, The Astrophysical Journal Letters 825 (2016) L11

  40. [40]

    Ghisellini and P

    G. Ghisellini and P. Madau, On the origin of the gamma-ray emission in blazars, Monthly Notices of the Royal Astronomical Society 280 (1996) 67

  41. [41]

    Chodorowski, Zdziarski and Sikora, Reaction rate and energy-loss rate fo photopair production by relativistic nuclei,

  42. [42]

    Dermer and G

    C.D. Dermer and G. Menon, High Energy Radiation from Black Holes: Gamma Rays, Cosmic Rays, and Neutrinos (2009)

  43. [43]

    Coppi and R.D

    P.S. Coppi and R.D. Blandford, Reaction rates and energy distributions for elementary processes in relativistic pair plasmas, Monthly Notices of the Royal Astronomical Society 245 (1990) 453

  44. [44]

    Ghisellini, Radiative Processes in High Energy Astrophysics, vol

    G. Ghisellini, Radiative Processes in High Energy Astrophysics, vol. 873 (2013), 10.1007/978-3-319-00612-3

  45. [45]

    Franceschini, G

    A. Franceschini, G. Rodighiero and M. Vaccari, The extragalactic optical-infrared background radiations, their time evolution and the cosmic photon-photon opacity, Astronomy & Astrophysics 487 (2008) 837

  46. [46]

    Mrk 421 as a case study for TeV and X-ray variability in leptohadronic models

    A. Mastichiadis, M. Petropoulou and S. Dimitrakoudis, Mrk 421 as a case study for TeV and X-ray variability in leptohadronic models, Monthly Notices of the Royal Astronomical Society 434 (2013) 2684 [1304.2957]

  47. [47]

    Time-dependent modelling of PKS 2155-304 in a low state

    M. Petropoulou, Time-dependent modelling of PKS 2155-304 in a low state, Astronomy & Astrophysics 571 (2014) A83 [1409.3739]

  48. [48]

    Polkas, M

    M. Polkas, M. Petropoulou, G. Vasilopoulos, A. Mastichiadis, C.M. Urry, P. Coppi et al., A numerical study of long-term multiwavelength blazar variability, Monthly Notices of the Royal Astronomical Society 505 (2021) 6103 [2105.07030]

  49. [49]

    Hadronic supercriticality as a trigger for GRB emission

    M. Petropoulou, S. Dimitrakoudis, A. Mastichiadis and D. Giannios, Hadronic supercriticality as a trigger for γ-ray burst emission, Monthly Notices of the Royal Astronomical Society 444 (2014) 2186 [1407.2915]

  50. [50]

    Mastichiadis, I

    A. Mastichiadis, I. Florou, E. Kefala, S.S. Boula and M. Petropoulou, A roadmap to hadronic supercriticalities: a comprehensive study of the parameter space for high-energy astrophysical sources, Monthly Notices of the Royal Astronomical Society 495 (2020) 2458 [2003.06956]. – 30 –

  51. [51]

    Y . Chen, Q. Gu, J. Fan, X. Yu, N. Ding, D. Xiong et al.,Curvature of the Spectral Energy Distribution, Compton Dominance, and Synchrotron Peak Frequency in Jetted Active Galactic Nuclei, The Astrophysical Journal 944 (2023) 157

  52. [52]

    The power of relativistic jets is larger than the luminosity of their accretion disks

    G. Ghisellini, F. Tavecchio, L. Maraschi, A. Celotti and T. Sbarrato, The power of relativistic jets is larger than the luminosity of their accretion disks, Natur 515 (2014) 376 [1411.5368]

  53. [53]

    Hadronic models of blazars require a change of the accretion paradigm

    A.A. Zdziarski and M. Bottcher, Hadronic models of blazars require a change of the accretion paradigm., Monthly Notices of the Royal Astronomical Society 450 (2015) L21 [1501.06124]

  54. [54]

    Petropoulou and C.D

    M. Petropoulou and C.D. Dermer, PROPERTIES OF BLAZAR JETS DEFINED BY AN ECONOMY OF POWER, The Astrophysical Journal 825 (2016) L11

  55. [55]

    Magnetically Arrested Disk: An Energetically Efficient Accretion Flow

    R. Narayan, I.V . Igumenshchev and M.A. Abramowicz,Magnetically Arrested Disk: an Energetically Efficient Accretion Flow, PASJ 55 (2003) L69 [astro-ph/0305029]

  56. [56]

    Efficient Generation of Jets from Magnetically Arrested Accretion on a Rapidly Spinning Black Hole

    A. Tchekhovskoy, R. Narayan and J.C. McKinney, Efficient generation of jets from magnetically arrested accretion on a rapidly spinning black hole, Monthly Notices of the Royal Astronomical Society 418 (2011) L79 [1108.0412]

  57. [57]

    Blandford and A

    R.D. Blandford and A. Königl, Relativistic jets as compact radio sources., The Astrophysical Journal 232 (1979) 34

  58. [58]

    Magnetic Driving of Relativistic Outflows in Active Galactic Nuclei. I. Interpretation of Parsec-Scale Accelerations

    N. Vlahakis and A. Königl, Magnetic Driving of Relativistic Outflows in Active Galactic Nuclei. I. Interpretation of Parsec-Scale Accelerations, The Astrophysical Journal 605 (2004) 656 [astro-ph/0310747]

  59. [59]

    Lynds, A.N

    C.R. Lynds, A.N. Stockton and W.C. Livingston, New Spectroscope Observations of Quasi-Stellar Sources., The Astrophysical Journal 142 (1965) 1667

  60. [60]

    Minute-Timescale >100 MeV gamma-ray variability during the giant outburst of quasar 3C 279 observed by Fermi-LAT in 2015 June

    M. Ackermann, R. Anantua, K. Asano, L. Baldini, G. Barbiellini, D. Bastieri et al., Minute-timescale >100 MeV γ-Ray Variability during the Giant Outburst of Quasar 3C 279 Observed by Fermi-LAT in 2015 June, The Astrophysical Journal Letters 824 (2016) L20 [1605.05324]

  61. [61]

    A hadronic minute-scale GeV flare from quasar 3C 279?

    M. Petropoulou, K. Nalewajko, M. Hayashida and A. Mastichiadis, A hadronic minute-scale GeV flare from quasar 3C 279?, Monthly Notices of the Royal Astronomical Society 467 (2017) L16 [1612.05699]

  62. [62]

    Elitzur, The obscuring torus in AGN, NewAR 50 (2006) 728

    M. Elitzur, The obscuring torus in AGN, NewAR 50 (2006) 728

  63. [63]

    Roychowdhury, E.T

    A. Roychowdhury, E.T. Meyer, M. Georganopoulos, P. Breiding and M. Petropoulou, Circumnuclear Dust in AP Librae and the Source of Its VHE Emission, The Astrophysical Journal 924 (2022) 57

  64. [64]

    emcee: The MCMC Hammer

    D. Foreman-Mackey, D.W. Hogg, D. Lang and J. Goodman, emcee: The MCMC Hammer, PASP 125 (2013) 306 [1202.3665]. A An empirical function for the Bethe-Heitler pair production spectrum Our goal is to provide an analytical function for the Bethe-Heitler pair production spectrum that could be easily implemented in numerical codes. Our function is benchmarked w...