Cosmology

The epoch nobody saw moves the tilt

Between the end of inflation and the hot universe that made the light elements lies an interval nothing has observed, in which the energy of the inflaton became radiation. How long that took changes how many e-folds before the end the observed scales left — by as many as fourteen — and that moves every model's predicted spectral index by more than the measurement's uncertainty. A potential is never tested by the tilt alone; a potential and a reheating history are tested together.

Assumes Inflation and Inflation.

Every prediction of inflation that has been compared with the sky is evaluated at one moment: when the pivot scale, a comoving wavelength of about 125 megaparsecs, was stretched past the Hubble radius. How far before the end of inflation that moment came is a number, N*, and the tilt reads as a count of it once a potential is granted — 56.5 for a mass term, 54.9 for the Starobinsky plateau.

The count is not a property of the potential. It depends on how much the universe expanded after inflation ended, and part of that expansion happened during an interval nobody has seen: the time it took the energy stored in the inflaton to become the hot plasma of ordinary particles. The essay on the mechanism described that step, reheating, as the least constrained part of the whole account, with a final temperature free across more than fifteen orders of magnitude. That is true. What this essay is about is that the freedom is not invisible. It leaves a mark on the one number the sky measures best.

Starobinsky: e-folds before the end, against how reheating went. N, the number of e-folds between the pivot scale leaving the Hubble radius and the end of inflation, for the Starobinsky potential, against the temperature at which reheating finished, for 3 equations of state during it. All the lines meet on the right at instant reheating, 2.6 × 10¹⁵ GeV, where N = 55.6. The left edge is 5 MeV, below which nucleosynthesis would not have happened. w = 0, oscillating field: 42.0 at 5 MeV; w = ⅓, like radiation: 55.6 at 5 MeV; w = 1, kination: 69.0 at 5 MeV. The reason is how far the universe stretches while the energy density falls: an oscillating field dilutes like matter, as a⁻³, so for the same fall in density it expands further than radiation would, more of the growth of today's scales happens after inflation, and fewer e-folds of inflation are needed to put them where they are. A stiff epoch with w = 1 dilutes as a⁻⁶, stretches less, and needs more. A radiation-like epoch changes nothing. None of this epoch has been observed; the lines are the arithmetic of energy and entropy, and the spread between them is how much an unobserved history moves a quantity the spectral index depends on.
Fig. 1 N*, the e-folds between the pivot scale’s exit and the end of inflation, for the Starobinsky potential, against the temperature at which reheating finished, for three equations of state during it. All three meet on the right at instant reheating, 2.6 × 10¹⁵ GeV, where N* = 55.6. At the left edge, 5 MeV, the lowest temperature nucleosynthesis allows, an oscillating field with w = 0 gives 42.0, a radiation-like w = ⅓ still gives 55.6, and a stiff w = 1 gives 69.0.

Why a count depends on what came after

The logic is geometric and short. A comoving scale leaves the Hubble radius when its wavenumber equals aHaH. The pivot scale’s size today is known — it is 222 times smaller than the present Hubble radius — and so is the expansion rate during inflation, to within the factor the tensor ratio leaves open. What connects the two is the total factor by which the universe has expanded since that scale left, and that factor is a sum of pieces laid end to end: N* e-folds of inflation after the exit, whatever expansion happened during reheating, and the expansion of the radiation and matter eras since.

The last of those is not free. Once the plasma is in thermal equilibrium its entropy is conserved, so its temperature falls in inverse proportion to the scale factor, and the expansion from reheating to today is fixed by how hot reheating left the universe compared with the 2.7 kelvin of the microwave background that is its cooled remnant. The expansion during reheating is fixed by how the energy density fell while it lasted, which depends on the equation of state.

