Slow-Roll Inflation

The slow-roll approximation is the backbone of inflationary cosmology — it reduces the full nonlinear inflaton dynamics to a tractable system, enabling analytic predictions for the duration of inflation, the spectrum of perturbations, and the connection between microphysics and observation.

1. Introduction — Why Slow-Roll?

The simplest physical picture of inflation is a ball rolling slowly down a gentle hillside. The "ball" is the inflaton field $\phi$, and the "hill" is the potential energy function $V(\phi)$. While the ball creeps along the nearly flat portion of the potential, the vacuum energy $V(\phi)$ dominates the energy budget of the universe. Because $V$ is approximately constant, the Friedmann equation gives a nearly constant Hubble parameter:

$H^2 \approx \frac{V(\phi)}{3M_P^2}$

and the scale factor grows quasi-exponentially: $a(t) \propto e^{Ht}$. The essential requirement is that the inflaton's kinetic energy remain negligible compared to its potential energy throughout inflation. This demands two conditions:

  • Flatness of the potential: The slope $V'(\phi)$ must be small enough that the field does not pick up too much kinetic energy.
  • Persistence: The curvature $V''(\phi)$ must also be small so that the flatness condition is not rapidly violated — inflation must last long enough to solve the horizon and flatness problems (at least $N \sim 50{-}60$ e-folds).

The slow-roll approximation formalizes these requirements via dimensionless parameters. When they are much less than unity, the second-order Klein-Gordon equation becomes first-order and all observables are expressed directly in terms of the potential and its derivatives.

Physical Intuition

Think of a marble on a slightly tilted table coated with honey. The Hubble drag $3H\dot{\phi}$ prevents acceleration; the marble moves at a terminal velocity balancing the slope $V'$ against friction $3H$. This is the slow-roll regime: $3H\dot{\phi} \approx -V'(\phi)$.

2. Derivation: Potential Slow-Roll Parameters

We start from the full equations governing an FRW universe with a single minimally-coupled scalar field:

$H^2 = \frac{1}{3M_P^2}\left(\frac{1}{2}\dot{\phi}^2 + V(\phi)\right)$

$\ddot{\phi} + 3H\dot{\phi} + V'(\phi) = 0$

The slow-roll approximation makes two assumptions:

$\text{(SR1):} \quad \frac{1}{2}\dot{\phi}^2 \ll V(\phi)$

$\text{(SR2):} \quad |\ddot{\phi}| \ll |3H\dot{\phi}|$

Under SR1, the Friedmann equation simplifies to $H^2 \approx V/(3M_P^2)$. Under SR2, the Klein-Gordon equation becomes:

$3H\dot{\phi} \approx -V'(\phi)$

We can now solve for $\dot{\phi}$:

$\dot{\phi} \approx -\frac{V'}{3H} \approx -\frac{V' M_P}{\sqrt{3V}}$

Defining $\epsilon_V$

Condition SR1 requires $\frac{1}{2}\dot{\phi}^2 \ll V$. Substituting the slow-roll expression for $\dot{\phi}$:

$\frac{1}{2}\dot{\phi}^2 \approx \frac{V'^2}{18H^2} \approx \frac{V'^2 M_P^2}{6V}$

The condition $\frac{1}{2}\dot{\phi}^2 \ll V$ then becomes:

