{"id":"91c4df63-bd1c-4b1f-96f2-4722c045ed86","arxiv_id":"2509.07848","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Gaussian process reconstruction of the cosmic distance duality parameter eta(z) from BAO, galaxy clusters, supernovae, and quasars finds consistency with eta=1 at the 2-sigma level out to z about 2.33.","lead":"This paper combines four kinds of cosmic distance measurements to check whether the standard distance duality relation still holds at high redshifts. The authors find no significant violation within current uncertainties.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quasar D_L values are likely model-dependent and the 1σ ΛCDM filter makes the high-z CDDR test circular, so the null result is not an independent validation.","rationale":"The reader identified the 1σ ΛCDM filtering as the weakest assumption, which is correct and explicitly self-admitted in the Conclusions. My concern extends one level deeper: the quasar distance moduli themselves are likely not cosmology-independent, so the filter is only the visible part of a more fundamental circularity. The low-redshift combinations (η_l and η_d) and the standard GP pipeline provide some independent support, and the paper is not unsalvageable, so a conditional acceptance is appropriate. The reader's CONDITIONAL verdict remains the right call; the conditional should require the authors to demonstrate (or correct for) the model dependence of the QSO D_L, or to restrict the abstract's high-redshift claim accordingly. My agreement is partial because the reader focused on the filter as the main issue, whereas I see the calibration of the quasar distances as the underlying load-bearing problem.","tokens_in":15595,"tokens_out":9944,"duration_ms":96744,"concrete_test":"Inspect the Lusso et al. (2020) data release to determine whether the quasar distance moduli were derived from a joint fit assuming a cosmological model such as flat ΛCDM. If they were, redo the high-z CDDR analysis using the raw X-ray and UV fluxes, fitting the relation parameters (α, β) simultaneously with a free CDDR-parameter function η(z), and do not apply the 1σ ΛCDM filter. If the resulting η(z) at z≳1.5 remains consistent with unity within 2σ, the filter is not the cause of the null result; if it deviates, the paper's central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's high-redshift claim rests on quasar luminosity distances from the Lusso et al. (2020) catalog. The manuscript does not state whether these distance moduli are model-independent, and the standard Risaliti-Lusso method calibrates the X-ray/UV relation by fitting the Hubble diagram under an assumed cosmology, so the quoted D_L values typically inherit that model dependence. Independently, Section III A explicitly filters quasars by requiring their distance moduli to lie within the 1σ dispersion of SNIa residuals around flat ΛCDM (upper panel, Fig. 2). Because flat ΛCDM satisfies the CDDR, this selection removes the high-redshift objects that would carry a CDDR-violation signal. After filtering, the binned D_L at z > 1.5 are dominated by the surviving quasars, so the GP reconstruction of η(z) is built from data preselected to agree with ΛCDM. The conclusion that η(z) is consistent with unity at 2σ at high z is therefore a consequence of the selection, not an independent test. The authors acknowledge the filter's model dependence in the Conclusions, but the abstract presents the result as supporting the CDDR, which overstates the analysis's power. The BAO D_A calibration with r_d = 147.05 Mpc from Planck ΛCDM adds a secondary but real model dependence on the angular-distance side.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper tests the cosmic distance duality relation (CDDR), η(z) = D_L(z)/[(1+z)^2 D_A(z)], by combining angular diameter distances from transverse BAO and galaxy clusters with luminosity distances from Pantheon+ SNe Ia and a quasar catalog. Quasar distances are filtered by retaining only objects whose residuals relative to flat ΛCDM lie within the 1σ dispersion of the SNIa residuals, then binned; Gaussian Process regression is used to reconstruct η(z). Three combinations are considered: GC+SNIa+QSO at low z, BAO+SNIa+QSO at high z, and an angular-only GC+BAO combination. All reconstructions are reported to be consistent with η=1 at 2σ, and the authors conclude that the CDDR is supported up to z≈2.33.","tokens_in":15889,"tokens_out":2641,"duration_ms":26238,"significance":"If the high-redshift test were independent, extending CDDR verification to z≈2.3 with quasars would be a valuable complement to SNIa-only analyses. The paper uses standard Gaussian Process machinery, provides a clear data summary in Table I, and explicitly acknowledges in the Conclusions that the quasar filter assumes ΛCDM. These are strengths. However, the central high-redshift claim is weakened by the model-dependent pre-selection of the quasar sample and by the unstated cosmology dependence of the quasar distance moduli. The paper is therefore more a consistency check within ΛCDM than an independent test of the CDDR, and the abstract overstates the evidential weight of the result.","major_comments":[{"comment":"The quasar filtering procedure selects only quasars whose distance-modulus residuals lie within the 1σ dispersion of SNIa residuals around flat ΛCDM. Because flat ΛCDM satisfies the CDDR by construction, this filter preferentially removes quasars that would carry a CDDR-violation signal, and the surviving high-z D_L sample is preselected to agree with η=1. Consequently, the high-redshift η_h reconstruction in Fig. 4 cannot be presented as an independent confirmation of the CDDR; it is partly a consequence of the selection. The Conclusions admit the model dependence, but the Abstract's claim that the results 'support its validity even at high redshifts' is not justified by the analysis as presented. A concrete remedy would be to repeat the analysis without the ΛCDM-based filter (e.g., using robust binning or a heavier-tailed likelihood) and to show the redshift-dependent retention fraction, or to inject a simulated η≠1 and demonstrate that the pipeline can recover it.","section":"III A and Fig. 2"},{"comment":"The manuscript does not state whether the quasar distance moduli taken from the Lusso et al. catalog are cosmology-independent. The standard Risaliti-Lusso method calibrates the L_X–L_UV relation by fitting the Hubble diagram under an assumed cosmological model, so the quoted D_L values generally inherit that model dependence. If the catalog distances are already ΛCDM-calibrated, the subsequent 1σ ΛCDM filter compounds the circularity. The authors should either demonstrate that the adopted D_L values are insensitive to the calibration cosmology or explicitly restrict their conclusions to 'consistency with ΛCDM-based quasar distances.'","section":"II D"},{"comment":"The angular-only combination η_d, constructed from GC and BAO D_A measurements alone, does not involve D_L and therefore cannot test the CDDR as defined in Eq. (15). It tests only the mutual consistency of two angular-diameter-distance datasets. The text and figure caption should state this limitation explicitly; otherwise the reader may infer that η_d provides evidence about the duality relation, which it does not.","section":"IV, Fig. 4 (η_d)"}],"minor_comments":[{"comment":"The binning uncertainty in Eq. (11) is the standard deviation of the binned values, not the standard error of the mean, and individual measurement errors are not propagated into the binned points. The authors should justify this choice or use inverse-variance weighting, since it directly affects the width of the GP confidence intervals.","section":"III A, Eq. (11)"},{"comment":"The caption says the red dashed lines represent the '1σ dispersion of the SNIa distribution,' but the plotted lines appear to be constant offsets in μ. Please clarify whether the dispersion is computed globally or as a function of redshift, and state the numerical value used.","section":"Fig. 2 caption"},{"comment":"Table I lists 2195 quasars, while the filtered sample contains 1160 quasars. Adding the filtered count and the retention fraction to the table or the text would improve transparency.","section":"II D and Table I"},{"comment":"There are minor typographical and grammatical issues, including 'FlAT-ΛCDM' and 'e Ωk' in Section IV, and the phrase 'Enabling access to the robustness of our analyses regarding dependence on priors' in the Conclusions is unclear and should be rewritten.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The core issue is not that the paper is sloppy; the authors explicitly acknowledge the ΛCDM dependence of their filter. The problem is that the abstract and the phrase 'supporting its validity even at high redshifts' go beyond what the filtered-QSO analysis can establish. If the authors reframe the high-z result as a consistency test within ΛCDM and add robustness checks without the filter, the paper could be publishable. I would also encourage the editor to ensure the quasar D_L calibration issue is fully addressed in the revision, since it is central to the claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know upfront: the abstract overstates the high-redshift result. The quasar sample is filtered by requiring each quasar's distance modulus to lie within the 1-sigma scatter of SNe around flat LambdaCDM. Since flat LambdaCDM satisfies the CDDR by construction, the high-z eta(z) reconstruction is largely built from data that already agrees with the relation being tested. That is a real circularity, and the paper admits it in the conclusions, though the abstract does not. The angular-only eta_d combination (GC vs BAO) avoids that problem, but it does not test the distance duality in any meaningful sense; it is just a cross-check of two angular-diameter distance catalogs.\n\nWhat is actually new: combining Pantheon+ SNe, Lusso et al. quasars, DESI DR2 transverse BAO, and De Filippis clusters into one GP reconstruction, with a binning pipeline, is a new data combination. The GP machinery is standard and competently executed. The low-redshift eta_l (clusters vs SNe+QSO over 0.023<z<0.784) is a reasonable consistency test, and the BAO-based eta_h at intermediate z has some value independent of the quasar filter if you treat it as a systematic check rather than a discovery-level test. The authors also deserve credit for explicitly flagging the LambdaCDM dependence of the filtering step, something many such papers gloss over.\n\nThe soft spots are proportionate to the claim. The central high-z null result is not an independent validation of the CDDR; removing the filter would be the obvious robustness check, and varying the 1-sigma threshold would quantify how much the conclusion is driven by selection. The Lusso QSO distance moduli are also not stated to be model-independent; the usual Risaliti-Lusso calibration fits the X-ray/UV relation against an assumed cosmology, so the D_L values likely inherit that dependence. A short paragraph addressing this would be enough. Finally, the data availability line says \"on request,\" which is weak for a paper whose pipeline is otherwise reproducible in principle.\n\nBottom line: this is a solid consistency check, not a decisive high-z test. A serious referee should see it; with revisions, including an unfiltered quasar analysis, a threshold variation, and a clearer statement of what is actually being tested, it could be a useful paper. I would not cite it as evidence for the CDDR at high z in its current form, but I would bring it to reading group and I would referee it myself.","headline":"The high-redshift CDDR claim is weakened by a LambdaCDM-based quasar cut, but the paper is a competent GP consistency check worth engaging with.","tokens_in":707,"tokens_out":881,"would_cite":false,"duration_ms":22094,"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":"The cosmic distance duality relation remains consistent with η = 1 within 2σ across 0.023 ≤ z ≤ 2.33 when tested with Gaussian Process reconstruction.","keywords":["cosmic distance duality relation","Gaussian process regression","non-parametric reconstruction","angular diameter distance","luminosity distance","quasars","type Ia supernovae","baryon acoustic oscillations"],"falsifier":"Re-run the $\\eta(z)$ reconstruction without the $\\Lambda$CDM-based quasar filter, or with a model-independent robust cut, and compare the binned quasar-only $\\eta$ values at $z > 1$ with the paper's filtered result; a mean deviation above $2\\sigma$ in the unfiltered sample would contradict the claim that the CDDR holds at high redshifts.","tokens_in":15384,"feed_emoji":"📏","tokens_out":8971,"duration_ms":72782,"temperature":0.7,"pith_summary":"This paper tries to establish whether the cosmic distance duality relation—the identity $D_L(z) = (1+z)^2 D_A(z)$ linking luminosity distance to angular diameter distance—continues to hold at redshifts beyond the supernova range. It combines angular diameter distances from transverse baryon acoustic oscillations and galaxy clusters with luminosity distances from Type Ia supernovae and quasars, then reconstructs the ratio $\\eta(z) = D_L/[(1+z)^2 D_A]$ with Gaussian Process regression. The central result is that $\\eta(z)$ stays within $2\\sigma$ of unity across $0.023 \\le z \\le 2.33$ in every dataset pairing tested. If correct, this supports the three assumptions behind the relation—metric gravity, null geodesics, and photon-number conservation—at high redshifts where quasars become the only standard candles.","feed_headline":"Distance-duality relation holds to z=2.33 at 2σ","feed_subtitle":"BAO, cluster, supernova, and quasar distances agree with η = 1 out to z ≈ 2.33","key_machinery":"The central object is the duality parameter $\\eta(z) = D_L(z)/[(1+z)^2 D_A(z)]$, with $\\eta = 1$ marking exact CDDR validity. The argument is carried by a Gaussian Process regression with a squared-exponential kernel, which reconstructs $D_A(z)$ and $D_L(z)$ from sparse, unevenly sampled measurements and propagates $1\\sigma$ and $2\\sigma$ bands into $\\eta(z)$. For galaxy clusters, the paper uses the relation $D_A^{\\rm cluster}(z) = \\eta^2(z) D_A(z)$ to rewrite the test as $\\eta(z) = D_A^{\\rm cluster}(z)(1+z)^2 / D_L(z)$, avoiding circular use of the relation itself. Quasar luminosity distances are first filtered against the $1\\sigma$ dispersion of supernova residuals around flat $\\Lambda$CDM and then binned into 25 points ($\\Delta z = 0.1$) to suppress their large intrinsic scatter.","core_discovery":"On the paper's own terms, the discovery is that no statistically significant deviation from the cosmic distance duality relation appears when non-parametric Gaussian Process reconstruction is applied to combined BAO/galaxy-cluster angular diameter distances and supernova/quasar luminosity distances. Three reconstructions are performed: a low-redshift pairing of galaxy clusters with supernovae plus quasars, a high-redshift pairing of BAO with supernovae plus quasars, and a pure angular-distance cross-check of clusters against BAO. In all three, $\\eta(z) = D_L(z)/[(1+z)^2 D_A(z)]$ is consistent with 1 at the $2\\sigma$ level over the full interval $0.023 \\le z \\le 2.33$. The paper therefore concludes that the CDDR remains valid under current observations even at high redshifts, while stressing that the quasar selection step explicitly assumes flat $\\Lambda$CDM.","pith_inferences":["The $\\Lambda$CDM-based quasar filter may bias the test toward $\\eta = 1$: if the CDDR is violated at high redshift, the excluded quasars would preferentially carry the deviation. Repeating the analysis with a model-independent outlier cut would show whether the consistency survives.","Gaussian Process smoothing has a characteristic correlation length; a narrow or oscillatory deviation in $\\eta(z)$ shorter than that length would be smeared into the $2\\sigma$ band. A redshift-bin-by-bin $\\chi^2$ test against unity would probe that regime.","Treating the quasar filter threshold as a free parameter and marginalizing over it would convert the present conditional test into a genuinely model-independent statement.","If future high-redshift samples reduce quasar distance errors below the current intrinsic scatter, the same pipeline could detect a redshift-dependent $\\eta(z)$ at the few-percent level."],"forward_implications":["Exotic photon-number-violating physics, such as axion-like photon mixing, is not needed to explain the distance measurements out to $z \\approx 2.33$.","Filtered and binned quasar distances can serve as high-redshift standard candles, extending cosmography beyond the Type Ia supernova horizon.","BAO-derived angular diameter distances provide the tightest high-redshift constraint, while galaxy clusters act as a low-redshift consistency check.","Improved quasar luminosity-distance precision would sharpen the same test from a $2\\sigma$ consistency into a few-percent constraint on $\\eta(z)$."],"supporting_citations":[{"why":"Defines the cosmic distance duality relation $D_L = (1+z)^2 D_A$ that the whole test targets.","marker":"[6]"},{"why":"Supplies the Type Ia supernova luminosity distances used as the low-redshift standard-candle anchor.","marker":"[66]"},{"why":"Supplies the 2421-quasar catalog whose X-ray/UV luminosity relation extends $D_L$ measurements to high redshift.","marker":"[69]"},{"why":"Provides the 25 galaxy-cluster angular diameter distances from SZE/X-ray observations under an oblate-spheroid model.","marker":"[63]"},{"why":"Provides 14 transverse BAO angular-scale measurements that anchor $D_A$ at intermediate redshifts.","marker":"[50]"},{"why":"Provides six $D_M/r_d$ measurements that extend the BAO-based $D_A$ reconstruction to $z=2.33$.","marker":"[48]"},{"why":"Establishes the Gaussian Process regression method used to reconstruct $D_A(z)$, $D_L(z)$, and $\\eta(z)$ non-parametrically.","marker":"[72]"},{"why":"Derives the general SZE/X-ray cluster distance expression that motivates writing $D_A^{\\rm cluster} = \\eta^2 D_A$ in the CDDR test.","marker":"[59]"}],"fun_headline_variants":["Distance-duality law holds to z=2.33 in GP reconstruction","No significant CDDR violation up to z=2.33","High-redshift distances obey cosmic duality within 2σ","GP analysis: distance-duality intact from 0.023 to 2.33"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes flat $\\Lambda$CDM is correct when it rejects quasars whose distance-modulus residuals fall outside the $1\\sigma$ supernova band; if the CDDR is actually violated at high redshift, that filter removes the very objects carrying the signal and biases $\\eta(z)$ toward 1.","fun_headline_variants_meta":{"raw":{"variants":["Distance-duality law holds to z=2.33 in GP reconstruction","No significant CDDR violation up to z=2.33","High-redshift distances obey cosmic duality within 2σ","GP analysis: distance-duality intact from 0.023 to 2.33"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000698,"raw_usage":{"total_tokens":3121,"prompt_tokens":879,"completion_tokens":2242,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":2165}},"tokens_in":495,"tokens_out":2242,"duration_ms":16523,"temperature":1.0,"reasoning_tokens":2165,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:08:59.398301+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the $\\eta(z)$ reconstruction without the $\\Lambda$CDM-based quasar filter, or with a model-independent robust cut, and compare the binned quasar-only $\\eta$ values at $z > 1$ with the paper's filtered result; a mean deviation above $2\\sigma$ in the unfiltered sample would contradict the claim that the CDDR holds at high redshifts.","supporting_citations":[{"cited_title":"Baryon acoustic oscillations in thin redshift shells from BOSS DR12 and eBOSS DR16 galaxies","cited_arxiv_id":"2112.10000","evidence_quote":"Provides 14 transverse BAO angular-scale measurements that anchor $D_A$ at intermediate redshifts."}],"review_version":2}