The numbers are worth having in hand. With instant reheating, the universe has expanded by about 1027.610^{27.6} — 63.5 e-folds — since inflation ended: the ratio of the reheating temperature to today’s, corrected for the particle species that have since annihilated and warmed the photons. Add the 55.6 e-folds of inflation that followed the pivot scale’s exit and the scale has grown by about 119 e-folds, a factor of 105210^{52}, from a size far smaller than a proton to 125 megaparsecs. The split between the two halves of that total is what the rest of this essay is about, because the total is fixed and the halves are not.

Anything that makes the universe expand more after inflation leaves less for inflation to do, and N* falls.

An oscillating field dilutes like matter, its energy density falling as a3a^{-3}, where radiation falls as a4a^{-4}. For the same fall in density from the end of inflation to the reheating temperature, matter-like dilution requires a larger expansion than radiation would. A long, matter-like reheating therefore stretches the universe further than an instant one, today’s scales grew more after inflation, and fewer e-folds of inflation were needed to put them where they are. A stiff fluid, w=1w = 1, dilutes as a6a^{-6}, stretches less, and needs more. A radiation-like epoch is indistinguishable from none, which is why its line in the opening figure is flat.

The ledger

Written out, N* is a sum of five terms, and it is worth seeing their sizes, because the familiar “fifty to sixty” hides how it is made.

N* = 55.6, as a difference of two much larger numbers. The number of e-folds between the pivot scale leaving the Hubble radius and the end of inflation, for the Starobinsky potential with instant reheating, built up term by term. The expansion rate during inflation over the rate today contributes 126.7; the pivot scale's size against today's Hubble radius takes away 5.4; the fall in temperature from reheating to the microwave background takes away 64.6; the particle species that annihilated since take away 1.1; and reheating itself takes away 0.0. The total is 55.6. Every term but the last is either measured or fixed by the Standard Model, and the last is the one nothing constrains — it is zero here only if reheating was instantaneous or radiation-like. That a number quoted as "fifty to sixty" is the small remainder of 127 minus about 71 is why so modest a spread in it carries so much of the argument.
Fig. 2 N* for the Starobinsky potential with instant reheating, built up term by term. The expansion rate during inflation over the rate today, ln(H*/H₀), contributes 126.7; the pivot scale’s size against today’s Hubble radius takes away 5.4; the fall in temperature from reheating to the microwave background takes away 64.6; the species that annihilated since take away 1.1; and reheating itself takes away nothing, because it was instantaneous. The total is 55.6.

The largest term is the logarithm of a ratio of expansion rates: HH during inflation was about 105510^{55} times its present value. The second largest, with the opposite sign, is the logarithm of a ratio of temperatures — from the 2.6×10152.6\times10^{15} GeV at which instant reheating would leave the plasma to the 2.3×10132.3\times10^{-13} GeV of the microwave background today. The number of e-folds is what remains when one is subtracted from the other. It is a small difference of two large logarithms, and that is why an unobserved term of only a few units in the ledger moves the answer by an amount that matters.

The ledger also shows which uncertainties do not matter. The expansion rate during inflation enters through a logarithm, so the whole range of energy scales left open by the tensor bound — a factor of ten in rr, a factor of three in HH — shifts N* by little more than one. The disagreement over the present expansion rate, 67 against 73 kilometres a second per megaparsec, shifts it by a tenth. The species term assumes the particle content of the Standard Model when reheating ended; doubling the number of particle species, as a supersymmetric spectrum roughly would, shifts it by a quarter. The reheating term alone is unconstrained, and it can be tens.

