arxiv: 2605.10728 · v1 · submitted 2026-05-11 · 🌌 astro-ph.CO

Recognition: no theorem link

Numerical Study of Some Generalizations of the Starobinsky Inflationary Model

Jordan Zambrano , \'Angel Rodr\'iguez , Melisa Alvear , Clara Rojas

Authors on Pith no claims yet

Pith reviewed 2026-05-12 04:42 UTC · model grok-4.3

classification 🌌 astro-ph.CO

keywords Starobinsky inflationalpha-Starobinsky modelpower-law inflationPlanck 2018numerical power spectraslow-roll parametersinflationary observablescosmic microwave background

0 comments

The pith

Numerical scans show that alpha and power-law alpha extensions of the Starobinsky model fit Planck 2018 data for specific parameter choices.

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

The paper examines three inflationary scenarios that generalize the classic Starobinsky model: the alpha version, a power-law version, and their combination. For each, the authors derive slow-roll formulas, compute the scalar and tensor power spectra numerically across parameter ranges, and map the resulting observables onto the planes used by cosmologists to compare with data. They find that the first and third scenarios can be brought into agreement with Planck 2018 measurements when parameters are chosen appropriately, while still recovering the original Starobinsky behavior at certain limits. This matters because it identifies concrete ways to adjust the shape of the inflationary potential without losing consistency with current observations of the cosmic microwave background.

Core claim

By deriving the relevant slow-roll expressions and numerically evaluating the scalar and tensor power spectra over the corresponding parameter spaces, the analysis shows that the alpha-Starobinsky model and the power-law alpha-Starobinsky model are favored by Planck 2018 observations for certain choices of parameters, while each scenario reproduces the standard Starobinsky limit when parameters are set appropriately.

What carries the argument

Numerical computation of scalar and tensor power spectra from slow-roll expressions evaluated across the parameter space of the three generalized Starobinsky models.

If this is right

The alpha-Starobinsky model produces inflationary observables consistent with Planck 2018 for appropriate parameter values.
The power-law alpha-Starobinsky model is likewise favored by the same data for suitable choices of its parameters.
All three models recover the standard Starobinsky inflationary predictions when their extra parameters are set to recover the original potential.
Comparative contour plots in the (n_s, A_s) and (r, n_s) planes quantify the regions of parameter space compatible with current observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The numerical approach could be extended to compute higher-order statistics such as the bispectrum to test whether the favored parameter regions also satisfy future non-Gaussianity bounds.
Because the models remain viable only inside limited parameter windows, they provide concrete targets for next-generation CMB experiments that aim to tighten constraints on the tensor-to-scalar ratio.
The same scan technique could be applied to other single-field potentials to map which functional forms remain viable once Planck-level precision is imposed.

Load-bearing premise

The slow-roll approximations remain accurate across the full scanned parameter space and the numerical power spectra capture the dynamics without higher-order corrections or initial-condition sensitivities.

What would settle it

A future measurement of the tensor-to-scalar ratio that lies systematically outside the contour regions obtained for the favored parameter choices in the (r, n_s) plane would falsify the claim that those parameter choices are consistent with observations.

read the original abstract

In this work, we perform a numerical study of three Starobinsky--type inflationary scenarios: the $\alpha$--Starobinsky inflationary model, the power--law Starobinsky inflationary model, and the power--law $\alpha$--Starobinsky inflationary model. For an appropriate choice of parameters, each scenario reproduces the standard Starobinsky limit. For each case, we derive the relevant slow--roll expressions in order to compute numerically the scalar and tensor power spectra over the corresponding parameter space and evaluate the associated inflationary observables. Finally, we provide a comparative analysis in the $(n_\sca,A_\sca)$ and $(r,n_\sca)$ parameter spaces using contour plots. Our results indicate that, for certain choices of parameters, the $\alpha$--Starobinsky model and the power--law $\alpha$--Starobinsky model are favored by \textit{Planck} 2018 observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This is a standard numerical parameter scan of known Starobinsky generalizations that finds some values fit Planck contours but skips checks on whether slow-roll holds at those points.

