Viscous cosmological fluids and large-scale structure

A Abebe; BG Mbewe; RR Mekuria; S Sahlu

REVIEW 4 major objections 6 minor 19 references

Viscous cosmological fluids and large-scale structure

T0 review · 4 major / 6 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 A viscous interacting dark-fluid model with a ΛCDM-like background predicts that matter density perturbations develop singularities and growing oscillations, so bound structures disintegrate in the late universe.

desk verdict The claimed late-time disintegration of structure is likely a singularity in an ill-defined perturbation variable, not a physical prediction. read the letter →

arxiv 2412.02276 v2 pith:TENM6SG4 submitted 2024-12-03 gr-qc

BG Mbewe , RR Mekuria , S Sahlu , A Abebe This is my paper · ORCID

classification gr-qc

keywords viscousdarkfluidsinteractingenergy1+3covariantformalismcosmologicalperturbationslarge-scalestructureformationMCMCparameterestimationTypeIasupernovaedisintegration

topics Dark Energy

open problems Dark 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

The paper sets out a cosmological model in which dark matter and dark energy exchange energy through a viscous interaction term, with a dark-energy equation of state tuned so the background expansion mimics ΛCDM. It tests that background against Type Ia supernova distances and studies linear density perturbations in a dust-dominated frame. Its central claim is that the perturbation equations, solved with the fitted parameters, make the dust density contrast blow up near redshift 2 and oscillate with growing amplitude at short wavelengths; the model therefore predicts that bound structures in the late universe rip apart. The significance is that this is a concrete, testable difference from ΛCDM: the same model that fits expansion history can fail on structure formation, so large-scale structure data could discriminate between the two.

What carries the argument

The load-bearing mechanism is the linear perturbation system (Eq. 10) for the dust density contrast Δ_d, the dark-energy density contrast Δ_Λ, and the volume-expansion gradient Z, built in the 1+3 covariant formalism (splitting spacetime into a fundamental time direction and spatial hypersurfaces). The interaction Q=δHρ_d appears directly in the Δ_d and Δ_Λ evolution equations, and the dark-energy equation of state is reduced to p_Λ=(A0−1−ζ0)ρ_Λ. Solving that system numerically with MCMC-best-fit parameters gives the singularity near z≈2 and the growing oscillations that the paper reads as a disintegration of bound structures.

What would settle it

Re-derive Eq. (10) from the full 1+3 covariant equations in a gauge-invariant formulation (for example using a comoving curvature perturbation or density contrast on uniform-density slices) and check whether the singularity near z≈2 and the growing oscillations persist; if they vanish, the disintegration is an artifact of the chosen variables. A complementary check would compare the predicted late-time decay or blow-up of clustering with observed growth data (e.g. fσ8) at z≲2, where ΛCDM shows continued growth and the VIDF model predicts disintegration.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that a viscous interacting dark-fluid (VIDF) universe, despite being constructed to reproduce the ΛCDM background under certain parameter choices, does not reproduce ΛCDM at the level of structure formation. In a dust-matter-dominated frame, the linear dust density contrast Δ_d develops a singularity near z≈2 for long wavelengths and growing oscillations for short wavelengths, depending on the MCMC parameter values. The paper interprets this as the disintegration — the 'rip' — of large-scale bound structures at late times, and attributes it mainly to the dark-fluid interaction term Q=δHρ_d.

Load-bearing premise

The load-bearing premise is that the linearized perturbation equations (Eq. 10), whose derivation is not included in the paper and is instead left to an external link in a footnote, are the correct and complete evolution equations for the VIDF model; if those equations contain sign errors, gauge artifacts, or unjustified simplifications, the predicted late-time disintegration would not follow.

Editorial extensions

If this is right

If the VIDF model is correct, bound structures such as galaxy clusters and filaments would not persist into the late universe; the matter distribution would be torn apart by growing oscillations and singularities in the density contrast.
The interaction between dark energy and dark matter is identified as the cause of this behavior, so the model provides a signature to distinguish interacting dark-sector models from ΛCDM using large-scale structure rather than only the expansion history.
The MCMC-fit parameters that agree with supernova data still produce the disintegration, meaning a good background fit does not guarantee viable structure formation.
Because the radiation-dust equality and dust-dark-energy equality occur at different redshifts than in ΛCDM, the model also predicts detectable shifts in the epochs of matter-radiation equality and acceleration.
The paper's own conclusion states that more data sets (BAO, CMB, OHD, R22) would be needed to test how well the model fits observations.

