INTEGRAL/ISGRI post 2024-periastron view of PSR B1259-63

A. Kuzin; A. Santangelo; B. van Soelen; D. Malyshev; M. Chernyakova

REVIEW 3 major objections 5 minor 28 references

INTEGRAL/ISGRI post 2024-periastron view of PSR B1259-63

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A single power law describes PSR B1259-63's 0.3-300 keV X-ray spectrum after the 2024 periastron passage.

desk verdict New 2024 X-ray data tighten the PSR B1259-63 spectrum, but the missing cross-normalization in the joint fits probably makes the quoted errors optimistic. read the letter →

arxiv 2505.21474 v1 pith:K4ZQYQL6 submitted 2025-05-27 astro-ph.HE

A. Kuzin , D. Malyshev , M. Chernyakova , B. van Soelen , A. Santangelo This is my paper

classification astro-ph.HE

keywords PSRB1259-63gamma-raybinariesX-rayspectroscopyINTEGRAL/ISGRISwift/XRTpulsarwindsynchrotronemissioninverseComptonscattering

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper establishes that the non-thermal X-ray emission from the gamma-ray binary PSR B1259-63, observed 17-33 days after the June 2024 periastron passage, is a single absorbed power law from 0.3 to 300 keV with photon index $\Gamma=1.42\pm0.03$. A joint analysis of soft and hard X-ray data excludes, at 95% confidence, a spectral break below about 30 keV for index changes larger than 0.3 and an exponential cutoff below 150 keV. Interpreting the X-rays as synchrotron radiation in a one-zone leptonic model, the paper concludes the radiating electrons have a break energy above $E_{\rm b,e}>0.8$ TeV and a cutoff above $E_{\rm cutoff,e}>1.7$ TeV. These results tighten the link between the observed X-ray spectrum and particle acceleration in the pulsar wind shock.

What carries the argument

The analysis rests on a joint spectral fit of soft X-ray (0.3-10 keV) and hard X-ray (30-300 keV) observations taken in July-August 2024, with nine quasi-simultaneous pairs around the periastron passage. Stacking datasets 4-9 (17-33 days after periastron) gives the tight power-law index $\Gamma=1.42\pm0.03$. The curved alternatives are tested with a broken power law (Equation 1) and an exponential cutoff power law (Equation 2), and the $\chi^2$ difference maps provide the break and cutoff lower limits. The connection to the electron distribution is made through the synchrotron peak relation $E_{\rm photon}\approx52$ keV $(B/1\,{\rm G})(E_e/1\,{\rm TeV})^2$, with the magnetic field $B\approx2$ G fixed by a one-zone leptonic fit that combines the X-ray spectrum with TeV data and uses super-exponential cutoff electron spectra with index $\beta=1$ or $\beta=3$.

What would settle it

A future, higher-sensitivity observation of PSR B1259-63 in the 30-300 keV band that finds a spectral break below about 30 keV with $\Delta\Gamma>0.3$, or an exponential cutoff below 150 keV, would refute the single-power-law claim; an independent determination that the emission-zone magnetic field is much weaker than 2 G would invalidate the translated electron limits.

Watch

Extended reading notes

Core claim

The central claim is that, in the 2024 post-periastron phase, the 0.3-300 keV spectrum of PSR B1259-63 is well described by an absorbed power law with photon index $\Gamma=1.42\pm0.03$ and hydrogen column density $N_H=(0.86^{+0.12}_{-0.09})\times10^{22}$ cm$^{-2}$, and that no spectral break or cutoff is required in this band. For the broken power-law model, break energies below $E_b\approx10$ keV are excluded for $\Delta\Gamma>0.1$, and break energies below $E_b\approx30$ keV for $\Delta\Gamma>0.3$; for an exponential cutoff model, $E_{\rm cutoff}<150$ keV is excluded at 95% confidence. If the emission is synchrotron from a single electron population with an inverse-Compton component matching the TeV spectrum, the magnetic field in the emission zone is about 2 G and the electron distribution has break and cutoff energies above $E_{\rm b,e}>0.8$ TeV and $E_{\rm cutoff,e}>1.7$ TeV. The paper thus claims that particle acceleration in the intrabinary shock produces electrons up to the TeV range with no cooling break in the X-ray band.

