Pith. sign in

REVIEW 2 major objections 4 minor 19 references

Observation of VER J2019+368 with the SST-1M stereoscopic system

T0 review · 2 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read The SST-1M stereoscopic system detects the extended very-high-energy source VER J2019+368 at 6.5 sigma post-trial and measures its power-law spectrum with index 2.44 up to 100 TeV.

desk verdict First stereoscopic SST-1M detection of an extended source, with a genuinely useful off-axis acceptance curve; the 6.5σ detection is credible, but the spectrum rests on a single-atmosphere MC and should be treated as preliminary. read the letter →

arxiv 2507.16408 v1 pith:IAOU5EDI submitted 2025-07-22 astro-ph.HE astro-ph.IM

classification astro-ph.HEastro-ph.IM
keywords VERJ2019+368SST-1MimagingatmosphericCherenkovtelescopesvery-high-energygammaraysextendedgamma-raysourcespulsarwindnebulaGalacticPeVatroncandidates
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

The paper reports the first detection of the very-high-energy gamma-ray source VER J2019+368 with the two SST-1M Cherenkov telescopes, from 97 hours of good-quality stereoscopic data collected between April and November 2024. The source shows up at 6.5 sigma post-trial significance in an extended-source search, with a best-fit morphology described by an asymmetric Gaussian and a spectrum that follows a power law with index 2.44 ± 0.13 up to about 100 TeV. The authors stress that the key instrumental property is a nearly flat gamma-ray acceptance out to roughly 2.5 degrees off-axis, which is what allows a small telescope with a 9-degree field of view to study extended galactic sources. If the result holds, it would demonstrate that compact Cherenkov telescopes can contribute to the follow-up of the multi-TeV sources discovered by large air-shower arrays and to searches for Galactic PeVatrons.

What carries the argument

The load-bearing instrument is the stereoscopic pair of SST-1M telescopes, whose Davies-Cotton optics (a segmented spherical reflector that keeps off-axis images compact) and 9-degree field of view give an integral gamma-ray acceptance after cuts that stays almost flat out to about 2.5 degrees, dropping only about 10 percent; this flatness is checked with dedicated Crab Nebula observations at offsets from 0.7 to 3.06 degrees. That flat acceptance is what allows the analysis to recover a source with roughly 0.4-degree extension. On the statistical side, the central machinery is a 0.25-degree extended convolution kernel combined with ring-background significance maps, with post-trial significances calibrated by running the identical analysis on 3,000,000 Monte Carlo background skymaps. The spectral result comes from a reflected-background 1D fit with a power law, followed by one-zone inverse-Compton and pion-decay modeling.

What would settle it

Re-run the identical 97-hour dataset through two independent instrument response sets produced from summer and winter atmospheric profiles and recompute the spectral energy distribution; if the flux normalization moves by more than the quoted 8% or the post-trial significance falls below 5 sigma, the fixed-atmosphere correction is the weak link.

Watch

Extended reading notes

Core claim

The central discovery is that the unknown source VER J2019+368 is clearly visible to the SST-1M stereoscopic system as an extended emitter: the strongest excess is found with 8.1 $\sigma$ pre-trial local significance for a 0.25-degree extended convolution kernel, corresponding to 6.5 $\sigma$ after accounting for the trial factor in 3,000,000 simulated background skymaps, while the point-like search yields 4.8 $\sigma$ post-trial. The best-fitting spatial model is an asymmetric Gaussian with $\sigma_{\mathrm{long}} = 0.42^\circ \pm 0.06^\circ$ and $\sigma_{\mathrm{lat}} = 0.16^\circ \pm 0.03^\circ$, centered at $\alpha = 304.81^\circ \pm 0.07^\circ$, $\delta = 36.71^\circ \pm 0.04^\circ$, consistent with the extension reported by earlier imaging-air-Cherenkov observations. A 1D spectral analysis in that region gives $\Gamma = 2.44 \pm 0.13$ and $\phi_0 = (3.17 \pm 0.41) \times 10^{-14}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{TeV}^{-1}$ at $E_0 = 7\,\mathrm{TeV}$, with flux points extending to 100 TeV. The authors do not confirm the two-source substructure suggested earlier and cannot reject an energy-independent morphology, but they show the integrated spectrum is compatible with either inverse-Compton radiation from electrons accelerated to about 200 TeV or pion decay from protons accelerated to about 1 PeV.