Reading between the lines

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

Editorial inference: the singularity near z≈2 may be a gauge artifact of the chosen scalar density variables; a gauge-invariant perturbation analysis would be needed to confirm that the disintegration is physical rather than an artifact of the slicing.
Editorial inference: the model's negative dark-energy density at some epochs violates energy conditions; this may be the underlying driver of the instability, and a version that imposes energy conditions could be tested to see whether the rip disappears.
Editorial inference: growth-rate data such as fσ8 at low redshift would provide a sharper test than the supernova distances used here, since the model predicts decaying or singular clustering rather than continued growth.

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.

Desk Editor's Note

read the letter

The paper is a proceedings-style cosmology paper that builds a viscous interacting dark-energy (VIDF) model, fits it to SNIa with MCMC, and claims the model predicts late-time disintegration of large-scale structure. The one genuinely new item is the linear perturbation system in Eq. (10) for this specific model and the resulting plots. That claim, though, does not survive contact with the paper's own equations.

The background treatment is standard but has typos and sign ambiguities, and the MCMC setup is underdescribed (no likelihood, no best-fit table). The perturbation derivation is not in the paper; it is relegated to a Google Drive folder. More importantly, the singularity in the dust density contrast at z≈2 is very likely an artifact of the perturbation variables. Their Eq. (5) explicitly allows ΩΛ to become negative, and they admit the model "violates all the energy conditions by having a dark energy density that is negative." The variable Δ_Λ is defined with 1/ρ_Λ in Eq. (9), and Eq. (10) has denominators ΩΛ. When ΩΛ crosses zero, Δ_Λ and the coupled system blow up. So the "rip" is a zero-crossing in the chosen variable, not a demonstrated physical instability. The paper offers no invariant criterion, no comparison to ΛCDM growth, and no check that the singularity is coordinate-independent. That's a load-bearing flaw, not a minor one.

Credit where due: the background equations are at least written out, and the MCMC plots are there, even if the parameters are not tabulated. The authors honestly call the results preliminary. But the central claim rests on an omitted derivation and a likely variable singularity.

I would not send this to a serious referee. It needs the derivation written out, the variables regularized or an invariant diagnostic used, and a genuine comparison with ΛCDM before it can support the disintegration claim. If the authors fix those, the model question might be worth revisiting; as it stands, the paper is a good example of why perturbation variables need to be checked for singularities.

Recommendation: desk reject, with the option to resubmit if the perturbation analysis is redone.

Referee Report

4 major / 6 minor

Summary. The manuscript studies a viscous interacting dark-fluid (VIDF) cosmological model with an inhomogeneous dark-energy equation of state, fits its background parameters to Type Ia supernova data via MCMC, and then uses a 1+3 covariant perturbation formalism to argue that linear density perturbations exhibit singularities and growing oscillations at late times, leading to the claim that the model predicts the disintegration of bound large-scale structures. The background solutions are compared with ΛCDM, and the MCMC results are used to evaluate the perturbation equations at the best-fit and boundary parameter values.

Significance. If established, the central claim would be significant: it would indicate that the VIDF model is ruled out by structure-formation considerations and would illustrate how interacting viscous dark fluids can produce unphysical late-time behavior. The paper also provides a concrete background-plus-perturbation framework that could be checked against other data sets. However, the presentation is not self-contained: the perturbation equations are not derived in the text, the MCMC setup is not described, and the interpretation of the perturbation singularities is not checked against variable artifacts. These issues currently prevent the central conclusion from being evaluated.

major comments (4)

