{"id":"17c2c01a-32f4-44c0-b75a-4a86b33d9a47","arxiv_id":"2603.17511","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Serendipitous VERITAS observations do not detect NGC 4278, but upper limits bound the quasi-quiet TeV index below 3.13 and, with Fermi-LAT, bracket the SED peak between 100 GeV and 2 TeV within an example hadronic-corona model.","lead":"VERITAS, an array of gamma-ray telescopes, saw no high-energy gamma rays from the nearby galaxy NGC 4278 during and around a claimed outburst seen by the Chinese observatory LHAASO, and set flux limits that constrain the galaxy's spectrum. The paper combines those limits with Fermi satellite data to suggest the emission peaks between 100 GeV and 2 TeV, and explores a model in which protons accelerated above the black hole produce the gamma rays and a faint, currently undetect","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quasi-quiet Γ<3.13 constraint rests on unverified assumption that 6.53 h of VERITAS observations (MJD 59280–59697) sample the time-averaged LHAASO quasi-quiet state; if the quiet-state flux varies over the year, the upper limit is not representative.","rationale":"We read the paper in good faith. The empirical analysis (VERITAS non-detection, upper limit, cross-check with independent package) is sound, and the authors are transparent about the example-level nature of the corona model and about the assumption needed for the Γ<3.13 constraint. The single most load-bearing concern is the quasi-quiet averaging assumption, exactly as the reader identified. It is load-bearing because the paper's strongest quantitative claim — that VERITAS rules out indices >=3.13 in the LHAASO quasi-quiet state — depends on the 6.53 h VERITAS exposure representing the time-averaged flux of a state that spans ~2.7 years (MJD 59278–60248). The VERITAS coverage (MJD 59280–59697) is temporally limited, and the source is known to vary by a factor ~4 between active and quiet states; variability within the quiet state is plausible and would bias the comparison in an unknown direction. We considered alternatives: the log-parabolic fit's treatment of the upper limit as a zero-flux point is non-standard, but the resulting peak range (100 GeV–2 TeV) is so broad that the abstract's 'indicate a peak' is a weak claim; the private-communication flux point is explicitly not used in fits; and the corona model is labeled as an example. Thus the index constraint remains the pivotal new result. The proposed test — time-matching LHAASO flux to the VERITAS windows — is feasible with public data and would settle whether the assumption holds. If it fails, the conclusion should be downgraded to a non-detection report plus an illustrative SED; if it passes, the conditional acceptance is justified. The reader's CONDITIONAL verdict is appropriate, and no verdict change is needed.","tokens_in":15505,"tokens_out":8136,"duration_ms":80666,"concrete_test":"Extract from the public LHAASO data the quasi-quiet light curve restricted to the exact time intervals of VERITAS quasi-quiet observations (MJD 59280–59449 and 59589–59697). Compute the exposure-weighted mean flux and its uncertainty in these windows and compare it with the full quasi-quiet average from Cao et al. (2024b) used for the LHAASO spectrum. If the two agree within ~1σ, the Γ<3.13 constraint stands; if they differ by a factor ≳1.5, the VERITAS upper limit is not representative of the quasi-quiet state and the constraint must be re-derived (or relaxed) using a time-matched LHAASO spectrum.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In §5 the paper states: 'Under the assumption that the VERITAS observations probe the average flux during this state, the power-law index cannot exceed 3.13 at 95% confidence.' This assumption is load-bearing because the LHAASO quasi-quiet state spans MJD 59278–60248, while the VERITAS quasi-quiet exposure (6.53 h) is confined to MJD 59280–59697: it covers the pre-active period and only ~3.5 months after the active window, not the full ~2.7-year quasi-quiet baseline. If the source's quiet-state emission decays, flares, or otherwise varies on month-to-year timescales (the active/quiet contrast is already a factor ~4), the mean flux during the VERITAS windows need not equal the time-averaged flux used to define the LHAASO quasi-quiet spectrum. An upper limit constrains the average flux over the VERITAS exposure only; if the true state average is higher (e.g., due to a flare missed by VERITAS), the derived Γ<3.13 would be artificially strong. The authors flag the assumption honestly, but it remains the central empirical claim's weakest link. The SED-peak range (100 GeV–2 TeV) is broad and based on the active period, and the corona model is explicitly illustrative, so the index constraint is the most consequential new result at risk.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports VERITAS observations of the LLAGN NGC 4278 that are contemporaneous with the LHAASO detection of 1LHAASO J1219+2915. NGC 4278 is not detected; 95% CL flux upper limits are derived for active and quasi-quiet states (Table 1). Combining the quasi-quiet VERITAS upper limit with LHAASO and Fermi-LAT data, the authors infer that the power-law index of the quasi-quiet VHE spectrum must be Gamma < 3.13 (if VERITAS samples the state average), and that the SED peak lies between 100 GeV and 2 TeV. They then present an illustrative hadronic corona model with p-gamma interactions that can simultaneously produce hard X-rays and TeV gamma rays, with neutrinos slightly below IceCube sensitivity.","tokens_in":15898,"tokens_out":4769,"duration_ms":49794,"significance":"The observational core is standard and reliable: the analysis uses reflected-region background estimation, Li&Ma significances, and is cross-checked with an independent analysis package. The upper limits at ~280 GeV are more constraining than Fermi-LAT for this source and are a useful new data point. The authors are careful to label their corona model as illustrative and to state the main assumptions. If the temporal averaging assumption is justified, the Gamma < 3.13 constraint and the SED peak bracket would be the first meaningful VHE spectral constraints for a quasi-quiet LLAGN and would favor models in which the VHE component peaks below the LHAASO band.","major_comments":[{"comment":"Section 5 states 'Under the assumption that the VERITAS observations probe the average flux during this state, the power-law index cannot exceed 3.13...' This assumption is load-bearing: the VERITAS quasi-quiet exposure (6.53 h over MJD 59280–59697) covers only the pre-active part and ~3.5 months after the active window, while the LHAASO quasi-quiet state spans MJD 59278–60248. The source shows factor ~4 active/quiet contrast, so unresolved variability within the quiet state could make the VERITAS sample unrepresentative of the time-averaged LHAASO spectrum. The constraint may be artificially strong if the state-average flux is higher. Please either justify representativeness with LHAASO light-curve data restricted to the VERITAS windows or explicitly downgrade the claim to a time-averaged upper limit that is not directly comparable to the LHAASO state average.","section":"§5 (Spectral Energy Distribution)"},{"comment":"The log-parabolic fit treats the VERITAS upper limit by 'assuming that it represents a flux of zero, with an uncertainty that is half of the flux upper limit.' This ad hoc replacement can bias the fitted curvature and hence the reported 100 GeV–2 TeV peak range. It is also unclear whether the active-period VERITAS limit (Table 1, 1e-7) or quasi-quiet limit (8.3e-9) was used. A proper likelihood (e.g., including the upper limit as a censored data point) and a sensitivity test should be reported before the peak bracket is used in the conclusions.","section":"§5 log-parabolic fit"},{"comment":"The text says Fermi-LAT 'provides spectral fluxes even when they are of low significance, and these were included in the fit,' yet Table 2 reports only 95% upper limits. If low-significance flux points are used as measurements, the fit may be biased; if only upper limits are used, the procedure should be stated. Please clarify and, if fluxes are used, report their values and significances.","section":"§5, Table 2"}],"minor_comments":[{"comment":"The abstract's statement that 'The flux upper limits constrain the photon spectrum of the quasi-quiet period' omits the temporal-representativeness caveat made explicit in §5. The abstract should carry the same conditionality as the body text.","section":"Abstract"},{"comment":"Typo: 'should not not be limited by the infall time' should be 'should not be limited by the infall time.' Also 'implicitely' should be 'implicitly.'","section":"§6.1"},{"comment":"Consider adding a scale bar or a clear angular-size annotation; the point spread function is shown, but the source separation from the 3.1σ fluctuation at 0.15° is easier to judge with an explicit angular scale.","section":"Figure 1"}],"recommendation":"major_revision","confidential_remarks":"The empirical upper limits and the honest presentation of the corona model are the paper's main strengths. The central new claim, Gamma < 3.13, rests on the unsupported assumption that 6.53 h of VERITAS data average over the quasi-quiet state. I believe the authors can address this with public LHAASO light-curve data restricted to the VERITAS time windows; if not, the claim should be softened. The log-parabolic fit's ad hoc use of the upper limit also needs revision before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a clean, honest observational paper. VERITAS serendipitously covered NGC 4278 during the LHAASO epoch, saw nothing, and the resulting upper limits actually do some work: they push the quasi-quiet VHE index below the LHAASO value and, together with Fermi-LAT, bracket a spectral peak somewhere in the 100 GeV–2 TeV range. That is a real, citable result.\n\nThe analysis itself is standard and reproducible: Li&Ma significance, cross-check with an independent package, effective areas per pointing offset. The authors also deserve credit for saying plainly that the log-parabolic fit treats the VERITAS upper limit as a zero-flux point, and that the corona model is an example, not a prediction. That is more transparency than most papers in this field.\n\nNow the soft spots, in order of seriousness.\n\nFirst, the Γ < 3.13 constraint is conditional on the VERITAS quasi-quiet exposure (6.53 h scattered over about a year) sampling the time-averaged LHAASO quasi-quiet state, which spans 2.7 years. The paper states this assumption explicitly, so the reader is warned, but it is still load-bearing: if the quiet-state flux drifts or flares on timescales longer than the VERITAS windows, the limit constrains only the average over those windows. I don't think this kills the result—the source's contrast between active and quiet is only a factor of ~4, and the VERITAS windows do bracket the active period—but it does mean the constraint is not as clean as the abstract implies.\n\nSecond, the SED peak bracket (100 GeV–2 TeV) is broad, and it comes from a fit where the VERITAS upper limit is injected as a pseudo-detection of zero. The authors are honest about this, and the bracket is what it is, but it would be easy for a reader to over-interpret 'peak between 100 GeV and 2 TeV' as a sharper statement than the data support.\n\nThird, the corona model: it is a plausible interpretation, but the proton spectrum is chosen by hand to reproduce the same SED it is meant to explain, and the neutrino expectation is normalized to the gamma-ray flux rather than derived from first principles. The paper labels this correctly, so I don't count it as a flaw, just as a limitation on what the model section adds.\n\nMinor: the SED figure includes a flux point from a private communication (Bronzini 2024). The authors explain why they don't use it, but citing a private comm in a published figure is awkward; that point should either be removed or the analysis should be public by the time this is refereed.\n\nBottom line: this is a solid, moderately significant observational contribution with an illustrative model attached. It deserves a serious referee. The main things a referee should probe are the sampling assumption behind Γ < 3.13 and whether the log-parabolic peak bracket is robust to the zero-flux pseudo-point treatment. Neither is a deal-breaker.","headline":"Solid VERITAS upper limits on NGC 4278, but the headline index constraint rests on an unverified sampling assumption and the corona model is explicitly illustrative—worth a careful referee, not a desk reject.","tokens_in":16859,"tokens_out":2348,"would_cite":true,"duration_ms":24782,"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":"Serendipitous VERITAS upper limits on the low-luminosity AGN NGC 4278 place its very-high-energy gamma-ray peak between 100 GeV and 2 TeV and make a coronal proton model a viable explanation.","keywords":["gamma-ray astronomy","low-luminosity active galactic nuclei","NGC 4278","very-high-energy gamma rays","flux upper limits","hadronic corona model","spectral energy distribution","neutrino predictions"],"falsifier":"A longer, deeper observation with an imaging Cherenkov telescope that detects NGC 4278 in the quasi-quiet state at a flux above the reported upper limit would falsify the index constraint; alternatively, if a future SED fit places the gamma-ray peak outside 100 GeV–2 TeV, or if a stacked search of similar low-luminosity AGN finds no neutrino signal matching the hadronic prediction, the coronal p-gamma interpretation would need revision.","tokens_in":15365,"feed_emoji":"🔭","tokens_out":4182,"duration_ms":39114,"temperature":0.7,"pith_summary":"The paper argues that the nearby low-luminosity active galactic nucleus NGC 4278, recently detected in TeV gamma rays by LHAASO, has its gamma-ray peak at a few hundred GeV rather than in the 1–20 TeV band. Using serendipitous VERITAS observations that overlap the LHAASO campaign, the authors derive a 95% confidence upper limit on the quasi-quiet flux that, combined with Fermi-LAT upper limits, constrains the power-law index to below 3.13 and locates the SED peak anywhere between 100 GeV and 2 TeV. They then show that a hadronic corona model, in which protons accelerated in the AGN corona interact with X-ray photons to produce pions, can naturally explain both the hard X-rays and the TeV gamma rays without violating energetics. If correct, this identifies a new mechanism for transient TeV emission from low-luminosity AGN and predicts a neutrino signal slightly below IceCube's current sensitivity.","feed_headline":"VERITAS upper limits place NGC 4278's gamma-ray peak at 0.1–2 TeV","feed_subtitle":"A quiet low-luminosity AGN's TeV emission is best explained by protons in its corona, not a jet.","key_machinery":"The argument rests on two pieces. First, the VERITAS flux upper limit of 8.3e-9 TeV^-1 m^-2 s^-1 at 680 GeV in the quasi-quiet state, obtained from 6.53 hours of serendipitous observations, which is used in a log-parabolic fit to the LHAASO and Fermi-LAT data to bracket the SED peak. Second, the analytic corona model: protons accelerated to ~20 TeV interact with ~10 keV synchrotron X-ray photons through the Delta resonance, producing pions that decay into gamma rays, neutrinos, and electron-positron pairs; the pairs cool by synchrotron emission and re-inject hard X-rays, providing a self-amplifying target photon field.","core_discovery":"The central claim is that the quasi-quiet very-high-energy spectrum of NGC 4278 is harder than a simple extrapolation of LHAASO's power law allows, with the spectral index constrained to Gamma < 3.13 at 95% confidence, and that the peak of the gamma-ray spectral energy distribution lies between 100 GeV and 2 TeV. The paper further claims that this peak can be produced by photohadronic interactions of ~20 TeV protons with ~10 keV X-ray photons in the AGN corona, without the need for relativistic jet boosting, and that the secondary electrons from pion decay naturally produce the observed hard X-ray component.","pith_inferences":["The same coronal p-gamma mechanism might apply to other low-luminosity AGN with hard X-ray peaks, making the class a population of faint, transient TeV emitters worth targeting in future surveys.","If the quasi-quiet state is not actually steady at the level probed by the 6.53 hours of VERITAS exposure, the index constraint and peak location could shift; a dedicated monitoring campaign could test this.","The paper's assumption that the two LHAASO quasi-quiet points share the active-state curvature is fragile; a future observation with more quasi-quiet flux points could confirm or alter the peak bracket.","The model's implication that the neutrino flux at ~100 TeV is comparable to that of NGC 1068 suggests that high-energy neutrino telescopes might distinguish hadronic corona models from leptonic jet models by looking for a spectral crossover."],"forward_implications":["The VHE gamma-ray peak of NGC 4278 lies in a band where both VERITAS and future Cherenkov telescopes are sensitive, so a firm detection or deeper upper limit can directly test the predicted index constraint.","The hadronic corona model predicts a neutrino spectrum that, though below IceCube's steady sensitivity, could be probed with stacked low-luminosity AGN samples or future neutrino telescopes.","If the peak is genuinely at a few hundred GeV, LHAASO's detection in the 1–20 TeV band requires a relatively hard tail, implying a specific high-energy cutoff in the proton spectrum.","The model provides a self-consistent explanation for the historically puzzling hard X-ray spectrum of NGC 4278 as synchrotron radiation from secondary pairs rather than intrinsic jet emission."],"fun_headline_variants":["VERITAS non-detection pins NGC 4278's TeV peak to 0.1–2 TeV","Corona protons, not jet, explain NGC 4278's TeV peak","Neutrino flux from NGC 4278's corona just below IceCube's reach","No VERITAS detection, but NGC 4278's TeV peak is now pinpointed"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The 6.53 hours of VERITAS exposure spread over roughly a year are assumed to represent the average flux of the entire quasi-quiet state, so that the upper limit can be combined with the time-averaged LHAASO spectrum to constrain the spectral index.","fun_headline_variants_meta":{"raw":{"variants":["VERITAS non-detection pins NGC 4278's TeV peak to 0.1–2 TeV","Corona protons, not jet, explain NGC 4278's TeV peak","Neutrino flux from NGC 4278's corona just below IceCube's reach","No VERITAS detection, but NGC 4278's TeV peak is now pinpointed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00119,"raw_usage":{"total_tokens":4801,"prompt_tokens":851,"completion_tokens":3950,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":3863}},"tokens_in":595,"tokens_out":3950,"duration_ms":26348,"temperature":1.0,"reasoning_tokens":3863,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T05:46:53.302925+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A longer, deeper observation with an imaging Cherenkov telescope that detects NGC 4278 in the quasi-quiet state at a flux above the reported upper limit would falsify the index constraint; alternatively, if a future SED fit places the gamma-ray peak outside 100 GeV–2 TeV, or if a stacked search of similar low-luminosity AGN finds no neutrino signal matching the hadronic prediction, the coronal p-gamma interpretation would need revision.","supporting_citations":[],"review_version":1}