{"id":"e16c7d04-1ae8-4072-82ee-bd11c8ca590f","arxiv_id":"2509.09624","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Constraints from Planck, DESI DR2, and supernovae favor LambdaCDM over a coupled quintessence model with warm dark matter, despite a large DIC improvement for the interacting model.","lead":"This paper constrains a quintessence dark energy model that interacts with dark matter, using Planck, DESI DR2, and three supernova datasets. The key result is a split verdict: a fit statistic prefers the interacting model with warm dark matter, but Bayesian evidence still favors the standard LambdaCDM model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The WDM perturbation equation (10) omits velocity-divergence and pressure-gradient terms needed for a fluid with ω_dm≠0; the reported 5σ detection of ω_dm may be an artifact of this incomplete fluid treatment.","rationale":"The paper's central quantitative claim is that Bayesian evidence favors ΛCDM over the coupled quintessence model, with the secondary claim being a 5σ detection of a non-zero dark matter equation of state ω_dm. The reader identified the barotropic-fluid treatment of WDM as the weakest assumption. I agree, but I sharpen the concern: even within a barotropic description, the perturbation equation Eq. (10) is internally inconsistent for ω_dm ≠ 0 because it omits the velocity-divergence and pressure-gradient terms that must appear for a fluid with c_s^2 = ω_dm. This is not merely a microphysical approximation; it is a concrete technical error that affects the predicted CMB spectra. If the implementation in the modified CAMB uses Eq. (10) as written, the constraints on ω_dm are unreliable. If the implementation actually includes the full equations, then the published equation is misleading, and the code should be released to verify. Either way, the 5σ detection and the DIC-based preference for CQ+WDM are not solid. The Bayesian evidence result favoring ΛCDM is less likely to be overturned because it is very strong for CQ+CDM and moderate for CQ+WDM, but a correct WDM perturbation treatment could change the WDM evidence. The reader's CONDITIONAL verdict remains appropriate, conditional on clarifying the perturbation implementation and either fixing it or demonstrating that the missing terms do not affect the conclusions. I therefore recommend UNCHANGED.","tokens_in":17214,"tokens_out":13939,"duration_ms":147933,"concrete_test":"Re-implement the DM component in the modified CAMB using the full coupled fluid perturbation equations (including θ_c and the c_s^2 k^2 δ/(1+w) term) with c_s^2 = ω_dm, and re-run the Pl+DESI+DESY5 analysis. If the 68% interval on ω_dm shifts by more than 2σ from 0.0025 ± 0.0005, or the Δχ²_MAP improvement over ΛCDM changes by more than 10%, then the paper's WDM detection and its evidence comparison are not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section II, after Eq. (10), the paper models warm dark matter as a barotropic fluid with c_s^2 = ω_dm, but the perturbation equation actually used, δ'_c = -(1+ω_dm) h'/2 + m φ' δ_c, drops the velocity-divergence term (θ_c) and the pressure-gradient term that are required for a fluid with a non-zero sound speed. In the synchronous gauge, the standard Ma–Bertschinger equations for constant w and c_s^2 include δ' = -(1+w)(θ + h'/2) and θ' = -H(1-3w)θ + c_s^2 k^2 δ/(1+w) - k^2 σ (plus coupling source terms). Setting θ_c = 0 is inconsistent with the Euler equation for non-zero c_s^2 at finite k. Thus the CMB power spectra and the inferred constraints on ω_dm are not those of a self-consistent warm fluid. The reported '5σ' detection of ω_dm and the large Δχ² improvements that drive the DIC preference could therefore be artifacts of an incomplete perturbation treatment, casting doubt on the secondary claim and on the validity of the CQ+WDM comparison itself.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript constrains an exponential-potential quintessence field coupled to dark matter through F(φ)=F0 e^{mφ}, considering both cold dark matter and a constant-EoS 'warm' dark matter fluid. The analysis uses Planck CMB likelihoods (including PR4 lensing), DESI DR2 BAO, and three SNe compilations (PantheonPlus, Union3, DESY5), implemented in a modified CAMB with Cobaya. The paper reports parameter constraints (Table II), goodness-of-fit and DIC comparisons (Table III), and Bayesian evidence ratios from MCEvidence (Table IV). The central quantitative results are: by DIC, the coupled-quintessence+warm-dark-matter model is strongly preferred over ΛCDM for combined datasets (ΔDIC ≈ −14 to −20), whereas by Bayesian evidence ΛCDM is preferred (lnB between −2.05 and −5.86 for WDM, and between −10.7 and −12.5 for CDM). The authors also report a ~5σ nonzero ω_dm for Pl+DESI+DESY5. The discrepancy between DIC and Bayesian evidence is acknowledged and discussed transparently.","tokens_in":17600,"tokens_out":9373,"duration_ms":104484,"significance":"If the warm-dark-matter modeling were complete, this would be a timely and useful test of a theoretically motivated interacting dark energy model against the latest CMB, BAO, and SN data, and the explicit comparison of DIC versus Bayesian evidence is a valuable cautionary exercise. The computational setup (modified CAMB, Cobaya, MCEvidence, GetDist) is standard, and the paper is transparent about priors, convergence criteria, and the model-selection discrepancy. However, the perturbation equations used for the warm dark matter fluid are incomplete (see major comments), so the quantitative WDM claims — including the 5σ ω_dm detection and the DIC preference for CQ+WDM — are not yet supported. The cold-dark-matter part is less affected, but the advertised WDM results need substantive revision.","major_comments":[{"comment":"The perturbation equation used for the WDM fluid is incomplete. For a barotropic fluid with ω_dm ≠ 0 and c_s^2 = ω_dm, the synchronous-gauge density and velocity perturbations obey δ'_c = −(1+ω_dm)(θ_c + h'/2) + coupling terms and θ'_c = −H(1−3ω_dm)θ_c + c_s^2 k^2 δ_c/(1+ω_dm) + coupling terms. Equation (10) drops the (1+ω_dm)θ_c term and provides no Euler equation for θ_c. Setting θ_c = 0 is only valid for pressureless CDM, not for a finite-sound-speed fluid; pressure gradients source θ_c at finite k, which then feeds back into δ_c and the CMB spectra. The Δχ², ΔDIC, and ω_dm constraints in Tables II–III and the reported 5σ detection in Section IV may therefore be artifacts of the incomplete fluid treatment. The full coupled fluid equations, or a proper WDM treatment with free-streaming, should be implemented and the WDM analysis redone.","section":"Section II, Eq. (10)"},{"comment":"The statement 'In this gauge, the velocity perturbation vanishes [108]' mischaracterizes the synchronous gauge of Ma & Bertschinger; in that gauge only pressureless CDM can be comoving, while a fluid with ω_dm ≠ 0 and nonzero sound speed generally has a nonvanishing velocity divergence. Thus the citation to [108] does not justify omitting θ_c. Relatedly, the interpretation of the constrained constant EoS ω_dm as 'warm dark matter' overreaches: a barotropic fluid with constant ω_dm is not a warm dark matter particle model and ignores free-streaming and viscosity. The authors should either adopt a microphysically motivated WDM description or clearly frame the result as constraints on a phenomenological dark matter fluid.","section":"Section II, after Eq. (10); Section IV"}],"minor_comments":[{"comment":"The equation is attributed to Ref. [109], but the equations in that reference include the θ_c terms; please verify that the reproduction is accurate or clarify which simplified limit is intended.","section":"Section II, Eq. (10)"},{"comment":"Typos and notation: 'Chaplyn gas' should be 'Chaplygin gas'; 'CamSpeclikelihood' should be 'CamSpec'; the symbol for the dark matter equation of state appears as both ω_dm and w_dm; please harmonize.","section":"Throughout"},{"comment":"For the Pl dataset the posterior for λ hits the lower boundary of the prior; this is mentioned but the robustness of the constraints under an extended prior (e.g., negative λ) is not discussed beyond the phantom-divide argument. A short comment on the effect of the chosen prior range would be useful.","section":"Section IV, Table II"},{"comment":"The paper states initial conditions are chosen 'in accordance with the values for the same obtained from a ΛCDM description' at lna=−7. A brief justification of the consistency of these initial conditions with the coupled model's Friedmann constraint would improve reproducibility.","section":"Section III"}],"recommendation":"major_revision","confidential_remarks":"The paper is a standard parameter-constraint analysis, and the main advertised WDM result rests on an incomplete perturbation treatment. I do not see a novelty or scope problem, but the WDM section needs substantive rework. The honest reporting of the DIC/Bayes-factor discrepancy is a strength; however, the 5σ framing in the abstract-adjacent text may be overinterpreted given that the model as a whole is disfavored by Bayesian evidence. The cold-dark-matter analysis appears sound and could be publishable even if the WDM part were removed, but as submitted the manuscript needs major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a solid constraints paper that re-fits a known coupled quintessence model to current data (full Planck, DESI DR2 BAO, three SNe compilations) and reports the model comparison honestly. The central claim—Bayesian evidence still favors ΛCDM—is believable. The DIC vs. Bayes discrepancy is handled transparently, and the MCMC setup is standard. The novelty is modest: the model is from the authors' earlier paper, and the main addition is a free ω_dm plus updated likelihoods.\n\nThe real problem is the WDM perturbation treatment. The paper models warm dark matter as a barotropic fluid with c_s^2 = ω_dm, but eq. (10) uses δ'_c = -(1+ω_dm) h'/2 + m φ' δ_c and simply drops the velocity-divergence term. There is no Euler equation. Setting θ_c = 0 is only consistent for CDM; for a fluid with nonzero sound speed, the pressure-gradient term in the Euler equation forces θ_c to evolve. I think the stress-test is right: the CMB spectra and the inferred ω_dm constraints are not those of a self-consistent warm fluid. The reported 5σ detection of ω_dm and the large Δχ² improvements that drive the DIC preference should be regarded as artifacts until this is fixed with the full Ma–Bertschinger fluid equations or a proper WDM implementation.\n\nMinor issues: eq. (16) has a term +Ω_φ Ω_r that looks like a missing factor of 3; λ sits at the prior boundary and should be flagged as prior-dominated. The modified CAMB and chains are not public, so the results are not reproducible as-is.\n\nThe Bayesian model-comparison conclusion is less affected—the evidence penalty is large enough that ΛCDM probably still wins—but the WDM row of Table IV would likely shift once perturbations are corrected. I'd be careful citing the WDM detection until that's addressed.\n\nBottom line: worth sending to peer review, but the revision needs to fix the fluid perturbations, release the code, and temper the 5σ language. The audience is people constraining interacting dark energy with DESI and SNe; I'd bring it to reading group with the perturbation caveat on the table.","headline":"Competent re-analysis of a known coupled quintessence model with DESI DR2 and full Planck; the Bayesian-evidence conclusion favoring ΛCDM is credible, but the WDM perturbation treatment is incomplete and the 5σ ω_dm detection is likely an artifact.","tokens_in":18063,"tokens_out":7477,"would_cite":false,"duration_ms":86820,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","85A40"],"pacs":["95.35.+d","95.36.+x","98.80.-k"],"model":"deepseek-v4-flash","headline":"Coupled quintessence with warm dark matter improves the fit to CMB, BAO, and supernovae data, but Bayesian model comparison still prefers LambdaCDM.","keywords":["coupled quintessence","dark energy-dark matter interaction","warm dark matter","Bayesian model selection","DESI DR2","Planck CMB","supernovae","dark matter equation of state"],"falsifier":"Compute the full linear matter power spectrum with a proper warm-dark-matter transfer function (free-streaming with finite velocity dispersion) for the same parameters; if the chi-square improvement over LambdaCDM disappears or the inferred omega_dm becomes consistent with zero, the paper's warmness detection is an artifact of the fluid approximation. Alternatively, a direct independent measurement of omega_dm from 21-cm or Lyman-alpha forest data at high significance would settle the question.","tokens_in":17084,"feed_emoji":"🌌","tokens_out":7025,"duration_ms":71333,"temperature":0.7,"pith_summary":"The paper asks whether a dark energy scalar field that exchanges energy with dark matter (coupled quintessence) can survive contact with the latest cosmological data: Planck CMB, DESI DR2 baryon acoustic oscillations, and three independent supernova compilations. Its main result is that when the dark matter is allowed to be warm (a small, non-zero equation of state), the model fits the data much better than LambdaCDM does, improving the maximum-likelihood chi-square by up to about 24 units and yielding a 4-5 sigma hint of a warm component with omega_dm ~ 0.0025. Yet when model complexity is penalized via Bayesian evidence, LambdaCDM is preferred in every dataset combination, with the coupled-cold-dark-matter version decisively disfavored and the warm version only moderately or weakly disfavored. A secondary point is that the simpler DIC criterion would have declared the warm coupled model a strong winner, illustrating how the choice of model-selection statistic changes the verdict.","feed_headline":"Bayesian evidence says LambdaCDM beats coupled quintessence","feed_subtitle":"Warm dark matter improves the fit but the extra parameters are not worth their Bayesian cost.","key_machinery":"The argument runs through a dynamical system for a canonical scalar field with exponential potential V proportional to e^{-lambda phi}, coupled to dark matter via the interaction Q_nu = F_{,phi} rho_dm grad_nu phi with F(phi) = F0 e^{m phi}. In the synchronous gauge the dark matter density perturbation obeys delta_c' = -(1+omega_dm) h'/2 + m phi' delta_c, which translates the coupling into a modified growth rate. The warm dark matter enters as a barotropic fluid with constant equation of state omega_dm and adiabatic sound speed c_s^2 = omega_dm. Solving these background and perturbation equations inside a Boltzmann solver and sampling the 8-9 parameters lets the paper compute chi^2, DIC, and","core_discovery":"On the paper's own terms, the central discovery is quantitative: across all dataset combinations considered, the log Bayes factor lnB for the coupled quintessence model relative to LambdaCDM is negative, ranging from -2.05 (Planck+DESI+DESY5, warm dark matter) to -12.5 (Planck+DESI, cold dark matter). LambdaCDM is therefore preferred, though the warm-dark-matter version comes close to being competitive. The paper also reports a positive dark matter equation of state, omega_dm = 0.0025 +/- 0.0005 at about 5-sigma for Planck+DESI+DESY5, and shows that warmth is what drives the improved fit: adding the WDM parameter lowers chi^2_MAP by roughly 24 relative to LambdaCDM while CDM versions of the","pith_inferences":["A microphysical warm-dark-matter implementation with free-streaming would likely change both the sound speed and the small-scale power spectrum, so the 4-5 sigma omega_dm detection should be re-tested with a non-fluid treatment before being interpreted as physical warmness.","The paper's prior keeps lambda non-negative so that the field stays quintessent; allowing the phantom side might remove the prior wall and change the Bayesian comparison, which is a natural extension not covered by the paper.","The same datasets could be used to score the alternative interaction forms the paper lists as future work; given the Bayes penalty for extra parameters, only interactions that deliver a large chi-square gain per parameter would plausibly beat LambdaCDM."],"forward_implications":["If the Bayesian evidence is taken at face value, the combined CMB+BAO+SN data do not require a dark-energy-dark-matter interaction; LambdaCDM remains the simplest viable description.","The warm dark matter equation of state is consistently positive and away from zero at 4-5 sigma for combined datasets, which is the strongest internal hint of a non-cold component, even though the full model is disfavored.","The cold dark matter branch of the coupled model is decisively disfavored once BAO data are added, ruling out substantial interaction strengths for CDM.","The DIC-versus-Bayes disagreement in the WDM branch implies that claims of a preferred interacting model should be checked with Bayesian evidence, since the two criteria disagree sharply."],"fun_headline_variants":["LambdaCDM wins Bayesian showdown with coupled quintessence","Bayes says no to quintessence, even with warm dark matter","Coupled quintessence loses to LambdaCDM in Bayesian analysis","Warm dark matter helps quintessence but Bayes still favors LambdaCDM","Positive dark matter EoS hints at warmth, but LambdaCDM preferred"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The warm dark matter is modeled as a perfect barotropic fluid with constant equation of state and adiabatic sound speed c_s^2 = omega_dm, so the reported nonzero warmness and fit improvement assume that dark matter behaves as such a fluid rather than as a free-streaming particle.","fun_headline_variants_meta":{"raw":{"variants":["LambdaCDM wins Bayesian showdown with coupled quintessence","Bayes says no to quintessence, even with warm dark matter","Coupled quintessence loses to LambdaCDM in Bayesian analysis","Warm dark matter helps quintessence but Bayes still favors LambdaCDM","Positive dark matter EoS hints at warmth, but LambdaCDM preferred"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000262,"raw_usage":{"total_tokens":1394,"prompt_tokens":664,"completion_tokens":730,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":408,"completion_tokens_details":{"reasoning_tokens":647}},"tokens_in":408,"tokens_out":730,"duration_ms":7475,"temperature":1.0,"reasoning_tokens":647,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T18:44:40.006948+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full linear matter power spectrum with a proper warm-dark-matter transfer function (free-streaming with finite velocity dispersion) for the same parameters; if the chi-square improvement over LambdaCDM disappears or the inferred omega_dm becomes consistent with zero, the paper's warmness detection is an artifact of the fluid approximation. Alternatively, a direct independent measurement of omega_dm from 21-cm or Lyman-alpha forest data at high significance would settle the question.","supporting_citations":[],"review_version":1}