pith. sign in

arxiv: 1906.09782 · v1 · pith:ONAFKMW3new · submitted 2019-06-24 · 🌌 astro-ph.HE · astro-ph.GA

Connection between optical and radio/millimeter flares in blazar OJ287

S.J Qian This is my paper

Pith reviewed 2026-05-25 17:28 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords OJ287blazar flaresoptical variationsradio variationssuperluminal knotshelical motionlighthouse effectrelativistic jets
0
0 comments X

The pith

Optical and radio light curves of blazar OJ287 decompose into 36 symmetric elementary flares produced by helical motion of superluminal knots.

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

The paper examines the connection between optical and radio/millimeter variations in the blazar OJ287 during its 1995.8-1996.1 outburst. It demonstrates that the multi-wavelength light curves can be decomposed into 36 individual flares, each with a symmetric profile. These flares arise from the lighthouse effect due to the helical motion of superluminal knots in the relativistic jet. This interpretation links the emissions across wavelengths and suggests helical knot motion is common in the inner jet regions formed near a supermassive black hole system. The work provides a framework for understanding the periodic optical outbursts observed in this source.

Core claim

It is shown that the multi-wavelength light curves at optical V-band and radio/mm wavelengths (37, 22, 14.5 and 8 GHz) can be decomposed into 36 individual elementary flares, each of which has a symmetric profile. The elementary flares can be understood to be produced through lighthouse effect due to the helical motion of corresponding superluminal optical/radio knots. Helical motion of superluminal knots should be prevailing in the inner regions of its relativistic jet formed in the magnetosphere of the putative supermassive black hole/accretion disk system. A comprehensive and compatible framework for understanding the entire phenomena in OJ287 is described.

What carries the argument

Lighthouse effect from helical motion of superluminal knots, which produces symmetric flare profiles in the light curves.

If this is right

  • The helical motion of superluminal knots prevails in the inner regions of the relativistic jet.
  • This decomposition explains the multi-wavelength connections observed in the 1995.8-1996.1 outburst.
  • The model is compatible with the quasi-periodic double-peaked optical outbursts with a 12-year cycle.
  • Emissions at optical V-band and radio/mm frequencies arise from the same knot motions in the jet.

Where Pith is reading between the lines

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

  • Similar decompositions into symmetric flares could be tested in other blazars to check for widespread helical jet structures.
  • Future monitoring of OJ287 outbursts could verify whether the 36-flare pattern repeats consistently.
  • The approach implies that jet emission is dominated by discrete knot passages rather than smooth continuous processes.

Load-bearing premise

The observed light curves can be accurately represented as the sum of exactly 36 independent symmetric elementary flares with no significant overlap or contributions from other emission mechanisms.

What would settle it

A new outburst observation where the light curves at these wavelengths cannot be fit by 36 symmetric flares or show asymmetric profiles would challenge the decomposition.

Figures

Figures reproduced from arXiv: 1906.09782 by S.J Qian.

