Pith. sign in

REVIEW 3 major objections 4 minor 26 references

Inflationary observables in $F(R)$ gravity

T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A quartic f(R) gravity model with fixed cubic and quartic couplings produces a spectral index and tensor-to-scalar ratio consistent with Planck and BICEP/Keck, and two of its solutions favor r below 0.0015.

desk verdict The quartic f(R) paper has an internal error in its spectral index formula that invalidates its headline low-r claim; the hand-picked parameters only compound the problem. read the letter →

arxiv 2412.19030 v2 pith:I6NWINAA submitted 2024-12-26 gr-qc hep-ph

Enrique Diaz , Oscar Meza-Aldama This is my paper · ORCID
classification gr-qchep-ph PACS 98.80.Cq
keywords f(R)gravityinflationarycosmologytensor-to-scalarratiospectralindexslow-rollapproximationEinsteinframequarticcurvaturecorrectionsCOBEnormalization
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that adding cubic and quartic Ricci-scalar terms to the usual $R+R^2$ inflation action produces a model that survives current data. With the ratios $x=\beta/\alpha=2\times10^{-8}$ and $y=\gamma/\alpha=10^{-8}$, the model's spectral index $n_s$ and tensor-to-scalar ratio $r$ fall inside the Planck 2018 and BICEP/Keck 2021 bounds. Two of the four mathematical solutions for the Einstein-frame scalar field are compatible with these bounds; they predict $0.0005

What carries the argument

The load-bearing mechanism is the Einstein-frame mapping of $f(R)$ gravity: an auxiliary field $\phi$ with $f'(\phi)=\sigma=e^{\sqrt{2/3}s}$ converts the quartic action into a minimally coupled scalar with potential $V_E(s)=(\phi f'(\phi)-f(\phi))/(2f'(\phi)^2)$. The e-folds integral $N=\sqrt{3/2}\int (V/(\sigma V'))\,d\sigma$ becomes a quartic polynomial in $\sigma$, so $\sigma(N)$ has four roots; selecting the two roots that keep the slow-roll parameters within observed bounds gives the paper's predictions. The named object is the quartic inversion of the e-fold relation, whose two accepted branches are called $\sigma_1$ and $\sigma_2$.

What would settle it

A future CMB B-mode experiment that measures $r>0.0015$ at 95% confidence would falsify the $\sigma_1$ branch's predicted window, and $r>0.00015$ would falsify the $\sigma_2$ branch. Recomputing $n_s$ and $r$ at other values of $x$ and $y$ (for example, $x=y=10^{-7}$) and finding that the agreement with Planck disappears would show that the central result is tied to the chosen parameter point rather than to the structure of the quartic model.

Watch

Extended reading notes

Core claim

The central discovery is that a quartic polynomial $f(R)=R+\tfrac12\alpha R^2+\tfrac13\beta R^3+\tfrac14\gamma R^4$, after conformal transformation to the Einstein frame, yields slow-roll observables that match current bounds for a specific parameter choice. The field redefinition $\sigma=e^{\sqrt{2/3}s}=1+\alpha\phi+\beta\phi^2+\gamma\phi^3$ converts the potential into a form where the number of e-folds is a quartic polynomial in $\sigma$; inverting it gives four branches $\sigma_i(N)$, two of which produce $n_s$ and $r$ within Planck and BICEP/Keck constraints at $N=60$. The compatible branches predict very low tensor-to-scalar ratios, and the COBE normalization fixes $\alpha\approx1.187\times10^9$. The results are quantitatively similar to earlier $R^3$ models, suggesting the quartic term does not dramatically alter the inflationary predictions.

Load-bearing premise

The load-bearing premise is the un-derived parameter choice $x=\beta/\alpha=2\times10^{-8}$ and $y=\gamma/\alpha=10^{-8}$; the two low-tensor solutions are computed for that point, and nothing in the paper derives or observationally motivates it.

Editorial extensions

