Recognition: unknown
Graviton Production from Inflaton Condensate: Boltzmann vs Bogoliubov
Pith reviewed 2026-05-10 14:57 UTC · model grok-4.3
The pith
For steeper inflaton potentials, Bogoliubov formalism captures sizable graviton production from non-adiabatic transitions missed by Boltzmann methods.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For inflaton potentials steeper than quadratic the Bogoliubov formalism naturally incorporates a sizable graviton contribution arising from the non-adiabatic transition between inflation and reheating; this contribution is absent in the Boltzmann description and remains important over a broad range of momenta.
What carries the argument
Bogoliubov coefficients obtained by matching mode functions across the abrupt change in effective frequency at the onset of reheating.
If this is right
- The total graviton energy density receives a non-perturbative addition whose scaling with n can be read from the derived analytic expressions.
- Boltzmann equations remain accurate only for quadratic potentials or for very short-wavelength modes.
- Analytic approximations in both frameworks clarify how the spectrum shape depends on the potential index n.
- Reheating models with n greater than 2 must employ the full Bogoliubov treatment to obtain correct gravitational-wave predictions.
Where Pith is reading between the lines
- Primordial gravitational-wave searches could distinguish inflaton potential shapes near the minimum through the low-momentum part of the spectrum.
- Analogous non-adiabatic contributions may appear in the production of other light degrees of freedom during the same transition.
- The regime boundary between perturbative and non-perturbative production can be mapped by varying the duration of the transition rather than assuming an instantaneous switch.
Load-bearing premise
The inflaton is treated as a homogeneous classical oscillating background whose potential near the minimum changes suddenly when inflation ends.
What would settle it
A direct numerical integration of the mode equations for an n=4 potential that yields a graviton spectrum identical to the Boltzmann prediction with no excess at intermediate momenta.
read the original abstract
We study graviton production from an oscillating inflaton condensate during reheating by systematically comparing Boltzmann and Bogoliubov descriptions for inflaton potentials of the form $V(\phi)\propto\phi^n$ around the minimum. The Bogoliubov framework provides a unified description of graviton production, capturing both perturbative and non-perturbative effects across short and long wavelengths, whereas the Boltzmann approach is restricted to perturbative production at short wavelengths. For the quadratic case ($n=2$), we find that the two approaches yield identical graviton spectra at short wavelengths, indicating that the Boltzmann treatments fully captures perturbative gravitational production in this regime. For steeper potentials ($n>2$), however, we identify a sizable contribution arising from the non-adiabatic transition between inflation and reheating. This component is naturally incorporated in the Bogoliubov formalism but absent in the Boltzmann description, and we show that it is important over a broad range of momenta. We derive analytic approximations within both frameworks that clarify the physical origin and scaling behavior of the spectrum. Our results delineate the regime of validity of Boltzmann approaches and show that, for steeper inflaton potentials, graviton production is governed by non-adiabatic transition dynamics for which the Bogoliubov formalism provides the most appropriate description.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript compares Boltzmann and Bogoliubov treatments of graviton production from an oscillating inflaton condensate with V(φ) ∝ φ^n near the minimum. For n=2 the two frameworks agree on the short-wavelength spectrum. For n>2 the Bogoliubov approach includes a non-adiabatic contribution from the inflation-to-reheating transition that is absent from the Boltzmann collision term and is reported to be sizable over a broad momentum range. Analytic approximations are derived in both frameworks to explain the scaling.
Significance. If the central comparison holds, the work usefully delineates the regime of validity of perturbative Boltzmann methods for gravitational particle production during reheating and shows that non-adiabatic dynamics at the end of inflation can dominate for steeper potentials. This has direct implications for the spectrum of primordial gravitational waves and for the choice of computational framework in early-universe cosmology.
major comments (2)
- [§3] §3 (Bogoliubov formalism) and the transition modeling paragraph: the non-adiabatic contribution for n>2 is obtained by matching mode functions across an instantaneous jump in the effective frequency when the condensate begins oscillating under V(φ)∝φ^n. No comparison is shown with a continuous transition of finite duration set by the inflaton mass scale; such a smoothing would suppress the Bogoliubov mixing for k ≲ a H_end and could shrink or eliminate the reported broad-momentum window.
- [Results] Results section (analytic approximations and spectra): the abstract states that the non-adiabatic term is 'sizable' and 'important over a broad range of momenta' for n>2, yet the manuscript supplies neither explicit numerical spectra with error bands nor a quantitative estimate of the size of the sudden-approximation error. Without these, the load-bearing claim that the excess is non-perturbative and broad cannot be verified.
minor comments (2)
- [Notation] The notation for the effective frequency ω_k(t) is used in both frameworks but never tabulated side-by-side, making it difficult to see exactly which term is omitted in the Boltzmann collision integral.
- [Figures] Figure captions for the momentum spectra should explicitly state the value of n and the range of k shown, and whether the curves include only the perturbative piece or the full Bogoliubov result.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and will revise the paper accordingly to strengthen the presentation while preserving the central results.
read point-by-point responses
-
Referee: §3 (Bogoliubov formalism) and the transition modeling paragraph: the non-adiabatic contribution for n>2 is obtained by matching mode functions across an instantaneous jump in the effective frequency when the condensate begins oscillating under V(φ)∝φ^n. No comparison is shown with a continuous transition of finite duration set by the inflaton mass scale; such a smoothing would suppress the Bogoliubov mixing for k ≲ a H_end and could shrink or eliminate the reported broad-momentum window.
Authors: The instantaneous matching is the standard sudden-transition approximation used throughout the literature on reheating and gravitational particle production, as the inflaton begins coherent oscillations on a timescale ~1/m_φ that is short compared with the Hubble time. For the momentum window where we report the non-adiabatic term to be sizable (extending well above a H_end), the violation of adiabaticity remains robust even when a finite-duration smoothing is considered, because the relevant modes oscillate many times during the transition. Nevertheless, we agree that an explicit estimate of the finite-duration correction is useful. In the revised manuscript we will add a short discussion quantifying the suppression for k ≲ a H_end and confirming that the broad-momentum excess for n>2 survives. revision: yes
-
Referee: Results section (analytic approximations and spectra): the abstract states that the non-adiabatic term is 'sizable' and 'important over a broad range of momenta' for n>2, yet the manuscript supplies neither explicit numerical spectra with error bands nor a quantitative estimate of the size of the sudden-approximation error. Without these, the load-bearing claim that the excess is non-perturbative and broad cannot be verified.
Authors: The analytic approximations derived in Sections 4 and 5 already supply explicit scaling relations that quantify the relative size of the non-adiabatic term versus the perturbative Boltzmann contribution; these scalings demonstrate that the non-adiabatic piece dominates for n>2 over a wide interval of comoving momenta. The numerical Bogoliubov spectra are obtained by direct integration of the mode equations. To make the evidence more transparent we will (i) overlay error bands on the numerical spectra that reflect integration tolerances and (ii) add a paragraph estimating the sudden-approximation error via the ratio of the transition timescale to the mode frequency, thereby providing the quantitative support requested. revision: yes
Circularity Check
No significant circularity; frameworks compared independently
full rationale
The paper derives graviton spectra separately in the Boltzmann collision-integral approach and the Bogoliubov mode-function approach, then compares them. For n=2 the two agree at short wavelengths by explicit calculation; for n>2 the extra non-adiabatic piece appears only in the Bogoliubov matching across the inflation-to-reheating transition. Neither result is obtained by fitting a parameter to the target spectrum, nor by renaming a known result, nor by a load-bearing self-citation. The sudden-transition modeling choice is an approximation whose accuracy can be tested externally, but it does not render the derivation self-referential. The central claim therefore remains independent of its inputs.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Bogoliubov transformation for bosonic mode functions in a time-varying background
- domain assumption Boltzmann equation applies to perturbative graviton production at short wavelengths
Forward citations
Cited by 1 Pith paper
-
A Unified Bogoliubov Approach to Primordial Gravitational Waves: From Inflation to Reheating
An improved Bogoliubov numerical method computes the full primordial GW spectrum from inflation to reheating and shows that inflaton anharmonicity imprints distinctive features at high frequencies.
Reference graph
Works this paper leans on
-
[1]
Starobinsky,A New Type of Isotropic Cosmological Models Without Singularity,Phys
A.A. Starobinsky,A New Type of Isotropic Cosmological Models Without Singularity,Phys. Lett. B91(1980) 99
1980
-
[2]
Guth,The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems,Phys
A.H. Guth,The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems,Phys. Rev. D23(1981) 347
1981
-
[3]
Linde,A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems,Phys
A.D. Linde,A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems,Phys. Lett. B108(1982) 389
1982
-
[4]
Albrecht and P.J
A. Albrecht and P.J. Steinhardt,Cosmology for Grand Unified Theories with Radiatively Induced Symmetry Breaking,Phys. Rev. Lett.48(1982) 1220
1982
-
[5]
Starobinsky,Spectrum of relict gravitational radiation and the early state of the universe,JETP Lett.30(1979) 682
A.A. Starobinsky,Spectrum of relict gravitational radiation and the early state of the universe,JETP Lett.30(1979) 682
1979
-
[6]
Allen,The Stochastic Gravity Wave Background in Inflationary Universe Models, Phys
B. Allen,The Stochastic Gravity Wave Background in Inflationary Universe Models, Phys. Rev. D37(1988) 2078
1988
-
[7]
Sahni,The Energy Density of Relic Gravity Waves From Inflation,Phys
V. Sahni,The Energy Density of Relic Gravity Waves From Inflation,Phys. Rev. D 42(1990) 453
1990
-
[8]
M.S. Turner, M.J. White and J.E. Lidsey,Tensor perturbations in inflationary models as a probe of cosmology,Phys. Rev. D48(1993) 4613 [astro-ph/9306029]
-
[9]
Cosmological Backgrounds of Gravitational Waves,
C. Caprini and D.G. Figueroa,Cosmological Backgrounds of Gravitational Waves, Class. Quant. Grav.35(2018) 163001 [1801.04268]
-
[10]
Primordial backgrounds of relic gravitons,
M. Giovannini,Primordial backgrounds of relic gravitons,Prog. Part. Nucl. Phys. 112(2020) 103774 [1912.07065]
-
[11]
Traschen and R.H
J.H. Traschen and R.H. Brandenberger,Particle Production During Out-of-equilibrium Phase Transitions,Phys. Rev. D42(1990) 2491
1990
-
[12]
Y. Shtanov, J.H. Traschen and R.H. Brandenberger,Universe reheating after inflation,Phys. Rev. D51(1995) 5438 [hep-ph/9407247]
- [13]
-
[14]
R. Brandenberger and J. Traschen,Inflationary Reheating to Preheating - A Personal Account,2602.16963
- [15]
- [16]
- [17]
-
[18]
M. Alsarraj and R. Brandenberger,Moduli and graviton production during moduli stabilization,JCAP09(2021) 008 [2103.07684]
- [19]
-
[20]
Giovannini,Ultra-high frequency spikes of relic gravitons,Class
M. Giovannini,Ultra-high frequency spikes of relic gravitons,Class. Quant. Grav.42 (2025) 175021 [2503.12175]
-
[21]
E.W. Kolb and A.J. Long,Cosmological gravitational particle production and its implications for cosmological relics,Rev. Mod. Phys.96(2024) 045005 [2312.09042]
- [22]
-
[23]
Y. Xu,Ultra-high frequency gravitational waves from scattering, Bremsstrahlung and decay during reheating,JHEP10(2024) 174 [2407.03256]
-
[24]
N. Bernal, Q.-f. Wu, X.-J. Xu and Y. Xu,Pre-thermalized gravitational waves, JHEP08(2025) 125 [2503.10756]
- [25]
-
[26]
K. Kaneta, S.M. Lee and K.-y. Oda,Boltzmann or Bogoliubov? Approaches compared in gravitational particle production,JCAP09(2022) 018 [2206.10929]
-
[27]
E. Basso, D.J.H. Chung, E.W. Kolb and A.J. Long,Quantum interference in gravitational particle production,JHEP12(2022) 108 [2209.01713]
-
[28]
A. Chakraborty, S. Clery, M.R. Haque, D. Maity and Y. Mambrini,Generalizing the Bogoliubov vs Boltzmann approaches in gravitational production,Phys. Rev. D112 (2025) 043511 [2503.21877]
-
[29]
N. Bhusal, M.E.C. M., M.A.G. Garcia, A.G. Menkara and M. Pierre,Fermion reheating with a quartic inflaton potential,2512.16203
-
[30]
Creation of spin-3/2 dark matter via cosmological gravitational particle production
E.W. Kolb, A.J. Long, E. McDonough and J. Wang,Creation of spin-3/2 dark matter via cosmological gravitational particle production,2512.16976
work page internal anchor Pith review Pith/arXiv arXiv
-
[31]
High Frequency Spectrum of Primordial Gravitational Waves
K. Mudrunka and K. Nakayama,High Frequency Spectrum of Primordial Gravitational Waves,2601.00378
work page internal anchor Pith review Pith/arXiv arXiv
-
[32]
Harris, K.J
C.R. Harris, K.J. Millman, S.J. van der Walt, R. Gommers, P. Virtanen, D. Cournapeau et al.,Array programming with NumPy,Nature585(2020) 357
2020
-
[33]
Hunter,Matplotlib: A 2d graphics environment,Computing in Science & Engineering9(2007) 90
J.D. Hunter,Matplotlib: A 2d graphics environment,Computing in Science & Engineering9(2007) 90
2007
-
[34]
Rackauckas and Q
C. Rackauckas and Q. Nie,Differentialequations.jl–a performant and feature-rich ecosystem for solving differential equations in julia,Journal of Open Research Software5(2017) 15
2017
-
[35]
Ellis,Ti k z-feynman: Feynman diagrams with ti k z,Computer Physics Communications210(2017) 103–123
J.P. Ellis,Ti k z-feynman: Feynman diagrams with ti k z,Computer Physics Communications210(2017) 103–123. – 39 –
2017
-
[36]
Shtabovenko, R
V. Shtabovenko, R. Mertig and F. Orellana,Feyncalc 10: Do multiloop integrals dream of computer codes?,Computer Physics Communications306(2025) 109357
2025
-
[37]
Universality Class in Conformal Inflation
R. Kallosh and A. Linde,Universality Class in Conformal Inflation,JCAP07 (2013) 002 [1306.5220]
work page Pith review arXiv 2013
-
[38]
Non-minimal Inflationary Attractors,
R. Kallosh and A. Linde,Non-minimal Inflationary Attractors,JCAP10(2013) 033 [1307.7938]
-
[39]
M. Drees and Y. Xu,Small field polynomial inflation: reheating, radiative stability and lower bound,JCAP09(2021) 012 [2104.03977]
-
[40]
M. Drees and Y. Xu,Large field polynomial inflation: parameter space, predictions and (double) eternal nature,JCAP12(2022) 005 [2209.07545]
-
[41]
Turner,Coherent Scalar Field Oscillations in an Expanding Universe,Phys
M.S. Turner,Coherent Scalar Field Oscillations in an Expanding Universe,Phys. Rev. D28(1983) 1243
1983
- [42]
-
[43]
D.G. Figueroa, A. Florio, F. Torrenti and W. Valkenburg,The art of simulating the early Universe – Part I: Integration techniques and canonical cases,JCAP04(2021) 035 [2006.15122]
-
[44]
E.W. Kolb and M.S. Turner,The Early Universe, vol. 69, Taylor and Francis (5, 2019), 10.1201/9780429492860
-
[45]
B.R. Holstein,Graviton Physics,Am. J. Phys.74(2006) 1002 [gr-qc/0607045]
-
[46]
Garcia, K
M.A. Garcia, K. Kaneta, Y. Mambrini and K.A. Olive,Inflaton oscillations and post-inflationary reheating,Journal of Cosmology and Astroparticle Physics2021 (2021) 012
2021
-
[47]
Gorbunov and V.A
D.S. Gorbunov and V.A. Rubakov,Introduction to the Theory of the Early Universe: Cosmological Perturbations and Inflationary Theory, World Scientific, Singapore ; Hackensack, N.J (2011)
2011
-
[48]
de Garcia Maia,Spectrum and energy density of relic gravitons in flat Robertson-Walker universes,Phys
M.R. de Garcia Maia,Spectrum and energy density of relic gravitons in flat Robertson-Walker universes,Phys. Rev. D48(1993) 647
1993
-
[49]
K. Ichikawa, T. Suyama, T. Takahashi and M. Yamaguchi,Primordial Curvature Fluctuation and Its Non-Gaussianity in Models with Modulated Reheating,Phys. Rev. D78(2008) 063545 [0807.3988]
- [50]
-
[51]
German,On theα-attractor T-models,JCAP09(2021) 017 [2105.05426]
G. German,On theα-attractor T-models,JCAP09(2021) 017 [2105.05426]
-
[52]
N. Aggarwal et al.,Challenges and opportunities of gravitational-wave searches above 10 kHz,Living Rev. Rel.28(2025) 10 [2501.11723]. – 40 –
- [53]
-
[54]
P. de Aquino, K. Hagiwara, Q. Li and F. Maltoni,Simulating graviton production at hadron colliders,JHEP06(2011) 132 [1101.5499]
-
[55]
B. Barman, N. Bernal, Y. Xu and ´O. Zapata,Gravitational wave from graviton Bremsstrahlung during reheating,JCAP05(2023) 019 [2301.11345]. [57]Atacama Cosmology Telescopecollaboration,The Atacama Cosmology Telescope: DR6 power spectra, likelihoods andΛCDM parameters,JCAP11(2025) 062 [2503.14452]. [58]BICEP, Keckcollaboration,Improved Constraints on Primor...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.