REVIEW 3 major objections 4 minor 24 references
Kinetic Equilibration after Preheating
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Kinetic equilibration after preheating is almost instantaneous on the Hubble time scale.
desk verdict Clever analytic transfer from heavy-ion physics to preheating, but the 'fast kinetic equilibration' claim depends on an unstated lower bound on the χ^4 coupling; referee it with that condition requested. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the exact solution of the relativistic Boltzmann equation for $2\to 2$ number-conserving scatterings in a Friedmann-Lemaître-Robertson-Walker metric, obtained in Ref. [17]. For the special initial spectrum $f_k(0) = \lambda (256/243)(k/T_0) \exp(-4k/(3T_0))$, the distribution at later times is $f_k(\tau) = \lambda e^{-k/[K(\tau)T_0]} K^{-4}(\tau) [4K(\tau) - 3 + (k/(K(\tau)T_0))(1 - K(\tau))]$, with dimensionless time $\tau = \int_{t_0/l_0}^{t/l_0} dt'\, a^{-3}(t')$ and $K(\tau) = 1 - (1/4)e^{-\tau/6}$. The function $K(\tau)$ carries the relaxation: as $\tau$ grows, $K \to 1$ and $f_k$ approaches a classical thermal distribution, with relaxation essentially complete at $\tau \approx 6$. This machinery turns the question of how fast thermalization occurs into a ratio of two timescales — the mean-free time $l_0$ (set by the $\chi$ self-coupling $h$ and preheating parameters) versus the Hubble time $t_0$ — with the exact solution supplying the $O(1)$ factor $N$ that converts $l_0$ into the equilibration duration.
What would settle it
A numerical lattice simulation of broad-resonance preheating that extracts the produced spectrum $f_k$ at the moment the $\chi$ modes become relativistic, and compares it quantitatively with Eq. (3.9) — checking both the infrared plateau and the exponential cutoff near $k \sim m_\phi$ — would test the central approximation. If the true spectrum differs markedly in shape, the exact solution does not apply and the inferred timescale $\delta t \approx N l_0$, along with the e-fold prediction $\Delta N_{\rm pr}$, would need to be revised.
Extended reading notes
Core claim
Using an exact solution of the number-conserving Boltzmann equation for a homogeneous, isotropic expanding universe, the paper shows that the distribution of $\chi$ quanta created by broad-resonance preheating relaxes to kinetic equilibrium on a timescale $\delta t \approx N l_0$, where $l_0$ is the mean-free length fixed by the $\chi$ self-interaction and $N$ is of order one (the exact solution reaches equilibrium at $\tau \approx 6$). Because the ratio of the Hubble time to $l_0$ at the moment the $\chi$ particles become relativistic is large, $t_0/l_0 \approx 7.4\times 10^{15} h^2 (10^{-3}/g)^4 (0.1 M_{\rm Pl}/A)$, kinetic equilibration is effectively instantaneous on the Hubble time scale. Consequently the total number of e-folds from the end of inflation to the completion of kinetic equilibration is $\Delta N_{\rm pr} \approx 4.38 + \ln(g/10^{-3}) + \ln(A/0.1 M_{\rm Pl}) - \ln(m_\phi/(5\times 10^{-6} M_{\rm Pl}))$, and the reheating temperature in this channel is bounded below by $T_r \approx 8.5\times 10^{13}\,{\rm GeV} \times (10^{-3}/g)^{3/4} (m_\phi/(5\times 10^{-6} M_{\rm Pl}))^{5/4} (0.1 M_{\rm Pl}/A)^{1/4}$. The same analysis, the authors argue, applies to tachyonic preheating, which also excites all modes below a critical wavenumber.
Load-bearing premise
The result depends on the assumption that the particle spectrum produced by broad-resonance preheating is well approximated by the special initial distribution $f_k(0) = \lambda (256/243)(k/T_0) \exp(-4k/(3T_0))$, for which the exact Boltzmann solution applies; the paper calls this 'qualitatively similar' and 'a good approximation' but gives no quantitative comparison to an actual preheating spectrum.
Editorial extensions
If this is right
- The reheating temperature in chaotic-inflation models with a $\chi^4$ self-interaction is bounded below by $T_r \approx 8.5\times 10^{13}$ GeV for the fiducial values $g=10^{-3}$, $m_\phi = 5\times 10^{-6} M_{\rm Pl}$, $A = 0.1 M_{\rm Pl}$, much higher than naive perturbative estimates, because kinetic equilibration completes before expansion dilutes the energy.
- The e-fold count from the end of inflation to kinetic equilibration, $\Delta N_{\rm pr} \approx 4.38 + \ln(g/10^{-3}) + \ln(A/0.1 M_{\rm Pl}) - \ln(m_\phi/(5\times 10^{-6} M_{\rm Pl}))$, is a concrete prediction that can feed into observational constraints on inflationary models through the duration of reheating.
- The entropy density of the $\chi$ field reaches its kinetic-equilibrium value on the same short timescale, so relic abundances such as dark matter candidates can in principle be set soon after preheating.
- Including additional particle species or interactions can only shorten the equilibration time, so the timescale computed here is an upper bound for kinetic equilibration in these models.
- The analysis carries over to tachyonic preheating, since that mechanism likewise excites all modes below a critical wavenumber.
Reading between the lines
- If kinetic equilibration is truly near-instantaneous, the radiation-dominated era after preheating begins with a nearly thermal spectrum at the high temperature $T_r$, which would suppress the abundance of massive relics whose production depends on the maximum post-inflationary temperature — a testable consequence in models with gravitinos or other superpartner relics.
- The same exact-solution technique could be exported to other cosmological out-of-equilibrium settings with similar initial spectra, such as gravitational particle production during inflation or thermalization of a hidden dark sector, where the timescale ratio $t_0/l_0$ would give immediate estimates.
- A lattice simulation of broad-resonance preheating that extracts the actual $f_k$ at the moment $\chi$ becomes relativistic, then feeds it into the full Boltzmann equation, would quantify how much the true spectrum deviates from Eq. (3.9); if the relaxation time departs from $\tau \approx 6$ substantially, the fast-equilibration conclusion would need revision.
- Because chemical equilibration and the transfer of remnant inflaton energy are not addressed here, the total thermalization of the universe could take longer than the kinetic timescale computed in this paper.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies kinetic equilibration of particles produced by broad-resonance preheating after inflation. The authors argue that the initial distribution produced by preheating matches the special initial condition of an exact solution of the relativistic Boltzmann equation in an FLRW background obtained in Ref. [17]. Using this solution, they compute the equilibration time scale in terms of the mean free path l0 and compare it with the Hubble time t0. Their central result is Eq. (4.6), t0/l0 ≈ 7.4×10^15 h^2 (10^-3/g)^4 (0.1 M_Pl/A), from which they conclude that kinetic equilibration is much faster than the Hubble expansion. They also derive the number of e-folds from the end of inflation to the onset of kinetic equilibration, Eq. (4.9), and estimate a reheating temperature, Eq. (4.13).
Significance. If the result holds, it would challenge the common perturbative expectation that post-inflationary thermalization takes much longer than a Hubble time, and it would give a concrete, falsifiable prediction for the e-fold duration of the reheating phase. The paper's use of an exact solution of the Boltzmann equation is a strength, and the parametric expressions in Eqs. (4.6) and (4.8) are simple and testable. However, the argument contains a load-bearing gap in the parameter range for the self-coupling h, and the applicability of the exact solution to preheating rests on an unquantified assertion about the shape of the preheating spectrum. These issues must be addressed before the central conclusion can be accepted.
major comments (3)
- [Sec. IV, Eqs. (4.6)-(4.8)] The conclusion that l0 << t0 is not valid for the full parameter range admitted by the paper. Eq. (4.6) gives t0/l0 = 7.4e15 h^2 (1e-3/g)^4 (0.1 M_Pl/A), but the paper only imposes an upper bound on h, h ≲ 8π^2. For example, at g = 1e-3, A = 0.1 M_Pl and h = 1e-8, Eq. (4.6) gives t0/l0 ≈ 0.74, so the mean free path is comparable to the Hubble time. For smaller h, l0 exceeds t0 and Eq. (4.8) (which neglects expansion) is no longer a valid estimate. Since h is a free parameter and the paper provides no lower bound, the central claim of Sec. V that kinetic equilibration is 'fast on the Hubble time scale' does not follow. The authors must either restrict h to a range where Eq. (4.6) gives t0/l0 >> 1, or provide a different argument for the equilibration time that remains valid when l0 is not small compared with t0.
- [Sec. III, between Eqs. (3.9) and (3.13)] The exact solution in Eq. (3.10) applies only to the special initial condition in Eq. (3.9), yet the paper asserts that this distribution is 'qualitatively similar' to the broad-resonance preheating spectrum and 'a good approximation' without providing any quantitative comparison, fit, or error estimate. In fact, the broad-resonance spectrum is approximately flat for k below the cutoff k_c, whereas Eq. (3.9) vanishes linearly as k -> 0, so the claimed similarity is not self-evident. Because the equilibration time scale τ ~ O(1) is a property of this particular exact solution, the entire analysis depends on this unsupported assumption. The authors should provide a direct comparison with the actual preheating spectrum, for example from lattice simulations, and quantify how the equilibration time changes for deviations from Eq. (3.9).
- [Sec. IV, Eq. (4.2)] The central quantitative result depends on the total cross-section σtot, which is taken from the literature with an enhancement factor: σtot ≃ σ0 k^3/(48π^2 nχ), where σ0 = h^2/(16π k^2). This enhancement and the numerical factor are not derived or checked in the present paper. Since l0 = 1/(σtot nχ) and Eq. (4.6) is proportional to h^2, any error or ambiguity in this cross-section propagates directly into the main claim. The authors should either provide a derivation of this medium-enhanced cross-section or at least quantify the uncertainty and show that the conclusion is robust.
minor comments (4)
- [Sec. IV, Eq. (4.5)] The identification t0 = tr = 1/[2H(tr)] is suspicious: during the matter-dominated period tp < t < tr, the cosmic time is related to the Hubble rate by t = 2/(3H), not 1/(2H). The numerical factor in Eq. (4.5) should be checked; the difference is only O(1), but it is an internal inconsistency.
- [Sec. III, Eq. (3.9)] The parameter T0 is stated to be '∼ mφ' but its relation to the preheating parameters g, A, mφ, and k_c is never specified. Since the exact solution depends on T0 through the initial distribution, the authors should define T0 or show that the final results are independent of it.
- [Sec. V, last paragraph] The claim that the analysis 'applies similarly to the case of tachyonic preheating' is asserted without any justification. If the tachyonic spectrum has a different shape, the same initial-condition concern raised in Major Comment 2 applies here as well.
- [Throughout] There are several typographical errors, such as 'themalization' in Sec. IV and 'Lematre' in the author affiliation; a careful proofread is needed.
Circularity Check
No significant circularity: the exact Boltzmann solution is imported from the independent Ref. [17], and Eq. (4.6) is derived from scattering rates rather than fitted to the fast-equilibration conclusion.
full rationale
The paper's central claim—fast kinetic equilibration after preheating—rests on an exact solution of the relativistic Boltzmann equation in FLRW spacetime taken from Bazow et al. (Ref. [17]), a group with no author overlap with the present paper. The special initial distribution Eq. (3.9) is chosen because it admits that exact solution; the paper explicitly labels it an approximation to the broad-resonance spectrum ('qualitatively similar', 'good approximation') rather than a fit to the target result. The equilibration timescale is then controlled by l0 = 1/(sigma_tot n_chi), with n_chi taken from standard preheating back-reaction estimates (Ref. [8]) and sigma_tot from the quartic self-interaction cross-section; Eq. (4.6) computes t0/l0 as a function of h, g, and A. No parameter is adjusted to make Eq. (4.8) or Eq. (4.9) come out; those are consequences of the computed ratio. The dependence of t0/l0 on h^2 means the 'fast' conclusion requires h large enough relative to g^4 A/M_Pl, and the paper's upper bound h ≲ 8*pi^2 alone does not guarantee that; this is a parameter-validity gap (a correctness risk), not circularity. Self-citations to Brandenberger's reviews and to preheating papers are contextual and not load-bearing; the exact-solution step is independent. Scope limits (chemical equilibration, inflaton remnant energy, and the unquantified match of Eq. (3.9) to the actual preheating spectrum) are stated caveats or modeling assumptions, not hidden circular assumptions. No derivation step reduces by construction to its input, so the circularity score is 0.
Assumptions & free parameters
free parameters (6)
- g =
benchmark 1e-3, allowed range 10^-4 to 10^-3
- h =
unspecified, only h ≲ 8 pi^2
- A =
0.1 M_Pl in benchmark
- m_phi =
5e-6 M_Pl in benchmark
- T0 =
~ m_phi
- N =
integer, 1 to 6
assumptions (6)
- ad hoc to paper The broad-resonance preheating spectrum is well approximated by Eq. (3.9).
- domain assumption The exact FLRW Boltzmann solution of Ref. [17] with number-conserving 2 to 2 scatterings applies to the chi gas for t > t_r.
- domain assumption The total cross-section is independent of particle momentum.
- domain assumption Between the end of preheating and t_r the universe is matter-dominated, rho proportional to a^-3.
- domain assumption The enhanced cross-section Eq. (4.2), with sigma_tot equal to sigma0/(48 pi^2) n_chi/k^3, is correct.
- domain assumption Backreaction stops preheating when n_chi(t_p) equals m_phi^2 A/g.
Cite this review
Pith. "Pith review of Kinetic Equilibration after Preheating." pith.science (2026). https://pith.science/paper/UI2LT5X7
@misc{pith2026190809866,
author = {Pith},
title = {Pith review of: Kinetic Equilibration after Preheating},
year = {2026},
howpublished = {\url{https://pith.science/paper/UI2LT5X7}},
note = {Machine review of arXiv:1908.09866}
}
read the original abstract
We study thermal equilibration after preheating in inflationary cosmology, which is an important step towards a comprehensive understanding of cosmic thermal history. By noticing that the problem is parallel to thermalization after a relativistic heavy ion collision, we make use of the methods developed in this context and that seek for an analytical approach to the Boltzmann equation. In particular, an exact solution for number-conserving scatterings is available for the distribution function in a Friedmann-Lema\^{i}tre-Robertson-Walker metric and can be utilized for the spectral evolution of kinetic equilibration process after preheating. We find that thermal equilibration is almost instantaneous on the time scale of the Hubble time. We also make an explicit prediction for the duration (the number of e-folds of expansion) required for this process of thermal equilibration to complete following the end of inflation.
Reference graph
Works this paper leans on
-
[8]
Theory and applications of Mathieu equations
McLachlan, N. W. “Theory and applications of Mathieu equations.” Clarendon, Oxford (1947)
work page 1947
-
[17]
J. Berges, S. Borsanyi and C. Wetterich, “Prether- malization,” Phys. Rev. Lett. 93, 142002 (2004) doi:10.1103/PhysRevLett.93.142002 [hep-ph/0403234]; J. Berges, K. Boguslavski, S. Schlichting and R. Venu- gopalan, “Turbulent thermalization process in heavy-ion collisions at ultrarelativistic energies,” Phys. Rev. D 89, no. 7, 074011 (2014) doi:10.1103/Ph...
arXiv 2004
-
[1]
The reduced Planck mass is denoted by MPl ≡ 1/ √ 8πG ≃ 2.4× 1018GeV. We assume a spatially flat Friedmann-Lemaˆ ıtre-Robertson-Walker (FLRW) metric, ds2 =−dt2 +a(t)2dx2, (1.1) where t is the physical time, x are the comoving spatial coordinates and a(t) is the scale factor, which describes the expansion of space. The Hubble expansion rate is H(t)≡ ˙a a, (1...
work page 2017
-
[2]
The Inflationary Universe: A Possible So- lution To The Horizon And Flatness Problems,
A. H. Guth, “The Inflationary Universe: A Possible So- lution To The Horizon And Flatness Problems,” Phys. Rev. D 23, 347 (1981); R. Brout, F. Englert and E. Gunzig, “The Creation Of The Universe As A Quantum Phenomenon, ” Annals Phys. 115, 78 (1978); K. Sato, “First Order Phase Transition of A Vacuum And Expansion Of The Universe, ” Mon. Not. Roy. Astron....
work page 1981
-
[3]
Warm In- flation and its Microphysical Basis,
A. Berera, I. G. Moss and R. O. Ramos, “Warm In- flation and its Microphysical Basis,” Rept. Prog. Phys. 72, 026901 (2009) doi:10.1088/0034-4885/72/2/026901 [arXiv:0808.1855 [hep-ph]]
arXiv 2009
-
[4]
Baryon Asymmetry in Inflationary Universe,
A. D. Dolgov and A. D. Linde, “Baryon Asymmetry in Inflationary Universe,” Phys. Lett. 116B, 329 (1982). doi:10.1016/0370-2693(82)90292-1
-
[5]
Particle Pro- duction in the New Inflationary Cosmology,
L. F. Abbott, E. Farhi and M. B. Wise, “Particle Pro- duction in the New Inflationary Cosmology,” Phys. Lett. 117B, 29 (1982). doi:10.1016/0370-2693(82)90867-X
-
[6]
Par- ticle Production During Out-of-equilibrium Phase Transitions,
J. H. Traschen and R. H. Brandenberger, “Par- ticle Production During Out-of-equilibrium Phase Transitions,” Phys. Rev. D 42, 2491 (1990). doi:10.1103/PhysRevD.42.2491
Show all 24 references
-
[7]
On Particle Creation By A Time Dependent Scalar Field,
A. D. Dolgov and D. P. Kirilova, “On Particle Creation By A Time Dependent Scalar Field,” Sov. J. Nucl. Phys. 51, 172 (1990) [Yad. Fiz. 51, 273 (1990)]
1990
-
[9]
Re- heating after inflation,
L. Kofman, A. D. Linde and A. A. Starobinsky, “Re- heating after inflation,” Phys. Rev. Lett. 73, 3195 (1994) doi:10.1103/PhysRevLett.73.3195 [hep-th/9405187]
1994 arXiv
-
[10]
Universe reheating after inflation,
Y. Shtanov, J. H. Traschen and R. H. Brandenberger, “Universe reheating after inflation,” Phys. Rev. D 51, 5438 (1995) doi:10.1103/PhysRevD.51.5438 [hep- ph/9407247]
1995
-
[11]
To- wards the theory of reheating after inflation,
L. Kofman, A. D. Linde and A. A. Starobinsky, “To- wards the theory of reheating after inflation,” Phys. Rev. D 56, 3258 (1997) doi:10.1103/PhysRevD.56.3258 [hep- ph/9704452]
1997
-
[12]
Reheating in Inflationary Cosmology: 7 Theory and Applications,
R. Allahverdi, R. Brandenberger, F. Y. Cyr-Racine and A. Mazumdar, “Reheating in Inflationary Cosmology: 7 Theory and Applications,” Ann. Rev. Nucl. Part. Sci. 60, 27 (2010) doi:10.1146/annurev.nucl.012809.104511 [arXiv:1001.2600 [hep-th]]
2010
-
[13]
Nonperturbative Dynamics Of Reheat- ing After Inflation: A Review,
M. A. Amin, M. P. Hertzberg, D. I. Kaiser and J. Karouby, “Nonperturbative Dynamics Of Reheat- ing After Inflation: A Review,” Int. J. Mod. Phys. D 24, 1530003 (2014) doi:10.1142/S0218271815300037 [arXiv:1410.3808 [hep-ph]]
2014 arXiv
-
[14]
Inflationary Fluctuations, Entropy Generation and Baryogenesis,
J. R. Ellis, K. Enqvist, D. V. Nanopoulos and K. A. Olive, “Inflationary Fluctuations, Entropy Generation and Baryogenesis,” Phys. Lett. B 191, 343 (1987). doi:10.1016/0370-2693(87)90620-4; S. Davidson and S. Sarkar, “Thermalization after inflation,” JHEP 0011, 012 (2000) doi:10...
1987 arXiv
-
[15]
Relativistic tur- bulence: A Long way from preheating to equi- librium,
R. Micha and I. I. Tkachev, “Relativistic tur- bulence: A Long way from preheating to equi- librium,” Phys. Rev. Lett. 90, 121301 (2003) doi:10.1103/PhysRevLett.90.121301 [hep-ph/0210202]; R. Micha and I. I. Tkachev, “Turbulent ther- malization,” Phys. Rev. D 70, 043538 (2004)...
2003 arXiv
-
[16]
Parametric resonance in quan- tum field theory,
J. Berges and J. Serreau, “Parametric resonance in quan- tum field theory,” Phys. Rev. Lett. 91, 111601 (2003) doi:10.1103/PhysRevLett.91.111601 [hep-ph/0208070]; J. Berges, D. Gelfand and J. Pruschke, “Quantum theory of fermion production after inflation,” Phys. Rev. Lett. 107,...
2003 arXiv
-
[18]
Nonlinear dynamics from the rela- tivistic Boltzmann equation in the Friedmann-Lematre- Robertson-Walker spacetime,
D. Bazow, G. S. Denicol, U. Heinz, M. Martinez and J. Noronha, “Nonlinear dynamics from the rela- tivistic Boltzmann equation in the Friedmann-Lematre- Robertson-Walker spacetime,” Phys. Rev. D 94, no. 12, 125006 (2016) doi:10.1103/PhysRevD.94.125006 [arXiv:1607.05245 [hep-ph]]
2016 arXiv
-
[19]
Chaotic Inflation,
A. D. Linde, “Chaotic Inflation,” Phys. Lett. 129B, 177 (1983). doi:10.1016/0370-2693(83)90837-7
1983 doi
-
[20]
Theory of Cosmological Perturbations
V.F. Mukhanov, H.A. Feldman and R.H. Brandenberger, “Theory of Cosmological Perturbations” Physics Reports 215, 203 (1992)
1992
-
[21]
Lectures on the theory of cosmo- logical perturbations
R. H. Brandenberger, “Lectures on the theory of cosmo- logical perturbations” Lect. Notes Phys. 646, 127 (2004) [arXiv:hep-th/0306071]
2004 arXiv
-
[22]
An Introduction to quantum field theory,
M. E. Peskin and D. V. Schroeder, “An Introduction to quantum field theory,” (Addison-Wesley, NY, 1995)
1995
-
[23]
Reheating constraints to inflationary models,
L. Dai, M. Kamionkowski and J. Wang, “Reheating constraints to inflationary models,” Phys. Rev. Lett. 113, 041302 (2014) doi:10.1103/PhysRevLett.113.041302 [arXiv:1404.6704 [astro-ph.CO]]
2014 arXiv
-
[24]
Dynamics of symmetry breaking and tachyonic preheating,
G. N. Felder, J. Garcia-Bellido, P. B. Greene, L. Kofman, A. D. Linde and I. Tkachev, “Dynamics of symmetry breaking and tachyonic preheating,” Phys. Rev. Lett.87, 011601 (2001) doi:10.1103/PhysRevLett.87.011601 [hep- ph/0012142]; E. J. Copeland, S. Pascoli and A. Rajantie, “D...
2001
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.