{"id":"7f6deabb-7410-4690-909f-4e5af7e0a507","arxiv_id":"2607.07202","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":10,"one_line_summary":"Simultaneous NICER+NuSTAR spectra of four TeV blazars show that 6 spectra require a blackbody component interpreted as accretion disk emission during low-flux states, plus Gaussian features at 1.4–1.7 keV in Mrk 421.","lead":"This paper analyzes X-ray spectra of four well-known TeV blazars using simultaneous NICER and NuSTAR observations, finding that 6 of 13 spectra require an additional blackbody component beyond the standard log-parabolic model. The authors interpret this as accretion disk emission becoming visible when the jet is in a low-flux state.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The bbodyrad-to-disk identification is unsupported: no significance test, and fitted temperatures vary by 10× across epochs of the same source without discussion.","rationale":"The reader correctly identified the load-bearing concern: the physical identification of bbodyrad with disk emission is the weak point, and the three sub-concerns (temperature inconsistency, physical implausibility, no significance test) are all valid. I agree with the reader's CONDITIONAL verdict. The paper's data reduction and analysis of the 7 cleanly-fit spectra are sound, and the flux-state context from MAXI/GSC light curves is a reasonable circumstantial argument. But the central novel claim — disk detection in BL Lac objects — requires more than improved χ²_red and a circumstantial flux-state argument. An F-test is the minimum bar; addressing the 10× temperature discrepancy is essential. Without these, the claim should not be accepted as established. The verdict should remain CONDITIONAL pending these checks. I note that the paper is honest about the Gaussian features possibly being instrumental, which is appropriate, but the same skepticism should have been applied to the bbodyrad component.","tokens_in":31708,"tokens_out":2166,"duration_ms":203330,"concrete_test":"For each of the 6 spectra where bbodyrad was added, perform an F-test comparing constant*TBabs*logpar against constant*TBabs*(bbodyrad+logpar). Report the F-statistic and p-value. If p > 0.05 for any spectrum, the added component is not statistically justified and the disk-contribution claim fails for that epoch. Additionally, refit the 4 Mrk 421 spectra with kT frozen to a physically motivated AGN disk temperature (~0.05 keV for a 10^9 M_sun BH); if the fit quality does not degrade, the bbodyrad component was not constraining a real disk temperature.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim — disk contribution in 6 spectra — rests entirely on interpreting an added bbodyrad component as accretion disk emission. This identification is fragile on two grounds. First, no statistical test (F-test, Bayesian evidence, or even Δχ² with degrees of freedom) is reported to establish that the bbodyrad component is a statistically significant improvement over the absorbed LP model alone. The χ²_red values change modestly (e.g., Mrk 501 MJD 59662: 0.945 → 0.906; MJD 59665: 1.10 → 1.01), and without a formal test it is unclear whether this reflects a real spectral component or overfitting soft-band residuals. Second, the fitted disk temperatures are internally inconsistent for the same source: Mrk 421 MJD 60077 gives kT = 1.51 keV while MJDs 60289/60296/60298 give kT ≈ 0.15 keV — a factor of 10 difference that is never discussed. For a ~10^9 M_sun black hole, the Shakura-Sunyaev inner disk temperature should be ~10–50 eV; both 1.51 keV and 0.15 keV are high, but 1.51 keV is in the regime of stellar-mass BH X-ray binaries, not AGN. The wild temperature variation suggests the bbodyrad component is absorbing different soft-band systematics across epochs rather than tracing a stable physical disk. The authors themselves acknowledge that the Gaussian features at 1.42–1.70 keV may be instrumental (Al-Kα at ~1.56 keV), which raises the broader question of whether NICER soft-band calibration residuals — not disk emission — drive the need for the bbodyrad component.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript presents X-ray spectral analysis of four classical TeV blazars (Mrk 421, Mrk 501, PG 1553+113, PKS 2155-304) using 13 simultaneous NICER and NuSTAR observations. Seven spectra are well fitted by an absorbed log-parabolic (LP) model. For four Mrk 421 and two Mrk 501 spectra, the authors add a bbodyrad component (and, for Mrk 421, a Gaussian line) to achieve acceptable fits. Using MAXI/GSC light curves, the authors argue these six epochs correspond to moderate-to-low or low flux states, and conclude that the bbodyrad component represents accretion disk emission visible when jet emission is diminished. The simultaneous NICER+NuSTAR data for PG 1553+113 and PKS 2155-304 are fitted with the LP model (plus an edge for PG 1553+113).","tokens_in":32043,"tokens_out":1116,"duration_ms":166910,"significance":"The paper provides a useful compilation of simultaneous soft/hard X-ray spectral fits for these well-studied sources, and the NICER+NuSTAR combination provides good broadband coverage. The hypothesis that disk emission may become detectable in BL Lacs during low-flux states is interesting and testable. However, the central claim of disk contribution rests on an identification that is not yet adequately supported.","major_comments":[{"comment":"§3, §5, Table 2: The identification of the bbodyrad component as accretion disk emission is the load-bearing claim of the paper, but it is not adequately justified. (1) No statistical test (F-test, Bayesian evidence, or explicit Δχ²/ΔDOF comparison) is reported to establish that the bbodyrad component is a statistically significant improvement over the absorbed LP model alone. The χ²_red changes are modest (e.g., Mrk 501 MJD 59662: 0.945→0.906; MJD 59665: 1.10→1.01), and it is unclear whether this reflects a real spectral component or fitting of soft-band residuals. (2) The fitted temperatures are internally inconsistent for the same source: Mrk 421 MJD 60077 gives kT = 1.51 keV while MJDs 60289/60296/60298 give kT ≈ 0.15 keV — a factor of 10 difference that is never discussed. For a ~10^9 M_sun black hole, the expected inner disk temperature is ~10–50 eV; both values are high, but 1.51-","section":null}],"minor_comments":[{"comment":"§4.1: The Gaussian features at 1.42–1.70 keV are acknowledged as possibly instrumental (Al-Kα at ~1.56 keV, Si-Kα at ~1.84 keV) or due to background contamination. Given that these features appear only in the four Mrk 421 spectra where bbodyrad is also added, the authors should discuss whether NICER soft-band calibration residuals — rather than disk emission — could be driving the need for additional components.","section":null},{"comment":"Table 2: The bbodyrad normalization for Mrk 421 MJDs 60289/60296/60298 is extremely small (of order 10^-4), with large asymmetric errors. The physical meaning of these values should be clarified.","section":null},{"comment":"§2.1.1: The choice of SCORPEON background model over 3C50 and Space Weather is justified, but it would help to state whether the SCORPEON model was fitted jointly with the source model in XSPEC or whether a fixed background spectrum was used.","section":null},{"comment":"Figure 3 and Figure 5: The MAXI/GSC flux values during the observational epochs (especially for Mrk 501 MJDs 59662 and 59665, where MAXI data were unavailable) are inferred from neighboring observations. This should be stated more explicitly in the text, not only in the figure captions.","section":null},{"comment":"§4.3: The edge feature at ~1.87 keV in PG 1553+113 is attributed to silicon. Is this an instrumental edge in NICER? Please clarify whether this is astrophysical or instrumental.","section":null},{"comment":"The manuscript would benefit from a brief discussion of whether the bbodyrad temperatures for Mrk 501 (kT ≈ 1.71–1.75 keV) are physically plausible for disk emission from a ~10^9 M_sun BH, or whether an alternative soft-band component (e.g., a steep power law or calibration residual) might be more appropriate.","section":null}],"recommendation":"major_revision","confidential_remarks":"The paper presents a reasonable observational effort but the interpretation overreaches. The bbodyrad-to-disk identification is the central claim and it is not supported by the data as presented: no significance test, implausible and inconsistent temperatures, and possible NICER calibration systematics in the soft band. The authors should either strengthen the claim with proper statistical tests and physical justification, or reframe the conclusion as a soft-band excess of uncertain origin. The Mrk 501 disk claim is presented as novel (no previous studies claiming disk contribution), which increases the burden of evidence."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"The headline: the central claim — accretion disk emission detected in 6 spectra of Mrk 421 and Mrk 501 during low-flux states — rests on adding a bbodyrad component to improve soft-band fits, but the physical identification is not well-supported and the fitted temperatures are internally inconsistent without discussion.","headline":"The bbodyrad-to-disk identification is unsupported: no significance test, and fitted temperatures vary by 10× across epochs of the same source without discussion.","tokens_in":32860,"tokens_out":139,"would_cite":false,"duration_ms":112337,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Quiet blazars reveal accretion disk signatures in X-ray spectra","keywords":[],"falsifier":"If the same six spectra can be equally well fitted by alternative soft-band models — such as a broken power law, a more complex synchrotron curvature, or a calibration artifact correction — without invoking a blackbody, the disk-contribution interpretation would lose its primary support.","tokens_in":31924,"feed_emoji":"🛰️","tokens_out":1171,"duration_ms":157083,"temperature":0.7,"pith_summary":"The paper examines 13 simultaneous X-ray spectra from two satellites (NICER and NuSTAR) covering four well-known blazars — Mrk 421, Mrk 501, PG 1553+113, and PKS 2155-304. Blazars are galaxies with supermassive black holes that launch jets of particles aimed nearly at Earth; their X-ray emission is normally dominated by synchrotron radiation from these jets and is well described by a curved spectral shape called the log-parabolic model. The authors find that for seven of the thirteen spectra this standard model works. But for four observations of Mrk 421 and two of Mrk 501, the standard jet-only model leaves systematic residuals in the soft X-ray band (below about 4 keV). Adding a blackbody component — a thermal spectrum characteristic of a hot accretion disk — improves the fits. Crucially, the authors show from independent long-term monitoring data that during all six of these epochs the blazars were in moderate to low flux states, meaning the jet was relatively dim. The argument is that when the jet fades, the thermal glow of the accretion disk, normally swamped by jet emission, becomes detectable. This is the first reported evidence of accretion disk contribution in the X-ray spectra of Mrk 501, and it extends earlier findings for Mrk 421. The paper also notes an unexplained Gaussian emission line near 1.4–1.7 keV in four Mrk 421 spectra, which the authors attribute to possible instrumental effects or background contamination rather than a physical origin in the source.","feed_headline":"Quiet blazars reveal accretion disk signatures in X-ray spectra","feed_subtitle":"When jet emission dims, thermal disk light emerges in six observations — the first such detection for Mrk 501.","key_machinery":"The central mechanism is the relative flux balance between two emission components: nonthermal synchrotron radiation from the relativistic jet (modeled as a log-parabola) and thermal radiation from the accretion disk (modeled as a blackbody, bbodyrad in XSPEC). The argument rests on the idea that the jet is Doppler-boosted and variable, while the disk is steadier. During high-flux states the jet dominates and the disk is invisible; during low-flux states the jet dims enough for the disk's soft X-ray thermal component to emerge as a residual in spectral fits. The authors use the 2–20 keV MAXI/GSC light curve as an independent flux-state diagnostic to establish that the six epochs requiring a黑","core_discovery":"When two classical TeV blazars (Mrk 421 and Mrk 501) are in moderate to low X-ray flux states, their spectra cannot be described by the standard jet-only log-parabolic model alone. Adding a blackbody component — interpreted as thermal emission from the accretion disk — produces satisfactory fits in six cases, suggesting that disk emission becomes visible when jet synchrotron radiation dims. This is the first such claim for Mrk 501.","pith_inferences":["The large discrepancy in fitted disk temperatures for Mrk 421 (1.51 keV for one epoch versus ~0.15 keV for three others) suggests the blackbody component may be absorbing different unmodeled spectral structure rather than tracing a single physical disk. If so, the disk-contribution interpretation is fragile.","A blackbody temperature of 1.51 keV for a black hole of ~10^9 solar masses is physically implausible under standard thin-disk models, where the inner disk temperature should be orders of magnitude lower. This raises the possibility that the added component is fitting residual curvature or calibration artifacts rather than disk emission.","Without a formal statistical test (such as an F-test or Bayesian model comparison), it is unclear whether the improvement from adding the blackbody component is statistically significant or reflects overfitting noise in the soft band."],"forward_implications":["If accretion disk emission is genuinely detectable in BL Lac objects during low states, systematic monitoring of blazar low-flux epochs could become a new method for estimating inner disk temperatures and, by extension, black hole masses and accretion rates in sources where the disk is otherwise invisible.","The finding could refine blazar unified models: BL Lacs are traditionally classified as having radiatively inefficient, geometrically thick accretion disks with negligible disk emission. Detectable thermal disk components would challenge that picture or suggest a transition between disk states.","The unexplained Gaussian line near 1.5 keV in Mrk 421, if confirmed as astrophysical rather than instrumental, could point to a previously unrecognized emission process in blazar X-ray spectra during low states."],"fun_headline_variants":["Low-flux TeV blazars reveal accretion disk in X-ray spectra","Dimmed jets expose disk contribution in TeV blazar X-rays","Quiet TeV blazars break jet-only model, show disk emission","Accretion disk signatures emerge in low-flux blazar X-rays","Faint blazar states reveal thermal disk light in X-ray band"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The blackbody component added to improve spectral fits is assumed to represent genuine thermal emission from the accretion disk. This identification rests on the coincidence that the epochs requiring the extra component occur during low flux states, but the fitted temperatures are internally inconsistent across epochs and, for one case, physically implausible for the inferred black hole mass. No statistical significance test is reported to distinguish a real disk signal from软","fun_headline_variants_meta":{"raw":{"variants":["Low-flux TeV blazars reveal accretion disk in X-ray spectra","Dimmed jets expose disk contribution in TeV blazar X-rays","Quiet TeV blazars break jet-only model, show disk emission","Accretion disk signatures emerge in low-flux blazar X-rays","Faint blazar states reveal thermal disk light in X-ray band","Low-flux blazars depart from standard jet-only X-ray model","Disk emission detected when TeV blazar jets dim","TeV blazars at low flux show disk emission in X-ray spectra"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":2089,"prompt_tokens":545,"completion_tokens":1544,"prompt_tokens_details":null},"tokens_in":545,"tokens_out":1544,"duration_ms":62754,"temperature":1.0,"reasoning_tokens":1227,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-09T17:17:54.915696+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the same six spectra can be equally well fitted by alternative soft-band models — such as a broken power law, a more complex synchrotron curvature, or a calibration artifact correction — without invoking a blackbody, the disk-contribution interpretation would lose its primary support.","supporting_citations":[],"review_version":1}