Bayesian analysis of $\alpha$-Starobinsky model with Planck, ACT and DESI data

TL;DR

This paper introduces a Bayesian analysis method for the $oldsymbol{ extalpha}$-Starobinsky inflationary model using multiple cosmological datasets, directly constraining primordial observables to improve accuracy and assess model viability.

Contribution

It develops an alternative sampling approach that directly constrains primordial observables, avoiding slow-roll approximation errors, and applies it to evaluate the $oldsymbol{ extalpha}$-Starobinsky model with current data.

Findings

Pure Starobinsky model requires >60 e-folds of inflation.

Allowing $oldsymbol{ extalpha}$>1 fits data better.

ACT DR6 lensing data has minimal impact on constraints.

Abstract

We present a joint Bayesian analysis to impose constraints on the generalized $\alpha$-Starobinsky inflationary model using the high-precision cosmological datasets: Planck, CMB lensing from ACT DR6, and Baryon Acoustic Oscillations (BAO) from DESI DR2. For the parameter inference, we introduce an alternative sampling approach. Rather than imposing priors on the cosmological parameters of the inflationary potential $(V_0, \, \alpha, \, N_*)$, we place priors directly on the primordial physical observables $(A_s,\, n_s,\, r)$ through using analytical slow-roll consistency relations, our pipeline internally maps these sampled observables to the corresponding $\alpha$-Starobinsky parameters. These values are then passed to a modified version of $\tt{CLASS}$, which solves the exact inflationary dynamics fully numerically. This pipeline ensures that the final reported posteriors for the observables are computed exactly, completely free from the slow-roll approximation. Applying this methodology, we explore the viability of the $\alpha$-Starobinsky model. We show that, when the full combined dataset is considered, the pure Starobinsky model (i.e., the canonical limit $\alpha = 1$) faces an apparent discrepancy: it requires a large number of $e$-folds of inflation after horizon crossing ($N_* > 60$) due to the shift in the scalar spectral index, $n_s$. In contrast, allowing the deformation parameter $\alpha$ as a free parameter yields a clear $1\sigma$ preference for $\log_{10} \alpha > 0$ present across all datasets. By favoring a broader inflationary plateau, the $\alpha$-Starobinsky model elegantly reconciles theoretically sound expansion histories with empirical data. Notably, we also show that the addition of ACT DR6 lensing data introduces no significant impact on these primordial constraints, confirming that our robust posteriors are primarily driven by Planck and DESI measurements.