If this is right

  • If the model is correct, the tensor-to-scalar ratio should be found between $0.0005$ and $0.0015$ on the $\sigma_1$ branch, or below $0.00015$ on the $\sigma_2$ branch, when the spectral index is within the Planck 2018 value.
  • The fixed value $\alpha\approx1.187\times10^9$ from the COBE normalization gives a concrete scalar amplitude prediction that future CMB measurements can check.
  • Because the results are quantitatively similar to earlier $R^3$ models, adding the quartic term is a mild deformation: it does not rescue a ruled-out model, nor does it spoil an allowed one.
  • The compatible branches already sit below the expected $0.001$ sensitivity of next-generation CMB experiments, so a detection of $r$ in this range would support the quartic model, while a null at the $0.001$ level would begin to disfavor the $\sigma_1$ branch.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The values $x=2\times10^{-8}$ and $y=10^{-8}$ are chosen without derivation; the low-$r$ predictions are conditional on this single parameter point. A systematic scan over $x$ and $y$ would show whether the agreement with Planck and BICEP/Keck is a broad region or a fine-tuned slice.
  • If the model is taken literally, the two branches $\sigma_1$ and $\sigma_2$ give mutually exclusive predictions for $r$ (one near $0.001$, the other far below). Distinguishing them observationally would require a measurement with sensitivity below $0.00015$, beyond currently planned single-experiment sensitivity; a joint analysis or an additional observable might break the degeneracy.
  • The same analytic machinery could be applied to quintic or higher polynomial $f(R)$, but the e-fold relation would no longer be quartic in $\sigma$, so the exact four-root inversion would need to be replaced by numerical root finding. Observing whether the low-$r$ property survives that generalization would test whether it is intrinsic to higher-order curvature corrections or specific to the quart
  • The COBE normalization is used only to fix $\alpha$ for the chosen $x$ and $y$; one could instead use the Planck measurement of the scalar amplitude to derive a joint posterior over $(\alpha,x,y)$ and see whether the values that fit the spectral index also reproduce the amplitude.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the quartic f(R) model f(R)=R+(1/2)αR^2+(1/3)βR^3+(1/4)γR^4, transforms it to the Einstein frame via the conformal transformation (4), computes the slow-roll parameters (17)–(18), inverts the e-fold expression (22) to obtain two roots σ1 and σ2, fixes x=β/α=2×10^-8 and y=γ/α=10^-8, and evaluates the spectral index n_s and tensor-to-scalar ratio r at N=60. It claims that these observables satisfy Planck/BICEP bounds and that σ1 and σ2 favor very low r values, and it uses a COBE normalization to fix α.

Significance. If the reported predictions were correct, the model would be a concrete quartic f(R) extension of Starobinsky inflation that is testable with LiteBIRD. The paper does not provide machine-checked proofs, code, or a parameter-free derivation; its main quantitative results are obtained from analytic expressions evaluated at hand-selected parameter values. Because the central observable n_s is computed with an incorrect slow-roll formula, and the parameter choice is not derived from any principle, the claimed agreement with Planck/BICEP is not established.