N* = 50.7, as a difference of two much larger numbers. The number of e-folds between the pivot scale leaving the Hubble radius and the end of inflation, for the Starobinsky potential with reheating ending at 1.0 × 10⁹ GeV and w = 0.00, built up term by term. The expansion rate during inflation over the rate today contributes 126.7; the pivot scale's size against today's Hubble radius takes away 5.4; the fall in temperature from reheating to the microwave background takes away 49.8; the particle species that annihilated since take away 1.1; and reheating itself takes away 19.8. The total is 50.7. Every term but the last is either measured or fixed by the Standard Model, and the last is the one nothing constrains — it is zero here only if reheating was instantaneous or radiation-like. That a number quoted as "fifty to sixty" is the small remainder of 127 minus about 76 is why so modest a spread in it carries so much of the argument.
Fig. 3 The same ledger with a matter-like reheating that ends at 10⁹ GeV. The expansion-rate term is unchanged at 126.7 and the pivot term at 5.4; the cooling term shrinks to 49.8, because the plasma started cooler; and reheating itself now takes away 19.8, the extra expansion of a matter-like epoch lasting from the end of inflation to that temperature. The net change is five e-folds, and N* falls to 50.7.

The two bars that change in that figure change in opposite directions, and the net effect is the difference between how a matter-like fluid and radiation dilute. That is why the shift is so much smaller than either bar: fifteen of the twenty e-folds of extra expansion during reheating are simply expansion that radiation would have done anyway, over the same fall in temperature. Only the excess counts.

Which way the tilt moves

For a given potential, fewer e-folds means the observed scales left closer to the end of inflation, where the slow-roll parameters are larger. Every potential considered in these essays has a spectral index that falls as the end approaches, so a smaller N* makes the prediction redder. For the plateau the index is close to 12/N1 - 2/N, so a fall of N* from 55.6 to 42 takes it from about 0.964 to about 0.952 — a shift of three times the present uncertainty, from a history that emitted nothing that survives.

Starobinsky: the predicted tilt against the reheating temperature, for w = 0. The spectral index the Starobinsky potential predicts, against the temperature at which reheating ended, for an equation of state w = 0.00 while it lasted. Instant reheating, on the right at 2.6 × 10¹⁵ GeV, gives 0.9653; a reheating that dragged on to 5 MeV gives 0.9545. The shaded band is today's measurement, 0.9649 ± 0.0042, with its two-sigma extent lighter; going down in temperature from instant reheating, the prediction crosses the lower one-sigma edge at about 4.3 × 10⁶ GeV and the lower two-sigma edge at about 1.9 GeV. The dashed lines are a hypothetical measurement with an uncertainty of 0.002, the level planned for the next generation of polarisation experiments, drawn around today's central value — which is an assumption about where the value will land and is labelled as one. At that precision the index would be a thermometer for an epoch that emitted nothing that survives: not a good one, since the temperature enters through a logarithm, but a thermometer where there is currently none.
Fig. 4 The spectral index the Starobinsky potential predicts, against the temperature at which a matter-like reheating ended. Instant reheating gives 0.9653, almost exactly the measured 0.9649; a reheating that dragged on to 5 MeV gives 0.9545. Going down in temperature the prediction leaves the one-sigma band at about 4.3 × 10⁶ GeV and the two-sigma band at about 1.9 GeV. The dashed lines are a hypothetical future uncertainty of ±0.002 drawn around today’s central value.

That is the sense in which the tilt is a thermometer, and also the sense in which it is a poor one. The temperature enters through a logarithm, so the prediction changes by a thousandth for every factor of about thirty in temperature, and the present measurement says only that a matter-like reheating in this model probably finished above a few GeV — which is not a demanding requirement, since the light elements need it finished above a few MeV anyway. The forecast lines show why the argument is worth making now. At an uncertainty of 0.002, the same curve would cross the edge of the band many decades higher in temperature, and a statement about an epoch that emitted nothing that survives would become a limit with teeth.

A stiff reheating runs the other way, and it shows that the direction is not a convention of the calculation.