Load-bearing premise

The electron break and cutoff limits follow from the photon limits only if the X-rays are synchrotron emission from one electron population whose inverse-Compton emission also produces the TeV spectrum with a magnetic field of about 2 G; if the magnetic field or the TeV spectrum is different, the quoted particle energies do not follow.

Editorial extensions

If this is right

The hard and soft X-ray bands are described by one power law, so the synchrotron component spans from 0.3 keV to at least 300 keV without a spectral break.
The parent electron population must carry a single power-law distribution up to electron energies of about a TeV, with break and cutoff energies at $E_{\rm b,e}>0.8$ TeV and $E_{\rm cutoff,e}>1.7$ TeV under the assumed field.
A magnetic field of about 2 G in the emission zone requires either an emission site very close to the pulsar or a significant stellar magnetic field contribution, both of which can be tested with orbital-phase modeling.
The gradual hardening of the photon index from about 1.6 to 1.35 over the monitored interval, while not statistically significant, matches the expectation from earlier periastron campaigns and indicates that the effective emission geometry changes as the pulsar recedes from the star.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The limits on the electron distribution are sensitive to the assumed magnetic field: the same photon break limit would correspond to a higher particle energy if the field were weaker, so an independent measurement of $B$ would move the quoted $E_{\rm b,e}$ and $E_{\rm cutoff,e}$ bounds.
The paper's one-zone interpretation assumes the X-ray and TeV emitting regions are the same electron population; if they are spatially separate, the magnetic field inferred from the combined SED would not apply cleanly to the X-ray-emitting electrons.
Because the TeV data used to fix the field come from a previous periastron, a contemporaneous TeV observation in a future orbit would test whether the magnetic field and electron cutoff are stable across periastron passages or vary with the orbital cycle.
The featureless power law from 0.3 to 300 keV implies that, if the one-zone picture is right, radiation losses do not imprint curvature in the X-ray band; a future detection of curvature would localize the cooling break and constrain the acceleration-site density.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Colleague,

The paper gives you a new, well-measured hard X-ray spectrum for PSR B1259-63 after the 2024 periastron, with a photon index of 1.42 ± 0.03 and 2-sigma lower limits of ~30 keV for a spectral break (ΔΓ > 0.3) and ~150 keV for a cutoff. The joint XRT-ISGRI analysis is careful: they handle the nearby 2RXP J130159.6-635806, fit N_H and Γ with proper error contours, and check broken/cutoff models with Δχ² maps. That part is solid and is a real step beyond prior work.

