pith. sign in

arxiv: 2604.11821 · v1 · submitted 2026-04-11 · 🌀 gr-qc

Polytropic f(Q) cosmology and its implications for the H₀ tension

Pith reviewed 2026-05-10 16:36 UTC · model grok-4.3

classification 🌀 gr-qc
keywords polytropic equation of statef(Q) gravityH0 tensionBayesian parameter estimationdeceleration parameterstatefinder diagnosticslate-time cosmologyexact solutions
0
0 comments X

The pith

A polytropic equation of state paired with power-law f(Q) gravity yields exact late-time solutions whose Bayesian fits address the H0 tension.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper seeks to model the dark energy component driving cosmic acceleration without presupposing a fixed fluid form. Instead it adopts a polytropic equation of state whose two free parameters are left adjustable and combines this with a recently proposed power-law f(Q) modification of gravity. Exact cosmological solutions are derived from the resulting Friedmann equations. These solutions are then confronted with observational data through a Bayesian analysis that employs the emcee sampler to obtain posterior constraints on all model parameters. From the constrained parameters the authors compute the present-day deceleration parameter, the statefinder trajectory, and the implied value of the Hubble constant, thereby reporting the status of the H0 tension within this framework.

Core claim

By adopting a polytropic equation of state with free parameters in a power-law f(Q) cosmological framework, exact solutions for the scale factor and other quantities are obtained. Rigorous Bayesian parameter estimation from data then yields constraints on the model parameters, allowing interpretation of the deceleration parameter and statefinder diagnostics. The model thereby provides a specific prediction for the current value of H0 and its implications for the observed tension between early- and late-time measurements.

What carries the argument

The polytropic equation of state with free parameters combined with a power-law form of f(Q) in symmetric teleparallel gravity, which permits exact integration of the background equations for late-time dynamics.

If this is right

  • The constrained polytropic index and other parameters furnish concrete predictions for the present deceleration parameter.
  • Statefinder trajectories generated from the posterior can be compared directly with LambdaCDM to test distinguishability.
  • The posterior on H0 quantifies how much of the early-versus-late discrepancy is absorbed by the model.
  • The same framework supplies a ready template for fitting future supernova or BAO datasets without additional tuning.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The polytropic approach could be ported to other modified-gravity settings such as f(T) or f(R) to compare their H0 posteriors.
  • Including early-universe physics within the same polytropic f(Q) ansatz might reveal whether the model can simultaneously satisfy CMB constraints.
  • Future surveys with percent-level H0 precision will tighten the allowed range on the polytropic exponent and thereby test the model's viability.

Load-bearing premise

The power-law assumption for f(Q) together with the polytropic equation of state with free parameters is sufficient to describe the late-time universe without needing extra adjustments or fine-tuning.

What would settle it

An independent, high-precision measurement of the Hubble constant lying well outside the one-sigma posterior interval obtained from the Bayesian fit of this model would falsify the claim that the polytropic f(Q) construction adequately accounts for the H0 tension.

Figures

Figures reproduced from arXiv: 2604.11821 by Raja Solanki.