read the letter

The main thing to know is that the paper scans parameters in the alpha-Starobinsky model, the power-law Starobinsky model, and their combination. It derives the slow-roll expressions, computes the scalar and tensor spectra numerically over the space, and plots the results in the usual (ns, As) and (r, ns) planes. For certain parameter choices the points fall inside the Planck 2018 contours, and the authors note that each case can recover the original Starobinsky limit when the extra parameters are set to one or zero. The contour comparisons are the clearest part of the work and give a direct visual sense of how the three scenarios overlap with data relative to each other. That is useful if you are tracking how these potentials behave under current constraints. The soft spots are more important than the plots suggest. The claim that the alpha and power-law alpha versions are favored rests on the slow-roll formulas being accurate, yet the paper gives no evidence that epsilon and eta remain much less than one at horizon exit for the selected points, nor does it mention any cross-check with full Mukhanov-Sasaki integration. Without that, the contour comparison can be misleading. The exercise is also post-hoc fitting: parameters are tuned to reproduce Starobinsky and then scanned until the observables land inside the allowed region. This is common in the field but does not constitute a new prediction. These model extensions already appear in earlier literature, so the numerical exercise itself is the only addition. This kind of paper is mainly for readers who maintain broad catalogs of inflationary models and their observational status. It does not introduce new analytic results, resolve open questions, or provide independently falsifiable predictions. Given the missing verification on the approximation and the low novelty, I would not send it to peer review.

Referee Report

2 major / 2 minor

Summary. The paper performs a numerical study of three Starobinsky-type inflationary scenarios (the α-Starobinsky model, the power-law Starobinsky model, and the power-law α-Starobinsky model). For appropriate parameter choices each recovers the standard Starobinsky limit. Slow-roll expressions for the scalar and tensor power spectra are derived and evaluated numerically over the respective parameter spaces; the resulting observables (n_s, A_s, r) are then compared to Planck 2018 constraints via contour plots in the (n_s, A_s) and (r, n_s) planes. The central claim is that, for certain parameter values, the α-Starobinsky and power-law α-Starobinsky models are favored by the data.

Significance. If the slow-roll regime is explicitly verified to hold for the parameter points that lie inside the observational contours and if the numerical evaluation of the spectra is cross-checked against the full Mukhanov-Sasaki equations, the work would supply a systematic comparative map of how the extra parameters in these generalizations shift the predictions relative to the original Starobinsky model. Such maps are useful for model-building and for guiding future CMB constraints, provided the underlying approximations remain valid.

major comments (2)

[Abstract] Abstract: The assertion that the α–Starobinsky and power–law α–Starobinsky models are favored by Planck 2018 for certain parameter choices rests on the numerically computed (n_s, A_s) and (r, n_s) values falling inside the observational contours. The abstract gives no indication that the slow-roll conditions (ε ≪ 1 and |η| ≪ 1 at horizon exit, remaining small over the relevant ~50–60 e-folds) were verified for those specific points, nor that a full numerical integration of the Mukhanov-Sasaki equations was performed as a cross-check. Without this verification the contour comparison is not demonstrably reliable.
[Abstract] Abstract and numerical-study description: The parameters are chosen so that each scenario reproduces the standard Starobinsky limit, after which the parameter space is scanned and values are selected that place the predictions inside the Planck contours. This procedure risks post-hoc fitting rather than an a-priori test; the manuscript should clarify whether any of the reported “favored” regions were identified before the observational comparison or whether they emerge only after the scan.

minor comments (2)

The notation n_ sca and A_ sca (and similarly for the tensor quantities) should be written consistently as n_s and A_s throughout; the current subscript is unclear on first reading.
[Abstract] The abstract states that all three models are studied but the final comparative analysis and the “favored” claim are stated only for the α–Starobinsky and power–law α–Starobinsky cases. Clarify whether the pure power-law Starobinsky model is also shown in the contour plots and, if not, why it is omitted from the final assessment.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable suggestions. We address the major comments point by point below, indicating the changes we will implement in the revised manuscript.