Load-bearing premise

The analysis assumes that one fixed-atmosphere Monte Carlo production, corrected only by scaling event image sizes (Hillas intensities) above 400 photoelectrons by a single cosmic-ray rate ratio, describes the real atmospheric conditions over the whole April–November 2024 campaign; if that scaling does not capture the actual variation, the flux, spectral index, and detection significance could shift beyond the quoted systematics.

Editorial extensions

If this is right

  • VER J2019+368 is established as an extended TeV source within reach of small-aperture Cherenkov telescopes, not only the large arrays that discovered it.
  • The near-flat acceptance out to about 2.5 degrees means SST-1M can observe extended galactic sources and poorly localized transients without the usual strong off-axis penalty.
  • The power-law spectrum measured up to 100 TeV gives an independent spectral constraint that models of the pulsar wind nebula and of possible PeV proton acceleration must match.
  • The failure to reject energy-independent morphology with the current data sets a sensitivity benchmark: a longer campaign would be needed to detect any energy-dependent extension of the source.

Reading between the lines

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

  • A direct way to stress-test the result is to split the campaign by month and re-fit the spectrum with per-season atmospheric corrections; if the flux normalization moves by more than the quoted 8%, the single fixed-atmosphere Monte Carlo is the dominant systematic.
  • If the inverse-Compton interpretation with electrons up to about 200 TeV is correct, the same electron population should produce a spatially extended X-ray synchrotron counterpart whose gradient could be checked against existing pulsar wind nebula X-ray maps.
  • The one-zone modeling cannot separate leptonic from hadronic emission, so a joint spectral fit covering the SST-1M and other TeV observatories' energy ranges over one common integration region would be the next decisive step.
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.

Referee Report

2 major / 4 minor

Summary. This ICRC 2025 proceeding reports the first stereoscopic observations of the VER J2019+368 region with the SST-1M system, using 97 hours of good-quality data taken between April and November 2024. The paper validates the telescope's off-axis acceptance with Crab observations, finds a main excess with 8.1 sigma pre-trial local significance and 6.5 sigma post-trial significance in an extended-source search, fits an asymmetric Gaussian morphology, and derives a power-law spectrum with index Gamma = 2.44 ± 0.13 and normalization phi0 = (3.17 ± 0.41) × 10^-14 cm^-2 s^-1 TeV^-1 at E0 = 7 TeV. It also presents simple one-zone IC and pion-decay models that can plausibly describe the SED.

Significance. If the detection and spectral results hold, this work demonstrates that the SST-1M system, with its large field of view, can detect and characterize extended Galactic gamma-ray sources, which is a valuable capability for follow-up of LHAASO sources. The off-axis acceptance study with Crab data is a concrete strength, and the paper is appropriately cautious in labeling the results as preliminary. The simple modeling provides context but does not add strong new physical constraints. The main significance is instrumental and phenomenological, with the detection itself being the central claim.

major comments (2)
  1. [Sec. 3] The post-trial significance is computed using 3,000,000 background-only skymaps generated with IRFs for a single 30 deg zenith angle, while the actual dataset covers 12-60 deg. If the acceptance and background rate vary with zenith and season, the simulated trial distribution may not match the data, so the quoted 6.5 sigma post-trial significance could shift. The authors should either generate trial maps with a representative distribution of zenith angles or demonstrate that the result is insensitive to the assumed zenith.
  2. [Sec. 4.2] The spectral analysis adopts the asymmetric Gaussian shape derived from the same dataset in Sec. 4.1 as the signal region. This is not an independent region and the spectral parameters may be biased by the morphology fit, especially since the two-source morphology suggested by VERITAS is not confirmed. The authors should either use an external region definition, propagate the morphology uncertainties into the spectral fit, or explicitly discuss this a posteriori selection effect.
