{"id":"0bdbe6b3-7b97-4053-8371-ea013ec30f62","arxiv_id":"2412.19665","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":13,"one_line_summary":"A homogeneous sample of 41 reverberation-mapped AGNs is standardizable via Hα and Hβ radius-luminosity relations, but cosmological constraints are weak.","lead":"Astronomers show that 41 active galactic nuclei with consistently measured Hα and Hβ time delays can be standardized through four radius-luminosity relations, giving distances that do not depend on the assumed cosmological model. The work supports the idea that using a homogeneous sample is the key to developing AGNs as a reliable cosmological probe.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Peculiar-velocity corrections enter every luminosity through D_L with no propagated uncertainty; a systematic error in the NED corrections would bias the R-L parameters that underpin the standardizability claim.","rationale":"The paper's central claim is modest: the 41-source homogeneous sample is standardizable in the sense that the four R-L relation parameter sets are consistent across six cosmological models, while the cosmological constraints themselves are admittedly weak. The analysis is a reasonable extension of the authors' program, and the use of a homogeneous ICCF lag set is a genuine improvement over earlier heterogeneous compilations. The model-independence test is weakened by the low redshift leverage of the sample, but the joint H(z)+BAO fits and the consistency of R-L parameters across those fits provide some additional support. The reader's weakest-assumption identification is correct and is the most load-bearing concern: peculiar-velocity corrections are applied to the very sources that anchor the low-luminosity end of the relations, and the absence of a propagated uncertainty leaves the R-L intercept and slope vulnerable to a systematic redshift error. Because the reader's CONDITIONAL verdict already reflects this issue, my stress test does not move the verdict. A secondary, non-central issue is the contradictory sentence in the Introduction ('we do find significant constraints against the standard model' after stating consistency with flat LambdaCDM), which should be corrected but does not affect the quantitative claims.","tokens_in":43793,"tokens_out":11682,"duration_ms":454411,"concrete_test":"Re-run the flat LambdaCDM and nonflat XCDM fits under two perturbations: (a) add a 300 km/s peculiar-velocity uncertainty in quadrature to sigma_tot, or draw Monte Carlo realizations of the NED correction; and (b) truncate the sample to z > 0.05 and refit. Compare gamma, beta, sigma_int, and the joint H(z)+BAO posteriors. A >1 sigma shift in any R-L parameter, or a systematic trend in beta with the z-cut, would show that the standardizability claim is not robust to peculiar-velocity systematics; no significant shift would retire the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest step is the treatment of peculiar velocities in Sec. III A. The NED Velocity Correction Calculator is applied to all 41 redshifts, and the corrected z enters D_L in Eq. (9) for every R-L and cosmological fit, but no uncertainty is assigned to the correction and none is propagated into sigma_tot in Eq. (12). For the roughly 13 sources with z < 0.05, a 300 km/s peculiar-velocity error changes log D_L by about 0.04-0.1 dex, comparable to or larger than the flux errors and a significant fraction of the fitted intrinsic scatter (~0.25 dex). These low-z sources populate the low-luminosity end of each R-L relation and strongly set the intercept and slope (e.g., flat LambdaCDM gamma = 0.571 +/- 0.055, beta = 2.336 in Eq. 15). A coherent error in the NED corrections, such as an imperfect local-flow model, would bias D_L for these anchor points and thereby shift beta and gamma and bias distance moduli, without necessarily altering the apparent model-independence used to claim standardizability.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper tests whether 41 low-redshift reverberation-mapped AGNs can be used as standardizable distance indicators. It simultaneously fits the parameters of four radius-luminosity relations (Hα and Hβ time delays, combined with either broad Hα luminosity or monochromatic 5100 Å luminosity) and the parameters of six cosmological models, using a Gaussian likelihood with an intrinsic scatter term. The authors report that the R-L parameters are stable across all six models, that the slopes are steeper than the simple photoionization value of 0.5 and steeper than those of earlier larger Hβ samples, and that the cosmological constraints from the 41 AGNs are weak but consistent with H(z)+BAO data. They interpret the sample homogeneity, particularly the consistent ICCF lag determinations, as the reason the results differ from their earlier 118-object Hβ analysis.","tokens_in":44100,"tokens_out":5517,"duration_ms":54668,"significance":"If the central standardizability claim holds, this is a useful step toward developing low-redshift RM AGNs as a cosmological probe, complementing higher-redshift Mg II and C IV samples. The paper is methodologically transparent: the likelihood, priors, MCMC setup, and model-comparison criteria are clearly specified, and the simultaneous fitting of R-L and cosmological parameters is the appropriate route around the circularity of assuming a cosmology when computing luminosities. The comparison with the independent 118-object sample is a genuine strength. However, the evidence for standardizability is currently limited by two issues: the low-redshift peculiar-velocity corrections enter every luminosity distance without propagated uncertainty, and the data have so little cosmological constraining power that model-independence of the R-L parameters is a weak test. The physical interpretation of the steeper slopes also relies on Eddington-ratio effects that are only marginally detected in the sample itself.","major_comments":[{"comment":"The NED Velocity Correction Calculator is applied to all 41 redshifts, and the corrected redshift enters the luminosity distance in Eq. (9) and hence every R-L and cosmological fit, but no uncertainty is assigned to the peculiar-velocity correction and none is propagated into sigma_tot in Eq. (12). For the eleven sources with z < 0.05 (e.g., NGC 4151, NGC 6814, Mrk 142, Arp 151, NGC 5548), a peculiar-velocity error of 300 km/s changes log D_L by roughly 0.02–0.08 dex, which is comparable to or larger than the flux errors and a substantial fraction of the fitted intrinsic scatter (sigma_int ~ 0.25 dex). These low-z, low-luminosity sources anchor the low-luminosity end of each R-L relation and therefore strongly influence the intercept and slope. A coherent error in the NED model, such as an imperfect local-flow correction, would bias beta and gamma and thus bias the distance moduli without necessarily changing the apparent model-independence of the R-L parameters. The authors should propagate a peculiar-velocity uncertainty through Eq. (12), or at minimum perform a sensitivity test in which all low-z redshifts are shifted by a plausible systematic velocity and show that the values in Tables IV and VI are stable within their quoted uncertainties.","section":"Sec. III A, Table I, Eq. (12)"},{"comment":"The main evidence for standardizability is that the R-L parameters vary by at most 0.22 sigma in gamma and 0.50 sigma in beta across the six cosmological models (Table VII). This is a weak test because the 41 AGNs by themselves barely constrain cosmology: in flat LambdaCDM the H-alpha broad data give only a 1-sigma lower limit Omega_m0 > 0.424, and the 2-sigma limits are set by the priors (Table IV). When the data have essentially no cosmological constraining power, the fitted R-L parameters are expected to be nearly model-independent whether or not the sources are truly standardizable. To make the standardizability claim load-bearing, the authors should calibrate the expected spread in beta and gamma under a null hypothesis, for example by running mock catalogs of non-standardizable sources drawn from a single cosmology, or by showing that the constraints tighten in a meaningful way when H(z)+BAO data are added. Absent such a test, the statement that the sources are 'standardizable' should be softened to consistency with standardizability.","section":"Sec. V, Tables IV and VII"},{"comment":"The interpretation that the steeper R-L slopes compared with the 118-object H-beta sample arise from the absence of high-accreting sources is not strongly supported by the internal diagnostic presented in Fig. 6. The anti-correlation between Delta R and Eddington ratio is weak and statistically marginal for the H-alpha relations (Spearman rho = -0.10, p = 0.55 for R_Halpha-L_Halpha; rho = -0.18, p = 0.26 for R_Halpha-L_5100), and only the H-beta relations show more suggestive correlations (p = 0.073 and p = 0.011). As the authors note, lambda_Edd scales as L/R_BLR while Delta R scales approximately as R_BLR/L^0.6, so the correlation may be partly a self-correlation. Since the slope difference versus earlier samples is one of the main physical conclusions, the manuscript should either fit an explicit Eddington-ratio term in the R-L relation to test the proposed interpretation, or clearly state that the self-correlation caveat prevents a definitive physical conclusion.","section":"Sec. VI A, Fig. 6"}],"minor_comments":[{"comment":"The caption of Fig. 4 and the corresponding discussion refer to 'Tables IV and IV'; the second table should be Table VI.","section":"Sec. V, Fig. 4 caption and text"},{"comment":"There is a typo in the sentence describing the comparison with Cho et al.: 'wehavecorected' should read 'we have corrected'.","section":"Sec. VI A"},{"comment":"The asymmetric treatment of the time-delay uncertainties is described in prose but not in the equation; defining sigma_log_tau explicitly as a piecewise function of the residual sign would make the likelihood unambiguous.","section":"Eq. (12)"},{"comment":"In the luminosity-luminosity fits, the luminosity distance is computed using posterior mean cosmological parameters, so cosmological-parameter uncertainty is not propagated in that step; this should be stated explicitly as a limitation.","section":"Sec. IV, Eqs. (13)-(14)"}],"recommendation":"major_revision","confidential_remarks":"I do not see grounds for rejection; the peculiar-velocity and weak-test issues are fixable with additional sensitivity analyses and a more careful interpretation. The paper would be strengthened by a null-distribution simulation for the model-independence test and by explicitly acknowledging that the standardizability claim is currently based on a low-power dataset. The physical interpretation of the slope difference should also be framed more cautiously given the marginal Eddington-ratio correlations reported in Fig. 6."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, read this if you work on AGN cosmology. The useful bit is the sample construction and the negative result framing: a homogeneous 41-object Hα/Hβ reverberation sample gives R–L parameters stable across six cosmological models, and the earlier ≥2σ disagreement with H(z)+BAO from the 118-object sample goes away. That is a legitimate new application of the simultaneous R–L/cosmology fitting method, not just a rerun: they correct redshifts for peculiar velocities, propagate asymmetric lag errors, and let cosmology float instead of fixing it. The paper is honest about the weak cosmological constraints—the 2σ limits are prior-dominated—so the standardizability claim rests mainly on the model-independence of the R–L parameters, which is a real but modest test.\n\nThe soft spots are proportionate. The peculiar-velocity corrections are the weakest link: the NED calculator is applied to every object, corrected z enters D_L in Eq. (9), and no uncertainty on the correction is propagated into σ_tot. For z<0.05 objects, a 300 km/s error gives a 0.04–0.1 dex shift in log D_L, comparable to the intrinsic scatter these low-luminosity anchors help set. This does not invalidate the paper, but it needs a robustness check—rerun with a few hundred km/s peculiar-velocity dispersion added, or drop the z<0.02 objects. If the results survive, the claim is considerably stronger. Also, the text in Sec. I appears to say \"we do find significant constraints against the standard model\" while the abstract and results say the opposite; that sentence needs fixing.\n\nOne more caution: the steeper-slope explanation (sample lacking high-accreting sources) is plausible but not established; the Eddington-ratio correlation is weak and self-correlation-prone, as the authors themselves note. Fine as a suggestion, not a result.\n\nOverall, the methods are clear, the likelihood and priors are specified, and the comparison with earlier work is fair. I’d cite it, and I’d send it to a referee—it deserves a serious look. The main request should be a sensitivity analysis on peculiar velocities. The paper's own flag that a larger homogeneous sample is needed is the right summary.","headline":"A careful, honest extension of the RM-AGN cosmological program: the homogeneous 41-object sample is a real step, but weak cosmology and unpropagated peculiar-velocity systematics keep it from being a decisive validation.","tokens_in":44653,"tokens_out":2231,"would_cite":true,"duration_ms":24569,"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":"A homogeneous sample of 41 low-redshift AGNs becomes standardizable through four radius–luminosity relations.","keywords":["reverberation mapping","active galactic nuclei","radius-luminosity relation","broad-line region time delays","H-alpha emission","H-beta emission","cosmological distance indicators","AGN standardizability"],"falsifier":"Recompute the four radius–luminosity fits after removing all sources with $z < 0.05$, or after replacing the adopted peculiar-velocity corrections with an independent local-flow velocity field; if the slopes or intercepts shift by more than the quoted uncertainties, or the cosmological consistency with H(z)+BAO data disappears, the standardizability claim would be falsified.","tokens_in":43594,"feed_emoji":"🔭","tokens_out":7194,"duration_ms":64014,"temperature":0.7,"pith_summary":"This paper argues that 41 low-redshift active galactic nuclei with reverberation-mapped Hα and Hβ time delays can serve as standardizable distance indicators. The basis is four radius–luminosity relations, pairing each line delay with either the broad Hα luminosity or the monochromatic 5100 Å luminosity, whose fitted slope and intercept stay essentially unchanged across six different cosmological models. That stability is the paper's criterion for standardizability, because it shows the calibration does not secretly depend on the assumed expansion history. The cosmological constraints from these 41 sources are weak but consistent within 2σ with those from H(z)+BAO data, in contrast to an earlier, larger but less homogeneous Hβ sample that disagreed by more than 2σ. If the claim holds, consistent time-lag measurement is the key to turning reverberation-mapped AGNs into a usable cosmological probe.","feed_headline":"A uniform method makes 41 AGNs standardizable","feed_subtitle":"Their radius–luminosity relations hold across six cosmologies, clearing a path to a new distance probe.","key_machinery":"The load-bearing object is the radius–luminosity relation, $\\log(\\tau/\\mathrm{day}) = \\beta + \\gamma \\log(L/10^{44}\\,\\mathrm{erg\\,s^{-1}})$, where $\\tau$ is the measured rest-frame line delay and $L$ is the luminosity inferred from the observed flux through $L = 4\\pi D_L^2 F$. Because the delay is measured and the flux is observed, a calibrated relation turns each AGN into a distance estimate whose only free ingredients are $\\beta$, $\\gamma$, and the assumed cosmology. The paper fits $\\beta$, $\\gamma$, an intrinsic scatter parameter, and cosmological parameters simultaneously in a joint likelihood; if the relation parameters are insensitive to the cosmological model, the dataset is declared standardizable. This machinery is applied to four combinations of Hα or Hβ delay with broad Hα or monochromatic 5100 Å luminosity, using measured fluxes directly rather than cosmology-dependent luminosities.","core_discovery":"The authors claim that, for a sample of 41 Type 1 AGNs whose Hα and Hβ lags were all determined with the interpolated cross-correlation function, the four radius–luminosity relations of the form $\\log(\\tau/\\mathrm{day}) = \\beta + \\gamma \\log(L/10^{44}\\,\\mathrm{erg\\,s^{-1}})$ have slopes $\\gamma$ near 0.54–0.62 and intercepts larger for Hα than for Hβ. These parameters do not shift by more than about half a $\\sigma$ when the assumed cosmology changes among flat and nonflat ΛCDM, XCDM, and ϕCDM models, which the authors take as evidence that the sources are standardizable. The measured slopes are steeper than the 0.5 predicted by a simple photoionization model and steeper than slopes from larger Hβ samples, a difference the authors attribute to the absence of high-accreting sources in their low-redshift sample. The cosmological parameters derived from these AGNs are only weakly constrained but are consistent within 2σ with those from more established probes, unlike the earlier 118-source Hβ analysis, whose constraints were in greater tension.","pith_inferences":["If peculiar-velocity systematics can be controlled, the calibrated radius–luminosity intercepts could provide absolute distance measurements at $z \\lesssim 0.5$, offering a low-redshift anchor independent of supernovae.","Combining this low-redshift Hα/Hβ calibration with higher-redshift Mg II and C IV reverberation samples could cover a much wider redshift range with one consistent method; the paper notes the current sample is too small for such a joint fit now.","The observed negative correlation between radius–luminosity residuals and Eddington ratio suggests a testable extension: explicitly including accretion rate as a second calibration parameter might remove the steep-slope offset and identify accretion physics, not sample selection, as the driver.","The mild anti-correlation between Hα equivalent width and continuum luminosity implies that line-based and continuum-based luminosities are not perfectly interchangeable, so future calibrations may need separate slope parameters for the two luminosity definitions."],"forward_implications":["Larger homogeneous reverberation-mapping samples would yield radius–luminosity calibrations independent of the assumed cosmology, allowing low-redshift AGNs to be appended to the distance ladder.","The radius–luminosity parameters change by at most about 0.8σ when H(z)+BAO data are added, so combining these AGNs with established probes is safe and shifts cosmological constraints by only about 0.1σ.","Because Hβ lags paired with Hα luminosity show the lowest intrinsic scatter, future monitoring campaigns should prioritize homogeneous Hβ lag measurement to tighten the calibration.","The steeper-than-0.5 slopes imply that the simple photoionization scaling $R \\propto L^{0.5}$ is incomplete, and samples dominated by low-Eddington-ratio sources will systematically overestimate the slope.","The earlier disagreement with standard cosmological probes is attributed to sample inhomogeneity rather than to the radius–luminosity method itself, so consistent lag determination is the enabling condition for the probe."],"supporting_citations":[{"why":"Supplies the homogeneous 41-source Hα and Hβ lag sample with interpolated-cross-correlation-function time delays and luminosity measurements.","marker":"[9]"},{"why":"Provides the earlier larger 118-source Hβ sample whose cosmological constraints disagreed with established probes, the baseline this paper contrasts with its homogeneous sample.","marker":"[8]"},{"why":"Defines the interpolated cross-correlation function technique whose uniform application is the stated source of sample homogeneity.","marker":"[87, 88]"},{"why":"Gives the likelihood function with an intrinsic scatter term used to fit the radius–luminosity relations.","marker":"[90]"},{"why":"Supplies the H(z)+BAO dataset and likelihood used as the better-established cosmological probe for joint fits and consistency checks.","marker":"[89]"},{"why":"Provides a larger 240-object Hβ-based radius–luminosity comparison with a shallower slope that this paper attributes to accretion-rate selection.","marker":"[80]"},{"why":"States the simple photoionization prediction of slope 0.5 that the measured steep slopes are compared against.","marker":"[12, 93]"},{"why":"Documents the shortening of broad-line-region time delays in high-accreting sources, which the paper uses to explain its steeper slope.","marker":"[15–17]"}],"fun_headline_variants":["41 AGNs standardized with Hα and Hβ lags","AGN R-L relations stable across six cosmologies","Uniform lag method yields standardizable AGNs","Homogeneous AGN sample eases cosmological tension","Steeper R-L slopes still standardize 41 AGNs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The peculiar-velocity corrections applied to the low-redshift sources are assumed accurate enough that the resulting luminosity distances are unbiased, and the correction's own uncertainty is not propagated into the analysis.","fun_headline_variants_meta":{"raw":{"variants":["41 AGNs standardized with Hα and Hβ lags","AGN R-L relations stable across six cosmologies","Uniform lag method yields standardizable AGNs","Homogeneous AGN sample eases cosmological tension","Steeper R-L slopes still standardize 41 AGNs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001173,"raw_usage":{"total_tokens":4920,"prompt_tokens":1086,"completion_tokens":3834,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":702,"completion_tokens_details":{"reasoning_tokens":3757}},"tokens_in":702,"tokens_out":3834,"duration_ms":25909,"temperature":1.0,"reasoning_tokens":3757,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:00:38.763098+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the four radius–luminosity fits after removing all sources with $z < 0.05$, or after replacing the adopted peculiar-velocity corrections with an independent local-flow velocity field; if the slopes or intercepts shift by more than the quoted uncertainties, or the cosmological consistency with H(z)+BAO data disappears, the standardizability claim would be falsified.","supporting_citations":[{"cited_title":"Cao and B","cited_arxiv_id":null,"evidence_quote":"Supplies the H(z)+BAO dataset and likelihood used as the better-established cosmological probe for joint fits and consistency checks."},{"cited_title":"Woo, et al., Astrophys","cited_arxiv_id":null,"evidence_quote":"Provides a larger 240-object Hβ-based radius–luminosity comparison with a shallower slope that this paper attributes to accretion-rate selection."}],"review_version":1}