Evading the BBN bound with a soft stiff period

Cosmic inflation is the leading theory to explain early Universe history and structure formation. Non-oscillatory inflation is a class of models which can naturally introduce a post-inflationary stiff period of the Universe’s evolution which boosts the signal of primordial gravitational waves (GWs), making it possible to observe them in forthcoming GW experiments. However, this pushes the GW energy density high enough to destabilise the process of Big Bang Nucleosynthesis (BBN). This problem can be overcome by “softening” the stiff period, so that the field is gradually tending towards freefall from a frozen start. Here, we consider a modified hybrid inflation model where the stiff period is driven by the waterfall field, allowing the barotropic parameter of the Universe to vary, so that it does not violate the BBN constraint but produces a characteristic gravitational wave spectrum soon to be observable.

💡 Research Summary

The paper addresses a long‑standing tension in early‑Universe cosmology: a post‑inflationary stiff (or “kination”) phase can dramatically amplify the stochastic gravitational‑wave background (SGWB), potentially bringing it within reach of upcoming detectors such as LISA, DECIGO, or BBO. However, the same amplification often leads to an integrated GW energy density that exceeds the bound imposed by Big‑Bang Nucleosynthesis (BBN), usually expressed as a limit on the effective number of extra neutrino species, ΔN_eff. Traditional approaches either avoid a stiff phase altogether or assume a constant equation‑of‑state w in the range 1/3 < w < 1, which requires a fine‑tuned balance between kinetic and potential energies and still struggles to satisfy the ΔN_eff constraint while producing an observable signal.

To overcome this, the authors propose a “soft” stiff period in which the barotropic parameter w evolves smoothly from –1 (field frozen) to +1 (free‑fall). The concrete realization is a modified hybrid inflation model. The inflaton σ drives standard slow‑roll inflation with a Coleman–Weinberg‑type potential, while a second scalar – the waterfall field ϕ – is initially held at the origin by a quadratic interaction term ΔV = ½ g²σ²ϕ². After σ reaches the end of inflation and oscillates around zero, the interaction shuts off, allowing ϕ to roll down a double‑exponential potential
 V(ϕ) = V₀ exp