major comments (3)
  1. [Section II, Eq. (11)] With the definitions (8)–(9), the standard slow-roll expression is n_s = 1 − 6ϵ + 2η, not n_s = 1 − 2ϵ − η. In the Starobinsky limit at N=60, the paper's own expressions give ϵ ≈ 1.4×10^-4 and η ≈ −0.014, so Eq. (11) yields n_s ≈ 1.013 whereas the correct formula gives n_s ≈ 0.972. All reported n_s values, the figures, and the abstract's claim of consistency with Planck are therefore based on an invalid formula. This error is internal to the paper's own definitions and is not a matter of convention.
  2. [Section II, Eqs. (4)–(7)] The conformal transformation (4) is written as g_{μν} → (1/f'(φ)) g_{μν}, but with the scalar field defined as s = √(3/2) ln f'(φ) in Eq. (6), the Einstein-frame metric should be g_{μν} → f'(φ) g_{μν}. The written version is the inverse, and since V_E(s) is never displayed, the reader cannot check whether the slow-roll expressions (17)–(18) are actually derived from a correctly transformed potential. This inconsistency needs to be resolved and the potential V_E(σ) shown explicitly.
  3. [Section IV and Appendix A] The values x=2×10^-8 and y=10^-8 are introduced without derivation, motivation, or a scan over allowed values. The low tensor-to-scalar ratios reported for σ1 and σ2 (0.0005<r_σ1<0.0015 and r_σ2<0.00015 in the abstract) are direct outputs of this choice; they are therefore fitted numbers, not predictions. A claim that the model favors very low r requires a derivation of x and y from some condition, for example from a consistency relation or a prior over parameters.
minor comments (4)
  1. [Introduction and Eq. (13)] The coefficient of βR^3 is written as 1/2 in the introduction but as 1/3 in Eq. (13); this inconsistency in the definition of the model should be fixed.
  2. [Section IV, Eq. (23)] Eq. (23) states the COBE normalization V/ϵ = 0.0274, but the standard amplitude A_s = V/(24π^2 ϵ) ≈ 2.1×10^-9 implies V/ϵ ≈ 5×10^-7. The quoted normalization (and hence the derived value α = 1.187×10^9 in Eq. (24)) needs to be justified or corrected; as it stands, this step is unexplained.
  3. [Section III and Appendix A] Eq. (16) and Eq. (25) are obtained by 'ignoring higher order terms in x and y' without an estimate of the introduced error; given the tiny values of x and y this is likely acceptable, but the approximation should be quantified.
  4. [Figures 1 and 2] The figures are impossible to reproduce from the text: the caption of Fig. 1 says 'we explored different values of 0<α≤1' while Fig. 2 is described as generated with the SpaceMath package [12], but the precise axes, parameter grids, and root-selection criteria are not specified.

Circularity Check

1 steps flagged · score 6.0 of 10

Central Planck/BICEP compatibility and low-r claims reduce to hand-set x,y and post-hoc rejection of two non-compliant roots.

  1. fitted input called prediction [Section IV, paragraph after Eq. (22) and before Eq. (23)]
    "Throughout this paper we use the values x = 2 × 10−8 and y = 10 −8. Since (22) is a quartic polynomial, we can invert it to get σ in terms of N ... Two of the solutions did not satisfy current bounds for the measured slow roll parameters, hence we only work with two of the roots which are presented in the appendix."

    The abstract's claim that the quartic f(R) model satisfies Planck/BICEP bounds and that σ1,σ2 favor very low r is not an independent prediction. The ratios x=β/α and y=γ/α are assigned by hand without derivation, and any algebraic root that failed the target bounds is explicitly removed before presenting the surviving observables. Thus the surviving σ1,σ2 branches are retained precisely because they satisfy the experimental bounds; the reported compatibility and low r values are a selection effect rather than a consequence forced by the f(R) action. This is the pattern of a tuned input or post-selected output being presented as a phenomenological prediction.

full rationale

The core calculation—transforming the quartic f(R) action to the Einstein frame, deriving the slow-roll parameters, and integrating the e-fold relation—is self-contained algebra, and the COBE normalization fixing α=1.187×10^9 is standard amplitude matching rather than circularity. No load-bearing self-citation was found: the cited [5] and [13] are not used to justify the central result, and the paper is not invoking a uniqueness theorem or smuggling in an ansatz via citation. However, the paper's central phenomenological claim does reduce to its own inputs. The parameters x and y are simply asserted ('we use the values...'), and the text then states that two of the four inversion branches 'did not satisfy current bounds... hence we only work with two of the roots.' That procedure guarantees that the presented branches satisfy Planck/BICEP; reporting their low tensor-to-scalar ratios as 'favoring' LiteBIRD sensitivity is therefore a post-selection effect, not an independent prediction. This is partial circularity of the fitted/tuned-input kind, scored 6. The separate and serious Eq. (11) problem—defining n_s ≡ 1 − 2ϵ − η instead of the standard slow-roll n_s = 1 − 6ϵ + 2η—is a correctness and internal-consistency failure that invalidates the Planck comparison, but it is not itself input-output circularity, so it does not further raise the circularity score.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The central claim relies on two hand-picked parameters (x and y), a COBE-fitted scale α, and several standard approximations. No new entities are introduced, but the free parameters are tuned to match the data, and the domain assumptions (N=60, σ≫σe) are not scrutinized.

free parameters (3)
  • x = β/α = 2×10⁻⁸
    Chosen by hand in Section IV; no derivation from first principles or independent data. It controls the shape of the potential and hence the low-r prediction.
  • y = γ/α = 10⁻⁸
    Chosen by hand in Section IV; no derivation. It modifies the quartic term and, together with x, sets the predicted observables.
  • α = 1.187×10⁹
    Fixed by the COBE normalization of the power spectrum (Eq. 23), i.e., fitted to the measured amplitude of perturbations. It is not predicted independently.
assumptions (4)
  • standard math Slow-roll approximation and the standard slow-roll definitions of ns and r (Eqs. 8–12)
    Used throughout to relate the potential to observables; standard in inflationary cosmology.
  • standard math The conformal transformation to the Einstein frame yields a minimally coupled scalar with potential (7)
    Standard method in f(R) gravity, though the transformation in Eq. (4) appears misstated.
  • domain assumption N = 60 e-folds at horizon exit
    Standard cosmological assumption to solve the flatness problem; the results are sensitive to this choice.
  • ad hoc to paper The inversion (16) keeps only up to linear terms in x and y, and the integral (22) assumes σ ≫ σe
    These approximations are not quantitatively justified and affect the accuracy of the final expressions; they are introduced to make the calculation tractable.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Inflationary observables in $F(R)$ gravity." pith.science (2026). https://pith.science/paper/I6NWINAA

@misc{pith2026241219030,
  author       = {Pith},
  title        = {Pith review of: Inflationary observables in $F(R)$ gravity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/I6NWINAA}},
  note         = {Machine review of arXiv:2412.19030}
}
abstract