[Section 4, footnote 1] The perturbation system in Eq. (10) is the sole basis for the paper's central claim, yet its derivation is not present in the manuscript; footnote 1 refers to an unversioned Google Drive folder. Because the signs, gauge choices, and approximations in these equations cannot be verified from the paper, and because an error in any of them would directly change the predicted structure growth, this omission is load-bearing and not acceptable for a self-contained journal submission.
[Section 4, Eqs. (9) and (10), Fig. 6] The reported singularity at z≈2 is structurally tied to the perturbation variables becoming singular. The paper states in Section 2 that the VIDF model has a negative dark-energy density (Fig. 1), so ΩΛ crosses zero; the variable ΔΛ = (a/ρΛ)∇_a ρΛ and the coefficients δΩ_d/Ω_Λ in Eq. (10) are singular at that crossing. The authors interpret the resulting divergence and oscillation in Δ_d as 'disintegration of bound structures,' but this is likely a coordinate or variable artifact rather than a physical instability of matter clustering. No gauge-invariant or regular variable is used to confirm that the effect survives a change of perturbation variables.
[Section 3] The MCMC results are not reproducible from the text. The paper does not specify the likelihood function, the SNIa sample, the number of data points, the priors, or convergence diagnostics, and it delegates the methodology to Ref. [13]. In addition, radiation is fixed by hand because the simulation fails to constrain it; the impact of this ad hoc treatment on the fitted background parameters and on the subsequent perturbation study is not assessed.
[Section 4] The perturbations are said to be studied in a 'dust-matter-dominated frame,' while the background includes radiation (Eq. (1)). The paper gives no justification for dropping radiation in the perturbation equations while keeping it in the background; this could modify the perturbation evolution at precisely the redshifts where the claimed singularity appears.

minor comments (6)

[Section 2, Eq. (7)] There is a typographical error in Eq. (7): '3(A0 − ζ0)1 + z)3(A0−ζ0)' should presumably read '3(A0 − ζ0)(1 + z)3(A0−ζ0)'.
[Section 2, Eq. (5)] The expression for ΩΛ in Eq. (5) has ambiguous bracket structures from the typesetting, making the intended formula difficult to parse; please rewrite it with clear notation.
[Fig. 6 caption] The caption of Fig. 6 refers to the 'VIDE model,' which is likely a typo for 'VIDF model'; please make the naming consistent throughout.
[Abstract and Section 5] The abstract uses 'disintegration' while Section 5 uses 'rip'; please define these terms operationally and use a consistent phrase.
[Section 2] The paper should either state the parameter ranges for which ΩΛ remains positive or explicitly discuss whether a negative dark-energy density is physically acceptable; currently the violation of all energy conditions is mentioned only in passing.
[Footnote 1] The phrase 'proceedings paper' in footnote 1 suggests this is a conference contribution; if the manuscript is intended for a journal, it should be made fully self-contained rather than relying on an external folder.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the predicted late-time disintegration is an output of the perturbation system evaluated at MCMC-fitted parameters, not a quantity used in the fit or defined in terms of the conclusion.

full rationale

The paper's derivation chain is linear and non-circular. Background conservation equations (1) with the chosen interaction Q=δHρ_d and reduced dark-energy equation of state (4) yield the analytic solutions (5)-(6); MCMC fits the free parameters to SNIa distance moduli via Eq. (8); the perturbation system (10) is then integrated at those fitted central values to produce the density-contrast plots and the late-time 'rip' conclusion. Nothing in the text indicates that the perturbation result fed back into the parameter fit, and no perturbation variable is defined in terms of the claimed disintegration. The heavy reliance on the authors' prior works ([13], [14], [17]) and the unversioned Google Drive link in footnote 1 are verifiability limitations, not circular reductions, because those sources supply methodology and formalism rather than the target conclusion. A separate, non-circular correctness concern should be noted: the paper itself states (Section 2, discussion of Fig. 1) that the model has negative dark-energy density, so ΩΛ crosses zero; the variable ΔΛ ≡ (a/ρΛ)∇_aρΛ in Eq. (9) and the explicit ΩΛ denominators in Eq. (10) are singular there. The claimed singularity near z≈2 may therefore be a variable/background artifact rather than a physically independent prediction, but this is an internal-consistency issue, not a case of the result being equivalent to its inputs by construction.

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

The model is a phenomenological extension of ΛCDM: it adds an interaction Q=δHρ_d and a simplified inhomogeneous dark-energy equation of state, relying on a standard FLRW background and the 1+3 covariant perturbation formalism. No new particle, force, or dimension is introduced. The free parameters are the MCMC-fitted interaction, equation-of-state, density, and Hubble parameters; Ωr0 is fixed by hand. A non-standard load-bearing assumption is that the full derivations are valid even though they are available only through an external Google Drive link.

free parameters (6)

