Pith. sign in

REVIEW 4 major objections 9 minor 20 references

Interacting dark energy models

T0 review · 4 major / 9 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Dark energy feeding dark matter at a small constant rate can keep the two densities comparable over cosmic time.

desk verdict The MCMC fits look fine internally, but the paper's own equations kill its headline claim about the coincidence problem. read the letter →

arxiv 2412.09024 v1 pith:TIN4TO4R submitted 2024-12-12 astro-ph.CO gr-qc

Robert Rugg , Shambel Sahlu , Amare Abebe This is my paper · ORCID
classification astro-ph.COgr-qc PACS 95.35.+d95.36.+x98.80.-k
keywords interactingdarkenergycoincidenceproblemmatter-darkcouplingPantheonsupernovaecosmicchronometersbaryonacousticoscillationsMarkovchainMonteCarloHubbleconstant
topics Dark EnergyDark Matter
open problems Dark MatterDark Energy
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to establish that a small energy transfer between dark matter and dark energy can repair the coincidence problem, the puzzling fact that the two densities are comparable only in the present epoch. Two linear interaction models are fitted to the Pantheon supernova sample and a compiled set of Hubble-parameter measurements. In the $\eta$ model, where the exchange rate is proportional to the dark-energy density, the combined data give $\eta = 0.042 \pm 0.023$ and a Hubble constant near 72 km/s/Mpc, with density histories that keep dark matter and dark energy at the same order of magnitude over late and future times. If this result holds, an interacting dark-energy model would offer a phenomenological alternative to $\Lambda$CDM that avoids the fine-tuning of equal densities today.

What carries the argument

The load-bearing object is the interaction term $Q$ inserted into the fluid continuity equations. Truncating a Taylor expansion of $Q(H\rho_m, H\rho_{\mathrm{de}})$ to first order gives two one-parameter families, $Q = 3\eta H\rho_{\mathrm{de}}$ and $Q = 3\beta H\rho_m$, which yield exact analytic solutions for the fractional densities $\Omega_m(z)$ and $\Omega_{\mathrm{de}}(z)$. The argument then turns on Markov chain Monte Carlo fits of the two models to the Pantheon and OHD data, and on comparing the predicted late-time density evolutions, where the $\eta$ model's curves stay close while the $\beta$ model's do not.

What would settle it

Measure the ratio $\rho_{\mathrm{de}}/\rho_m$ over redshifts $0 \lesssim z \lesssim 2$ using tomographic weak lensing and baryon acoustic oscillations. The $\eta$ model with $\eta \approx 0.042$ predicts this ratio levels off and stays of order one at late times, whereas $\Lambda$CDM predicts it grows as $(1+z)^{-3}$; observing the $\Lambda$CDM-like growth would falsify the paper's central claim.

Watch

Extended reading notes

Core claim

The paper's central claim is that the $\eta$ model---the interaction $Q = 3\eta H\rho_{\mathrm{de}}$ with a small positive coupling---effectively addresses the coincidence problem. With the Pantheon and OHD data combined, the best fit is $\eta = 0.042 \pm 0.023$, $H_0 = 72.024 \pm 0.020$ km/s/Mpc, and $\Omega_m = 0.276 \pm 0.012$; with this coupling, the dark-energy and dark-matter densities evolve so that they remain of the same order of magnitude instead of dark energy taking over completely. The paper also finds that Pantheon data alone prefer the opposite sign of the interaction (matter flowing into dark energy), that the $\beta$ model fails to produce the clear matter-dominated epoch needed for structure formation, and that information criteria do not decisively favour any model on the combined data.

Load-bearing premise

Everything rests on treating the interaction as exactly proportional to $H$ times one of the two energy densities with a constant coupling, which is only the first term of a Taylor expansion; if higher-order or time-dependent terms matter, the fitted couplings and the density evolutions built from them are biased.

Editorial extensions

