{"id":"955ad44f-6b07-4673-88dd-4f183cccf664","arxiv_id":"2508.15187","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"The interacting polytropic dark energy model is claimed to fit late-time cosmological data, but the central analytic solution for the matter density is algebraically inconsistent with the model's continuity equation.","lead":"A dark energy model with a polytropic equation of state and matter interaction is fitted to Hubble, BAO, DESI, and supernova data, with reported parameters favoring energy flow from dark energy to matter. The central matter-density solution does not satisfy the model's own continuity equation, so the fitted constraints and the viability claim are unsupported.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Matter density solution Eq. (17) violates continuity Eq. (12) except for unphysical ρd0=-(1+K)ρm0; H(z) Eq. (21) and all constraints rest on this error.","rationale":"The single most load-bearing issue is the incorrect matter density solution. The paper's own equations show that Eq. (17) cannot satisfy Eq. (12) unless an unphysical relation between ρd0, ρm0, and K holds. The exact solution includes an additional term that is omitted. Because the Hubble parameter, the MCMC likelihoods, and all derived cosmological diagnostics depend directly on Eq. (21), the central claim of a viable interacting polytropic dark energy model is unsupported. The reader identified precisely this flaw, and the rejection is warranted. The concern is not about external consensus or stylistic issues; it is a mechanical mathematical inconsistency within the paper's derivation.","tokens_in":13954,"tokens_out":8504,"duration_ms":76101,"concrete_test":"Symbolically solve Eq. (15) with integrating factor (1+z)^{-3} and initial condition ρm(0)=ρm0. Verify the exact solution is ρm = (1+z)^3[ρm0 + (ρd0/(1+K))(1-g^η)] with g = (1+K)(1+z)^{-3} - K, and substitute it into Eq. (12) to confirm it satisfies the ODE. Then substitute the paper's Eq. (17) and evaluate the residual at z=0 and z=1 using the best-fit parameters from Table 1 (e.g., K=-0.93, η=-0.22, Ωd0/Ωm0≈2.7). The residual is nonzero and of order ρm0, demonstrating that Eq. (17) is not a solution.","verdict_should_be":"REJECT","load_bearing_attack":"The model's central H(z) (Eq. 21) is built on the matter density solution Eq. (17) and dark energy solution Eq. (18). Substituting Eq. (17) into the modified conservation equation (Eq. 12) gives, after canceling common factors, the condition -3η(1+K)ρm0 g^{η-1} = 3ηρd0 g^{η-1}, i.e., ρd0 = -(1+K)ρm0, where g = (1+K)(1+z)^{-3} - K. This is not satisfied by the best-fit parameters: for the Hubble77+BAO26 fit, K ≈ -0.93, so 1+K ≈ 0.07, while ρd0/ρm0 = Ωd0/Ωm0 ≈ 2.7; the relation fails by roughly a factor of 40 in magnitude. The general solution of the linear ODE Eq. (15), obtained with the integrating factor (1+z)^{-3}, is ρm = (1+z)^3 [ρm0 + (ρd0/(1+K))(1 - g^η)]. The paper's Eq. (17) drops the term (ρd0/(1+K))(1 - g^η), which is non-negligible at z>0 (e.g., at z=1 with the first dataset it is a substantial fraction of ρm). Consequently Eq. (21) is not the Hubble parameter of the model for the fitted parameters; the MCMC constraints, ω_eff, q, statefinder, and age diagnostics all inherit this inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a spatially flat FLRW universe containing pressureless matter and a polytropic dark energy fluid with equation of state p_d = α ρ_d^{1+1/β}, coupled through Q = 3ηHρ_d. After imposing β = 1 - η, the authors derive matter and dark energy densities (Eqs. 17 and 18), form the Hubble parameter (Eq. 21), and fit the parameters H0, Ωd0, η, and K to three joint data combinations (Hubble77+BAO26, Hubble77+Pantheon+, Hubble77+BAO26+DESI DR2) via MCMC. They report best-fit values, show plots of effective energy density, pressure, equation-of-state, deceleration, statefinder, age, and the interaction term, and conclude that the interacting polytropic dark energy model is a viable candidate for late-time cosmic acceleration.","tokens_in":14305,"tokens_out":6316,"duration_ms":68964,"significance":"If the derivation were correct, the paper would provide useful constraints on an analytically tractable interacting polytropic dark energy model, and the use of multiple public cosmological datasets is appropriate. However, the central matter-density solution is not a solution of the modified conservation equation, so the Hubble parameter used in the likelihood is not the model's true Hubble parameter. All parameter constraints and diagnostic conclusions inherit this error. The paper does not currently support its viability claim.","major_comments":[{"comment":"Eq. (17) does not satisfy the modified conservation equation Eq. (12). Setting g(z) = (1+K)(1+z)^{-3} - K, substituting ρm = ρm0(1+z)^3 g^η and ρd = ρd0 g^{η-1} into Eq. (12) gives -3η(1+K)ρm0 g^{η-1} on the left and 3ηρd0 g^{η-1} on the right, so consistency requires ρd0 = -(1+K)ρm0. This is incompatible with the best fits in Table 1: with Ωd0 ≈ 0.73, ρd0/ρm0 ≈ 2.7, and K ≈ -0.93, so 1+K ≈ 0.07, the relation fails by roughly a factor of 40. Direct integration of Eq. (15) gives ρm = (1+z)^3 [ρm0 + (ρd0/(1+K))(1 - g^η)], with a term that Eq. (17) drops. Thus Eq. (21) is not the Hubble parameter of the interacting model, and the MCMC constraints and all derived diagnostics are invalid. Note also that Eq. (16) is miswritten: d/dz(ρm/(1+z)^3) equals the integrand, not its integral.","section":"Sec. 3, Eqs. (12), (17), (21)"},{"comment":"The diagnostic results are presented as validations, but they are deterministic outputs of Eq. (21) evaluated at the same best-fit parameters used for the fit. The effective EoS, deceleration parameter, density parameters, statefinder pair, and age are algebraic consequences of the fitted H(z); Figures 5-10 and Table 2 therefore restate the fit rather than independently testing the model. These should be presented as derived predictions with appropriate caveats, not as observational confirmation.","section":"Secs. 5-8"},{"comment":"The BAO likelihood is ambiguous and appears to double-count DESI DR2. The text says the 26 BAO measurements include the most recent DESI DR2 results (Refs. [34-37]), yet one of the joint datasets is labeled 'Hubble77+BAO26+DESI DR2'. If the 26-point BAO sample already contains the DESI DR2 points, then the third dataset enters those points twice, which would bias the reported constraints and make the comparison across dataset combinations unreliable.","section":"Sec. 4, BAO likelihood"}],"minor_comments":[{"comment":"As noted above, Eq. (16) incorrectly writes d/dz(ρm/(1+z)^3) as an integral; the integral sign on the right should be removed. This typo accompanies the more serious integration error in Eq. (17).","section":"Eq. (16)"},{"comment":"The interaction term is written as Q = 3ηHρ_d in the abstract and throughout the derivation, but Sec. 9 states Q = ηHρ_d. This inconsistency should be fixed.","section":"Abstract and Sec. 9"},{"comment":"The equation-of-state parameter is written ambiguously as ω_eff = ρ_eff / P_eff rather than ω_eff = P_eff / ρ_eff. The later formulas use the standard convention, so the displayed relation should be corrected.","section":"Eq. (7)"},{"comment":"The pressure expression in Eq. (34) is typeset in a way that obscures its derivation and contains apparent dimensional inconsistencies. A clean, dimensionally correct expression should be provided.","section":"Eq. (34)"},{"comment":"The condition β = 1 - η is imposed to make the integrals tractable; this is an additional model assumption that reduces the generality of the polytropic EoS. It should be stated explicitly as such in the model setup, not introduced as a mere integration trick.","section":"General"}],"recommendation":"reject","confidential_remarks":"The central derivation error is load-bearing and cannot be fixed by rewriting; the fitted Hubble parameter is not that of the stated model. The diagnostic sections would also need reframing even after a correct derivation. I see no path to acceptance within the current scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's core analytic result—the matter density solution (Eq. 17)—is not a solution of the model's own conservation equation (Eq. 12). Substituting it forces ρ_d0 = -(1+K)ρ_m0, which the best-fit parameters violate by a factor of ~40. So Eq. 21, the MCMC constraints, and every diagnostic built on them are invalid.\n\nWhat the paper does well: it is clearly organized, the interacting polytropic setup is standard, and the MCMC pipeline with three data combinations (including DESI DR2) follows normal practice. The special case β = 1-η is a legitimate simplifying choice, and the writing is straightforward. But none of that rescues the central derivation.\n\nThe problem occurs at Eq. (16): the derivative of ρ_m/(1+z)^3 is set equal to an integral, and the integration constant is then dropped. The actual solution is ρ_m = (1+z)^3[ρ_m0 + (ρ_d0/(1+K))(1 - g^η)], not their Eq. (17). The omitted term is non-negligible at z>0 (at z=1 it is a substantial fraction of ρ_m). This is not a minor slip; it changes the Hubble parameter. Additional issues: the 'predictions' in Sections 5-7 are deterministic outputs of the fitted parameters, so they are consistency checks at best, not independent tests. The β=1-η restriction is not physically motivated, and derived quantities lack error bars. These secondary problems would matter if the derivation were sound.\n\nWho is this for? A reader wanting a worked example of an interacting dark energy fit might skim it, but the core result is wrong. The citation pattern is unremarkable—some self-citations, but they are relevant. If you referee it, the substitution check alone is enough for rejection. I would not send it to peer review.","headline":"The paper's central matter-density solution violates its own conservation equation, invalidating the Hubble parameter and all downstream constraints.","tokens_in":14848,"tokens_out":5373,"would_cite":false,"duration_ms":48509,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that an interacting polytropic dark energy fluid, with a linear coupling to dark matter, fits the combined cosmic expansion data and serves as a viable alternative to ΛCDM for late-time acceleration.","keywords":["interacting dark energy","polytropic equation of state","Hubble parameter constraints","MCMC parameter estimation","quintessence","statefinder diagnostic","DESI DR2 BAO","Pantheon+ supernovae"],"falsifier":"Check the model's own matter conservation equation, Eq. (12), at z = 0 using the published best-fit values. Consistency requires (1 + K)ρ_m0 = −ρ_d0; the fits give (1 + K) ≈ 0.07 (since K ≈ −0.93) and ρ_d0/ρ_m0 = Ω_d0/Ω_m0 ≈ 0.73/0.27 ≈ 2.7, so the identity fails by a large factor. The same check at several redshifts in 0 < z < 2 would settle whether the derived H(z) actually satisfies the coupled equations.","tokens_in":13809,"feed_emoji":"🌌","tokens_out":8954,"duration_ms":100611,"temperature":0.7,"pith_summary":"The paper tries to establish that dark energy can be described as a polytropic fluid—a fluid whose pressure is a power of its own density—and that when this fluid interacts linearly with dark matter, the resulting model fits the observed expansion history. The authors derive an analytic Hubble parameter under the condition β = 1 − η, then fit its free parameters to joint Hubble, BAO, Pantheon+, and DESI DR2 data with Markov Chain Monte Carlo. Their best fits give H0 ≈ 67.8–69.2 km/s/Mpc, Ω_d0 ≈ 0.73, and a negative coupling η, meaning energy currently flows from dark energy into dark matter. With these values the model reproduces ΛCDM-like expansion, yields a quintessence equation of state, shifts the deceleration parameter from positive to negative around z ≈ 0.77, and gives a present age near 14 Gyr. A sympathetic reader would care because the model offers a single fluid, rather than a cosmological constant, as the engine of late-time acceleration, and it can be tested with the same joint datasets.","feed_headline":"Interacting polytropic dark energy fits Hubble, BAO, and DESI data","feed_subtitle":"MCMC fits give Ω_d0 ≈ 0.73, negative coupling η, and a quintessence equation of state today.","key_machinery":"The machinery is the polytropic equation of state p_d = α ρ_d^{1+1/β} together with the linear interaction Q = 3ηHρ_d. The crucial step is the substitution β = 1 − η, which turns the coupled conservation equations into closed-form expressions for ρ_m(z), ρ_d(z), and hence the Hubble parameter H(z) of Eq. (21). That single formula carries the whole analysis: it is fed into chi-square fits against Hubble, BAO, Pantheon+, and DESI DR2 data, and then differentiated to yield the effective pressure, equation-of-state parameter, deceleration parameter, statefinder pair, and cosmic age.","core_discovery":"The paper's central claim is that a dark-energy fluid obeying the polytropic equation of state p_d = α ρ_d^{1+1/β}, coupled to pressureless matter by Q = 3ηHρ_d, yields a closed-form Hubble parameter (Eq. 21) that fits the combined observational data. The best-fit parameters are H0 ≈ 67.8–69.2 km/s/Mpc, Ω_d0 ≈ 0.73, K ≈ −0.9, and η < 0, with the precise η depending on the dataset. With these values the model reproduces the expansion history of ΛCDM, keeps the energy density positive and pressure negative, places the effective equation of state in the quintessence range, and predicts a transition from deceleration to acceleration near z ≈ 0.77 and a present age near 14 Gyr. The paper therefor","pith_inferences":["A direct follow-up would be to compute the linear growth rate of matter fluctuations from Eq. (21); because the fitted coupling is negative, the growth should deviate from ΛCDM in a way that galaxy clustering surveys could detect independently of the expansion history.","Since K is defined through α, ρ_d0, and η, the best-fit K implies a specific polytropic constant α; checking that implied constant against independent constraints on the dark-energy equation of state would connect the statistical fit to the fluid interpretation.","The analytic frame could be extended to a redshift-dependent coupling η(z) or a different polytropic index; the statefinder trajectories approaching the ΛCDM fixed point indicate that late-time, low-redshift data will be where such extensions are most easily distinguished."],"forward_implications":["If the central claim holds, late-time acceleration can be driven by an interacting polytropic fluid with no cosmological constant, and the coincidence between ρ_d and ρ_m can be adjusted by the coupling strength η.","The negative η obtained from all three data combinations means energy currently flows from dark energy into dark matter; this direction of transfer is itself testable through galaxy growth-rate measurements.","The predicted present deceleration parameter q0 ≈ −0.42 to −0.50 and effective equation of state ω_eff ≈ −0.61 to −0.67 put the model firmly in the quintessence class, distinguishing it from phantom dark-energy candidates.","Because H0 comes out near 68 km/s/Mpc, the model tracks Planck's H0 rather than the higher local value; this makes it a candidate that preserves CMB-calibrated expansion while fitting lower-redshift probes."],"supporting_citations":[{"why":"Supplies the Planck 2018 H0 baseline against which the fitted H0 values are compared.","marker":"[8]"},{"why":"Investigates the polytropic gas model using cosmological observations, providing the direct predecessor for the present model.","marker":"[24]"},{"why":"Explores polytropic dark energy for late-time cosmic acceleration through an observational approach, motivating the equation of state used here.","marker":"[25]"},{"why":"Compiles the 77 Hubble parameter measurements and the Pantheon+ 1701-supernova sample used in the chi-square fits.","marker":"[32]"},{"why":"Provides the 26 BAO data points and the DESI DR2 BAO list used in the joint observational constraints.","marker":"[33]"},{"why":"Supplies the DESI DR2 BAO measurements included in the third joint dataset combination.","marker":"[34]"},{"why":"Documents energy transfer from dark energy to matter, cited to support the negative η result.","marker":"[38]"}],"fun_headline_variants":["Polytropic dark energy fits DESI+BAO+Hubble data","Negative coupling in polytropic dark energy passes DESI","Polytropic dark energy yields quintessence and 14 Gyr age","Interacting polytropic dark energy matches ΛCDM expansion","DESI data support interacting polytropic dark energy"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The derivation assumes the density formulas it writes solve the coupled conservation equations at every redshift with nothing left over from the integration; if a leftover term is present, the fitted Hubble law is not actually a solution of the model.","fun_headline_variants_meta":{"raw":{"variants":["Polytropic dark energy fits DESI+BAO+Hubble data","Negative coupling in polytropic dark energy passes DESI","Polytropic dark energy yields quintessence and 14 Gyr age","Interacting polytropic dark energy matches ΛCDM expansion","DESI data support interacting polytropic dark energy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000533,"raw_usage":{"total_tokens":2516,"prompt_tokens":977,"completion_tokens":1539,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":721,"completion_tokens_details":{"reasoning_tokens":1451}},"tokens_in":721,"tokens_out":1539,"duration_ms":12858,"temperature":1.0,"reasoning_tokens":1451,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T18:04:48.221313+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Check the model's own matter conservation equation, Eq. (12), at z = 0 using the published best-fit values. Consistency requires (1 + K)ρ_m0 = −ρ_d0; the fits give (1 + K) ≈ 0.07 (since K ≈ −0.93) and ρ_d0/ρ_m0 = Ω_d0/Ω_m0 ≈ 0.73/0.27 ≈ 2.7, so the identity fails by a large factor. The same check at several redshifts in 0 < z < 2 would settle whether the derived H(z) actually satisfies the coupled equations.","supporting_citations":[{"cited_title":"Astrophys","cited_arxiv_id":null,"evidence_quote":"Supplies the Planck 2018 H0 baseline against which the fitted H0 values are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Investigates the polytropic gas model using cosmological observations, providing the direct predecessor for the present model."},{"cited_title":"C, 48, 115110 (2024)","cited_arxiv_id":null,"evidence_quote":"Explores polytropic dark energy for late-time cosmic acceleration through an observational approach, motivating the equation of state used here."},{"cited_title":"Singh, H","cited_arxiv_id":null,"evidence_quote":"Compiles the 77 Hubble parameter measurements and the Pantheon+ 1701-supernova sample used in the chi-square fits."},{"cited_title":"Samaddar, S","cited_arxiv_id":null,"evidence_quote":"Provides the 26 BAO data points and the DESI DR2 BAO list used in the joint observational constraints."},{"cited_title":"Roy, Gen","cited_arxiv_id":null,"evidence_quote":"Documents energy transfer from dark energy to matter, cited to support the negative η result."}],"review_version":1}