The main thing to probe is the cross-normalization. The text says the fit uses a single normalization for the PL; it does not mention a free constant between XRT and ISGRI. Since the two instruments observe up to a day apart and the source flux swings by about a factor of two in that span, a 10–20% relative offset (within ISGRI's calibration uncertainty) would be absorbed into the slope and can shift Γ by more than 0.03. That effectively raises the systematic floor on the headline number and on the break/cutoff limits. The paper needs to either include a cross-normalization, or explicitly argue it is negligible.

Two smaller points. They stack only datasets 4–9, calling them 'high-quality' with no quantitative threshold; a sentence on why that subset is justified would help. And the electron break/cutoff limits assume a magnetic field in the conversion formula B ≈ 1 G (they later fit B ≈ 2 G) and combine 2024 X-ray data with 2021 H.E.S.S. TeV data. The electron parameters are model-translated, not measured; this is fine if labeled clearly, but the B dependence should not be buried.

None of this kills the central power-law result, which is robust. The concerns are about error bars and the model-dependent translation. I'd send it to a referee. The referee should ask for a cross-normalization statement, a clearer dataset-selection rule, and a more explicit caveat on the electron limits. With those, the paper will be worth citing.

Referee Report

3 major / 5 minor

Summary. This paper analyzes INTEGRAL/ISGRI and quasi-simultaneous Swift/XRT observations of the gamma-ray binary PSR B1259-63/LS 2883 obtained 9–33 days after the 2024 periastron. The authors perform joint spectral fits with an absorbed power law, finding that the 0.3–300 keV spectrum of datasets 4–9 is well described by a photon index Gamma = 1.42 ± 0.03. They then use Delta-chi-squared maps to limit spectral curvature, excluding a break above about 30 keV for Delta-Gamma > 0.3 and an exponential cutoff below about 150 keV at 95% confidence. Interpreting these limits in a one-zone synchrotron/inverse-Compton framework and combining the X-ray data with 2021 H.E.S.S. TeV spectra, they infer electron spectral parameters and a magnetic field B about 2 G.

Significance. If the X-ray results are robust, the paper provides the hardest X-ray spectral coverage of this system after a periastron passage and yields a clean constraint on the location of any spectral break or cutoff, which is valuable for particle acceleration models. The analysis is careful in several respects: the contribution of the nearby source 2RXP J130159.6-635806 is modeled using contemporaneous NuSTAR data, the confidence contours are obtained with Delta-chi-squared rather than simple 1-sigma errors, and the SED modeling is explicitly separated from the direct X-ray spectroscopy. The main caveats are the apparent omission of a cross-normalization constant in the joint fits and an internal inconsistency in the electron limits, both of which are fixable in revision.

major comments (3)

[§2, Figs. 2–3] The joint XRT+ISGRI spectral fits appear to use a single power-law normalization with no cross-normalization constant between the two instruments. Because the XRT and ISGRI observations in a given dataset are separated by up to one day (Table 1) and the source flux varies by roughly a factor of two on these timescales (Fig. 2, top panel), any relative offset between the two instruments at the 10–20% level (typical of ISGRI absolute calibration) will be absorbed into the photon index and can shift Gamma by an amount comparable to the quoted ±0.03. The break and cutoff limits in Fig. 4 are likewise sensitive to the relative normalization of the 30–80 keV and 80–300 keV ISGRI bins with respect to the XRT spectrum. I ask the authors to redo the fits including a multiplicative cross-normalization constant (e.g., a *const in XSPEC) and to report its best-fit value and the resulting changes to Gamma and the Fig. 4 contours.
[§3, after Eq. (3), and Eqs. (4)–(5)] The translation of the X-ray break/cutoff limits into electron spectral parameters uses E_ph ≈ 52 keV (B/1 G)(E_e/1 TeV)^2, but the quoted values E_b,e > 0.8 TeV and E_cutoff,e > 1.7 TeV are obtained only if B = 1 G is assumed. Later in the same section the SED fit yields B ≈ 2 G. With B = 2 G, the same photon limits translate to E_b,e > 0.54 TeV and E_cutoff,e > 1.2 TeV. The abstract and conclusions should either adopt a consistent B and propagate its uncertainty, or explicitly state that the quoted electron limits depend on the assumed B value and are not direct measurements.
[§3, Fig. 5] The one-zone SED fit uses the 2021 H.E.S.S. TeV spectrum to represent the 2024 TeV state, while the X-ray data are from 2024. Because the paper itself notes that the X-ray light curve differs between periastron passages and that the previous 2017 modeling yielded B five to ten times weaker, the derived B ≈ 2 G and the consequent electron limits are conditional on an untested assumption. I ask the authors to state this assumption more prominently and to discuss how period-to-period variations would affect B and the electron limits.

minor comments (5)

[§3 and §4] The text refers to the stacked time interval as '18–32 days' in one place and '17–33 days' in the abstract and Table 1; please make the numbers consistent.
[§4] The phrase 'the cutoff upper limit E_cutoff ≳ 27 TeV reported by H.E.S.S.' is contradictory: the symbol ≳ denotes a lower limit, so rephrase as 'lower limit'.
[Eq. (4)] The numerical coefficient in Eq. (4) is not dimensionally transparent; please specify the normalization of the magnetization parameter sigma and of the distance r_PE so that the G unit follows explicitly.
[§2] There is a typo: 'lies well within the ISGRI PSR for our target' should read 'ISGRI PSF'.
[§2] The choice of only one ISGRI energy bin in the 80–300 keV range limits sensitivity at the high-energy end; this should be acknowledged as a caveat in the discussion of the break/cutoff limits.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the X-ray spectral results are direct fits to the data, and the electron-parameter translation is model-dependent but is not fed back into the spectral derivation.

full rationale

The paper's central quantitative claims (Gamma = 1.42 +/- 0.03; E_b > 30 keV for DeltaGamma > 0.3; E_cutoff > 150 keV at 2 sigma) are obtained by direct XSPEC fits to Swift/XRT and INTEGRAL/ISGRI spectra, with the observational data coming from outside the paper. No parameter is fitted to the target result and then renamed a prediction. The broken-PL and cutoff-PL scans in Fig. 4 are independent model comparisons against the same data, and fixing Gamma_1 = 1.42 to the joint best-fit value is a standard, explicitly verified choice. The later translation to electron spectral parameters (p = 2Gamma - 1; E_b,e > 0.8 TeV; E_cutoff,e > 1.7 TeV) is a derived application of the standard synchrotron relation with an assumed magnetic field; the B ~ 2 G value from the Naima SED fit is itself a fit to the X-ray data plus an independent external TeV dataset, and the electron limits are not fed back into the X-ray spectral analysis. Self-citations (Chernyakova et al. 2021, 2024, 2025) supply prior light-curve context, XRT data reduction, and Fermi spectra; they are observational inputs or prior context, not unverified uniqueness theorems or ansatze used to force the present conclusion. The H.E.S.S. 2021 TeV spectrum is an external, independent dataset. Thus no circular step is present; the stated model dependencies (one-zone assumption, magnetic field choice) are scientific assumptions, not cases where the paper's input is equivalent to its output by construction.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central X-ray spectral measurements require only the standard free parameters N_H, Gamma, and normalization. The model-dependent electron constraints additionally assume the one-zone leptonic scenario and rely on the magnetic field B fitted to the SED, with the 2021 TeV data serving as the TeV anchor.

free parameters (5)

Interstellar hydrogen column density N_H = 0.86 (+0.12, -0.09) x 10^22 cm^-2 for stacked datasets 4-9; individual values in Table 1
Fitted to the X-ray absorption in each spectrum; part of the measured result, not hidden.
Photon index Gamma = 1.42 +/- 0.03 (stacked datasets 4-9)
Main fitted spectral slope; is itself the central measurement.
Magnetic field B in emission zone = ~2 G
Fitted with Naima to the combined X-ray and 2021 TeV SED; carries the model-dependent electron constraints.
Electron distribution normalization = not quoted
Normalization of the electron spectrum in the one-zone model.
Super-exponential cutoff index beta = 1 or 3 (two fixed choices)
Chosen by hand, not fit; affects the sharpness of the cutoff but not the constraint derivation.

assumptions (5)

domain assumption X-rays are synchrotron and TeV gamma rays are inverse Compton from a single electron population (one-zone leptonic model)
Used to translate the photon break/cutoff limits into electron energy limits and to fit B in Section 3.
domain assumption The electron distribution follows a super-exponential cutoff power law, Eq. (3), with beta = 1 or 3
The shape is assumed, not derived; other shapes would alter the inferred limits.
ad hoc to paper The 2021 H.E.S.S. TeV spectrum represents the 2024 TeV state
The SED fit combines 2021 TeV data with 2024 X-ray data; known orbit-to-orbit variability makes this uncertain.
domain assumption Gamma-gamma absorption is negligible in the TeV band
Stated explicitly in Section 3; it affects the TeV model fit but not the X-ray constraints.
domain assumption The contaminating source 2RXP J130159.6-635806 is accurately described by the NuSTAR model of Salganik et al. (2025)
Used to subtract the nearby source contribution in ISGRI spectra.

how reviews work

0 comments

Cite this review

Pith. "Pith review of INTEGRAL/ISGRI post 2024-periastron view of PSR B1259-63." pith.science (2026). https://pith.science/paper/K4ZQYQL6

@misc{pith2026250521474,
  author       = {Pith},
  title        = {Pith review of: INTEGRAL/ISGRI post 2024-periastron view of PSR B1259-63},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/K4ZQYQL6}},
  note         = {Machine review of arXiv:2505.21474}
}