If this is right

  • The fitted $\eta$ model predicts that dark matter and dark energy densities stay within the same order of magnitude into the future, removing the need to explain why they are comparable precisely today.
  • With the combined data, the model pins $H_0$ to about 72.0 km/s/Mpc with a formal uncertainty near 0.02, sitting above the CMB-inferred value and in line with local distance-ladder measurements.
  • The $\beta$ model, with the interaction proportional to matter density, is disfavoured because its best-fit density histories lack a clear matter-dominated era, the epoch required for galaxies and large-scale structure to form.
  • The sign of the interaction depends on the data set: Pantheon alone favours matter-to-dark-energy flow, while OHD and the combined data favour the opposite, so the direction of energy transfer is not settled.
  • The information criteria AIC and BIC do not uniformly prefer any model on the combined data, so the interacting models are competitive with $\Lambda$CDM rather than decisively better.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the small positive coupling is real, it should leave a measurable signature in the growth of cosmic structure: the growth rate $f\sigma_8$ at $z \lesssim 1$ would be shifted relative to $\Lambda$CDM, a test the paper motivates but does not perform.
  • The sign flip between Pantheon and OHD suggests the fitted coupling may be absorbing systematics in distance or age calibrations rather than a physical interaction; cross-calibrating the two data sets with independent anchors could reveal which.
  • Because the interaction is derived from a first-order Taylor expansion, these data cannot yet distinguish a constant coupling from one that slowly varies; extending Hubble-parameter measurements to higher redshift would test whether the linear form holds.
  • If independent probes confirm $H_0 \approx 72$ with sub-percent precision, an interacting dark-energy sector with a positive coupling becomes a concrete alternative in the Hubble-tension debate, though the paper does not claim to solve that tension.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 9 minor

Summary. The paper studies two linear interacting dark energy models, labelled β and η, in which the interaction term is respectively Q = 3βHρ_m and Q = 3ηHρ_de. The authors perform MCMC fits to the Pantheon supernova sample (1048 points) and to a combined OHD dataset (26 BAO + 31 cosmic-chronometer measurements), obtaining parameter constraints for Ω_m, the coupling, and H0. They then use the best-fit parameters to integrate the energy-density evolution and discuss the coincidence problem, arguing that the η model with a small positive coupling 'effectively addressed' the coincidence problem by keeping the two densities comparable in magnitude at late times. The paper also reports AIC/BIC model comparisons and notes tensions between the SNIa-only and OHD-only results, including negative couplings that yield negative energy densities.

Significance. If the coincidence-problem claim were correct, the η model would be an interesting minimal alternative to ΛCDM: a single constant coupling would keep dark matter and dark energy densities in a fixed, order-unity ratio asymptotically, removing the fine-tuning of the present density ratio. The paper does provide useful, reproducible-looking MCMC constraints and writes down explicit analytic solutions, which is a strength because the central claims can be checked directly from the equations. However, that check fails: the asymptotic ratio is η/(1−η), which for the reported best fit is about 0.044, not of order unity. The plotted 'densities' are not proper fractional densities, and the abstract contains an error in the quoted Hubble-constant uncertainties. The central phenomenological conclusion is therefore not supported by the paper's own equations.

