{"id":"0ff0bb93-3117-46ec-9ad0-dac146ac3884","arxiv_id":"2608.07654","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Whether the data prefer dynamical dark energy depends on the assumed early-universe values of the sound horizon and matter density; early solutions to the Hubble tension move to a nexus where the preference drops below 0.5 sigma.","lead":"This analysis of galaxy and supernova data shows that the claimed evidence for dark energy changing over time depends on assumptions about the early universe. Models that solve the Hubble tension shift two key parameters into a region where that evidence all but disappears, connecting two big cosmological debates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sigma-reduction map is built on a Gaussian-Mahalanobis approximation that is never validated; if the CPL posteriors are non-Gaussian or unrepresentative, the nexus and the claimed evidence drop could shift materially.","rationale":"I read the paper as a conditional-evidence argument: for each possible (Ω_m, h_rd) the BAO+SN data yield some preference for dynamical dark energy, and early-universe Hubble-tension models sit where that preference is weak. The logic is sound and the CMB-independent confirmation that the nexus lies near the early-model region is interesting. The weak point is the map itself: it is produced by an analytic Gaussian approximation rather than by full likelihood evaluation, and the authors explicitly rely on 'near-Gaussian' contours without showing convergence or diagnostics. If the approximation is inaccurate, the quoted drop from >2σ to 1.3σ and the exact nexus position may be off, though the qualitative conclusion may survive. I further note that the Chebyshev reconstruction is a single-point check, not a full map, so it cannot rescue the grid-wide approximation. The paper is worthy of conditional acceptance; the authors should either validate Gaussianity at representative points or replace the approximation with full posterior computations. The 'necessarily' in the conclusions is an overreach given only three early-model classes are plotted, but that is secondary to the map validation issue.","tokens_in":7665,"tokens_out":14534,"duration_ms":140096,"concrete_test":"At three representative points in Fig. 1 — the nexus center, the ΛCDM CMB+BAO+SN best fit, and one point on the 2σ BAO contour — re-run the BAO+SN likelihood with MCMC over (w0, wa) for fixed (Ω_m, h_rd), using the same data and priors. Compute the actual posterior density at (-1,0) (e.g., by Savage-Dickey density ratio or by the fraction of samples within a small radius) and convert to a significance; compare with the Mahalanobis value from Eqs. (A1)-(A2). If any point differs by more than ~0.3σ, the Gaussian approximation fails and the green contours must be recomputed. As a second check, re-grid the map with a 3-parameter principal-component or Chebyshev dark-energy parameterization at a coarse grid; if the nexus moves by more than the grid spacing in (Ω_m, h_rd), the claim is parametrization-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim — that early-universe solutions push the DDE preference from >2σ to ~1.3σ — depends entirely on the green iso-contours in Fig. 1. These contours are computed at each grid point by fitting a 2-parameter CPL model to BAO+SN data, taking the mean and covariance of the (w0, wa) posterior, and then converting the Mahalanobis distance to (-1,0) into a sigma level via Eqs. (A1)-(A2), which assume the posterior is a 2D Gaussian. The supplementary states this is justified by 'the near-Gaussian nature of the contours throughout the parameter space,' but no diagnostics are shown. At fixed (Ω_m, h_rd) values far from the BAO+SN best fit, or where w0-wa degeneracies are strong, posteriors can be significantly non-Gaussian; a skewed or heavy-tailed posterior can make the Mahalanobis distance a poor proxy for the true evidence ratio. Additionally, the CPL form cannot capture rapid or oscillatory dark-energy features, so the map may underestimate or overestimate preferences in a way that moves the nexus. The Chebyshev reconstruction of Fig. 2 provides an independent check, but only at a single fixed point and only for Union3, not a grid-wide map. Because the paper's headline numbers and the location of the nexus are read directly off this map, the absent Gaussianity validation is the most load-bearing unexamined assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This letter argues that the apparent evidence for dynamical dark energy (DDE) from BAO and supernova data is not model-independent but is conditioned on the assumed early-universe parameters H0 rd and Ωm. The authors show geometrically that BAO and uncalibrated SN data constrain only the expansion history E(z) and the product H0 rd (Eq. 1), while the CMB angular scale provides a nearly model-independent degeneracy line in (Ωm, H0 rd). They then divide the (Ωm, H0 rd) plane into a grid, fit a CPL dark-energy model to BAO+SN at each point, and convert the distance of the posterior from w0=-1, wa=0 into a significance using the Mahalanobis distance (Eqs. A1-A2). The resulting map has a minimum ('nexus') at high H0 rd and low Ωm, close to the region preferred by early-universe Hubble-tension solutions such as varying me, WZDR, and NEDE. They conclude that early-universe solutions generically reduce the DDE preference from >2σ to about 1.3σ, and that claims of phantom crossing or thawing DDE are premature.","tokens_in":8002,"tokens_out":8922,"duration_ms":80846,"significance":"If the central claim survives scrutiny, it has a significant impact on the interpretation of DESI and SN results, reframing DDE evidence as conditional on early cosmology. The geometric decomposition in Eqs. (1)-(2) is elegant, internally consistent, and largely assumption-light; the use of three SN catalogs and an independent Chebyshev reconstruction are commendable robustness checks. The paper also makes a concrete, falsifiable prediction: future SN and BAO data should show a lower Ωm if the true model has high H0 rd. The main quantitative result, however, rests on an unvalidated Gaussian approximation, so the significance values and the exact nexus location are not yet established to the standard required.","major_comments":[{"comment":"The significance map in Fig. 1 is computed from the Mahalanobis distance of the CPL posterior to (-1,0), assuming that each (w0, wa) posterior is a two-dimensional Gaussian. The supplementary text asserts that contours are 'typically close to Gaussian' and 'near-Gaussian nature ... throughout the parameter space', but no diagnostic is shown. This is load-bearing because the paper's headline numbers (reduction from about 2.5σ to 1.3σ, and the nexus position at hrd ≈ 103 Mpc, Ωm ≈ 0.285) are read from this map. At grid points far from the BAO+SN best fit, or where the w0-wa degeneracy is strong, the posterior can be skewed or heavy-tailed, in which case the Mahalanobis distance to a single point is not the evidence ratio against a cosmological constant. I request a validation of the Gaussian assumption at several representative grid points, for example by comparing the Mahalanobis significance with the full parameter-shift metric of Ref. [11] or with an evidence-ratio calculation, and by over-plotting the actual 2D posterior contours.","section":"Supplementary material, Eqs. (A1)-(A2)"},{"comment":"The Chebyshev-polynomial check is presented as confirming the main claim, but it is performed at a single point in parameter space: hrd is fixed to 102.4 Mpc (the yellow star in Fig. 1) and only Union3 SN data are used in Fig. 2, while Fig. 6 compares only this fixed hrd with a free hrd. This does not test the shape of the green iso-contours across the (Ωm, hrd) plane, nor the location of the nexus. To support the statement 'We confirm this with a full reconstruction', a grid of Chebyshev fits (or at least a handful of representative points) would be needed. Without that, the claim of robustness to the CPL parametrization is only demonstrated at the nexus itself.","section":"Fig. 2 and Fig. 6"},{"comment":"The conclusion that 'Models that solve the Hubble tension necessarily push the investigation to trade a >2σ preference for DDE for an 0.5σ increase in the SN tension' is stronger than what the analysis demonstrates. The evidence consists of three specific model classes (varying me, WZDR, NEDE) whose best-fit positions are taken from Ref. [11], plus a heuristic argument citing Ref. [17] that typical early-universe solutions increase hrd and lower Ωm. This supports 'generically' or 'the models considered here', but not 'necessarily'. If there exist early-universe Hubble-tension solutions that do not move toward the nexus (e.g., models that simultaneously raise Ωm h²), the claim would be false. Please either provide a more rigorous proof of the universality or soften the wording.","section":"Conclusions, final paragraph"}],"minor_comments":[{"comment":"There is a typo: 'redhsift' should be 'redshift'.","section":"Introduction, first sentence"},{"comment":"'Mahalabonis' is a misspelling of 'Mahalanobis'.","section":"Supplementary material, first paragraph"},{"comment":"The text claims the preference is 'below 0.5σ' in the nexus, but the green iso-contours in Fig. 1 only label 1σ, 2σ, 3σ, 4σ, and 5σ; please add a 0.5σ contour or state the numerical minimum significance from the grid.","section":"Fig. 1 and surrounding text"},{"comment":"The notation for the chi-squared cumulative distribution is unclear ('Cχ2d' in one place, 'F χ2 d' in another); please use a single, unambiguous symbol such as F_{χ²_d} throughout.","section":"Eq. (A2)"},{"comment":"The reference for Brieden, Gil-Marín, and Verde appears malformed ('JCAP12(12)'); please correct the journal volume, issue, or article number.","section":"Ref. [24]"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of a cosmology journal and the central idea is timely. The main risk is the unvalidated Gaussian-Mahalanobis pipeline; if the authors provide the requested diagnostics, I would be willing to reconsider. The reliance on the authors' own Ref. [11] for model star positions is acceptable given the explicit cross-reference, but an independent check of those positions would strengthen the paper. There is also a slight mismatch between the abstract's 'generically' and the conclusions' 'necessarily' that should be harmonized."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's core point is right: the apparent DESI preference for dynamical dark energy (DDE) is not a model-independent fact. The authors build a significance map in (Ω_m, h_rd) space for CPL dark energy vs. ΛCDM, fitted to BAO+SN data, and show that early-universe Hubble-tension solutions (which push toward higher h_rd and lower Ω_m) land in a 'nexus' where the preference drops from >2σ to ~1.3σ. That is a genuinely useful framing, and I think it survives scrutiny.\n\nWhat's new: the explicit construction of the sigma iso-contours on this plane, and the observation that the nexus sits right on the CMB angular constraint even though no CMB data were used. The argument is robust to the choice of SN catalog (Pantheon+, Union3, DES Dovekie all give similar results), and the Chebyshev reconstruction at the fixed point gives an independent check. The geometric decomposition in Eqs. (1)-(2) is clean and the reasoning is transparent. This is a well-written letter.\n\nSoft spots, in decreasing order of concern. First, the sigma map is computed using a Mahalanobis distance from the CPL posterior to (-1,0), assuming the posterior is a 2D Gaussian. The authors assert near-Gaussianity 'throughout the parameter space' but show no diagnostics. If the w0-wa contours are skewed or multimodal near the nexus, the exact sigma levels and even the nexus location could shift. The conclusion is qualitative, so I doubt it flips, but the headline numbers deserve a check. Second, the nexus location is quoted without error bars; a dense grid plus bootstrap would help. Third, no code or parameter files are released, which slows independent verification. Fourth, the conclusion uses 'necessarily' in a couple of places where the argument actually supports 'generically' – an early-universe solution that does not move h_rd or Ω_m much would not weaken DDE evidence.\n\nThe paper is squarely aimed at people interpreting DESI results and working on the Hubble tension. It does not oversell; it explicitly says the nexus is not a silver bullet and trades DDE preference for increased SN tension. This deserves a serious referee. A referee should ask for Gaussianity diagnostics, error bars on the nexus, and a code release, but the central claim is sound.","headline":"A clean, likely-correct argument that the DESI dynamical dark energy preference is conditioned on early-universe assumptions; the main weakness is an unvalidated Gaussianity assumption in the significance map.","tokens_in":8502,"tokens_out":2473,"would_cite":true,"duration_ms":20083,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that the evidence for dynamical dark energy is conditioned on early-universe cosmology: Hubble-tension solutions shift the inferred $H_0 r_d$ and $\\Omega_m$ into a region where the preference drops from roughly 2.5σ to…","keywords":["dynamical dark energy","Hubble tension","baryon acoustic oscillations","supernovae","equation of state","CPL parametrization","sound horizon","early universe cosmology"],"falsifier":"A reader could test the paper's central claim by computing the full parameter-shift metric or nested-sampling evidence ratio for CPL versus $\\Lambda$CDM at every grid point instead of using the Gaussian Mahalanobis conversion; if the significance at the nexus rises above about 1σ, or the minimum moves noticeably, the claimed conditioning on early cosmology weakens.","tokens_in":7467,"feed_emoji":"🔭","tokens_out":14168,"duration_ms":115975,"temperature":0.7,"pith_summary":"This paper argues that the widely quoted preference for dynamical dark energy from baryon acoustic oscillations and supernovae is not a model-independent fact: it depends on what one assumes about the early universe, specifically on the inferred product of the Hubble constant and the sound horizon, $H_0 r_d$, and the matter density $\\Omega_m$. Using a grid of fixed values for these two parameters, the authors fit a two-parameter dark energy equation of state (the Chevallier-Polarski-Linder form) to BAO and supernova data and convert the distance from a cosmological constant into a significance. The significance is smallest near a region they call the 'nexus,' centered near $h r_d \\approx 103\\,\\mathrm{Mpc}$ and $\\Omega_m \\approx 0.285$ for the DES Dovekie sample. Models that resolve the Hubble tension with early-universe physics generically shift the inferred parameters toward this nexus, weakening the preference for dynamical dark energy from roughly 2.5σ to 1.3σ while slightly increasing a residual tension with supernovae. If correct, claims of phantom crossing or thawing dark energy are premature until early cosmology is pinned down.","feed_headline":"Early-universe fixes cut dark energy evidence to 1.3 sigma","feed_subtitle":"BAO and supernovae alone put the preference near 2.5 sigma; shifting to Hubble-tension solutions drops it to 1.3 sigma.","key_machinery":"The central object is the geometric degeneracy encoded in $H(z)=H_0 E(z)$: BAO data measure $\\beta(H_0 r_d)/f(z_i)$ and $(H_0 r_d)E(z_i)$, while uncalibrated supernovae measure $(1+z)^2/f(z)$, so these probes constrain only the product $H_0 r_d$ and the shape of $E(z)$. The CMB angular scale contributes the nearly model-independent band $\\theta_{\\rm CMB}=\\beta_{\\rm drag}\\beta(H_0 r_d)/f(z_{\\rm CMB})$. The paper scans a grid in ($H_0 r_d,\\Omega_m$); at each point it fits the Chevallier-Polarski-Linder parameters $w_0,w_a$ with $w(a)=w_0+w_a(1-a)$, and converts the covariance-weighted Mahalanobis distance from the cosmological-constant point ($w_0=-1,w_a=0$) into a $\\sigma$ level through a $\\chi^2_2$ distribution. The 'nexus' is the local minimum of this $\\sigma$ map. A fourth-order Chebyshev expansion of the dark energy density serves as an independent reconstruction check.","core_discovery":"The central discovery is that the level of evidence for dynamical dark energy, measured by fitting the CPL parametrization $w(a)=w_0+w_a(1-a)$ to BAO and uncalibrated supernova data, is not a single number but a field over the early-universe parameters ($H_0 r_d,\\Omega_m$). Only the product $H_0 r_d$ and the expansion shape $E(z)$ are accessible to these late-time probes, so the same data can support 'cosmological constant' or 'dynamical dark energy' depending on where the early universe places the model. The CMB angular scale provides an almost model-independent ridge in this plane, and the minimum of the dynamical-dark-energy preference, the 'nexus,' sits on that ridge without having used CMB data. Early-universe solutions to the Hubble tension, which move the joint fit to higher $H_0 r_d$ and lower $\\Omega_m$, push the fit toward that nexus and reduce the significance from about 2.5σ to 1.3σ, making the apparent preference depend on the assumed early cosmology. A fourth-order Chebyshev reconstruction of the dark energy density confirms the same weakening near the nexus.","pith_inferences":["A natural extension is that any late-time reconstruction of dark energy, not just CPL, will have its evidence reweighted by early-universe assumptions, since BAO and supernovae only fix $H_0 r_d$ and $E(z)$.","The paper's logic implies that an independent, model-free measurement of $\\Omega_m$, for example from galaxy clustering or CMB lensing, would be more decisive for settling dynamical dark energy than more low-redshift supernovae.","If future data push $\\Omega_m$ down toward the nexus, the same geometric tension now read as 'dynamical dark energy' could instead be read as early-universe physics plus a residual supernova offset.","A stress test the authors did not run is to compute full Bayesian evidence at each grid point; doing so would show whether the quantitative sigma levels, not just the qualitative nexus, survive non-Gaussian contours."],"forward_implications":["The reported preference for dynamical dark energy carries an implicit early-universe prior: with early-universe Hubble-tension solutions, it drops from roughly 2.5σ to 1.3σ.","A phantom crossing or thawing behavior inferred from BAO and supernovae is not evidence for new physics until the early-universe parameters $H_0 r_d$ and $\\Omega_m$ are independently fixed.","The nexus lies on the CMB angular constraint even though no CMB data enter its construction, linking the weak-dark-energy region to the geometric CMB ridge.","The residual low-redshift tension with supernovae grows by about 0.5σ when early-universe models move to the nexus, so the Hubble-tension fix and the dark-energy preference trade against each other.","Future measurements of $\\Omega_m$ from galaxy full-shape clustering and Lyman-$\\alpha$ data will identify where the true cosmology sits on the CMB constraint, and therefore how much dark-energy evidence survives."],"supporting_citations":[{"why":"The DESI free-form reconstruction that this paper reinterprets; it sets the sound horizon through a $\\Lambda$CDM-based early-time expression, one of the priors that conditions the evidence.","marker":"[1]"},{"why":"Supplies the fourth-order Chebyshev polynomial expansion of the dark energy density used to cross-check the CPL-based significance result.","marker":"[2]"},{"why":"A recent dynamical-dark-energy analysis that imposes a $\\Lambda$CDM prior on $r_d$, illustrating the early-universe assumption whose effects the paper traces.","marker":"[5]"},{"why":"Provides the best-fit ($h r_d,\\Omega_m$) values and the roughly 1.3σ dynamical-dark-energy preference for early-universe models, the comparison point for the nexus calculation.","marker":"[11]"},{"why":"Defines the phase-shift factor $\\beta$ and the template-fit correction entering the BAO and CMB geometric relations.","marker":"[12]"},{"why":"Supplies the DESI DR2 BAO constraints that define the BAO contour in ($\\Omega_m,h r_d$) space.","marker":"[13]"},{"why":"Explains why early-universe Hubble-tension solutions that raise $h$ without a proportional rise in $\\Omega_m h^2$ end up at lower $\\Omega_m$ and higher $h r_d$.","marker":"[17]"},{"why":"Provides the Mahalanobis distance used to convert CPL ($w_0,w_a$) contour offsets into a sigma-level preference.","marker":"[20]"},{"why":"Supplies the DES Y5 Dovekie uncalibrated supernova sample used for the main significance map.","marker":"[21]"},{"why":"Supplies the Union3 supernova sample used for the robustness check and for the Chebyshev reconstruction in Figure 2.","marker":"[23]"}],"fun_headline_variants":["Early cosmology slices dark energy evidence to 1.3σ","Hubble fixes trim dark energy preference to 1.3σ","Dynamical dark energy support rides on early universe","Early universe dictates strength of dark energy evidence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes the dark-energy fit contours are nearly bell-shaped at every grid point, so the Mahalanobis distance can be read as a sigma level; the paper states this is a good approximation but does not verify it at every point, and a failure there would shift the reported significance and the exact nexus location.","fun_headline_variants_meta":{"raw":{"variants":["Early cosmology slices dark energy evidence to 1.3σ","Hubble fixes trim dark energy preference to 1.3σ","Dynamical dark energy support rides on early universe","Early universe dictates strength of dark energy evidence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000613,"raw_usage":{"total_tokens":2836,"prompt_tokens":917,"completion_tokens":1919,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":1854}},"tokens_in":533,"tokens_out":1919,"duration_ms":15163,"temperature":1.0,"reasoning_tokens":1854,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:26:40.351166+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A reader could test the paper's central claim by computing the full parameter-shift metric or nested-sampling evidence ratio for CPL versus $\\Lambda$CDM at every grid point instead of using the Gaussian Mahalanobis conversion; if the significance at the nexus rises above about 1σ, or the minimum moves noticeably, the claimed conditioning on early cosmology weakens.","supporting_citations":[{"cited_title":"Model-Independent Indication for a Localized Anomaly in the Late-Time Expansion History","cited_arxiv_id":"2607.13009","evidence_quote":"A recent dynamical-dark-energy analysis that imposes a $\\Lambda$CDM prior on $r_d$, illustrating the early-universe assumption whose effects the paper traces."},{"cited_title":"The $H_0$ world cup. II. A comprehensive competition between proposed Hubble tension solutions","cited_arxiv_id":"2607.13283","evidence_quote":"Provides the best-fit ($h r_d,\\Omega_m$) values and the roughly 1.3σ dynamical-dark-energy preference for early-universe models, the comparison point for the nexus calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Mahalanobis distance used to convert CPL ($w_0,w_a$) contour offsets into a sigma-level preference."},{"cited_title":"Exploring New Frontiers in Cosmology","cited_arxiv_id":null,"evidence_quote":"Supplies the DES Y5 Dovekie uncalibrated supernova sample used for the main significance map."}],"review_version":1}