{"id":"79271439-d74d-4ae2-bcd9-ec0d0193023e","arxiv_id":"2505.04160","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"Two z>1 FSRQ blazars that point at IceCube neutrino events are modeled as efficient PeV neutrino emitters, with muon neutrino detection probabilities of about 2% and 0.8%.","lead":"This paper models two very distant blazars, each linked to an IceCube neutrino alert, with a jet model that includes protons, and computes the chance that either source produces a detectable neutrino during its gamma-ray flare. It finds probabilities of about 2% and 0.8%, and argues that despite the low odds these are efficient PeV neutrino emitters, with X-rays the better electromagnetic tracer.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The BLR-location assumption (r'_in ≈ r'_BLR) is the load-bearing pillar: outside the BLR, Figure 6 shows fpγ drops by orders of magnitude, so the 2%/0.8% neutrino probabilities and the 'efficient PeV emitter' claim collapse; a re-fit with r'_in outside the BLR would test this.","rationale":"The reader's conditional verdict is well aligned with my reading. The paper is transparent about its weaknesses (by-eye fitting, non-simultaneous data, assumed NVSS X-ray luminosity, emission-duration ambiguity), and the strongest quantitative claims are honestly presented as scenario-dependent estimates. The single most decisive vulnerability is the r'_in ≈ r'_BLR assumption: it is the reason BLR photons dominate pγ interactions, and it is the reason the sources can be called efficient PeV neutrino emitters with moderate baryon loading. The paper itself provides the counter-evidence: short-term variability alone is insufficient, and known FSRQs (PKS B1424-418, PKS 1502+106) have their emission regions argued to be outside the BLR. Figure 6 makes the quantitative consequence explicit: outside the BLR, fpγ drops by orders of magnitude, so the neutrino luminosity at fixed proton luminosity collapses and the quoted detection probabilities become much smaller. Since the SED fit is by-eye and non-simultaneous, it cannot uniquely determine r'_in. Thus the central claim is conditional on a location that the data do not force. This does not invalidate the paper; it means the probabilities should be read as scenario-dependent estimates rather than measurements. The proposed re-fit with r'_in outside the BLR directly tests whether the assumption is required. I agree with the reader's weakest_assumption and see no reason to change the CONDITIONAL verdict.","tokens_in":30990,"tokens_out":4789,"duration_ms":47811,"concrete_test":"Re-fit the GB6 J2113+1121 SED with r'_in fixed outside the BLR (e.g., at the dust torus radius), re-optimizing the nine free parameters (R, B, δ, Le,inj, γe,break, γe,max, ne,1, ne,2, Lp,inj) against the same multi-band data, and recompute the IceCube muon-neutrino count from Eq. 23. If a comparably acceptable by-eye SED exists with PDet dropping below ~0.1% (or requiring L_p,inj above the Eddington limit by multiple orders), the BLR-location assumption is confirmed as load-bearing and the 2%/0.8% probabilities should be weakened; if no such fit exists, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is the placement of the emission region at the BLR boundary, r'_in ≈ r'_BLR (Section 2.3, Table 1). This is what makes BLR photons the dominant pγ target and produces the claimed 'efficient PeV neutrino emitter' status. Section 4 openly states that short-term variability alone does not prove this location, and that other FSRQs (PKS B1424-418, PKS 1502+106) are argued to emit beyond the BLR. Figure 6 then shows that if the blob sits outside the BLR, fpγ falls by orders of magnitude over the relevant proton energies. With the same proton luminosity, the predicted neutrino fluxes, the 2% and 0.8% IceCube detection probabilities, and the cascade-dominated X-ray/GeV components all drop correspondingly; the neutrino-to-gamma luminosity ratio Yνγ ≈ 0.5–0.7 would no longer follow. The by-eye, non-simultaneous SED fit (Section 2.3) does not uniquely constrain r'_in, and for NVSS the cascade flux that limits the neutrino flux is set by an assumed X-ray luminosity L_X ≈ 2×10^46 erg/s rather than by data. The quoted probabilities are therefore scenario-dependent predictions, not robust measurements.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies a one-zone leptohadronic model with a BLR external photon field to two high-redshift FSRQs, GB6 J2113+1121 and NVSS J171822+423948, which have been temporally and spatially associated with IceCube neutrino events. From by-eye SED fits during their gamma-ray flares, the authors derive neutrino spectra, IceCube single-muon detection probabilities of about 2% and 0.8%, baryon loading factors of a few hundred, and all-flavor neutrino-to-gamma luminosity ratios Y_nu_gamma ~ 0.5-0.7. They conclude that the sources are efficient PeV neutrino emitters and that hadronic cascade synchrotron emission dominates the X-ray and GeV bands. The paper also compares these sources with other neutrino candidate blazars and discusses how the location of the emission region relative to the BLR affects neutrino production efficiency.","tokens_in":31254,"tokens_out":4171,"duration_ms":41436,"significance":"If the emission region is indeed at the BLR boundary, the paper provides a concrete multimessenger interpretation for two distant FSRQ-neutrino associations and makes a falsifiable prediction that X-rays, rather than gamma-rays, are the better electromagnetic tracer of neutrino output. The modeling is transparent: the authors explicitly state that the fit is by eye, that the multi-band data are not simultaneous, that the NVSS X-ray luminosity is assumed rather than measured, and that the BLR-location assumption is not uniquely required by variability. These candor and the use of a standard, well-documented leptohadronic framework are strengths. The quantitative conclusions are nevertheless scenario-dependent because the neutrino normalization is set by the free proton luminosity and because the central claim hinges on the placement of the emission region at r'_in ~ r'_BLR.","major_comments":[{"comment":"The assumption r'_in ~ r'_BLR is load-bearing for the central claim. The authors note in Section 4 that short-term variability alone does not prove this location and that other FSRQs (PKS B1424-418, PKS 1502+106) are argued to emit beyond the BLR. Figure 6 shows that if the blob is outside the BLR, fp_gamma drops by orders of magnitude over the relevant proton energies; with the same proton luminosity, the neutrino fluxes, the 2% and 0.8% detection probabilities, and the cascade-dominated X-ray/GeV components would all drop correspondingly. A quantitative refit with r'_in outside the BLR, or at least a likelihood-based comparison of allowed locations, is needed before the conclusion that these sources are efficient PeV neutrino emitters can be considered robust.","section":"Section 4, Section 2.3, Table 1, Figure 6"},{"comment":"The fitting procedure is by-eye, with nine free parameters adjusted without a likelihood function and without reported uncertainties. The text itself states that the data are not simultaneous, 'which weakens the use of a steady-state solution.' Consequently, the quoted detection probabilities are point estimates whose statistical meaning is unclear. The authors should provide a sensitivity scan or alternative fits showing how P_Det and the cascade contributions vary within a plausible parameter range; as it stands, the abstract's phrase 'strongly suggests' is not commensurate with the fitting procedure.","section":"Section 2.3, Table 1"},{"comment":"For NVSS J171822+423948 there is no X-ray data, and the model instead assumes L_X ~ 2x10^46 erg/s. Since the text states that the X-ray flux 'strictly constrains the pair cascades emission, thereby limiting the neutrino flux,' the 0.8% detection probability is partly an input assumption rather than an output prediction. The authors should show how the neutrino flux and detection probability depend on the assumed X-ray luminosity or on the X-ray upper limit shown in Figure 4; otherwise the quoted number is not a robust prediction for this source.","section":"Section 3, Table 1, Figure 4"},{"comment":"The neutrino normalization is set by Lp,inj, which is a free parameter that the authors tune to 'maximize the neutrino flux while explaining the electromagnetic spectrum.' The resulting Y_nu_gamma ~ 0.5-0.7 and the detection probabilities are therefore scenario-maximized outputs of a tuned model, not parameter-free predictions. The paper should explicitly label them as such and soften the abstract's claim that the investigation 'strongly suggests these sources are efficient PeV neutrino emitters.'","section":"Section 2.3, Section 3, Eq. (27)"}],"minor_comments":[{"comment":"The text gives the energy of IceCube-201221A as 0.174 TeV, but the same event is described elsewhere as ~0.2 PeV; this appears to be a typo and should read 0.174 PeV.","section":"Section 3"},{"comment":"The caption lists 'GB6 J2113+1121' twice in the legend; presumably one entry should be 'NVSS J171822+423948.'","section":"Figure 8 caption"},{"comment":"The text says 'we consider eta = 103' where the intended value is eta = 10^3; this should be typeset consistently as 10^3 in the text and in the discussion of Eq. (22).","section":"Section 2.3, item 5"},{"comment":"The notation in Eq. (23) mixes the integration variable epsilon_nu^obs_mu with the differential d epsilon_nu^obs_mu; using a single symbol for the energy and the differential would improve readability.","section":"Equation (23)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the authors are unusually candid about their modeling limitations. The main issues are not internal inconsistencies but the gap between the conditional nature of the BLR-location assumption and the strength of the abstract's claim, plus the absence of sensitivity or uncertainty quantification. These are fixable with additional tests and reframing, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is concrete SED and neutrino modeling for GB6 J2113+1121 and NVSS J171822+423948, the two z>1 FSRQs with reported IceCube associations. The paper applies the standard one-zone leptohadronic machinery, so the novelty is the application, not the method. It does that carefully: the model is described in enough detail to be transparent, and the authors disclose the by-eye fitting, the non-simultaneous multi-band data, the assumed NVSS X-ray luminosity, and the emission-duration ambiguity. I credit them for that honesty.\n\nThe main quantitative claims—2% and 0.8% single-muon-neutrino detection probabilities, baryon loading ~300–400, Y_nu_gamma ~0.5–0.7—are not parameter-free predictions. The proton luminosity Lp,inj is a free parameter tuned to maximize neutrinos while matching the SED, so the detection probabilities are outputs of a tuned model. Worse, the whole neutrino efficiency story relies on placing the emission region at the BLR boundary (r'_in ~ r'_BLR). The paper itself acknowledges that short-term variability doesn't prove that location and that several FSRQs are argued to emit beyond the BLR, and Figure 6 shows that moving outside the BLR drops fpγ by orders of magnitude. At that point the 'efficient PeV emitter' claim collapses under the same proton luminosity. So the stress-test note is right: the BLR-location assumption is load-bearing.\n\nThat said, the soft spots are mostly disclosed by the authors rather than hidden. The by-eye fit and lack of uncertainties are real weaknesses, but they are stated plainly. The NVSS X-ray luminosity is assumed, not measured—that is a genuine limitation that the authors flag. The non-simultaneous data weakens the steady-state solution, as they admit.\n\nThe paper also does something useful beyond the headline numbers: it shows BLR photon absorption creates a gamma-ray cutoff near 5–10 GeV and argues that X-rays are a better neutrino tracer. That is a testable, physically sensible consequence. The comparison with other neutrino-associated blazars and the jet-power/Eddington discussion are reasonable and well-grounded in prior work.\n\nOverall: this is a solid, unflashy modeling paper with an honest presentation. The central interpretation is plausible but not uniquely determined by the data. It deserves a serious referee. My recommendation is to send it to review, with the expectation that the quoted probabilities will be presented as scenario-dependent estimates, not robust measurements. A constructive referee request would be a re-fit with the emission region outside the BLR, which the paper's own Figure 6 argues would be unproductive for neutrinos—but that still needs to be shown.","headline":"A transparent, standard leptohadronic model applied to two new z>1 FSRQ–neutrino correlations; the 2%/0.8% detection probabilities are scenario-dependent because the BLR-location assumption carries the whole neutrino efficiency claim.","tokens_in":31920,"tokens_out":1286,"would_cite":false,"duration_ms":15676,"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 paper argues that GB6 J2113+1121 and NVSS J171822+423948 are efficient PeV neutrino emitters, with cascade synchrotron radiation dominating their X-ray and GeV emission.","keywords":["neutrino astronomy","blazars","flat-spectrum radio quasars","leptohadronic model","p-gamma interactions","IceCube alerts","high-redshift blazars","multimessenger astronomy"],"falsifier":"Take a pointed X-ray observation of NVSS J171822+423948 during a gamma-ray flare: the model predicts a cascade-dominated X-ray component peaking near 0.5 keV with a flux tied to the hadronic luminosity, whereas a purely leptonic interpretation would give a different spectral shape and normalization. A spectrum showing primary-electron synchrotron or SSC dominance instead of cascade synchrotron, or radio-core or variability evidence placing the emission region beyond the BLR, would falsify the central claim.","tokens_in":30653,"feed_emoji":"🔭","tokens_out":11494,"duration_ms":96691,"temperature":0.7,"pith_summary":"The paper tries to establish that two very distant flat-spectrum radio quasars (FSRQs), GB6 J2113+1121 at redshift 1.3 and NVSS J171822+423948 at redshift 2.7, are genuinely efficient PeV neutrino emitters, although their distance makes individual IceCube detections unlikely (about 2% and 0.8% per gamma-ray flare). It builds a one-zone leptohadronic model in which protons in the jet collide with broad-line region (BLR) photons to produce pions and neutrinos, and fits the multiwavelength emission during the flares. A sympathetic reader would care because it would mean that high-redshift blazar jets in the early universe accelerate hadrons to extreme energies and that the X-ray band, not the gamma-ray band, is the most direct electromagnetic tracer of that neutrino production. The paper also argues that the hadronic cascade synchrotron emission dominates the X-ray and GeV emission of both sources, and predicts future IceCube detections of neutrinos from FSRQs with redshifts greater than 1.","feed_headline":"Two far-off blazars are efficient PeV neutrino emitters, model says","feed_subtitle":"Fits to IceCube-associated flares put neutrino production in the broad-line region and X-rays as the best tracer.","key_machinery":"The central object is the one-zone leptohadronic emission region: a spherical blob of radius $R \\sim 1$-$3\\times10^{16}$ cm with magnetic field $B \\sim 6$ G and $2.5$ G, moving with Doppler factor $\\delta \\sim 30$ and $20$, and filled with co-accelerated electron and proton distributions. The argument is carried by $p\\gamma$ pion production on BLR photons, whose energy density in the jet frame is boosted as $u_{\\rm BLR} = \\Gamma^2 L'_{\\rm BLR}/(4\\pi (r'_{\\rm BLR})^2 c [1+(r'_{\\rm in}/r'_{\\rm BLR})^3])$. Hadronic cascades are computed iteratively: $\\pi^0$-decay gamma rays and high-energy pair synchrotron photons annihilate on soft photons to produce successive generations of electron-positron pairs, and the synchrotron radiation of those pairs is what dominates the X-ray band. The $p\\gamma$ efficiency $f_{p\\gamma} \\approx t_{\\rm dyn}/t_{p\\gamma}$ converts proton luminosity into neutrino luminosity through $\\epsilon_\\nu L_{\\epsilon_\\nu} \\approx (3/8) f_{p\\gamma} \\epsilon_p L_{\\epsilon_p}$, and the BLR photon field is what raises this efficiency to the level the paper calls 'efficient' PeV neutrino emission.","core_discovery":"Under the assumption that the emission region sits at the BLR boundary ($r'_{\\rm in} \\approx r'_{\\rm BLR}$, about 0.1 pc for GB6 J2113+1121 and 0.2 pc for NVSS J171822+423948), the BLR photon field acts as the dominant target for $p\\gamma$ interactions. The model reproduces the observed spectral energy distribution with proton injection luminosities about two orders of magnitude above electron injection, neutrino spectra peaking near 2-3 PeV in the observer frame, and single-muon-neutrino detection probabilities of about 2% and 0.8% over the roughly 1-2 year gamma-ray flares. The central discovery is that electromagnetic cascades from pion decay and Bethe-Heitler pairs, radiated as synchrotron emission by the secondary pairs, dominate the X-ray band (around 0.5-1 keV) and contribute significantly above about 5 GeV, while BLR photons absorb gamma rays above roughly 5-10 GeV. The paper concludes that for these sources the X-ray flux is a better neutrino tracer than the gamma-ray flux.","pith_inferences":["One extension the paper leaves implicit is that quasi-simultaneous X-ray pointed observations of gamma-ray-alerted FSRQs would be an efficient way to select neutrino flares, because the cascade-dominated X-ray flux directly limits the expected neutrino flux.","A testable consequence of the BLR-absorption picture is that gamma-ray-selected samples may systematically miss neutrino emitters; a search for IceCube neutrinos from X-ray-selected, gamma-faint high-redshift FSRQs would probe this selection bias.","The same BLR-boundary assumption could be applied to other luminous FSRQs with measured disk luminosities to produce a ranked list of neutrino candidates, a ranking that stacking IceCube alerts on X-ray-selected flares could falsify.","If super-Eddington high-redshift FSRQ jets are common in the early universe, the cumulative contribution of such sources to the diffuse IceCube neutrino flux could be larger than the small per-source detection probabilities suggest; stacking analyses of flaring redshifts-greater-than-1 FSRQs would test this."],"forward_implications":["If the modeling is right, both sources are efficient PeV neutrino emitters, and their per-flare single-muon-neutrino detection probabilities are about 2% and 0.8%, low mainly because of their high redshifts.","X-rays, not gamma rays, are the most promising electromagnetic messenger for these two sources: hadronic cascade synchrotron dominates the X-ray band while BLR absorption suppresses gamma rays above about 5-10 GeV.","IceCube should accumulate further neutrino alerts from flat-spectrum radio quasars at redshifts greater than 1 during flaring states, with neutrino spectra peaking near 2-3 PeV.","A spectral break in the gamma-ray spectrum at a few GeV can be used as a marker that a flaring FSRQ is in a neutrino-productive (BLR-near) environment.","The required jet powers exceed the Eddington luminosities by factors of about 4-10 during the flares, implying that temporary super-Eddington jet power may be an ingredient for IceCube neutrino production."],"supporting_citations":[{"why":"Established the temporal and spatial correlation between IceCube-191001A and the GB6 J2113+1121 flare, and supplied its disk luminosity, black hole mass, and multiwavelength data.","marker":"Liao et al. (2022)"},{"why":"Established the IceCube-201221A association with NVSS J171822+423948 and supplied its disk luminosity, black hole mass, and light curves.","marker":"Jiang et al. (2024)"},{"why":"Provided the one-zone leptohadronic modeling approach and the TXS 0506+056 comparison values used for baryon loading and neutrino luminosity.","marker":"Keivani et al. (2018)"},{"why":"Supplied the iterative electromagnetic cascade scheme used to compute successive generations of pairs and their radiation.","marker":"Cerruti et al. (2015)"},{"why":"Provided the methods for computing the spectra of secondary particles from p-gamma and Bethe-Heitler processes.","marker":"Kelner & Aharonian (2008)"},{"why":"Gave the BLR characteristic radius and luminosity scaling used to set the BLR radius and energy density.","marker":"Ghisellini & Tavecchio (2008)"},{"why":"Provided the formula for the BLR radiation energy density in the jet comoving frame.","marker":"Hayashida et al. (2012)"},{"why":"Provided the p-gamma energy-loss timescale used to compute the neutrino production efficiency.","marker":"Waxman & Bahcall (1997)"},{"why":"Justified the choice of low proton acceleration efficiency and the link between maximum proton energy and the sub-PeV neutrino flux.","marker":"Oikonomou et al. (2019)"}],"fun_headline_variants":["Distant blazars seen as PeV neutrino factories","X-rays, not gamma rays, betray far blazar neutrinos","Two ancient quasars can churn out PeV neutrinos","At z>2, X-rays reveal blazar neutrinos"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the emission region sits at the boundary of the broad-line region, about 0.1-0.2 pc from the black hole, so BLR photons are the main target for proton-photon collisions; if the same flares instead originated outside the BLR, the paper's own comparison shows the neutrino production efficiency dropping by orders of magnitude and the claim that these are efficient PeV neutrino emitters collapsing.","fun_headline_variants_meta":{"raw":{"variants":["Distant blazars seen as PeV neutrino factories","X-rays, not gamma rays, betray far blazar neutrinos","Two ancient quasars can churn out PeV neutrinos","At z>2, X-rays reveal blazar neutrinos"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000715,"raw_usage":{"total_tokens":3305,"prompt_tokens":1130,"completion_tokens":2175,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":746,"completion_tokens_details":{"reasoning_tokens":2104}},"tokens_in":746,"tokens_out":2175,"duration_ms":16906,"temperature":1.0,"reasoning_tokens":2104,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:36:48.386644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a pointed X-ray observation of NVSS J171822+423948 during a gamma-ray flare: the model predicts a cascade-dominated X-ray component peaking near 0.5 keV with a flux tied to the hadronic luminosity, whereas a purely leptonic interpretation would give a different spectral shape and normalization. A spectrum showing primary-electron synchrotron or SSC dominance instead of cascade synchrotron, or radio-core or variability evidence placing the emission region beyond the BLR, would falsify the central claim.","supporting_citations":[{"cited_title":"2019, , 489, 4347, 10.1093/mnras/stz2246","cited_arxiv_id":null,"evidence_quote":"Justified the choice of low proton acceleration efficiency and the link between maximum proton energy and the sub-PeV neutrino flux."}],"review_version":1}