δ (interaction coupling) = not reported in text; MCMC posterior shown in Fig. 5
Dimensionless coupling in Q=δHρ_d (Eq. 2); enters background and perturbation equations; constrained by SNIa MCMC but central value not tabulated.
A0 (EoS coefficient) = not reported in text
Appears in dark-energy pressure p_Λ=(A0-1-ζ0)ρ_Λ (Eq. 4); fitted by MCMC; the ΛCDM limit requires A0=ζ0.
ζ0 (viscosity parameter) = not reported in text
Coefficient of the viscous term in p_Λ; fitted by MCMC; appears in all background and perturbation equations.
Ωd0 (present dust density) = not reported; described as fairly good agreement with Planck 2018
Present fractional dust density; constrained by SNIa MCMC; used in the Hubble parameter and density-contrast equations.
H0 (Hubble constant) = not reported; described as outside 1σ of Planck 2018
Present expansion rate; constrained by SNIa distance modulus data; needed for the perturbation equations.
Ωr0 (radiation density) = fixed constant; value not stated
Radiation fraction is treated as constant because MCMC cannot resolve it (Section 3); it enters the background energy densities.

assumptions (6)

domain assumption The universe is a flat FLRW spacetime with 8πG=c=1 units.
Eq. (1) sets the background geometry and units for all subsequent equations.
domain assumption The dark-sector interaction is Q=δHρ_d with constant δ.
Taken from Brevik [16] (Eq. 2); a phenomenological choice not derived from a microphysical model.
domain assumption The dark-energy equation of state has the inhomogeneous form p_Λ=A0ρ_Λ^α - ρ_Λ - ζ0ρ_Λ0(ρ_Λ/ρ_Λ0)^m, with α=m=1.
Equation (3) is a model ansatz; assuming α=m=1 reduces it to Eq. (4).
domain assumption Linear perturbations can be analyzed in a dust-matter-dominated frame using the 1+3 covariant formalism, with radiation neglected in the perturbations.
Section 4 states the frame; the perturbation equations (10) do not include radiation variables.
domain assumption The MCMC likelihood for Type Ia supernovae follows the methodology of ref. [13] and is applicable to this model.
Section 3 does not present the likelihood or dataset; it refers to earlier work by the same group.
ad hoc to paper An unversioned Google Drive folder contains the complete derivations of Eqs. (5)-(10).
Footnote 1 sends readers to an external folder; the derivations are not in the manuscript.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Viscous cosmological fluids and large-scale structure." pith.science (2026). https://pith.science/paper/TENM6SG4

@misc{pith2026241202276,
  author       = {Pith},
  title        = {Pith review of: Viscous cosmological fluids and large-scale structure},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TENM6SG4}},
  note         = {Machine review of arXiv:2412.02276}
}

abstract

In this paper, we study the viscous fluid cosmological model that when certain conditions are invoked mimics the $\Lambda$CDM model. The background equations governing the evolution of viscous interacting fluids in a multifluid system are derived. The Markov Chain Monte Carlo (MCMC) simulation is applied to constrain the best-fit cosmological parameters with Supernova Type 1a data. In addition, linear cosmological perturbations are investigated in a dust-matter-dominated frame using a $1+3$ covariant formalism approach. It is evident from the perturbation results obtained that the model predicts the disintegration of bound structures of large-scale structures in the late-time universe.

Figures

Figures reproduced from arXiv: 2412.02276 by the authors.

**Figure 1.** The evolution of energy density parameter for a multi-fluid universe. H =H0 Ωr0 (1 + z) 4 − (1 + z) 3(A0−ζ0) + Ωd0 3(1 + ζ0 − A0) 3(1 + ζ0 − A0) − δ (1 + z) 3−δ − (1 + z) 3(A0−ζ0) + (1 + z) 3(A0−ζ0) 1 2 . (6) The deceleration parameter can be obtained from Eq. (6) and it yields as: q(z) = 1 2h 2 Ωr0 4(1 + z) 4 − 3(A0 − ζ0)1 + z) 3(A0−ζ0) + Ωd0 3(1 + ζ0 − A0) 3(1 + ζ0 − A0) − δ (3 − δ)(1 + z) 3… view at source ↗

