Pith. Sign in

REVIEW

Towards the Theory of Reheating After Inflation

High Energy Physics - Phenomenology · Submitted 30 Apr 1997 · 46 references

Pith's verdict: The paper that made preheating calculable: the parabolic-scattering recursion is sound and reproduces the numerics, with soft spots confined to the last oscillation or two and to an acknowledged μ calibration.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/9704452 (version 2)
Pith IDWCW7NAI6, Pith's permanent record of this paper and its review
CitationAlexei Starobinsky, Andrei Linde, Lev Kofman. "Towards the Theory of Reheating After Inflation." arXiv:hep-ph/9704452 (1997). https://pith.science/paper/WCW7NAI6
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Reheating after inflation is not a slow one-particle-at-a-time decay but a stochastic resonance, and the authors claim it can be computed analytically even though each growing mode jumps across thousands of resonance bands within a single振荡

The reading

After inflation ends, the inflaton field oscillates about the minimum of its potential and dumps its energy into other fields. The paper asks how that transfer actually happens in the simplest model, a massive inflaton $m^2varphi^2/2$ coupled to a second scalar $chi$ through $g^2varphi^2chi^2$. The authors argue that the answer is not perturbation theory and not the textbook narrow-resonance picture either: because the oscillation amplitude and frequency change quickly as the universe expands, the phase of each $chi$ mode at the moment the inflaton crosses zero is effectively random, so the standard Mathieu-equation band structure does not apply. They call the resulting process stochastic resonance, show that it still amplifies occupation numbers exponentially, and build an analytic method that treats each zero-crossing as a scattering event off a parabolic potential. With that tool they follow the whole history, including the backreaction of the produced particles and their rescattering, and extract when preheating stops, how many particles and how much fluctuation it leaves behind, and whether it can make particles heavier than the inflaton itself.

What carries the argument

The load-bearing device is the reduction of the mode equation near each zero of the inflaton to a universal form. Writing $X_k = a^{3/2}chi_k$ removes the Hubble friction, and near $varphi(t_j)=0$ the frequency becomes $k^2/a^2 + g^2Phi^2 m^2 (t-t_j)^2$, so with $tau = k_*(t-t_j)$ and $kappa = k/(a k_*)$, where $k_* = sqrt{g m Phi}$, one gets $d^2X/dtau^2 + (kappa^2+tau^2)X = 0$. That equation has exact parabolic-cylinder solutions, hence exact reflection and transmission amplitudes $R_k$ and $D_k$ (Eqs. 51–52), which convert into a two-by-two recursion for the Bogoliubov coefficients across one kick. Iterating the recursion gives the change in occupation number per kick, $n_k^{j+1} = e^{-pi kappa^2} + (1+2e^{-pi kappa^2})n_k^j - 2 e^{-pi kappa^2/2}sqrt{1+e^{-pi kappa^2}}sqrt{n_k^j(1+n_k^j)}sinTheta_j$. The sine term is the interference between consecutive kicks; when the phases $Theta_j$ are random, the whole band structure dissolves and what remains is a broad, phase-insensitive amplification with a resonance width $k lesssim k_*/sqrt{pi}$. This is what lets the authors replace a Mathieu stability chart, which fails here, with a single growth index and then integrate over momenta.

Extended reading notes

Core claim

The central claim is that broad parametric resonance in an expanding universe with a quadratic inflaton potential is a stochastic process, and that it is nevertheless tractable analytically. In the usual treatment each mode sits in a fixed instability band of the Mathieu equation and grows at a fixed rate. Here, the authors point out, the parameter $q = g^2Phi^2/4m^2$ falls as $t^{-2}$, so within a single oscillation a mode sweeps through order $sqrt{q}$ different bands, and the phase accumulated between successive zero-crossings of $varphi$ changes by far more than $pi$. The kicks therefore arrive with uncorrelated phases, and the occupation number $n_k$ can fall as well as rise from one kick to the next — something impossible for classical particles in a fixed background. Averaged over the randomness, the growth is still exponential, with an effective index $mu$ of order $0.1$–$0.2$ (the authors quote $mu approx 0.13$). From this they derive the duration of the first stage of preheating, the amplitude $Phi_1$ and $q_1$ at its end, the criterion $g lesssim 3times10^{-4}$ for preheating to finish before backreaction matters, the growth of the effective inflaton mass as $g n_chi/|varphi|$ (a $|varphi|$ term in the potential, not $varphi^2$), and finally a condition $g m varphi Phi gtrsim 4M^2$ under which $chi$ particles of mass $M gg m$ are copiously produced.

Load-bearing premise

Every exponential estimate rests on replacing the wildly fluctuating, phase-dependent growth index of Eq. (81) by one effective constant, $mu approx 0.13$, and then integrating $n_chi propto e^{2mu m t}$ over the entire first stage.

Editorial extensions

