{"id":"31ee24a7-d3fd-4aa6-b756-5f09804e2424","arxiv_id":"2501.14137","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"For 23 high-redshift blazars, one-zone SED modeling favors the molecular torus over the broad-line region as the source of seed photons for gamma-ray emission.","lead":"This paper studies 30 distant blazars (beyond redshift 2.5) using 15 years of gamma-ray data and models their light output to ask where the gamma-ray emission is produced. It finds that infrared seed photons from the dusty torus explain the data better than photons from the broad-line region, suggesting the emission occurs far from the black hole.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MT-over-BLR chi-square ranking is not robust: the paper's own R_MT definition is internally inconsistent, changing the MT seed-photon energy density by a factor ~6.","rationale":"The reader's weakest assumption was the fixed environmental parameters (τ_BLR, τ_MT, T_MT, and the R_BLR/R_MT scaling). My stress-test sharpens this: the paper itself contains a factor-2.5 conflict in the definition of R_MT, which translates to a factor of roughly 6 in the MT seed-photon energy density. Because the MT-versus-BLR chi-square differences are often below 3 and no significance test is reported, this internal inconsistency sits directly under the central claim. The concrete re-fit test would settle whether the discrepancy is a harmless typo or a model-definition error that could flip the conclusion. I do not recommend rejecting the paper outright; a conditional acceptance with the re-fit as a requirement matches the reader's verdict. I only partially agree with the reader's framing because the specific R_MT inconsistency is a stronger, checkable flaw than the general worry about parameter choices, and the reader did not identify this particular conflict. If the re-fit preserves the MT preference, the central claim would survive this attack, though the broader lack of model-selection statistics and the poor absolute fit quality would remain weaknesses.","tokens_in":29234,"tokens_out":15224,"duration_ms":130360,"concrete_test":"Re-run the JetSet fits for all 23 sources twice: once with the R_MT values currently in Table 1, and once with R_MT = 10^18 sqrt(L_disk/10^45) cm as printed in the Table 1 note (i.e., divide the current R_MT values by 2.5), keeping all other settings, parameters, and fitting procedures fixed. Compare the sign, magnitude, and statistical significance of Δχ^2 = χ^2_BLR − χ^2_MT in the two runs. If the MT preference disappears or becomes non-significant in the second run, the central claim is not supported as stated; if it persists, the inconsistency is an inconsequential typo.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the MT scenario is systematically better than the BLR rests on the chi-square comparison in Table 3, which depends on the MT seed-photon energy density u'_MT ∝ τ_MT L_disk / (R_MT^2). The paper is internally inconsistent about R_MT: the Table 1 note states R_MT = 10^18 sqrt(L_disk/10^45) cm, but the Table 1 R_MT values (in units of 10^19 cm) are 2.5 times larger, i.e., R_MT = 2.5×10^18 sqrt(L_disk/10^45) cm. This factor changes u'_MT by about 6.25. Many of the reported gamma-ray chi-square differences are small (e.g., J1510.1: 2.95 vs 1.86; J1127.4: 1.21 vs 0.61; J2313.9: 2.3 vs 1.76), so a factor-of-6 change in the seed photon field can plausibly flip the ranking. No Δχ^2 significance or model-selection statistic is reported, so the abstract's 'systematically better' claim is asserted without a formal test. The inference that the emission region lies outside the BLR is therefore not robust until the fitting is repeated with a consistent R_MT (or equivalently with τ_MT scaled by 1/6.25).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies 30 Fermi-LAT-detected blazars at z>2.5 using 15 years of gamma-ray data. It fits the gamma-ray spectra with power-law, power-law-with-cutoff, and log-parabola models, finding significant curvature for most sources. For 23 sources with multiwavelength SEDs, the authors apply a one-zone leptonic jet model using JetSet with emcee optimization, placing the dissipation region either at the broad-line region (Rdiss=RBLR) or at the molecular torus (Rdiss=RMT). The central claim is that the MT scenario gives systematically lower gamma-ray-band chi-square values than the BLR scenario, implying that the gamma-ray emission region lies outside the BLR. The paper also reports jet power, magnetic field, energy densities, and correlations among derived parameters, including a claimed accretion-jet connection via the Pjet-Ldisk correlation.","tokens_in":29533,"tokens_out":12453,"duration_ms":102557,"significance":"If the MT-over-BLR result is robust, it is astrophysically valuable: it locates the gamma-ray dissipation zone at ~10^18 cm from the central engine for z>2.5 blazars, constraining external-Compton seed photon fields in the early Universe. The sample is sizable and the use of a public SED modeling code with MCMC optimization is appropriate. The paper provides complete parameter tables and per-source SED plots, which is helpful for reproducibility. However, the central claim currently rests on a raw chi-square comparison without a significance test, and the MT radius used in the fits is internally inconsistent with the stated scaling. These issues must be addressed before the main conclusion can be accepted.","major_comments":[{"comment":"The MT radius is defined inconsistently. The Table 1 note gives R_MT ≃ 10^18 sqrt(L_disk/10^45) cm, but the tabulated R_MT values (in units of 10^19 cm) are 2.5 times larger: for J1510.1+5702, log L_disk=46.80 gives 0.79×10^19 cm from the formula, while the table lists 1.99×10^19 cm. Since the MT seed-photon energy density scales as u'_MT ∝ τ_MT L_disk / R_MT^2, this changes u'_MT by a factor of about 6.25 between the stated scaling and the values actually used. The MT-vs-BLR chi-square ranking in Table 3 depends directly on u'_MT, and several reported differences are small (e.g., for J1510.1+5702 the gamma-band chi-square is 2.95 for BLR and 1.86 for MT; for J1127.4+5648 it is 1.21 vs 0.61; for J2313.9-4501 it is 2.3 vs 1.76). The fits should be redone with a single, consistent R_MT, and the sensitivity of the ranking to this choice should be reported.","section":"§3.1 and Table 1"},{"comment":"The abstract's claim that MT modeling is 'systematically better' is not supported by a statistical test. Table 3 reports only the gamma-band chi-square and the all-band reduced chi-square, with no Delta chi-square, number of gamma-band degrees of freedom, p-value, or information criterion. The raw differences are small for several sources, and the reduced chi-square values are large in both scenarios (e.g., J1510.1+5702: 12.73 vs 14.94; J0746.4+2546: 29.31 vs 27.80; J1344.2-1723: 26.38 vs 16.07), indicating poor absolute fits. In J1510.1+5702 the gamma-band chi-square favors MT while the all-band reduced chi-square favors BLR, so the two metrics conflict. The authors should report Delta chi-square with the appropriate degrees of freedom, a model-selection statistic, and a discussion of what the poor absolute reduced chi-square values imply for the relative comparison.","section":"§3.2 and Table 3"},{"comment":"Several of the reported spectral fits are unphysical. Best-fit photon indices are negative in a number of cases, e.g., alpha_PL = -1.09 +/- 0.18 for J0539.6+1432, alpha_PL = -1.05 +/- 3.06 for J2320.8-0823, alpha_PL = -1.79 +/- 0.36 for J0224.9+1843, and alpha_PLC = -1.19 +/- 0.33 for J2015.4+6556. A negative photon index implies a spectrum rising with energy, which is not expected in the Fermi-LAT band and likely indicates fit convergence problems or very low signal-to-noise. These fits are nevertheless included in the sample statistics and should be repaired, rejected, or explicitly discussed as nonphysical or upper-limit cases.","section":"Table 2"},{"comment":"The correlations used to support the accretion-jet connection are at risk of being partly induced by construction. L_disk enters the model as an input to the external photon fields that determine the IC emission, and the derived jet power includes the radiation power Pr; therefore a positive Pjet-Ldisk correlation (and the nu_IC-Ldisk anticorrelation) may reflect the input-output structure of the fitting procedure rather than an independent physical relation. A partial-correlation analysis, or a test in which L_disk is removed from the seed-photon inputs while other parameters are refit, is needed before these correlations can be used as evidence.","section":"§4.2 and Figure 4"}],"minor_comments":[{"comment":"Some source names in the appendix figures are inconsistent with Table 1, e.g., 'J10224.9+1843' should be 'J0224.9+1843' and 'J0914.2+4127' should be 'J0912.2+4127'.","section":"Figure 12 and Appendix B"},{"comment":"The paper should state explicitly that the correlation analysis uses the Rdiss=RMT parameter set, since both scenarios are fitted in Table 3 and the choice affects the reported correlations.","section":"§3.3 and §4.1"},{"comment":"The fixed environmental parameters (tau_BLR=0.1, tau_MT=0.2, T_MT=10^3 K, and R proportional to sqrt(L_disk)) are load-bearing for the MT/BLR comparison; a robustness test varying these values, or at least a quantitative discussion of the resulting systematic uncertainty, should be added.","section":"§3.1"},{"comment":"There are numerous typographical artifacts in the text (e.g., 'B la˙ zejowski' in the introduction and apparent table-splitting artifacts in the extraction); the manuscript should be carefully proofread before resubmission.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The MT-over-BLR claim is interesting and potentially important, but the internal inconsistency in R_MT and the absence of any significance test mean the central result is not yet established. If the authors rerun the fits with a consistent R_MT and report a proper model comparison, the paper could become acceptable. The correlation analysis also needs to be framed more cautiously because of the input-output coupling with L_disk."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful core is the sample study: 30 high-redshift blazars, 15 years of Fermi-LAT data, careful PL/PLC/LP spectral fits, and one-zone leptonic SED modeling for 23 of them with a genuine BLR-versus-MT dissipation-region comparison. The sample and spectral fits are a reasonable resource, and the qualitative conclusion that the torus provides the seed photons matches earlier work (Sahakyan et al. 2020; Wu et al. 2024). The paper is not circular; the MT/BLR comparison is a real model comparison.\n\nThe problems are in the evidence for the central claim. No significance statistic is reported; the paper states MT is 'systematically better' based on chi-square, but the reduced values are often >10, so we are ranking two poor fits. The stress-test inconsistency is real: Table 1's R_MT values are 2.5 times the stated formula, changing the MT seed-photon energy density by about a factor of 6. Several MT-vs-BLR chi-square differences are small, so that factor can easily flip the ranking. The paper does not say which R_MT was used in the fits. Third, some fitted photon indices are negative and some errors are enormous; the authors flag them, but they still feed into the correlations. And the Pjet-Ldisk correlation is partly constructed because Ldisk sets the external photon field in the same model that gives Pjet.\n\nWhat holds up: the sample selection is transparent, the analysis is standard, and the direction of the result is plausible. It just isn't demonstrated yet.\n\nI'd send it to a referee, not desk reject it, but I'd expect heavy revision: fix the R_MT inconsistency, report a formal model-selection statistic (delta-chi^2 or AIC/BIC), vary the environmental parameters (tau_BLR, tau_MT, T_MT), and show the ranking survives. The paper is for people working on Fermi blazars and jet dissipation; it would be a decent resource once those checks are in.","headline":"Useful sample of z>2.5 blazars with a plausible but unproven MT-over-BLR conclusion; the chi-square ranking is not robust as presented.","tokens_in":30103,"tokens_out":4304,"would_cite":false,"duration_ms":36402,"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 gamma-ray-emitting regions of z>2.5 blazars are most likely located outside the broad-line region, near the molecular torus, based on systematically lower chi-square values in one-zone leptonic SED fits.","keywords":["high-redshift blazars","gamma-ray spectra","Fermi-LAT","spectral energy distribution","one-zone leptonic model","external Compton scattering","molecular torus","broad-line region"],"falsifier":"Re-fit the same 23 sources while varying the BLR reprocessing fraction (tau_BLR) from 0.01 to 0.3 and the torus temperature (T_MT) from 500 to 2000 K; if a substantial number of sources then show lower chi-square for the BLR scenario, the claim that the emitting region lies outside the BLR is not robust.","tokens_in":28993,"feed_emoji":"🔭","tokens_out":4993,"duration_ms":40430,"temperature":0.7,"pith_summary":"This paper studies 30 blazars with redshifts above 2.5 using 15 years of gamma-ray observations. It finds that most of their gamma-ray spectra are curved, and it fits the broadband spectral energy distributions of 23 of them with a one-zone leptonic emission model. The central claim is that the gamma-ray-emitting region is more likely located near the molecular torus than inside the broad-line region, because the torus scenario gives systematically lower chi-square values. If true, this locates the dissipation zone at about $10^{18}$ cm from the supermassive black hole and means infrared seed photons, not UV line photons, dominate the external Compton emission. The paper also reports that these high-redshift blazars have higher jet powers and disk luminosities than lower-redshift counterparts, with a positive jet power–disk luminosity correlation supporting an accretion–jet connection.","feed_headline":"Distant blazars shine from outside their broad-line regions","feed_subtitle":"A 23-source SED comparison finds infrared, not UV, seed photons dominate the jet's gamma-ray output.","key_machinery":"The central object is the one-zone leptonic emission model with a log-parabolic electron energy distribution, in which the high-energy peak is produced by external Compton scattering of seed photons from either the broad-line region (UV) or the molecular torus (IR). The argument is carried by comparing the chi-square of fits with the dissipation region placed at the BLR radius versus the torus radius, with both radii set by the scaling R proportional to the square root of the accretion disk luminosity. The lower chi-square of the torus scenario is what locates the emission outside the BLR.","core_discovery":"For 23 high-redshift blazars (z>2.5) with multiwavelength data, fitting a one-zone leptonic model with either broad-line-region or molecular-torus seed photons yields systematically lower chi-square values for the torus scenario in the gamma-ray band, so the gamma-ray-emitting region is most likely located outside the broad-line region, at distances of order $10^{18}$ cm from the central engine. The gamma-ray spectra of most of the 30 sources show significant curvature, modeled better by a power law with an exponential cutoff or a log-parabola than by a simple power law. High-redshift blazars exhibit higher gamma-ray luminosities, softer spectral indices, higher jet powers, and higher accretion disk luminosities than their low-redshift counterparts, and the data support a positive correlation between jet power and disk luminosity.","pith_inferences":["If the torus location is correct, the gamma-ray emitting region lies beyond the zone probed by UV broad emission lines, so gamma-ray variability timescales may be longer than the BLR light-crossing time.","The chi-square ranking depends on fixed environmental parameters, so a systematic variation of the BLR and torus reprocessing fractions could change individual source classifications even if the sample-average conclusion holds.","The same one-zone comparison could be applied to lower-redshift flat-spectrum radio quasars to test whether an external location of the gamma-ray zone is universal or specific to the high-luminosity, high-redshift regime.","Because the EBL absorption correction significantly affects gamma-ray fluxes above 10^25 Hz, improved EBL models would directly alter the inferred intrinsic spectra and the chi-square comparison between the two seed-photon scenarios."],"forward_implications":["The gamma-ray dissipation zone in high-redshift blazars sits at distances of order 10^18 cm from the central engine, outside the broad-line region.","External Compton emission in these sources is dominated by infrared seed photons from the molecular torus, so the high-energy peak of the SED is set by the torus temperature and energy density.","High-redshift blazars have jet powers that are systematically larger than their accretion disk luminosities, supporting an accretion–jet connection in the early universe.","Blazars with higher disk luminosities tend to have lower IC peak frequencies, implying stronger external photon fields and more efficient cooling of the radiating electrons.","Most high-redshift blazar gamma-ray spectra show significant curvature, so simple power-law fits are insufficient and the electron energy distribution is better described by a log-parabola."],"supporting_citations":[{"why":"Supplies the scaling relations R_BLR and R_MT proportional to the square root of disk luminosity, which set the dissipation distances used in both scenarios.","marker":"Ghisellini & Tavecchio 2008"},{"why":"Earlier study of high-redshift blazars whose conclusion of radiation region near the molecular torus is cited as consistent with the present finding.","marker":"Sahakyan et al. 2020"},{"why":"Recent study of high-redshift blazars that also supports the torus location of the gamma-ray emitting region.","marker":"Wu et al. 2024"},{"why":"Provides the extragalactic background light model used to convert intrinsic gamma-ray fluxes to observed fluxes in the SED fits.","marker":"Finke et al. 2010"},{"why":"Provides the one-zone leptonic radiative model used to compute the theoretical SEDs in the fitting procedure.","marker":"Tramacere et al. 2011"}],"fun_headline_variants":["Blazars' gamma rays originate outside their broad-line regions","Infrared seed photons dominate distant blazar gamma-ray sites","High-redshift blazars' gamma-ray zone sits beyond the broad-line region","Blazar gamma-ray emission located outside the broad-line region","Infrared photons win over UV for high-redshift blazar gamma rays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison assumes fixed environmental parameters: a BLR reprocessing fraction of 0.1, a torus reprocessing fraction of 0.2, a torus temperature of 1000 K, and radii scaling with the square root of disk luminosity; if these values are wrong, the seed photon energy densities change and the chi-square ranking could flip.","fun_headline_variants_meta":{"raw":{"variants":["Blazars' gamma rays originate outside their broad-line regions","Infrared seed photons dominate distant blazar gamma-ray sites","High-redshift blazars' gamma-ray zone sits beyond the broad-line region","Blazar gamma-ray emission located outside the broad-line region","Infrared photons win over UV for high-redshift blazar gamma rays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000548,"raw_usage":{"total_tokens":2661,"prompt_tokens":1031,"completion_tokens":1630,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":1540}},"tokens_in":647,"tokens_out":1630,"duration_ms":11220,"temperature":1.0,"reasoning_tokens":1540,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:21:07.086043+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the same 23 sources while varying the BLR reprocessing fraction (tau_BLR) from 0.01 to 0.3 and the torus temperature (T_MT) from 500 to 2000 K; if a substantial number of sources then show lower chi-square for the BLR scenario, the claim that the emitting region lies outside the BLR is not robust.","supporting_citations":[],"review_version":1}