A Minimal and Stable Vacuum Bounce in Exponential $f(R)$ Gravity

Abstract

We investigate the realization of a nonsingular cosmological bounce in metric $f(R)$ gravity using a controlled exponential deformation of the Starobinsky $R^{2}$ model. Adopting a smooth Gaussian-type bouncing scale factor, we first demonstrate a no-go result showing that a positive-curvature vacuum bounce cannot be supported by the model $f(R)=R+αR^{2}(1-e^{-R/R_b})$ alone. We then show that a minimal extension obtained by introducing a constant term restores the bounce exactly, with the constant fixed algebraically by the bounce condition. A systematic parameter-space scan is performed to identify regions free of ghost and tachyonic instabilities. Working in the Einstein frame, we study the evolution of scalar and tensor perturbations across the bounce and show that both remain finite and well behaved. Our results establish a minimal, perturbatively stable realization of a vacuum bounce in $f(R)$ gravity that goes beyond background-level constructions.

Figures (6)

• Figure 1: \ref{['fig:1']} The Gaussian-type bouncing scale factor $a(t)=a_b e^{\lambda t^2}$, \ref{['fig:2']} Phase portrait in the $(a,\dot a)$ plane, showing that the trajectory passes smoothly through $\dot a=0$, confirming the nonsingular bounce; \ref{['fig:3']} Phase portrait in the $(H,\dot H)$ plane. The straight-line trajectory reflects the fact that $\dot H=2\lambda$ is constant, corresponding to a maximally regular bounce with positive acceleration at $H=0$.
• Figure 2: Curvature evolution associated with the Gaussian bouncing background $a(t)=a_b e^{\lambda t^2}$. Panel (a) shows the time dependence of the Ricci scalar $R(t)=12\lambda+48\lambda^2 t^2$, which remains finite and positive throughout the evolution and attains its minimum value at the bounce point $t=0$. Panel (b) displays the phase portrait in the $(R,\dot R)$ plane, where $\dot R(t)=96\lambda^2 t$. The trajectory crosses $\dot R=0$ smoothly at finite $R$, indicating that the bounce is free of curvature singularities and that the curvature dynamics are regular across the transition from contraction to expansion.
• Figure 3: Parameter-space scan for the extended model $f(R)=R-2\Lambda+\alpha R^2\!\left(1-e^{-R/R_b}\right)$ in terms of the dimensionless parameters $(\bar{\alpha},\bar{R}_b)$ defined in Eq. (\ref{['eq:dimensionless_defs']}). For each grid point, $\Lambda$ is fixed by the bounce condition via Eq. (\ref{['eq:fR_extended']}), and the ghost- and tachyon-free requirements $f_R(\bar{R}(t))>0$ and $f_{RR}(\bar{R}(t))>0$ are imposed throughout the bounce interval $t\in[-T,T]$ as in Eqs. (\ref{['eq:ghost_condition']})--(\ref{['eq:tachyon_condition']}), with $\bar{R}(t)=1+4\lambda t^2$ Eq. (\ref{['eq:Rbar_t']}). Panel \ref{['fig:4']} displays the resulting viability classification over the scan grid, while Panel \ref{['fig:5']} presents the same outcome in an alternative view for readability.
• Figure 4: Einstein-frame description of the vacuum bounce for the extended model $f(R)=R-2\Lambda+\alpha R^{2}\!\left(1-e^{-R/R_b}\right)$. Panel \ref{['fig:9']} shows the scalaron potential $V(\phi)=\frac{M_{\textrm{P}l}^2}{2}\,(R f_R-f)/f_R^{2}$ as a function of the canonically normalized scalar field $\phi=\sqrt{3/2}\,M_{\textrm{P}l}\ln f_R$. Panel \ref{['fig:10']} displays the time evolution of $\phi(t)$ along the Jordan-frame bouncing background. The scalar field remains finite and evolves smoothly across the bounce, reaching a turning point near $t=0$, which signals the occurrence of the bounce in the Einstein frame.
• Figure 5: Tensor perturbations across the vacuum bounce in the Einstein frame. Panel \ref{['fig:11']} shows the effective friction term $\Gamma(t)=3\tilde{H}(t)$ governing the evolution of tensor modes, computed from the background quantities constructed in Sec. \ref{['sec:perturbations']}. The friction term remains finite and continuous throughout the bounce. Panel \ref{['fig:12']} displays the numerical evolution of tensor perturbation modes $\tilde{h}_k(\tilde{t})$ for several representative comoving wavenumbers $k$. All modes remain finite and well behaved across the bounce, with no evidence of divergent or exponentially growing behavior.