major comments (4)
  1. [§3.3 and Table 1] The claim that the η model keeps dark matter and dark energy densities 'of the same magnitude in future times' is contradicted by the analytic solutions in Table 1. With w_de = −1, the η-model solutions give Ω_m/Ω_de → η/(1−η) as z → −1 for 0 < η < 1. Using the combined OHD+SNIa fit η = 0.042, this asymptotic ratio is approximately 0.044, meaning dark energy dominates by a factor of about 23. The model does not render the two densities comparable; it only replaces the ΛCDM limit of zero with a small nonzero constant. A ratio of 0.044 is not a resolution of the coincidence problem, and the conclusion in §3.3 and in the final paragraph of §4 is therefore unsupported.
  2. [Table 1 and Figure 2] The quantities labelled Ω_m and Ω_de in Table 1 are not the usual fractional energy densities, because they do not sum to 1 at epochs other than today. For the combined-fit parameters (η = 0.042, Ω_m0 = 0.276, Ω_de0 = 0.724), evaluating the Table 1 expressions at z = 1 gives Ω_m ≈ 1.99 and Ω_de ≈ 0.79, with a sum of about 2.78. These are densities scaled by the present critical density, not by the instantaneous critical density. Consequently, Figure 2, presented as the 'scaled evolution of densities', is not a valid coincidence-problem diagnostic. The authors should plot ρ_m/ρ_de directly, or normalize the densities by the instantaneous critical density, and then reassess their qualitative conclusions.
  3. [Abstract and Table 2] The abstract quotes Pantheon-only H0 uncertainties of ±0.003 and ±0.004 km/s/Mpc, while Table 2 lists 72.158 ± 0.275 and 72.371 ± 0.355 km/s/Mpc for the same fits. The abstract understates the errors by two orders of magnitude. As written, this is an internal inconsistency in a headline number and must be corrected before the paper can be considered further.
  4. [§3.2 and Table 2] The SNIa-only best fits give negative couplings (β = −0.203 ± 0.131, η = −0.265 ± 0.144) that, as the authors acknowledge in §3.3, produce unphysical negative energy densities. Presenting these fits in Table 2 alongside the physically viable OHD and combined results, without marking them as unphysical, obscures the analysis. The paper should either restrict the parameter space to the region where the energy densities remain positive for all relevant redshifts, or explicitly exclude the SNIa-only fits from the physical discussion and from the claimed 'tension' between datasets.
minor comments (9)
  1. [§3.1] The likelihood and MCMC setup are not fully specified: the authors do not state the priors on the parameters, the number of walkers and steps, the burn-in length, or the convergence criteria. This information is needed for reproducibility.
  2. [§3.2] The sentence 'the width of the Gaussian distribution for the Hubble constant includes values from ±60 to ±70' is unclear; presumably it means the 1σ allowed range spans values from about 60 to 70 km/s/Mpc.
  3. [Table 1] The equations in Table 1 are poorly typeset, especially the β-model expression for Ω_DE, which appears garbled. Please use clear notation and define all quantities (e.g., z, w, Ω_i,0) in the caption.
  4. [§3.2] Reference [17] is a preprint, and the comparison with Planck 2018 is made only qualitatively; please provide the quoted H0 value and the error bar so the reader can judge the claimed consistency.
  5. [Abstract] The phrase 'with one having a large error margin' is vague; specify which dataset or fit has the large uncertainty.
  6. [§2] The footnote marker after 'variation of parameter method' appears misplaced, and the plural 'parameters' is the standard terminology.
  7. [§3.2] The name 'Plank' should be 'Planck'.
  8. [Figure 2 caption] The caption says 'a = 0, represent the modern universe as seen today', which is confusing; the present epoch corresponds to scale factor a = 1 (or redshift z = 0).
  9. [General] No mention is made of code or data availability; adding a reproducibility statement would strengthen the paper.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the coincidence-problem discussion is a posterior model consequence, and any issue with it is a correctness matter, not a circularity.

full rationale

The paper's derivation chain is: (i) assume a linear interaction Q = 3ηHρ_de or Q = 3βHρ_m via a Taylor-expansion ansatz; (ii) solve the continuity equations analytically (Table 1); (iii) fit the parameters η, β, Ω_m, H0 to Pantheon and OHD data with MCMC; (iv) plot the resulting density evolutions and discuss the coincidence problem. The coincidence-problem claim is a derived posterior statement: the density ratio evolution is computed from the fitted parameters, not used as a fit target. This is a standard model consequence, not a fitted input renamed as a prediction; the fit to distance moduli and Hubble parameters does not by construction force the late-time density ratio to any particular value. The interaction ansatz itself is an input assumption, not justified by the conclusion, so its physical adequacy is a modeling question, not circularity. The self-citations present ([10], [14], [16]) are not load-bearing for the central claim: the interaction formalism is traced to Campo et al. [9], the emcee package to Foreman-Mackey et al. [15], and the Hubble-parameter data are external observations. The skeptic's observation that the η-model asymptotic ratio is η/(1−η) ≈ 0.044, and that the Table 1 expressions are not properly normalized fractional densities, is a serious scientific objection to the paper's conclusions, but it is a correctness/validity concern, not an instance of the derivation reducing to its own inputs. Therefore, under the provided circularity criteria, no circular step is exhibited.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central claims rest on four assumptions: the FLRW background, the linear form of Q with constant couplings, the fixed equations of state, and the inclusion of baryons in the matter sector. The free parameters beta, eta, Omega_m0, and H0 are all fitted to the data, so the density evolution shown in Figure 2 is a consequence of the fit rather than an independent prediction.