Figure 1
Figure 1. Figure 1: FIG. 1. The 1 [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The 1 [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The 1 [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Plot showing the normalized expansion rate [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Covariance matrices representing the relationship among the constrained model parameters corresponding to the dataset [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Plot shows the comparison of the deceleration parameter behavior corresponding to the median values from the posterior [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. The effective equation of state parameter behavior corresponding to the median values from the posterior distributions [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The effective dark energy equation of state parameter behavior corresponding to the median values from the posterior [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Behavior of the statefinder parameters corresponding to the median values from the posterior distributions using three [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Plot showing the comparison among the different [PITH_FULL_IMAGE:figures/full_fig_p019_10.png] view at source ↗
read the original abstract

Understanding the late-time cosmic phenomenon of the universe commonly referred to as the dark energy problem, which is one of the prominent tension in the field of theoretical as well as observational cosmology. In this work, we attempt to analyze the nature of the missing fluid of the universe. In order to do so, we employ a poly tropic equation of state consisting of free parameters rather assuming directly a particular form of the fluid. In addition, for the background geometry we consider a $f(Q)$ cosmology exhibiting power-law assumption, which is recently proposed and found to be attractive in the study of late-time cosmology. We find exact cosmological solution along with a rigorous data analysis, utilizing the Bayesian statistics approach and the emcee ensemble sampler, to find the parameter constraints and then we interpret the parameters of physical interests such as deceleration and statefinder parameter. Also, we present status of the $H_0$ tension predicted by our polytropic $f(Q)$ cosmological model.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript derives exact FLRW solutions in f(Q) gravity under a power-law ansatz f(Q) = α(-Q)^n combined with a polytropic equation of state p = K ρ^γ (two free parameters), performs Bayesian parameter estimation via the emcee sampler on cosmological datasets, extracts physical quantities including the deceleration and statefinder parameters, and reports the resulting status of the H0 tension.

Significance. If the exact solutions are free of derivation gaps and the posterior constraints demonstrate stable H0 values that reduce tension without the polytropic index and power-law index absorbing all mismatches via fine-tuning, the work would provide a concrete modified-gravity alternative for late-time acceleration. The use of ensemble sampling and explicit exact solutions aids reproducibility and verifiability.

major comments (2)
  1. [§4] §4 (Bayesian constraints and H0 discussion): the status of the H0 tension is evaluated after fitting the polytropic index, amplitude K, and power-law index n to the same datasets that define the tension; without unfitted predictions, cross-validation on independent benchmarks, or external priors, the reported alleviation is effectively a fitted outcome rather than a model prediction. This is load-bearing for the central claim.
  2. [§3] §3 (exact solutions): the power-law f(Q) fixes the non-metricity contribution at all redshifts, so consistency of the effective w_DE(z) trajectory with data across the full likelihood range must be shown explicitly; any mismatch can be absorbed only by tuning the two free parameters, risking that tension reduction is an artifact of the restricted functional family.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We respond point by point to the major concerns, offering clarifications on the methodology and indicating revisions where they will strengthen the presentation without altering the core results.

read point-by-point responses
  1. Referee: [§4] §4 (Bayesian constraints and H0 discussion): the status of the H0 tension is evaluated after fitting the polytropic index, amplitude K, and power-law index n to the same datasets that define the tension; without unfitted predictions, cross-validation on independent benchmarks, or external priors, the reported alleviation is effectively a fitted outcome rather than a model prediction. This is load-bearing for the central claim.

    Authors: We agree that the reported status of the H0 tension follows from a Bayesian fit of the free parameters (polytropic index, K, and n) to the same datasets used to quantify the tension. This is the standard procedure for testing modified-gravity models against observations: the exact FLRW solutions fix the functional form of the expansion history, after which the posterior constrains the parameters and yields the implied H0 value. The use of multiple independent probes (Pantheon+, BAO, CMB, etc.) already provides an internal cross-check on consistency. We do not claim an unfitted a-priori prediction; rather, the exact solutions demonstrate that a single functional family can simultaneously accommodate a higher H0 while remaining compatible with the full dataset. No additional unfitted forecasts or external priors are required for the claim as stated, and we therefore do not revise the manuscript on this point. revision: no

  2. Referee: [§3] §3 (exact solutions): the power-law f(Q) fixes the non-metricity contribution at all redshifts, so consistency of the effective w_DE(z) trajectory with data across the full likelihood range must be shown explicitly; any mismatch can be absorbed only by tuning the two free parameters, risking that tension reduction is an artifact of the restricted functional family.

    Authors: We concur that an explicit demonstration of the effective dark-energy equation-of-state trajectory w_DE(z) is valuable. The power-law ansatz for f(Q) does fix the non-metricity scalar at every redshift once the scale factor is known, and the polytropic EOS then determines the remaining dynamics. In the revised manuscript we will add a figure (or panel) displaying w_DE(z) evaluated along the posterior samples, together with a brief discussion of its redshift dependence and comparison to the observational constraints. This will make transparent that the model remains physically consistent across the likelihood range and will clarify the extent to which parameter tuning is responsible for the reported H0 values. The choice of power-law f(Q) is retained for its analytic solvability, but the added plot will address the referee’s concern about possible artifacts. revision: partial

Circularity Check

1 steps flagged

H0 tension status obtained directly from Bayesian fits of polytropic f(Q) parameters to the same data

specific steps
  1. fitted input called prediction [Abstract]
    "we present status of the H0 tension predicted by our polytropic f(Q) cosmological model."

    The model parameters (including those controlling late-time acceleration and H0) are fitted to the data using Bayesian statistics and emcee; the H0 tension value is then read off the posterior. This makes the quoted 'prediction' of tension status statistically equivalent to the input fit rather than an independent derivation.

full rationale

The paper derives exact FLRW solutions under the assumed power-law f(Q) and polytropic EOS, then constrains the two free parameters via emcee MCMC on observational data. The reported H0 tension status is extracted from the resulting posterior; no unfitted, out-of-sample prediction or external benchmark is shown. This reduces the central claim to a fitted outcome by construction. No other circular steps (self-citation chains, imported uniqueness theorems, or renamed ansatze) are evident from the provided text.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the power-law ansatz for f(Q), the polytropic form of the equation of state, and the validity of the background FLRW metric in f(Q) gravity; these are introduced without independent derivation in the abstract.

free parameters (2)
  • polytropic index and amplitude
    Two free parameters in the polytropic EOS are left open and fitted to data.
  • power-law index of f(Q)
    The exponent in the assumed power-law form of f(Q) is treated as a free parameter.
axioms (2)
  • domain assumption Power-law form for f(Q) is sufficient to describe late-time cosmology
    Invoked to obtain exact solutions; no justification supplied in abstract.
  • domain assumption Polytropic EOS adequately models the missing fluid
    Chosen instead of a specific dark-energy form; parameters fitted rather than derived.

pith-pipeline@v0.9.0 · 5462 in / 1458 out tokens · 24366 ms · 2026-05-10T16:36:41.103245+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages

  1. [1]

    Cosmic Chronometers The Cosmic Chronometer dataset providesH(z)measurements derived from observations of massive galaxies, hav- ing the redshift interval 0.07≤z≤2.41. These measurements are obtained using the differential age technique, which estimates the Hubble parameter based on the age difference between galaxies at slightly different redshifts [46]. ...

  2. [2]

    SNIa datasets The Pantheon+SH0ES data points sample represents a comprehensive compilations of Type Ia supernovae observa- tions to date, significantly enhancing the statistical power of cosmological analyses based on distance-redshift mea- surements. It extends the original Pantheon sample by incorporating improved calibrations, light-curve fits, and add...

  3. [3]

    CMB We employ the compressed CMB likelihood as the dynamicalf(Q)dark energy model along with the polytropic EoS primarily affect the late-time expansion history of the Universe and mainly alter the geometrical aspects of the CMB. The full CMB temperature polarization power spectrum contains several small scale non-geometric features, such as the low-ℓpowe...

  4. [4]

    BAO (DESI DR2) We consider the Baryon Acoustic Oscillation (BAO) measurements from the DESI Data Release 2 (DR2) sample, which provide one of the most precise geometric probes of the late-time Universe. BAO originate from sound waves propagating in the tightly coupled photon-baryon plasma before recombination, leaving a characteristic comoving scale at th...

  5. [5]

    (48) The absolute difference between the two determinations is, ∆H ΛCDM 0 =73.57−69.45=4.12 km s −1 Mpc−1

    Hubble T ension in theΛCDM Model Using theΛCDM model constraints obtained in Table (II), we have HΛCDM 0 (CC+SN) =73.57 +0.49 −0.36 km s−1 Mpc−1, (47) HΛCDM 0 (BAO+CMB) =69.45 +0.57 −0.48 km s−1 Mpc−1. (48) The absolute difference between the two determinations is, ∆H ΛCDM 0 =73.57−69.45=4.12 km s −1 Mpc−1. (49) The combined uncertainty is, σcomb = q (0.3...

  6. [6]

    The corresponding Hubble constant estimates are, H f(Q) 0 (CC+SN) =72.41 +0.56 −0.60 km s−1 Mpc−1, (52) H f(Q) 0 (BAO+CMB) =68.19 +0.37 −0.34 km s−1 Mpc−1

    Hubble T ension in the Polytropic f(Q)Model We now repeat the same analysis for the polytropicf(Q)cosmology using the constraints reported in Table (I). The corresponding Hubble constant estimates are, H f(Q) 0 (CC+SN) =72.41 +0.56 −0.60 km s−1 Mpc−1, (52) H f(Q) 0 (BAO+CMB) =68.19 +0.37 −0.34 km s−1 Mpc−1. (53) The difference between the late and early t...

  7. [7]

    In fact, the level of internal tension between the CC+SN and BAO+CMB datasets indicates a mild relieve but comparable to that found in theΛCDM model

    Physical Interpretation The above results indicate that, although the polytropicf(Q)model provides an excellent fit to late-time cosmo- logical data and allows for a dynamically effective dark energy sector, it does not fully reconcile the discrepancy between early and late Universe determinations of the Hubble constant. In fact, the level of internal ten...

  8. [8]

    Hubble, M

    E. Hubble, M. L. Humason,Astrophys. J.74, 43 (1931)

  9. [9]

    Astrophys.641, A6 (2020)

    Planck Collaboration,Astron. Astrophys.641, A6 (2020)

  10. [10]

    Aiola et al.,J

    S. Aiola et al.,J. Cosmol. Astropart. Phys.12, 047 (2020)

  11. [11]

    Dutcher et al.,Phys

    D. Dutcher et al.,Phys. Rev. D104, 022003 (2021)

  12. [12]

    E. D. Valentino et al.,Class. Quantum Grav.38, 153001 (2021)

  13. [13]

    E. D. Valentino et al.,Phys. Dark Univ.49, 101965 (2025)

  14. [14]

    A. G. Riess et al.,Astrophys. J.876, 85 (2019)

  15. [15]

    A. G. Riess et al.,Astrophys. J. Lett.908, L6 (2021)

  16. [16]

    J. B. Jim ´enez, L. Heisenberg and T.S. Koivisto,Universe5, 173 (2019)

  17. [17]

    J. M. Nester and H.-J. Yo,Chin. J. Phys.37, 113 (1999)

  18. [18]

    J. B. Jim ´enez, L. Heisenberg, and T. S. Koivisto,JCAP08, 039 (2018)

  19. [19]

    R ¨unkla and O

    M. R ¨unkla and O. Vilson,Phys. Rev. D98, 084034 (2018),

  20. [20]

    J. B. Jim ´enez et al.,Phys. Rev. D101, 103507 (2020)

  21. [21]

    Blixt et al.,Phys

    D. Blixt et al.,Phys. Rev. D109, 044061 (2024)

  22. [22]

    J. B. Jim ´enez and T.S. Koivisto,Int. J. Geom. Meth. Mod. Phys.19, 2250108 (2022)

  23. [23]

    Adak,Int

    M. Adak,Int. J. Geom. Meth. Mod. Phys.15, 1850198 (2018)

  24. [24]

    Tomonari and S

    K. Tomonari and S. Bahamonde,Eur. Phys. J. C84, 349 (2024)

  25. [25]

    D Ambrosio et al.,Phys

    F. D Ambrosio et al.,Phys. Rev. D105, 024042 (2022)

  26. [26]

    Jimenez, L

    J.B. Jimenez, L. Heisenberg and T. Koivisto,Phys. Rev. D98, 044048 (2018)

  27. [27]

    Quiros,Phys

    I. Quiros,Phys. Rev. D105, 104060 (2022)

  28. [28]

    Lazkoz et al.,Phys

    R. Lazkoz et al.,Phys. Rev. D100, 104027 (2019)

  29. [29]

    Khyllep et al.,Phys

    W. Khyllep et al.,Phys. Rev. D107, 044022 (2023)

  30. [30]

    Frusciante,Phys

    N. Frusciante,Phys. Rev. D103, 044021 (2021)

  31. [31]

    F. K. Anagnostopoulos et al.,Eur. Phys. J. C83, 58, (2023)

  32. [32]

    D ´Agostino, R

    R. D ´Agostino, R. C. Nunes,Phys. Rev. D106, 124053 (2022)

  33. [33]

    T. B. Goncalves, L. Atayde, N. Frusciante,Phys. Rev. D109, 084003 (2024)

  34. [34]

    Yang et al.,Mon

    Y. Yang et al.,Mon. Not. Roy. Astron. Soc.533, 2232-2241 (2024)

  35. [35]

    G. N. Gadbail, S. Mandal, P . K. Sahoo,Astrophys. J972, 172 (2024)

  36. [36]

    Nojiri, S

    S. Nojiri, S. D. Odintsov,Phys. Dark Univ.45, 101538 (2024)

  37. [37]

    Sahlu, A

    S. Sahlu, A. d. l. Cruz-Dombriz, A. Abebe,Mon. Not. Roy. Astron. Soc.539, 690-703 (2025)

  38. [38]

    Wang et al.,Astrophys

    Q. Wang et al.,Astrophys. J974, 7 (2024)

  39. [39]

    I. R. Vasquez, A. Oliveros,Gen. Relativ. Grav.57, 67 (2025)

  40. [40]

    Karmakar, S

    P . Karmakar, S. Haridasu, A. Nishizawa,arXiv, arXiv:2508.21054 (2025). 22

  41. [41]

    R. F. Tooper,Astrophys. J.140, 434 (1964)

  42. [42]

    P . H. Chavanis,Eur. Phys. J. Plus130, 130 (2015)

  43. [43]

    Mukhopadhyay and S

    U. Mukhopadhyay and S. Ray,Mod. Phys. Lett. A23, 3187 (2008)

  44. [44]

    Moradpour and M

    H. Moradpour and M. T. M. Sabet,Can. J. Phys.94, 1 (2016)

  45. [45]

    R. H. Lin and X. H. Zhai,Phys. Rev. D103, 124001 (2021)

  46. [46]

    J. C. N. D. Araujo and H. G. M. Fortes,arXiv, arXiv:2407.08884 (2024)

  47. [47]

    Harko,Phys

    T. Harko,Phys. Rev. D98, 084043 (2018)

  48. [48]

    C. G. B ¨ohmer and T. Harko,JCAP2007, 025 (2007)

  49. [49]

    A. Y. Kamenshchik, U. Moschella and V . Pasquier,Phys. Lett. B511, 265 (2001)

  50. [50]

    M. C. Bento, O. Bertolami and A. A. Sen,Phys. Rev. D66, 043507 (2002)

  51. [51]

    Shabani, A

    H. Shabani, A. De, and Tee-How Loo,Eur. Phys. J. C83, 535 (2023)

  52. [52]

    D. F. Mackey et al.,Publ. Astron. Soc. Pac.125, 306(2013)

  53. [53]

    H. Yu, B. Ratra, and F.-Y. Wang,Astrophys. J.856, 3 (2018)

  54. [54]

    Solanki et al.,Phys

    R. Solanki et al.,Phys. Dark Univ.32, 100820 (2021)

  55. [55]

    Brout et al.,Astrophys

    D. Brout et al.,Astrophys. J.938, 110 (2022)

  56. [56]

    Eric Aubourg et al.,Phys. Rev. D92, 123516 (2015)

  57. [57]

    Astrophys.607, A95 (2017)

    Planck Collaboration,Astron. Astrophys.607, A95 (2017)

  58. [58]

    L. A. Escamilla et al.,Phys.Dark Univ.46, 101652 (2024)

  59. [59]

    Y. Wang, P . Mukherjee,Phys. Rev. D76, 103533 (2007)

  60. [60]

    DESI Collaboration,Phys. Rev. D112, 083515 (2025)

  61. [61]

    Betoule et al.,Astron

    M. Betoule et al.,Astron. Astrophys.568, A22 (2014)

  62. [62]

    Sahni et al.,JETP Letters77, 201-206 (2003)

    V . Sahni et al.,JETP Letters77, 201-206 (2003)

  63. [63]

    Alam et al.,Mon

    U. Alam et al.,Mon. Not. R. Astron. Soc.344, 1057-1074 (2003)