Figure 1
Figure 1. Figure 1: Left panel: Scenario for a precessing jet and helical motion. The straight line denotes the precessing jet axis (projected on the plane of the sky) which is defined by precession phase ω=–2.0 rad. The helix indicates the trajectory of the optical knot moving around the jet axis in the perfect collimation zone for the elementary flare in 1995.98 (peaking at JD2450075). The corresponding Lorentz and Doppler … view at source ↗
Figure 2
Figure 2. Figure 2: Modeled spectrum of the elementary optical/radio flare in 1995.98 (peaking at JD2450075; filled black circles). In order to simulate the radio light-curve the spectral indexes between the optical and radio frequencies are selected to be: α37,v≡α(37GHz,V-band)=0.51, α22,v≡α(22GHz,V-band)=0.49, α15,v≡α(14.5GHz,V-band)=0.48 and α8,v≡ α(8GHz,V￾band)=0.45, respectively. The K-, U- and UV-band flux densities (op… view at source ↗
Figure 3
Figure 3. Figure 3: Quasi-periodic optical outburst in 1995.84: optical V-band light curve (upper panel) and radio/millimeter light curves observed at 37/22 GHz (middle panel) and 14.5/8 GHz (bottom panel). Zero point of time=Julian date 2450000 (1995.7735). 7 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Only five data-points are available and four el [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: Model simulation of the optical light curve for the 1995.8 optical outburst event (1995.84–1996.10): the light curve is displayed in three segments, separately. The model parameters are listed in [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Model simulation of the 37 GHz light curve during the outburst event between 1995.84 and 1996.10 (JD2450024– 2540118). The model parameters are listed in [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Model simulation of the 14.5 GHz light curve. The model parameters are listed in [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
read the original abstract

Blazar OJ287 is a unique source in which optical outbursts with double-peak structure have been observed quasi-periodically with a cycle of 12yr. It may be one of the best candidates for searching supermassive black hole binaries. We investigate the connection between its optical and radio/millimeter variations and interpret the emissions in terms of relativistic jet models. Specifically, we make a detailed analysis and model simulation of the optical and radio/mm light curves for the outburst during the period of 1995.8--1996.1. It is shown that the multi-wavelength light curves at optical V-band and radio/mm wavelengths (37, 22, 14.5 and 8 GHz) can be decomposed into 36 individual elementary flares, each of which has a symmetric profile. The elementary flares can be understood to be produced through lighthouse effect due to the helical motion of corresponding superluminal optical/radio knots. Helical motion of superluminal knots should be prevailing in the inner regions of its relativistic jet formed in the magnetosphere of the putative supermassive black hole/accretion disk system. A comprehensive and compatible framework for understanding the entire phenomena in OJ287 is described.

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 analyzes the 1995.8--1996.1 outburst in blazar OJ287, claiming that the optical V-band and radio/mm light curves (37, 22, 14.5, and 8 GHz) can be decomposed into exactly 36 individual elementary flares, each with a symmetric profile. These are interpreted as arising via the lighthouse effect from helical motion of superluminal knots, providing a framework linking the flares to jet structure in a putative supermassive black hole binary system.

Significance. If the decomposition is demonstrated to be unique and robust with quantitative support, the result would supply a concrete multi-wavelength model for helical jet motions in OJ287 and strengthen connections between periodic optical outbursts and radio/mm emission. The elementary-flare decomposition approach could inform jet studies in other blazars if the fitting methodology is made reproducible.

major comments (2)
  1. [Abstract / model simulation] Abstract and model-simulation description: the central claim that the light curves 'can be decomposed into 36 individual elementary flares' supplies no fitting procedure, degrees of freedom, goodness-of-fit metrics, or justification for selecting precisely 36 flares rather than other counts. Without these, the decomposition cannot be assessed for uniqueness or overfitting.
  2. [Interpretation of elementary flares] Helical-motion interpretation: the inference that the flares result from lighthouse-effect helical knot motion rests entirely on the untested premise that the observed curves are the sum of exactly 36 independent symmetric components with negligible overlap or other emission mechanisms. No comparison to alternative models (asymmetric profiles, different flare numbers, or additional components) is described, so the physical conclusion is not shown to follow uniquely from the data.
minor comments (1)
  1. Notation for the five frequencies and the 1995.8--1996.1 interval is clear in the abstract but should be repeated with explicit units and reference epochs when the light-curve figures are introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments. We address each major point below and will revise the manuscript accordingly to provide the requested details on methodology and model comparisons.

read point-by-point responses

  1. Referee: [Abstract / model simulation] Abstract and model-simulation description: the central claim that the light curves 'can be decomposed into 36 individual elementary flares' supplies no fitting procedure, degrees of freedom, goodness-of-fit metrics, or justification for selecting precisely 36 flares rather than other counts. Without these, the decomposition cannot be assessed for uniqueness or overfitting.

    Authors: We agree that the original manuscript does not describe the decomposition procedure in sufficient detail. The 36 flares were identified by matching symmetric components to distinct peaks visible across the optical V-band and radio/mm (37, 22, 14.5, 8 GHz) light curves during 1995.8--1996.1, with the count determined by the number of such features that could be consistently aligned at all frequencies. In the revised version we will add an explicit description of the identification criteria, the effective degrees of freedom, and quantitative goodness-of-fit measures (e.g., reduced chi-squared) to allow evaluation of uniqueness and overfitting. revision: yes

  2. Referee: [Interpretation of elementary flares] Helical-motion interpretation: the inference that the flares result from lighthouse-effect helical knot motion rests entirely on the untested premise that the observed curves are the sum of exactly 36 independent symmetric components with negligible overlap or other emission mechanisms. No comparison to alternative models (asymmetric profiles, different flare numbers, or additional components) is described, so the physical conclusion is not shown to follow uniquely from the data.

    Authors: The lighthouse-effect interpretation follows from the observed symmetry of the components and their simultaneous appearance at optical and radio frequencies, consistent with helical trajectories of superluminal knots. We acknowledge that the manuscript does not present explicit comparisons with asymmetric profiles or alternative flare counts. The revised manuscript will include a dedicated discussion comparing the symmetric 36-component model against asymmetric and alternative-count models, showing why the former provides a more consistent description of the multi-frequency data. revision: yes

Circularity Check

0 steps flagged

No significant circularity; standard model fit to multi-band light curves

full rationale

The paper performs a model simulation in which the 1995.8--1996.1 light curves are decomposed into 36 symmetric elementary flares whose parameters are adjusted to reproduce the observed data at five frequencies. This decomposition is presented as an empirical representation that is then interpreted as consistent with the lighthouse effect from helical knot motion. No step claims a first-principles derivation of the flare count or profiles independent of the fit; the number 36 and the symmetric shapes are outputs of the fitting procedure itself. No self-citation chain, uniqueness theorem, or ansatz imported from prior work is invoked to force the result. The analysis is therefore a conventional forward-modeling exercise whose central claim (that such a decomposition is possible and compatible with helical motion) does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim depends on fitting multiple parameters to match the light-curve shapes and on domain assumptions about flare symmetry and production mechanism; no independent evidence or external benchmarks are provided in the abstract.

free parameters (2)
  • Number of elementary flares = 36
    Set to 36 to reproduce the observed light-curve structure
  • Amplitude, duration, and peak time for each of the 36 flares
    Individually adjusted to fit the multi-wavelength data points
axioms (2)
  • domain assumption Each elementary flare has a symmetric profile
    Invoked as the basis for decomposing the light curves in the 1995.8--1996.1 analysis
  • domain assumption Flares arise from lighthouse effect of helically moving superluminal knots
    Used to interpret the decomposed flares as physical jet features

pith-pipeline@v0.9.0 · 5735 in / 1624 out tokens · 50859 ms · 2026-05-25T17:28:30.552175+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

57 extracted references · 57 canonical work pages · 3 internal anchors

  1. [1]

    Ackermann M., Ajello M., Allafort A., et al., 2011, ApJ 743, 1 71

    work page 2011

  2. [2]

    2012, ApJ, 7 47, 63

    Agudo, I., Marscher, A.P., Jorstad, S.G., et al. 2012, ApJ, 7 47, 63

    work page 2012

  3. [3]

    1982, MNRAS, 199, 969

    Allen, D.A., W ard, M.J., Hyland, A.R. 1982, MNRAS, 199, 969

    work page 1982

  4. [4]

    1994, in: W orkshop on Intensive Monitoring of OJ287, Tuorla Obs

    Aller, M.D., Aller, H.D., Hughes, P.A. 1994, in: W orkshop on Intensive Monitoring of OJ287, Tuorla Obs. Rep. Informo 174 , ed. M.R. Kidger & L.O. Takalo (Turku: Univ. Turku), 60

    work page 1994

  5. [5]

    Fermi Meets Jansky - AGN in Radio and Gamma-rays

    Aller, M.F., Hughes, P.A., Aller, H.D. 2010, in “Fermi Meets Jansky - AGN in Radio and Gamma-rays”, Eds.: Savolainen, T., Ros, E. ,

    work page 2010

  6. [6]

    2014, ApJ, 791 , 53

    Aller, M.F., Hughes, P.A., Aller, H.D., et al. 2014, ApJ, 791 , 53

    work page 2014

  7. [7]

    1996, ApJ, 467, L77

    Artymowicz, P., Lubow, S.H. 1996, ApJ, 467, L77

    work page 1996

  8. [8]

    1998, in: Theory of Black Hole Accretion Disk s, ed

    Artymowicz, P. 1998, in: Theory of Black Hole Accretion Disk s, ed. M.A. Abramowicz, G. Bj¨ ornsson, J.E. Pringle, p202

    work page 1998

  9. [9]

    http://altamira.asu.cas.cz/iblwg/d ata/oj287, 2006

    Basta, M., Hudec, L. http://altamira.asu.cas.cz/iblwg/d ata/oj287, 2006

    work page 2006

  10. [10]

    2010, Physics-Uspekhi, 53, 1199 Bj¨ ornsson, C.I

    Beskin, V.C. 2010, Physics-Uspekhi, 53, 1199 Bj¨ ornsson, C.I. 1982, ApJ, 260, 855 Bj¨ ornsson, C.I. & Blumenthal, G.R. 1982, ApJ, 259, 802

    work page 2010

  11. [11]

    1977, MNRAS, 179, 433

    Blandford, R.D., & Znajek, R.L. 1977, MNRAS, 179, 433

    work page 1977

  12. [12]

    Blandford, R.D., K¨ onigl A., 1979, ApJ 232, 34

    work page 1979

  13. [13]

    1982, MNRAS, 199, 883

    Blandford, R.D., & Payne, D.G. 1982, MNRAS, 199, 883

    work page 1982

  14. [14]

    2018, 4 78, 3199

    Britzen, S., Fendt, C., Witzel, G., Qian, S.J., et al. 2018, 4 78, 3199

    work page 2018

  15. [15]

    1989, A pJ, 340, 150

    Brown, L.M.J., Robson, E.I., Gear, W.K., Smith, M.G. 1989, A pJ, 340, 150

    work page 1989

  16. [16]

    1990, Reviews in Modern Astronomy, Vol.3, 234

    Camenzind, M. 1990, Reviews in Modern Astronomy, Vol.3, 234

    work page 1990

  17. [17]

    & Krockenberger, M

    Camenzind, M. & Krockenberger, M. 1992, A&A, 255, 59

    work page 1992

  18. [18]

    2017, Galaxies, 5, 12 D’Arcangelo, F.D., Marscher, A.P., Jorstad, S.D., et al

    Cohen, M.H. 2017, Galaxies, 5, 12 D’Arcangelo, F.D., Marscher, A.P., Jorstad, S.D., et al. 20 09, ApJ, 697, 985

    work page 2017

  19. [19]

    2018, ApJ, 866 , article id

    Dey, L., Valtonen, M.J., Gopakumar, A., et al. 2018, ApJ, 866 , article id. 11D

    work page 2018

  20. [20]

    2012, Scie nce, 338, 355 D’Orazio, D.L., Haiman, Z., & Macfadyen A

    Doeleman, S.S., Fish, V.L., Schenck, D.E., et al. 2012, Scie nce, 338, 355 D’Orazio, D.L., Haiman, Z., & Macfadyen A. 2013, MNRAS, 436, 2997

    work page 2012

  21. [21]

    1916, Sitzungberichte der K¨ oniglich Presssi chen Akademie der Wissenshafte (SPA W, Berlin), 688

    Einstein, A. 1916, Sitzungberichte der K¨ oniglich Presssi chen Akademie der Wissenshafte (SPA W, Berlin), 688

    work page 1916

  22. [22]

    1918, Sitzungberichte der K¨ oniglich Presssichen Akademi der Wissenshafte (SPA W, Berlin), 154

    Einstein, A. 1918, Sitzungberichte der K¨ oniglich Presssichen Akademi der Wissenshafte (SPA W, Berlin), 154

    work page 1918

  23. [23]

    2016, MNRAS, 458 , 1127

    Gupta, A.C., Agarwal, A., Mishra, A., et al. 2016, MNRAS, 458 , 1127

    work page 2016

  24. [24]

    Hartman, R.C., Bertsch, D.L., Bloom, S.D., et al.1999, ApJS , 123, 79

    work page 1999

  25. [25]

    & Ho, L.C

    Hayasaki, K., Mineshige, S. & Ho, L.C. 2008, ApJ, 682, 1134 16 S.J. Qian: Connection between optical and radio/millimete r flares in blazar OJ287

    work page 2008

  26. [26]

    2017 , A&A, 597, 80

    Hodgson, J.A., Krichbaum, T.P., Marscher, A.P., et al. 2017 , A&A, 597, 80

    work page 2017

  27. [27]

    Distance measures in cosmology

    Hogg, D.W. 1999, astro-ph/9905116

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  28. [28]

    1984, MNR AS, 211, 497

    Holmes, P.A., Brand, P.W.J.L., Impey, C.D.,et al. 1984, MNR AS, 211, 497

    work page 1984

  29. [29]

    1988, A&A, 190, L8 K¨ onigl, A

    Kikuchi, S., Inoue, M., Mikami, Y., et al. 1988, A&A, 190, L8 K¨ onigl, A. & Choudhuri, A.R. 1985, ApJ, 289, 173

    work page 1988

  30. [30]

    2009, ApJS, 180 , 330

    Komatsu, E., Dunkley, J., Nolta, M.R., et al. 2009, ApJS, 180 , 330

    work page 2009

  31. [31]

    1996, ApJ, 460, 207

    Lehto, H.J., & Valtonen, M.J. 1996, ApJ, 460, 207

    work page 1996

  32. [32]

    1992, ApJ, 394, 459

    Li, Z.Y., Chiueh, T.H., Begelman, M.C. 1992, ApJ, 394, 459

    work page 1992

  33. [33]

    2001, Science, Vol.291, 8 4

    Meier, D.L., Koide, S., Uchida, Y. 2001, Science, Vol.291, 8 4

    work page 2001

  34. [34]

    2013, EPJ W eb of Conference 61, 01001

    Meier, D.L. 2013, EPJ W eb of Conference 61, 01001

    work page 2013

  35. [35]

    & Vovk, Ie

    Neronov, A. & Vovk, Ie. 2011, MNRAS, 412, 1389

    work page 2011

  36. [36]

    1991, Acta As tron

    Qian, S.J., Witzel, A., Krichbaum, T.P., et al. 1991, Acta As tron. Sin., 32, 369 (english translation: in Chin. Astro. Astroph ys., 16, 137 (1992))

    work page 1991

  37. [37]

    2013, Res

    Qian, S.J. 2013, Res. Astron. Astrophys., 13, 783

    work page 2013

  38. [38]

    2015, Res

    Qian, S.J. 2015, Res. Astron. Astrophy., 15, 687

    work page 2015

  39. [39]

    2017, A&A, 604, A9 0

    Qian, S.J., Britzen, S., Witzel, A., et al. 2017, A&A, 604, A9 0

    work page 2017

  40. [40]

    A tentative double-jet model for blazar OJ287

    Qian, S.J. 2018, arXiv-1811.11514

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  41. [41]

    Model simulation of optical light curves for blazar OJ287

    Qian, S.J. 2019, arXiv-1904.03357 Schramm K.-J., Borgeest U., Camenzind M.,et al., 1993, A&A 2 78,391

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  42. [42]

    2012, ApJ , 749, 118 Sillanp¨ a¨ a, A., Teerikorpi, P., Haarala, S., et al

    Shi, J.M., Krolik, J.H., Loubow, S.H., Hawley J.F. 2012, ApJ , 749, 118 Sillanp¨ a¨ a, A., Teerikorpi, P., Haarala, S., et al. 1985, A&A, 147, 67 Sillanp¨ a¨ a, A., Haarala, S., Valtonen, M.J., et al. 1988, ApJ, 325, 628 Sillanp¨ a¨ a, A., Takalo, L.O., Pursimo, T., et al. 1996a, A&A, 315, L13 Sillanp¨ a¨ a, A., Takalo, L.O., Pursimo, T., et al. 1996b, AS...

    work page 2012

  43. [43]

    2003, ApJS, 14 8, 175

    Spergel, D.N., Verde, L., Peiris, H.V., et al. 2003, ApJS, 14 8, 175

    work page 2003

  44. [44]

    1997, ApJ, 484, 1 80

    Sundelius, B., W ade, M., Lehto, H.J., et al. 1997, ApJ, 484, 1 80

    work page 1997

  45. [45]

    2013, MNRAS, 434, 2275

    Tanaka, T.L. 2013, MNRAS, 434, 2275

    work page 2013

  46. [46]

    1999, Ap J, 520, 627

    Tateyama, C.E., Kingham, K.A., Kayfmann, P., et al. 1999, Ap J, 520, 627

    work page 1999

  47. [47]

    1979, Astron

    Usher, P.D. 1979, Astron. J., 84, 1253

    work page 1979

  48. [48]

    1987, A&A, 184, 57

    Valtaoja, L., Sillanp¨ a¨ a, A., & Valtaoja, E. 1987, A&A, 184, 57

    work page 1987

  49. [49]

    2000, A pJ, 531, 744

    Valtaoja, E., Terasranta, H., Tornikoski, M., et al. 2000, A pJ, 531, 744

    work page 2000

  50. [50]

    2006, ApJ, 646, 36

    Valtonen, M.J., Lehto, H.J., Sillanp¨ a¨ a, A., et al. 2006, ApJ, 646, 36

    work page 2006

  51. [51]

    2007, ApJ, 659, 1074

    Valtonen, M.J. 2007, ApJ, 659, 1074

    work page 2007

  52. [52]

    & Ciprini, S

    Valtonen, M.J. & Ciprini, S. 2012, Mem.S.A.It., 83, 219

    work page 2012

  53. [53]

    2018, in: Gravitat ional W aves Astrophysics: Early Results from Gravitational W ave Searches and Electromagnetic Counterparts, Proceedings o f IAU

    Valtonen, M.J., Dey, L., Hudec, R., et al. 2018, in: Gravitat ional W aves Astrophysics: Early Results from Gravitational W ave Searches and Electromagnetic Counterparts, Proceedings o f IAU

    work page 2018

  54. [54]

    2017, Galaxies, 5, 83

    Valtonen, M.L., Zola, S., Jermak, H., et al. 2017, Galaxies, 5, 83

    work page 2017

  55. [55]

    1998, MN RAS 293, L13

    Villata, M., Raiteri, C.M., Sillanp¨ a¨ a, A., et al. 1998, MN RAS 293, L13

    work page 1998

  56. [56]

    2010, MNRAS 402 , 2087

    Villforth, C., Nilsson, K., Heidt, J., et al. 2010, MNRAS 402 , 2087

    work page 2010

  57. [57]

    2004, ApJ, 605, 656 W agner S.J., Camenzind M., Dreissigacker O., et al., 1995, A &A 298, 688 17

    Vlahakis, N., & K¨ onigl, A. 2004, ApJ, 605, 656 W agner S.J., Camenzind M., Dreissigacker O., et al., 1995, A &A 298, 688 17

    work page 2004