{"id":"a43d1110-fa59-483f-8844-a1b5b8977619","arxiv_id":"1909.02360","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A closed-form hypergeometric and Taylor treatment of dark-energy and Cardassian distance integrals is fitted to SN and GRB data, but the key integral identity and the data treatment are flawed.","lead":"This paper rewrites the Hubble-radius integral for wCDM, evolving dark energy, and Cardassian cosmologies as hypergeometric functions and Taylor series, then fits them to supernova and gamma-ray burst distance moduli. The printed central substitution appears algebraically wrong and the fits ignore published covariance matrices, so the reported cosmological parameters are not reliable.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (12) has an incorrect t exponent in the change of variable, and Eq. (13) fails a direct derivative test at ΩM=1; the analytic distance modulus and all Section 6 fits are built on this formula.","rationale":"The reader's weakest assumption identified the change of variable in Eq. (12) and the validity of Eq. (13) as the hinge of the paper. Independent derivation confirms the substitution error: the printed Eq. (12) is off by a factor t^{1/3}. Moreover, Eq. (13) itself cannot be the antiderivative, even if Eq. (12) were corrected, because its ΩM=1 limit has the wrong z-scaling: a constant times (1+z)^{-3/2} differentiates to a constant times (1+z)^{-5/2}, not (1+z)^{-3/2}. This is an internal algebraic inconsistency, not a matter of outside consensus or subtle statistics. Since the analytic distance modulus and the fits in Tables 1-3 all depend on Eq. (13)/(14), the central claim of the paper is unsupported. The statistical criticisms (ignored covariance matrices) and the novelty question are real but secondary; they do not need to be reached once the core formula fails. The reader's REJECT verdict is therefore correct and remains unchanged.","tokens_in":10823,"tokens_out":10584,"duration_ms":94634,"concrete_test":"Set ΩM=1 in Eq. (13), where the hypergeometric argument vanishes, and substitute t=(1+z)^3. The right-hand side becomes a constant times (1+z)^{-3/2}; differentiating with respect to z gives a constant times (1+z)^{-5/2}, while Eq. (10) with ΩM=1 is identically (1+z)^{-3/2}. Since (d/dz)(1+z)^{-3/2} ∝ (1+z)^{-5/2} cannot equal (1+z)^{-3/2} for all z, Eq. (13) is provably not the indefinite integral of Eq. (10). A numerical check at z=1 confirms the ratio is (1/2)(1+z)^{-1}=1/4, not a constant. This single derivative test settles that the central analytic construction is incorrect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (13)/(14) supplies the exact indefinite integral of the wCDM Hubble-distance integrand Eq. (10), with Eqs. (15), (34) and (36) then giving exact luminosity distances and distance moduli. This fails already at the change of variable. With t=(1+z)^3, one has dz=(1/3)t^{-2/3}dt and (1+z)^{3+3w}=t^{1+w}, so the integrand in t is (1/3)t^{-7/6}[ΩM+(1-ΩM)t^w]^{-1/2}. The printed Eq. (12) instead contains t^{2/3} inside the square root, giving (1/3)t^{-5/6}[ΩM+(1-ΩM)t^w]^{-1/2}, a factor t^{1/3}=(1+z) too large. More seriously, Eq. (13) is not an antiderivative of either Eq. (12) or Eq. (10). Setting ΩM=1 makes the hypergeometric argument 0, so Eq. (13) reduces to a constant times (1+z)^{-3/2}; differentiating gives a constant times (1+z)^{-5/2}, whereas Eq. (10) at ΩM=1 is (1+z)^{-3/2}. The mismatch is z-dependent, so no choice of normalization or of the regularized-hypergeometric Gamma factor can repair it. All later results—Eqs. (14), (15), (34)–(36), the explicit distance-modulus formulas, and the Section 6 fits—inherit this error. The advertised analytic hypergeometric solution is therefore unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims to derive exact analytical expressions, in terms of the regularized Gauss hypergeometric function, for the comoving distance integral in flat wCDM cosmology (Eq. (13)) and in flat Cardassian cosmology (Eq. (31)). On this basis it proposes analytic luminosity distances and distance moduli (Eqs. (34)-(36) and (42)-(43)), derives seventh-order Taylor approximations for the constant and variable equation-of-state cases, and fits the Union 2.1, JLA, and Union 2.1 + Hymnium GRB samples to determine H0, Omega_M, and w (or w0, w1, or n) using a Levenberg-Marquardt chi-square minimization.","tokens_in":11124,"tokens_out":3736,"duration_ms":33637,"significance":"If the central hypergeometric identities were correct, the paper would provide a compact and useful closed form for the Hubble-distance integral and would make the full set of distance-modulus fits straightforward. However, the derivation fails at the first step: the change of variable in Eq. (12) is algebraically wrong, and Eq. (13) is not an antiderivative of the original integrand even in the Omega_M=1 limit. Since the exact distance-modulus formulas in Section 6 are built on this unsupported expression, the advertised analytical solution is not established. The paper also contains useful, carefully listed Taylor coefficients (Appendices A and C), but they do not rescue the main claim. Given the fundamental nature of the error, the significance of the paper in its present form is minimal.","major_comments":[{"comment":"The change of variable 1+z = t^{1/3} gives dz = (1/3) t^{-2/3} dt and (1+z)^{3+3w} = t^{1+w}, so the integrand in Eq. (11) becomes (1/3) t^{-7/6} [Omega_M + (1-Omega_M) t^w]^{-1/2}. The printed Eq. (12), with t^{2/3} inside the square root, produces instead (1/3) t^{-5/6} [Omega_M + (1-Omega_M) t^w]^{-1/2}, which is too large by a factor t^{1/3}. This error is inherited by every later expression that derives from Eq. (12).","section":"§3.1, Eq. (12)"},{"comment":"Eq. (13) is not an antiderivative of Eq. (10). Setting Omega_M=1 makes the hypergeometric argument vanish, reducing Eq. (13) to a constant times t^{-1/2}; differentiating with respect to z yields a term proportional to (1+z)^{-5/2}, whereas Eq. (10) at Omega_M=1 is (1+z)^{-3/2}. The mismatch is z-dependent, so it cannot be removed by any choice of normalization or by the Gamma-function prefactor of the regularized hypergeometric function. Consequently Eqs. (14), (15), (34), and (36), and all fits based on them, are unsupported.","section":"§3.1, Eq. (13)"},{"comment":"The Cardassian change of variable has the same defect. With t=(1+z)^3, the correct integrand is (1/3) t^{-7/6} [Omega_M + (1-Omega_M) t^{n-1}]^{-1/2}, but Eq. (30) contains an extra t^{2/3} factor inside the square root and an incorrect term structure. Eq. (31) therefore does not follow from Eq. (28), and the Cardassian distance modulus in Eq. (43) lacks a valid derivation.","section":"§5, Eq. (30)"},{"comment":"The chi-square in Eq. (7) uses only the diagonal uncertainties and ignores the full covariance matrix of the supernova samples, particularly for JLA, so the quoted parameter errors and the reported Q values are statistically incomplete. In addition, Table 2 shows an internal inconsistency: the wCDM Taylor approximation gives Omega_M=0.133 and w=-0.709, far from the hypergeometric result Omega_M=0.293 and w=-0.996, even though both fits are meant to describe the same wCDM model. This discrepancy indicates a breakdown in the Taylor approximation that is not discussed.","section":"§6, Tables 1-3 and Eq. (7)"}],"minor_comments":[{"comment":"In the sentence 'In the case of Ω_K = we have the flat case' the value of Ω_K is missing; presumably Ω_K=0.","section":"§2.1"},{"comment":"Eq. (44) is malformed: the factors '×' and '+ 3.4146' appear without clear placement, and Eq. (45) contains the confusing expression '1/ln(10) 25 ln(10)'. These formulas should be rederived and re-typeset.","section":"§6, Eq. (44)"},{"comment":"The cross-reference 'equation (reficardz)' in the conclusion is broken and should refer to Eq. (32).","section":"§7"},{"comment":"The paper states that Eq. (13) uses the regularized hypergeometric function, but the text and Appendix B denote it by the same symbol 2F1 used for the ordinary hypergeometric function; this notational ambiguity is confusing because the Gamma prefactor changes the Omega_M=1 limit.","section":"§3.1, Eq. (13) and Appendix B"}],"recommendation":"reject","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper claims new analytic hypergeometric expressions for the Hubble-distance integral in wCDM and Cardassian cosmologies, plus Taylor approximations and fits to SN/GRB samples. The central formula is wrong. Set ΩM=1 in Eq. (13): the hypergeometric argument vanishes, leaving −(1/3)(1+z)^{-3/2}. Its derivative with respect to z is (1/2)(1+z)^{-5/2}, but the integrand in Eq. (10) at ΩM=1 is (1+z)^{-3/2}. So Eq. (13) is not an antiderivative of Eq. (10) (or of Eq. (12), which is actually the correct substitution as far as I can read the typesetting). This is not a normalization slip; the z-dependence is wrong. Everything built on Eq. (13)—Eqs. (14), (15), the distance-modulus formulas, and the “Hypergeometric solution” fits in Tables 1–3—inherits the error.\n\nThe Taylor expansions themselves appear routine and the coefficients in the appendices are carefully generated. The paper is also honest that the hypergeometric dependence for wCDM was already recognized in Ref. [15]; the Cardassian version is a direct analogue, so novelty is low. The statistical treatment has separate problems: the χ² ignores off-diagonal covariance in JLA, and Table 2 shows a Taylor fit (ΩM=0.133, w=−0.709) that is inconsistent with the hypergeometric fit at nearly the same χ², which suggests one of the two is badly broken.\n\nWho is this for? Someone wanting a fast closed-form distance modulus for wCDM, but since the closed form is wrong, the paper would only mislead. The flaw is elementary and demonstrable, so I would not send this to a full referee cycle; a desk reject is appropriate. If the author can fix the antiderivative—and there is a correct hypergeometric form in the literature—then there might be a niche result for Cardassian, but the current version is not salvageable as is.","headline":"The paper's advertised hypergeometric solution for the wCDM Hubble-distance integral is not an antiderivative; a simple Ω_M=1 check kills the central claim, so the paper should not be published.","tokens_in":11741,"tokens_out":6718,"would_cite":false,"duration_ms":60408,"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":"The paper derives closed-form expressions for the Hubble-distance integral in wCDM and Cardassian cosmologies using the regularized Gaussian hypergeometric function, enabling analytic distance moduli and parameter fits to supernova and…","keywords":["dark energy","wCDM cosmology","Cardassian cosmology","hypergeometric function","distance modulus","Type Ia supernovae","equation of state","gamma-ray bursts"],"falsifier":"Differentiate the closed-form expression in Eq. (14) numerically with respect to $z$ at several points (for example $z=0.2,0.5,1.0$ with $\\Omega_M=0.3$, $w=-1$) and subtract the integrand $d_H(z;\\Omega_M,w)$ from Eq. (10); if the residual exceeds $10^{-6}$, the hypergeometric antiderivative is not a true antiderivative.","tokens_in":10506,"feed_emoji":"🔭","tokens_out":11430,"duration_ms":100620,"temperature":0.7,"pith_summary":"The paper claims that the integral of the Hubble radius in wCDM cosmology—and in flat Cardassian cosmology—can be expressed exactly in terms of a regularized Gaussian hypergeometric function, replacing numerical integration for the luminosity distance and distance modulus. This yields analytic distance moduli that the paper fits to Type Ia supernova samples (Union 2.1, JLA, and the Hymnium gamma-ray-burst sample) to obtain cosmological parameters such as $H_0 \\approx 70\\,\\mathrm{km\\,s^{-1}\\,Mpc^{-1}}$, $\\Omega_M \\approx 0.28$, and $w \\approx -1$. The paper also derives seventh-order Taylor polynomial approximations about $z=0$ and $z=1$, giving a fast approximate route for the variable-equation-of-state case. The central step is a substitution $t=(1+z)^3$ that turns the integrand into a binomial radical whose antiderivative is hypergeometric.","feed_headline":"Hypergeometric formula solves dark-energy distance modulus","feed_subtitle":"Closed-form integral lets supernova and GRB samples yield H0, Omega_M, and w without numerics.","key_machinery":"The central object is the regularized Gaussian hypergeometric function $_2\\tilde{F}_1(a,b;c;z)$, the standard $_2F_1$ normalized by Gamma functions, which appears as the exact antiderivative of the reciprocal of the Hubble function after the substitution $t=(1+z)^3$. This identity is the load-bearing result of the paper: it turns the Hubble-distance integral into an expression that can be evaluated pointwise, from which the luminosity distance and distance modulus follow directly, and it also yields the two Taylor approximations (order 7 about $z=0$ and order 2 about $z=1$) by series-expanding the integrand before integrating. The same substitution produces the Cardassian antiderivative in Eq. (31).","core_discovery":"The central discovery is an explicit antiderivative for the Hubble-distance integrand in a flat wCDM universe, written with the regularized Gaussian hypergeometric function $_2\\tilde{F}_1(a,b;c;z)$ in Eq. (13), and an analogous expression for the flat Cardassian model in Eq. (31). From these, the luminosity distance $d_L = \\frac{c}{H_0}(1+z)F(z;\\Omega_M,w)$ and the distance modulus $(m-M)=25+5\\log_{10}(d_L)$ follow analytically, letting the paper obtain best-fit parameters from Type Ia supernova data via the Levenberg-Marquardt method, for example $H_0\\approx70\\,\\mathrm{km\\,s^{-1}\\,Mpc^{-1}}$, $\\Omega_M\\approx0.28$, and $w\\approx-1$ for wCDM, and $n\\approx-0.08$ for Cardassian cosmology.","pith_inferences":["Because the distance modulus is now an analytic function of the parameters, its gradients with respect to $H_0$, $\\Omega_M$, and $w$ could be computed in closed form, which would accelerate and stabilize chi-square fits relative to numerical differentiation.","Special values of $w$ may reduce the hypergeometric form to elementary or elliptic functions, giving compact distance-modulus formulas for regimes such as $w=-1/3$ or $w=0$.","The same $t=(1+z)^3$ substitution should apply to other models in which the Hubble function is a binomial radical in $t$, such as extended dark-energy or curvature models, although the paper does not treat those cases."],"forward_implications":["The distance modulus in wCDM and Cardassian cosmologies can be evaluated in closed form, avoiding numerical quadrature of the Hubble function.","From the Union 2.1, JLA, and Hymnium samples, the paper derives cosmological parameters such as $H_0\\approx70\\,\\mathrm{km\\,s^{-1}\\,Mpc^{-1}}$, $\\Omega_M\\approx0.28$, $w\\approx-1$, with $\\chi^2_{\\mathrm{red}}\\approx0.85\\text{--}0.98$ depending on the sample.","The seventh-order Taylor approximation about $z=0$ and the second-order approximation about $z=1$ reproduce the exact integral to within a small percent error, with the crossover at $z\\approx0.58$.","For the variable equation of state $w(z)=w_0+w_1 z/(1+z)$, the Taylor approximation extends the analytic treatment to that case."],"supporting_citations":[{"why":"Defines the wCDM Hubble radius (Eq. 3.4) that the paper integrates.","marker":"[13]"},{"why":"Supplies the luminosity-distance formulas for wCDM and Cardassian models that the paper extends to closed form.","marker":"[14]"},{"why":"Introduces the flat Cardassian expansion model whose Hubble radius is integrated in Section 5.","marker":"[4]"},{"why":"Recognized the hypergeometric dependence of the integral but did not develop it, motivating the new solution.","marker":"[15]"},{"why":"Standard reference for the Gaussian hypergeometric function and the regularized form used in Eq. (13).","marker":"[21]"},{"why":"Standard reference for the hypergeometric function's analytical properties and series definition.","marker":"[25]"},{"why":"Provides the Union 2.1 compilation of 580 supernovae used in the fits of Tables 1 and 3.","marker":"[9]"},{"why":"Provides the JLA compilation of 740 supernovae used in the fits of Table 2.","marker":"[10]"},{"why":"Provides the calibrated Hymnium gamma-ray-burst sample that extends the distance-modulus fits to high redshift.","marker":"[11]"},{"why":"Supplies the Levenberg-Marquardt subroutine (MRQMIN) used to minimize chi-square in the parameter fits.","marker":"[16]"}],"fun_headline_variants":["Hypergeometric formula delivers exact distance modulus","Exact dark-energy distance modulus via hypergeometric","One hypergeometric solution for distance modulus in modern cosmology","Closed-form distance modulus from hypergeometric functions","Hypergeometric distance modulus for dark energy and Cardassian"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the change of variable in Eq. (12) is algebraically correct, since the printed exponent appears to insert an extra factor $t^{1/3}$, which would make the hypergeometric antiderivative fail to differentiate back to the Hubble function.","fun_headline_variants_meta":{"raw":{"variants":["Hypergeometric formula delivers exact distance modulus","Exact dark-energy distance modulus via hypergeometric","One hypergeometric solution for distance modulus in modern cosmology","Closed-form distance modulus from hypergeometric functions","Hypergeometric distance modulus for dark energy and Cardassian"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001118,"raw_usage":{"total_tokens":4608,"prompt_tokens":857,"completion_tokens":3751,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":473,"completion_tokens_details":{"reasoning_tokens":3675}},"tokens_in":473,"tokens_out":3751,"duration_ms":27960,"temperature":1.0,"reasoning_tokens":3675,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:55:13.414882+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Differentiate the closed-form expression in Eq. (14) numerically with respect to $z$ at several points (for example $z=0.2,0.5,1.0$ with $\\Omega_M=0.3$, $w=-1$) and subtract the integrand $d_H(z;\\Omega_M,w)$ from Eq. (10); if the residual exceeds $10^{-6}$, the hypergeometric antiderivative is not a true antiderivative.","supporting_citations":[{"cited_title":"Utilizing the Updated Gamma-Ray Bursts and Type Ia Supernovae to Constrain the Cardassian Expansion Model and Dark Energy","cited_arxiv_id":"1504.02308","evidence_quote":"Supplies the luminosity-distance formulas for wCDM and Cardassian models that the paper extends to closed form."},{"cited_title":"Analytical Approach for the Determination of the Luminosity Distance in a Flat Universe with Dark Energy","cited_arxiv_id":"1003.0483","evidence_quote":"Recognized the hypergeometric dependence of the integral but did not develop it, motivating the new solution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Standard reference for the Gaussian hypergeometric function and the regularized form used in Eq. (13)."},{"cited_title":"(Cambridge: Cambridge University Press","cited_arxiv_id":null,"evidence_quote":"Standard reference for the hypergeometric function's analytical properties and series definition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Union 2.1 compilation of 580 supernovae used in the fits of Tables 1 and 3."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the JLA compilation of 740 supernovae used in the fits of Table 2."},{"cited_title":"Observational Constraints on Cosmological Models with the Updated Long Gamma-Ray Bursts","cited_arxiv_id":"1004.4951","evidence_quote":"Provides the calibrated Hymnium gamma-ray-burst sample that extends the distance-modulus fits to high redshift."},{"cited_title":"The Art of Scientific Computing (Cambridge, UK: Cambridge University Press)","cited_arxiv_id":null,"evidence_quote":"Supplies the Levenberg-Marquardt subroutine (MRQMIN) used to minimize chi-square in the parameter fits."}],"review_version":1}