read point-by-point responses

Referee: [Abstract] Abstract: The assertion that the α–Starobinsky and power–law α–Starobinsky models are favored by Planck 2018 for certain parameter choices rests on the numerically computed (n_s, A_s) and (r, n_s) values falling inside the observational contours. The abstract gives no indication that the slow-roll conditions (ε ≪ 1 and |η| ≪ 1 at horizon exit, remaining small over the relevant ~50–60 e-folds) were verified for those specific points, nor that a full numerical integration of the Mukhanov-Sasaki equations was performed as a cross-check. Without this verification the contour comparison is not demonstrably reliable.

Authors: We agree with the referee that explicit verification of the slow-roll conditions for the parameter points inside the observational contours is necessary to ensure the reliability of our results. In the original manuscript, we derived the slow-roll expressions and used them to compute the observables numerically, but we did not include explicit checks of ε and η for the selected points. We will revise the manuscript to include such verifications, for example by plotting or tabulating the values of ε and η at horizon exit for the points that lie within the Planck contours, confirming that they satisfy ε ≪ 1 and |η| ≪ 1. Regarding the cross-check with the full Mukhanov-Sasaki equations, our study focuses on the slow-roll approximation as is common in such numerical explorations of inflationary models. Performing a full numerical integration would require significant additional computational effort and is outside the scope of this work. However, we will add a discussion noting the expected accuracy of the slow-roll approximation for these models based on existing literature, and we will state that the results should be interpreted within the slow-roll regime. revision: partial
Referee: [Abstract] Abstract and numerical-study description: The parameters are chosen so that each scenario reproduces the standard Starobinsky limit, after which the parameter space is scanned and values are selected that place the predictions inside the Planck contours. This procedure risks post-hoc fitting rather than an a-priori test; the manuscript should clarify whether any of the reported “favored” regions were identified before the observational comparison or whether they emerge only after the scan.

Authors: We appreciate the referee highlighting this potential issue. The parameter ranges were selected to include the standard Starobinsky limit as a special case, which is a standard practice when studying generalizations of known models. The numerical scan was then performed over these ranges to compute the observables for various parameter combinations. The regions reported as consistent with Planck 2018 data are those that emerged from the scan and satisfy the observational constraints. No parameters were adjusted specifically to fit the data after seeing the results; the scan was done independently of the comparison. In the revised version, we will modify the abstract and the description of the numerical study to explicitly state the procedure: first, parameters are fixed to recover the Starobinsky limit, then the space is scanned, and finally the results are compared to data. We will change the wording from 'favored by' to 'consistent with' or 'compatible with' to avoid any implication of post-hoc selection. revision: yes

Circularity Check

0 steps flagged

No significant circularity; standard parameter scan against external data

full rationale

The paper derives standard slow-roll expressions for three parameterized generalizations of the Starobinsky model, numerically evaluates the resulting (n_s, A_s) and (r, n_s) values over the scanned parameter space, and compares those values to the external Planck 2018 contours. The statement that certain parameter choices are 'favored' is a direct outcome of this external comparison rather than any self-referential reduction. No equation is defined in terms of its own output, no fitted parameter is relabeled as a prediction, and no load-bearing step collapses to a self-citation or ansatz imported from the authors' prior work. The reproduction of the Starobinsky limit by appropriate parameter choice is an explicit modeling choice, not a hidden circularity.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the slow-roll approximation being sufficient, the validity of the chosen parameter ranges, and the accuracy of the numerical power-spectrum evaluation. No new entities are introduced; the alpha and power-law indices function as free parameters tuned to data.

free parameters (2)

alpha
Dimensionless parameter controlling the shape of the inflationary potential; values are scanned to match observations.
power-law index
Exponent introduced in the power-law and combined models; chosen to reproduce the Starobinsky limit and fit data.

axioms (1)