Starobinsky: the predicted tilt against the reheating temperature, for w = 1. The spectral index the Starobinsky potential predicts, against the temperature at which reheating ended, for an equation of state w = 1.00 while it lasted. Instant reheating, on the right at 2.6 × 10¹⁵ GeV, gives 0.9653; a reheating that dragged on to 5 MeV gives 0.9719. The shaded band is today's measurement, 0.9649 ± 0.0042, with its two-sigma extent lighter; going down in temperature from instant reheating, the prediction crosses the upper one-sigma edge at about 1.7 × 10⁶ GeV. The dashed lines are a hypothetical measurement with an uncertainty of 0.002, the level planned for the next generation of polarisation experiments, drawn around today's central value — which is an assumption about where the value will land and is labelled as one. At that precision the index would be a thermometer for an epoch that emitted nothing that survives: not a good one, since the temperature enters through a logarithm, but a thermometer where there is currently none.
Fig. 5 The same potential with a reheating whose equation of state is w = 1, the stiff fluid of a field whose energy is almost all kinetic. The prediction now moves blue as the reheating temperature falls: from 0.9653 at instant reheating to 0.9719 at 5 MeV, crossing the upper one-sigma edge at about 1.7 × 10⁶ GeV. A kinetic-dominated epoch stretches the universe less than radiation, so more e-folds of inflation are needed and the observed scales left further from the end.

The minimum decides the direction

Whether reheating was matter-like, radiation-like or stiff sounds like one more free parameter, and for most of its duration it is not free at all. It is fixed by the shape of the potential at its minimum, by an argument that is older than cosmology.

After inflation ends the field oscillates about the bottom of its potential, and it oscillates quickly: once the slow-roll conditions fail, the curvature of the minimum exceeds the expansion rate, so the field completes many swings in a single expansion time. The equation of state that governs the expansion is then the average over a swing, set by how the kinetic and potential energies share out over a cycle. For a potential that rises as a power of the field near its minimum, VϕnV \propto \phi^n, the virial theorem fixes that share: the time-averaged kinetic energy is n/2n/2 times the time-averaged potential energy. The pressure is the kinetic energy minus the potential, the density is their sum, and so

w=n2n+2.\langle w \rangle = \frac{n-2}{n+2}.

A quadratic minimum, n=2n = 2, gives w=0w = 0: the oscillating field has no average pressure and dilutes like matter, exactly as a pendulum’s kinetic and potential energies average to equal shares. A quartic minimum, n=4n = 4, gives w=1/3w = 1/3, the equation of state of radiation, with no particles involved at all. A minimum rising as the sixth power gives w=1/2w = 1/2, stiffer than radiation, and its delayed reheating would push the tilt blue. The whole family is one line of the same theorem that weighs a cluster of galaxies from the speeds of its members.

The Starobinsky potential, whatever it looks like on its plateau, is quadratic at its minimum. So once its inflation ends and before its energy has been transferred, it dilutes like matter, and the direction in which its prediction moves is not a choice: a delayed reheating can only redden it. A potential with a quartic minimum would behave like radiation from the start, and its prediction would hardly depend on how long reheating took. Dark energy’s equation of state is the same number measured for a field today; here it is fixed by the curvature of the bottom of a potential nobody has measured.

What the oscillation does not fix is when it stops. The field decays into particles at a rate set by its couplings to them, and the reheating temperature is set by that rate. A feeble coupling, of gravitational strength, can prolong the matter-like phase until the universe has cooled to a few GeV or less. A strong one can end it within a few swings, and not by ordinary decay. A field oscillating in time makes the effective mass of anything coupled to it oscillate too, and a mode whose mass is modulated at twice its own natural frequency grows exponentially, the way a swing is pumped by a child rising and crouching twice per cycle. That parametric resonance, worked out for inflation by Lev Kofman, Andrei Linde and Alexei Starobinsky in 1994, can transfer most of the field’s energy into particles long before any perturbative decay would have started — at which point the expansion stops being matter-like, but the resulting plasma is far from thermal and takes longer still to become the equilibrium the ledger assumes.

A potential and a history, tested together