If this is right

  • For $g lesssim 3times10^{-4}$ preheating ends while backreaction is still negligible, and the perturbative reheating calculation remains the correct description of the end of the process.
  • For larger $g$ a second stage follows in which the produced $chi$ particles raise the inflaton's effective frequency, shortening each oscillation and making particle production faster rather than slower, until $q sim 1/4$.
  • The variance of $chi$ can reach $sqrt{langle chi^2 rangle} sim 10^{16}$ GeV, orders of magnitude above any thermal value at the same energy density, which can restore symmetry nonthermally and seed topological defects.
  • Particles with mass $M gg m$ can be produced in quantity provided $g m varphi Phi gtrsim 4M^2$; producing $M = 10m$ requires $g gtrsim 0.065$, or about $0.036$ with the relaxed criterion.
  • At the end of resonance the $chi$ particles outnumber nothing: rescattering leaves roughly ten times more $varphi$ particles than $chi$ particles, so the reheating temperature is still set by the later perturbative decay of the $varphi$ remnants, not by preheating.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If $n_k$ really obeys a log-normal distribution, the volume-averaged energy density is dominated by rare patches with anomalously large occupation rather than by the mean, so preheating would be strongly inhomogeneous on sub-horizon scales; that would be a source of non-Gaussian density perturbations and possibly of primordial black holes, which the paper only gestures at.
  • The parabolic-kick formalism depends only on the behaviour of $varphi$ near its zero, so it should transfer directly to any potential where the inflaton crosses zero with a linear zero, giving a general map from $V(varphi)$ to a growth index without ever writing down a Mathieu equation.
  • Because $mu$ is a random variable, the honest object is a Fokker–Planck equation for the distribution of $ln n_k$; replacing it by a single mean $mu$ will systematically underestimate the upper tail, which would push the superheavy-particle mass bound of Eq. (144) upward.
  • Soft $varphi$ modes with $k ll k_*$ are indistinguishable from the coherent condensate over the timescales that matter, so 'number of inflaton particles' is not a clean observable during preheating; field variances, not occupation numbers, are probably the right bookkeeping.
Original abstract

Reheating after inflation occurs due to particle production by the oscillating inflaton field. In this paper we describe the perturbative approach to reheating, and then concentrate on effects beyond the perturbation theory. They are related to the stage of parametric resonance called preheating. It may occur in an expanding universe if the initial amplitude of oscillations of the inflaton field is large enough. We investigate a simple model of a massive inflaton field coupled to another scalar field X. Parametric resonance in this model is very broad. It occurs in a very unusual stochastic manner, which is different from the parametric resonance in the case when the expansion of the universe is neglected. Quantum fields interacting with the oscillating inflaton field experience a series of kicks which occur with phases uncorrelated to each other. We call this process stochastic resonance. We develop the theory of preheating taking into account the expansion of the universe and backreaction of produced particles, including the effects of rescattering. The process of preheating can be divided into several distinct stages. At the first stage the backreaction of created particles is not important. At the second stage backreaction increases the frequency of oscillations of the inflaton field, which makes the process even more efficient than before. Then the effects related to scattering of X-particles terminate the resonance. We calculate the density of X-particles and their quantum fluctuations with all backreaction effects taken into account. This allows us to find the range of masses and coupling constants for which one has efficient preheating. In particular, under certain conditions this process may produce particles with a mass much greater than the mass of the inflaton field.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 11,802 input tokens, $0.0005 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.16 · human participants 0.01 · controlled comparison 0.09 · randomized assignment 0.01 · states quantitative result 0.08 · compares with prior results 0.41 · mentions code or data availability 0.01.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