abstract

PSR B1259-63/LS 2883 is a well-studied gamma-ray binary hosting a pulsar in a 3.4-year eccentric orbit around a Be-type star. Its non-thermal emission spans from radio to TeV energies, exhibiting a significant increase near the periastron passage. This paper is dedicated to the analysis of INTEGRAL observations of the system following its last periastron passage in June 2024. We aim to study the spectral evolution of this gamma-ray binary in the soft (0.3-10 keV) and hard (30-300 keV) X-ray energy bands. We performed a joint analysis of the data taken by INTEGRAL/ISGRI in July-August 2024 and quasi-simultaneous Swift/XRT observations. The spectrum of the system in the 0.3-300 keV band is well described by an absorbed power law with a photon index of $\Gamma=1.42\pm 0.03$. We place constraints on potential spectral curvature, limiting the break energy $E_\mathrm{b}>30$ keV for $\Delta\Gamma>0.3$ and cutoff energy $E_\mathrm{cutoff}>150$ keV at 95% confidence level. For one-zone leptonic emission models, these values correspond to electron distribution spectral parameters of $E_\mathrm{b,e}>0.8$ TeV and $E_\mathrm{cutoff,e}>1.7$ TeV, consistent with previous constraints derived by H.E.S.S.

Figures

Figures reproduced from arXiv: 2505.21474 by the authors.