minor comments (4)
  1. [Sec. 2] The text states that the measured Crab rates in Fig. 1 confirm the MC acceptance, but it does not report the statistical significance of the agreement or the systematic uncertainties on the points; a brief statement would strengthen the validation claim.
  2. [Sec. 4.1] The energy-dependent morphology test reports ΔTS = 11.5 and says H0 cannot be rejected, but no p-value or threshold is given; stating the resulting p-value would aid interpretation.
  3. [Sec. 4.2] The SED comparison in Fig. 4 notes that integration regions are not the same across observatories; this caveat is useful, but the text could also explicitly mention that the SST-1M spectral region may blend multiple VERITAS sources, which is relevant for the physical interpretation.
  4. [References] References [10] and [15] are cited as arXiv preprints or software repositories; for a proceedings paper this is acceptable, but the authors should ensure the final version includes journal or DOI identifiers where available.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the detection and spectrum are empirical measurements calibrated with MC and validated against Crab; self-citations are methodological, not load-bearing.

full rationale

This paper is an observational analysis, not a derivation whose conclusion is built into an input assumption. The main claims (8.1 sigma pre-trial excess, 6.5 sigma post-trial extended-source significance, power-law index 2.44 +/- 0.13, and flux normalization) are obtained by applying a Monte-Carlo-based reconstruction and instrument response functions to recorded stereoscopic data, with the MC serving as calibration rather than as the target quantity. The one internal reuse, taking the asymmetric Gaussian morphology fitted in Sec. 4.1 as the fixed signal region for the 1D spectral fit in Sec. 4.2, is a standard and explicitly stated analysis choice; the spectral index and normalization are free parameters and are not algebraically forced by the morphology, the IRFs, or the exclusion regions. Self-citations to [10], [15], and [16] document the sst1mpipe pipeline and random-forest classifier, and the paper shows external validation with Crab Nebula observations and simulated gammas in Sec. 2 (Fig. 1). The acknowledged limitation in Sec. 3, that only a single fixed-atmosphere MC production is used with a single Hillas-intensity scaling for atmospheric differences, is a systematic-risk caveat rather than a circular step: the recorded data are not constructed from the MC, and the authors explicitly assign additional 8% and 5% systematics. The post-trial significance computed from 3,000,000 background MC skymaps at 30 deg zenith is a trials-factor estimate, not a derivation of the observed excess from the background model. No quoted equation or fitted parameter is identical by construction to an input, so no circularity step is identified.

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

The paper introduces no new physical entities, but its central observational claims rest on standard IACT assumptions: MC fidelity, background symmetry, and the adequacy of a single fixed-atmosphere simulation with a rate-ratio correction. The morphological and spectral fit parameters are free parameters, and the SED modeling adds fitted particle populations.

free parameters (7)
  • Power-law spectral index Gamma = 2.44 +/- 0.13
    Free parameter in the power-law fit to the SED in Sec 4.2; used to characterize the source spectrum up to 100 TeV.
  • Flux normalization phi0 = (3.17 +/- 0.41) x 10^-14 cm^-2 s^-1 TeV^-1
    Normalization fitted at reference energy 7 TeV in the same spectral fit.
  • Asymmetric Gaussian extensions = sigma_long = 0.42 +/- 0.06 deg, sigma_lat = 0.16 +/- 0.03 deg
    Spatial model fitted to the significance map in Sec 4.1 and then adopted as the spectral signal region.
  • Source position = RA = 304.81 +/- 0.07 deg, Dec = 36.71 +/- 0.04 deg
    Best-fitting centroid of the asymmetric Gaussian morphology in Sec 4.1.
  • Maximum electron energy in IC model = about 200 TeV
    Fitted via naima one-zone inverse Compton modeling of the SED in Sec 4.2.
  • Maximum proton energy in pion decay model = about 1 PeV
    Fitted via naima one-zone pion decay modeling of the SED in Sec 4.2.
  • Particle spectral indices in SED models = about 3 for electrons, about 2.4 for protons
    Injected particle spectra in the naima fits; the paper reports them as physically plausible values.
