REVIEW 3 major objections 5 minor 4 references
A 1.7-day optical quasi-periodic oscillation in the quasar PKS 0805-07 passes high-confidence tests against red noise, pointing to a compact jet or disk structure.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 18:57 UTC pith:62DZD2NY
load-bearing objection A plausible but overclaimed TESS QPO candidate in PKS 0805-07; the 9-day window selection likely inflates the significance, but the paper deserves a careful referee. the 3 major comments →
A Possible Short-Timescale Optical Quasi-Periodic Oscillation in PKS\,0805-07 from High-Cadence TESS Observations
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The core claim is that a quasi-periodic modulation at P ≈ 1.7 d (f = 0.5968 ± 0.0351 d^-1) is present in the TESS Sector 34 optical light curve of PKS 0805-07, with the Lomb-Scargle peak exceeding the 99.99% confidence level and the independent weighted wavelet Z-transform recovering a consistent timescale (f = 0.5977 ± 0.0437 d^-1) at about the 99.9% confidence level. The wavelet map shows the signal is localized in time rather than persistent throughout the 25-day sector. The authors interpret the modulation as a transient quasi-periodic feature spanning approximately five cycles, superimposed on stochastic jet variability. They argue that a disk-based hotspot orbiting near the ISCO would
What carries the argument
The analysis rests on two complementary time-series tools: the Lomb-Scargle periodogram (a standard method for detecting periodic signals in unevenly sampled data) and the weighted wavelet Z-transform (WWZ, a time-frequency decomposition that reveals whether a periodic signal is persistent or localized in time). Significance is evaluated by generating 20,000 (for LSP) and 10,000 (for WWZ) Monte Carlo light curves that reproduce both the observed power spectral density and the flux distribution, following the Emmanoulopoulos approach, to set red-noise confidence levels. The period uncertainty is derived by fitting a Gaussian to the peak. A physical interpretation uses the ISCO orbital period–
Load-bearing premise
The 9-day window over which the periodicity is measured was chosen without regard to where a signal might appear; if the window was selected after seeing the peak, the Monte Carlo significance levels are systematically inflated by the number of trial windows and sectors examined.
What would settle it
Compute the Lomb-Scargle periodogram over the full 25-day TESS Sector 34 baseline (not just the 9-day window) and also over the other TESS sectors (7, 61, 88, 99) using the same pipeline, and count how many independent red-noise light curves produce a peak exceeding 99.99% significance in any 9-day window. If the 1.7-day peak disappears in the full-baseline analysis or appears in fewer than ~1 in 10^4 red-noise simulations after a look-elsewhere correction, the candidate oscillation would be a false positive.
If this is right
- If the 1.7-day oscillation is real, it would demonstrate that day-scale optical quasi-periodic modulation can occur in a jet-dominated FSRQ, extending the range of known blazar QPO timescales.
- The transient, five-cycle nature suggests the mechanism is episodic rather than a stable disk mode, favoring jet kink instabilities or short-lived disk hotspots over long-lived binary or precession scenarios.
- The inferred black hole mass (if the disk-hotspot interpretation holds) of ~10^8–10^9 M_sun is consistent with typical FSRQs and could provide an independent mass estimate for PKS 0805-07.
- The detection would motivate high-cadence multiwavelength follow-up (X-ray, radio, polarimetric) to test for correlated variability and to distinguish between disk-based and jet-based origins.
- If confirmed to recur in other TESS sectors, the signal could be linked to the longer-timescale gamma-ray periodicities previously reported for this source, suggesting a common geometric or physical origin.
Where Pith is reading between the lines
- The quoted confidence levels do not account for the fact that the 9-day analysis window was selected from a 25-day sector after inspecting the light curve; a look-elsewhere correction over possible window positions would reduce the effective significance, possibly below the 99.9% level.
- Because the signal spans only five cycles, the paper's own reference to Vaughan et al. (2016) implies that such a pattern can occasionally be produced by red noise; a formal false-alarm probability computed over the entire sector (rather than a pre-chosen window) would be a stronger test.
- A testable prediction of the kink-instability scenario is that the optical modulation should be accompanied by correlated variability in the polarized flux or in the X-ray/radio bands, with a timescale set by the jet's transverse size and Doppler factor.
- The same analysis applied to the other available TESS sectors (7, 61, 88, 99) would provide a direct test of recurrence; if the 1.7-day signal appears only in Sector 34, the transient interpretation is supported but the statistical evidence is weaker than claimed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the TESS Sector 34 optical light curve of the flat-spectrum radio quasar PKS 0805−07 over a 9-day sub-interval (BTJD 2231.4047–2240.4047; MJD 59230.90–59239.90). Using Lomb–Scargle periodogram (LSP) and weighted wavelet Z-transform (WWZ), it reports a dominant quasi-periodic modulation at f ≈ 0.597 d⁻¹ (P ≈ 1.7 d), with LSP significance exceeding 99.99% and WWZ significance ~99.9% based on Monte Carlo red-noise simulations. The modulation is claimed to be transient, spanning ~5 cycles, and is interpreted as either a disk hotspot near the ISCO (yielding MBH ~ 1.1×10⁸–7.2×10⁸ M⊙) or magnetohydrodynamic kink instabilities in the jet. The abstract and §5 conclude that the signal is consistent with a compact, short-lived structure embedded in stochastic jet variability.
Significance. If the detection and significance claims are correct, this would add a useful example of a short-timescale optical QPO in a blazar observed with TESS, complementing previous TESS-based QPO searches. The paper employs standard complementary period-search tools (LSP, WWZ) and Monte Carlo simulations following Emmanoulopoulos et al. (2013), which are appropriate for red-noise-dominated blazar light curves. The use of the public QUAVER pipeline and the reporting of Baluev FAP values are also positive features. However, the central significance claim rests on treating the 9-day analysis window as fixed in advance, while §2.2 indicates it was selected for the QPO investigation, and other TESS sectors exist but were not analyzed. This look-elsewhere effect is not accounted for in the Monte Carlo simulations, so the reported 99.99%/99.9% confidence levels are likely optimistic. The physical interpretations are illustrative but do not independently strengthen the detection.
major comments (3)
- [§2.2, §4] The analysis is restricted to the 9-day interval BTJD 2231.4047–2240.4047, which §2.2 states was 'selected from the 1.5 hr binned light curve for the QPO investigation.' The Monte Carlo significance estimates in §4 generate red-noise light curves over this same fixed window, so they do not account for the number of trials involved in selecting this window from the full 25-day sector, nor for the fact that the source is also observed in Sectors 7, 61, 88, and 99 (§2.1). If the window or sector was chosen after inspecting the data, the effective false-alarm probability is multiplied by the number of independent sub-windows/sectors searched. A modest trial factor of ~10–20 would raise the nominal Baluev FAP of 1.07×10⁻⁴ (which corresponds to 99.989%, not 99.99%) above 10⁻³, undermining the central claim that the signal 'exceeds the 99.99% confidence level.' The authors should either pre-spe
- [§3.1, Figure 2 caption] The text in §3.1 and §5 states that the LSP peak 'exceeds the 99.99% confidence level,' but the caption of Figure 2 says it 'exceeds the 99.9% confidence level derived from 2×10^4 Monte Carlo simulations.' These are inconsistent. Additionally, the Baluev FAP of 1.07×10⁻⁴ corresponds to 99.989%, not 99.99%. This numerical inconsistency must be resolved, and the significance level used in the abstract and conclusions must match the actual simulation output and the FAP calculation.
- [§3.2, §4] The claim that the WWZ 'shows that the signal is temporally localized rather than persistent across the full light curve' (§3.2 and abstract) is not supported by the analysis: the WWZ is computed only over the pre-selected 9-day interval, so localization within that interval cannot establish non-persistence over the full 25-day sector. To support the transient interpretation, the authors should present a WWZ map (or equivalent time-frequency analysis) over the entire Sector 34 light curve, or at least analyze other sub-windows of the same sector. As written, the transient claim is an artifact of the analysis window.
minor comments (5)
- [§2.2, Figure 1] The caption of Figure 1 says the red points are the hybrid method and the blue points are the Simple PCA (SPO) light curve, but the text says the SPO reduction was adopted for the timing analysis. Clarify which curve is used and ensure the label colors match the description.
- [§4, Figure 3] Figure 3 caption: '104 Monte Carlo simulations' should read '10^4 Monte Carlo simulations.' Also, the caption mentions dashed cyan and red curves for 99.9% and 99.99% confidence levels, but the text only claims ~99.9% for the WWZ; clarify the levels shown and which is relevant.
- [§3.2, Eq. (1)-(2)] The WWZ equations are written in integral form, but the implementation is necessarily discrete. A brief statement of the discrete weighted projection used would improve reproducibility.
- [§5] The 'disk-based hotspot' mass estimate spans MBH ≈ 1.1×10⁸ M⊙ (Schwarzschild) to 7.2×10⁸ M⊙ (maximal Kerr), a factor of ~7. The abstract quotes only the upper value (7.2×10⁸ M⊙) as the implication, which overstates the constraint. Either quote the full range or justify why the Kerr value is preferred.
- [§5] The kink-instability scenario uses R_KI = 10^16–10^17 cm to obtain T_obs ≈ 1.6–16 days, then selects R_KI ≈ 10^16 cm to match P ≈ 1.7 d. This is a fit to the observation, not a prediction. Consider framing this as an illustrative consistency check rather than as evidence for the model.
Circularity Check
No circular derivation; the QPO detection is validated against Monte Carlo red noise, with only a statistical window-selection caveat.
full rationale
The central detection chain is not circular. The LSP peak at f=0.596841 d^-1 and the WWZ peak at f=0.597719 d^-1 are both compared with Monte Carlo red-noise simulations that reproduce the observed PSD and PDF with the same sampling, so the claimed significance is not obtained by fitting the detected frequency to the noise model. The WWZ and LSP agreement is same-data cross-method consistency, not an independent confirmation from new data, but that is not circularity. Self-citations (Akbar et al. 2024, 2025; Akbar 2026; Nazir et al. 2026) are motivational or methodological, not load-bearing for the detection. The physical interpretations are also non-circular consistency arguments: the ISCO hotspot scenario derives M_BH from the observed period, and the kink scenario adopts a range of R_KI and notes that P lies at the lower end; neither is presented as an ab initio prediction of P. The main caveat is the selection of the 9-day interval (BTJD 2231.4047-2240.4047) 'selected from the 1.5 hr binned light curve for the QPO investigation' while other TESS sectors were not analyzed; if the window or sector was chosen after seeing the signal, the quoted 99.99%/99.9% confidence levels would need a trial-factor correction. That is a statistical validity concern, not a definitional or self-citation circularity. The paper itself acknowledges that the significance is marginal and that the limited number of cycles is insufficient for firm confirmation.
Axiom & Free-Parameter Ledger
free parameters (3)
- Red-noise PSD slope β and normalization A =
not reported
- Analysis-window boundaries BTJD 2231.4047-2240.4047 =
9 days
- Kink emitting-region size R_KI =
10^16-10^17 cm (effectively ~10^16 cm)
axioms (4)
- domain assumption Underlying variability is stationary red noise described by a single power-law PSD and a fixed PDF over the 9-day window
- domain assumption The optical light curve is jet-dominated synchrotron emission, so red noise is the appropriate null model
- standard math The observed period corresponds to the orbital timescale at the ISCO, giving M_BH = 3.23e4 P (r^1.5 + a)(1+z)
- domain assumption Representative blazar kink-instability parameters: Doppler factor δ=15 and transverse propagation speed ⟨v_tr⟩≈0.16c
read the original abstract
We present a timing analysis of the high-cadence optical light curve of the high-redshift flat-spectrum radio quasar PKS\,0805$-$07 obtained during \textit{TESS} Sector~34 (MJD~=~59230.90--59239.90). We search for short-timescale quasi-periodic oscillations (QPOs) using complementary time-series techniques, including the Lomb--Scargle periodogram (LSP) and the weighted wavelet $Z$-transform (WWZ), and evaluate their significance against red-noise variability using Monte Carlo simulations. The LSP reveals a dominant modulation at $f \approx 0.6\,\mathrm{d^{-1}}$ ($P \approx 1.7$\,d) exceeding the $99.99\%$ confidence level, while the WWZ independently recovers a consistent timescale at the $\sim$99.9\% level and shows that the signal is temporally localized rather than persistent. The modulation spans $\sim$5 coherent cycles, indicating a transient feature. However, the significance only marginally exceeds the highest WWZ confidence threshold, and the limited number of cycles is insufficient for a firm confirmation. In a disk-based scenario, orbital motion of a hotspot near the innermost stable circular orbit implies a black hole mass of $M_{\rm BH} \sim 7.2 \times 10^{8}\,M_\odot$, consistent with typical FSRQ values. Alternatively, magnetohydrodynamic kink instabilities in the relativistic jet can naturally produce day-scale variability for standard blazar parameters and account for its transient character. Compact SMBH binary and tidal disruption event scenarios are disfavoured by the jet-dominated nature of the source, the optical rather than X-ray character of the modulation, and a black hole mass well above that of confirmed QPE hosts. The candidate modulation is consistent with a compact, short-lived structure embedded within stochastic jet variability, and high-cadence multiwavelength monitoring will be essential to confirm its recurrence and constrain its physical origin.
Figures
Reference graph
Works this paper leans on
-
[2004]
The Importance of Discovering a 3:2 Twin-Peak Quasi-periodic Oscillation in an Ultraluminous X-Ray Source, or How to Solve the Puzzle of Intermediate-Mass Black Holes. Astrophys. J. Lett. 609, L63–L65. doi:10.1086/422810,arXiv:astro-ph/0402012. Akbar, S., 2026. A multi-technique search for year-scaleg-ray quasi-periodic modulation in the high-redshift FSR...
Pith/arXiv arXiv 2026
-
[2017]
Gamma-ray and optical oscillations of 0716+714, MRK 421, and BL Lacertae. Astron. Astrophys. 600, A132. doi:10.1051/0004-6361/ 201630288,arXiv:1701.04454. Scargle, J.D., 1982. Studies in astronomical time series analysis. ii- statistical aspects of spectral analysis of unevenly spaced data. As- trophysical Journal, Part 1, vol. 263, Dec. 15, 1982, p. 835-...
Pith/arXiv arXiv 1982
-
[2021]
Detection of quasi-periodic oscillations from the blazar PKS 1510-089. Journal of Astrophysics and Astronomy 42, 92. doi:10. 1007/s12036-021-09773-9. Liska, M., Hesp, C., Tchekhovskoy, A., Ingram, A., van der Klis, M., Markoff, S., 2018. Formation of precessing jets by tilted black hole discs in 3D general relativistic MHD simulations. Mon. Not. R. Astron...
Pith/arXiv arXiv 2018
-
[3152]
ˇCemelji´ c, M., Yang, H., Yuan, F., Shang, H., 2022
doi:10.1093/mnras/stw1412,arXiv:1606.02620. ˇCemelji´ c, M., Yang, H., Yuan, F., Shang, H., 2022. Formation of Episodic Jets and Associated Flares from Black Hole Accretion Sys- tems. Astrophys. J. 933, 55. doi:10.3847/1538-4357/ac70cc, arXiv:2205.10531. Wagoner, R.V ., Silbergleit, A.S., Ortega-Rodríguez, M., 2001. “Stable” Quasi-periodic Oscillations an...
Pith/arXiv arXiv 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.