{"id":"22a083c3-95d0-47df-a407-12d30d9df639","arxiv_id":"2502.08366","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":17,"one_line_summary":"New optical reverberation lags in 3C 273 show an accretion disc 2-7 times larger than thin-disc predictions, with extrapolated dusty-disc radii of 100-200 light-days matching the broad-line region.","lead":"Astronomers measured how long it takes light from different parts of the quasar 3C 273's accretion disc to reach Earth, using seven optical filters over three years. The observed delays imply the disc is about 2 to 7 times larger than standard thin-disc theory predicts, and it may extend far enough to host the dusty gas that forms the broad-line region.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Factor ~2–7 disc-size excess is not robust to the adopted black-hole mass: Li et al. 2022 mass changes predicted lags by ~1.56, dropping the X=4.96 discrepancy to ~1.2–1.6.","rationale":"I read the paper as a genuine first RM campaign on 3C273 with two independent lag estimators, and the FVG and structure-function checks give independent support for disc-dominated variability. The most central quantitative assertion, however, is the factor ~2–7 excess over thin-disc predictions. That factor is computed against Eq. 5, which is sensitive to M_BH and mdot_Edd. The authors' statement that the Li et al. mass change 'roughly cancels' is not correct under the paper's own scaling: for fixed L_acc, tau0 ∝ M^{1/3}, so the predicted lags shift by ~56%. Recomputing the Table 5 ratios with M_BH = 1.15e9 M_sun would change the X=4.96 comparison from a clear excess to roughly 1.2–1.6, which is within or near the uncertainties for several bands. Even with X=2.49, the claimed range shrinks from ~4–7 to ~3–4. This is not a disagreement with external consensus; it is a sensitivity of the paper's own formula to a value it cites, and it should be settled by a one-line recalculation. The near-IR/BLR extrapolation flagged by the reader is also fragile, but it is downstream of the size-problem claim: the mass sensitivity is the more load-bearing issue for the paper's principal result. The reader's CONDITIONAL verdict already includes a request to quantify black-hole mass and X uncertainties, so my read does not change that verdict; it sharpens the reason the condition matters.","tokens_in":29347,"tokens_out":10855,"duration_ms":121257,"concrete_test":"Recompute the Table 5 predicted lags using Eq. 5 with M_BH = 1.15e9 M_sun and mdot_Edd = 0.31 (using L_acc = 4.5e46 erg/s and L_Edd = 1.45e47 erg/s), keeping X and lambda0 fixed, and multiply the predicted differential lags by 1.56. Compare the resulting predicted lags with the Year3 PyROA and Javelin lags for g, V, r, and z. If the X=4.96 ratios fall to ~1.2–1.6 while the X=2.49 ratios fall to ~3–4, the paper should report a mass- and X-dependent excess of roughly 1–4 rather than the stated factor ~2–7, and the size-problem conclusion should be qualified accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central size-problem claim rests on the predicted differential lags from Eq. 5, which is evaluated in Section 4.1 with M_BH = 3e8 M_sun and mdot_Edd = 1.2. The paper dismisses the Li et al. (2022) mass M_BH = 1.15e9 M_sun by saying the higher Eddington luminosity and lower Eddington ratio roughly cancel. For fixed L_acc, mdot_Edd = L_acc/L_Edd scales as M^{-1}, so Eq. 5 gives tau0 ∝ M^{2/3} mdot_Edd^{1/3} ∝ M^{1/3}. Switching to 1.15e9 M_sun changes tau0, and hence all predicted differential lags, by a factor (1.15e9 / 3e8)^{1/3} ≈ 1.56. Applying this scaling to the Year3 lags in Table 5 reduces the X=4.96 discrepancy from the quoted ~2–3 to roughly 1.2–1.6 (g, V, r, z), i.e. effectively no size problem in the traditional Wien convention, and reduces the X=2.49 discrepancy from ~4–7 to roughly 3–4. Thus the abstract's 'factor ~2–7' and the claim that 3C273 joins the size-problem sample are conditional on an adopted black-hole mass that the authors themselves cite but do not propagate. This is the most load-bearing assumption because it directly sets the denominator of the observed/predicted ratio that defines the size problem.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first dedicated accretion disc reverberation mapping campaign of the quasar 3C 273, using seven optical bands from Las Cumbres Observatory over four observing seasons. Inter-band lags are measured with two independent codes, Javelin and PyROA, and the main analysis focuses on the third season, which has the largest variability amplitude. The claimed findings are: (i) the two codes give mutually consistent lags; (ii) the lags exceed thin-disc predictions by factors of ~2-7, placing 3C 273 in the 'accretion disc size problem' class; (iii) power-law fits tau ~ lambda^beta and nu f_nu ~ nu^beta are consistent with beta = 4/3; (iv) a flat disc with a steep rim can reproduce both the lags and the variable SED; and (v) extrapolating the optical lags to near-infrared wavelengths gives 100-200 light-day radii, matching the BLR radius and supporting a dusty-disc/BLR connection. Flux variation gradient and structure function analyses are used to argue that the variability is disc-dominated.","tokens_in":29757,"tokens_out":11007,"duration_ms":121537,"significance":"If the lags are robust, this is a valuable addition to the accretion disc RM sample because 3C 273 is a high-luminosity, near-Eddington quasar, whereas the 'disc size problem' has mostly been studied in lower-luminosity AGN. The paper has clear strengths: it presents season-resolved lag tables and MCMC corner plots for both codes, it uses a contemporaneous spectrum to estimate broad-line contamination, it checks the disc-variability assumption with two independent methods, and it makes an explicit, testable prediction for a near-infrared RM campaign. However, the headline size-problem claim is sensitive to the adopted black hole mass and the chosen Wien factor convention, and the near-infrared/BLR inference is an extrapolation based on the same optical lag data. These issues are correctable but affect the central conclusions as currently stated.","major_comments":[{"comment":"The mass dependence of the predicted lags is misstated. Equation (5) gives tau0 proportional to M^(2/3) mdot^(1/3), and since mdot = L_acc/L_Edd with L_Edd proportional to M, for fixed L_acc one has mdot proportional to M^(-1) and hence tau0 proportional to M^(1/3). The claim in Section 4.1 that increasing the black hole mass to the Li et al. (2022) value of 1.15e9 M_sun makes the changes 'roughly cancel' is therefore not correct. Using M = 1.15e9 M_sun instead of 3e8 M_sun raises the predicted lags by (1.15/0.3)^(1/3) ~ 1.56. Applying this to the Year 3 lags in Table 5 lowers the X = 4.96 observed/predicted ratios from ~2-3 to roughly 1.0-1.9 (for example, the PyROA g-band ratio drops from ~3.0 to ~1.9, V from ~2.0 to ~1.3, and r from ~2.4 to ~1.5). The abstract's 'factor ~2-7' and the statement that 3C 273 joins the size-problem sample depend on the adopted mass, so the analysis should propagate the mass uncertainty and report the discrepancy under both the GRAVITY and Li et al. masses and both X conventions.","section":"Section 4.3, Table 3"},{"comment":"The decision to base the remaining analysis on Year 3 alone was made after inspecting the season-by-season results. The PyROA lags in Years 1 and 2 do not show the size excess claimed for Year 3: for example, Year 2 gives z = 6.2 days against an X = 4.96 prediction of 15.7 days, and Year 1 gives r = 5.0 days against 6.3 days. Only Year 3 shows the large lags used in Tables 5 and 8. Javelin lags are more stable across seasons, but the quoted uncertainties are internal to a single season and do not include season-to-season scatter. The paper also notes in Section 4.3 that the variability amplitude falls below the ~10% threshold recommended for RM. The authors should either report the all-season comparison in the size-problem analysis or justify the Year 3 choice with a clear, objective variability-based criterion rather than selecting the season that produces the desired lag pattern.","section":"Section 5.4, Table 9"},{"comment":"The near-infrared 'predictions' in Table 9 are not independent measurements or independent theoretical predictions. The absolute lags are computed by taking tau0 = 19.7 days, which is the average of values back-calculated from the observed optical differential lags in Table 5, and then extrapolating a fitted power law to 1.2-2.2 microns. This assumes that tau ~ lambda^beta continues unchanged beyond the observed bands. That assumption is directly challenged by the steep-rim model in Section 5.3, which predicts the rim response to become increasingly important at redder wavelengths. If the lag spectrum flattens or steepens beyond the i/z bands, the inferred 100-200 light-day dusty disc radii and the BLR connection would not follow. The BLR/FRADO conclusion should therefore be presented as a model-dependent consistency check rather than a measured result, and the abstract's causal 'therefore' should be softened accordingly.","section":"Section 5.4, Table 9"}],"minor_comments":[{"comment":"The PyROA Year 4 lags (e.g., g = 52.9 days, z = 99.0 days) are close to or at the 100-day upper prior. The text calls these 'unrealistically large' but should state explicitly that they are prior-dominated, since this is part of the justification for excluding Year 4.","section":"Section 4.3, Table 3"},{"comment":"The phrase 'Bayesian Information Criterion loss function' is imprecise; BIC is a model comparison criterion, not a loss function in the usual sense. Consider rewording to 'Bayesian Information Criterion penalty' or similar.","section":"Section 4.2.2"},{"comment":"Table 9 quotes integer values for the extrapolated absolute lags without uncertainties. Since both beta and tau0 carry uncertainties, the JHK predictions should include propagated errors or at least a sensitivity range.","section":"Section 5.4, Table 9"},{"comment":"The text reports k = 160 as the best-fit disc shape index, but the corner plot shows that log k flattens beyond log k ~ 2. The parameter k is therefore effectively a lower limit on the rim steepness, and the paper should describe it as such rather than as a tightly constrained index.","section":"Section 5.3, Figure B4"},{"comment":"The outlier-rejection threshold chi^2 > 70 is presented without justification or a sensitivity test. A brief demonstration that the fitted lags are stable under different thresholds would strengthen confidence in the cleaned light curves.","section":"Appendix A"},{"comment":"The structure function likelihood treats the N(N-1)/2 magnitude pairs as independent even though they share the same light curve. The quoted parameter uncertainties may therefore be underestimated; this should be acknowledged or tested with a bootstrap.","section":"Section 4.5, Eq. (16)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the first dedicated accretion disc RM campaign for 3C273, and it is a genuinely useful dataset: high-cadence LCO photometry over four seasons, two independent fitting codes, and a sensible effort to check variability origin via FVG and structure functions. The Year 3 lags from Javelin and PyROA agree, the power-law index is consistent with 4/3, and the detection of the disc size problem in a luminous, near-Eddington quasar would be a meaningful extension of the sample. Credit where due: the paper is honest about the low variability amplitudes, the seasonal gaps, the trend correction, and the caveats on the steep-rim model. The processed tables and the description of the outlier rejection are clear enough to follow.\n\nThe soft spots are real. The most load-bearing is the black-hole mass. The paper dismisses Li et al. (2022) with “the two changes roughly cancel,” but that is simply wrong. For fixed L_acc, mdot_Edd goes as M^{-1}, so tau0 from Eq. 5 scales as M^{2/3} mdot_Edd^{1/3} ∝ M^{1/3}. Switching from 3e8 to 1.15e9 M_sun changes every predicted lag by a factor of 1.56. Applied to the Year 3 measured lags, the X=4.96 discrepancy drops from the quoted ~2–3 to roughly 1.2–1.6 — essentially no size problem in the traditional Wien convention. The X=2.49 discrepancy shrinks to ~3–4. So the abstract’s “factor ~2–7” and the claim that 3C273 joins the size-problem sample are conditional on an adopted mass that the paper itself cites but does not propagate. This needs to be quantified, not hand-waved.\n\nThe season selection is also a concern: the final analysis rests on one of four seasons, chosen after inspecting the results. The authors give a defensible reason — Year 3 has the largest F_var and the cleanest fits — but robustness to this choice is not demonstrated. Excluding u and i from the power-law fits is reasonable given Balmer contamination, but the i-band anomaly is dismissed rather quickly. The near-IR extrapolation in Section 5.4 is not an independent measurement: it uses tau0 = 19.7 days, which is back-calculated from the observed optical lags, and it assumes the power law continues to 2 microns. No error bars are given for Table 9. If the lag spectrum flattens or steepens beyond the observed bands, the 100–200 light-day radii and the BLR connection do not follow.\n\nProcessed data are only available on request. For a result this mass-sensitive, releasing the light curves and fitting scripts should be a condition.\n\nThis paper deserves a serious referee, but it needs major revision before acceptance. The mass dependence alone could flip the central claim; the near-IR story needs to be framed as a speculative extrapolation with uncertainties. I would not cite it as a clean size-problem data point until the mass question is resolved.","headline":"First dedicated disc RM for 3C273, but the size-problem claim hinges on a black-hole mass choice the paper waves away, and the near-IR/BLR story is an extrapolation, not a measurement.","tokens_in":30363,"tokens_out":2131,"would_cite":true,"duration_ms":24687,"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":"Reverberation mapping of quasar 3C 273 shows its accretion disc is 2–7 times larger than thin-disc theory predicts.","keywords":["accretion disc","reverberation mapping","quasar 3C 273","inter-band continuum lags","accretion disc size problem","thin disc model","broad line region","dusty outflow"],"falsifier":"A dedicated near-infrared reverberation campaign measuring J-, H-, and K-band lags in 3C 273: if those lags do not fall in the 100–200 day range, or if the lag spectrum flattens or steepens beyond the optical bands, the claim that the disc extends to dust-forming radii and feeds the BLR is ruled out. A shorter-term test is to measure an absolute B-band lag independently (for example via X-ray or UV to optical cross-correlation) and check whether tau_0 = 19.7 days holds.","tokens_in":29114,"feed_emoji":"🔭","tokens_out":5642,"duration_ms":55340,"temperature":0.7,"pith_summary":"This paper reports the first dedicated accretion-disc reverberation-mapping campaign on the quasar 3C 273, using four seasons of high-cadence light curves in seven optical filters from the Las Cumbres Observatory. The inter-band time delays, recovered independently with the Javelin and PyROA fitting algorithms, agree with each other and place the emitting regions 2–7 times farther out than the Shakura–Sunyaev thin-disc model predicts, making 3C 273 one of the brightest, highest-luminosity objects yet to show the 'accretion disc size problem.' The slope of the lag–wavelength relation still matches the thin-disc exponent β = 4/3, and a 'flat disc with a steep rim' geometry can simultaneously reproduce the lags and the spectral energy distribution of the variations. Extrapolating the optical lag spectrum into the near-infrared gives radii of 100–200 light-days for putative dust-forming disc regions, consistent with the measured broad-line-region size and the rim's outer edge, suggesting the disc may extend far enough to be dusty and that the BLR could form in a dusty outflow.","feed_headline":"Quasar 3C273's accretion disc is 2-7 times too big","feed_subtitle":"Two independent reverberation-mapping fits agree, and the disc may reach dusty radii where the broad-line region forms.","key_machinery":"The central object is the lag–wavelength (reverberation) spectrum tau($\\lambda$), the light-travel delay between continuum variations at different photometric bands; it is measured with Javelin (a damped-random-walk Gaussian-process fit with a top-hat transfer function) and PyROA (a running-optimal-average fit with a delta-function transfer function), and cross-checked against thin-disc predictions tau ∝ (X $\\lambda$)^(4/3) with Wien factor X = 4.96 or 2.49. The argument then turns on two extensions of this spectrum: a 'flat disc with steep rim' model, a finite-height power-law disc H(r) = H_out (r/r_out)^k with k > 1 irradiated by a lamp-post, which fits both the lags and the variable SED, and a power-law extrapolation from the measured optical lags to the roughly 1000 K dust-forming region at J, H, and K band wavelengths.","core_discovery":"3C 273's accretion disc is a factor of about 2–7 larger than predicted by the standard geometrically thin, optically thick disc model, based on inter-band continuum lags measured with two independent reverberation-mapping codes. The lag spectrum follows tau ∝ $\\lambda$^$\\beta$ with $\\beta$ consistent with 4/3, so the disc matches the thin-disc temperature profile in shape but not in absolute size. A disc with a flat interior and a steep irradiated rim, parametrised following Starkey et al., fits both the observed lags and the variable-flux SED, placing the outer rim at roughly 120–150 light-days. Extrapolating the measured optical power law to the about 1000 K dust-sublimation region yields near-infrared lags of about 100–200 days, which match the BLR radius measured by near-infrared interferometry and the rim's outer edge; the paper therefore argues that the disc in 3C 273 may extend into dust-forming territory, so the broad-line region could emerge from a dusty disc wind.","pith_inferences":["The paper notes that a dedicated near-infrared RM campaign on 3C 273 has been completed by some of the authors with results forthcoming; if those measured JHK lags fall outside 100–200 days, the dusty-disc and BLR-connection conclusions would be directly falsified.","The absolute B-band reference lag of 19.7 days is itself back-calculated from the relative optical lags and the thin-disc normalisation; an independent absolute lag measurement, for example from X-ray or UV to optical cross-correlation, would test whether the extrapolated radii are an artefact of that zero-point.","Because the structure-function analysis finds a decorrelation timescale above roughly 300 days, monitoring longer than a decade is needed to distinguish the apparent ~3-year quasi-periodic trend from red noise; a confirmed periodicity would strengthen the disc-dominated variability interpretation.","The dusty-disc idea implies a continuous transition from accretion disc to torus, predicting that the hot-dust radius from near-infrared RM (about 400–900 light-days) is physically connected to the ~100–200 light-day disc rim, a relation testable with joint optical and near-infrared RM campaigns."],"forward_implications":["3C 273 becomes a high-luminosity, near-Eddington data point in the 'accretion disc size problem', showing that the discrepancy is not confined to low-luminosity AGN.","The consistency of Javelin and PyROA lags, despite very different assumptions about the driving variability, suggests the measured disc sizes are robust to the choice of variability model for this data set.","If the steep-rim geometry is correct, the outer disc rim at roughly 120–150 light-days should shine at about 5000 K and produce a small near-infrared excess that can be searched for in the SED.","If the dusty-disc extension is right, the inner edge of the BLR in 3C 273 should coincide with the dust-sublimation temperature, providing a testable site for dusty-outflow BLR formation models.","For a rim at radius r_out, the mean rim lag scales as (r_out/c)(1 + (2/3) sin i), giving about 120 days face-on and about 180 days at inclination 45 degrees, so future near-infrared RM can also constrain the disc inclination."],"supporting_citations":[{"why":"Supplies the black hole mass of 3e8 solar masses and the BLR radius of 145 +/- 35 light-days that the extrapolated lags are compared against.","marker":"Gravity Collaboration et al. 2018"},{"why":"Provides the thin-disc lag normalisation (Wien factor X = 2.49 and the tau_0 scaling relation) used to predict the optical lags.","marker":"Fausnaugh et al. 2016"},{"why":"Supplies the wavelength-dependent Wien factors X_Mudd for the LCO filters and the prior restrictions (decorrelation timescale and lag range) adopted in the Javelin fits.","marker":"Mudd et al. 2018"},{"why":"Provides the 'flat disc with steep rim' model that simultaneously fits the lags and the SED, including the rim-lag formula used here.","marker":"Starkey et al. 2023"},{"why":"Provides the FRADO mechanism by which the BLR emerges from dusty disc outflows at about 1000 K, which the extrapolated radii are used to test.","marker":"Czerny & Hryniewicz 2011"},{"why":"Supplies the Javelin algorithm (damped random walk with top-hat transfer function) used for one set of lag estimates.","marker":"Zu et al. 2011"},{"why":"Supplies the PyROA algorithm (running optimal average with delta-function transfer function) used for the other set of lag estimates.","marker":"Donnan et al. 2021"},{"why":"Supplies the data-reduction pipeline and the flux-flux and structure-function analysis methods applied to the LCO light curves.","marker":"Hernández Santisteban et al. 2020"},{"why":"Defines the 'accretion disc size problem' and reviews the RM method, framing the paper's motivation and placing 3C 273 in context.","marker":"Cackett et al. 2021"}],"fun_headline_variants":["3C 273 disc 2-7x bigger than thin-disc prediction","Quasar 3C273's disc lags hint at extended dusty rim","Reverberation mapping reveals oversized disc in 3C 273","3C 273 disc may extend to dusty BLR region"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dusty-disc and BLR conclusions assume that the lag–wavelength relation measured in the optical continues unchanged as the same power law out to about 2 microns, and that the back-calculated B-band reference lag of 19.7 days is the correct absolute zero-point.","fun_headline_variants_meta":{"raw":{"variants":["3C 273 disc 2-7x bigger than thin-disc prediction","Quasar 3C273's disc lags hint at extended dusty rim","Reverberation mapping reveals oversized disc in 3C 273","3C 273 disc may extend to dusty BLR region"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000354,"raw_usage":{"total_tokens":2000,"prompt_tokens":1096,"completion_tokens":904,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":712,"completion_tokens_details":{"reasoning_tokens":825}},"tokens_in":712,"tokens_out":904,"duration_ms":7916,"temperature":1.0,"reasoning_tokens":825,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T05:24:11.208561+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A dedicated near-infrared reverberation campaign measuring J-, H-, and K-band lags in 3C 273: if those lags do not fall in the 100–200 day range, or if the lag spectrum flattens or steepens beyond the optical bands, the claim that the disc extends to dust-forming radii and feeds the BLR is ruled out. A shorter-term test is to measure an absolute B-band lag independently (for example via X-ray or UV to optical cross-correlation) and check whether tau_0 = 19.7 days holds.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the wavelength-dependent Wien factors X_Mudd for the LCO filters and the prior restrictions (decorrelation timescale and lag range) adopted in the Javelin fits."},{"cited_title":"M., Bentz M","cited_arxiv_id":null,"evidence_quote":"Defines the 'accretion disc size problem' and reviews the RM method, framing the paper's motivation and placing 3C 273 in context."}],"review_version":1}