assumptions (5)
  • domain assumption Air-shower and detector Monte Carlo simulations accurately model the Cherenkov signal and are adequate for random forest reconstruction and IRF generation.
    Sec 3 uses a single fixed-atmosphere MC production for all 97 hours of data, with only a rate-ratio intensity scaling to correct data-MC differences.
  • domain assumption The wobble and ring background method assumes radially symmetric acceptance of the camera and that residual background after exclusions follows the modeled distribution.
    Sec 4.1 uses ring background with a 0.3 deg exclusion radius for known VHE sources; the authors test two variants but the assumption is intrinsic to the significance estimate.
  • ad hoc to paper The cosmic-ray rate ratio computed for events with Hillas intensities above 400 p.e. is a sufficient correction for all atmospheric and instrument variations over the campaign.
    Sec 3 introduces this scaling to match MC to data; if incomplete, flux and spectral index systematics exceed the stated 8% and 5%.
  • domain assumption The gamma and hadron separation classifier trained on Monte Carlo generalizes to real data.
    Sec 3 optimizes an energy-dependent gammaness cut on MC and applies it to data.
  • domain assumption One-zone emission models with simple power-law particle distributions can describe the gamma-ray emission of this complex region.
    Sec 4.2 fits naima IC and pion decay models and interprets the maximum energies; the authors themselves caution that multi-wavelength information and 3D analysis are needed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Observation of VER J2019+368 with the SST-1M stereoscopic system." pith.science (2026). https://pith.science/paper/IAOU5EDI

@misc{pith2026250716408,
  author       = {Pith},
  title        = {Pith review of: Observation of VER J2019+368 with the SST-1M stereoscopic system},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IAOU5EDI}},
  note         = {Machine review of arXiv:2507.16408}
}
read the original abstract

The Single-Mirror Small Size Telescope (SST-1M) is a small Cherenkov telescope designed to detect gamma rays with energies more than about 1 TeV. The optical design of the SST-1M follows the Davies-Cotton concept to ensure good off-axis performance. In 2022, two SST-1M telescope prototypes were installed in Ondrejov, Czech Republic, and stereoscopic observations of astrophysical gamma-ray sources have been performed since then. VER J2019+368 is an unidentified very-high-energy (VHE) gamma-ray source, surrounded by several gamma-ray point-like and diffuse sources, together with their multi-wavelength counterparts. VHE emission was discovered by MILAGRO in 2012, followed by VERITAS observation, which revealed the complex morphology of the source. Recently, the LHAASO observatory detected photons with multi-TeV energies, opening up the possibility of particle acceleration up to PeV energies. In this contribution, we present preliminary results of the first observing campaign of the VER J2019+368 region, performed with SST-1M from April to November 2024. We present the data analysis, focusing on the morphological and spectroscopic study of the region. We also present the off-axis performance of SST-1M in the context of the prospects for detecting extended galactic gamma-ray sources. As one of the brightest and hardest sources in the LHAASO catalog, VER J2019+368 is an ideal candidate for testing the capabilities of the SST-1M, with its large field of view, to detect extended gamma-ray sources.

Figures

Figures reproduced from arXiv: 2507.16408 by the authors.

