REVIEW 4 major objections 5 minor 2 cited by
X-ray reverberation modelling of the continuum, optical/UV time-lags in quasars
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The UV/optical time lags of 128 quasars are fully consistent with X-ray reverberation, provided the X-ray corona sits above about 40 gravitational radii.
desk verdict First real application of the K21 model to quasar continuum lags, but the 'h > 40 Rg' lower limit is not supported by the paper's own 1σ error bars. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the K21b analytic time-lag model, built on the K21a disc response functions. These functions compute, at each disc radius, the extra time-variable UV/optical flux produced when X-rays from a point-like corona are absorbed by the disc, including relativistic light travel and bending, the disc ionisation state, and a radial-dependent ratio of external (X-ray) to internal (accretion) heating. Given a black hole mass, accretion rate, 2-10 keV luminosity, spin, colour-correction factor (the standard adjustment for electron scattering in the disc atmosphere), and corona height $h$, the model predicts the lag between two wavelengths; the authors fix all parameters from the literature and leave $h$ as the only free parameter, scanning 5 to 80 $R_{\rm g}$ (gravitational radii). A later colour-correction prescription is used to adapt the model to $f_{\rm col}$ values of 1 and 1.7 as well as the original 2.4.
What would settle it
A direct measurement of the X-ray corona height in one of these quasars — for example from an iron-K X-ray reverberation lag that locates the corona, or from microlensing size measurements — that placed the corona below about 20 gravitational radii would contradict the best-fit models, which require about 40-50 Rg. Alternatively, if future time-lag data for quasars above $10^9\,M_\odot$ (where the analytic relations are extrapolated) require heights below 40 Rg to fit, the paper's central claim would be falsified.
Extended reading notes
Core claim
On the paper's own terms, the discovery is that the continuum UV/optical time lags of luminous quasars can be explained by thermal reverberation of a standard Novikov-Thorne accretion disc, with no need for anomalously large disc radii. The mean rest-frame lag spectrum of 128 quasars is well fitted by the K21b analytic time-lag relations for both zero and maximal black hole spin, as long as the X-ray corona height is above about 40 gravitational radii. For a non-spinning black hole with colour correction 1.7, the best-fit height is 42 Rg ($\chi^2=11.4$ for 5 degrees of freedom); for a maximally spinning black hole with colour correction 2.4, it is 50 Rg ($\chi^2=11.3$). The residuals of individual quasars are consistent with a zero-mean Gaussian with the expected scatter, and the paper notes the same corona-height range independently explains the half-light radii of microlensed quasars as supporting evidence. The paper states this is the first time quasar continuum UV/optical time lags have been shown to be consistent with X-ray reverberation of standard discs.
Load-bearing premise
The load-bearing premise is that the analytic lag relations, originally calibrated for black holes at or below about $10^8\,M_\odot$ with accretion rate below 0.5, stay accurate for the quasars in this sample, which reach $10^9\,M_\odot$ and Eddington ratios up to 1; if they lose accuracy there, the inferred corona heights and the good fits would both change.
Editorial extensions
If this is right
- The observed UV/optical quasar lags do not require disc radii larger than standard thin-disc predictions; X-ray heating of a Novikov-Thorne disc accounts for them.
- The X-ray corona in quasars is typically located at or above about 40 gravitational radii, matching the height independently inferred from microlensed quasar half-light radii.
- Continuum time-lag measurements can constrain corona height, and the best-fitting height is insensitive to whether the black hole is non-spinning or maximally spinning.
- Earlier models that reported factor-of-2-3 lag discrepancies likely failed because they assumed a constant ratio of X-ray to internal heating across the disc, rather than computing it radius by radius.
- If the flat-disc assumption is relaxed, smaller corona heights may fit equally well, so the quoted heights are upper-side estimates.
Reading between the lines
- A natural extension the paper does not pursue is to split the sample by black hole mass or Eddington ratio and re-fit the corona height in each bin; a trend would reveal that height is not universal.
- If the best-fit heights are taken at face value, the X-ray corona sits well above the disc plane in luminous quasars, which bears on models of the disc-corona interface and on jet-launching geometry; the paper does not discuss this connection.
- Future wide-field time-domain surveys with longer, denser UV/optical light curves could apply the same stacking method in several redshift bins, and simultaneous far-UV monitoring could help break the spin-colour degeneracy the paper identifies.
- Because the flat-disc assumption tends to push the inferred height upward, a finite disc scale height modelled with radiative transfer could bring the best-fit heights down to around 20 Rg; this is a testable extension of the same model rather than a claim of the paper.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper compiles archival UV/optical continuum time-lag measurements for 128 quasars, computes a rest-frame mean time-lag spectrum in six wavelength bins, and fits it with the analytical X-ray reverberation models of Kammoun et al. (2021a,b), extended to different colour-correction factors by Kammoun et al. (2023). For each of six model variants (spin a*=0 or 1 with fcol = 1, 1.7, 2.4), the only free parameter is the X-ray corona height h. The authors report that only the models with (a*=0, fcol=1.7) and (a*=1, fcol=2.4) can fit the data, with best-fit heights h=42^{+34}_{-23} R_g and h=50^{+15}_{-12} R_g respectively, both with chi^2_min/dof = 11.4/5 and p_null = 0.04. The central claim is that the observed quasar time lags are fully consistent with the X-ray reverberation hypothesis, provided the corona height is larger than about 40 R_g, and that this is consistent with microlensing disc-size constraints.
Significance. If substantiated, the result would be a notable step: it would show that standard Novikov-Thorne discs illuminated by a compact X-ray corona can reproduce the UV/optical time lags in luminous quasars, potentially resolving the apparent discrepancy between observed and predicted disc sizes without invoking non-standard disc physics. The study uses a large archival sample and a physically motivated model that includes relativistic effects and radial-dependent external heating, and it explicitly connects to independent microlensing constraints. However, the statistical support is marginal (p_null = 0.04, below the conventional 0.05 threshold), and the headline lower limit on the corona height is not supported by the confidence intervals reported in Table 1. The analysis also relies on an arbitrary error-inflation factor and on the applicability of K21b relations outside their calibrated parameter range.
major comments (4)
- [Abstract; Sect. 5.2; Table 1] The claim that 'the corona height is larger than ~40 Rg' is not supported by the best-fit uncertainties reported in Table 1. For model M01.7 the best fit is h=42^{+34}_{-23} Rg, so h≈19 Rg is within the 1σ lower interval; for model M12.4 h=50^{+15}_{-12} Rg, so h≈38 Rg is within 1σ. Thus the data are consistent with heights well below 40 Rg, and the statements in the abstract and in Sect. 5.2 ('the corona height should be equal or larger than 40 Rg') are overstatements. Please quote the confidence intervals and rephrase the conclusion as a best-fit range rather than a firm lower limit.
- [Sect. 4; Table 1] The reported p_null = 0.04 for the two acceptable models is below the conventional 0.05 threshold for rejecting the null hypothesis. The text nevertheless describes these fits as fitting 'well' and later concludes that the data are 'fully consistent' with the model. This is an internal inconsistency in the statistical interpretation. The authors should acknowledge that the fits are marginal at the 5% level, report the exact p-values, and discuss how the conclusion depends on the choice of significance threshold.
- [Sect. 4] The error-inflation factor of 1.1 applied to the observed mean time-lag errors is ad hoc and is not a propagation of the stated uncertainties on MBH, Lbol, and L2-10keV. Since the chi^2 values are near the rejection threshold, the choice of this factor directly affects whether a model is deemed acceptable. The authors should either propagate the input parameter uncertainties through the model to obtain model time-lag uncertainties, or at minimum perform a sensitivity analysis varying the inflation factor and reporting how chi^2_min and p_null change.
- [Sect. 2.2] The K21b analytic time-lag relations were developed for AGN with MBH <~ 1e8 Msun and accretion rate <~ 0.5, but the final sample extends to MBH < 1e9 Msun and lambda_Edd <= 1. The text states that the analysis was repeated on smaller subsamples with MBH < 5e8 and 2e8 Msun and that the results were 'almost identical', but no numerical results are provided. This robustness test is load-bearing for the extrapolation claim; please present the best-fit heights, chi^2 values, and p-values for these restricted samples, or justify why the K21b relations remain valid in the extrapolated regime.
minor comments (5)
- [Abstract] The abstract states that the model assumes 'the measured BH mass, accretion rate and X-ray luminosity', but the 2-10 keV luminosity is not measured; it is estimated from the Lusso et al. (2012) relation. Please clarify this in the text.
- [Section 5] The 'Conclusions' heading is followed by no text. A concise conclusions paragraph summarizing the main findings and their caveats should be added.
- [Table 1] For models M01 and M11 the best-fit height is listed as 80 Rg, which is the upper boundary of the searched range. Please state explicitly whether the chi^2 minimum is at the boundary and whether the fit was attempted beyond 80 Rg.
- [Sect. 4] The description 'we increased the error of \bar t_lag,obs by a factor of 1.1' is ambiguous: it should state whether the factor is applied to the standard error of the mean or to the individual time-lag errors before averaging.
- [Sect. 4; Fig. 7] The Kolmogorov–Smirnov test result is quoted as 'pnull = 0.15'. Please define this quantity explicitly and ensure the notation is consistent with the p_null used for the chi^2 fits.
Circularity Check
No significant circularity: the corona height is a fitted parameter and the consistency claim is a goodness-of-fit against independent archival data.
full rationale
The paper's central derivation chain is not circular. The observed mean time-lags spectrum is constructed from archival measurements (Jiang et al. 2017; Homayouni et al. 2019; Guo et al. 2022; Jha et al. 2022), which are external to the authors. The model time lags are computed from the K21a/K21b relativistic response functions and the Kammoun et al. (2023) colour-correction prescription; these are previously published physical calculations, not fits to the quasar time lags used here. The only free parameter is the X-ray corona height h, scanned over 5-80 Rg, with χ² minimized against the six binned mean lags. Thus the statement that X-ray reverberation is 'fully consistent' with the observed lags is a goodness-of-fit result with one free parameter, not an input recycled as an output. Heavy self-citation to K21a/b, Kammoun et al., Dovčiak et al., and Papadakis et al. (2022) is present, but it is not load-bearing in a circular way: the K21b equations are used as an independent model, and the Papadakis et al. microlensing comparison is a consistency check across separate observables (disc half-light radii vs. continuum time lags). No equation in the paper defines the target result in terms of itself, and no fitted parameter is renamed as a prediction. The skeptic's concern that the h > 40 Rg abstract claim is not supported by the 1σ intervals (e.g., h = 42^{+34}_{-23} Rg allowing h ≈ 19 Rg) is a statistical calibration or overstatement issue, not a circularity of the derivation chain.
Assumptions & free parameters
free parameters (2)
- X-ray corona height h =
42 Rg (a*=0, fcol=1.7); 50 Rg (a*=1, fcol=2.4)
- Error inflation factor =
1.1
assumptions (5)
- domain assumption K21b analytic time-lag relations remain valid for MBH up to 1e9 Msun and lambdaEdd up to 1, despite being developed for MBH <= 1e8 Msun and accretion rate < 0.5.
- domain assumption X-ray corona is a point source on the BH rotation axis at height h, emitting isotropically in its rest frame.
- domain assumption The accretion disc is a standard, geometrically thin, flat Novikov-Thorne disc, with a single colour-correction factor applied to the whole disc.
- domain assumption 2-10 keV luminosity for every source follows the Lusso et al. (2012) relation, log(Lbol/L2-10keV) = 0.75 log(lambdaEdd) + 2.13.
- domain assumption All quasars in the sample share the same corona height, black hole spin, and colour-correction factor.
Cite this review
Pith. "Pith review of X-ray reverberation modelling of the continuum, optical/UV time-lags in quasars." pith.science (2026). https://pith.science/paper/H72H3HZT
@misc{pith2026241109681,
author = {Pith},
title = {Pith review of: X-ray reverberation modelling of the continuum, optical/UV time-lags in quasars},
year = {2026},
howpublished = {\url{https://pith.science/paper/H72H3HZT}},
note = {Machine review of arXiv:2411.09681}
}
abstract
Context: Extensive, multi-wavelength monitoring campaigns of nearby and higher redshift active galactic nuclei (AGN) have shown that the UV/optical variations are well correlated with time delays which increase with increasing wavelength. Such behaviour is expected in the context of the X-ray thermal reverberation of the accretion disc in AGN. Aims: Our main objective is to use time-lag measurements of luminous AGN and fit them with sophisticated X-ray reverberation time-lags models. In this way we can investigate whether X-ray reverberation can indeed explain the observed continuum time lags, and whether time-lag measurements can be used to measure physical parameters such as the X-ray corona height and the spin of the black hole (BH) in these systems. Methods: We use archival time-lag measurements for quasars from different surveys, and we compute their rest frame, mean time-lags spectrum. We fit the data with analytical X-ray reverberation models, using $\chi^2$ statistics, and fitting for both maximal and non spinning BHs, for various colour correction values and X-ray corona heights. Results: We found that X-ray reverberation can explain very well the observed time lags, assuming the measured BH mass, accretion rate and X-ray luminosity of the quasars in the sample. The model agrees well with the data both for non-rotating and maximally rotating BHs, as long as the corona height is larger than $\sim 40$ gravitational radii. This is in agreement with previous results which showed that X-ray reverberation can also explain the disc radius in micro-lensed quasars, for the same corona heights. The corona height we measure depends on the model assumption of a perfectly flat disc. More realistic disc models may result in lower heights for the X-ray corona.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 2 Pith papers
-
X-ray disc reverberation modelling of the X-ray/UV/optical spectral/timing properties of Fairall 9
X-ray reprocessing by the accretion disc in Fairall 9 simultaneously explains the SED, UV/optical power spectra, and interband time lags only if the corona is powered externally rather than by accretion.
-
Continuum optical-UV and X-ray variability of AGN: current results and future challenges
A comprehensive review of AGN continuum variability from optical/UV to X-rays, with no new data.
Reference graph
Works this paper leans on
-
[1]
M., Chiang, C.-Y ., McHardy, I., et al
Cackett, E. M., Chiang, C.-Y ., McHardy, I., et al. 2018, ApJ, 857, 53
2018
-
[2]
M., Gelbord, J., Li, Y .-R., et al
Cackett, E. M., Gelbord, J., Li, Y .-R., et al. 2020, ApJ, 896, 1
2020
-
[3]
M., Horne, K., & Winkler, H
Cackett, E. M., Horne, K., & Winkler, H. 2007, MNRAS, 380, 669
2007
-
[4]
2019, Nature As- tronomy, 3, 251
Chelouche, D., Pozo Nuñez, F., & Kaspi, S. 2019, Nature As- tronomy, 3, 251
work page 2019
-
[5]
Davis, S. W. & El-Abd, S. 2019, ApJ, 874, 23
2019
-
[6]
W., Jin, C., Blaes, O., & Ward, M
Done, C., Davis, S. W., Jin, C., Blaes, O., & Ward, M. 2012, MNRAS, 420, 1848
2012
-
[7]
Donnan, F. R., Hernández Santisteban, J. V ., Horne, K., et al. 2023, MNRAS, 523, 545 Dovˇciak, M., Papadakis, I. E., Kammoun, E. S., & Zhang, W. 2022, A&A, 661, A135
work page 2023
- [8]
Show all 57 references
-
[9]
2019, ApJ, 870, 123
Edelson, R., Gelbord, J., Cackett, E., et al. 2019, ApJ, 870, 123
2019
-
[10]
M., Horne, K., et al
Edelson, R., Gelbord, J. M., Horne, K., et al. 2015, ApJ, 806, 129
2015
-
[11]
M., Denney, K
Fausnaugh, M. M., Denney, K. D., Barth, A. J., et al. 2016, ApJ, 821, 56
2016
-
[12]
E., et al
Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748
1996
-
[13]
Gaskell, C. M. & Peterson, B. M. 1987, ApJS, 65, 1
1987
-
[14]
Gaskell, C. M. & Sparke, L. S. 1986, ApJ, 305, 175
1986
-
[15]
J., & Wang, S
Guo, H., Barth, A. J., & Wang, S. 2022, ApJ, 940, 20
2022
-
[16]
B., Sarrouh, G
Hall, P. B., Sarrouh, G. T., & Horne, K. 2018, ApJ, 854, 93
2018
-
[17]
R., Millman, K
Harris, C. R., Millman, K. J., Van Der Walt, S. J., et al. 2020, Nature, 585, 357 Hernández Santisteban, J. V ., Edelson, R., Horne, K., et al. 2020, MNRAS, 498, 5399
2020
-
[18]
R., Trump, J
Homayouni, Y ., Sturm, M. R., Trump, J. R., et al. 2022, ApJ, 926, 225
2022
-
[19]
R., Grier, C
Homayouni, Y ., Trump, J. R., Grier, C. J., et al. 2019, ApJ, 880, 126
2019
-
[20]
Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90
2007
-
[21]
K., Joshi, R., Chand, H., et al
Jha, V . K., Joshi, R., Chand, H., et al. 2022, MNRAS, 511, 3005
2022
-
[22]
J., Greene, J
Jiang, Y .-F., Green, P. J., Greene, J. E., et al. 2017, ApJ, 836, 186
2017
-
[23]
E., Dov ˇciak, M., & Panagiotou, C
Kammoun, E., Papadakis, I. E., Dov ˇciak, M., & Panagiotou, C. 2024, A&A, 686, A69
2024
-
[24]
S., Papadakis, I
Kammoun, E. S., Papadakis, I. E., & Dov ˇciak, M. 2019, ApJ, 879, L24
2019
-
[25]
2023, MNRAS, 526, 138
Panagiotou, C. 2023, MNRAS, 526, 138
2023
-
[26]
J., Cackett, E
Kara, E., Barth, A. J., Cackett, E. M., et al. 2023, ApJ, 947, 62
2023
-
[27]
A., et al
Kara, E., Mehdipour, M., Kriss, G. A., et al. 2021, ApJ, 922, 151
2021
-
[28]
& Nayakshin, S
Kazanas, D. & Nayakshin, S. 2001, ApJ, 550, 655
2001
-
[29]
C., Bechtold, J., & Siemiginowska, A
Kelly, B. C., Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895
2009
-
[30]
C., Becker, A
Kelly, B. C., Becker, A. C., Sobolewska, M., Siemiginowska, A., & Uttley, P. 2014, ApJ, 788, 33
2014
-
[31]
Korista, K. T. & Goad, M. R. 2001, ApJ, 553, 695
2001
-
[32]
Korista, K. T. & Goad, M. R. 2019, MNRAS, 489, 5284
2019
-
[33]
J., Trakhtenbrot, B., Ricci, C., et al
Koss, M. J., Trakhtenbrot, B., Ricci, C., et al. 2022, ApJS, 261, 1
2022
-
[34]
C., Papadakis, I
Kumari, K., Dewangan, G. C., Papadakis, I. E., et al. 2023, MN- RAS, 521, 4109
2023
-
[35]
C., Papadakis, I
Kumari, K., Dewangan, G. C., Papadakis, I. E., & Singh, K. P. 2024, MNRAS, 527, 5668
2024
-
[36]
R., Korista, K
Lawther, D., Goad, M. R., Korista, K. T., Ulrich, O., & Vester- gaard, M. 2018, MNRAS, 481, 533
2018
-
[37]
D., et al
Lusso, E., Comastri, A., Simmons, B. D., et al. 2012, MNRAS, 425, 623
2012
-
[38]
M., Cameron, D
McHardy, I. M., Cameron, D. T., Dwelly, T., et al. 2014, MN- RAS, 444, 1469
2014
-
[39]
A., Cackett, E
Miller, J. A., Cackett, E. M., Goad, M. R., et al. 2023, ApJ, 953, 137
2023
-
[40]
2018, ApJ, 862, 123
Mudd, D., Martini, P., Zu, Y ., et al. 2018, ApJ, 862, 123
2018
-
[41]
2022, MNRAS, 509, 2637
Netzer, H. 2022, MNRAS, 509, 2637
2022
-
[42]
Novikov, I. D. & Thorne, K. S. 1973, 343
1973
-
[43]
2022, ApJ, 935, 93
Panagiotou, C., Papadakis, I., Kara, E., Kammoun, E., & Dovˇciak, M. 2022, ApJ, 935, 93
2022
-
[44]
E., Kammoun, E
Panagiotou, C., Papadakis, I. E., Kammoun, E. S., & Dov ˇciak, M. 2020, MNRAS, 499, 1998
2020
-
[45]
E., Dov ˇciak, M., & Kammoun, E
Papadakis, I. E., Dov ˇciak, M., & Kammoun, E. S. 2022, A&A, 666, A11
2022
-
[46]
Peterson, B. M. 2004, in The Interplay Among Black Holes, Stars and ISM in Galactic Nuclei, ed. T. Storchi-Bergmann, L. C. Ho, & H. R. Schmitt, V ol. 222, 15–20
2004
-
[47]
Bauer, F. E. 2024, A&A, 686, A286
2024
-
[48]
T., Lacy, M., Storrie-Lombardi, L
Richards, G. T., Lacy, M., Storrie-Lombardi, L. J., et al. 2006, ApJS, 166, 470
2006
-
[49]
R., Fabian, A
Ross, R. R., Fabian, A. C., & Mineshige, S. 1992, MNRAS, 258, 189
1992
-
[50]
A., Huang, J., Horne, K., & Lin, D
Starkey, D. A., Huang, J., Horne, K., & Lin, D. N. C. 2023, MN- RAS, 519, 2754 Article number, page 9 of 10 A&A proofs: manuscript no. timelags_no_bold
2023
-
[51]
L., Stubbs, C
Tonry, J. L., Stubbs, C. W., Lykke, K. R., et al. 2012, ApJ, 750, 99
2012
-
[52]
M., McHardy, I., Cackett, E
Vincentelli, F. M., McHardy, I., Cackett, E. M., et al. 2021, MN- RAS, 504, 4337
2021
-
[53]
E., et al
Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261
2020
-
[54]
White, R. J. & Peterson, B. M. 1994, PASP, 106, 879
1994
-
[55]
M., et al
Yu, Z., Martini, P., Davis, T. M., et al. 2020, ApJS, 246, 16
2020
-
[56]
S., Kozłowski, S., & Peterson, B
Zu, Y ., Kochanek, C. S., Kozłowski, S., & Peterson, B. M. 2016, ApJ, 819, 122
2016
-
[57]
S., & Peterson, B
Zu, Y ., Kochanek, C. S., & Peterson, B. M. 2011, ApJ, 735, 80 Article number, page 10 of 10
2011
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.