{"id":"deb777e6-36eb-4d66-a4f6-a4a9fa3c7174","arxiv_id":"2608.07878","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"An exact analytic formula for ln(H(z)/H0) in the Lambert W cosmic fluid model, obtained by evaluating log-integral antiderivatives of the Lambert W function.","lead":"This paper derives a closed-form analytic expression for the Hubble factor ratio H(z)/H0 in the Lambert W cosmic fluid model, replacing the numerical integration used in a 2026 fitting paper. A generalist should read it because it turns a numerical step in analyzing a cosmological equation of state into an exact formula, making model exploration faster and more transparent.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The antiderivative algebra is internally correct, but the closed form inherits the unverified Dubey et al. integral (7)/(8); if that Friedmann-derived formula is mistaken, the analytic H/H0 is the exact solution of the wrong integral.","rationale":"I verified the internal integration steps: eqs (13)-(16), (20)-(25), and the Ω simplifications all check out, so the antiderivative of the claimed integrand is correct. The remaining risk is the provenance of eqs (7)-(8) from [4]: if those are not the correct Friedmann consequence of eq (1), the analytic sum of (31)-(33) is an exact integral of the wrong expression. The reader's weakest_assumption names the same point. The missing figure in §3 is a presentation defect rather than a mathematical one; it weakens the empirical check but does not alter the derivation. Since the reader already issued CONDITIONAL and my analysis does not move the verdict, UNCHANGED is appropriate.","tokens_in":6123,"tokens_out":7719,"duration_ms":85499,"concrete_test":"Derive eq (8) independently from the flat Friedmann equation and the Lambert W EoS (1): use H² ∝ ρ, d ln ρ/d ln a = -3[1+θ1 ln W(a)+θ2 W(a)³], and integrate from a=1 to a=a1. Check that the result matches eq (8) exactly, including the 3/2 factor, the '+1' term, and the direction of integration. If it differs in any factor or sign, the analytic formula (31)-(33) does not represent H/H0 for the LW model. A quick consistency check is θ1=θ2=0, for which the formula should reduce to (1+z1)^(3/2).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytic derivation in §2 is correct: substituting y=W(a), d ln a = (1+1/y)dy, gives L(y) and I3(y) as stated, and the constants W(1)=Ω, ln Ω=-Ω are handled properly. The load-bearing assumption is external: eqs (7) and (8), taken from eqn (14) of Dubey et al. [4], are the purported Friedmann-equation consequence of the EoS (1). The paper never derives them. Specifically, the factor 3/2, the '+1' inside the integrand, and the orientation of the limits a1→1 must follow from d ln ρ/d ln a = -3(1+w(a)) and H² ∝ ρ. If [4]'s formula has a missing term, a wrong factor, or reversed bounds, then the sum of (31)-(33) is exactly the antiderivative of the wrong integrand. The graphical check described in §3 does not appear in the manuscript, so the claimed agreement with numerical integration is unverifiable from the text. Thus the correctness of the central claim is conditional on [4].","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":6367,"tokens_out":5427,"duration_ms":62439,"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":[{"comment":"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.","section":"§3"},{"comment":"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.","section":"§2, eqs. (7)-(8)"}],"minor_comments":[{"comment":"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.","section":"§2, eqs. (7)-(8)"},{"comment":"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.","section":"§2, eqs. (20)-(25)"},{"comment":"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.","section":"§3"},{"comment":"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.","section":"§2, eq. (30)"}],"recommendation":"major_revision","confidential_remarks":"This is a short technical note rather than a full research article. The central derivation is correct, but the missing verification figure and the unstated dependence on an externally supplied integrand are both significant for a journal reading. If the journal publishes concise technical notes, the paper could be suitable after the requested revisions; otherwise the editor may wish to consider whether the scope and length match the journal's usual contributions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does exactly what it says: it evaluates the ln(H/H0) integral from Dubey et al. in closed form, and the math is right. The substitution y=W(x), the handling of d[ln x]=(1+1/y)dy, the antiderivatives for W^3 and ln W, and the simplification at W(1)=Omega all check out. For this specific Lambert W cosmic fluid model, the resulting analytic formula appears to be new; the cited literature used numerical integration. That is a real, if modest, contribution.\n\nSoft spots, in proportion: the source integral, equations (7)-(8), is taken on faith from Dubey et al. [4]. The paper never derives it from the Friedmann equations, so if the 3/2 factor, the '+1', or the orientation of limits in [4] is wrong, this formula is the exact antiderivative of the wrong integrand. The authors should state this conditional dependence explicitly rather than leave it implicit. Second, the verification figure promised in Section 3 does not appear in the manuscript, so the claimed agreement with numerical integration is currently unverifiable from the text. That is a real publication gap. The Python snippet is also incomplete as printed, though that is minor.\n\nThese are addressable. The derivation itself is self-contained, non-circular, and formally reproducible. No fitting, no invented parameters. The contribution is a convenience formula, not new physics, but it is a convenience that people using the Lambert W model will likely want.\n\nWho it is for: researchers exploring Lambert W cosmological models, or anyone wanting a closed-form check on numerical integrals of this type. It deserves a serious referee: the math is correct and the paper is honest, and the missing figure plus an explicit caveat about equations (7)-(8) can be fixed in revision.","headline":"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.","tokens_in":6836,"tokens_out":1669,"would_cite":true,"duration_ms":19599,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["Lambert W function","cosmic fluid","Hubble parameter","equation of state","analytic integration","redshift","cosmic chronometers","logarithmic integral"],"falsifier":"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).","tokens_in":5951,"feed_emoji":"🌌","tokens_out":9312,"duration_ms":95718,"temperature":0.7,"pith_summary":"The paper's claim is that ln(H(z1)/H0), the logarithm of the Hubble factor relative to its present-day value, has an exact closed-form expression in the Lambert W cosmic fluid model. The expression is the sum of the three definite integrals (31), (32), and (33), so H(z1)/H0 is obtained by exponentiating that sum. If the claim is right, numerical quadrature of the Hubble integral is unnecessary for exploring this model, and the model's parameter dependence becomes explicit. The paper supports the claim by comparing the analytic formula with an independently computed numerical integration of the same integral and with cosmic-chronometer data used in the recent fit.","feed_headline":"Cosmic Hubble ratio becomes an exact closed form","feed_subtitle":"For the Lambert W cosmic fluid, ln(H/H0) becomes the sum of three exact integrals, so numerical integration drops out.","key_machinery":"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).","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"It supplies the integral formula for ln(H/H0) and the fitted θ values that the analytic formula is meant to replace.","marker":"[4]"},{"why":"It introduces the Lambert W equation of state whose Hubble factor the paper integrates.","marker":"[1]"},{"why":"It supplies the logarithmic-integral techniques used for the change of variable y=W(x).","marker":"[8]"},{"why":"It provides the 32 cosmic chronometer measurements used in the graphical comparison.","marker":"[5]"}],"fun_headline_variants":["Exact form for cosmic Hubble ratio in Lambert W fluid","Lambert W fluid: H(z)/H0 becomes closed-form","Analytic integral replaces numerics for Lambert W cosmology","Cosmic expansion formula: closed-form Lambert W solution","Exact H(z)/H0 via Lambert W antiderivative"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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)=Ω.","fun_headline_variants_meta":{"raw":{"variants":["Exact form for cosmic Hubble ratio in Lambert W fluid","Lambert W fluid: H(z)/H0 becomes closed-form","Analytic integral replaces numerics for Lambert W cosmology","Cosmic expansion formula: closed-form Lambert W solution","Exact H(z)/H0 via Lambert W antiderivative"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000538,"raw_usage":{"total_tokens":2536,"prompt_tokens":856,"completion_tokens":1680,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":472,"completion_tokens_details":{"reasoning_tokens":1597}},"tokens_in":472,"tokens_out":1680,"duration_ms":11438,"temperature":1.0,"reasoning_tokens":1597,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:43:17.667689+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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).","supporting_citations":[{"cited_title":"The LambertWequation of state in light of DESI BAO","cited_arxiv_id":null,"evidence_quote":"It supplies the integral formula for ln(H/H0) and the fitted θ values that the analytic formula is meant to replace."},{"cited_title":"The Lambert $W$ function: A newcomer in the Cosmology class?","cited_arxiv_id":"1908.11712","evidence_quote":"It introduces the Lambert W equation of state whose Hubble factor the paper integrates."},{"cited_title":"[end] ASBH-20260807 Page 13 07-Aug-2026","cited_arxiv_id":null,"evidence_quote":"It supplies the logarithmic-integral techniques used for the change of variable y=W(x)."}],"review_version":1}