Figure 1
Figure 1. Acceptance of SST-1M as a function of the offset angle. The lines represent the integral rate of simulated gammas with Crab Nebula spectrum in the signal region after analysis cuts as a function of the offset at 30◦ zenith angle. The points show a measured rate of gamma-like excess events resulting from a preliminary analysis of the Crab Nebula observation [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Maps of local significance for the VER J2019+368 region as seen with the SST-1M stereoscopic system. Point-like (0.1 ◦ , left) and extended (0.25◦ , right) integration radii were used to produce the skymaps. Left skymap shows the multi-wavelength counterparts, and right skymap shows the VHE sources in the region. The radii of the VHE sources correspond to their extensions as reported in [5, 6, 18]. The contours in b… view at source ↗
Figure 3
Figure 3. Left: The profile of VER J2019+368 emission along its main axis. The excess counts (blue distribution) are calculated from the skymap constructed using the point-like integration region. The red line represents the best-fitting Gaussian, and the vertical lines mark the coordinates of the three VERITAS sources along the profile. Right: The map of local significance between 2.5 TeV and 100 TeV showing the best-fitting… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Left: SED of the VER J2019+368 region as measured with SST-1M stereoscopic system, compared with results of different observatories (we note, however, that the integration regions are not the same). Right: Modeling of the SST-1M SED assuming inverse Compton (IC) and pi…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

19 extracted references · 11 canonical work pages

  1. [1]

    J.W.T.Hessels,M.S.E.Roberts, S.M.Ransom,V.M.Kaspi,R.W.Romani,C.Y.Ngetal.,ApJ612(2004)389 [astro-ph/0403632]

  2. [2]

    Rings and Jets around PSR J2021+3651: the `Dragonfly Nebula'

    A. Van Etten, R. W. Romani and C. Y. Ng,ApJ680(2008) 1417 [0802.3707]

  3. [3]

    E. Aliu, T. Aune, B. Behera, M. Beilicke, W. Benbow, K. Berger et al.,ApJ788(2014) 78 [1404.1841]

  4. [4]

    A. A. Abdo, U. Abeysekara, B. T. Allen, T. Aune, D. Berley, E. Bonamente et al.,The Astrophysical Journal753 (2012) 159

  5. [5]

    A. U. Abeysekara, A. Archer, T. Aune, W. Benbow, R. Bird, R. Brose et al.,The Astrophysical Journal861(2018) 134

  6. [6]

    Z. Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao et al.,ApJS271(2024) 25 [2305.17030]

  7. [7]

    J. Woo, H. An, J. D. Gelfand, C. J. Hailey, K. Mori, R. Mukherjee et al.,ApJ954(2023) 9 [2306.07347]

  8. [8]

    Optical observations of PSR J2021+3651 in the Dragonfly Nebula with the GTC

    A. Kirichenko, A. Danilenko, P. Shternin, Y. Shibanov, E. Ryspaeva, D. Zyuzin et al.,ApJ802(2015) 17 [1501.04594]

Show all 19 references
  1. [9]

    Alispach, A

    C. Alispach, A. Araudo, M. Balbo, V. Beshley, A. Biland, J. Blažek et al.,JCAP2025(2025) 047 [2409.11310]

  2. [10]

    Alispach, A

    C. Alispach, A. Araudo, M. Balbo, V. Beshley, J. Blažek, J. Borkowski et al.,2506.01733

  3. [11]

    Serrano et al., inICALEPCS, Kobe, Japan, 12th International Conference on Accelerator and Large Experimental Physics Control Systems, 2009

    J. Serrano et al., inICALEPCS, Kobe, Japan, 12th International Conference on Accelerator and Large Experimental Physics Control Systems, 2009

  4. [12]

    V. P. Fomin, A. A. Stepanian, R. C. Lamb, D. A. Lewis, M. Punch and T. C. Weekes,Astroparticle Physics2 (1994) 137

  5. [13]

    Aleksić, S

    J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, P. Bangale et al.,Astroparticle Physics72(2016) 76 [1409.5594]

  6. [14]

    Heller, J

    M. Heller, J. Schioppa, E., A. Porcelli, I. T. Pujadas, K. Zie ¸tara, D. d. Volpe et al.,European Physical Journal C 77(2017) 47 [1607.03412]

  7. [15]

    Jurysek et al., Feb., 2025

    J. Jurysek et al., Feb., 2025. 10.5281/zenodo.14808846

  8. [16]

    Jurysek, T

    J. Jurysek, T. Tavernier, V. Novotný, M. Heller, D. Mandat, M. Pech et al.,arXiv e-prints(2023) arXiv:2307.09799 [2307.09799]

  9. [17]

    Donath et al., Nov., 2021

    A. Donath et al., Nov., 2021. 10.5281/zenodo.5721467

  10. [18]

    A. U. Abeysekara, A. Albert, R. Alfaro, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu et al., PRL124(2020) 021102 [1909.08609]

  11. [19]

    Zabalza, in34th International Cosmic Ray Conference (ICRC2015), vol

    V. Zabalza, in34th International Cosmic Ray Conference (ICRC2015), vol. 34 ofInternational Cosmic Ray Conference, p. 922, July, 2015,1509.03319, DOI. 8 VER J2019+368 with SST-1MJ. Juryšek Full Authors List: SST-1M Collaboration C. Alispach1, A. Araudo2, M. Balbo1, V. Beshley3,...

Pith tools

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