{"id":"b9f58a54-ec0e-4f3a-8e67-f108885161f5","arxiv_id":"2601.05646","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A microphysically motivated Q∝ρ² dark-sector interaction is constrained by cosmological data, giving upper limits A<7.6×10⁻²⁵, B<0.048, and H₀=67.71±0.65 km/s/Mpc.","lead":"This paper proposes a new way to write the interaction between dark matter and dark energy, based on particle collisions rather than an arbitrary formula, and fits it to recent cosmological data. The best fit says the interaction is extremely weak and the Hubble constant is 67.7 km/s/Mpc, in line with early-Universe measurements.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boltzmann-to-fluid mapping fails for w≈−1 dark energy; B limit is not a physical cross-section bound","rationale":"The reader's conditional verdict is appropriate. The most robust result—that a Q∝ρ² interaction is tightly constrained by background data—likely survives. However, the headline cross-section bound depends on the microphysical mapping. For non-relativistic DM, ρ=mn is standard; for DE with w≈−1, it is not. Section III's assumption that ϕ is a light scalar or condensate does not establish a particle number or mass, and the appeal to Eq. (4) does not justify writing the reverse process as an independent Bρ_DE² term. Detailed balance further undermines the two-parameter independence. Therefore the B limit should not be described as a thermally-averaged annihilation cross-section per unit mass without an explicit field-theoretic derivation. Since this is exactly the reader's weakest-assumption concern, I agree with the conditional verdict and recommend no change: the paper should be accepted only if the cross-section interpretation is either rigorously derived or substantially softened.","tokens_in":11915,"tokens_out":8630,"duration_ms":99828,"concrete_test":"Build an explicit model, e.g. L = L_χ + ½(∂ϕ)² − V(ϕ) − λχ²ϕ² with V(ϕ)=½m_ϕ²ϕ², and derive the background energy-transfer Q from the full coupled field/Boltzmann equations without setting ρ_ϕ=m_ϕ n_ϕ. Check whether Q can be cast exactly as Eq. (6) with constant A and B, and whether B equals ⟨σ|v|⟩/m_ϕ for a w≈−1 solution. Also test detailed balance: if A and B are related by the underlying 2→2 process, the two-parameter fit is not microphysically reversible. Failure of this reduction confirms that the B limit is not a cross-section bound.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that A and B are bounds on thermally-averaged annihilation cross-sections per unit mass rests on Eq. (5), obtained from Eq. (4) via ρ=mn. This substitution is valid for non-relativistic DM, but the paper applies the same replacement to DE, which it models as a light scalar field or condensate with w≈−1 (Section III). For such a field, a particle number density n and a mass m entering ρ=mn are not defined: the energy density is dominated by the potential, ρ≈V(ϕ), and the particle-number current is not conserved. Excluding a cosmological constant because it lacks particle-like excitations does not by itself supply a particle interpretation for a w≈−1 scalar. Moreover, Eq. (6) does not follow from Eq. (4) as written: for one 2→2 process χχ↔ϕϕ, detailed balance relates forward and reverse rates, so A and B cannot be treated as independent, and the reverse creation term for DM should involve the DM equilibrium density, not ρ_DE² with an arbitrary coefficient. Thus the B-bound—and the mutual independence of A and B—cannot be interpreted as a physical dark-energy scattering cross-section limit. The background constraint B<0.048 may survive as a phenomenological bound, but the advertised cross-section conclusion for the dark-energy sector is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an interacting dark energy–dark matter model with interaction term Q = -A (H0/ρc,0) ρ_DM^2 + B (H0/ρc,0) ρ_DE^2, motivated by a Boltzmann collision term for the reversible process χχ ↔ φφ. The model is fit to Pantheon+, cosmic chronometers, DESI DR2 BAO, and Planck distance priors. The combined analysis yields H0 = 67.71 ± 0.65 km/s/Mpc and 95% upper limits A < 7.586×10^-25 and B < 0.048. The authors claim these limits translate into bounds on thermally-averaged annihilation cross-sections per unit mass, and that the model fits the data comparably to flat ΛCDM.","tokens_in":12311,"tokens_out":7973,"duration_ms":88236,"significance":"If the microphysical mapping were sound, this would be a useful step toward particle-physics-motivated interacting dark energy models, avoiding the H-dependent interaction terms that suffer from large-scale instabilities. The dataset combination is current and the paper is transparent about many numerical details, including the bimodality of the B posterior and per-dataset chi-square values. However, the central cross-section interpretation is not established: the Boltzmann-to-fluid mapping fails for a w ≈ -1 dark energy field, the forward and reverse rates are not independent, and no numerical cross-section limit is actually computed. The results are plausible as phenomenological background constraints on a ρ² interaction, but the advertised claim of 'first bounds on dark-sector scattering cross sections' overreaches.","major_comments":[{"comment":"Equation (6) is not a valid two-species reduction of the Boltzmann equation. Equation (5) uses ρ = m n, which is appropriate only for non-relativistic, dilute particles. Section III explicitly models dark energy as a light scalar field or condensate with w ≈ -1; for such a field the energy density is potential-dominated and there is no conserved particle number, so n_DE and m_DE entering Eq. (5) are undefined. The B term therefore has no particle-scattering interpretation, and the claimed dark-energy annihilation cross-section bound is unsupported. In addition, for a single χχ ↔ φφ process, forward and reverse rates are related by detailed balance; treating A and B as independent positive coefficients requires a microphysical justification that the paper does not provide.","section":"Section II–III, Eqs. (4)–(6)"},{"comment":"No cross-section bound is ever presented. The title and abstract advertise a limit on the thermally-averaged annihilation cross-section per unit mass, but the paper only reports constraints on the dimensionless coefficients A and B. If Eqs. (5)–(6) are accepted, the conversion is ⟨σv⟩/m_χ = A H0/ρc,0, yet no numerical value with units is given anywhere. The authors should either state the derived cross-section limits explicitly with units and error propagation, or revise the central claim to be purely a bound on the coefficients A and B.","section":"Title, Abstract, Section V, Conclusion"},{"comment":"The posterior for B is bimodal, with a best-fit peak at the edge of the sampled range. The paper itself describes a 'sharp tip at the end of the sample' and shows a two-peaked likelihood distribution. A 95% upper limit B < 0.048 obtained from such a posterior is not a stable summary statistic and is not robust evidence of a constraint. A profile-likelihood analysis or a prior-sensitivity check is needed before B can be quoted as a headline upper limit.","section":"Appendix A, Figures 2–3"},{"comment":"The reported reduced chi-square for the full combination is inconsistent: the main text states χ²ν = 0.8889, while Table III gives 0.8805. One of these values is wrong. Additionally, Table III shows that the IDE model fits BAO alone and BAO+PCMB substantially worse than ΛCDM (χ²ν = 1.5050 vs 1.1426 and 1.7362 vs 1.3379, respectively). This is acknowledged only in passing, yet it weakens the conclusion that the model 'fits the observational data with high precision.' The BAO tension should be quantified and discussed in the context of the combined fit.","section":"Section V and Appendix A, Table III"}],"minor_comments":[{"comment":"The text says 'log10 A < 7.586×10^-25', which is dimensionally inconsistent. It should read A < 7.586×10^-25, or equivalently log10 A < -24.12.","section":"Section V"},{"comment":"The header 'Top 5% Max' is unclear; the standard definition of the 95% upper limit should be stated explicitly.","section":"Table II"},{"comment":"The vertical axis is labeled 'Log Probability (ln)' but the text refers to lnL (log-likelihood); the labeling should be made consistent and precise.","section":"Figures 2–3"},{"comment":"The normalization by H0/ρc,0 is arbitrary from a microphysical viewpoint. The authors should state clearly that A and B are defined relative to this cosmological normalization, and that this choice affects the numerical values of the quoted bounds.","section":"Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The paper's observational analysis is usable and the DM-side A constraint can likely be salvaged, but the advertised cross-section bound for the dark-energy sector is not supported by the current derivation. A revision that restricts the cross-section claim to DM, properly derives the two-species Boltzmann reduction, computes the actual cross-section limits, and addresses the B posterior would be needed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth a look if you care about interacting dark energy, but the main advertised result — cross-section limits — is not actually delivered. What is new is the specific Boltzmann-motivated form Q = -A(H0/ρc,0)ρ_DM^2 + B(H0/ρc,0)ρ_DE^2 with independent DM and DE self-annihilation coefficients, and the first constraints on it from the Pantheon+ / CC / DESI DR2 / Planck-prior combination. Those constraints are honestly presented: the interaction is consistent with zero, H0 comes out 67.71±0.65, and the fit is essentially as good as ΛCDM. So as a phenomenological bound on a plausible dark-sector interaction, the paper is fine and probably correct in its main qualitative conclusion: any such interaction must be very small at late times.\n\nThe soft spots are real, and they are concentrated in the bridge from the fitted parameters to the physical cross-sections. Eq. (5) is obtained from the Boltzmann equation via ρ = m n, which is fine for non-relativistic DM but not for dark energy with w ≈ -1, whether you call it a light scalar or a condensate. For such a field there is no well-defined particle number density entering ρ = m n, so the B limit as a thermally-averaged cross-section is unsupported. The paper essentially acknowledges this by insisting DE must be dynamical, but that does not fix the substitution. There is also a detailed-balance problem: if the interaction is one 2→2 process χχ ↔ ϕϕ, the forward and reverse rates are related by equilibrium abundances, so treating A and B as independent coefficients with an arbitrary B is not microphysically motivated.\n\nThere are also smaller numerical issues. The reported chi-square for the full combination changes between the abstract (0.8889), the conclusions (0.8805), and the appendix table (0.8805); the B posterior is bimodal and edge-peaked; and the BAO and BAO+PCMB fits are noticeably worse than ΛCDM. None of these kill the paper as a phenomenological study, but they need cleaning up.\n\nWho is this for? Someone working on dark-sector interactions or H0 tension who wants a quick sense of what current background data say about a quadratic coupling. The paper deserves a serious referee, not because the cross-section claim is right, but because the model is new enough and the constraints are useful enough. My recommendation: send it to peer review, but the referee should insist the authors either compute the cross-section bounds properly for a well-defined particle model, or re-scope the paper as constraints on the dimensionless coefficients A and B only.","headline":"A useful phenomenological constraint on a quadratic dark-sector interaction, but the advertised cross-section bounds do not follow from the model's own equations, especially for dark energy.","tokens_in":12801,"tokens_out":1653,"would_cite":false,"duration_ms":20443,"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 paper claims that dark matter and dark energy can exchange energy through reversible two-body collisions, and that the resulting quadratic interaction fits all background cosmological data while yielding the first tight bounds on dark-s","keywords":["interacting dark energy","Boltzmann equation","dark matter annihilation","dark energy microphysics","Hubble constant tension","cosmological parameter constraints","baryon acoustic oscillations","CMB distance priors"],"falsifier":"Solve the first-order perturbation equations for this Q ∝ ρ^2 interaction: if dark-energy density perturbations grow without bound on large scales in the early universe—as happens for many Hubble-proportional interacting-dark-energy models—the model is ruled out even though it fits the background data.","tokens_in":11824,"feed_emoji":"🌌","tokens_out":7212,"duration_ms":71912,"temperature":0.7,"pith_summary":"This paper tries to put interacting dark energy on a microphysical footing rather than choosing the interaction by dimensional analysis. It proposes that dark matter and dark energy exchange energy through local, reversible 2-to-2 collisions, which at the background level gives an interaction proportional to the square of the energy densities rather than to the Hubble rate. Fitting the model to combined late-time and CMB distance-prior data returns a Hubble constant of 67.71 ± 0.65 km/s/Mpc and 95% upper limits on the two interaction coefficients, which the paper interprets as the first observational bounds on thermally averaged dark-sector annihilation cross sections per unit mass. If correct, any dark-matter annihilation into dark-energy states must be extremely weak, while the reverse process is only constrained at the few-percent level. This matters because it offers a collision-based alternative to the phenomenological interacting-dark-energy models that are prone to instabilities.","feed_headline":"Cosmic data bound dark matter–dark energy collisions","feed_subtitle":"A collision-based interaction model fits the data but forces the dark-matter coupling below 10^-24.","key_machinery":"The load-bearing object is the interaction term Q = -A (H0/ρc,0) ρ_DM^2 + B (H0/ρc,0) ρ_DE^2, obtained from the standard Boltzmann annihilation term by substituting ρ = m n for each dark-sector species. This coarse-grained term turns the dark-matter–dark-energy coupling into a local, number-changing collision process; A and B are dimensionless coefficients that the paper reads as thermally averaged cross-sections per unit mass. Supporting the term is an effective-field-theory description with quartic operators (such as χ^2 φ^2) that realize the reversible χ+χ ↔ φ+φ channel and exclude linear decay interactions.","core_discovery":"The central claim is that a bottom-up interacting-dark-energy model with interaction Q = -A (H0/ρc,0) ρ_DM^2 + B (H0/ρc,0) ρ_DE^2 is viable. The quadratic form follows from the Boltzmann collision term for a 2-to-2 process χ+χ ↔ φ+φ after writing number densities as mass densities. Using supernova distances, cosmic chronometers, baryon acoustic oscillations, and CMB distance priors, the model fits the data as well as flat ΛCDM, with H0 = 67.71 ± 0.65 km/s/Mpc and 95% upper limits A < 7.586 × 10^-25 and B < 0.048. The paper's key interpretive step is that these dimensionless coefficients, together with the expansion rate, translate directly into limits on the thermally averaged annihilation c","pith_inferences":["Editorial extension: if the interaction really is quadratic in densities, its effects scale strongly with redshift, so nucleosynthesis or recombination-era data may provide much stronger tests than the late-time background probes used here.","Editorial extension: the B-bound's lack of a particle interpretation is a genuine open question—if dark energy is a bare cosmological constant, the central cross-section claim reduces to a constraint on an effective coefficient, not a physical scattering rate.","Editorial extension: the same Boltzmann-derived interaction could be applied to dark matter–dark radiation or dark matter–dark matter self-interactions in the early universe, where the resulting relic-density changes would be observable."],"forward_implications":["Dark-matter self-annihilation into dark-energy states is forced to be extremely weak, with A below 7.6 × 10^-25, making such a process cosmologically negligible.","The dark-energy-side coupling is only bounded at B < 0.048, so the reverse process could still be active at the few-percent level.","The model prefers H0 ≈ 67.7 km/s/Mpc, aligning with early-universe estimates rather than the higher local values, so it does not resolve the Hubble tension by raising H0.","The interaction coefficients are strongly correlated with the sound horizon; CMB distance priors fix the sound horizon near 147.7 Mpc and thereby suppress the allowed interaction range.","The model's fit quality is nearly identical to flat ΛCDM, so any detectable signature must come from epochs or observables where the quadratic interaction term is larger, not from the late-time background expansion alone."],"fun_headline_variants":["Dark matter–dark energy scattering tightly bounded","Collision model sets strict dark sector bounds","Cosmic data pin down dark sector scattering","Dark matter–dark energy collisions face tight limits","Dark energy interactions constrained by data"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The cross-section interpretation rests on treating dark energy as a gas of non-relativistic particles with a definite mass, so that its energy density equals mass times number density; if dark energy is a cosmological constant or a field whose energy cannot be counted as particles, the B term and its bound lose their microphysical meaning.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter–dark energy scattering tightly bounded","Collision model sets strict dark sector bounds","Cosmic data pin down dark sector scattering","Dark matter–dark energy collisions face tight limits","Dark energy interactions constrained by data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000817,"raw_usage":{"total_tokens":3515,"prompt_tokens":945,"completion_tokens":2570,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":689,"completion_tokens_details":{"reasoning_tokens":2505}},"tokens_in":689,"tokens_out":2570,"duration_ms":19090,"temperature":1.0,"reasoning_tokens":2505,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T11:34:59.006373+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the first-order perturbation equations for this Q ∝ ρ^2 interaction: if dark-energy density perturbations grow without bound on large scales in the early universe—as happens for many Hubble-proportional interacting-dark-energy models—the model is ruled out even though it fits the background data.","supporting_citations":[],"review_version":1}