{"id":"4b87c5ec-88d8-4d09-9710-eef0a7dfcfca","arxiv_id":"2501.03388","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"Bulk viscous matter in a linear f(T) gravity model is fitted to cosmic data, but the velocity-dependent viscosity case is built on an incorrect Hubble parameter solution.","lead":"A cosmology paper fits a model where dark matter is treated as a viscous fluid in modified teleparallel gravity, and finds it can mimic cosmic acceleration. The model's main case, however, contains a mathematical inconsistency in the derived Hubble formula, so the headline results do not hold as stated.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (24) is not the solution of Equation (23) for Case I; all Case I results are therefore invalid.","rationale":"The reader's verdict is REJECT, and my analysis confirms the same load-bearing defect: the explicit Hubble rate used for Case I is not the solution of the model's own differential equation. This is not a matter of interpretation or of disagreement with a consensus assumption; it is an internal inconsistency in the manuscript's central derivation. The ODE (23) is linear and first-order in H, so its solution is uniquely determined up to the initial condition; any expression put forward as the solution must satisfy it identically. Eq. (24) fails this test whenever ζ1≠0, which is precisely Case I. Because the MCMC likelihood uses H(z) directly, the quoted posterior contours and all derived cosmological diagnostics for Case I are based on a function that does not follow from the field equations. The abstract's headline numbers (q0≈−0.49, H0=60.0) come from Case I and are therefore not credible. Case II is a correct but routine fit with no model-selection analysis, so it cannot by itself salvage the paper's central claim. Thus the reader's REJECT verdict stands. I see no additional concern that would change the verdict, and no reason to soften it: the error is in the core analytic result, not in a peripheral detail.","tokens_in":1047,"tokens_out":1032,"duration_ms":39468,"concrete_test":"Substitute Eq. (24) into Eq. (23) at, say, z=1 using the Case I best-fit values (H0=60.0, α=1.01, ζ0=40.1, ζ1=0.123); the left-hand side will not vanish, proving Eq. (24) is not a solution. The definitive check is to re-solve Eq. (23) analytically as H(z)=H0(1+z)^A + ζ0/(α−ζ1)[1−(1+z)^A] with A=3(α−ζ1)/(2α), then re-run the MCMC on the combined dataset with this corrected H(z); if the recovered best-fit parameters and derived q0, j0, and diagnostics shift from the quoted values, the Case I conclusions are unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central Case I claim—that a bulk viscous f(T) model with ζ1≠0 fits H(z)+Pantheon+BAO and gives q0≈−0.49—rests on Eq. (24). But Eq. (24) is not a solution of the model's own Eq. (23). Rewriting Eq. (23) as dH/d ln a + A H = B, with A = 3(α−ζ1)/(2α) and B = 3ζ0/(2α), the general solution is H(z) = H0(1+z)^A + (B/A)[1−(1+z)^A] = H0(1+z)^{3(α−ζ1)/(2α)} + ζ0/(α−ζ1)[1−(1+z)^{3(α−ζ1)/(2α)}]. The second term of Eq. (24) has this form, but the first term is H0(1+z)^{3/2}, with exponent 3/2 instead of A. For ζ1≠0, A≠3/2, so Eq. (24) does not satisfy Eq. (23); direct substitution leaves a residual (A−3/2)H0(1+z)^{3/2}. Every Case I result—best-fit parameters, the quoted H0=60.0, q0≈−0.49, j0≈0.68, and the statefinder and Om(z) diagnostics—is computed from this incorrect expression. Case II (ζ1=0) is internally consistent but is only a two-parameter fit without model comparison or error propagation; the paper's advertised 'path to cosmic acceleration' depends on the invalid Case I.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript considers a flat FLRW universe in f(T)=αT teleparallel gravity whose matter content is a pressureless fluid with bulk viscosity ζ(t)=ζ0+ζ1H. It solves the background Friedmann equation for two cases (ζ1≠0 and ζ1=0), fits the resulting H(z) to the combined H(z)+Pantheon++BAO dataset with an MCMC, and studies the deceleration, jerk, effective EoS, statefinder, and Om(z) diagnostics. The abstract reports accelerating solutions with q0≈−0.49 (Case I) and q0≈−0.32 (Case II) and concludes that bulk viscous matter in f(T) gravity can explain cosmic acceleration without dark energy.","tokens_in":16351,"tokens_out":9631,"duration_ms":90289,"significance":"The idea is not without interest: a single viscous matter component in a simple f(T) gravity could in principle mimic dark energy, and the paper uses standard cosmological datasets with a clearly stated likelihood. The analytical derivation is also checkable, which is a strength. However, the main announced result (Case I) is built on an incorrect solution of the model's own field equation, so none of the Case I constraints or diagnostics can be trusted. In addition, the paper's diagnostics are deterministic functions of the fitted parameters and are not independent tests. As it stands, the manuscript does not provide credible evidence for the claimed path to cosmic acceleration.","major_comments":[{"comment":"Equation (24) is not the solution of Eq. (23) for ζ1≠0. Writing Eq. (23) as dH/d ln a + A H = B with A=3(α−ζ1)/(2α) and B=3ζ0/(2α), the general solution is H(z)=H0(1+z)^A + (B/A)[1−(1+z)^A] = H0(1+z)^{3(α−ζ1)/(2α)} + ζ0/(α−ζ1)[1−(1+z)^{3(α−ζ1)/(2α)}]. Equation (24) instead uses H0(1+z)^{3/2} for the homogeneous term, so direct substitution leaves a residual (A−3/2)H0(1+z)^{3/2}. Since ζ1≠0, A≠3/2 and Eq. (24) does not satisfy Eq. (23). Every Case I result—the best-fit parameters in Fig. 1, H(z) in Fig. 3, q0≈−0.49, j0≈0.68, and the statefinder and Om(z) curves—is computed from this invalid expression and must be either recomputed with the correct solution or removed.","section":"III, Eqs. (23)–(24)"},{"comment":"The claimed confirmation of acceleration is partly circular. The deceleration, jerk, EoS, statefinder, and Om(z) are evaluated by inserting the best-fit parameters from Section IV into the same H(z) that was fitted to the data; they are not independent predictions. To substantiate the title's claim, the paper needs at least a model comparison (for example Δχ², AIC, or DIC relative to ΛCDM and to Case II) and a demonstration that the fitted model is statistically preferred. Without this, the statement that the model explains acceleration is largely a restatement of the fit.","section":"V–VII"}],"minor_comments":[{"comment":"The MCMC description should report the burn-in length, convergence diagnostics (such as Gelman-Rubin statistics), and the final χ² per degree of freedom; otherwise the quoted 1σ intervals cannot be independently checked.","section":"IV"},{"comment":"The text says 'the present-day values of the jerk parameter are approximately ω0≈−0.78 and ω0≈−0.55', but ω0 is the effective EoS parameter, not the jerk parameter; this sentence should be corrected.","section":"V.C and Conclusion"},{"comment":"The units of ζ0 and ζ1 should be reconciled with the use of H in km/s/Mpc; as written, ζ=ζ0+ζ1H mixes quantities whose physical dimensions differ unless H is converted to s⁻¹ or natural units are explicitly specified.","section":"III and IV"},{"comment":"The diagnostic curves are shown only for best-fit parameter values; including MCMC uncertainty bands would make the plots more informative and would also clarify how robust the quoted q0, j0, and ωeff values are.","section":"Figures 4–11"}],"recommendation":"reject","confidential_remarks":"The error in Eq. (24) is elementary and invalidates the primary Case I analysis, which is the advertised central result of the paper. In my view rejection is the appropriate outcome; if the authors correct the ODE solution and redo the full fitting, error propagation, and model comparison, a substantially revised manuscript could be considered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a routine viscous-cosmology fit in f(T) gravity, and the main advertised case is invalid. Equation (24) is not the solution of Eq. (23) when ζ1 ≠ 0. Solving Eq. (23) gives H(z) = H0(1+z)^(3(α−ζ1)/(2α)) + [ζ0/(α−ζ1)](1 − (1+z)^(3(α−ζ1)/(2α))). The paper drops the factor 3(α−ζ1)/(2α) from the first exponent and uses 3/2 instead. Direct substitution leaves a residual (A − 3/2)H0(1+z)^(3/2), with A = 3(α−ζ1)/(2α). Every Case I number in the paper — H0 = 60, α = 1.01, ζ0 = 40.1, ζ1 = 0.123, q0 ≈ −0.49, j0 ≈ 0.68, the statefinder and Om(z) plots — comes from that wrong expression. So the central claim falls.\n\nWhat the paper does do well: the f(T) setup is standard and the review of viscous cosmology is adequate. Case II (ζ1 = 0) is internally consistent, and the fit to H(z)+Pantheon+BAO using emcee is a routine but legitimate exercise. The transition redshift and q0 values for Case II are plausible. That is about the limit of the positive content.\n\nThe other soft spots are secondary but real. The diagnostics (q, j, ω_eff, statefinder, Om(z)) are deterministic functions of the fitted parameters, so calling them \"predictions\" overstates the case. There is no model comparison — no χ² for ΛCDM to see whether the viscous term is even preferred. The MCMC setup is stated but the priors, burn-in, and convergence checks are thin. And Case I pulls H0 to 60, which is hard to take seriously even before the math error.\n\nBottom line: this paper is not publishable as is. The error in Eq. (24) is load-bearing, not cosmetic. If the authors correct the exponent and refit, they get a different answer — and they need to include a comparison with ΛCDM and a discussion of why the viscous model is competitive. As submitted, I would not send it to a referee; it is a routine fit with a fatal flaw in its main case.","headline":"The Case I results rest on an H(z) that does not solve the paper's own Eq. (23), so the advertised path to cosmic acceleration is unsupported.","tokens_in":16893,"tokens_out":2789,"would_cite":false,"duration_ms":24122,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","83D05"],"pacs":["98.80.-k","04.50.Kd"],"model":"deepseek-v4-flash","headline":"Bulk viscous matter in $f(T)=\\alpha T$ gravity, with viscosity $\\zeta(t)=\\zeta_0+\\zeta_1 H$, fits combined cosmological data and can drive the observed transition to accelerated expansion without a dark-energy component.","keywords":["bulk viscosity","f(T) teleparallel gravity","cosmic acceleration","deceleration parameter","statefinder diagnostic","Om(z) diagnostic","Hubble parameter constraints","quintessence"],"falsifier":"Direct substitution of Eq. (24) into Eq. (23) fails for $\\zeta_1\\neq0$: the homogeneous term needs the exponent $3(\\alpha-\\zeta_1)/(2\\alpha)$, not $3/2$. Re-running the MCMC with the correctly integrated $H(z)$ would settle whether the quoted Case I best-fit parameters and $q_0\\approx-0.49$ survive.","tokens_in":15806,"feed_emoji":"🌌","tokens_out":10211,"duration_ms":86009,"temperature":0.7,"pith_summary":"This paper sets out to show that the universe's late-time acceleration can be produced by bulk viscosity in the matter fluid, within teleparallel modified gravity with $f(T)=\\alpha T$. It takes the viscosity coefficient to be $\\zeta(t)=\\zeta_0+\\zeta_1 H$, derives an explicit Hubble parameter $H(z)$ from the modified Friedmann equations, and fits that expression to the combined $H(z)+\\text{Pantheon}^{+}+\\text{BAO}$ dataset. With the best-fit parameters the deceleration parameter changes from positive to negative, so the expansion switches from decelerating to accelerating around redshift $z\\approx0.8$--$0.9$. The present-day values are $q_0\\approx-0.49$ (for $\\zeta_1\\neq0$) and $q_0\\approx-0.32$ (for $\\zeta_1=0$), which the paper reads as evidence that a bulk-viscous matter component in $f(T)$ gravity can replace dark energy.","feed_headline":"Data fits show bulk viscosity could drive cosmic acceleration","feed_subtitle":"A viscous-matter f(T) model fits H(z), Pantheon+ and BAO data with a late-time deceleration-to-acceleration transition.","key_machinery":"The object carrying the argument is the explicit Hubble parameter $H(z)$ obtained by integrating the modified Friedmann equation under the linear torsion model $f(T)=\\alpha T$ and the bulk-viscosity prescription $p_{\\mathrm{eff}}=p-3\\zeta(t)H$. That single function enters every subsequent diagnostic: the deceleration parameter $q$, jerk $j$, effective equation of state $\\omega_{\\mathrm{eff}}$, statefinder pair $(r,s)$, and $Om(z)$ are all computed from $H(z)$ or its derivatives, so the fit quality and the predicted cosmic acceleration depend entirely on this expression and on the two fitted viscosity parameters $\\zeta_0$ and $\\zeta_1$.","core_discovery":"The central claim is that a pressureless matter fluid with bulk viscosity $\\zeta(t)=\\zeta_0+\\zeta_1 H$, coupled through $f(T)=\\alpha T$ teleparallel gravity, yields the Hubble-rate solution $H(z)=H_0(1+z)^{3/2}+\\frac{\\zeta_0}{\\alpha-\\zeta_1}\\left[1-(1+z)^{\\frac{3(\\alpha-\\zeta_1)}{2\\alpha}}\\right]$, and that this solution, when fitted to $H(z)+\\text{Pantheon}^{+}+\\text{BAO}$ data, gives an accelerating universe without dark energy. For $\\zeta_1\\neq0$ the best fit is $H_0=60.0^{+2.0}_{-1.9}$ km/s/Mpc, $\\alpha=1.01^{+0.10}_{-0.098}$, $\\zeta_0=40.1^{+1.9}_{-2.0}$, and $\\zeta_1=0.123^{+0.093}_{-0.088}$; for $\\zeta_1=0$ it is $H_0=67.5^{+1.3}_{-1.3}$ km/s/Mpc, $\\alpha=0.94^{+0.14}_{-0.13}$, and $\\zeta_0=34.7^{+2.0}_{-2.0}$. The resulting jerk parameter, effective equation of state, statefinder trajectories, and $Om(z)$ diagnostic all indicate quintessence-like behavior, with the model approaching the $\\Lambda$CDM fixed point in the far future.","pith_inferences":["The same derivation could be repeated for nonlinear $f(T)$ models, such as adding a quadratic torsion term, where the torsion correction and the viscous pressure would interact; those extended models could be tested against the same $H(z)$, Pantheon+, and BAO data.","The gap between the inferred $H_0$ values in the two cases ($\\approx60$ versus $\\approx67.5$ km/s/Mpc) suggests that deciding whether viscosity scales with $H$ strongly affects the Hubble-constant estimate; a formal model-selection comparison between the two cases is a natural next step.","Because every quoted diagnostic is a functional of the single $H(z)$ expression, a reader can reconstruct the entire analysis from Eqs. (20)--(24) and the quoted best-fit parameters; no additional closure assumptions are needed to reproduce the figures."],"forward_implications":["If the model is correct, no cosmological constant or scalar field is needed: the negative effective pressure that accelerates the expansion comes from the bulk-viscosity term $-3\\zeta H$ in a universe dominated by ordinary matter.","The combined dataset fixes the parameters tightly enough to predict a transition redshift $z_{\\mathrm{tr}}\\approx0.90$ for $\\zeta_1\\neq0$ and $z_{\\mathrm{tr}}\\approx0.80$ for $\\zeta_1=0$, with present-day deceleration $q_0\\approx-0.49$ and $-0.32$.","The statefinder and $Om(z)$ trajectories lie in the quintessence region and converge to the $\\Lambda$CDM fixed point $(r,s)=(1,0)$ and the de Sitter point $(r,q)=(1,-1)$ in the future, so the model mimics $\\Lambda$CDM at late times while differing at intermediate redshifts.","The fitted $\\alpha$ values are close to $1$, so the model can be viewed as a small torsion-gravity correction to the standard Friedmann dynamics of a viscous fluid."],"supporting_citations":[{"why":"Supplies the bulk-viscosity ansatz $\\zeta(t)=\\zeta_0+\\zeta_1 H$ used throughout the model.","marker":"[74-76]"},{"why":"Provides the modified Friedmann equations for $f(T)$ gravity with a viscous fluid that the paper integrates.","marker":"[77]"},{"why":"The emcee MCMC package used to constrain the model parameters.","marker":"[78]"},{"why":"The 31 Hubble parameter measurements entering the $H(z)$ likelihood.","marker":"[79-85]"},{"why":"The Pantheon+ supernova sample and its covariance matrix used in the joint fit.","marker":"[87, 88]"},{"why":"The BAO measurements from SDSS, 6dFGS, and BOSS used in the joint fit.","marker":"[89-94]"},{"why":"Planck values used as the reference for $H_0$ and for the $\\Lambda$CDM comparison curves.","marker":"[95]"},{"why":"Defines the statefinder pair $(r,s)$ used to classify the model as quintessence-like.","marker":"[119, 120]"},{"why":"Defines the $Om(z)$ diagnostic whose negative slope is used to support quintessence behavior.","marker":"[121]"}],"fun_headline_variants":["Bulk viscosity in f(T) gravity drives cosmic acceleration","Viscous matter model fits data, no dark energy needed","f(T) gravity plus bulk viscosity yields quintessence-like expansion","Data favors viscous f(T) acceleration"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Hubble parameter expression in Eq. (24) is the correct integral of Eq. (23) when $\\zeta_1\\neq0$; if that integration is not correct, every Case I parameter constraint and diagnostic built on it has to be recomputed.","fun_headline_variants_meta":{"raw":{"variants":["Bulk viscosity in f(T) gravity drives cosmic acceleration","Viscous matter model fits data, no dark energy needed","f(T) gravity plus bulk viscosity yields quintessence-like expansion","Data favors viscous f(T) acceleration"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000683,"raw_usage":{"total_tokens":3261,"prompt_tokens":1268,"completion_tokens":1993,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":884,"completion_tokens_details":{"reasoning_tokens":1939}},"tokens_in":884,"tokens_out":1993,"duration_ms":13300,"temperature":1.0,"reasoning_tokens":1939,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:52:16.494212+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Direct substitution of Eq. (24) into Eq. (23) fails for $\\zeta_1\\neq0$: the homogeneous term needs the exponent $3(\\alpha-\\zeta_1)/(2\\alpha)$, not $3/2$. Re-running the MCMC with the correctly integrated $H(z)$ would settle whether the quoted Case I best-fit parameters and $q_0\\approx-0.49$ survive.","supporting_citations":[{"cited_title":"Avelino, U","cited_arxiv_id":null,"evidence_quote":"Provides the modified Friedmann equations for $f(T)$ gravity with a viscous fluid that the paper integrates."},{"cited_title":"Sharif and S","cited_arxiv_id":null,"evidence_quote":"The emcee MCMC package used to constrain the model parameters."},{"cited_title":"Giostri et al., J","cited_arxiv_id":null,"evidence_quote":"Planck values used as the reference for $H_0$ and for the $\\Lambda$CDM comparison curves."},{"cited_title":"alam et al., Mon","cited_arxiv_id":null,"evidence_quote":"Defines the $Om(z)$ diagnostic whose negative slope is used to support quintessence behavior."}],"review_version":1}