free parameters (4)
  • beta (beta model coupling) = OHD: 0.128 +/- 0.057; SNIa: -0.203 +/- 0.131; combined: 0.041 +/- 0.016
    Dimensionless coupling in Q=3*beta*H*rho_m fitted to data via MCMC; central to the derived density evolution.
  • eta (eta model coupling) = OHD: 0.243 +/- 0.151; SNIa: -0.265 +/- 0.144; combined: 0.042 +/- 0.023
    Dimensionless coupling in Q=3*eta*H*rho_de fitted to data; the coincidence-problem analysis depends on its fitted sign and magnitude.
  • Omega_m0 (current matter density parameter) = beta model OHD: 0.555 +/- 0.128; eta model OHD: 0.445 +/- 0.115; combined values near 0.28-0.30
    MCMC free parameter; determines the density evolution and the amount of dark matter in the models.
  • H0 (Hubble constant) = beta model OHD: 65.773 +/- 2.617; eta model OHD: 65.806 +/- 3.216; combined near 72.0
    Fitted to data; used in density evolution and model comparison.
assumptions (4)
  • domain assumption The universe is described by the FLRW metric with the standard continuity equations for matter, dark energy, and radiation.
    Section 2 uses the Friedmann background and equations (1)-(4) as the starting point.
  • ad hoc to paper The interaction Q is a function of H*rho_m and H*rho_de, expandable as a Taylor series, and is truncated to first order with constant coefficients.
    Section 2 introduces Q approx lambda_m*H*rho_m + lambda_de*H*rho_de and then Q=3*eta*H*rho_de or Q=3*beta*H*rho_m; no physical mechanism or scale dependence is given.
  • domain assumption Dark energy has equation of state w_de=-1 and matter has w_m=0.
    Footnote 1 in Section 2 fixes these equations of state.
  • domain assumption The interaction involves total matter (dark plus baryonic), not dark matter alone.
    Section 2 states 'the interaction shall take place between matter as a whole (dark matter and baryonic matter) and dark energy.'

how reviews work

0 comments
Cite this review

Pith. "Pith review of Interacting dark energy models." pith.science (2026). https://pith.science/paper/TIN4TO4R

@misc{pith2026241209024,
  author       = {Pith},
  title        = {Pith review of: Interacting dark energy models},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TIN4TO4R}},
  note         = {Machine review of arXiv:2412.09024}
}
read the original abstract

This work focuses on two linear interaction models between dark matter and dark energy, which are proposed as key factors in explaining cosmic history, specifically the latetime accelerating expansion of the universe. Both models are constrained using a Markov chain Monte Carlo analysis (MCMC) using different sets of observational data. The analysis was composed using the Pantheon data set, consisting of 1048 points of SNIa distance moduli measurements from the Pantheon analysis and the Observed Hubble Parameter (OHD) data set using Baryon acoustic Oscillation (BAO), consisting of 57 data points using distance and expansion rate measurement. Both models showed promising results with the OHD data (BAO), with a interaction that results in a higher dark matter content of 56% and 44%, and a Hubble parameter of 65.7+-3km/ s/Mpc and 65.8+-3km/ s/Mpc for the interaction dependent on dark matter and dark energy respectively. The pantheon data set however predicted a reverse interaction for both models which does not follow initial assumptions that were made. The pantheon data measured a dark matter content of 18% and 20% with a Hubble parameter of 72.1 +- 0.003km/ s/Mpc and 72.3 +- 0.004km /s/Mpc. The constrained results are used to revisit the coincidence problem and other problems in standard cosmology. The analysis provided a discrepancy between the different data sets with one having a large error margin.

