REVIEW 2 major objections 4 minor 8 references
Analytic Integration of the Lambert W Cosmic Fluid Model H(z) Formula
T0 review · 2 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The paper shows that the Lambert W cosmic fluid model's ln(H(z1)/H0) equals exactly the sum of equations (31), (32), and (33), so H/H0 is the exponential of that sum.
desk verdict Correct, modest, honestly presented analytic integration for the Lambert W cosmology, conditional on the correctness of Dubey et al.'s source integral and missing its verification figure. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing device is the change of variable y=W(x), where y=W(x) means x=y e^y, applied to logarithmic integrals of the form ∫ f(W(x)) d[ln x]. For a positive power p, the substitution turns ∫ W(x)^p d[ln x] into ∫ (y^p+$y^{{p-1}}$) dy. For ln W(x), it turns the measure d[ln x] into d[ln y]+dy, so the antiderivative becomes (ln W(x))^2/2 + W(x) ln W(x) - W(x). These two antiderivatives, together with ∫ d[ln x]=ln x, assemble into the three terms (31), (32), and (33) that the paper sums to evaluate ln(H(z1)/H0).
What would settle it
Compute the sum (31)+(32)+(33) for θ1=0.087, θ2=-3.36 at a grid of redshifts z1 in (0,2) and compare it with independent high-precision numerical quadrature of equation (7). Any difference beyond the quadrature tolerance would falsify the analytic formula; a secondary check is to differentiate the closed form and compare its derivative with the integrand of (7).
Extended reading notes
Core claim
The paper establishes that the integral defining ln(H(z1)/H0), equations (7) and (8), has a closed-form antiderivative in terms of the Lambert W function. Writing y=W(a), so that a=y e^y, turns the logarithmic measure d[ln a] into (1+1/y) dy, which makes ∫ W(a)^p d[ln a] equal to ∫ (y^p+$y^{{p-1}}$) dy and makes ∫ ln W(a) d[ln a] equal to ∫ u du + ∫ ln y dy with u=ln y. Evaluating the resulting antiderivatives from a1 to a0=1 gives the three terms (31), (32), and (33), and H(z1)/H0 is the exponential of their sum. The identity W(1)=Ω, with ln Ω=-Ω, compacts the present-time boundary terms and produces partial cancellations.
Load-bearing premise
The derivation takes the starting integral, equations (7) and (8), from the fitting paper [4] and does not rederive it from the Lambert W equation of state; if that integral is not the correct Friedmann consequence of the model, or if the flat-universe assumption fails, the closed form inherits the error. It also assumes the principal real branch of W(a) for a in (0,1], with W(1)=Ω.
Editorial extensions
If this is right
- Numerical quadrature of equations (7) and (8) is no longer required to get H(z1)/H0; the exponential of the sum (31)+(32)+(33) gives the same value in closed form.
- The formula separates cleanly into a θ1 term, a θ2 term, and a constant term, so the dependence of the Hubble factor on each parameter of the Lambert W model is visible directly.
- The closed form can be used as a drop-in replacement in fits of the model to Hubble data, removing one source of numerical error from the comparison with the recent fit.
- Derivatives of ln(H/H0) with respect to z1 and with respect to θ1 and θ2 become analytic, which should make parameter exploration faster and more stable.
Reading between the lines
- The same substitution y=W(x) works for any positive power p in ∫ W(x)^p d[ln x], so the method would extend immediately to equation-of-state variants that use other integer powers of W(a).
- Differentiating the closed form would give an analytic deceleration parameter q(z), including the redshift where the universe switches from deceleration to acceleration; the paper notes a future working paper on that transition.
- A re-fit of the model using the exact formula instead of numerical integration could shift the fitted θ values slightly, because the analytic form removes quadrature error that grows with redshift.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This working paper treats the Lambert W cosmic-fluid model of Saha and Bamba. It takes as its starting point the integral formula for ln(H(z1)/H0) that appears in a recent fitting paper by Dubey et al. [4], reproduced here as eqs. (7) and (8). The core of the paper is an analytic evaluation of that integral: the substitution y = W(a), with d ln a = (1 + 1/y) dy, is used to compute the three logarithmic antiderivatives I3(a), L(a), and K(a). These are combined into the definite integral from a1 to 1, giving the closed-form expression in eqs. (31)-(33), with the exponential yielding H(z)/H0. A Python code snippet for the formula is included in Section 3, where the text promises a graphical comparison with numerical integration and cosmic-chronometer data, although no figure appears in the manuscript.
Significance. If the integral in eqs. (7)-(8) is accepted as correct, the derivation is sound and the closed form is genuinely useful: it removes the need for numerical integration of that particular integral, and the substitution-based derivation is transparent and self-contained. The paper performs no fitting and introduces no free parameters; the θ values are taken from [4]. The derivation of I3 and L is correct, including the evaluation at W(1) = Ω and the use of ln Ω = -Ω. The main limitations are that the starting integrand is an unverified external input from [4], and the promised numerical/visual verification is absent from the manuscript. The novelty is modest, but the result is a clean technical contribution if the external integral is correct.
major comments (2)
- [§3] The verification promised at the start of Section 3 is not present in the manuscript. The text refers to 'the following figure' and describes a comparison with figure 3 of [4], including orange dots from numerical integration, green dots from the analytic formula, and blue ΛCDM curves, but no figure is included anywhere in the paper. Because the conclusion states that the analytic formula agrees with the numerical integration, this missing figure is load-bearing evidence for the verification claim. Please include the figure with a proper caption and legend, or explicitly state that the numerical verification is deferred to a later version.
- [§2, eqs. (7)-(8)] The integral that is integrated analytically is adapted from eqn. (14) of [4] and is not derived in this paper. The final H/H0 formula inherits every potential error in that integrand: the factor 3/2, the '+1' inside the bracketed integrand, and the orientation of the integration limits. The manuscript should either derive eqs. (7)-(8) from the Friedmann equation together with the equation of state (1), or state unambiguously that the analytic evaluation is conditional on the correctness of the Dubey et al. integral. This is the only load-bearing external assumption, but it is essential because the closed form in eqs. (31)-(33) is exact for the given integrand and not for a corrected version of it.
minor comments (4)
- [§2, eqs. (7)-(8)] The equality of the redshift integral and the scale-factor integral is asserted rather than shown; adding the one-line derivation using a = 1/(1+z) and da/a = -dz/(1+z) would make the orientation of the limits explicit and remove any ambiguity.
- [§2, eqs. (20)-(25)] The notation d[u] + d[y] in eq. (21) is informal. Since y = e^u, the separation d ln x = du + dy is exact, and stating this before splitting the integral would make the computation easier to follow.
- [§3] The Python code snippet is not self-contained: theta1, theta2, log, and exp are used without definitions or imports, and lambertw is only mentioned in a note. Please provide a complete runnable function, including the scipy.special import and the parameter values as arguments or explicit globals.
- [§2, eq. (30)] The symbol Ω for W(1) should be introduced as the Omega constant satisfying Ω e^Ω = 1, so that it is not confused with the matter density parameter Ωm used later in Section 3.
Circularity Check
No significant circularity: the analytic H(z) formula is an antiderivative evaluation of an explicitly attributed integral, with no fitted parameters and no load-bearing self-citation.
full rationale
The paper's claimed result, the sum of equations (31), (32), and (33), is derived from the starting integrals (7) and (8), which the paper explicitly states are 'adapted from eqn (14) of [4]' — an external, attributed formula. The derivation then proceeds by the substitution y=W(a), giving d[ln(a)] = (1+1/y)dy, and evaluates the three logarithmic integrals I3(a), L(a), and K(a) in closed form. This is a self-contained antiderivative computation; the target quantity ln(H/H0) is not used to define any input or parameter. The parameters θ1 and θ2 are taken from Dubey et al. [4] as inputs and appear linearly; they are not fitted to the claimed formula. The only load-bearing step external to the paper is the correctness of the Friedmann-equation consequence (7)/(8) from [4], but an inherited assumption is a dependency, not circularity: if the integral is mis-derived, the analytic expression is exactly the antiderivative of the wrong integrand, which is a correctness concern rather than a circularity. No load-bearing self-citation appears; the current authors have no overlap with [4], and reference [8] (Mező's book) supplies background techniques rather than the central result. There is no uniqueness theorem invoked, no ansatz smuggled in via citation, and no renaming of a known empirical pattern. Even the described-but-not-shown graphical verification would affect empirical validation, not circularity of the derivation. Therefore no circular step is present; score 0.
Assumptions & free parameters
assumptions (3)
- standard math Principal real branch of the Lambert W function and its properties, including W(x) exp(W(x)) = x, W(1)=Omega, and ln Omega = -Omega.
- domain assumption The integral formula for ln(H/H0) (equations 7 and 8) correctly follows from the Lambert W equation of state (equation 1) and the Friedmann equations.
- domain assumption The principal real branch of W is the physical branch for a in (0,1].
Cite this review
Pith. "Pith review of Analytic Integration of the Lambert W Cosmic Fluid Model H(z) Formula." pith.science (2026). https://pith.science/paper/XOOONUHE
@misc{pith2026260807878,
author = {Pith},
title = {Pith review of: Analytic Integration of the Lambert W Cosmic Fluid Model H(z) Formula},
year = {2026},
howpublished = {\url{https://pith.science/paper/XOOONUHE}},
note = {Machine review of arXiv:2608.07878}
}
abstract
The Lambert $W$ (LW) model for the cosmic fluid equation of state was proposed by S. Saha and K. Bamba in 2019-2020. A recent (early 2026) paper by Dubey, et al, carries out a new fit of the LW model to observational data, in order to estimate the model's parameters. That paper exhibits a formula for the logarithm of the relative Hubble factor at redshift $z$, $\ln(H(z)/H_0)$, which is based upon a numerical integration. That integral can be evaluated analytically, and we present the details in this working paper. The resulting analytic expression for $H(z)/H_0$ may be convenient for exploration of the LW model.
Reference graph
Works this paper leans on
-
[4]
The LambertWequation of state in light of DESI BAO
Vipin Chandra Dubey, Subhajit Saha and Abdulla Al Mamon, “The LambertWequation of state in light of DESI BAO”.New Astronomy, 2026, vol 128, art 102616. doi:10.1016/j.newast.2026.102616. Preprint available as arXiv:2601.20972v1
-
[1]
The Lambert $W$ function: A newcomer in the Cosmology class?
Subhajit Saha and Kazuharu Bamba, “The Lambert W Function: A Newcomer in the Cosmology Class?”,Z. Naturforsch., 2020, 75(1)a: pp 23-27. doi:10.1515/zna-2019-0240. Preprint available as arXiv:1908.11712v2
work page Pith review arXiv 2020
-
[2]
Growth of Perturbations using Lambert$W$ Equation of State
Manisha Banerjee, Sudipta Das, Abdulla Al Mamon, Subhajit Saha and Kazuharu Bamba, “Growth of Perturbations using Lam- bert W Equation of State”,Int. J. Geom. Meth. Mod. Phys., 2021, vol 18, no 09, art 2150139. doi:10.1142/S0219887821501395. Pre- print available as arXiv:2011.12058v1
work page Pith review arXiv 2021
-
[3]
Testing lam- bert W equation of state with observational hubble para- meter data
Abdulla Al Mamon and Subhajit Saha, “Testing lam- bert W equation of state with observational hubble para- meter data”,New Astronomy, 2021, vol 86, art 101567. doi:10.1016/j.newast.2020.101567. Preprint available as arXiv:2005.14061v2
arXiv 2021
-
[5]
Arianna Favale, Adri` a G´ omez-Valent, and Marina Migliac- cio, “Cosmic chronometers to calibrate the ladders and meas- ure the curvature of the Universe. A model-independent study,”Mon. Not. Roy. Astron. Soc., 2023, vol 523, pp 3406-3422. doi:10.1093/mnras/stad1621. Preprint available as arXiv:2301.09591v3
arXiv 2023
-
[6]
Elena Tomasetti, M. Moresco, N. Borghi, K. Jiao, and 5 other authors, “A new measurement of the expansion history of the Universe atz= 1.26 with cosmic chronometers in V ANDELS”,Astronomy and Astrophysics, 2023, vol 679, num- ber A96. doi:10.1051/0004-6361/202346992. Preprint available as arXiv:2305.16387v1. ASBH-20260807 Page 12 07-Aug-2026
arXiv 2023
-
[7]
Arianna Favale, Adri` a G´ omez-Valent, and Marina Migliaccio, “Revisiting model-independent constraints on spatial curvature and cosmic ladders calibration: updated and forecast ana- lyses,”Mon. Not. Roy. Astron. Soc., 2026, vol 546, no 4, art stag303. doi:10.1093/mnras/stag303. Preprint available as arXiv:2511.19332v2
arXiv 2026
-
[8]
[end] ASBH-20260807 Page 13 07-Aug-2026
Istv´ an Mez˝ o,The Lambert W Function: Its Generalizations and Applications, CRC Press, 2022. [end] ASBH-20260807 Page 13 07-Aug-2026
work page 2022
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.