{"id":"e43ecaaf-7a9a-4987-988e-5bc8d1544589","arxiv_id":"1908.08393","paper_version":1,"verdict":"UNVERDICTED","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The CTA collaboration outlines its scheduling and real-time analysis pipeline for following up gravitational-wave alerts with gamma-ray telescopes, supported by illustrative simulations but no derived detection rates.","lead":"This conference paper describes the planned gravitational-wave follow-up program of the Cherenkov Telescope Array and shows example simulations of gamma-ray counterpart searches. A generalist might read it to see how next-generation gamma-ray observatories plan to respond to neutron-star merger alerts, though the paper stops short of predicting detection rates.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 8σ pipeline validation may be inflated by simulating the VHE spectrum without EBL absorption: Sec. 3.2 extrapolates a 0.1–10 GeV power law to 10 TeV, yet the text never applies extragalactic background light attenuation for the 500 Mpc GWCOSMoS sample.","rationale":"The reader's UNVERDICTED verdict is appropriate: this ICRC contribution is a program description with an illustrative injection study, and the authors explicitly defer the GW-EM detection-rate derivation to future work (Section 4). I nevertheless found one load-bearing weak point in the strongest concrete claim, namely that the 8σ detection in Fig. 3 validates the low-latency follow-up pipeline. The validation depends on the Sec. 3.2 source model, which extrapolates a 0.1–10 GeV power law to 10 TeV without any stated EBL absorption, despite the mock catalog reaching 500 Mpc. At z≈0.1, EBL absorption suppresses multi-TeV photons by a large factor; ignoring it inflates the VHE flux and can turn a sub-threshold source into an 8σ detection. A single Gammapy rerun with EBL absorption would determine whether the 8σ claim survives. This is a concrete, testable concern rather than general skepticism about sGRB phenomenology. It does not change the UNVERDICTED status, because the paper does not present a quantitative science result, but it does mean the simulation evidence for the program's capabilities should be treated as provisional. The reader's weakest_assumption also centers on the Sec. 3.2 emission model, so there is partial agreement, but my concern targets a specific missing physical ingredient (EBL) rather than the broader phenomenological extrapolation.","tokens_in":6877,"tokens_out":7749,"duration_ms":71864,"concrete_test":"Repeat the Section 3.3 injection used for Fig. 3 in Gammapy, adding an EBL attenuation model (e.g., Franceschini et al. 2008 or Domínguez et al. 2011) to the Sec. 3.2 power-law spectrum at the actual redshift/distance of the injected GWCOSMoS event, and recompute the TS map. If the peak significance drops from 8σ to below the 5σ threshold used in Eq. 2.1, the simulation demonstration is not robust and the claim of validated low-latency detectability needs qualification. If the injected event is very nearby and the significance remains at or above 5σ, the concern is resolved for that event but should still be checked for the farthest events in the sample.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest concrete demonstration is the 8σ detection of an injected GRB in a Gammapy simulation (Sec. 3.3, Fig. 3). That demonstration inherits the source model of Sec. 3.2, where the 0.1–10 GeV Fermi-LAT behavior of GRB 090510 is extrapolated as a power law with photon index -2.1 up to 10 TeV. The text does not state that EBL attenuation is applied to this extrapolation, although the mock BNS sample extends to 500 Mpc (Sec. 3.1). At the farthest distances, the EBL optical depth at several TeV is of order unity or larger, so an unattenuated extrapolation can overproduce the VHE photon flux by a factor e^τ. If EBL was omitted, the 8σ excess in Fig. 3 is not a realistic detection for a distant sGRB, and the claimed validation of the low-latency follow-up pipeline is weaker than presented. This is a correctness risk in the central simulation, not merely a disagreement with current sGRB models; a rerun with EBL absorption would settle it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper describes the planned gravitational-wave (GW) follow-up program of the Cherenkov Telescope Array (CTA). It outlines an observation scheduling strategy that optimizes low-energy coverage, sky-probability coverage, and dynamic observation windows, and it describes a Real-Time Analysis (RTA) system that is claimed to issue science alerts with latencies below 30 seconds. The paper then presents simulations built from the public GWCOSMoS catalog of mock BNS mergers, a phenomenological model of very-high-energy (VHE) emission from short GRBs based on the Fermi-LAT behavior of GRB 090510, and Gammapy-based CTA observations using instrument response functions. The central demonstration is a simulated follow-up observation in which an injected GRB is reported as detected at 8-sigma significance (Section 3.3, Figure 3). The paper explicitly states that the derivation of GW-EM detection rates is left to future work.","tokens_in":7135,"tokens_out":4772,"duration_ms":51739,"significance":"If the program performs as described, the paper is a useful description of an important multimessenger follow-up capability. The use of public simulation tools and databases (GWCOSMoS, Gammapy, CTA IRFs) is a strength, and the simulated detection is an internally consistent proof-of-concept that the pipeline can find an injected source. However, the concrete 8-sigma claim rests on a single simulated observation and on several unstated or unvalidated modeling choices, so the validation is weaker than the abstract and Section 3.3 suggest. The paper does not yet provide the detection-rate estimates that would motivate the follow-up program quantitatively; it explicitly defers them. The contribution is therefore best read as a program description and a preliminary pipeline demonstration, not as a measured sensitivity statement.","major_comments":[{"comment":"The simulation extrapolates a 0.1-10 GeV power law with photon index -2.1 up to 10 TeV but does not state whether extragalactic background light (EBL) attenuation is applied. Since the GWCOSMoS sample extends to 500 Mpc and the EBL optical depth is of order unity above several TeV for such distances, an unattenuated extrapolation can overproduce the expected VHE flux and inflate the reported 8-sigma significance. The authors should either explicitly state that EBL absorption is included, with the adopted model, or rerun the simulation with EBL attenuation and report the resulting significance.","section":"Section 3.2, Figure 3"},{"comment":"The claimed 8-sigma detection is based on a single simulated observation of a single injected event, and the paper provides no test-statistic distribution, no background-only comparison, and no accounting for the number of trials or hot-spots searched. Without this statistical context, the significance cannot be interpreted and the statement that the pipeline is validated is not fully supported. At minimum, the authors should specify how the 8-sigma value is derived from the TS analysis and show that it is not an artifact of a particular noise realization or of multiple testing across the field of view.","section":"Section 3.3, Figure 3"},{"comment":"The simulated VHE signal is built entirely from a phenomenological model of one short GRB (GRB 090510): a t^2 rise, an onset at 3 s, a t^-1.4 decay, a photon index of -2.1, and a homogeneous jet with a 5-degree opening angle extrapolated to 10 TeV. Because this model is used to generate the injected source in the 8-sigma validation, the simulation is conditional on that unvalidated assumption. The paper should either present a bracketing study over plausible short-GRB VHE models or state explicitly that the detection claim holds only under this model, so that readers can judge the robustness of the pipeline demonstration.","section":"Section 3.2"}],"minor_comments":[{"comment":"The short GRB detected by MAGIC is named '1901114C' in the abstract but '190114C' in the Introduction; the correct name, GRB 190114C, should be used consistently.","section":"Abstract and Introduction"},{"comment":"The sentence 'In order to asses the proposed strategy' contains a typo: 'asses' should be 'assess'.","section":"Section 3"},{"comment":"The notation 'Fint_5σ' is used without a definition, and the phrase 'T obs → inf' in the text should read 'T_obs → ∞'.","section":"Section 2, Eq. (2.1)"},{"comment":"Figure 2 shows six observations of 15 seconds while Figure 3 shows a single observation with Tobs=2 s; the relationship between these two exposure choices should be stated explicitly so the reader does not infer an inconsistency.","section":"Figures 2 and 3"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings paper, and the standards for a peer-reviewed journal may be different from those for an ICRC contribution. The program description is reasonable and the use of public tools is commendable. However, the paper's most concrete quantitative claim, the 8-sigma detection, is not yet backed by a statistically complete analysis, and the EBL question is a genuine correctness risk rather than a stylistic issue. A revision that addresses these two points, or that explicitly reframes the result as a pipeline demonstration rather than a validated detection, would be appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick note on arXiv:1908.08393. It's an ICRC2019 proceeding describing the CTA gravitational-wave follow-up program. The honest summary: it's a status report, not a discovery paper. The scheduling strategy (low-zenith preference, probability-ranked tiling, galaxy-catalog convolving, dynamic windows) and the RTA science-alert pipeline are described transparently, and the GWCOSMoS-based simulation bank plus Gammapy analysis is a legitimate piece of reusable infrastructure. The 8σ detection in Fig. 3 is an injection test, not a real event, and the paper itself says the detection-rate derivation is future work (Sec. 4). So credit where it's due: the program description is coherent, the simulation setup is concrete, and the GWCOSMoS catalog is public and traceable.\n\nNow the soft spots. The stress-test concern lands: the simulated VHE spectrum is extrapolated as a power law to 10 TeV with no mention of EBL attenuation, while the mock BNS sample reaches 500 Mpc. At several TeV the EBL optical depth is order unity or larger, so an unabsorbed power law can overproduce the flux by e^τ. That makes the 8σ figure optimistic if EBL was indeed omitted. This is not fatal to the paper's main purpose—the pipeline still does its job on an injected source—but it is a real omission in the simulation description, and the authors should state explicitly whether EBL was applied. There's also only one simulated injection, no error bars on the TS, and no null-hypothesis comparison. The source model leans heavily on GRB 090510 with a fixed 5° homogeneous jet; that is acceptable as a phenomenological starting point, but it is not a validated sGRB template. Minor: the abstract says GRB 1901114C; it should be 190114C. The citation pattern is fine, building on Patricelli et al. and earlier IACT follow-up studies.\n\nBottom line: if you need a snapshot of CTA's GW follow-up plans, this is a clear, useful paper. It does not overclaim—the authors explicitly defer rates to future work. I'd send it to a referee for this conference proceeding, but I'd ask for a clarification of the EBL treatment and, ideally, a null-case example to calibrate the TS significance.","headline":"A clear CTA status report that describes a sensible GW follow-up program and a reusable simulation pipeline, but the example 8σ detection likely ignores EBL attenuation and should not be read as a realistic sensitivity claim.","tokens_in":7781,"tokens_out":1992,"would_cite":false,"duration_ms":22021,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The CTA gravitational-wave follow-up program can issue real-time science alerts in under 30 seconds and recover a modeled short-GRB counterpart at 8σ in simulated observations.","keywords":["gravitational waves","very-high-energy gamma rays","short gamma-ray bursts","Cherenkov Telescope Array","multi-messenger astronomy","real-time analysis","binary neutron star mergers","follow-up scheduling"],"falsifier":"A targeted search: take the first on-axis BNS-merger short GRB observed by CTA in follow-up mode and compare its VHE light curve and spectrum to the assumed template ($t^2$ rise, $t_{\\rm peak}=3$ s, $t^{-1.4}$ decay, photon index $-2.1$). A substantial mismatch—or a non-detection where the template predicts a 5σ detection—would show the phenomenological model does not describe real sources.","tokens_in":6686,"feed_emoji":"🔭","tokens_out":7732,"duration_ms":71754,"temperature":0.7,"pith_summary":"The paper presents the Cherenkov Telescope Array's program for following up gravitational-wave alerts, and argues that the program is fast enough and sensitive enough to catch very-high-energy (VHE) counterparts of binary neutron star mergers. Its short-term scheduler chooses observing windows that favour low-energy coverage, cover the most probable sky regions first, and adapt their duration to the predicted flux decay of a short gamma-ray burst. Its real-time analysis pipeline is said to issue science alerts to external observatories in under 30 seconds. In a full end-to-end simulation built from a public catalog of simulated mergers and a phenomenological model of short-GRB VHE emission, the pipeline recovers the injected source at 8σ in a 2-second observation. If the simulation reflects reality, CTA would be able to catch and quickly broadcast VHE emission from gravitational-wave counterparts, linking GW detections to the wider multi-wavelength transient picture.","feed_headline":"Simulated GRB in CTA follow-up is recovered at 8 sigma","feed_subtitle":"Real-time pipeline issues alerts in under 30 seconds, validating CTA's low-latency GW follow-up.","key_machinery":"The machinery that carries the argument is the follow-up scheduler plus real-time analysis pipeline, combined with a simulation chain that produces the expected VHE signal. The scheduler optimises three axes: low-zenith (low-energy-threshold) observations, probability-coverage ordering over the GW localisation map, and dynamic observation windows satisfying the condition that the time-integrated flux over the window equals the 5σ sensitivity. The real-time pipeline makes a run-by-run test-statistic fit and issues alerts in under 30 seconds. The simulation chain links a sample of synthetic BNS mergers to a phenomenological short-GRB VHE model ($t^2$ rise to $t_{\\rm peak}=3$ s, $t^{-1.4}$ decay, $E^{-2.1}$ spectrum to 10 TeV, homogeneous jet with a 5-degree opening angle, normalised by the isotropic-energy-to-GeV-luminosity correlation) and folds the signal through CTA instrument response functions. This machinery converts a GW alert into a concrete observation schedule and a detection decision, and produces the 8σ recovery that supports feasibility.","core_discovery":"The paper's central claim is that CTA can conduct a complete low-latency gravitational-wave follow-up: upon receiving an alert, a short-term scheduler computes observation windows optimised along three axes—low zenith angle to keep the energy threshold low, a probability-ordered tiling of the gravitational-wave localisation map (optionally weighted by a galaxy catalog), and dynamic window lengths set by the time needed to reach a 5σ detection given the assumed source decay. The accompanying real-time analysis pipeline is stated to detect sub-minute emission and to issue science alerts below 30 seconds. To test the strategy, the authors simulate BNS mergers with public waveforms and detection criteria, attach a VHE light curve to each: the flux rises as $t^2$ to a peak at 3 s, decays as $t^{-1.4}$, has a photon index of $-2.1$, is extrapolated to 10 TeV, normalised via a correlation between isotropic energy and GeV luminosity, and corrected for viewing angle assuming a homogeneous 5-degree jet. Simulating a CTA observation of one such event with instrument response functions and analysing it run-by-run with a test-statistic fit, they recover the modeled GRB at 8σ. They do not yet report joint detection rates—that is the stated next step—but the simulation infrastructure is explicitly built for that purpose.","pith_inferences":["If the true VHE emission from short GRBs is fainter or decays faster than the GRB 090510 template, the 8σ recovery and any future detection-rate numbers would degrade, but the pipeline's sub-30-second alert latency would remain valid; the conclusion that CTA can respond fast is more robust than the conclusion that it will detect.","The same scheduling machinery could be applied to other poorly localised transients, such as neutrino alerts, where probability-ordered tiling and low-energy prioritisation are also the limiting factors.","A direct way to benchmark the pipeline is to replay archival VHE observations of GRB 190114C through the same run-by-run test-statistic analysis and compare recovery significance and latency."],"forward_implications":["A BNS merger whose short GRB emits VHE light along the model's template can be caught by CTA within seconds of the alert, and the detection can be broadcast to other observatories in under 30 seconds.","The simulation bank, once detection rates are computed, will let the collaboration state how many GW events per year CTA should expect to observe at VHE, and how that number depends on source distance and orientation.","The run-by-run test-statistic analysis is appropriate for real-time transient searches, since it recovered a point source at 8σ in a 2-second exposure without needing a stacked analysis.","Favouring low-zenith angles and probability-ordered tiling increases the chance that the first observation windows cover the true counterpart."],"supporting_citations":[{"why":"Establishes the empirical basis for BNS mergers as short-GRB progenitors, motivating the follow-up program.","marker":"[1]"},{"why":"Provides the original simulated BNS merger sample and the GW detection criteria used for the mock catalog.","marker":"[17]"},{"why":"Supplies the follow-up scheduling methodology and dynamic observation-window estimates that the program extends.","marker":"[19]"},{"why":"Is cited for the real-time analysis pipeline that issues science alerts in under 30 seconds.","marker":"[24]"},{"why":"Provides the intrinsic isotropic-energy distribution for short GRBs used to assign Eiso to simulated mergers.","marker":"[38]"},{"why":"Supplies the Eiso-to-GeV-luminosity correlation used to normalise the simulated VHE light curve.","marker":"[39]"},{"why":"Gives the homogeneous-jet prescription used to correct the simulated emission for viewing angle.","marker":"[40]"},{"why":"Provides the instrument response functions from detailed Monte Carlo simulations used to simulate CTA observations.","marker":"[42]"},{"why":"Supplies the test-statistic root-finding analysis method used in the run-by-run real-time analysis.","marker":"[43]"}],"fun_headline_variants":["CTA's GW follow-up hits 8σ on simulated GRB","Low-latency CTA pipeline catches simulated GRB at 8σ","30-second alerts: CTA recovers simulated GRB at 8σ","Tiled follow-up recovers simulated GRB in CTA test","CTA's rapid response: simulated GRB detected at 8σ"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulation rests on the assumption that short-GRB very-high-energy emission follows the GRB 090510 template—the $t^2$ rise to a 3 s peak, the $t^{-1.4}$ decay, the $-2.1$ photon index, the extrapolation to 10 TeV, and the 5-degree homogeneous jet—and that real events behave this way.","fun_headline_variants_meta":{"raw":{"variants":["CTA's GW follow-up hits 8σ on simulated GRB","Low-latency CTA pipeline catches simulated GRB at 8σ","30-second alerts: CTA recovers simulated GRB at 8σ","Tiled follow-up recovers simulated GRB in CTA test","CTA's rapid response: simulated GRB detected at 8σ"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000324,"raw_usage":{"total_tokens":1871,"prompt_tokens":1051,"completion_tokens":820,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":725}},"tokens_in":667,"tokens_out":820,"duration_ms":8086,"temperature":1.0,"reasoning_tokens":725,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:39:59.167101+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A targeted search: take the first on-axis BNS-merger short GRB observed by CTA in follow-up mode and compare its VHE light curve and spectrum to the assumed template ($t^2$ rise, $t_{\\rm peak}=3$ s, $t^{-1.4}$ decay, photon index $-2.1$). A substantial mismatch—or a non-detection where the template predicts a 5σ detection—would show the phenomenological model does not describe real sources.","supporting_citations":[{"cited_title":"2017, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes the empirical basis for BNS mergers as short-GRB progenitors, motivating the follow-up program."},{"cited_title":"2016, JCAP, 11, 056","cited_arxiv_id":null,"evidence_quote":"Provides the original simulated BNS merger sample and the GW detection criteria used for the mock catalog."},{"cited_title":"2018, JCAP, 5, 056","cited_arxiv_id":null,"evidence_quote":"Supplies the follow-up scheduling methodology and dynamic observation-window estimates that the program extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is cited for the real-time analysis pipeline that issues science alerts in under 30 seconds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the intrinsic isotropic-energy distribution for short GRBs used to assign Eiso to simulated mergers."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Eiso-to-GeV-luminosity correlation used to normalise the simulated VHE light curve."},{"cited_title":"2002, ApJ, 570, L61","cited_arxiv_id":null,"evidence_quote":"Gives the homogeneous-jet prescription used to correct the simulated emission for viewing angle."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the instrument response functions from detailed Monte Carlo simulations used to simulate CTA observations."},{"cited_title":"A&A 495.3, pp","cited_arxiv_id":null,"evidence_quote":"Supplies the test-statistic root-finding analysis method used in the run-by-run real-time analysis."}],"review_version":1}