Figures

Figures reproduced from arXiv: 2412.09024 by the authors.

Figure 1
Figure 1. (left) Best fit model constraints of the η model with SNIa, OHD and a combination of both data sets. (right) Best fit model constraints of the β model with SNIa, OHD and a combination of both data sets. From this analysis, it can be observed that the η model was constrained with regards to the values following in the table below. Model Data set Ωm Coupling H0 (km/s)/M pc AIC BIC ΛCDM SNIa 0.28 ± 0.01 NA 71.85 ± 0.22… view at source ↗
Figure 2
Figure 2. (left)Scaled evolution of densities for the η model. (right) Scaled evolution of densities for the β model. Densities follow the same colour codes, where a = 0, represent the modern universe as seen today. Beyond this is considered into the future and behind the past. Referring to the left-hand side of [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 19 canonical work pages

  1. [1]

    Bertone G, Hooper D and Silk J 2005Physics Reports 405 279–390

    work page

  2. [2]

    Riess A G and et al 1998The Astronomical Journal 116 1009–1038

    work page

  3. [3]

    Velten H, vom Marttens R F and Zimdah W 2014The European Physics Journal 74 3160

    work page

  4. [4]

    Weinberg S 1989Rev. Mod. Phys. 61(1) 1–23

    work page

  5. [5]

    Chantavat T 2023Journal of Physics: Conference Series 2431

    work page

  6. [6]

    Randall L and Sundrum R 1999Physical Review Letters 83 4690–4693

    work page

  7. [7]

    Bahamonde S, Böhmer C G and Wright M 2015Physical Review D 92

    work page

  8. [8]

    Bludman S A and Roos M 2002Physical Review D 65

    work page

Show all 20 references
  1. [9]

    Campo S d, Herrera R and Pavón D 2009Journal of Cosmology and Astroparticle Physics 2009 020–020

    2009

  2. [10]

    van der Westhuizen M and Abebe A 2024Journal of Cosmology and Astroparticle Physics 2024 048

    2024

  3. [11]

    Marsh M D 2017Physical Review Letters 118

  4. [12]

    2018 The Astrophysical Journal 859 101

    Scolnic D Met al. 2018 The Astrophysical Journal 859 101

    2018

  5. [13]

    Farooq O and Ratra B 2013The Astrophysical Journal Letters 766 L7

  6. [14]

    Akalu S, Alfedeel A and Abebe A 2024Eur. Phys. J. C

  7. [15]

    Foreman-Mackey D, Hogg D W, Lang D and Goodman J 2013PASP 125 306–312 (Preprint 1202.3665)

    arXiv

  8. [16]

    Hough R T, Abebe A and Ferreira S E S 2019arXiv: General Relativity and Quantum Cosmology

  9. [17]

    Di Valentino E, Melchiorri A, Mena O and Vagnozzi S 2019 Interacting dark energy after the latest planck, des, and h0 measurements: an excellent solution to theh0 and cosmic shear tensions

    2019

  10. [18]

    Chakrabarti A and Ghosh J K 2011Philosophy of Statistics (Handbook of the Philosophy of Science vol 7) ed Bandyopadhyay P S and Forster M R (Amsterdam: North-Holland) pp 583–605

  11. [19]

    Salcido J, Bower R G, Barnes L A, Lewis G F, Elahi P J, Theuns T, Schaller M, Crain R A and Schaye J 2018 Monthly Notices of the Royal Astronomical Society 477 3744–3759

    2018

  12. [20]

    Moresco M and et al 2012Journal of Cosmology and Astroparticle Physics 2012 006–006

    2012

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.