$\frac{M_P^2}{2}\left(\frac{V'}{V}\right)^2 \ll 1$

This motivates the definition of the first potential slow-roll parameter:

$\epsilon_V \equiv \frac{M_P^2}{2}\left(\frac{V'}{V}\right)^2 \ll 1$

When $\epsilon_V \ll 1$, the kinetic energy is negligible compared to the potential energy, and inflation proceeds. Moreover, inflation ends when $\epsilon_V = 1$ — at this point the kinetic and potential energies become comparable, and the accelerated expansion ceases. More precisely, the condition for accelerated expansion ($\ddot{a} > 0$) is equivalent to $\epsilon_H < 1$, and in the slow-roll regime $\epsilon_H \approx \epsilon_V$.

Defining $\eta_V$

Condition SR2 requires $|\ddot{\phi}| \ll |3H\dot{\phi}|$. To see what this implies about the potential, differentiate the slow-roll equation $3H\dot{\phi} \approx -V'$with respect to time:

$3\dot{H}\dot{\phi} + 3H\ddot{\phi} \approx -V''\dot{\phi}$

Dropping the $\ddot{\phi}$ term (by assumption SR2) and the $\dot{H}$ term (which is of the same order as $\epsilon_V$ corrections), we find:

$\frac{\ddot{\phi}}{H\dot{\phi}} \approx -\frac{V''}{3H^2} + \frac{\dot{H}}{H^2} \approx -\frac{M_P^2 V''}{V} + \epsilon_V$

Thus $|\ddot{\phi}| \ll |H\dot{\phi}|$ requires both $\epsilon_V \ll 1$ and:

$\eta_V \equiv M_P^2 \frac{V''}{V} \quad \text{with} \quad |\eta_V| \ll 1$

The parameter $\eta_V$ measures the curvature of the potential relative to its height. While $\epsilon_V \ll 1$ ensures that inflation is occurring at a given moment,$|\eta_V| \ll 1$ ensures that the flatness condition persists — that is, inflation lasts long enough. If $|\eta_V|$ were of order unity, the potential would curve sharply and the field would quickly roll into a steep region, ending inflation prematurely.

Summary of Conditions

$\epsilon_V \ll 1$ — the potential is flat enough for inflation to occur (kinetic energy subdominant).

$|\eta_V| \ll 1$ — the potential stays flat long enough for inflation to last $N \gtrsim 50$ e-folds.

Inflation ends when $\epsilon_V(\phi_\text{end}) = 1$ (or equivalently $\epsilon_H = 1$).

3. Derivation: Hubble Slow-Roll Parameters

An alternative and in many ways more fundamental set of slow-roll parameters is defined directly in terms of the Hubble parameter and its time derivatives, without reference to a specific potential. These are called the Hubble slow-roll parameters (or Hubble flow parameters).

Defining $\epsilon_H$

The condition for inflation ($\ddot{a} > 0$) can be rewritten. From$a(t) = a_0 e^{\int H\, dt}$, we compute:

$\frac{\ddot{a}}{a} = \dot{H} + H^2$

So $\ddot{a} > 0$ requires $\dot{H} + H^2 > 0$, i.e.,$-\dot{H}/H^2 < 1$. This motivates:

$\epsilon_H \equiv -\frac{\dot{H}}{H^2}$

Inflation occurs when $\epsilon_H < 1$, and slow-roll inflation requires$\epsilon_H \ll 1$. Let us now express $\epsilon_H$ in terms of the inflaton field. From the Friedmann and Raychaudhuri equations:

$H^2 = \frac{1}{3M_P^2}\left(\frac{1}{2}\dot{\phi}^2 + V\right)$

$\dot{H} = -\frac{\dot{\phi}^2}{2M_P^2}$

The second equation uses $\dot{H} = -(1/2M_P^2)(\rho + p)$ with$\rho + p = \dot{\phi}^2$. Therefore:

$\epsilon_H = \frac{\dot{\phi}^2}{2M_P^2 H^2}$

This is exact — no approximation. It shows $\epsilon_H$ measures the ratio of kinetic to total energy.

Defining $\eta_H$

The second Hubble slow-roll parameter encodes the condition SR2 ($|\ddot{\phi}| \ll |3H\dot{\phi}|$):

$\eta_H \equiv -\frac{\ddot{\phi}}{H\dot{\phi}}$

When $|\eta_H| \ll 1$, the acceleration of the field is negligible compared to the Hubble friction term.

Relating Hubble and Potential Parameters

We now derive the connection between the two sets of parameters. In the slow-roll regime, using $H^2 \approx V/(3M_P^2)$ and $\dot{\phi} \approx -V'/(3H)$:

$\epsilon_H = \frac{\dot{\phi}^2}{2M_P^2 H^2} \approx \frac{V'^2/(9H^2)}{2M_P^2 H^2} = \frac{V'^2}{18 H^4 M_P^2}$

Substituting $H^2 \approx V/(3M_P^2)$ so that $H^4 \approx V^2/(9M_P^4)$:

$\epsilon_H \approx \frac{V'^2 \cdot 9M_P^4}{18 V^2 M_P^2} = \frac{M_P^2}{2}\left(\frac{V'}{V}\right)^2 = \epsilon_V$

$\epsilon_H \approx \epsilon_V \quad \text{(to lowest order in slow-roll)}$

For the second parameter, we differentiate $\dot{\phi} \approx -V'/(3H)$ with respect to time:

$\ddot{\phi} \approx -\frac{V''\dot{\phi}}{3H} + \frac{V'\dot{H}}{3H^2}$

Dividing by $-H\dot{\phi}$:

$\eta_H = -\frac{\ddot{\phi}}{H\dot{\phi}} \approx \frac{V''}{3H^2} - \frac{V'\dot{H}}{3H^3\dot{\phi}}$

Using $3H^2 \approx V/M_P^2$ for the first term and $\dot{H}/H^2 = -\epsilon_H \approx -\epsilon_V$ together with $V'/(3H\dot{\phi}) \approx -1$ for the second:

$\eta_H \approx \frac{M_P^2 V''}{V} - \epsilon_V = \eta_V - \epsilon_V$

$\eta_H \approx \eta_V - \epsilon_V$

The two sets agree at leading order; the distinction matters at second order, relevant for next-to-leading corrections to $n_s$ and its running.

4. Derivation: Number of e-Folds

The number of e-folds measures how much the universe expands during inflation. It is defined as:

$N \equiv \ln\frac{a_\text{end}}{a} = \int_t^{t_\text{end}} H\, dt$

To express this as an integral over the field, we change variable from $t$to $\phi$. Since $dt = d\phi/\dot{\phi}$:

$N = \int_{\phi}^{\phi_\text{end}} \frac{H}{\dot{\phi}}\, d\phi$

Note the sign: since the field rolls from large to small values (for typical potentials),$\dot{\phi} < 0$ and $\phi > \phi_\text{end}$, so we write with the convention that $N > 0$:

$N = \int_{\phi_\text{end}}^{\phi} \frac{H}{|\dot{\phi}|}\, d\phi$

Applying the slow-roll approximation $3H\dot{\phi} \approx -V'$, so$H/|\dot{\phi}| \approx 3H^2/|V'| \approx V/(M_P^2|V'|)$:

$N \approx \frac{1}{M_P^2}\int_{\phi_\text{end}}^{\phi} \frac{V}{V'}\, d\phi$

This is the master formula for the number of e-folds in slow-roll inflation. Given a potential, one computes $\phi_\text{end}$ from $\epsilon_V(\phi_\text{end}) = 1$, and then evaluates the integral to find how many e-folds of inflation occur for a given initial field value.

Example 1: Quadratic Potential $V = \frac{1}{2}m^2\phi^2$

Here $V' = m^2\phi$, so $V/V' = \phi/2$. The slow-roll parameter is:

$\epsilon_V = \frac{M_P^2}{2}\left(\frac{m^2\phi}{\frac{1}{2}m^2\phi^2}\right)^2 = \frac{2M_P^2}{\phi^2}$

Inflation ends when $\epsilon_V = 1$, giving $\phi_\text{end} = \sqrt{2}\, M_P$. The e-fold integral is:

$N = \frac{1}{M_P^2}\int_{\phi_\text{end}}^{\phi} \frac{\phi}{2}\, d\phi = \frac{\phi^2 - \phi_\text{end}^2}{4M_P^2} = \frac{\phi^2}{4M_P^2} - \frac{1}{2}$

For large $\phi$ (i.e., $\phi \gg \phi_\text{end}$):

$N \approx \frac{\phi^2}{4M_P^2} \quad \Longrightarrow \quad \phi_* \approx 2\sqrt{N}\, M_P$

For $N = 60$, the field value at horizon exit is $\phi_* \approx 15.5\, M_P$, deep in the super-Planckian regime. This is the well-known Lyth bound implication for large-field models.

Example 2: Quartic Potential $V = \frac{\lambda}{4}\phi^4$

Here $V' = \lambda\phi^3$, so $V/V' = \phi/4$. The slow-roll parameter:

$\epsilon_V = \frac{M_P^2}{2}\left(\frac{4}{\phi}\right)^2 = \frac{8M_P^2}{\phi^2}$

Inflation ends at $\phi_\text{end} = 2\sqrt{2}\,M_P$. The e-folds:

$N = \frac{1}{M_P^2}\int_{\phi_\text{end}}^{\phi}\frac{\phi}{4}\,d\phi = \frac{\phi^2 - \phi_\text{end}^2}{8M_P^2}$

$N \approx \frac{\phi^2}{8M_P^2} \quad \Longrightarrow \quad \phi_* \approx 2\sqrt{2N}\, M_P$

Example 3: Starobinsky Model $V = \Lambda^4\left(1 - e^{-\sqrt{2/3}\,\phi/M_P}\right)^2$

The Starobinsky potential arises from $R + R^2/(6M^2)$ gravity after a conformal transformation. Define $x = e^{-\sqrt{2/3}\,\phi/M_P}$ so $V = \Lambda^4(1-x)^2$. Then:

$V' = \Lambda^4 \cdot 2(1-x)\cdot\sqrt{\frac{2}{3}}\frac{x}{M_P}, \qquad \frac{V'}{V} = \frac{2x}{(1-x)}\sqrt{\frac{2}{3}}\frac{1}{M_P}$

$\epsilon_V = \frac{4}{3}\frac{x^2}{(1-x)^2}$

For $\phi \gg M_P$ (i.e., $x \ll 1$), $\epsilon_V \approx (4/3)x^2 \ll 1$, which is exponentially small. The e-fold integral gives:

$N \approx \frac{3}{4}e^{\sqrt{2/3}\,\phi/M_P} \approx \frac{3}{4x}$

So $x \approx 3/(4N)$, and $\epsilon_V \approx 3/(4N^2)$,$\eta_V \approx -1/N$. For $N = 55$ this gives $n_s \approx 1 - 2/N \approx 0.964$and $r \approx 12/N^2 \approx 0.004$ — in excellent agreement with Planck data.

5. Derivation: Attractor Behavior & Hamilton-Jacobi Formulation

A remarkable property of slow-roll inflation is that it is an attractor: regardless of initial conditions (provided the field starts in the slow-roll region), different trajectories converge to the same solution. The Hamilton-Jacobi formalism makes this manifest.

The Hamilton-Jacobi Equation

Instead of treating $\phi(t)$ as the dynamical variable, we use the inflaton field value as the time variable and express everything in terms of $H(\phi)$. From $\dot{H} = H'(\phi)\dot{\phi}$ and$\dot{H} = -\dot{\phi}^2/(2M_P^2)$, we obtain:

$\dot{\phi} = -2M_P^2 H'(\phi)$

This is an exact first-order equation for $\dot{\phi}$, replacing the second-order Klein-Gordon equation. Substituting into the Friedmann equation$H^2 = (\dot{\phi}^2/2 + V)/(3M_P^2)$:

$H^2 = \frac{1}{3M_P^2}\left(2M_P^4 H'^2 + V\right)$

Rearranging gives the Hamilton-Jacobi equation:

$\left[H'(\phi)\right]^2 - \frac{3}{2M_P^2}H(\phi)^2 = -\frac{V(\phi)}{2M_P^4}$

Given $V(\phi)$, this is a first-order ODE for $H(\phi)$. Different "initial conditions" correspond to different integration constants, i.e., different amounts of kinetic energy at a given field value.

Proving the Attractor Property

Let $H_0(\phi)$ be a particular solution and consider a perturbation$H(\phi) = H_0(\phi) + \delta H(\phi)$. Substituting into the Hamilton-Jacobi equation and linearizing:

$2H_0' \delta H' - \frac{3}{M_P^2}H_0 \delta H \approx 0$

This gives:

$\frac{\delta H'}{\delta H} = \frac{3H_0}{2M_P^2 H_0'}$

Using $\dot{\phi} = -2M_P^2 H_0'$ and the number of e-folds$dN = -Hd\phi/\dot{\phi} = H_0/(2M_P^2 H_0')d\phi$:

$\delta H \propto \exp\left(-3\int dN\right) = e^{-3N}$

The perturbation $\delta H$ decays as $e^{-3N}$ — exponentially fast in the number of e-folds. After just a few e-folds, all initial conditions converge to the same slow-roll trajectory. This is why the slow-roll attractor is so robust: the universe "forgets" its initial kinetic energy within the first handful of e-folds.

In terms of the slow-roll parameter, the Hubble parameter in the Hamilton-Jacobi formalism gives:

$\epsilon_H = 2M_P^2\left(\frac{H'}{H}\right)^2$

This provides a clean parametrization: the slow-roll hierarchy is simply$\epsilon_H \ll 1$, $\eta_H \ll 1$, and higher-order flow parameters are small. The Hamilton-Jacobi approach forms the basis for the inflationary flow equations used in Monte Carlo reconstructions of the inflaton potential.

6. Applications

Computing Observables for Specific Models

The slow-roll parameters connect directly to CMB observables through the relations:

$n_s = 1 - 6\epsilon_V + 2\eta_V$

$r = 16\epsilon_V$

$n_t = -2\epsilon_V$

$\alpha_s = \frac{dn_s}{d\ln k} = 16\epsilon_V\eta_V - 24\epsilon_V^2 - 2\xi_V^2$

where $\xi_V^2 = M_P^4 V' V'''/ V^2$ is the third slow-roll parameter governing the running of the spectral index.

$V = \frac{1}{2}m^2\phi^2$

$\epsilon_V = 2M_P^2/\phi^2 \approx 1/(2N)$

$\eta_V = 2M_P^2/\phi^2 \approx 1/(2N)$

$n_s \approx 1 - 2/N \approx 0.967$

$r \approx 8/N \approx 0.13$

Disfavored by Planck+BICEP (r too large)

Starobinsky

$\epsilon_V \approx 3/(4N^2)$

$\eta_V \approx -1/N$

$n_s \approx 1 - 2/N \approx 0.964$

$r \approx 12/N^2 \approx 0.004$

Excellent agreement with Planck data

Flow Equations

The inflationary flow equations, introduced by Kinney (2002), extend the slow-roll hierarchy to infinite order. Define parameters $ {}^\ell\lambda_H$ via:

$ {}^\ell\lambda_H = \left(2M_P^2\right)^\ell \frac{(H')^{\ell-1}}{H^\ell}\frac{d^{\ell+1}H}{d\phi^{\ell+1}}$

with $ {}^0\lambda_H = \epsilon_H$ and $ {}^1\lambda_H = \eta_H$. The flow equations form a coupled system of ODEs in the number of e-folds $N$. These can be integrated numerically to explore the space of inflationary models without committing to a specific potential.

Monte Carlo Reconstruction

By randomly sampling initial values of the flow parameters and evolving them, one can generate a large ensemble of inflationary models and map out the allowed region in the $(n_s, r)$ plane. This "Monte Carlo reconstruction" approach (Easther & Kinney, 2003) provides a model-independent way to test the slow-roll paradigm against CMB data. The results confirm that the observed values of$n_s$ and the upper bound on $r$ are naturally accommodated by the simplest slow-roll models, particularly plateau-type potentials like the Starobinsky model.

7. Historical Context

The slow-roll paradigm emerged from early attempts to fix the "graceful exit" problem of Guth's original inflationary model (1981), which relied on a first-order phase transition that produced an unacceptable inhomogeneous universe.

Linde (1982) — New & Chaotic Inflation

Linde proposed "new inflation" (slow roll from an unstable maximum) and later "chaotic inflation" (1983), where the inflaton starts at super-Planckian values with $V = m^2\phi^2/2$. The slow-roll conditions were implicit in these models.

Steinhardt & Turner (1984)

Formalized the slow-roll approximation, identifying conditions for the dynamics to reduce to first order and establishing the connection between slow-roll parameters and e-fold duration.

Lyth & Riotto Review (1999)

Systematized the formalism in Phys. Rep. 314, connecting slow-roll to particle physics model building, deriving perturbation spectra to second order, and providing the dictionary between potentials and observables still used today.

Further Developments

Hamilton-Jacobi formulation: Salopek & Bond (1990), Muslimov (1990). Systematic slow-roll hierarchy: Liddle, Parsons & Barrow (1994). Flow equations: Kinney (2002). Monte Carlo reconstruction: Easther & Kinney (2003). Baumann's TASI lectures (2009) and the Planck inflation papers (2013–2018) cemented slow-roll as the standard framework.

8. Python Simulation: Slow-Roll Analysis

The simulation below computes and compares the slow-roll parameters, e-fold counts, and field trajectories for three inflationary potentials: the quadratic model ($m^2\phi^2/2$), the Starobinsky model, and natural inflation. For each model it plots $\epsilon_V(\phi)$,$\eta_V(\phi)$, $N(\phi)$, and compares the slow-roll approximate trajectory against the exact numerical solution of the full Klein-Gordon + Friedmann system.

Slow-Roll Parameters, e-Folds, and Exact vs Approximate Trajectories

Python
script.py122 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server