domain assumption Slow-roll approximation holds throughout the relevant field range
Used to derive the expressions for the power spectra and observables.

pith-pipeline@v0.9.0 · 5466 in / 1426 out tokens · 31269 ms · 2026-05-12T04:42:29.066864+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

[1]

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

work page 1981
[2]

A. D. Linde. A new inflationary universe scenario: a possible solution of the horizon, flat- ness, homogeneity, isotropy and primordial monopole problems.Phys. Lett. B, 108:389, 1982

work page 1982
[3]

A. H. Guth and S. Pi. Quantum mechanics of the scalar field in the new inflationary universe.Phys. Rev. D, 32:1899, 1985

work page 1985
[4]

Martin, C

J. Martin, C. Ringeval, and V. Vennin. Encyclopaedia Inflationaris.Phys. Dark Univ., 5–6:75, 2014

work page 2014
[5]

V. F. Mukhanov and G.V. Chibisov. Quantum fluctuations and a nonsingular universe. Pis’ma Zh. Eksp. Teor. Fiz., 33:549, 1981

work page 1981
[6]

Jarosiket al.Seven-year wilkinson microwave anisotropy probe (WMAP*) observa- tions: Sky maps, systematic errors, and basic results.ApJS, 192:14, 2011

N. Jarosiket al.Seven-year wilkinson microwave anisotropy probe (WMAP*) observa- tions: Sky maps, systematic errors, and basic results.ApJS, 192:14, 2011

work page 2011
[7]

Akramiet al.Planck 2018 results

Y. Akramiet al.Planck 2018 results. X. Constraints on inflation.Astron. Astrophys., 641:A10, 2020

work page 2018
[8]

A. A. Starobinsky. A new type of isotropic cosmological models without singularity. Phys. Lett. B, 91:99, 1980

work page 1980
[9]

Saini and A

S. Saini and A. Nautiyal. Power lawα-starobinsky inflation.Phys. Rev. D, 113:043546, 2026

work page 2026
[10]

Bhargabi Saha and Malay K. Nandy. Theα-Attractor E-Model in Warm Inflation: Observational Viability from Planck 2018, 2025

work page 2018
[11]

Saini and A

S. Saini and A. Nautiyal. Observational constraints on power law Starobinsky inflation. Phys. Rev. D, 108:123505, 2023

work page 2023
[12]

S. D. Odintsov and V. K. Oikonomou. GeneralizedR p-attractor Cosmology in the Jordan and Einstein Frames: New Type of Attractors and Revisiting Standard Jordan FrameR p Inflation.Int. J. Mod. Phys. D, 32:2250135, 2023

work page 2023
[13]

Rodrigues-da Silva and L

G. Rodrigues-da Silva and L. G. Medeiros. Second–order corrections to Starobinsky inflation.EPJC, 83:1032, 2023

work page 2023
[14]

S. Meza, D. Altamirano, M. Z. Mughal, and Clara Rojas. Numerical analysis of the generalized Starobinsky inflationary model.Int. J. Mod. Phys. D, 30:2150062, 2021

work page 2021
[15]

D. D. Canko, I. D. Gialamas, and G. P. Kodaxis. A simpleF(R, ϕ) deformation of Starobinsky inflationary model.EPJC, 80:458, 2020

work page 2020
[16]

G. K. Chakravarty and S. Mohanty. Power law Starobinsky model of inflation from no–scale SUGRA.Phys. Lett. B, 746:242, 2015

work page 2015
[17]

Kallosh and A

R. Kallosh and A. Linde. Universality Class in Conformal Inflation.JCAP, 07:002, 2013

work page 2013
[18]

Kallosh, A

R. Kallosh, A. Linde, and D. Roest. Superconformal Inflationaryα-Attractors.JHEP, 11:198, 2013

work page 2013
[19]

Ellis, D

J. Ellis, D. V. Nanopoulos, and K. A. Olive. No-Scale Supergravity Realization of the Starobinsky Model of Inflation.Phys. Rev. Lett., 111:111301, 2013

