{"id":"821135a1-eb24-4a46-8f76-1e38e3edac13","arxiv_id":"2411.13379","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A six-parameter viscous interacting dark energy model is fitted to Union 2.1 supernovae, yielding h0=0.68, omega_de=-1.07, eta=0.05, but the underlying evolution equations for the interacting case are not presented.","lead":"This paper fits a model of dark energy with both bulk viscosity and dark matter interactions to supernova data, reporting best-fit values for six parameters. It presents derived expansion diagnostics, but the key equations for the interacting model are not shown, so the fit cannot be reproduced from the text.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (16), the only Hubble law entering the Union 2.1 likelihood, is not closed and is not the solution of the stated interacting-viscous system; the quoted MCMC parameters are therefore not attached to a well-defined model.","rationale":"The paper's central claim requires that the H(z) used in the Union 2.1 distance-modulus likelihood be the H(z) of the model defined by Eqs. (1)-(3) and (11). This condition fails at Eq. (16). The noninteracting VDE Hubble laws in Eqs. (12)-(15) are self-consistent and derivable, so the problem is not viscosity itself; it is the step to the interacting case. There, the densities rho_chi and rho_de are left undefined, and the functional form of Eq. (16) is not obtainable from the coupled linear system: the exact densities are superpositions of two eigen-solutions, and the DE redshift exponent is missing the (1+omega_de) factor present in Eq. (15). Since the fit returns omega_de=-1.07, this is not a negligible special-case coincidence. A full MCMC with an unspecified or incorrect H(z) therefore does not support the quoted parameters. I agree with the reader's weakest-assumption diagnosis; my check sharpens it by showing that even supplying evolution equations would not recover Eq. (16) as written. Additional weaknesses (no Lambda-CDM baseline comparison, dated Union 2.1 data only) are secondary; the Eq. (16) problem is sufficient to reject the central numerical claim. If the authors supply the correct derivation and rerun the fit, a CONDITIONAL resubmission would be appropriate, but as submitted the claim is not reproducible.","tokens_in":13360,"tokens_out":16413,"duration_ms":160499,"concrete_test":"Derive the exact solution of Eqs. (1)-(2) with Q=3H(lambda_de rho_de+lambda_chi rho_chi) and p_de=omega_de rho_de-3 eta H^2, including baryons and curvature, and compare it with Eq. (16). Concretely, use the quoted best fit h0=0.68, lambda_r=0.09, lambda_theta=0.41, omega_de=-1.07, Omega_k0=0.10, eta=0.05; numerically integrate the coupled ODEs from z=0 to z=1.5 with initial conditions satisfying E(0)=1; compute distance moduli from the resulting H(z). Re-run the Union 2.1 chi^2 at the quoted parameters and, if feasible, repeat the MCMC with this correct H(z). If the chi^2 changes by more than the typical data scatter or the best-fit moves outside the quoted 1-sigma intervals, the paper's numerical claim fails. Also check the redshift exponent by replacing (1+z)^{-3 eta} in Eq. (16) with (1+z)^{3(1+omega_de-eta)}.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is a parameter estimate from a likelihood, but the H(z) entering that likelihood is not defined by the paper. Equation (16) contains rho_chi and rho_de as functions of z with no evolution equations supplied; if they are meant as present densities, the dark-energy term is independent of omega_de, and if they are functions, the system is unclosed. Solving the model's own equations (1), (2), (3), (11) in a flat universe (ignoring baryons and curvature for the argument) gives the linear system d r_chi/dx = -3(1-lambda_chi) r_chi + 3 lambda_de r_de, d r_de/dx = -3(lambda_chi-eta) r_chi - 3(1+omega_de+lambda_de-eta) r_de, with x=ln a and r=rho/(3H0^2). The exact solution is a superposition of the eigenvalues of this 2x2 matrix, not the factorized form of Eq. (16), in which dark matter is bundled with baryons in a single w/(w-eta) term and dark energy is multiplied by (1+z)^{-3 eta}. Moreover, the dark-energy scaling in Eq. (16), (1+z)^{-3 eta}, disagrees with the noninteracting limit Eq. (15), whose last factor is (1+z)^{3(1+omega_de-eta)}; the two match only for special parameter values, not at the quoted best fit omega_de=-1.07. Thus the expansion history used for the distance moduli is not the expansion history of the model being claimed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a generalized interacting dark energy (IDE) model with a two-parameter coupling Q = 3H(λdeρde + λχρχ), parameterized by λχ = λr sin λθ and λde = λr cos λθ, in a bulk-viscous dark energy framework with ζ = ηH. It presents three viscous dark energy Hubble expressions, a generalized form, and then an 'interacting' Hubble expression (Eq. 16). Using that expression, the authors compute deceleration, jerk, snap, Om and Om3 diagnostics, and perform an MCMC fit to the Union 2.1 supernova dataset, reporting optimum parameters h0 = 0.68, λr = 0.09, λθ = 0.41, ωde = −1.07, Ωk0 = 0.10, η = 0.05. The central quantitative claim is that this interacting viscous model fits the data and yields mildly phantom dark energy with small dark-sector coupling.","tokens_in":13724,"tokens_out":8225,"duration_ms":86934,"significance":"If the model and fit were correctly derived and reproducible, the paper would offer a useful two-parameter unification of common IDE interaction forms and a first constraint on their combination with bulk viscosity. That would be of interest to the dark-energy phenomenology community. However, the central Hubble expression used in the likelihood is never derived or closed, the viscous-pressure definition is internally inconsistent, and the MCMC analysis reports no goodness-of-fit or model comparison. As presented, the paper does not deliver reliable constraints on the interaction or viscosity parameters. It also provides no code, data products, or machine-checked derivations, so the quantitative claims cannot be independently verified.","major_comments":[{"comment":"Equation (16), which is the only Hubble law entering the Union 2.1 likelihood, is not closed and does not follow from the stated model. The quantities ρχ and ρde appear without evolution equations or a specification of whether they are present-day densities or redshift-dependent functions; as written, the expression cannot be differentiated or integrated to distance moduli. Solving the continuity equations (1)-(2) with the interaction (3) and the viscous pressure (11) in a flat two-fluid system yields a coupled linear system whose solution is a superposition of two power laws, not the factored form of Eq. (16). Furthermore, Eq. (16) does not reduce to the noninteracting limit Eq. (15): the dark-energy term in Eq. (16) scales as (1+z)^{-3η} and lacks the factor (1+z)^{3(1+ωde)} present in Eq. (15), so the two expressions agree only for special parameter values, not for the quoted best fit ωde = −1.07. The reported parameters h0 = 0.68, λr = 0.09, λθ = 0.41, ωde = −1.07, Ωk0 = 0.10, η = 0.05 are therefore not attached to a well-defined model, and the fit cannot be reproduced or checked.","section":"III, Eq. (16); IV, MCMC fit"},{"comment":"The definition of the viscous pressure is internally inconsistent. Equation (9) states pde = ωdeρde = ζθ, while Eq. (11) gives pde = ωdeρde − 3ηH²; with ζ = ηH and θ = 3H these two expressions differ by the sign of the viscous term. The viscous contribution in Eq. (11) is the physically expected −ζθ, so Eq. (9) is not merely a notation issue but a contradictory definition of pde. Since η enters all subsequent Hubble expressions (12)-(16) through this pressure, the derivation of the three VDE models rests on an ill-defined equation.","section":"III, Eqs. (9) and (11)"},{"comment":"The MCMC fit is underspecified and its output is not assessed. Equation (16) contains Ωb0, which is not listed in Table III or given anywhere, and the value of Ωm0 used in the likelihood is never reported; without these values the distance moduli cannot be reproduced. In addition, the paper reports only the posterior peak and contour plots. It gives no χ², no log-evidence, no comparison with ΛCDM or with the noninteracting VDE models (12)-(14), and no goodness-of-fit statistic. The conclusion that the generalized interacting model is favored—or even that the fit is acceptable—is therefore unsupported by the presented evidence.","section":"IV, MCMC paragraph and Table III"},{"comment":"The three-point diagnostic is misdefined. With Om(z2,z1) = [h²(z2) − h²(z1)]/[(1+z2)³ − (1+z1)³] and Om(z3,z1) = [h²(z3) − h²(z1)]/[(1+z3)³ − (1+z1)³], the ratio Om(z2,z1)/Om(z3,z1) has h²(z3) − h²(z1) in the denominator, not h²(z3) − h²(z2) as written in Eq. (23). The plotted Om3 curves and the phantom/quintessence interpretation in Fig. 3(b) are therefore based on an incorrect formula.","section":"IV, Eq. (23)"}],"minor_comments":[{"comment":"The text says the three benchmark models are considered to investigate the effect of DM-DE interaction on brightness temperature, but the paper contains no brightness-temperature analysis; this appears to be a leftover from a different study.","section":"II, first paragraph"},{"comment":"The sentence 'the black and green solid lines represent the evolution of dark matter density parameter Ωde' is confusing: the notation Ωde is the dark energy density parameter, not dark matter, and the accompanying text should distinguish Ωde, Ωχ, Ωb, and Ωrad clearly.","section":"IV, Fig. 2(a) discussion"},{"comment":"The text says the topmost plots in the first, second, third, and fourth columns of Figure 7 represent six parameters (h0, λr, λθ, ωde, Ωk0, η) respectively, which is arithmetically inconsistent; the posterior peak for λθ is also quoted as both λθ = 0.41 and λθ ≈ 7π/10, and these two values are not compatible.","section":"IV, MCMC paragraph"},{"comment":"There is a stray '.' line immediately after the definition of the snap parameter, which should be removed.","section":"IV, after Eq. (19)"},{"comment":"The figures display Mathematica-style output artifacts (e.g., 'Out[27]=', 'Out[33]=', 'Out[39]=', 'Out[45]=' in the captions), and the λθ axis ticks are displayed in radians without a label; these should be cleaned for publication.","section":"Figures 5 and 6"}],"recommendation":"reject","confidential_remarks":"The central quantitative result is not reproducible because Eq. (16), the Hubble expression used in the likelihood, is not the solution of the model's own equations and is not closed. The internal inconsistency in the viscous-pressure definition and the incorrect Om3 formula further weaken the technical content. I see no presentation-level fixes that could rescue the central claim without redoing the model derivation and the MCMC analysis, which is why I recommend rejection rather than major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The reader's main concern is correct, and it is load-bearing. The central MCMC fit is built on Eq. (16), which still contains unsolved rho_chi and rho_de terms. The evolution equations the paper writes down for the interacting viscous system do not lead to that factorized H(z). The stress-test note has it right: solving the pair of conservation equations gives a 2x2 linear system whose solution is a superposition of eigenmodes, not the single dark-energy scaling (1+z)^{-3eta} shown in Eq. (16). There is also a direct internal inconsistency: the noninteracting limit Eq. (15) has a factor (1+z)^{3(1+omega_de-eta)} for dark energy, while Eq. (16) gives (1+z)^{-3eta}; these match only for special parameter values, not at the reported best fit. So the likelihood was evaluated with an expansion history that is not the expansion history of the claimed model. That alone breaks the paper's central quantitative claim.\n\nWhat is worth crediting: the paper does try to unify viscous dark energy and linear dark-sector interactions under a single trigonometric parametrization, and the exploratory figures convey how viscosity and coupling shift q, j, s, and the Om diagnostics. That is a reasonable phenomenological exercise. The bibliography covers the standard IDE and VDE literature, and the MCMC machinery is standard. The problem is not the ambition; it is that the model is not actually closed before the fit.\n\nOther soft spots, in roughly decreasing severity. Equation (9) conflates the barotropic pressure p = omega rho with the viscous stress zeta theta, which is not the correct Eckart relation. The diagnostics q, j, s, Om, Om3 are all recomputed from the same fitted H(z), so they are functions of the fitted parameters, not independent tests. There is no baseline comparison to Lambda-CDM or to the noninteracting viscous models, so we learn little about whether the interaction term is actually preferred. The posterior for lambda_theta is almost uniform (the paper itself says the 'optimum' is not significant), yet it is still listed as a headline parameter.\n\nWho is this for? A reader interested in the phenomenology of viscous dark energy with interactions might skim the early sections, but anyone wanting reliable constraints should not use these numbers. The paper deserves a serious referee only after the authors derive the correct H(z) for the interacting system, redo the fit, and add a Lambda-CDM comparison. In its current form I would not send it to peer review; I would desk-reject with an invitation to resubmit once the model is well-defined and reproducible.","headline":"The MCMC headline numbers rest on a Hubble law that the paper never derives and that does not follow from the model's own equations; until that is fixed, the central result is not reproducible.","tokens_in":14259,"tokens_out":3120,"would_cite":false,"duration_ms":37677,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.36.+x","95.35.+d","98.80.-k"],"model":"deepseek-v4-flash","headline":"The paper claims that a generalized viscous interacting dark energy model fits Union 2.1 supernova data and returns a mildly phantom equation of state near $\\omega_{\\rm de} \\approx -1.07$, along with small but nonzero dark-sector coupling…","keywords":["dark matter-dark energy interaction","bulk viscosity","viscous dark energy","interacting dark energy","cosmic acceleration","Union 2.1 supernovae","Markov Chain Monte Carlo","Om3 diagnostic"],"falsifier":"Recompute the distance moduli from equations (1), (2), (3), and (16) using the reported best-fit parameters and compare them with the Union 2.1 likelihood; if the resulting posteriors differ from Figure 7, the stated parameters are not the model's actual best fit. A simpler check is to solve the coupled density equations numerically with the reported parameters and verify that the Hubble parameter used in the analysis is recovered.","tokens_in":13111,"feed_emoji":"🌌","tokens_out":9457,"duration_ms":87100,"temperature":0.7,"pith_summary":"The paper tries to establish that the cosmic expansion history can be described by a dark energy fluid that is both bulk-viscous and interacting with dark matter, and that this combination is visible in supernova data. To do so it proposes a generalized interacting dark energy model whose coupling constants are packaged into a strength $\\lambda_r$ and a mixing angle $\\lambda_\\theta$, and it takes the dark energy viscosity to be proportional to the Hubble rate. Fitting the resulting Hubble parameter to the Union 2.1 Type Ia supernova sample with Markov Chain Monte Carlo, the paper reports best-fit values $h_0 = 0.68$, $\\lambda_r = 0.09$, $\\lambda_\\theta = 0.41$, $\\omega_{\\rm de} = -1.07$, $\\Omega_{k,0} = 0.10$, and $\\eta = 0.05$. A sympathetic reader would take this as evidence that the data mildly prefer a phantom-like dark energy with a small viscous coefficient and a nonzero dark-sector coupling. The paper also shows that viscosity and the interaction move the deceleration, jerk, snap, and Om/Om3 diagnostics away from the $\\Lambda$CDM behavior, with viscosity playing the most visible role.","feed_headline":"Viscous interacting dark energy fits supernova data","feed_subtitle":"Generalized coupling plus bulk viscosity returns h0 = 0.68 and equation of state -1.07.","key_machinery":"The central object is the generalized Hubble parameter for interacting viscous dark energy, written in equation (16) as a modification of the viscous dark energy expression (15), with the interaction entering through the transfer $Q = 3H(\\lambda_{\\rm de}\\rho_{\\rm de} + \\lambda_\\chi\\rho_\\chi)$. This quantity carries the argument because it is used both to compute the distance moduli compared with Union 2.1 data and to generate the evolution of density parameters and the higher-order diagnostics. The model's six free parameters are the present Hubble constant $h_0$, the coupling strength $\\lambda_r$, the coupling mixing angle $\\lambda_\\theta$, the dark energy equation of state $\\omega_{\\rm de}$, the curvature density $\\Omega_{k,0}$, and the dimensionless viscosity $\\eta$ defined through $\\zeta = \\eta H$. The diagnostics side uses the Om parameter and the three-point Om3 statistic as null tests against $\\Lambda$CDM.","core_discovery":"The central discovery the paper argues for is that a single generalized model, obtained by writing the dark-sector energy transfer as $Q = 3H(\\lambda_{\\rm de}\\rho_{\\rm de} + \\lambda_\\chi \\rho_\\chi)$ and parametrizing the couplings by $\\lambda_\\chi = \\lambda_r \\sin\\lambda_\\theta$ and $\\lambda_{\\rm de} = \\lambda_r \\cos\\lambda_\\theta$, can reproduce the observed distance-redshift relation once dark energy is given a bulk viscosity $\\zeta = \\eta H$. With the Union 2.1 data the MCMC fit returns $h_0 = 0.68$, $\\lambda_r = 0.09$, $\\lambda_\\theta = 0.41$, $\\omega_{\\rm de} = -1.07$, $\\Omega_{k,0} = 0.10$, and $\\eta = 0.05$. The paper presents this as a mildly phantom dark energy with small viscosity and nonzero coupling, and it claims that the resulting model produces characteristic deviations from $\\Lambda$CDM in the density-parameter evolution and in the $q$, $j$, $s$, and Om/Om3 diagnostics. It also notes that the curvature parameter has little dynamical effect and that lower interaction strengths are preferred.","pith_inferences":["A natural extension would be to re-run the same generalized model on the Pantheon+ or DES-SN samples; those data sets would show whether the mildly phantom value $\\omega_{\\rm de} \\approx -1.07$ is a Union 2.1 artifact or a stable feature of viscous interacting dark energy.","The missing density evolution can be supplied by numerically integrating the coupled ODEs; such an independent computation would either reproduce the reported best fit or expose the degree to which the printed Hubble parameter is approximate.","The parametrization in terms of $\\lambda_r$ and $\\lambda_\\theta$ invites a physical reading: the ratio $\\lambda_\\chi/\\lambda_{\\rm de} = \\tan\\lambda_\\theta$ determines whether dark matter or dark energy is the dominant receiver of the energy transfer, and the posterior's broad peak near $7\\pi/10$ suggests the data favor a transfer dominated by the dark energy coupling.","To break degeneracies among $\\omega_{\\rm de}$, $\\eta$, and $\\lambda_r$, a joint fit that adds CMB and BAO data would be needed; the present single-probe analysis cannot separate these effects."],"forward_implications":["If the best-fit parameters are right, the present acceleration is driven by a mildly phantom dark energy ($\\omega_{\\rm de} \\approx -1.07$) with a small bulk viscosity ($\\eta \\approx 0.05$) and a weak but nonzero dark-sector coupling ($\\lambda_r \\approx 0.09$).","The model predicts that increasing viscosity slows the late-time expansion, so precise measurements of $H(z)$ around $z \\lesssim 2$ could distinguish a viscous contribution from a pure cosmological constant.","The Om3 diagnostic, which is a constant unity for $\\Lambda$CDM, is expected to dip below unity at higher redshifts and approach quintessence-like values at low redshifts, a signature that can be searched for in independent expansion data.","Because lower $\\lambda_r$ values are preferred and the curvature posterior is nearly flat, the model anticipates that future distance data will tighten the dark-sector coupling while leaving $\\Omega_{k,0}$ poorly constrained."],"supporting_citations":[{"why":"Supplies the generalized viscous dark energy Hubble parameter expression that the paper modifies for interactions.","marker":"[32]"},{"why":"Gives the coupled dark matter–dark energy evolution equations and the energy transfer $Q$ used to build the generalized interaction.","marker":"[34]"},{"why":"One of the benchmark interacting dark energy models whose $Q = 3\\lambda H\\rho_{\\rm de}$ is reproduced by the generalized coupling.","marker":"[35]"},{"why":"Another benchmark model, $Q = 3\\lambda H\\rho_\\chi$, folded into the generalized coupling parametrization.","marker":"[36]"},{"why":"Defines the Om3 three-point diagnostic used to distinguish the model from $\\Lambda$CDM.","marker":"[33]"},{"why":"Provides the adaptive MCMC algorithm used for the Bayesian fit to Union 2.1.","marker":"[50]"},{"why":"Provides the adaptive Metropolis method that underlies the MCMC package used in the analysis.","marker":"[51]"},{"why":"Eckart's theory defines the bulk-viscous effective pressure that turns the dark energy into a viscous fluid.","marker":"[21]"}],"fun_headline_variants":["Coupled viscous dark energy yields h0=0.68 from supernovae","Phantom viscous coupling best explains supernova distances","Generalized interacting viscous DE fits Union2.1 with h0=0.68","Viscous dark sector with energy swap matches SNe Ia data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Hubble parameter actually used in the MCMC fit follows from the model's coupled dark matter and dark energy density equations; the paper never displays the evolution of $\\rho_\\chi$ and $\\rho_{\\rm de}$ that would turn equation (16) into concrete distance moduli, so the fit could in principle be disconnected from the stated physics.","fun_headline_variants_meta":{"raw":{"variants":["Coupled viscous dark energy yields h0=0.68 from supernovae","Phantom viscous coupling best explains supernova distances","Generalized interacting viscous DE fits Union2.1 with h0=0.68","Viscous dark sector with energy swap matches SNe Ia data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00069,"raw_usage":{"total_tokens":3122,"prompt_tokens":941,"completion_tokens":2181,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":2112}},"tokens_in":557,"tokens_out":2181,"duration_ms":16512,"temperature":1.0,"reasoning_tokens":2112,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:29:04.496557+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the distance moduli from equations (1), (2), (3), and (16) using the reported best-fit parameters and compare them with the Union 2.1 likelihood; if the resulting posteriors differ from Figure 7, the stated parameters are not the model's actual best fit. A simpler check is to solve the coupled density equations numerically with the reported parameters and verify that the Hubble parameter used in the analysis is recovered.","supporting_citations":[{"cited_title":"Global 21-cm brightness temperature in viscous dark energy models","cited_arxiv_id":"2207.09177","evidence_quote":"Supplies the generalized viscous dark energy Hubble parameter expression that the paper modifies for interactions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the coupled dark matter–dark energy evolution equations and the energy transfer $Q$ used to build the generalized interaction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"One of the benchmark interacting dark energy models whose $Q = 3\\lambda H\\rho_{\\rm de}$ is reproduced by the generalized coupling."},{"cited_title":"Di Valentino, A","cited_arxiv_id":null,"evidence_quote":"Another benchmark model, $Q = 3\\lambda H\\rho_\\chi$, folded into the generalized coupling parametrization."},{"cited_title":"Haario, M","cited_arxiv_id":null,"evidence_quote":"Provides the adaptive MCMC algorithm used for the Bayesian fit to Union 2.1."},{"cited_title":"Haario, E","cited_arxiv_id":null,"evidence_quote":"Provides the adaptive Metropolis method that underlies the MCMC package used in the analysis."},{"cited_title":"Eckart, The thermodynamics of irreversible processes","cited_arxiv_id":null,"evidence_quote":"Eckart's theory defines the bulk-viscous effective pressure that turns the dark energy into a viscous fluid."}],"review_version":1}