Reheating after inflation occurs due to particle production by the oscillating inflaton field. In this paper we describe the perturbative approach to reheating, and then concentrate on effects beyond the perturbation theory. They are related to the stage of parametric resonance called preheating. It may occur in an expanding universe if the initial amplitude of oscillations of the inflaton field is large enough. We investigate a simple model of a massive inflaton field coupled to another scalar field X. Parametric resonance in this model is very broad. It occurs in a very unusual stochastic manner, which is different from the parametric resonance in the case when the expansion of the universe is neglected. Quantum fields interacting with the oscillating inflaton field experience a series of kicks which occur with phases uncorrelated to each other. We call this process stochastic resonance. We develop the theory of preheating taking into account the expansion of the universe and backreaction of produced particles, including the effects of rescattering. The process of preheating can be divided into several distinct stages. At the first stage the backreaction of created particles is not important. At the second stage backreaction increases the frequency of oscillations of the inflaton field, which makes the process even more efficient than before. Then the effects related to scattering of X-particles terminate the resonance. We calculate the density of X-particles and their quantum fluctuations with all backreaction effects taken into account. This allows us to find the range of masses and coupling constants for which one has efficient preheating. In particular, under certain conditions this process may produce particles with a mass much greater than the mass of the inflaton field.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Reheating after inflation occurs due to particle production by the oscillating inflaton field. background 1.00 0.06 0.03 asserted 1.00
s1 In this paper we describe the perturbative approach to reheating, and then concentrate on effects beyond the perturbation theory. contribution 0.88 0.76 0.04 asserted 0.80
s2 They are related to the stage of parametric resonance called preheating. background 0.91 0.11 0.04 asserted 0.99
s3 It may occur in an expanding universe if the initial amplitude of oscillations of the inflaton field is large enough. background 0.91 0.11 0.04 hedged possibility 1.00
s4 We investigate a simple model of a massive inflaton field coupled to another scalar field X. method 0.85 0.76 0.02 asserted 0.97
s5 Parametric resonance in this model is very broad. result 0.97 0.72 0.05 asserted 1.00
s6 It occurs in a very unusual stochastic manner, which is different from the parametric resonance in the case when the expansion of the universe is neglected. result 0.96 0.75 0.09 asserted 1.00
s7 Quantum fields interacting with the oscillating inflaton field experience a series of kicks which occur with phases uncorrelated to each other. result 0.71 0.53 0.06 asserted 1.00
s8 We call this process stochastic resonance. interpretation 0.49 0.83 0.09 asserted 0.99
s9 We develop the theory of preheating taking into account the expansion of the universe and backreaction of produced particles, including the effects of rescattering. contribution 0.78 0.94 0.12 asserted 0.99
s10 The process of preheating can be divided into several distinct stages. result 0.51 0.50 0.05 asserted 1.00
s11 At the first stage the backreaction of created particles is not important. result 0.73 0.44 0.03 asserted 1.00
s12 At the second stage backreaction increases the frequency of oscillations of the inflaton field, which makes the process even more efficient than before. result 0.81 0.71 0.09 asserted 1.00
s13 Then the effects related to scattering of X-particles terminate the resonance. result 0.87 0.42 0.07 asserted 1.00
s14 We calculate the density of X-particles and their quantum fluctuations with all backreaction effects taken into account. contribution 0.51 0.96 0.07 asserted 1.00
s15 This allows us to find the range of masses and coupling constants for which one has efficient preheating. result 0.83 0.91 0.17 asserted 0.99
s16 In particular, under certain conditions this process may produce particles with a mass much greater than the mass of the inflaton field. result 0.98 0.82 0.26 hedged possibility 1.00

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 0.95.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1. Oscillations of the field
Figure 2. FIG. 2. Narrow parametric resonance for the field
Figure 3. In the semiclassical regime when the frequency ωk(t) is changing adiabatically, nk is a constant which coincides with an adiabatic invariant. To appreciate the useful￾ness of the introduction of the adiabatic invariant nk, one should compare the evolution of the modes χk with the evolution of the occupation numbers corresponding to each of these modes shown in Figs. 2 and 3. As we see, in the narrow resonance regime χk vigorously oscillates, whereas ln nk grows like a straight line. In the broad resonance regime the field amplification occurs near the points φ(t) = 0 where the proce
Figure 4. FIG. 4. Early stages of parametric resonance in the the
Figure 5. FIG. 5. The same process as in Fig. 4 during a longer
Figure 7. FIG. 7. The structure of the resonance bands for the
Figure 8. FIG. 8. The change of the comoving particle number
Figure 9. FIG. 9. The characteristic exponent
Figure 10. FIG. 10. The characteristic exponent
Figure 12. FIG. 12. Two diagrams for the polarization operator of

Showing 10 of the paper's 12 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:hep-ph/9704452 by the authors.

References
How later papers used it

755 later papers cite it

Evidence: substantial. 748 of the citing papers share no author with it, and 110 of those carry a confident label from Pith's citation classifier.

Pith read 40 of the 1,465 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 12 of 76 papers with this label · 1 by the same authors

By the same authors

How the literature uses it

668 arXiv papers cite this in 1,296 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 734 sentences without a confident label

Contests it 2 papers
Builds on it 3 papers
Uses its method or tool 70 papers
Uses its data or values 3 papers
Compares with it 20 papers
Cites it as support 11 papers
Cites it as background 343 papers
  • Studies of reheating in inflationary cosmology [ref] have taught us that in a dynamical space-time only the broad resonance regime with $q>1$ can lead to a significant enhancement of the generation of the $χ$ particles and of the corresponding fluctuations on super-Hubble scales.

    Preheating a bouncing universe · Broad Resonance with Weak Coupling · 0.84

  • Particle production can occur whenever the effective mass $m_eff$ of $χ$ becomes light [ref] , given that the coupling constant is chosen appropriately.

    Preheating after Multi-field Inflation · Introduction · 0.86

  • Given the shape of the potential ([ref]), very different (p)reheating mechanisms could in principle take place, from a tachyonic production in the region between the end of inflation and the inflection point [ref] , an instant preheating mechanism [ref] and a parametric resonance effect around the minimum of the potential [ref] .

    Preheating in the Standard Model with the Higgs-Inflaton coupled to gravity · Reheating in the Standard Model of particle physics · 0.80

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.