{"id":"3f5fe423-ec53-4c3a-8274-9a2662459cdd","arxiv_id":"1909.01998","paper_version":4,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"A Brans-Dicke-Lyra model is fitted to H(z) and supernova data, but the acceleration claim conflicts with the model's own field equations and the fit reduces to ΛCDM.","lead":"This paper fits a Brans-Dicke scalar field model in Lyra's geometry to Hubble and supernova data and reports a universe that switches from decelerating to accelerating expansion. The fitted parameters match Planck and WMAP, but the model's own equations conflict on the sign of the acceleration, and the model is observationally indistinguishable from standard dark energy cosmology.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (20) does not satisfy the scalar wave equation (7): with the paper's own parameters, Eq. (7) gives q0 ≈ +1.58 while Eq. (23) gives q0 = -0.58, so the fitted H(z) is not a solution of the model.","rationale":"The reader's REJECT verdict is confirmed by an independent check. The most load-bearing premise is not merely that the power-law φ is a solution of Eq. (15) (it does satisfy that particular combination identically), but that the resulting H(z) in Eq. (20) solves the full system (5)-(7). I derived the deceleration parameter implied by the scalar wave equation (7) under the same ansatz and fitted parameters. At z = 0, Eq. (7) requires q0 ≈ +1.58, while the paper's kinematic formula Eq. (23) gives q0 = -0.58. This two-unit discrepancy means Eq. (20) is not a solution of the stated field equations. Since the observational constraints, the transition redshift, and the inferred age are all computed from Eq. (20), the central claim that the model explains cosmic acceleration without a cosmological constant collapses. I agree with the reader that this is the weakest assumption, though I would phrase it as a failure to satisfy the full system rather than a failure of the ansatz to be 'general.' The verdict remains REJECT; no adjustment is needed.","tokens_in":10911,"tokens_out":14349,"duration_ms":121387,"concrete_test":"Perform an independent symbolic substitution of φ = φ0(a/a0)^{1/(ω+1)} and H(z) from Eq. (20) into the scalar wave equation (7) at z = 0 with Ωm0 = 0.28 and ω = 40000. The residual of Eq. (7) should be zero for a true solution. Directly compare the deceleration parameter implied by Eq. (7), q0 = n + 2 − 3Ωm0(ω+1)/(2ω+3) ≈ +1.58, with the paper's Eq. (23) value q0 = −0.58. A non-zero residual or a disagreement in q0 by O(1) demonstrates that Eq. (20) is not a solution of the stated Brans-Dicke-Lyra field equations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper constructs H(z) by inserting the ansatz φ = φ0(a/a0)^{1/(ω+1)} into the Friedmann equation (5) together with ρ_m ∝ a^{-3}, but it never verifies that this H(z) satisfies the remaining field equations, and it cannot. For dust, the scalar wave equation (7) is φ¨/φ + 3Hφ˙/φ = 8πρ/[(2ω+3)φc^2]. Substituting φ˙/φ = nH and φ¨/φ = nH˙ + n^2H^2 with n = 1/(ω+1), then using (5) to eliminate ρ, gives at z = 0: H˙/H^2 = (1/n)[3Ωm0/(2ω+3) − n(n+3)], after eliminating Ωβ via Ωβ = 1 + C − Ωm. With ω = 40000 and Ωm0 = 0.28, this yields q0 = −1 − H˙/H^2 ≈ +1.58. Equation (23), obtained kinematically from the same H(z), gives q0 = −0.58. The two deceleration parameters differ by more than two, so Eq. (20) cannot simultaneously satisfy Eqs. (5), (6), and (7). The claim that the power-law φ is the 'general solution' of Eq. (15) is also misleading: Eq. (15) is identically satisfied for arbitrary H(t) once φ ∝ a^n, so it does not select an expansion history; the subsequent fit uses only Eq. (5) and thereby drops the scalar-wave constraint.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a flat FLRW Brans-Dicke cosmological model formulated in Lyra's geometry with a time-like displacement vector. After writing the field equations (5)-(7), the authors adopt a power-law scalar field φ ∝ a^{1/(ω+1)}, derive the Hubble function (20), and fit its two free parameters H0 and Ωm0 to 46 H(z) points and an SN Ia Gold sample. From the best fit they compute a present deceleration q0=-0.58, a transition redshift zt=0.52, an age of 14.29 Gyr, and j0=0.80, concluding that the model yields late-time cosmic acceleration without a cosmological constant.","tokens_in":11293,"tokens_out":8835,"duration_ms":84992,"significance":"If the model were correct, it would provide an interesting example of an accelerating dust-filled universe without dark energy, with the Brans-Dicke scalar field playing the role usually assigned to a cosmological constant in Lyra geometry. The paper reports a standard χ² likelihood analysis with contour plots, which is a useful route to falsifiability. However, the central expansion history is not a solution of the field equations: the best-fit H(z) violates the scalar-wave equation. Because every subsequent observable (q0, age, jerk, transition redshift) is derived from this H(z), the claimed cosmological conclusions do not follow. The paper also correctly notes that the Lyra displacement vector β cannot act as Λ in Brans-Dicke theory, but this insight is not developed into a workable model.","major_comments":[{"comment":"The best-fit Hubble function is not a solution of the model's field equations. For the ansatz (16), one has φ˙/φ = nH with n = 1/(ω+1). With dust (p=0), the scalar-wave equation (7) becomes n(\\dot H/H²) + n(n+3) = 3Ωm/(2ω+3). Evaluating this at z=0 with ω=40000 and Ωm0=0.28 gives q0 = -1 - \\dot H/H² ≈ +1.58, whereas Eq. (23), obtained from the fitted H(z), gives q0 = -0.58. Both expressions describe the same model at the same point, so Eq. (20) cannot satisfy Eqs. (5)-(7) simultaneously; the fitted H(z) solves only the Friedmann equation used to construct it.","section":"§III, Eqs. (20) and (23); §II, Eq. (7)"},{"comment":"The statement that Eq. (16) is 'the general solution' of Eq. (15) is incorrect. Substituting φ ∝ a^{1/(ω+1)} makes Eq. (15) an identity for arbitrary a(t); it imposes no restriction on H(z). Consequently, the power-law ansatz does not select an expansion history, and the derivation of Eq. (20) uses only Eq. (5) together with the density-parameter decomposition. This is the structural reason why the scalar-wave constraint is lost and why the fitted H(z) is not a solution of the full system.","section":"§III, Eqs. (15) and (16)"},{"comment":"The observational validation is partly circular and the data handling is inconsistent. The deceleration parameter, transition redshift, age, and jerk are all deterministic functions of H0 and Ωm0 obtained from the same OHD/SNIa fits; their agreement with WMAP/Planck values is therefore a restatement of the fit rather than an independent prediction. In addition, Table I is described in the text as cosmic-chronometer data, but it includes BAO measurements (for example, entries cited to refs. [49], [51], [53], [55], [57]-[59]). Cosmic-chronometer and BAO data have different systematics and should not be combined in a single 'OHD' likelihood without modeling those systematics.","section":"§III, Table I and Eq. (21); §IV, Eqs. (23)-(27)"}],"minor_comments":[{"comment":"There are numerous typographical errors ('obtian', 'Brans-Dike', 'Plank collaboration', 'co-moving co-ordinate', 'generak', 'red-shift' as one word) that should be corrected in a thorough proofread.","section":"Throughout"},{"comment":"The equation cross-references are unreliable: the sentence before Eq. (16) refers to 'equation (18)' when it should refer to Eq. (15), and the sentence before Eq. (18) says 'Using equation (23) in equation (6)' when it should refer to Eq. (16).","section":"§III, Eq. (16) and Eq. (18)"},{"comment":"The text says the jerk parameter is 'graphed in Figure 5', but Figure 5 is the age plot H0(t0−t) versus z; the jerk plot is Figure 6.","section":"§IV.C"},{"comment":"The SN Ia sample is described only as 'Gold Sample and New Gold Sample' from ref. [60]; this is an old compilation, and the paper should state why more recent samples (e.g., Pantheon or DES) were not used.","section":"§III, SN Ia data"}],"recommendation":"reject","confidential_remarks":"The manuscript cites ref. [29], a Comment on a closely related Lyra-Brans-Dicke model, but does not discuss whether the objections raised there apply to the present derivation. Given the inconsistency between Eq. (20) and Eq. (7) demonstrated above, the editor may wish to ask the authors to reconcile these two papers. The observational claim also rests on a mixed OHD/BAO sample presented as cosmic chronometers only."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Jim,\n\nBottom line: this one does not survive contact with its own equations. The fitted H(z) is not a solution of the Brans-Dicke field equations in Lyra geometry. Substituting their power-law φ into the scalar wave equation (7) yields q0 ≈ +1.58 with their OHD best-fit parameters, while Eq. (23) from the same H(z) gives q0 = −0.58. That is not a small numerical quibble; it means the central object of the paper, the expansion history used for the χ² fit, does not satisfy the theory it claims to test. I checked the algebra and the stress-test note is right.\n\nThe one thing the paper does well is present a clean, reproducible data exercise. The 46 OHD points are tabulated, the χ² procedure is transparent, and they give contours and best-fit values. If all you want is an example of a routine modified-gravity fit to H(z) and SNIa, that part is fine. But the physics added by the model is not. The “general solution” for φ is really just an ansatz; Eq. (15) is identically satisfied for any H(t) once φ ∝ a^n, so it does not select the expansion history. The fit uses only the Friedmann equation (5) and silently drops the wave equation.\n\nThere are other, softer problems. With ω = 40000 the model reduces numerically to flat ΛCDM, so the fit is not evidence for Lyra geometry or Brans-Dicke physics. The parameters come without uncertainties, and the OHD sample mixes cosmic chronometer and BAO points despite the text calling all 46 uncorrelated CC measurements. The age and jerk quoted later are derived from the same best-fit H(z), so they are not independent predictions.\n\nGiven the load-bearing contradiction, I would not send this to peer review; it needs a desk reject. The paper is honest in tone and cites its sources, but the inconsistent derivation is fatal on its own terms.\n\nNot worth citing, and I wouldn’t bring it to reading group.\n\nBest,\n[Name]","headline":"The fitted H(z) is not a solution of the model's own field equations, so the central claim fails; the paper is mostly a routine data fit of a known form with no new physics.","tokens_in":11881,"tokens_out":4431,"would_cite":false,"duration_ms":39002,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.kd","98.80.-k","98.80.JK"],"model":"deepseek-v4-flash","headline":"A Brans-Dicke scalar field in Lyra's geometry drives cosmic acceleration without a cosmological constant.","keywords":["Brans-Dicke theory","Lyra geometry","scalar field cosmology","late-time acceleration","deceleration parameter","observational Hubble data","Type Ia supernovae","cosmological parameters"],"falsifier":"Substitute Eq. (20) and the power-law $\\varphi$ back into the Raychaudhuri equation (Eq. 6) and compute the deceleration parameter $q$ directly from the field equations; if it does not reproduce Eq. (23) for the fitted $\\Omega_{m0}\\approx 0.28$ and $\\omega = 40000$, the fitted expansion history is not a solution of the model as written.","tokens_in":10650,"feed_emoji":"🌌","tokens_out":15868,"duration_ms":136839,"temperature":0.7,"pith_summary":"This paper claims that a Brans-Dicke scalar field living in Lyra's geometry can reproduce the observed expansion history of a flat, homogeneous, isotropic universe without a cosmological constant or a separate dark-energy fluid. The authors solve the field equations for pressureless matter using a power-law scalar field, obtain a closed-form Hubble function, and fit it to 46 cosmic-chronometer $H(z)$ points plus a Type Ia supernova sample, finding $H_0 \\approx 67.44$ km/s/Mpc and $\\Omega_{m0}\\approx 0.28$ from the Hubble data. The same solution gives a deceleration-to-acceleration transition at $z_t\\approx 0.52$, a present deceleration $q_0\\approx -0.58$, an age of about 14.29 Gyr, and a jerk parameter near 0.80. If correct, the model would show that the observed late-time acceleration can be a dynamical scalar-field effect rather than the cosmological constant.","feed_headline":"A scalar field alone drives cosmic acceleration, no dark energy","feed_subtitle":"Fits Hubble and supernova data, flips deceleration at z≈0.52, and sets the universe's age at 14.29 Gyr.","key_machinery":"The carrying object is the power-law scalar field ansatz $\\varphi = \\varphi_0 (a/a_0)^{1/(\\omega+1)}$, taken together with the Lyra displacement vector $\\beta$ inside the effective density and pressure. This ansatz converts the two Friedmann-type equations into a single algebraic constraint, $\\Omega_m + \\Omega_\\beta = 1 + (5\\omega+6)/(6(\\omega+1)^2)$, which fixes the $\\beta$-energy fraction whenever the matter density and $\\omega$ are specified. Feeding that constraint into the Friedmann equation produces the closed-form Hubble function in Eq. (20); the rest of the paper is the cosmographic series derived from that function: $q = -1 + (1+z)H'(z)/H(z)$ for deceleration, the age integral $t_0 = \\int_0^\\infty dz / [(1+z)H(z)]$, and the jerk expression obtained from derivatives of $H(z)^2$.","core_discovery":"The paper's central claim is that the expansion of a pressureless universe obeying the Brans-Dicke field equations in Lyra's manifold is governed by $H(z) = H_0\\,[1+C]^{-1/2}\\,\\bigl[1+\\Omega_{m0}(1+z)^{(3\\omega+4)/(\\omega+1)} + C - \\Omega_{m0}\\bigr]^{1/2}$, where $C = (5\\omega+6)/(6(\\omega+1)^2)$, and that this history fits both the 46-point $H(z)$ sample and the Type Ia supernova sample used in the paper. Fitting yields $H_0=67.44$ km/s/Mpc with $\\Omega_{m0}=0.28$ from the Hubble data and $H_0=70.02$ km/s/Mpc with $\\Omega_{m0}=0.272$ from the supernova data. The same $H(z)$ is then differentiated to obtain a deceleration parameter that changes sign at $z_t\\approx 0.52$ and equals $q_0\\approx -0.58$ today, integrated to give an age of about 14.29 Gyr, and differentiated twice to give a present jerk $j_0\\approx 0.80$. The authors read these numbers as showing that the model describes a universe that decelerates early, accelerates now, and needs no cosmological constant, with the acceleration supplied by the dynamics of the Brans-Dicke scalar field rather than by Lyra's displacement vector.","pith_inferences":["Not stated in the paper: at $\\omega \\approx 40000$ the correction $C = (5\\omega+6)/(6(\\omega+1)^2)$ is tiny, so the fitted $H(z)$ closely tracks flat $\\Lambda$CDM with a rescaled matter density; the successful fit may be largely a test of how close the ansatz is to $\\Lambda$CDM, not of the Lyra-specific terms.","Not stated in the paper: since the displacement vector has effective equation of state $+1$ in the matter-free limit, the acceleration must be carried by the scalar field; a high-precision measurement of the late-time effective equation of state would separate these two contributions.","Not stated in the paper: applying the same power-law scalar-field ansatz to anisotropic Bianchi models would give the Lyra displacement field a testable shear or anisotropy signature that the isotropic fit cannot reveal."],"forward_implications":["If Eq. (20) is the correct expansion history, the universe's transition from deceleration to acceleration occurs at $z_t\\approx 0.52$, so independent high-redshift probes of the expansion rate should see deceleration above that redshift.","The model yields a present age of about 14.29 Gyr, which is consistent with independent CMB-based age estimates and removes the need for a cosmological constant to stretch the age.","The fitted matter density $\\Omega_{m0}\\approx 0.28$ and Hubble constant $H_0\\approx 67.44$ km/s/Mpc from the Hubble data fall inside the ranges reported by CMB experiments, so the model is not immediately excluded by these basic expansion probes.","Because the predicted present jerk $j_0\\approx 0.80$ differs from the $\\Lambda$CDM value $j=1$, future cosmographic measurements of higher-order expansion can distinguish this scalar-field model from vacuum-energy cosmology."],"supporting_citations":[{"why":"Introduces Lyra's gauge geometry and the displacement vector field used in the modified field equations.","marker":"[15]"},{"why":"Supplies the scalar-tensor field equations that the model joins with Lyra's geometry.","marker":"[16]"},{"why":"Gives the combined Brans-Dicke field equations in Lyra's manifold that the paper takes as its starting point.","marker":"[31]"},{"why":"Provides the solar-system bound $\\omega \\approx 40000$ used in the numerical evaluations of $q_0$, age, and jerk.","marker":"[36, 37]"},{"why":"Contributes cosmic-chronometer $H(z)$ measurements to the 46-point sample used in the Hubble-rate fit.","marker":"[47]"},{"why":"Supplies additional cosmic-chronometer $H(z)$ measurements used in the same observational Hubble data fit.","marker":"[48]"},{"why":"Provides low-error $H(z)$ data used in the observational Hubble data fit.","marker":"[52]"},{"why":"Provides the Gold and New Gold Type Ia supernova sample used for the second parameter fit.","marker":"[60]"},{"why":"Gives the CMB-based reference values for the Hubble constant and age that the model's results are compared against.","marker":"[70]"},{"why":"Provides independent CMB-based cosmological parameter estimates used to check the model's $H_0$, matter density, and age.","marker":"[71]"}],"fun_headline_variants":["Cosmic acceleration without dark energy, just a scalar field","Brans-Dicke scalar field fits supernova and Hubble data, no Λ","Scalar field alone explains cosmic acceleration, data agrees","No dark energy needed: Brans-Dicke scalar field drives expansion","Universe accelerates on scalar field alone, matches observations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the assumed power-law form of the scalar field really solves the full field equations, so the fitted $H(z)$ is a true solution of the theory and not just a convenient curve; if that premise gives way, the derived expansion history and its parameters no longer describe the model.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic acceleration without dark energy, just a scalar field","Brans-Dicke scalar field fits supernova and Hubble data, no Λ","Scalar field alone explains cosmic acceleration, data agrees","No dark energy needed: Brans-Dicke scalar field drives expansion","Universe accelerates on scalar field alone, matches observations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000546,"raw_usage":{"total_tokens":2642,"prompt_tokens":1011,"completion_tokens":1631,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":1546}},"tokens_in":627,"tokens_out":1631,"duration_ms":11524,"temperature":1.0,"reasoning_tokens":1546,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:05:50.710009+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Substitute Eq. (20) and the power-law $\\varphi$ back into the Raychaudhuri equation (Eq. 6) and compute the deceleration parameter $q$ directly from the field equations; if it does not reproduce Eq. (23) for the fitted $\\Omega_{m0}\\approx 0.28$ and $\\omega = 40000$, the fitted expansion history is not a solution of the model as written.","supporting_citations":[{"cited_title":"Lyra, Ubereine Modifikation der Riemannschen Ge- ometrie, Math","cited_arxiv_id":null,"evidence_quote":"Introduces Lyra's gauge geometry and the displacement vector field used in the modified field equations."},{"cited_title":"Brans, R","cited_arxiv_id":null,"evidence_quote":"Supplies the scalar-tensor field equations that the model joins with Lyra's geometry."},{"cited_title":"galaxy differ- ential age","cited_arxiv_id":null,"evidence_quote":"Gives the combined Brans-Dicke field equations in Lyra's manifold that the paper takes as its starting point."},{"cited_title":"Zhang et al, Four new observational H(z) data from lu- minous red galaxies in the Sloan Digital Sky Survey data release seven, Res","cited_arxiv_id":null,"evidence_quote":"Contributes cosmic-chronometer $H(z)$ measurements to the 46-point sample used in the Hubble-rate fit."},{"cited_title":"Simon, L","cited_arxiv_id":null,"evidence_quote":"Supplies additional cosmic-chronometer $H(z)$ measurements used in the same observational Hubble data fit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides low-error $H(z)$ data used in the observational Hubble data fit."},{"cited_title":"Delubac et al, Which fundamental constants for cos- mic microwave background and baryon-acoustic oscilla- tion?, Astron & Astrophys 584 (2015) A69","cited_arxiv_id":null,"evidence_quote":"Provides the Gold and New Gold Type Ia supernova sample used for the second parameter fit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the CMB-based reference values for the Hubble constant and age that the model's results are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides independent CMB-based cosmological parameter estimates used to check the model's $H_0$, matter density, and age."}],"review_version":1}