The consequence for comparing models with the sky is that a model is not a point, or even the short track in e-folds drawn in the plane of tilt and tensor ratio. It is a line parameterised by its reheating temperature.

How far each prediction moves with the reheating temperature. Predictions in the plane of spectral index and tensor-to-scalar ratio, each drawn as a track in the temperature at which reheating finishes rather than as a point, for an equation of state w = 0.00 during it. One end of each track is instant reheating and the other is 5 MeV, the lowest temperature at which nucleosynthesis still works. The vertical band is the measured index, 0.9649 ± 0.0042 with its two-sigma extent, and the shaded region above 0.036 is excluded at 95 per cent by the B-mode polarisation limit — drawn as two independent limits, which the published constraint is not quite: the real likelihood is a correlated contour, and it is somewhat tighter than this box in the corner where the tilt is high. Starobinsky: nₛ 0.9653 and r 3.4e-3 with instant reheating at 2.6 × 10¹⁵ GeV, nₛ 0.9545 and r 5.8e-3 if reheating drags on to 5 MeV (N from 55.6 to 42.0); φ^⅔: nₛ 0.9764 and r 4.7e-2 with instant reheating at 4.2 × 10¹⁵ GeV, nₛ 0.9689 and r 6.2e-2 if reheating drags on to 5 MeV (N from 56.4 to 42.7); natural, f = 7: nₛ 0.9614 and r 7.3e-2 with instant reheating at 2.7 × 10¹⁵ GeV, nₛ 0.9515 and r 1.1e-1 if reheating drags on to 5 MeV (N* from 57.1 to 43.5).
Fig. 6 Three potentials drawn as tracks in their reheating temperature, for a matter-like reheating, from instant reheating to 5 MeV. The plateau moves from ns=0.9653n_s = 0.9653, r = 0.0034 to ns=0.9545n_s = 0.9545, r = 0.0058 as N* falls from 55.6 to 42.0; φ^⅔ from 0.9764 and 0.047 to 0.9689 and 0.062; natural inflation with f = 7 from 0.9614 and 0.073 to 0.9515 and 0.11. Every track runs down and to the left: fewer e-folds, a redder tilt and a larger tensor ratio.

Two things in that figure are not obvious from any single number. The first is that a long reheating raises the tensor ratio as well as reddening the tilt, because the scales left where the slope was steeper. A model that sits just below the tensor bound with instant reheating can cross it if its reheating was slow, and a model rescued in the tilt by a slow reheating pays for the rescue in tensors.

The second is that the tracks are long compared with the band. The plateau’s track spans almost three standard deviations of the tilt. Quoting a potential as “consistent with the data” without saying how it reheated is therefore a statement about a segment of that track, usually the instant end, and the consistency is a joint property of the potential and an assumption.

Natural inflation shows the ordinary case. With instant reheating its index, 0.9614, is already near the lower edge of the band, and a matter-like reheating carries it out of one sigma below about 7×10137\times10^{13} GeV and out of two sigma below about 2×1052\times10^{5} GeV. For that potential the tilt prefers a quick reheating, and its tensor ratio of 0.073 is above the bound whatever the reheating did, so the tilt’s preference is academic. The power law with an exponent of two-thirds shows the trade in its sharpest form.

φ^⅔: the predicted tilt against the reheating temperature, for w = 0. The spectral index the φ^⅔ potential predicts, against the temperature at which reheating ended, for an equation of state w = 0.00 while it lasted. Instant reheating, on the right at 4.2 × 10¹⁵ GeV, gives 0.9764; a reheating that dragged on to 5 MeV gives 0.9689. The shaded band is today's measurement, 0.9649 ± 0.0042, with its two-sigma extent lighter; going down in temperature from instant reheating, the prediction crosses the upper two-sigma edge at about 9.6 × 10⁶ GeV and the upper one-sigma edge at about 13 MeV. The dashed lines are a hypothetical measurement with an uncertainty of 0.002, the level planned for the next generation of polarisation experiments, drawn around today's central value — which is an assumption about where the value will land and is labelled as one. At that precision the index would be a thermometer for an epoch that emitted nothing that survives: not a good one, since the temperature enters through a logarithm, but a thermometer where there is currently none.
Fig. 7 The spectral index of the φ^⅔ potential against the temperature at which a matter-like reheating ended. With instant reheating it is 0.9764, nearly three standard deviations too blue. Delaying reheating reddens it, and going down in temperature it crosses the upper two-sigma edge at about 9.6 × 10⁶ GeV and the upper one-sigma edge only at about 13 MeV, barely above the floor nucleosynthesis sets, where it reaches 0.9689.