We present phenomenological signatures for a modified gravity model f(R), constructed with linear, quadratic, cubic and quartic terms. The obtained signatures satisfy current phenomenological bounds reported by PLANCK and BICEP3. Furthermore, two of the model solutions $\sigma_1$ and $\sigma_2$ seem to favor a much lower value for the tensor-to-scalar ratio $0.0005<r_{\sigma_1}<0.0015$ and $r_{\sigma_2}<0.00015$ than the current reported experimental bounds. The results we obtained are quantitatively similar to those presented in previous studies for $R^3$ models.

Figures

Figures reproduced from arXiv: 2412.19030 by the authors.

Figure 1
Figure 1. FIG. 1: The blue line corresponds to the solution [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Fits corresponding to solution [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 13 canonical work pages

  1. [14]

    Arroyo-Ure˜ na, R

    M.A. Arroyo-Ure˜ na, R. Gait´ an and T.A. Valencia-P´ erez, SpaceMath version 1.0, a Mathematica package for beyond the Standard Model parameter space searches, Rev. Mex. Fis. E 19 (2022) 020206 [arXiv:2008.00564] [INSPIRE]

    arXiv 2022

  2. [1]

    This plot was generated using the package presented in [12]

    We also present a scatter plot of allowed values for r and ns in figure 2. This plot was generated using the package presented in [12]. The number of e-folds is N ≡ Z s se 1√ 2ϵ ds; (19) upon changing variables from s to σ σ = e √ 2 3 s, (20) we get N = r 3 2 Z σ σe V (σ) σV ′(σ) dσ, (21) and assuming σ ≫ σe, we obtain the expression N = r 3 2   (−36y +...

    work page

  3. [2]

    in which we added higher-order Ricci scalar terms as opposed to just having the linear R term given in Einstein gravity. Several f (R) models, such as the Starobinsky model f (R) = R + αR2, have been extensively studied and are known to produce predictions consistent with the latest observational data from the Planck satellite, particularly the spectral i...

    work page Pith review arXiv 2025

  4. [3]

    Guth, Inflationary universe: A possible solution to the horizon and flatness problems , Phys

    A. Guth, Inflationary universe: A possible solution to the horizon and flatness problems , Phys. Rev. D, 23 (1981) pg.347

    work page 1981

  5. [4]

    De Felice and S

    A. De Felice and S. Tsujikawa, f(R) theories, Living Rev. Rel. 13 (2010) 3 [1002.4928]

    arXiv 2010

  6. [5]

    Starobinsky, A New Type of Isotropic Cosmological Models Without Singularity Phys

    A.A. Starobinsky, A New Type of Isotropic Cosmological Models Without Singularity Phys. Lett. B91 99 (1980)

    work page 1980

  7. [6]

    Hodges, G.R

    H.M. Hodges, G.R. Blumenthal, L.A. Kofman and J.R. Primack, Nonstandard Primordial Fluctuations From a Polynomial Inflaton Potential , Nucl. Phys. B 335 (1990) 197

    work page 1990

  8. [7]

    D ´ ıaz E and Meza-Aldama O 2020 J. Phys. Commun. 4 125011

    work page 2020

Show all 26 references
  1. [8]

    S. D. Odintsov and V. K. Oikonomou, Phys. Rev. D 99 (2019) no.6, 064049 doi:10.1103/PhysRevD.99.064049 [arXiv:1901.05363 [gr-qc]]

    2019 arXiv

  2. [9]

    Inflation (2023)

    Ellis, J.; Wands, D. Inflation (2023). arXiv 2023, arXiv:2312.13238v1

    2023

  3. [10]

    Planck collaboration, Planck 2018 results. X. Constraints on inflation, Astron. Astrophys. 641 (2020) A10

    2020

  4. [11]

    Cosmology in Scalar-Tensor Gravity, 2004

    Faraoni, V. Cosmology in Scalar-Tensor Gravity, 2004. Springer

    2004

  5. [12]

    P. A. R. Ade et al. [BICEP and Keck], Phys. Rev. Lett. 127, no.15, 151301 (2021)

    2021

  6. [13]

    Dodelson S, 2003, Modern Cosmology, (NewYork:Academic)

    2003

  7. [15]

    Beyond the Starobinsky model for inflation

    Dhong Yeon Cheong, Hyun Min Lee, and Seong Chan Park. Beyond the Starobinsky model for inflation. Phys. Lett. B, 805:135453, 2020

    2020

  8. [16]

    The inflationary scenario in the f (R) gravity model with a R4 term

    Sahazada Aziz, Sohan Kumar Jha, Anisur Rahaman. The inflationary scenario in the f (R) gravity model with a R4 term. 2021 Class. Quantum Grav. 38 225008, arXiv:2211.08159 [gr-qc]

    2021 arXiv

  9. [17]

    Nojiri and S

    S. Nojiri and S. D. Odintsov, Phys. Rept. 505 (2011), 59-144 doi:10.1016/j.physrep.2011.04.001 [arXiv:1011.0544 [gr-qc]]

    2011 arXiv

  10. [18]

    Nojiri, S

    S. Nojiri, S. D. Odintsov and V. K. Oikonomou, Phys. Rept. 692 (2017), 1-104 doi:10.1016/j.physrep.2017.06.001 [arXiv:1705.11098 [gr-qc]]

    2017 arXiv

  11. [19]

    S. D. Odintsov, V. K. Oikonomou, I. Giannakoudi, F. P. Fronimos and E. C. Lymperiadou, Symmetry 15 (2023) no.9, 1701 doi:10.3390/sym15091701 [arXiv:2307.16308 [gr-qc]]

    2023 arXiv

  12. [20]

    Fuskeland et al

    U. Fuskeland et al. (LiteBIRD), Astron. Astrophys. 676, A42 (2023), arXiv:2302.05228 [astro-ph.CO]. 6

    2023 arXiv

  13. [21]

    S. D. Odintsov, V. K. Oikonomou and L. Sebastiani, Nucl. Phys. B923 (2017), 608-632 doi:10.1016/j.nuclphysb.2017.08.018 [arXiv:1708.08346 [gr-qc]]

    2017 arXiv

  14. [22]

    Elizalde, S

    E. Elizalde, S. D. Odintsov, L. Sebastiani and R. Myrzakulov, Nucl. Phys. B 921 (2017), 411-435 doi:10.1016/j.nuclphysb.2017.06.003 [arXiv:1706.01879 [gr-qc]]

    2017 arXiv

  15. [23]

    Myrzakulov, S

    R. Myrzakulov, S. Odintsov and L. Sebastiani, Phys. Rev. D 91 (2015) no.8, 083529 doi:10.1103/PhysRevD.91.083529 [arXiv:1412.1073 [gr-qc]]

    2015 arXiv

  16. [24]

    S. D. Odintsov and V. K. Oikonomou, Phys. Rev. D 104 (2021) no.12, 124065 doi:10.1103/PhysRevD.104.124065 [arXiv:2112.06269 [gr-qc]]

    2021 arXiv

  17. [25]

    S. D. Odintsov and V. K. Oikonomou, Phys. Lett. B 833 (2022), 137353 doi:10.1016/j.physletb.2022.137353 [arXiv:2206.06024 [gr-qc]]

    2022

  18. [26]

    Sebastiani, G

    L. Sebastiani, G. Cognola, R. Myrzakulov, S. D. Odintsov and S. Zerbini, Phys. Rev. D 89 (2014) no.2, 023518 doi:10.1103/PhysRevD.89.023518 [arXiv:1311.0744 [gr-qc]]

    2014 arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.