**Figure 2.** The evolution of Hubble parameter [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 5.** Best fit parameters of VIDF model using MCMC. In [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗

Figures from the paper (1 more)

**Figure 6.** Figure 6: The density contrast of VIDE model as a function of redshift. noticed when considering the lower boundary values obtained from the simulation. However, in Fig. 6b, there exist no singularity which was experienced by the model in Figs. 6a and 6c and we thus see the grow…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 18 canonical work pages

[13]

Hough R, Sahlu S, Sami H, Elmardi M, Swart A M and Abebe A 2021 arXiv preprint arXiv:2112.11695

work page Pith review arXiv 2021
[1]

Bamba K, Capozziello S, Nojiri S and Odintsov S D 2012 Astrophysics and Space Science342 155–228

work page 2012
[2]

Copeland E J, Sami M and Tsujikawa S 2006 International Journal of Modern Physics D15 1753–1935

work page 2006
[3]

Shi K, Huang Y and Lu T 2012 Monthly Notices of the Royal Astronomical Society426 2452–2462

work page 2012
[4]

Fabris J C, Goncalves S and Ribeiro R d S 2006 General Relativity and Gravitation38 495–506

work page 2006
[5]

Bolotin Y L, Kostenko A, Lemets O A and Yerokhin D A 2015 International Journal of Modern Physics D 24 1530007

work page 2015
[6]

Chimento L P, Jakubi A S and Pav´ on D 2000 Physical Review D62 063508

work page 2000
[7]

Astorga-Moreno J, Chagoya J, Flores-Urbina J and Garc ´ ıa-Aspeitia M A 2019 Journal of Cosmology and Astroparticle Physics 2019 005

work page 2019

Show all 19 references

[8]

Garcia-Aspeitia M A, Hernandez-Almada A, Magana J, Amante M H, Motta V and Mart ´ ınez-Robles C 2018 Physical Review D97 101301

2018
[9]

Garc ´ ıa-Aspeitia M A, Mart ´ ınez-Robles C, Hern´ andez-Almada A, Maga˜ na J and Motta V 2019Physical Review D 99 123525
[10]

Hernandez-Almada A, Maga˜ na J, Garc ´ ıa-Aspeitia M A and Motta V 2019The European Physical Journal C 79 1–9
[11]

Abebe A, Abdelwahab M, De la Cruz-Dombriz A and Dunsby P K 2012 Classical and quantum gravity29 135011

2012
[12]

Ellis G F, Maartens R and MacCallum M A 2012 Relativistic cosmology(Cambridge University Press)

2012
[14]

Sahlu S, Ntahompagaze J, Elmardi M and Abebe A 2019 The European Physical Journal C79 1–31

2019
[15]

Brevik I, Obukhov V and Timoshkin A 2015 Astrophysics and Space Science355 399–403

2015
[16]

Brevik I, Myrzakulov R, Nojiri S and Odintsov S D 2012 Physical Review D 86 ISSN 1550-2368 URL http://dx.doi.org/10.1103/PhysRevD.86.063007

2012 doi
[17]

van der Westhuizen M A and Abebe A 2023 Interacting dark energy: clarifying the cosmological implications and viability conditions ( Preprint 2302.11949)

2023 arXiv
[18]

Chanda A, Roy B C, Bamba K and Paul B C 2023 arXiv preprint arXiv:2309.09158

2023 arXiv
[19]

Benetti M, Borges H, Pigozzo C, Carneiro S and Alcaniz J 2021 Journal of Cosmology and Astroparticle Physics 2021 014

2021

Pith tools

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

[1] [13]

Hough R, Sahlu S, Sami H, Elmardi M, Swart A M and Abebe A 2021 arXiv preprint arXiv:2112.11695

work page Pith review arXiv 2021

[2] [1]

Bamba K, Capozziello S, Nojiri S and Odintsov S D 2012 Astrophysics and Space Science342 155–228

work page 2012

[3] [2]

Copeland E J, Sami M and Tsujikawa S 2006 International Journal of Modern Physics D15 1753–1935

work page 2006

[4] [3]

Shi K, Huang Y and Lu T 2012 Monthly Notices of the Royal Astronomical Society426 2452–2462

work page 2012

[5] [4]

Fabris J C, Goncalves S and Ribeiro R d S 2006 General Relativity and Gravitation38 495–506

work page 2006

[6] [5]

Bolotin Y L, Kostenko A, Lemets O A and Yerokhin D A 2015 International Journal of Modern Physics D 24 1530007

work page 2015

[7] [6]

Chimento L P, Jakubi A S and Pav´ on D 2000 Physical Review D62 063508

work page 2000

[8] [7]

Astorga-Moreno J, Chagoya J, Flores-Urbina J and Garc ´ ıa-Aspeitia M A 2019 Journal of Cosmology and Astroparticle Physics 2019 005

work page 2019