work page 2013
[20]

D. Roest. Universality classes of inflation.JCAP, 01:007, 2014

work page 2014
[21]

Cecotti and R

S. Cecotti and R. Kallosh. Cosmological Attractor Models and Higher Curvature Su- pergravity.JHEP, 05:114, 2014

work page 2014
[22]

Odintsov, and Vasilis K

Shin’ichi Nojiri, Sergei D. Odintsov, and Vasilis K. Oikonomou. Modified gravity theo- ries on a nutshell: Inflation, bounce and late–time evolution.arXiv: General Relativity and Quantum Cosmology, 2017

work page 2017
[23]

Saini and A

S. Saini and A. Nautiyal. Observational Constraints onα–Starobinsky Inflation.Grav- itation & Cosmology, 31:517, 2025. 22 Jordan Zambrano et al

work page 2025
[24]

Kehagias, A

A. Kehagias, A. M. Dizgah, and A. Riotto. Remarks on the Starobinsky model of inflation and its descendants.Phys. Rev. D, 89:043527, 2014

work page 2014
[25]

Ferrara, R

S. Ferrara, R. Kallosh, A. Linde, and M. Porrati. Minimal supergravity models of inflation.Phys. Rev. D, 88:085038, 2013

work page 2013
[26]

S.D Odintsov and V. K. Oikonomou. Inflationaryα-attractors fromF(R) gravity.Phys. Rev. D, 94:124026, 2016

work page 2016
[27]

Frolovsky and S

D. Frolovsky and S. V. Ketov. One-loop corrections to the E-typeα-attractor models of inflation and primordial black hole production.Phys. Rev. D, 111:083533, 2025

work page 2025
[28]

Ellis, M

J. Ellis, M. A. G. Garcia, K. A. Olive, and S. Verner. Constraints on Attractor Models of Inflation and Reheating from Planck, BICEP/Keck, ACT DR6, and SPT-3G Data. 2510.18656, 2025

work page arXiv 2025
[29]

Alho and C

A. Alho and C. Uggla. Inflationaryα-attractor cosmology: A global dynamical systems perspective.Phys. Rev. D, 95:083517, 2017

work page 2017
[30]

K. I. Maeda. Towards the Einstein-Hilbert action via conformal transformation.Phys. Rev. D, 39:3159, 1989

work page 1989
[31]

Odintsov

Shin’ichi Nojiri and Sergei D. Odintsov. Unified cosmic history in modified gravity: From F(R) theory to Lorentz non–invariant models.Phys. Rep., 505:59, 2011

work page 2011
[32]

Muller, H

V. Muller, H. J. Schmidt, and A. A. Starobinsky. Power-law inflation as an attractor solution for inhomogeneous cosmological models.Class. Quantum Grav., 7:1163, 1990

work page 1990
[33]

Sebastiani, G

L. Sebastiani, G. Cognola, R. Myrzakulov, S. D. Odintsov, and S. Zerbini. Nearly starobinsky inflation from modified gravity.Phys. Rev. D, 89, 2014

work page 2014
[34]

J. O. Gong Y. F. Cai and S. Pi. Inflation beyond T-models and primordial B-modes. Phys. Lett. B, 738:20–24, 2014

work page 2014
[35]

A. R. Liddle, and D. H. Lyth.Cosmological inflation and large–scale structure. Cam- bridge university press, 2000

work page 2000
[36]

Mathematica, version 14.3.0.0

Wolfram Research, Inc. Mathematica, version 14.3.0.0. https://www.wolfram.com/mathematica/, 2025

work page 2025
[37]

V. F. Mukhanov, H. A. Feldman, and R. H. Brandenberger. Theory of Cosmological Perturbations.Phys. Rep., 215:203, 1992

work page 1992
[38]

S. D. Odintsov, and V. K. Oikonomou. A power–law inflation tail for the standardR 2– inflation and the Trans–Planckian censorship conjecture.Phys. Lett. B, 865:139458, 2025

