{"id":"ca0463b1-5e60-4437-b330-fec4903afd90","arxiv_id":"2504.19523","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A Bianchi type I bulk-viscous model with pressure -3ζH^2, fitted to OHD, Union 2.1, Pantheon, and BAO data, yields H0 estimates near 66.9 and 74.2 km/s/Mpc and late-time quintessence-like acceleration.","lead":"This paper fits a Bianchi type I universe with a bulk viscous fluid to Hubble, supernova, and BAO data, and reports best-fit values of H0, a shear parameter, and a viscosity coefficient. A generalist might read it to see whether a simple viscous-fluid model can mimic late-time cosmic acceleration, but internal data inconsistencies prevent the claimed constraints from being trusted.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"BAO likelihood and combined-data χ² are internally inconsistent; the headline H0 constraints are unsupported as presented.","rationale":"Good-faith reading: the paper derives a consistent LRS Bianchi I solution—Eq. (14) follows algebraically from Eq. (12)—and it attempts to fit that solution to standard public data. That is a legitimate way to constrain a phenomenological viscous model, and the derivation of H(z) itself is not the main weakness. The decisive problem is that the numerical pipeline described in the manuscript cannot be reproduced as written. The BAO section first gives one data vector and then evaluates a different vector in Eq. (24), while the combined-data likelihood in Eq. (26) adds Pantheon to a fit labeled OHD+BAO and the four-set combination stacks Union 2.1 and Pantheon without handling overlap. These are internal inconsistencies, not merely choices outside current consensus. They directly determine the fitted parameters in Table 5 and therefore the derived age, q0, transition redshift, and statefinder claims. The reader's weakest_assumption about the viscous pressure law −3ζH² is a fair scientific caveat, but it is less load-bearing than the likelihood corruption: the pressure law is explicitly assumed, whereas the headline numbers are claimed to be estimates from the data. I agree with the REJECT verdict, with a partial disagreement on where the weakest point lies. A corrected BAO likelihood and combined-χ² rule could potentially rescue the analysis, so the rejection is based on the analysis as presented rather than on the model being impossible.","tokens_in":19940,"tokens_out":7551,"duration_ms":74032,"concrete_test":"Recompute the Table 5 constraints from the manuscript's own equations: (i) replace the vector in Eq. (24) with either the zBAO/dz/σ list from §III.D using dz=rs(z*)/Dv(z), or the original Giostri et al. data with its correct covariance; (ii) set χ²_OHD+BAO = χ²_OHD + χ²_BAO only, and for the four-set fit use OHD+Pan or OHD+Union plus BAO, not both SN compilations; (iii) rerun the same MCMC with these corrected likelihoods. If the resulting H0/l/ζ posteriors differ from Table 5 by more than the quoted 1σ errors, the headline constraints are artifacts of the corrupted likelihood rather than the data.","verdict_should_be":"REJECT","load_bearing_attack":"The headline constraints in Table 5 rest on a BAO likelihood whose data vector, covariance, and point count do not match each other, and on a combined-χ² rule that does not correspond to the named data sets. In §III.D the authors list a six-point BAO sample with dz(zBAO)=[30.95,17.55,10.11,8.44,6.69,5.45] and σ=[1.46,0.60,0.37,0.67,0.33,0.31], and define dz=rs(z*)/Dv in Eq. (22). Equation (24) then evaluates a different vector, dA(z*)/Dv − [30.84,10.33,6.72,8.41,6.66,5.43], at redshifts {0.106,0.35,0.57,0.44,0.60,0.73}; neither the redshifts nor the values match §III.D, and dA(z*) is not the sound-horizon ratio defined in Eq. (22). The abstract and text call the BAO set '5 redshifts' while six points are used with a 6×6 covariance. In §III.E, Eq. (26) defines the 'OHD+BAO' likelihood as χ²_OHD + χ²_PAN + χ²_BAO, silently adding Pantheon to a fit labeled OHD+BAO, while the 'OHD+Pan+BAO+Union' fit stacks Union 2.1 and Pantheon without accounting for overlapping SNIa or the paper's own 680/580/1048 dataset-size contradictions. Since H0, l, and ζ in Table 5 are direct outputs of these likelihoods, the reported H0 values 66.912 and 74.216, and all derived age, transition-redshift, and statefinder conclusions evaluated at those parameters, are not supported by the analysis as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs an LRS Bianchi type-I cosmological model with a bulk viscous fluid whose effective pressure is p_eff = -3ζH^2 and p = 0 for dust, solves the field equations to obtain H(z) in Eq. (14), and fits the three free parameters H0, l, ζ to OHD, Union 2.1, Pantheon, and BAO data, individually and in combinations. The headline results are H0 = 66.912^{+0.497}_{-0.501} and 74.216^{+0.150}_{-0.148} km/s/Mpc for the combined fits, together with derived ages, transition redshifts, statefinder, Om, and ω-ω' diagnostics, leading to the claim that the model behaves like quintessence and approaches ΛCDM.","tokens_in":20296,"tokens_out":8547,"duration_ms":74693,"significance":"The paper's strength is an exact anisotropic solution with a simple viscous-pressure ansatz and an explicit MCMC/χ² estimation procedure. If the statistical analysis were valid, the model would be a useful worked example of a bulk-viscous mechanism for late-time acceleration with a deceleration-to-acceleration transition. However, the statistical basis for the central claims is internally inconsistent: the BAO likelihood is not reproducible from the stated equations, the Union 2.1 sample size changes across sections, and the combined likelihood does not match the data sets named in Table 5. The late-time diagnostics are also functions of the same fitted parameters rather than independent tests. The headline H0 constraints and the quintessence/ΛCDM conclusion are therefore unsupported as presented.","major_comments":[{"comment":"The BAO likelihood used for Table 4 and for the combined fits is not consistently defined. The text lists a six-point sample with zBAO = [0.106, 0.2, 0.35, 0.44, 0.6, 0.73] and dz = [30.95, 17.55, 10.11, 8.44, 6.69, 5.45], but Eq. (24) evaluates X at the redshifts [0.106, 0.35, 0.57, 0.44, 0.60, 0.73] against the reference values [30.84, 10.33, 6.72, 8.41, 6.66, 5.43]; neither the points nor the values match those in §III.D. In addition, Eq. (22) defines the distance ratio as dz(z) = rs(z*)/Dv(z), while Eq. (24) uses dA(z*)/Dv(z), which is a different quantity. Because Eq. (23) defines χ²_BAO = X^T C^{-1} X with this X, the BAO constraints in Tables 4 and 5 are not reproducible from the stated ingredients.","section":"§III.D, Eqs. (22)–(25)"},{"comment":"The sample sizes used in the fits are mutually inconsistent. The Union 2.1 compilation is described as 650 data sets in §III.B, as 580 in Eq. (19) and the abstract, and as 518 in the caption of Fig. 2b and in the corner-plot text. For BAO, the abstract and conclusion say '5 redshifts', while §III.D and Eqs. (23)–(25) use six points with a 6×6 covariance matrix. These discrepancies change the χ² normalization and degrees of freedom and therefore affect the parameter values and uncertainties reported in Tables 1–5.","section":"§III.B and §III.D"},{"comment":"The combined likelihood in Eq. (26), labeled χ²_OHD+BAO, is defined as χ²_OHD + χ²_PAN + χ²_BAO, so it silently adds the Pantheon sample to a fit that Table 5 and the abstract call OHD+BAO. The 'OHD+Pan+BAO+Union' fit stacks Union 2.1 and Pantheon without accounting for the overlap between SNIa compilations or for the size contradictions noted above. Since the headline values H0 = 66.912 and H0 = 74.216 in Table 5, and all derived quantities evaluated at those parameters in Sections IV–X, are direct outputs of these likelihoods, the central quantitative claims are unsupported by the analysis as written.","section":"§III.E, Eq. (26)"},{"comment":"The late-time diagnostics are not independent checks of the model. The quantities q0, zt, t0, the statefinder pair (r,s), Om(z), and the ω-ω' trajectory are all evaluated at the fitted H0, l, and ζ through Eq. (14), so the statement in the abstract and conclusion that the model 'behaves like quintessence and approaches ΛCDM' is a restatement of the fitted viscosity parameter rather than a falsifiable prediction. No comparison with a ΛCDM fit to the same data or model-selection statistics (Δχ², AIC, or BIC) is provided, so the claim of good agreement with observations in Sec. XI is not quantitatively supported.","section":"§VI–§X, Eqs. (36)–(46)"},{"comment":"Equation (34) reports the present density for the H0 = 74.216 fit as ρ0 = 4.06629×10^{-31} g/cm^3. Evaluating Eq. (33) at z = 0 gives ρ0/ρc = 1 - l^2/9 ≈ 0.99, which with H0 = 74.216 km/s/Mpc yields ρ0 ≈ 1.03×10^{-29} g/cm^3; the reported value is too small by a factor of about 25. The corresponding value for the H0 = 66.912 fit in Eq. (35) is consistent with the formula, so this appears to be a numerical error in one of the two reported densities.","section":"§V, Eqs. (33)–(35)"}],"minor_comments":[{"comment":"The Data Availability statement, 'No data was used for the research described in the article,' is inconsistent with the extensive use of the OHD, Union 2.1, Pantheon, and BAO data sets in Sec. III.","section":"Data Availability"},{"comment":"The abstract and conclusion contain a stray closing parenthesis in '66.912^{+0.497}_{-0.501}) Km/s/Mpc', and the same H0 expression appears without consistent units in Table 5.","section":"Abstract and Sec. XI"},{"comment":"The caption of Fig. 5 and the legends of Figs. 8–14 use different names for the same combined fits ('OHD+BAO', 'OHD+BAO+Union2.1', 'Pantheon+OHD+BAO+UNION2.1'), making it hard to associate the curves with Table 5.","section":"Figs. 5, 8–14"},{"comment":"The paragraph in Sec. III.A describing the least-squares, χ², and MCMC procedure is repeated verbatim a few lines later in the same subsection; one copy should be deleted.","section":"§III.A"}],"recommendation":"reject","confidential_remarks":"The number and type of inconsistencies (dataset sizes, likelihood definitions, combined-data labels, and a factor-of-25 error in a derived density) suggest that the paper needs a full reanalysis rather than a routine revision. I would not recommend inviting a revision unless the authors redo the fits with clearly specified likelihoods and report model-selection statistics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The central H(z) formula, Eq. (14), is a legitimate solution of the stated field equations, and the authors do confront the model with public data and run a standard MCMC pipeline. That is the most that can be said in its favor, because the statistical analysis as written is not internally consistent, and the headline H0 values rest on likelihoods that do not match the data they claim to use.\n\nThe BAO section is the clearest problem. The paper lists a six-point sample with redshift vector [0.106, 0.2, 0.35, 0.44, 0.6, 0.73] and dz values [30.95, 17.55, 10.11, 8.44, 6.69, 5.45], but Eq. (24) evaluates a different vector at a different redshift set, including 0.57 instead of 0.2 and values like 30.84 and 10.33. The abstract and conclusion call it a 5-redshift sample while six points and a 6x6 covariance are used. Eq. (26) defines the 'OHD+BAO' chi-square as OHD plus Pantheon plus BAO, silently adding a supernova sample to a fit that is labeled otherwise. The Union 2.1 sample size appears as 680, 650, and 580 in different places. These are not cosmetic typos; H0, l, and zeta in Table 5 are direct outputs of these likelihoods, so the quoted 66.912 and 74.216 values are not reproducible from the information given.\n\nThere is also a legitimate circularity concern, though I would phrase it carefully. The late-time diagnostics (q0, zt, t0, statefinder, Om, omega-omega') are all functions of the fitted H0, l, and zeta, so the statement that the model 'behaves like quintessence and approaches Lambda CDM' is not an independent test; it is a restatement of the assumed viscous pressure law -3 zeta H^2 and the fitted parameters. The model's acceleration is put in by hand in that sense, not derived from new physics.\n\nThe data availability statement, saying 'No data was used,' is contradicted by the paper’s own use of OHD, Union 2.1, Pantheon, and BAO catalogs. No code or processed likelihoods are provided, which makes the internal inconsistencies harder to sort out.\n\nWhom is this for? Readers working on Bianchi-I viscous cosmologies might find the derived H(z) useful as a template, but the observational constraints section needs substantial correction before the numerical results can be trusted. I would not cite the current version. It deserves a serious referee in the sense that the modeling is coherent and the data-analysis problems are fixable; a desk rejection would be premature. But my recommendation is reject in current form, with an invitation to resubmit after the likelihoods, dataset sizes, and BAO vector are corrected and the code or exact likelihoods are made available.","headline":"A standard Bianchi-I viscous model with a consistent Hubble solution, but the data analysis is internally inconsistent and the headline H0 constraints are not supported as presented.","tokens_in":20881,"tokens_out":1793,"would_cite":false,"duration_ms":19477,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k"],"model":"deepseek-v4-flash","headline":"This paper claims that an axially symmetric Bianchi type-I universe with pressureless dust and bulk viscous pressure −3ζH² fits the observed late-time acceleration, giving H0 = 66.9 km/s/Mpc and 74.2 km/s/Mpc from two combined data sets.","keywords":["LRS Bianchi type-I model","bulk viscous fluid","barotropic viscous pressure","late-time cosmic acceleration","Hubble constant","Markov Chain Monte Carlo","statefinder diagnostic","quintessence"],"falsifier":"Use the best-fit Hubble law with, say, $H_0 = 74.216$, $l = 0.276$, $\\zeta = 0.661$ to predict distance moduli over $0<z<2.3$, and compare them with a supernova sample not used in the fit; a systematic deviation beyond the reported uncertainties would show that the viscous law cannot alone carry the acceleration. Alternatively, measure the deceleration-to-acceleration transition redshift from an independent cosmic-chronometer sample: the model predicts $z_t \\approx 1.17$ for the $H_0 = 74.216$ fit and $z_t \\approx 2.10$ for the $H_0 = 66.912$ fit, so a measured transition outside both ranges would contradict the constant-$\\zeta$ picture.","tokens_in":19668,"feed_emoji":"🌌","tokens_out":13968,"duration_ms":117266,"temperature":0.7,"pith_summary":"The paper tries to establish that the observed late-time acceleration of the universe can be produced by bulk viscosity alone, without a cosmological constant or scalar field, inside an axially symmetric Bianchi type-I spacetime. It assumes the cosmic fluid is pressureless dust with an additional barotropic viscous pressure of the form $-3\\zeta H^2$, and uses that assumption to derive an exact Hubble law $H(z)$ whose three parameters ($H_0$, anisotropy strength $l$, and viscosity coefficient $\\zeta$) are fitted to Hubble, supernova, and BAO data. The headline results are best-fit values of $H_0 = 66.912^{+0.497}_{-0.501}$ km/s/Mpc from the OHD+BAO combination and $H_0 = 74.216^{+0.150}_{-0.148}$ km/s/Mpc from OHD+Pan+BAO+Union. The paper further reports a deceleration-to-acceleration transition at redshift $z_t \\approx 1.17$ or $2.10$, and statefinder and Om diagnostics that place the late-time behavior in the quintessence region, approaching $\\Lambda$CDM. If right, the model offers a purely fluid-dynamical explanation for dark energy and a concrete set of parameters that future surveys can test.","feed_headline":"Bulk viscosity alone can drive cosmic acceleration, model says","feed_subtitle":"Fits to supernova, Hubble, and BAO data put the late-time expansion near ΛCDM with quintessence behavior.","key_machinery":"The load-bearing object is the effective pressure $p_{\\rm eff} = p - 3\\zeta H^2$, with dust pressure $p=0$ and a constant viscosity coefficient $\\zeta$. Because $H^2$ is always positive, the viscous term is always a negative pressure whose magnitude grows with the expansion rate; that is what drives the acceleration. The second ingredient is the anisotropy parameter $l = c_1/(H_0 a_0^3)$, which encodes the shear between the axial and transverse scale factors of the LRS Bianchi type-I metric. Together they reduce the field equations to a single ordinary differential equation for $H(z)$, whose closed-form solution is the function fitted to all data sets. The machinery converts a transport coefficient into a geometric expansion history.","core_discovery":"Within Einstein's field equations, the paper derives an exact Hubble parameter of the form $H(z) = H_0 (1+z)^{3/2} (1+z)^{-3\\zeta/2} \\sqrt{l^2((1+z)^{3\\zeta+3}-1)+9(\\zeta+1)}/(3\\sqrt{\\zeta+1})$. Fitting this law to the four data sets---46 Hubble parameter measurements, the Union 2.1 supernova compilation, the 1048-point Pantheon apparent-magnitude sample, and six BAO distance-ratio points---the paper reports best-fit parameters for $H_0$, $l$, and $\\zeta$ from each set and from two combined sets. It then computes the deceleration parameter $q(z)$, finds the sign change from deceleration to acceleration, and uses the statefinder pair $(r,s)$, the jerk parameter, the Om diagnostic, and the $\\omega$--$\\omega'$ plane to classify the dark-energy behavior. The central discovery claim is that the viscous term $-3\\zeta H^2$ with constant $\\zeta$ is sufficient to make the model transition from deceleration to acceleration and to approach $\\Lambda$CDM at late times, with fitted $H_0$ values that lie on both sides of the current Hubble tension.","pith_inferences":["Because $-3\\zeta H^2$ is assumed rather than derived, the same expansion history could be reinterpreted as an effective dark-energy equation of state; deriving $\\zeta$ from kinetic theory or from an underlying field would turn the fit into a prediction.","The two combined-data $H_0$ values differ by roughly 7 km/s/Mpc, so a single viscosity law may not resolve the Hubble tension by itself; fitting all four data sets simultaneously with the full covariance matrix would show whether one parameter set can satisfy everything.","A natural extension is to use the fitted $H(z)$ to compute the sound-horizon scale or CMB distance priors, which the paper fixes rather than fits; those would be independent cross-checks of the viscosity parameters.","The age estimates (12.4 and 15.4 Gyr) straddle the usually quoted 13.7 Gyr; comparing the predicted age-redshift relation with the oldest observed objects would be a direct test."],"forward_implications":["Late-time acceleration can be obtained with no cosmological constant: a single constant bulk-viscosity coefficient suffices.","The fitted values give $H_0 \\approx 66.9$ km/s/Mpc from OHD+BAO and $H_0 \\approx 74.2$ km/s/Mpc from OHD+Pan+BAO+Union, so the model can be tuned toward either side of the Hubble-tension range.","The predicted transition redshift $z_t \\approx 1.17$ or $2.10$ and present deceleration parameter $q_0 \\approx -0.48$ or $-0.39$ are direct observational targets for cosmic-chronometer and redshift-drift surveys.","The statefinder trajectory ending near $(r,s)=(1,0)$ and the Om slope classify the model as quintessence-like at late times, so data requiring a phantom crossing would falsify this viscosity mechanism."],"supporting_citations":[{"why":"supplies the 46-point Hubble-parameter (OHD) data set used for the first fit and as part of the combined fits.","marker":"[41]"},{"why":"supplies the Union 2.1 supernova distance-modulus compilation used in the SNIa fits and in combined fits.","marker":"[47]"},{"why":"supplies the 1048-point Pantheon apparent-magnitude sample used in the SNIa fits and combined fits.","marker":"[48]"},{"why":"supplies the BAO distance-ratio data and inverse covariance matrix that define the BAO chi-square used in the combined fits.","marker":"[70]"},{"why":"provides the CMB-based Hubble-constant value used as the lower-side comparison for the fitted H0.","marker":"[62]"},{"why":"provides the Cepheid-calibrated local Hubble-constant value used as the higher-side comparison for the fitted H0.","marker":"[61]"},{"why":"defines the statefinder pair (r,s) and the classification regions used to conclude the model behaves like quintessence.","marker":"[72]"},{"why":"defines the (r,s) and (r,q) planes and the ΛCDM fixed point used to say the model approaches ΛCDM.","marker":"[76]"}],"fun_headline_variants":["Bulk viscosity alone drives late-time cosmic acceleration","Viscous fluid model fits data, approaches ΛCDM","Bianchi-I viscous model matches supernova and BAO data","Viscous dark energy mimics quintessence in fit","One viscous term explains acceleration, model shows"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire late-time acceleration rests on the assumed viscous pressure law $p_{\\rm eff} = -3\\zeta H^2$ with $p=0$ and a constant viscosity coefficient $\\zeta$; if that negative-pressure form is wrong, or if $\\zeta$ varies with time or density, the reported $H_0$ values, transition redshifts, and the approach to $\\Lambda$CDM inherit the error rather than following from independent physics.","fun_headline_variants_meta":{"raw":{"variants":["Bulk viscosity alone drives late-time cosmic acceleration","Viscous fluid model fits data, approaches ΛCDM","Bianchi-I viscous model matches supernova and BAO data","Viscous dark energy mimics quintessence in fit","One viscous term explains acceleration, model shows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000182,"raw_usage":{"total_tokens":1422,"prompt_tokens":1166,"completion_tokens":256,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":782,"completion_tokens_details":{"reasoning_tokens":177}},"tokens_in":782,"tokens_out":256,"duration_ms":3233,"temperature":1.0,"reasoning_tokens":177,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:52:15.499010+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use the best-fit Hubble law with, say, $H_0 = 74.216$, $l = 0.276$, $\\zeta = 0.661$ to predict distance moduli over $0<z<2.3$, and compare them with a supernova sample not used in the fit; a systematic deviation beyond the reported uncertainties would show that the viscous law cannot alone carry the acceleration. Alternatively, measure the deceleration-to-acceleration transition redshift from an independent cosmic-chronometer sample: the model predicts $z_t \\approx 1.17$ for the $H_0 = 74.216$ fit and $z_t \\approx 2.10$ for the $H_0 = 66.912$ fit, so a measured transition outside both ranges would contradict the constant-$\\zeta$ picture.","supporting_citations":[{"cited_title":"Ren and X.-H","cited_arxiv_id":null,"evidence_quote":"supplies the 46-point Hubble-parameter (OHD) data set used for the first fit and as part of the combined fits."},{"cited_title":"Koussour, Abdelghani Errehymy, O","cited_arxiv_id":null,"evidence_quote":"supplies the Union 2.1 supernova distance-modulus compilation used in the SNIa fits and in combined fits."},{"cited_title":"Suzuki, et al., Astrophy","cited_arxiv_id":null,"evidence_quote":"supplies the 1048-point Pantheon apparent-magnitude sample used in the SNIa fits and combined fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the BAO distance-ratio data and inverse covariance matrix that define the BAO chi-square used in the combined fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the (r,s) and (r,q) planes and the ΛCDM fixed point used to say the model approaches ΛCDM."}],"review_version":1}