{"id":"e835430a-89b2-4ee1-9afe-9fb006aa693c","arxiv_id":"2501.09108","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A single-zone model with added proton-proton collisions suggests W Comae and 1ES 1959+650 need hadronic emission, and CTAO sensitivity forecasts show all four BL Lacs are detectable.","lead":"This paper models four BL Lac blazars as potential sources of cosmic rays, neutrinos, and gamma rays, using a single-zone lepto-hadronic model with added proton-proton interactions. It then forecasts that the future CTAO observatory could detect gamma-ray emission from all four sources, and that two of them may require hadronic processes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that W Comae and 1ES 1959+650 require hadronic emission rests on adopting the maximum allowed proton luminosity from Rodrigues et al. (2024) without a statistical comparison; lowering Lp or refitting the leptonic component could remove the need.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing issue: the analysis uses maximum proton luminosities from Rodrigues et al. (2024) as if they were measured model parameters, and it lacks a statistical comparison between leptonic and hadronic models. My stress-test agrees and adds specificity: the comparison is visual, the electron parameters are also inherited, and the CTAO detectability analysis in Section 3 is independent of this hadronic claim because it is based on observed gamma-ray fluxes. Since the reader already issued a CONDITIONAL verdict requiring clarification of Lp and a goodness-of-fit assessment, my read does not change the verdict; the paper should be revised to fit Lp (or label the predictions as upper limits) and to provide a statistical test of the hadronic component before the claim can be accepted.","tokens_in":23976,"tokens_out":5597,"duration_ms":57774,"concrete_test":"For W Comae and 1ES 1959+650, fit the multiwavelength SED with a purely leptonic SSC model (electron luminosity, B, and possibly a second zone as free parameters) using the same data set and a chi-square or likelihood statistic. Then add protons with Lp as a free parameter and compare via AIC/BIC. If the leptonic-only best fit is acceptable (e.g., p-value > 0.01) or the delta-AIC for adding hadrons is < 10, the paper's claim that a leptonic model is insufficient fails. As a second check, rerun the published model with Lp = 0.1 times its adopted maximum and verify whether the TeV excess above the leptonic component remains significant.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central conclusion (§4) that a purely leptonic model is insufficient for W Comae and 1ES 1959+650, and that these sources are therefore effective cosmic-ray emitters, rests on the proton luminosity values adopted in §2. Specifically, Table 1 lists log10 Lmax_p with upper-limit symbols ('<') for several sources, and the text imposes only the constraint Lp/Le < 10^3. The model is evidently run with these maximum allowed proton luminosities, so the hadronic gamma-ray and neutrino fluxes are upper limits, not values inferred from the SED. The evidence for hadronic emission is a visual comparison in Figs. 1–2 between the 'Primary e− Syn+IC' curve and VERITAS/HESS data, with no chi-square, likelihood, or information criterion reported for the SED fits. Moreover, the electron luminosity and other parameters are also inherited from Rodrigues et al. (2024) and are not refit here; a purely leptonic model with different parameters (e.g., higher electron luminosity, different B, or an additional zone) might account for the TeV points. Thus the claimed need for hadronic processes is underdetermined by the analysis as presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the open-source AM3 lepto-hadronic one-zone code to four BL Lacs (W Comae, 1ES 1959+650, PKS 2005-489, PKS 2155-304), adopting model parameters from Rodrigues et al. (2024), extending the framework to include pp interactions, and simulating CTAO detectability with Gammapy and prod5 instrument response functions. The authors conclude in Section 4 that a purely leptonic model is insufficient for W Comae and 1ES 1959+650, that hadronic processes contribute substantially at high energies, and that all four sources could be effective cosmic-ray emitters detectable by CTAO in 50-hour exposures.","tokens_in":24248,"tokens_out":5429,"duration_ms":51652,"significance":"The CTAO visibility and sensitivity analysis in Section 3 is concrete and reproducible: it uses public IRFs, Gammapy, specific array configurations, and reports AIC-based model selection with fit parameters in Tables 6 and 7. The choice to use the open-source AM3 code and to state model parameters explicitly is a strength. However, the paper's central astrophysical claim, that hadronic emission is required for W Comae and 1ES 1959+650, is not established by the analysis as presented. The hadronic component is normalized to the maximum proton luminosity allowed by the prior constraint, and the comparison with TeV data is visual only. If reframed as an upper-limit study, the model predictions and CTAO forecasts would be a useful contribution, but the current wording overreaches the evidence.","major_comments":[{"comment":"The hadronic normalization is taken at the upper limit, not determined by a fit. Table 1 lists log10 Lmax_p with '<' symbols for all four sources, and Section 2 states only the constraint Lp/Le < 10^3. The hadronic gamma-ray and neutrino fluxes in Figs. 1–4 therefore scale linearly with a parameter that is set to its maximum allowed value, so these curves are upper limits, not measured or fitted components. Section 4's conclusion that a purely leptonic model is insufficient for W Comae and 1ES 1959+650 depends entirely on this maximal normalization. The authors should either fit Lp as a free parameter and report uncertainties, or explicitly relabel all hadronic curves and neutrino fluxes as upper limits and soften the conclusion accordingly.","section":"Section 2, Table 1"},{"comment":"The evidence against the purely leptonic model is a visual comparison. No goodness-of-fit statistic (chi-square, likelihood, or AIC) is computed for the SED fits in Section 2; AIC appears only in Section 3.1 for the CTAO source-model selection. Visual insufficiency is especially fragile here because the leptonic parameters (electron luminosity, B, R, Gamma) are inherited from Rodrigues et al. (2024) rather than refit in this work, and a different leptonic parameter set (for example a higher electron luminosity or a different magnetic field) might describe the VERITAS and HESS points without hadronic contributions. The claim that hadronic interactions 'significantly improve the fit' requires a quantitative model comparison.","section":"Section 2, Figs. 1–2"},{"comment":"The CTAO detectability forecasts are based on the log-parabola spectrum fitted to Fermi and VTSCat data in Section 3.1 (Table 7), not on the lepto-hadronic SED from Section 2. Consequently, the statement in Section 4 that CTAO observations will test these sources as cosmic-ray emitters is not directly connected to the hadronic model: the simulated detection would occur even for a purely leptonic source with the same observed spectrum. The authors should clarify which model is being simulated and state explicitly what a CTAO detection or non-detection would imply for the proton luminosity and the hadronic scenario.","section":"Section 3.2, Figs. 6, 8, 10, 12"}],"minor_comments":[{"comment":"The definition of deltaAICm in the table caption uses BICm (1 - BICm/AICPL) while the text and header refer to AIC; please correct the notation and also align the label 'BPL' with Equation (4), which is called the exponential cutoff power law (ECPL).","section":"Section 3.1, Table 6"},{"comment":"The phrase 'the ashes indicating the source is not visible' should be 'the dashes indicating the source is not visible'.","section":"Section 3.2, Table 11 caption"},{"comment":"The caption references 'see Table 5' for the 1ES 1959+650 counterparts, but the correct table is Table 3.","section":"Figure 8 caption"},{"comment":"The sentence 'Figure 7 illustrates the energy spectra for the W Comae region' should refer to the 1ES 1959+650 region.","section":"Section 3.2, 1ES 1959+650 paragraph"},{"comment":"The caption contains a duplicated sentence: 'The black curve represents the total cumulative SED... The black curve depicts the total cumulative SED...' Please remove the duplicate.","section":"Figure 1 caption"},{"comment":"The text refers to 'AM31 software' but the cited reference [67] and the project name are 'AM3'; please fix this typo.","section":"Section 2, AM3 description"},{"comment":"Several declaration items ('Ethics approval and consent to participate', 'Consent for publication', 'Materials availability') are present but empty; please either complete them or remove the placeholders.","section":"Declarations"},{"comment":"The title 'Blazars Jets' should be 'Blazar Jets' for grammatical correctness.","section":"Title and abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper's main risk is the gap between the maximal-Lp modeling in Section 2 and the strong conclusion in Section 4. The CTAO forecast section is the strongest part and could form the core of a revised paper; the hadronic-necessity claim should be reframed as an upper-limit or prospect study unless a proper fit is performed. The manuscript also contains several labeling errors that suggest a careful proofreading pass."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the CTAO forecast part is a solid, reproducible piece of work; the hadronic-emission claim is not yet supported by the analysis.\n\nThe genuinely new pieces are the addition of pp interactions in AM3 to the Rodrigues et al. (2024) model for these four BL Lacs, and a careful CTAO visibility/sensitivity study using prod5 IRFs and Gammapy. That gives observers a concrete target list and exposure recommendation, which is useful. The use of AIC for the catalog spectral models is also fine.\n\nThe problem is Section 4. The claim that W Comae and 1ES 1959+650 require hadronic emission rests on two things: a visual comparison of the electron synchrotron+SSC curve with VERITAS points, and a proton luminosity taken at its maximum from the prior Rodrigues fit. Table 1 gives Lmax_p as upper limits, and the only constraint stated is Lp/Le < 1e3. Running at that ceiling means the hadronic gamma-ray and neutrino fluxes are an upper-limit scenario, not an inference from the SED. No chi-square, likelihood, or AIC is reported for the SED fits, and the electron-side parameters are also inherited. A lower Lp or a refitted electron component could plausibly do the job. I would also want to see the predicted neutrino fluxes checked against IceCube stacking limits for these sources.\n\nMinor issues: duplicated caption text and a mislabeled figure (the 1ES 1959+650 sensitivity plot is described as W Comae). Easy to fix, but sloppy.\n\nBottom line: the paper deserves a serious referee. The CTAO section alone is worth publishing as a target-selection study. Ask for fit statistics, an explicit rationale for Lp, a lower-Lp test or leptonic refit, and an IceCube consistency check; conditional acceptance after that.","headline":"Solid CTAO forecast, but the hadronic claim rests on upper-limit proton luminosity and no fit statistics; worth a serious referee.","tokens_in":24804,"tokens_out":5599,"would_cite":false,"duration_ms":52934,"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":"W Comae and 1ES 1959+650 show TeV gamma-ray excesses that electron-only emission cannot explain, pointing to proton interactions and making all four BL Lacs promising neutrino and cosmic-ray sources.","keywords":["blazars","cosmic rays","gamma rays","neutrinos","lepto-hadronic model","BL Lacertae objects","CTAO","multimessenger astronomy"],"falsifier":"A 50-hour CTAO observation of W Comae or 1ES 1959+650 that finds the TeV spectrum matches the electron-only prediction with no excess around $10^{12}$–$10^{14}$ eV, or a stacked neutrino-telescope search over these four directions that sees no events at the modeled all-flavor fluxes near $10^{14}$–$10^{16}$ eV, would falsify the central claim.","tokens_in":23757,"feed_emoji":"🔭","tokens_out":11230,"duration_ms":106970,"temperature":0.7,"pith_summary":"This paper argues that four BL Lacertae blazars — W Comae, 1ES 1959+650, PKS 2005-489, and PKS 2155-304 — can account for their multiwavelength emission only if their jets contain both electrons and protons in a single compact blob. For W Comae and 1ES 1959+650, the electron-only component under-predicts the observed TeV gamma-ray flux, so the fit needs pion production from proton-photon and proton-proton collisions; the other two sources are mostly leptonic but show hints of hadronic emission at the highest energies. If correct, these blazar jets accelerate protons to very high energies and could be sources of cosmic rays and of neutrinos peaking near $10^{14}$–$10^{16}$ eV. The paper then shows that the Cherenkov Telescope Array Observatory should detect all four sources in 50 hours, which would put this hadronic picture to a direct test.","feed_headline":"Two blazars need hadronic beams to explain their TeV light","feed_subtitle":"A lepto-hadronic fit says W Comae and 1ES 1959+650 emit proton-driven gamma rays; CTAO can test it in 50 hours.","key_machinery":"The load-bearing machinery is the one-zone lepto-hadronic radiation model, in which all radiation comes from a single spherical emission region containing power-law populations of electrons and protons. The calculation solves time-dependent coupled kinetic equations for electrons, positrons, protons, neutrons, photons, neutrinos, pions, and muons, combining primary electron synchrotron and synchrotron self-Compton emission with proton synchrotron and inverse-Compton emission, photon-pair production, and neutral-pion decay from both proton-photon and proton-proton interactions; the proton-proton channel is added on top of the earlier framework. On the observability side, the CTAO projection uses a simultaneous likelihood fit over multi-GeV-to-TeV catalog data, an extragalactic-background-light absorption model, and a 1D ON/OFF analysis with 50-hour instrument response functions for the North and South arrays at several zenith angles.","core_discovery":"The central claim is that a single-zone lepto-hadronic model — a spherical blob with radius $R \\sim 10^{16.7}$–$10^{17.3}$ cm, magnetic field $B \\sim 0.7$–$1.4$ G, bulk Lorentz factor $\\Gamma \\sim 4$–$12$, and equal power-law indices for electrons and protons — reproduces the gamma-ray-to-radio spectral energy distributions of all four BL Lacs. For W Comae and 1ES 1959+650, the high-energy tail around $10^{10}$–$10^{14}$ eV cannot be matched by electron synchrotron and inverse Compton emission alone; the fit requires hadronic channels, chiefly neutral-pion decay from proton-photon ($p\\gamma \\to \\pi^0 \\to \\gamma\\gamma$) and proton-proton ($pp \\to \\pi^0 \\to \\gamma\\gamma$) interactions. The corresponding all-flavor neutrino fluxes peak near $10^{14}$–$10^{16}$ eV, with maximum values around $10^{-14.4}$ to $10^{-11.3}$ erg cm$^{-2}$ s$^{-1}$ for the four sources. The paper concludes that these objects could be effective cosmic-ray emitters and that CTAO observations will be able to test this interpretation.","pith_inferences":["A testable extension is to compare the predicted neutrino fluxes directly with stacked neutrino-telescope limits: a non-detection at the modeled fluxes would force the proton luminosity down and weaken the hadronic interpretation.","The TeV excess in W Comae and 1ES 1959+650 could also be explained if those data points come from a separate flare zone with different physical conditions; a time-resolved, multi-zone fit would settle whether the single blob is truly hadronic.","The same fitting procedure could be applied to other gamma-ray-selected blazars, prioritizing objects whose existing TeV data already lie above the electron-only envelope; those would be the best candidates for the next generation of gamma-ray and neutrino instruments."],"forward_implications":["W Comae and 1ES 1959+650 would have a substantial hadronic component in their TeV emission, making them the two most promising neutrino targets among the four.","CTAO should detect all four blazars in 50-hour exposures, with 1ES 1959+650's modeled spectrum extending to about 30 TeV.","The predicted all-flavor neutrino fluxes near $10^{14}$–$10^{16}$ eV give neutrino telescopes a concrete flux level to test with point-source or stacked searches.","A positive detection at multi-TeV energies would strengthen the case that BL Lac jets accelerate nuclei and contribute to the cosmic-ray flux, rather than emitting only leptonic radiation."],"supporting_citations":[{"why":"supplies the base lepto-hadronic one-zone model, the four source parameters, and the maximum neutrino flux values this paper extends with proton-proton interactions.","marker":"[24]"},{"why":"provides the multiwavelength data compilation that the modeled spectral energy distributions are fitted against.","marker":"[33]"},{"why":"the open-source code that solves the time-dependent coupled equations for all particle species in the radiation model.","marker":"[67]"},{"why":"the hadronic interaction generator used to compute pion production in the added proton-proton channel.","marker":"[69]"},{"why":"supplies the very-high-energy gamma-ray measurements of W Comae that drive the conclusion that a purely leptonic model is insufficient.","marker":"[70]"},{"why":"provides the extragalactic-background-light optical depth used to convert intrinsic spectra into observed spectra.","marker":"[87]"},{"why":"supplies the simultaneous-likelihood and ON/OFF sensitivity methodology used for the CTAO detection projections.","marker":"[74]"},{"why":"the gamma-ray analysis software used to run the likelihood fits and sensitivity calculations.","marker":"[76]"},{"why":"provides the CTAO instrument response functions and sensitivity thresholds assumed for the 50-hour observations.","marker":"[90]"}],"fun_headline_variants":["Two blazars need hadronic beams to explain their TeV light","Proton-driven gamma rays fit TeV tails in two blazars","CTAO can test hadronic model for W Comae and 1ES 1959+650","Blazar TeV emission hints proton acceleration, CTAO to verify"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The hadronic gamma-ray and neutrino fluxes are set by the maximum proton luminosity, which the paper adopts as an upper-limit parameter from the earlier model rather than as a measured quantity; if the true proton luminosity is lower, or if the TeV emission comes from a separate flare zone, the claimed need for hadronic processes disappears.","fun_headline_variants_meta":{"raw":{"variants":["Two blazars need hadronic beams to explain their TeV light","Proton-driven gamma rays fit TeV tails in two blazars","CTAO can test hadronic model for W Comae and 1ES 1959+650","Blazar TeV emission hints proton acceleration, CTAO to verify"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1389,"prompt_tokens":985,"completion_tokens":404,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":601,"completion_tokens_details":{"reasoning_tokens":319}},"tokens_in":601,"tokens_out":404,"duration_ms":4690,"temperature":1.0,"reasoning_tokens":319,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:10:58.739633+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A 50-hour CTAO observation of W Comae or 1ES 1959+650 that finds the TeV spectrum matches the electron-only prediction with no excess around $10^{12}$–$10^{14}$ eV, or a stacked neutrino-telescope search over these four directions that sees no events at the modeled all-flavor fluxes near $10^{14}$–$10^{16}$ eV, would falsify the central claim.","supporting_citations":[{"cited_title":"A gamma-ray study of galactic PeVatron candidates LHAASO J1825-1326 and LHAASO J1839-0545","cited_arxiv_id":null,"evidence_quote":"supplies the simultaneous-likelihood and ON/OFF sensitivity methodology used for the CTAO detection projections."},{"cited_title":"Gammapy - A prototype for the CTA science tools; 2017","cited_arxiv_id":null,"evidence_quote":"the gamma-ray analysis software used to run the likelihood fits and sensitivity calculations."}],"review_version":1}