work page 2025
[39]

S. D. Odintsov, V. K. Oikonomou, E. I. Manouri, and A. T. Papadopoulos. Solving the trans–Planckian censorship problem with a power–law tail inR 2 inflation: A dynamical system approach.Nucl. Phys. B, 1020:117114, 2025

work page 2025
[40]

Martin and R

J. Martin and R. H. Brandenberger. Trans–Planckian issues for inflationary cosmology. Class. Quantum Grav., 30:113001, 2013

work page 2013
[41]

Habib, A

S. Habib, A. Heinen, K. Heitmann and G. Jungman. Inflationary Perturbations and Precision Cosmology.Phys. Rev. D, 71:043518, 2005

work page 2005
[42]

Zambrano, M

J. Zambrano, M. Agama, M. Garz´ on, W. Br¨ amer–Escamilla, C. Rojas, and Te´ ofilo Vargas. Quintessential inflation studied through semiclassical methods.Int. J. Mod. Phys. D, 35:2550993, 2026

work page 2026
[43]

Giar` e, S

W. Giar` e, S. Pan, E. Di Valentino, W. Yang, J. de Haro, and A. Melchiorri. Inflation- ary Potential as seen from Different Angles: Model Compatibility from Multiple CMB Missions.JCAP, 09:019, 2023

work page 2023
[44]

Das and R

S. Das and R. O. Ramos. Running and Running of the Running of the Scalar Spectral Index in Warm Inflation.Universe, 9:76, 2023

work page 2023
[45]

J. A. Vazquez, L. E. Padilla, and T. Matos. Inflationary cosmology: from theory to observations.Rev. Mex. Fis. E, 17:73, 2020

work page 2020
[46]

Finelliet al

F. Finelliet al. Exploring cosmic origins with CORE: Inflation.JCAP, 04:016, 2018

work page 2018
[47]

J. A. Vazquez, M. Bridges, Y–Z. Ma, and M. P. Hobson. Constraints on the tensor–to– scalar ratio for non–power–law models.JCAP, 08:001, 2013

work page 2013
[48]

ACT DR6 Insights on the Inflationary Attractor models and Reheating, 2025

Md Riajul Haque, Sourav Pal, and Debarun Paul. ACT DR6 Insights on the Inflationary Attractor models and Reheating, 2025

work page 2025
[49]

Burikham, T

P. Burikham, T. Chantavat, and P. Boonaom. Observational constraints on extended starobinsky and weyl gravity model of inflation.J. High Energy Astrophys., 43:178, 2024

work page 2024
[50]

Codello, J

A. Codello, J. Joergensen, F. Sannino, and O. Svendsen. Marginally deformed starobin- sky gravity.JHEP, 2015:050, 2015. Numerical Study of Some Generalizations of the Starobinsky Inflationary Model 23

work page 2015
[51]

Ben-Dayan, S

I. Ben-Dayan, S. Jing, M. Torabian, A. Westphal, and L. Zarate.R 2 logRquantum corrections and the inflationary observables.JCAP, 2014:005, 2014

work page 2014
[52]

J. G. Rodr´ ıguez, S. Santos da Costa, and J. S. Alcaniz. Observational constraints on α–attractor inflationary models with a Higgs–like potential.Phys. Lett. B, 815:136156, 2021

work page 2021
[53]

Adameet al.DESI 2024 III: baryon acoustic oscillations from galaxies and quasars

A.G. Adameet al.DESI 2024 III: baryon acoustic oscillations from galaxies and quasars. JCAP, 04:012, 2025

work page 2024
[54]

V. K. Oikonomou. Strong gravity effects onR 2–corrected single scalar field inflation and compatibility with the ACT data.Phys. Lett. B, 871:139972, 2025

work page 2025
[55]

S. D. Odintsov, and V. K. Oikonomou. Power–lawF(R) gravity as deformations to Starobinsky inflation in view of ACT.Phys. Lett. B, 870:139907, 2025

work page 2025