{"id":"06892e36-39a5-4d88-bf19-e7524a7a2dd3","arxiv_id":"2501.12585","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An exponential f(Q) gravity model, a small perturbation of ΛCDM, gives quintessence for one sign of its parameter b and phantom behavior for the other, while changing late-time structure growth.","lead":"This paper studies a modified gravity model in which dark energy arises from an exponential correction to the standard gravity action, and computes how it changes the universe's expansion, energy conditions, and the growth of galaxies. It finds the model can mimic both quintessence and phantom dark energy depending on a sign parameter, and that it slightly alters structure growth at late times.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Section 7 structure-growth predictions depend on Eq. (92), whose validity is not established for this model because f_QQ≠0 triggers the known strong-coupling/ghost pathology; the paper's own Sec. 8 acknowledges but does not resolve this.","rationale":"The central claim has two parts: background dark-energy phenomenology (quintessence/phantom, energy conditions, statefinders) and structure-growth predictions. The background part is supported by the small-b analytical solutions validated against numerical integration to ≤0.008%, and by the transparent two-parameter model. The structure-growth part, however, is the unique part that goes beyond the authors' earlier work [24] and is the basis for the claim that the model 'impacts the growth of structures in the universe.' That part depends on a single borrowed equation, Eq. (92), together with the substitution G_eff=G/f_Q. The known strong-coupling/ghost pathology of f(Q) theories (which the paper itself cites at [61]) directly threatens this substitution: if the scalar degree of freedom is strongly coupled, the linear perturbation equation is not derived from a well-defined quadratic action, and the quasistatic approximation is uncontrolled. The paper's own conclusion ('ghost fields may be propagating...') shows the authors are aware of the issue, but awareness without a concrete check leaves the growth predictions ungrounded. This is not a disagreement with consensus about, say, parameter values; it is a question of whether the perturbation equations used are the correct ones for this model. A conditional verdict is therefore appropriate: accept the background analysis, but require either a demonstration of ghost-freedom for the exponential model (e.g., by fitting the conditions of [62]) or an explicit softening of the structure-growth claims. The H0 inconsistency (71.99 vs 67.36) and the visual fσ8 comparison are secondary; they affect precision, not the validity of the perturbation framework.","tokens_in":26391,"tokens_out":9567,"duration_ms":99432,"concrete_test":"Derive the quadratic scalar perturbation action for f(Q)=Q+2Λ exp[-(bΛ/Q)^n] around the FLRW background of Sec. 4, using the kinetic-matrix formalism of [61]; evaluate the kinetic coefficient at z=0 for b=-0.1, n=1 (and, if needed, b=0.1, n=2). If the coefficient is zero or negative, strong coupling or a ghost is present and Eq. (92) is invalid, so the Section 7 fσ8 curves should be discarded. If the coefficient is positive and the mode decouples in the quasistatic limit, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim—that the model impacts the growth of structures through a modified effective gravitational constant (Section 7, Figs. 19–23)—rests entirely on Eq. (92), the sub-horizon growth equation with G_eff = G/f_Q taken from [48]. This equation presumes a standard quasistatic scalar perturbation regime. For f(Q) with f_QQ ≠ 0, which is the case here for every b ≠ 0 (Eqs. 94–95), the scalar perturbation around FLRW has been shown to be strongly coupled or ghost-like (Gomes et al., PRL 132, 141401 (2024) = [61]). The exponential model does not obviously belong to the ghost-free class of [62], and the paper never checks whether the conditions of that class hold for the parameter range b ∈ [−0.2, 0.2], n = 1, 2. If the scalar mode is strongly coupled, the linearized equation (92) is not a valid effective description even in the sub-horizon limit, so the computed δm(z), fg(z), γ(z), and fσ8(z) are not trustworthy predictions. The paper itself flags this in Sec. 8 ('ghost fields may be propagating...') but only cites the possibility, never demonstrating that the model avoids it. Hence the perturbation-based part of the central claim is conditionally supported at best.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes the exponential f(Q) model f(Q) = Q + 2Λ exp[-(bΛ/Q)^n] in a flat FLRW universe. It derives approximate analytic Hubble-parameter solutions for n = 1 and n = 2 using a small-b expansion around ΛCDM, checks these against a numerical integration of the Friedmann ODE, and then computes the effective dark-energy equation of state, deceleration and statefinder parameters, Om(z), energy conditions, and linear matter perturbations. The perturbation analysis uses a sub-horizon equation with G_eff = G/f_Q to obtain δ_m, the growth factor f_g, the growth index γ, and fσ8, which is compared visually with 23 observational data points. The authors conclude that b > 0 gives quintessence-like behavior and b < 0 gives phantom-like behavior, and that the model produces small late-time deviations from ΛCDM that affect structure growth.","tokens_in":26624,"tokens_out":21441,"duration_ms":201446,"significance":"If the background and perturbation results were fully established, the paper would provide a convenient closed-form approximate treatment of a two-parameter f(Q) model, including analytic H(z) expressions with a documented numerical check and a concrete fσ8 comparison. The paper is transparent in stating that statistical constraints and the ghost-field issue are left for future work. However, the perturbation sector, which carries the main new claim about structure growth, rests on an equation whose validity for f_QQ ≠ 0 is not demonstrated, and the numerical validation of the background solutions contains an apparent algebraic error. These issues make the current support for the central claims conditional rather than conclusive.","major_comments":[{"comment":"The linear growth equation (92) is taken from Ref. [48] in the quasistatic sub-horizon limit with G_eff = G/f_Q. For every b ≠ 0 in this model, f_Q ≠ 1 and f_QQ ≠ 0 (Eqs. 94–95), and Ref. [61] shows that generic f(Q) modifications develop strong coupling or ghosts in the scalar sector. The manuscript never checks whether the exponential model belongs to the ghost-free class of Ref. [62] over the parameter range used here (b ∈ [−0.2, 0.2], n = 1, 2). Since Eqs. (100) and (101) and all of the δ_m, f_g, γ, and fσ8 results in Section 7 are built on Eq. (92), the structure-growth part of the central claim is not established. Section 8 explicitly acknowledges the possible ghost propagation but does not resolve it; a demonstration of the validity of the quasistatic approximation for this model is needed before the perturbation results can be accepted.","section":"Section 7, Eqs. (92) and (100)"},{"comment":"The numerical ODE displayed in Eq. (58) does not follow from Eq. (25) for n = 1. With g' = bΛ/(6H^2)^2, the coefficient of H_dot in Eq. (25) is 1 + Λ e^g [24H^2(g'' + g'^2) + 2g'], which reduces to 1 + e^{-H0^2 ΩΛ b/(2H^2)} [−(3/2) b H0^4 ΩΛ^2/H^4 + (1/2) b^2 H0^6 ΩΛ^3/H^6]. Equation (58) instead contains the combination b (2/3) H0^4 ΩΛ^2/H^4 − b H0^4 ΩΛ^2/(2H^4), which has the opposite sign and a different magnitude for the linear-in-b term. Because this equation is the stated basis for the reported 0.008% agreement between the analytic and numerical H(z), the numerical validation needs to be corrected and re-run before the analytic solutions are relied on.","section":"Section 4.3, Eq. (58)"}],"minor_comments":[{"comment":"The bracket in Eq. (33) appears to have an algebraic typo: 12H^2 n (bΛ)^n (6H^2)^{-n-1} reduces to 2n [bΛ/(6H^2)]^n, which becomes 2n [H0^2 ΩΛ b/(2H^2)]^n after substituting Λ = 3H0^2 ΩΛ; the printed expression with H0^2 factors and (2H^2)^{-n-1} is off by a factor of 3. The subsequent n = 1 and n = 2 equations are consistent with the corrected form, so this should be fixed for consistency.","section":"Equation (33)"},{"comment":"The first-order coefficient δu1 in Eq. (42) should be −3 ΩΛ^2/(2ξ), not −7 ΩΛ^2/(2ξ). The final solution (44) is consistent with the corrected coefficient, so this is a typographical error that should be corrected.","section":"Equation (42)"},{"comment":"The numerical check uses H0(−0.1) = 71.99 km/s/Mpc, while the rest of the paper uses H0 = 67.36 km/s/Mpc (e.g., Figs. 2–17). Please clarify the fiducial H0 convention and how H0(b) in Eq. (45) is related to the Planck value, since this affects the interpretation of the numerical comparison.","section":"Section 4.3"},{"comment":"The statement that the model 'fits nicely' the 23 fσ8 data points is based on a visual comparison; since no χ² or likelihood is quoted, please either quantify the agreement or explicitly label the plot as illustrative. The text already says that a statistical constraint is future work, so this is a presentation issue rather than a fatal one.","section":"Section 7, Fig. 23"},{"comment":"Some figure captions are not in English (e.g., 'Evolución de δm...' in Fig. 19) and the Fig. 1 caption contains 'Der' in parentheses; please unify the language and correct the caption typos.","section":"Figures and captions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a phenomenological extension of the same authors' earlier work in Ref. [24]. Its main new contributions beyond that paper are the n = 2 analytic solution, the energy conditions, the statefinder diagnostics, and the perturbation analysis. The perturbation section is the part that would justify publication as more than a background study, and it is currently unresolved because of the strong-coupling/ghost issue. If the authors can correct Eq. (58), re-run the numerical validation, and add either a rigorous justification of Eq. (92) for this model or an explicit caveat that the growth predictions are conditional on the ghost-free status of the model, the paper could be suitable for publication after major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful, workmanlike extension of the authors' own exponential f(Q) model. The genuinely new piece is the n=2 analytic H(z) solution, and they check both analytic solutions against the full numerical ODE with a maximum relative error of about 0.008%. That check is real evidence. The background and energy-condition phenomenology are handled cleanly and the ΛCDM limit is correctly recovered. The paper deserves a referee, but not a clean acceptance.\n\nStrengths: the small-b expansion is done honestly, the cubic and quartic roots are treated without hand-waving, and the statefinder and Om(z) diagnostics are standard but competently executed. The quintessence-like (b>0) versus phantom-like (b<0) classification is a direct consequence of the model's form, and the authors do not oversell it as a discovery. They also explicitly flag that a full statistical constraint on b is future work, which is the right level of modesty.\n\nSoft spots, in order of seriousness. First and most important: Section 7's growth results all rest on Eq. (92), the sub-horizon perturbation equation with G_eff = G/f_Q from [48]. The paper never addresses whether that equation is valid for this model. Since f_QQ ≠ 0 for every b ≠ 0, the scalar perturbations around FLRW in f(Q) can be strongly coupled or ghostly (Gomes et al., PRL 132, 141401, which is their [61]). The conclusions mention that ghosts 'may be propagating' and cite a ghost-free class [62], but they never show this exponential model belongs to that class. So the fσ8 curves and growth index are conditional at best. If Eq. (92) is not the correct effective description, those numbers are not trustworthy predictions. This does not sink the background analysis, but it does mean the paper's most observationally flavored claim is not established.\n\nMinor issues: the H0 inconsistency (71.99 km/s/Mpc in the numerical check, 67.36 elsewhere) needs a comment. The fσ8 comparison with 23 data points is visual only; the authors admit a statistical constraint is future work, but the text's 'aligns nicely' is overstatement without a chi-square. Also, the model is engineered as a small perturbation around ΛCDM, so small deviations are built in, not discovered.\n\nBottom line: conditional accept. The background and analytic solution work is publishable after minor revision. The perturbation section should either demonstrate that the model evades the strong-coupling/ghost issue or be explicitly labeled as illustrative pending that question. I would send this to a serious referee, not desk reject.","headline":"Solid background extension of an exponential f(Q) model, but the structure-growth predictions rest on a perturbation equation whose validity for this model is not demonstrated.","tokens_in":786,"tokens_out":875,"would_cite":true,"duration_ms":43489,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83D05","83F05"],"pacs":["04.50.Kd","95.36.+x","98.80.-k"],"model":"deepseek-v4-flash","headline":"A two-parameter exponential modification of nonmetric gravity can mimic dark energy and change how cosmic structure grows, while staying close to ΛCDM.","keywords":["f(Q) gravity","nonmetricity","dark energy","quintessence","phantom","growth index","fσ8","energy conditions"],"falsifier":"Compute the full second-order action for $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$ in the parameter range used here and check for ghosts or strong coupling; if an instability appears, the sub-horizon growth equation is not trustworthy. Observationally, a measurement of $f\\sigma_8(z)$ at several redshifts with errors smaller than the model's percent-level deviations from $\\Lambda$CDM would distinguish the $b>0$ and $b<0$ branches.","tokens_in":26073,"feed_emoji":"🌌","tokens_out":7623,"duration_ms":75608,"temperature":0.7,"pith_summary":"The paper studies a modified theory of gravity in which the usual action is altered by an exponential term involving the non-metricity scalar $Q$. The model $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$ is a smooth small perturbation of $\\Lambda$CDM, and the paper derives analytic expressions for the Hubble parameter at late times for $n=1$ and $n=2$. Using these, it calculates the effective dark-energy equation of state, energy conditions, deceleration and statefinder parameters, and the growth of linear matter perturbations. The central conclusion is that the model is a workable dark-energy surrogate: it produces late-time accelerated expansion, behaves like quintessence for $b>0$ and like phantom energy for $b<0$, and its growth predictions, including $f\\sigma_8$, are consistent with observational data within the quoted errors. A sympathetic reader should care because the model offers a geometric origin for dark-energy phenomenology without new fields, and it makes a concrete prediction that growth observables deviate from $\\Lambda$CDM near the present.","feed_headline":"One exponential term in gravity mimics dark energy and alters growth","feed_subtitle":"A small deviation from ΛCDM changes how matter clusters and matches observed fσ8 data.","key_machinery":"The central object is the exponential nonmetricity function $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$, a smooth perturbative deformation of the $\\Lambda$CDM action. Expanding in the small parameter $b$ turns the Friedmann equation into an algebraic equation for $u=H^2/H_0^2$, which is solved order by order to give analytic $H(z)$ for $n=1$ and $n=2$. The same expansion controls the effective dark-energy density and pressure, and its derivative $f_Q$ enters the perturbation equation through $G_{\\rm eff}=G/f_Q$, which is what carries the model's influence on structure growth.","core_discovery":"The paper argues that the exponential nonmetricity model $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$, which reduces to $\\Lambda$CDM when $b=0$, is a viable late-time cosmology. With $|b|$ small it yields analytic Hubble parameters $H(z)$ for $n=1$ and $n=2$ that stay within about $0.008\\%$ of the numerical solution. The effective dark energy arising from the geometry then behaves as quintessence for $b>0$ ($w_{DE}>-1$) and as phantom-like for $b<0$ ($w_{DE}<-1$), while the effective gravitational constant $G_{\\rm eff}=G/f_Q$ changes near the present, altering the linear growth of matter perturbations, the growth index $\\gamma$, and the observable $f\\sigma_8$. The model's $f\\sigma_8$ curves fall within the error bars of 23 published growth measurements, and the authors take this as showing that the model reproduces the main effects attributed to dark energy.","pith_inferences":["The authors assign $b$ by hand in the growth plots; a full statistical fit to the $f\\sigma_8$ compilation could tighten the model and reveal whether the preferred branch is the quintessence or phantom side.","The perturbation-sector caveat the authors cite implies a decisive check: computing the full action-level kinetic structure for this $f(Q)$ could determine whether ghosts or strong coupling appear at the parameter values used, in which case the background and energy-condition results would survive but the growth results would need revision.","Because the model's effective gravitational constant changes with redshift, the same exponential ansatz could be tested with cosmic-shear or CMB-lensing statistics, where a time-varying gravitational strength leaves distinctive scale-dependent signatures."],"forward_implications":["If the model is correct, late-time cosmic acceleration can be produced by the geometry of nonmetric gravity alone, with no additional scalar field or separate dark-energy fluid.","The sign of the deviation parameter $b$ controls the effective dark-energy sector: $b>0$ gives quintessence-like behavior with $w_{DE}>-1$, while $b<0$ gives phantom-like behavior with $w_{DE}<-1$.","Because $G_{\\rm eff}=G/f_Q$ varies with redshift, structure growth is modified near the present; the model predicts different clustering compared with $\\Lambda$CDM for the two signs of $b$.","The $n=2$ variant differs from $\\Lambda$CDM only at order $b^2$ in the background, so its expansion history is nearly indistinguishable from $\\Lambda$CDM and growth data become the main way to detect it."],"supporting_citations":[{"why":"Introduces the exponential $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$ model and provides the MCMC constraints on $b$ from $H(z)$ data that this paper extends.","marker":"[24]"},{"why":"Supplies the perturbative expansion in $b$ used to solve the Friedmann equation analytically around the $\\Lambda$CDM solution.","marker":"[31]"},{"why":"Provides the sub-horizon linear matter perturbation equation whose source term is proportional to $G/f_Q$.","marker":"[48]"},{"why":"Gives the growth-index formulation $f_g=\\Omega_m^{\\gamma}$ used to compute $\\gamma(z)$ and compare with the $\\Lambda$CDM value $6/11$.","marker":"[49]"},{"why":"Provides the 'Gold-2017' compilation of growth-rate data that is included in the $f\\sigma_8$ comparison.","marker":"[57]"},{"why":"Tests the internal robustness of the growth-rate dataset used for the $f\\sigma_8$ comparison.","marker":"[56]"},{"why":"Raises the strong-coupling and ghost pathologies of $f(Q)$ modifications that the paper cites as a caveat for the perturbation-sector results.","marker":"[61]"},{"why":"Supplies the Planck values of $\\Omega_{m,0}$, $H_0$, and $\\sigma_8$ used in the numerical evaluations.","marker":"[32]"}],"fun_headline_variants":["Exponential gravity term mimics dark energy, alters growth","Tiny exponential tweak to gravity shifts cosmic structure growth","Modified gravity with exponential term matches fσ8 data","Exponential f(Q) model reproduces dark energy and growth effects"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The growth-of-structure part assumes the standard sub-horizon perturbation equation of $f(Q)$ gravity, with the gravitational constant replaced by $G/f_Q$, remains valid for this model; if the strong-coupling or ghost problems noted in the paper invalidate that equation, the growth predictions would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Exponential gravity term mimics dark energy, alters growth","Tiny exponential tweak to gravity shifts cosmic structure growth","Modified gravity with exponential term matches fσ8 data","Exponential f(Q) model reproduces dark energy and growth effects"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000741,"raw_usage":{"total_tokens":3314,"prompt_tokens":960,"completion_tokens":2354,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":2287}},"tokens_in":576,"tokens_out":2354,"duration_ms":18676,"temperature":1.0,"reasoning_tokens":2287,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T17:01:43.685207+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full second-order action for $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$ in the parameter range used here and check for ghosts or strong coupling; if an instability appears, the sub-horizon growth equation is not trustworthy. Observationally, a measurement of $f\\sigma_8(z)$ at several redshifts with errors smaller than the model's percent-level deviations from $\\Lambda$CDM would distinguish the $b>0$ and $b<0$ branches.","supporting_citations":[{"cited_title":"Oliveros and M","cited_arxiv_id":null,"evidence_quote":"Introduces the exponential $f(Q)=Q+2\\Lambda e^{-(b\\Lambda/Q)^n}$ model and provides the MCMC constraints on $b$ from $H(z)$ data that this paper extends."},{"cited_title":"Basilakos, S","cited_arxiv_id":null,"evidence_quote":"Supplies the perturbative expansion in $b$ used to solve the Friedmann equation analytically around the $\\Lambda$CDM solution."},{"cited_title":"Beltr´ an Jim´ enez, L","cited_arxiv_id":null,"evidence_quote":"Provides the sub-horizon linear matter perturbation equation whose source term is proportional to $G/f_Q$."},{"cited_title":"Khyllep, et al","cited_arxiv_id":null,"evidence_quote":"Gives the growth-index formulation $f_g=\\Omega_m^{\\gamma}$ used to compute $\\gamma(z)$ and compare with the $\\Lambda$CDM value $6/11$."},{"cited_title":"Nesseris, G","cited_arxiv_id":null,"evidence_quote":"Provides the 'Gold-2017' compilation of growth-rate data that is included in the $f\\sigma_8$ comparison."},{"cited_title":"Sagredo, S","cited_arxiv_id":null,"evidence_quote":"Tests the internal robustness of the growth-rate dataset used for the $f\\sigma_8$ comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Raises the strong-coupling and ghost pathologies of $f(Q)$ modifications that the paper cites as a caveat for the perturbation-sector results."},{"cited_title":"Aghanim et al","cited_arxiv_id":null,"evidence_quote":"Supplies the Planck values of $\\Omega_{m,0}$, $H_0$, and $\\sigma_8$ used in the numerical evaluations."}],"review_version":1}