{"id":"b14ba656-98b9-4ffc-8446-4479139e7fa3","arxiv_id":"2509.03159","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Using X-ray reverberation models, the authors reproduce the SDSS Stripe-82 quasar UV/optical power spectra at 1300-4000 A, favouring an accretion-powered corona, spin below 0.7, and corona height 20-60 Rg.","lead":"Astronomers tested whether X-ray light from a hot corona reflecting off the accretion disc can explain how quasars flicker in ultraviolet and optical light. The model fits the observed flickering well, and implies the corona is powered by the disc itself, with black hole spin below 0.7 and corona height 20 to 60 gravitational radii.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Predicted UV/optical PSD amplitude scales linearly with the assumed X-ray PSD normalization (Eq. 10); if the local-Seyfert A=6e-3 mdot^-0.8 is off by a factor of a few, the case-B rejection and h/spin inferences are not secure.","rationale":"The paper's central claim is that X-ray reverberation can fit the observed UV/optical PSDs, and that the best fits imply an accretion-powered corona, spin <0.7, and height 20-60 Rg. The model predictions are physically coherent, and the authors are transparent about their simplifications, but the chain from assumption to inference has one clearly weakest link: the X-ray PSD is not measured for these quasars. Because Eq. (10) is linear in PSD_X, the amplitude of every predicted UV/optical PSD is set by the adopted normalization A and the X-ray luminosity factor. The case-B rejection is an amplitude statement, and even a modest systematic error in A (well within the scatter of the local-Seyfert relation and the redshift/luminosity extrapolation) could change the conclusion. The reader's verdict correctly identifies this as the weakest assumption; my read does not move the verdict. I would keep the paper as CONDITIONALLY acceptable: the framework and fits are useful, but the headline physical inferences (especially the rejection of case B, and the h/spin values) should be treated as provisional until the X-ray PSD sensitivity is quantified or actual X-ray variability constraints for these quasars are incorporated.","tokens_in":21442,"tokens_out":11426,"duration_ms":137642,"concrete_test":"Re-run the global fits with the X-ray PSD normalization A free (one global or per-bin), and compare the best-fit A for case B to the Ponti et al. scatter; also check whether case-A h and spin move by >30% or change ranking.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claims are conditional on the assumed X-ray PSD. In Eq. (10), PSD_lambda,mod = |Gamma_norm|^2 * PSD_X * L_X^2, so the UV/optical PSD amplitude is directly proportional to the X-ray PSD normalization A from Eq. (8), which is taken from local Seyferts (Ponti et al. 2012). The bend frequency from McHardy et al. (2006) lies above the observed P24 frequency band for these quasars, so the shape is mostly set by the -1 low-frequency slope; the amplitude is the main lever. The rejection of the externally powered corona (Sect. 7.2) is explicitly an amplitude argument: |Gamma_norm|^2 is said to be more than 10 times lower than in case A. However, the full model also contains L_X^2; in case B, larger |Ltransf/Ldisc| increases L_X^2 while |Gamma_norm|^2 decreases, so the net model amplitude is not shown and could be less suppressed than claimed. Even setting that aside, if the true quasar X-ray PSD normalization is a factor of 2-3 higher (plausible given the scatter in the local-Seyfert relation and the extrapolation to z~1-2 quasars), the case-B model could become viable and the inferred Ltransf/Ldisc, corona height, and spin would shift. The authors acknowledge this in Sect. 8 but do not quantify the sensitivity of their headline physical parameters to this input.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes accretion-disc response/transfer functions for X-ray reprocessing using the KYNXiltr/KYNSED models (K23, Dovčiak et al. 2022), studies their dependence on wavelength, Ltransf/Ldisc, BH spin, and fcol, and then fits UV/optical power spectra predicted from an assumed X-ray bending-power-law PSD to the ensemble quasar PSDs of Petrecca et al. (2024). The claim is that X-ray reverberation can fit the observed quasar PSDs from 1300–4000 Å, but only for an accretion-powered corona, with best-fit parameters Ltransf/Ldisc ~ 0.7–0.8, h ~ 20–60 Rg, fcol ~ 2, and spin < 0.7. The central formalism is standard and the grid fitting is transparent, but the physical conclusions are conditional on an assumed X-ray PSD normalization calibrated on local Seyferts, and the statistical treatment of the P24 PSD points as independent data is questionable.","tokens_in":21832,"tokens_out":5279,"duration_ms":61557,"significance":"If the result holds, it would be an important step: the apparently universal UV/optical PSD shape of quasars could be interpreted physically as X-ray disc reprocessing, and the inferred corona height, spin, and coronal power source would be directly relevant to AGN accretion physics. The paper uses a state-of-the-art radiative-transfer code, computes transfer functions from physical inputs rather than fitting them, and is transparent about its assumptions and limitations. However, the headline constraints are only as secure as the adopted X-ray PSD, which is not measured for these quasars; the discrimination between the two coronal power scenarios is essentially an amplitude argument, and the statistical framework treats strongly correlated, derived PSD points as independent. The qualitative agreement is encouraging, but the quantitative parameter claims need additional robustness analysis.","major_comments":[{"comment":"The rejection of the externally powered corona is presented as an amplitude argument, but §7.2 compares only |Γ_norm|^2 and states it is “more than 10 times lower” in case B. The actual model PSD is PSD_{λ,mod} = |Γ_norm|^2 · PSD_X · L_X^2 (Eq. 10). In case B, increasing |Ltransf/Ldisc| increases L_X^2 while decreasing |Γ_norm|^2, so the net amplitude may be substantially less suppressed than the transfer-function amplitude alone. Please report the net product |Γ_norm|^2·L_X^2 over the allowed case-B grid, and show explicitly that even including L_X^2 the case-B model cannot reach the observed PSD amplitude.","section":"§7.2, Eq. (10)"},{"comment":"The X-ray PSD of the target quasars is unknown and is taken from local Seyfert scalings: A = 6×10^{-3} ṁ^{-0.8} (Eq. 8) and the McHardy et al. (2006) bend relation (Eq. 9). Since Eq. (10) is linear in PSD_X, the predicted UV/optical amplitude scales linearly with A. A factor-of-2–3 change in A—well within the scatter of Ponti et al. (2012) and the extrapolation from local Seyferts to z~1–2 quasars—rescales all model amplitudes and could shift the inferred Ltransf/Ldisc, h, fcol, and spin, and could alter the case-B conclusion. The authors acknowledge the unknown X-ray PSD in §8, but they do not quantify the sensitivity of the headline parameters to A. Please add a sensitivity test, e.g. refit with A varied by its 1σ scatter or by factors of 2–3.","section":"§5.1, Eqs. (7)–(10), §8"},{"comment":"The three rest-frame frequencies used in the χ² fit are not independent measurements: they derive from P24’s two-parameter power-law fit, log PSD(ν) = log PSD_amp + PSD_slope[log ν + 2.6], so the PSD values at the three frequencies are perfectly correlated functions of the fitted amplitude and slope. Treating them as 18 independent points per (MBH, ṁEdd) bin in Eq. (12) overstates the number of constraints and affects the reported χ² and the parameter uncertainties. Please fit directly in the (PSD_amp, PSD_slope) space, propagate the P24 covariance matrix, or otherwise justify the independence assumption.","section":"§6, Eq. (12)"}],"minor_comments":[{"comment":"The fcol grid contains only 1, 1.7, and 2.5, but the text reports a mean best-fit fcol ~ 2. Please clarify whether the reported value is an interpolation between grid points or a grid-edge/integrated estimate, and indicate the resulting systematic uncertainty.","section":"§7.1, Table 1"},{"comment":"The aliasing correction uses a single rest-frame sampling rate and a simplified cutoff at (1+z)/90 day^{-1}, although the P24 light curves are irregularly binned. The ad hoc 0.1 dex model error may not fully cover this simplification. A short discussion or test of the sensitivity of the fits to the aliasing prescription would strengthen the conclusions.","section":"§5.3"},{"comment":"The statement that the dependence of |Γ_norm|^2 on fcol “becomes weaker (and even non-existent in the case B corona)” is based on λ=3000 Å and the fiducial parameters; please state explicitly that this is parameter-dependent, since the text just before it describes a general fcol dependence.","section":"§3.4.5, Fig. 2"},{"comment":"There are numerous typographical issues: “di fferent”, “˙mEdd” rendering, “Sect. 8.1” equation references, and the use of “K21a/K23/P22” without expanding the first time in the abstract or introduction. A careful language and notation pass is needed.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A and the central approach is not circular: the transfer functions are computed from physical radiative-transfer inputs, not from the fitted PSDs. The main concerns are statistical and sensitivity-related: the case-B rejection argument in §7.2 omits the L_X^2 factor in Eq. (10), and the headline parameters depend linearly on an assumed X-ray PSD normalization. The third major comment (correlated P24 points) is also important for the credibility of the χ² values. These are fixable within the manuscript's scope, hence major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is fitting the K23 accretion-powered transfer functions to P24's Stripe-82 ensemble PSDs at six wavelengths, and using that to distinguish an accretion-powered corona (case A) from an externally powered one (case B). The case A/B distinction is new, and the resulting spin <0.7 and height 20-60 Rg are physically interesting. The method is transparent: standard convolution formalism, a clearly described grid, and a simple chi-square statistic. They also check competing contributions (intrinsic disc variability, Balmer continuum, TDEs) and find them subdominant. That is good practice, and the fits themselves look credible (chi2 around 110 for 135 dof).\n\nThe soft spot is exactly where the reader put it. The X-ray PSD is not measured; it comes from local Seyfert relations (Eqs. 7-9), and Eq. 10 makes the whole UV/optical PSD amplitude proportional to that normalization A. The case B rejection in Sect. 7.2 is an amplitude argument, and the spin/height constraints inherit the same dependence. The authors acknowledge this in Sect. 8 but never quantify how much A would have to change to alter the conclusions. A factor of 2-3 in A is plausible given scatter in the Ponti et al. relation and the extrapolation to z~1-2 quasars. The stress-test note about L_X^2 partly cancelling is less damaging than it first appears, because they actually fit case B with the full Eq. 10 rather than comparing |Gamma_norm|^2 alone; still, the full model amplitude is not shown, and sensitivity to A remains unquantified. Two smaller issues: the best-fit parameters come with no confidence intervals, and the three frequencies per wavelength are treated as independent even though they come from P24's two-parameter power-law fits and are correlated. The code is not shipped either, though that is a common frustration rather than a fatal one.\n\nOverall this is a solid, honest paper that deserves a serious referee. The referee should push for a sensitivity analysis on A and for confidence intervals on the fitted parameters. I would not desk-reject it; I would send it out with a request for conditional revision.","headline":"Worth reading: a transparent grid fit of X-ray reverberation to ensemble quasar PSDs, but the headline spin/height constraints and the case-A-vs-B rejection all rest on an adopted X-ray PSD that the paper does not stress-test.","tokens_in":22350,"tokens_out":4610,"would_cite":true,"duration_ms":49943,"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":"X-ray reverberation can explain quasar UV/optical flickering across 1300–4000 Å, provided the X-ray corona is powered by the accretion flow.","keywords":["accretion discs","active galactic nuclei","X-ray reverberation","quasar variability","power spectra","corona geometry","black hole spin","lamp-post model"],"falsifier":"Measure the X-ray power spectrum of a Stripe-82 quasar with the same black hole mass and accretion rate as the sample bins; if its break frequency or normalization deviates from the assumed Seyfert scaling relations by more than the adopted 0.1 dex uncertainty, the predicted UV/optical PSD amplitude would shift linearly and the spin and height constraints would no longer hold. A second test: search the UV/optical PSD of a lower-mass AGN for the predicted high-frequency flattening feature (around ν ≳ 0.1 day⁻¹ for M_BH = 8×10⁸ M☉); its absence would rule out the accretion-powered corona signatu","tokens_in":21350,"feed_emoji":"🔭","tokens_out":2565,"duration_ms":29451,"temperature":0.7,"pith_summary":"This paper tests whether the observed ultraviolet and optical variability of distant quasars can be produced entirely by X-ray illumination of the accretion disc—the X-ray reverberation scenario. Using a lamp-post corona model and observed ensemble power spectra of SDSS Stripe-82 quasars, the authors claim that reverberation alone fits the power spectra well at all wavelengths from 1300 Å to 4000 Å. The fits require that the corona be powered by energy extracted from the accretion disc, not by an external source; an externally powered corona under-predicts the variability amplitude by more than an order of magnitude. The best fits also imply black hole spins lower than about 0.7 and a compact X-ray source sitting 20–60 gravitational radii above the disc.","feed_headline":"Quasar UV/optical flickering traced to X-ray echoes off the disc","feed_subtitle":"Reverberation fits power spectra from 1300 to 4000 Å; best fits demand an accretion-powered corona and spin below 0.7.","key_machinery":"The disc transfer function |Γλ(ν)|², computed as the Fourier transform of the disc response function Ψλ(t), is the central object. It converts an assumed X-ray power spectrum into a predicted UV/optical power spectrum through PSDλ(ν) = |Γλ,norm(ν)|² · PSDX(ν) · L²_X,Edd. The response functions come from the KYNXiltr code, which incorporates general-relativistic light bending, disc reflection, and a distinction between accretion-powered (case A) and externally powered (case B) coronae. The transfer function's frequency shape—flat at low frequencies, bending down at high frequencies, with an additional flattening feature for case A—is what allows the model to match the observed power spectra.","core_discovery":"Under the assumption that each quasar's X-ray power spectrum follows a bending power law calibrated on local Seyferts, the authors show that the predicted UV/optical power spectrum, computed as the disc transfer function times the X-ray power spectrum, reproduces the observed quasar power spectra across six rest-frame wavelengths. The central quantitative result is the rejection of an externally powered corona: in that case the transfer function normalized to the disc flux is more than ten times too small to match the observed variability. Only an accretion-powered corona, with roughly 70–80 percent of the inner accretion power transferred to the X-ray source, matches the data. The best-fit","pith_inferences":["A natural extension is to test the model on individual quasars with measured X-ray light curves, rather than ensemble-averaged spectra, which would remove the main systematic uncertainty and directly check the predicted UV/optical lag and PSD amplitude.","The result that an externally powered corona under-predicts variability by more than an order of magnitude suggests that any successful alternative to the lamp-post model must also produce a large reprocessed-to-intrinsic disc flux ratio in the UV/optical bands.","The spin constraint (α* < 0.7) is conditional on the assumed X-ray PSD normalization; if quasar X-ray PSDs are systematically more variable than Seyfert scaling relations predict, the spin and height estimates would shift.","Because the model assumes a single fixed spin for all quasars in the sample, a more realistic treatment allowing a spin distribution would change the inferred parameters and likely improve the scatter in the residuals; this could be tested with mock quasar populations."],"forward_implications":["If correct, the observed UV/optical variability of quasars in the Stripe-82 sample is dominated by X-ray reprocessing rather than intrinsic disc fluctuations, since adding an intrinsic −1 slope power-law component does not significantly improve the fits.","The inferred corona height of 20–60 Rg and spin below 0.7 become physical constraints on the accretion geometry of luminous quasars, consistent with earlier micro-lensing disc-size and time-lag results.","The high-frequency flattening feature predicted for an accretion-powered corona offers a direct observational test: it should appear in the UV/optical PSDs of lower-mass AGNs at frequencies above the Poisson noise level.","The model reinforces the 'universal PSD shape' idea: after rescaling by the gravitational timescale, the predicted power spectra for different black hole masses converge to a common shape.","Future dense long-baseline monitoring (e.g., LSST, ULTRASAT, UVEX) could extend these power spectra over a wide frequency range and sharpen the constraints on corona height, spin, and the accretion-power transfer fraction."],"supporting_citations":[{"why":"Provides the observed ensemble UV/optical power spectra of SDSS Stripe-82 quasars that the model is fitted against.","marker":"P24"},{"why":"Supplies the KYNXiltr code used to compute disc response functions for accretion-powered and externally powered coronae.","marker":"K23"},{"why":"Defines the disc response function formalism and the baseline illuminated-disc model that K23 extends.","marker":"K21a"},{"why":"Establishes the method of computing disc transfer functions and modelling UV/optical PSDs via X-ray reverberation, applied here to quasars.","marker":"P22"},{"why":"Gives the empirical relation between X-ray PSD break frequency, black hole mass, and bolometric luminosity used to set the X-ray power spectrum shape.","marker":"McHardy et al. 2006"},{"why":"Provides the normalization of the assumed X-ray PSD as a function of Eddington ratio.","marker":"Ponti et al. 2012"},{"why":"Previous application of the X-ray reverberation model to the quasar micro-lensing disc size problem, cited as independent support for the geometry inferred here.","marker":"Papadakis et al. 2022"},{"why":"Shows that observed UV/optical time lags in a large quasar sample can be explained by X-ray reverberation, supporting the same scenario.","marker":"Langis et al. 2024"}],"fun_headline_variants":["X-ray echoes explain quasar UV/optical flickering from 1300 to 4000 Å","Quasar UV/optical variability traced to X-ray disc reverberation","X-ray reverberation model fits quasar power spectra, implies spin <0.7","Accretion-powered X-ray corona best matches quasar flickering"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The predicted UV/optical power spectrum is directly proportional to an assumed X-ray power spectrum that is calibrated from nearby Seyferts, not measured for these quasars; if the true quasar X-ray PSD differs in normalization or shape, the inferred corona height, spin, and rejection of the external corona would change.","fun_headline_variants_meta":{"raw":{"variants":["X-ray echoes explain quasar UV/optical flickering from 1300 to 4000 Å","Quasar UV/optical variability traced to X-ray disc reverberation","X-ray reverberation model fits quasar power spectra, implies spin <0.7","Accretion-powered X-ray corona best matches quasar flickering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000575,"raw_usage":{"total_tokens":2613,"prompt_tokens":868,"completion_tokens":1745,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":1671}},"tokens_in":612,"tokens_out":1745,"duration_ms":14471,"temperature":1.0,"reasoning_tokens":1671,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:07:05.134794+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the X-ray power spectrum of a Stripe-82 quasar with the same black hole mass and accretion rate as the sample bins; if its break frequency or normalization deviates from the assumed Seyfert scaling relations by more than the adopted 0.1 dex uncertainty, the predicted UV/optical PSD amplitude would shift linearly and the spin and height constraints would no longer hold. A second test: search the UV/optical PSD of a lower-mass AGN for the predicted high-frequency flattening feature (around ν ≳ 0.1 day⁻¹ for M_BH = 8×10⁸ M☉); its absence would rule out the accretion-powered corona signatu","supporting_citations":[{"cited_title":"2012, A&A, 542, A83","cited_arxiv_id":null,"evidence_quote":"Provides the normalization of the assumed X-ray PSD as a function of Eddington ratio."},{"cited_title":"E., Dovˇciak, M., & Kammoun, E","cited_arxiv_id":null,"evidence_quote":"Previous application of the X-ray reverberation model to the quasar micro-lensing disc size problem, cited as independent support for the geometry inferred here."},{"cited_title":"A., Papadakis, I","cited_arxiv_id":null,"evidence_quote":"Shows that observed UV/optical time lags in a large quasar sample can be explained by X-ray reverberation, supporting the same scenario."}],"review_version":1}