**Figure 3.** Chi-square (χ 2 ) maps in the Γ–NH parameter space. The color scale represents the joint fit χ 2 values, with regions below 3σ significance shown. White contours indicate the 1σ and 2σ confidence levels. Red contours correspond to χ 2 levels from XRT-only fits (see labels). The analysis includes data from datasets 4 − 9 only. We note that a magnetic field of the same order may originate from the optical star. We es… view at source ↗

**Figure 2.** Top: the 2024 light curve in soft X-rays (blue) and hard X [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: The difference in total χ 2 between a modified PL model and a simple PL model. Left: the modified model is a broken PL, with the map shown in coordinates of break energy Eb and index difference ∆Γ = Γ2 − Γ1. The first index was fixed at Γ1 = 1.42, while the second inde…

**Figure 5.** Figure 5: Top: SED of PSR B1259-63. Black: Swift/XRT stacked spectrum for the datasets 4–9. Orange: INTEGRAL/ISGRI stacked spectrum for the same datasets. Red bound: Fermi 2024 spectrum 0 − 20 days after periastron. Green: Fermi spectrum during gamma-ray flares: 19 − 77 days af…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

28 extracted references · 22 canonical work pages

[1]

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

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
[3]

A., Kelner , S

Aharonian , F. A., Kelner , S. R., & Prosekin , A. Y. 2010, , 82, 043002

work page 2010
[4]

N., Hill , J

Burrows , D. N., Hill , J. E., Nousek , J. A., et al. 2005, , 120, 165

2005
[5]

M., Takata , J., Yi , S

Chen , A. M., Takata , J., Yi , S. X., Yu , Y. W., & Cheng , K. S. 2019, , 627, A87

work page 2019
[6]

& Malyshev , D

Chernyakova , M. & Malyshev , D. 2020, in Multifrequency Behaviour of High Energy Cosmic Sources - XIII. 3-8 June 2019. Palermo, 45

work page 2020
[7]

2020, , 497, 648

Chernyakova , M., Malyshev , D., Mc Keague , S., et al. 2020, , 497, 648

work page 2020
[8]

2025, , 536, 247

Chernyakova , M., Malyshev , D., van Soelen , B., et al. 2025, , 536, 247

work page 2025

Show all 28 references

[9]

2024, , 528, 5231

Chernyakova , M., Malyshev , D., van Soelen , B., et al. 2024, , 528, 5231

2024
[10]

2021, Universe, 7, 242

Chernyakova , M., Malyshev , D., van Soelen , B., et al. 2021, Universe, 7, 242

2021
[11]

2015, , 454, 1358

Chernyakova , M., Neronov , A., van Soelen , B., et al. 2015, , 454, 1358

2015
[12]

2007, , 14, 171

De Becker , M. 2007, , 14, 171

2007
[13]

H. E. S. S. Collaboration , Aharonian , F., Ait Benkhali , F., et al. 2024, , 687, A219

2024
[14]

N., Lyne , A

Johnston , S., Manchester , R. N., Lyne , A. G., et al. 1992, , 387, L37

1992
[15]

Kennel , C. F. & Coroniti , F. V. 1984, , 283, 694

1984
[16]

A., Bogovalov , S

Khangulyan , D., Aharonian , F. A., Bogovalov , S. V., & Rib \'o , M. 2011, , 742, 98

2011
[17]

A., & Kelner , S

Khangulyan , D., Aharonian , F. A., & Kelner , S. R. 2014, , 783, 100

2014
[18]

Melatos , A., Johnston , S., & Melrose , D. B. 1995, , 275, 381

1995
[19]

Miller-Jones , J. C. A., Deller , A. T., Shannon , R. M., et al. 2018, , 479, 4849

2018
[20]

2011, , 732, L11

Negueruela , I., Rib \'o , M., Herrero , A., et al. 2011, , 732, L11

2011
[21]

S., Chernyakova , M., Malyshev , D., & Poutanen , J

Salganik , A., Tsygankov , S. S., Chernyakova , M., Malyshev , D., & Poutanen , J. 2025, arXiv e-prints, arXiv:2504.11263

2025 arXiv
[22]

E., Chernyakova , M., Rodriguez , J., et al

Shaw , S. E., Chernyakova , M., Rodriguez , J., et al. 2004, , 426, L33

2004
[23]

& Arons , J

Tavani , M. & Arons , J. 1997, , 477, 439

1997
[24]

Tavani , M., Arons , J., & Kaspi , V. M. 1994, , 433, L37

1994
[25]

P., & Kee , N

ud-Doula , A., Owocki , S. P., & Kee , N. D. 2018, , 478, 3049

2018
[26]

Usov , V. V. & Melrose , D. B. 1992, , 395, 575

1992
[27]

A., Petit , V., Grunhut , J., & Neiner , C

Wade , G. A., Petit , V., Grunhut , J., & Neiner , C. 2014, arXiv e-prints, arXiv:1411.6165

2014 arXiv
[28]

2015, in International Cosmic Ray Conference, Vol

Zabalza , V. 2015, in International Cosmic Ray Conference, Vol. 34, 34th International Cosmic Ray Conference (ICRC2015), 922

2015

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [1]

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

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

[2] [2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

[3] [3]

A., Kelner , S

Aharonian , F. A., Kelner , S. R., & Prosekin , A. Y. 2010, , 82, 043002

work page 2010

[4] [4]

N., Hill , J

Burrows , D. N., Hill , J. E., Nousek , J. A., et al. 2005, , 120, 165

2005

[5] [5]

M., Takata , J., Yi , S

Chen , A. M., Takata , J., Yi , S. X., Yu , Y. W., & Cheng , K. S. 2019, , 627, A87

work page 2019

[6] [6]

& Malyshev , D

Chernyakova , M. & Malyshev , D. 2020, in Multifrequency Behaviour of High Energy Cosmic Sources - XIII. 3-8 June 2019. Palermo, 45

work page 2020

[7] [7]

2020, , 497, 648

Chernyakova , M., Malyshev , D., Mc Keague , S., et al. 2020, , 497, 648

work page 2020

[8] [8]

2025, , 536, 247

Chernyakova , M., Malyshev , D., van Soelen , B., et al. 2025, , 536, 247

work page 2025