In the tilt alone, that potential is compatible with the data only if the universe took almost the longest time it could to reheat — a reheating that ended a few MeV above the point at which the helium abundance and the count of light particle species would have gone wrong. At that end of its track, though, its tensor ratio is 0.062, well above the bound. The slow reheating that rescues its tilt condemns it in tensors, and the plane shows both at once. The model is excluded not by either number but by the fact that no single reheating history satisfies both.

What is actually measured

Nothing about reheating itself has been observed, and nothing in this essay changes that. The inference runs entirely through arithmetic that is well tested at the other end: entropy conservation in an expanding plasma, which is the same physics that makes four light-element abundances fit one free parameter, and the relation between an equation of state and the rate at which density falls.

The two ends of the temperature axis are measured, and they are measured quite differently. The lower end, a few MeV, comes from the requirement that neutrinos had time to reach thermal equilibrium before nucleosynthesis; a later reheating would leave too few of them, which would change both the expansion rate during element formation and the damping of the microwave background’s smallest scales, and neither change is seen. The upper end is the energy density at the end of inflation, which for the plateau is fixed by the measured amplitude of the fluctuations once the potential is chosen.

Between those ends the assumption that matters most is that entropy was conserved after reheating. Any later injection of entropy — a heavy particle that decayed after the plasma formed, a phase transition that released latent heat — stretches the universe in the same way a slow reheating does and would be read, in this ledger, as one. The microwave background’s spectrum limits energy injected later than a few months after the beginning, because a blackbody that perfect could not survive it; it says nothing about injections in the first second.

What the ledger leaves out

The equation of state is held constant. A real reheating passes through stages: coherent oscillation, a resonant burst of particle production, a period of turbulence, and thermalisation, and its average ww drifts from zero towards one third as it goes. The figures bracket the possibilities with constant values rather than following any one history.

Instant reheating is a limit, not a scenario. The highest temperature drawn assumes all of the energy at the end of inflation became radiation at once, which no known coupling achieves exactly; the real curve starts slightly to the left of the point where the lines meet.

The species count is the Standard Model’s. If there are particles beyond it, the entropy they carried into the plasma changes the cooling term by a third of the logarithm of the extra count — a fraction of an e-fold for any plausible spectrum.

And some theories cap the temperature from above. In supersymmetric theories a reheating temperature much above 10910^9 GeV overproduces gravitinos, whose late decays would disrupt the light elements. That bound depends on a theory with no experimental support, and it is not drawn.

Still open: what the inflaton decayed into

The ledger needs only a temperature and an equation of state, and the tilt can in principle constrain both. What it cannot reach is the physics that set them: the couplings through which the inflaton gave its energy to the Standard Model, which decide whether the transfer was a resonant burst or a slow leak and whether it produced the particles that survive today. A second field that decayed later, like the curvaton that would leave the sky lopsided, would add its own chapter to the same ledger. How the potential energy of a field became every particle now in existence is still a question with a temperature attached and no mechanism.

About the same objects

Not linked from either essay — found by the objects both name.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

Big bang nucleosynthesisE-foldsEquation of stateHubble parameterInflatonReheatingScalar spectral indexSlow-rollTensor to scalar ratioVirial theorem