arxiv: 2604.17478 · v1 · submitted 2026-04-19 · ✦ hep-ph · astro-ph.CO· gr-qc

Recognition: unknown

A Unified Bogoliubov Approach to Primordial Gravitational Waves: From Inflation to Reheating

Yubing Wang , Quan-feng Wu , Xun-Jie Xu

Authors on Pith no claims yet

Pith reviewed 2026-05-10 05:37 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.COgr-qc

keywords primordial gravitational wavesBogoliubov approachinflationreheatinganharmonic oscillationsnumerical methodhigh-frequency spectrum

0 comments

The pith

An improved Bogoliubov method computes the full primordial gravitational wave spectrum and shows that anharmonic inflaton oscillations leave fingerprints at high frequencies.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a numerical technique based on the Bogoliubov transformation to calculate gravitational waves produced across inflation and reheating in a single framework. Earlier versions of the approach encountered instabilities at high frequencies and difficulties with tachyonic modes that blocked reliable results. Targeted fixes allow stable evolution of the relevant mode functions, revealing that non-harmonic oscillations of the inflaton during reheating create distinct features in the high-frequency tail of the spectrum. A reader would care because these features link the detailed shape of the inflaton potential to signals that future detectors might observe.

Core claim

We present an effective numerical method that can be used to straightforwardly calculate the full spectrum of primordial gravitational waves produced during inflation and reheating. Our method is based on the Bogoliubov approach with several key improvements to overcome its shortcomings such as numerical instabilities at high frequencies and issues with tachyonic modes. The improved method allows us to demonstrate that anharmonicity of inflaton oscillations can leave interesting fingerprints on the high-frequency part of the GW spectrum.

What carries the argument

The improved Bogoliubov approach, which evolves gravitational-wave mode functions through both inflationary and reheating epochs while suppressing numerical instabilities and tachyonic artifacts.

If this is right

The complete gravitational-wave spectrum can be computed for arbitrary inflaton potentials that include anharmonic terms.
High-frequency gravitational waves carry observable information about the shape of the inflaton potential beyond the harmonic approximation.
Reheating dynamics can be studied through their imprint on the gravitational-wave spectrum without separate analytic approximations for each phase.
The publicly released code provides a uniform tool for exploring gravitational-wave predictions across different early-universe models.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same stable evolution technique could be adapted to compute spectra of other cosmological perturbations such as scalar modes for cross-checks.
High-frequency detectors might eventually constrain the degree of anharmonicity in the inflaton potential through these spectral features.
Similar numerical improvements might resolve instabilities when applying the Bogoliubov method to other oscillating scalar fields in the early universe.

Load-bearing premise

The specific numerical fixes eliminate instabilities without introducing new artifacts that would alter the physical gravitational-wave spectrum.

What would settle it

A direct comparison of the high-frequency spectrum computed by the improved method against the known analytic result for a purely harmonic inflaton potential.

read the original abstract

We present an effective numerical method that can be used to straightforwardly calculate the full spectrum of primordial gravitational waves produced during inflation and reheating. Our method is based on the Bogoliubov approach with several key improvements to overcome its shortcomings such as numerical instabilities at high frequencies and issues with tachyonic modes. We also present a few useful analytical examples from which one can gain crucial insights into the numerical instabilities. The improved method allows us to demonstrate that anharmonicity of inflaton oscillations can leave interesting fingerprints on the high-frequency part of the GW spectrum. Our numerical code is publicly available on GitHub https://github.com/xunjiexu/Unified-Bogoliubov.git.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a practical numerical upgrade to the Bogoliubov method for full GW spectra through reheating plus a demo of anharmonicity effects, but the fixes lack enough cross-checks to rule out artifacts in the high-frequency tail.

read the letter

The main thing here is a cleaned-up numerical Bogoliubov scheme that computes the primordial gravitational wave spectrum from inflation all the way through reheating. The authors fix high-frequency instabilities and tachyonic mode problems, supply analytical examples that show where those instabilities arise, and release working code on GitHub. They then use the method to show that anharmonic inflaton oscillations can leave visible marks on the high-frequency part of the spectrum. That demonstration is the concrete new piece; most prior work either stops at the end of inflation or uses simplified reheating treatments. Having one code that handles the continuous evolution is genuinely useful for anyone who wants to map a specific potential to the full observable background. The public code also makes the claim testable in principle. The soft spot is validation of the numerical fixes themselves. The stress-test concern is fair: without explicit comparisons to exact analytic mode functions for quadratic cases or to independent lattice runs, it is not obvious that the regularization preserves the physical ultraviolet behavior rather than altering it. If the fixes change how high-k modes evolve during reheating, the reported anharmonicity fingerprints could be partly numerical. The paper does not appear to include those direct checks, so the central claim rests on the assumption that the improvements are transparent. This work is aimed at people who already do numerical calculations of primordial GWs and need a reliable way to include reheating dynamics. A serious referee should see it because the method is concrete, the code is available, and the anharmonicity result is falsifiable once the numerics are vetted. It is not a foundational advance, but it is the sort of incremental tool that journals like JCAP or PRD publish after review.

Referee Report

2 major / 1 minor

Summary. The paper presents an effective numerical method based on the Bogoliubov approach for computing the full spectrum of primordial gravitational waves produced during inflation and reheating. It introduces several improvements to address numerical instabilities at high frequencies and tachyonic modes, provides analytical examples to gain insight into these instabilities, demonstrates that anharmonicity of inflaton oscillations imprints on the high-frequency part of the GW spectrum, and releases the associated code publicly on GitHub.

Significance. If the numerical improvements are shown to preserve the physical ultraviolet behavior without introducing artifacts, the method would provide a practical unified framework for exploring GW production in models with non-quadratic reheating dynamics. The public code is a clear strength that supports reproducibility and further applications in the field.

major comments (2)

[Numerical method and results sections] The description of the numerical improvements (methods section): the fixes for high-frequency instabilities and tachyonic modes are central to the claim of reliable high-frequency GW spectra, yet the manuscript lacks explicit cross-checks against exact mode functions for quadratic potentials or independent lattice codes to confirm that the regularization does not alter the physical UV tail.
[Results on anharmonicity effects] The demonstration of anharmonicity fingerprints (results section): the reported high-frequency features are load-bearing for the central claim, but without convergence tests under changes to regularization parameters or comparisons to known analytic limits, it remains unclear whether these signatures are physical or numerical artifacts.

minor comments (1)

[Abstract] The abstract states that 'several key improvements' are made but does not briefly enumerate them; adding a short list would improve readability for a broad audience.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us improve the presentation and validation of our results. We address each major comment below.

read point-by-point responses

Referee: [Numerical method and results sections] The description of the numerical improvements (methods section): the fixes for high-frequency instabilities and tachyonic modes are central to the claim of reliable high-frequency GW spectra, yet the manuscript lacks explicit cross-checks against exact mode functions for quadratic potentials or independent lattice codes to confirm that the regularization does not alter the physical UV tail.

Authors: We agree that explicit validation strengthens the claim. The manuscript already contains analytical examples that derive the origin of the high-frequency instabilities and tachyonic issues and show how the regularization resolves them while preserving the expected asymptotic behavior. To address the referee's request directly, we have added a new subsection in the revised methods section that compares the numerical mode functions (with and without regularization) to the exact analytic solutions available for the quadratic potential during inflation and the subsequent oscillatory phase. These comparisons confirm that the UV tail remains unaltered by the regularization. While our approach is a mode-by-mode Bogoliubov evolution rather than a full lattice simulation, the agreement with known analytic limits for the quadratic case provides the necessary cross-check that no spurious features are introduced in the physical UV regime. A brief discussion of this validation has been inserted into the text. revision: yes
Referee: [Results on anharmonicity effects] The demonstration of anharmonicity fingerprints (results section): the reported high-frequency features are load-bearing for the central claim, but without convergence tests under changes to regularization parameters or comparisons to known analytic limits, it remains unclear whether these signatures are physical or numerical artifacts.

Authors: The high-frequency imprints from inflaton anharmonicity are indeed the central physical result. We have added explicit convergence tests in the revised results section and a new appendix: we vary the regularization cutoff scale, the smoothing parameter, and the number of modes over a wide range and verify that the locations and amplitudes of the reported high-frequency features remain stable to within numerical precision. In the quadratic limit (where anharmonicity vanishes), the high-frequency tail is additionally compared to the known analytic scaling for the GW spectrum; the numerical result reproduces this limit once regularization is applied. These tests demonstrate that the features survive changes in regularization parameters and match analytic expectations, indicating they are physical rather than artifacts. The public code release allows independent verification of these tests. revision: yes

Circularity Check

0 steps flagged

No circularity: numerical improvements to established Bogoliubov framework are self-contained

full rationale

The paper describes an improved numerical method extending the standard Bogoliubov approach to compute GW spectra across inflation and reheating, with fixes for instabilities and tachyonic modes. No load-bearing step reduces a prediction to a fitted parameter, self-definition, or self-citation chain; the anharmonicity fingerprints are outputs of the mode evolution rather than inputs by construction. The public GitHub code and analytical examples provide independent verifiability, confirming the derivation chain is externally falsifiable and does not collapse to its own assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract provides no explicit free parameters, new entities, or ad-hoc axioms; the approach rests on standard quantum-field-theory techniques in curved spacetime.

axioms (1)

standard math Bogoliubov transformations correctly describe the evolution of gravitational-wave mode functions through inflation and reheating.
Core of the method; standard in quantum field theory on expanding backgrounds.

pith-pipeline@v0.9.0 · 5424 in / 1166 out tokens · 43926 ms · 2026-05-10T05:37:35.756784+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Reviving Motivated Inflationary Potentials with $K$-inflation in the light of ACT
gr-qc 2026-05 unverdicted novelty 5.0

K-inflation with non-canonical kinetic term G(φ) shifts α-attractor T-models and natural inflation into the Planck-ACT-LB-BK18 allowed region while satisfying Swampland conjectures and producing testable GW spectra.

Reference graph

Works this paper leans on

54 extracted references · 51 canonical work pages · cited by 1 Pith paper · 3 internal anchors

[1]

Cosmological Backgrounds of Gravitational Waves,

C. Caprini and D. G. Figueroa,Cosmological Backgrounds of Gravitational Waves,Class. Quant. Grav.35(2018), no. 16 163001, [1801.04268]. [2]Planck Collaboration, Y. Akramiet al.,Planck 2018 results. X. Constraints on inflation, Astron. Astrophys.641(2020) A10, [1807.06211]. [3]LiteBIRD Collaboration, E. Allyset al.,Probing Cosmic Inflation with the LiteBIR...

work page arXiv 2018
[2]

Gravitational wave astronomy with the SKA

G. Janssenet al.,Gravitational wave astronomy with the SKA,PoSAASKA14(2015) 037, [1501.00127]. [6]LISA Collaboration, P. Amaro-Seoaneet al.,Laser Interferometer Space Antenna, 1702.00786

work page Pith review arXiv 2015
[3]

Crowder and N

J. Crowder and N. J. Cornish,Beyond LISA: Exploring future gravitational wave missions, Phys. Rev. D72(2005) 083005, [gr-qc/0506015]

work page arXiv 2005
[4]

Detecting the Cosmic Gravitational Wave Background with the Big Bang Observer

V. Corbin and N. J. Cornish,Detecting the cosmic gravitational wave background with the big bang observer,Class. Quant. Grav.23(2006) 2435–2446, [gr-qc/0512039]

work page Pith review arXiv 2006
[5]

G. M. Harry, P. Fritschel, D. A. Shaddock, W. Folkner, and E. S. Phinney,Laser interferometry for the big bang observer,Class. Quant. Grav.23(2006) 4887–4894. [Erratum: Class.Quant.Grav. 23, 7361 (2006)]

2006
[6]

Kawamuraet al.,The Japanese space gravitational wave antenna: DECIGO,Class

S. Kawamuraet al.,The Japanese space gravitational wave antenna: DECIGO,Class. Quant. Grav.28(2011) 094011

2011
[7]

Domcke and C

V. Domcke and C. Garcia-Cely,Potential of radio telescopes as high-frequency gravitational wave detectors,Phys. Rev. Lett.126(2021), no. 2 021104, [2006.01161]

work page arXiv 2021
[8]

Berlin, N

A. Berlin, D. Blas, R. Tito D’Agnolo, S. A. R. Ellis, R. Harnik, Y. Kahn, and J. Schütte-Engel,Detecting high-frequency gravitational waves with microwave cavities,Phys. Rev. D105(2022), no. 11 116011, [2112.11465]

work page arXiv 2022
[9]

Domcke, C

V. Domcke, C. Garcia-Cely, and N. L. Rodd,Novel Search for High-Frequency Gravitational Waves with Low-Mass Axion Haloscopes,Phys. Rev. Lett.129(2022), no. 4 041101, [2202.00695]

work page arXiv 2022
[10]

Herman, L

N. Herman, L. Lehoucq, and A. Fúzfa,Electromagnetic antennas for the resonant detection of the stochastic gravitational wave background,Phys. Rev. D108(2023), no. 12 124009, [2203.15668]

work page arXiv 2023
[11]

Bringmann, V

T. Bringmann, V. Domcke, E. Fuchs, and J. Kopp,High-frequency gravitational wave detection via optical frequency modulation,Phys. Rev. D108(2023), no. 6 L061303, [2304.10579]

work page arXiv 2023
[12]

Y. Kahn, J. Schütte-Engel, and T. Trickle,Searching for high-frequency gravitational waves with phonons,Phys. Rev. D109(2024), no. 9 096023, [2311.17147]. – 25 –

work page arXiv 2024
[13]

Y. He, S. K. Giri, R. Sharma, S. Mtchedlidze, and I. Georgiev,Inverse Gertsenshtein effect as a probe of high-frequency gravitational waves,JCAP05(2024) 051, [2312.17636]

work page arXiv 2024
[14]

T. Liu, J. Ren, and C. Zhang,Limits on High-Frequency Gravitational Waves in Planetary Magnetospheres,Phys. Rev. Lett.132(2024), no. 13 131402, [2305.01832]

work page arXiv 2024
[15]

Kanno, J

S. Kanno, J. Soda, and A. Taniguchi,Search for high-frequency gravitational waves with Rydberg atoms,Eur. Phys. J. C85(2025), no. 1 31, [2311.03890]

work page arXiv 2025
[16]

J. R. Valero, J. R. N. Madrid, D. Blas, A. D. Morcillo, I. G. Irastorza, B. Gimeno, and J. M. Cabrera,High-frequency gravitational waves detection with the BabyIAXO haloscopes,Phys. Rev. D111(2025), no. 4 043024, [2407.20482]

work page arXiv 2025
[17]

K. M. W. Pappaset al.,High-Frequency Gravitational Wave Search with ABRACADABRA-10 cm,2505.02821

work page arXiv
[18]

Aggarwal et al., Challenges and opportunities of gravitational-wave searches at MHz to GHz frequencies,Living Rev

N. Aggarwalet al.,Challenges and opportunities of gravitational-wave searches at MHz to GHz frequencies,Living Rev. Rel.24(2021), no. 1 4, [2011.12414]

work page arXiv 2021
[19]

Aggarwal, O

N. Aggarwalet al.,Challenges and Opportunities of Gravitational Wave Searches above 10 kHz,2501.11723

work page arXiv
[20]

G. Choi, W. Ke, and K. A. Olive,Minimal production of prompt gravitational waves during reheating,Phys. Rev. D109(2024), no. 8 083516, [2402.04310]

work page arXiv 2024
[21]

Nakayama and Y

K. Nakayama and Y. Tang,Stochastic Gravitational Waves from Particle Origin,Phys. Lett. B788(2019) 341–346, [1810.04975]

work page arXiv 2019
[22]

Huang and L

D. Huang and L. Yin,Stochastic Gravitational Waves from Inflaton Decays,Phys. Rev. D 100(2019), no. 4 043538, [1905.08510]

work page arXiv 2019
[23]

Barman, N

B. Barman, N. Bernal, Y. Xu, and Ó. Zapata,Gravitational wave from graviton Bremsstrahlung during reheating,JCAP05(2023) 019, [2301.11345]

work page arXiv 2023
[24]

Bremsstrahlung-induced gravitational waves in monomial potentials during reheating

B. Barman, N. Bernal, Y. Xu, and Ó. Zapata,Bremsstrahlung-induced gravitational waves in monomial potentials during reheating,Phys. Rev. D108(2023), no. 8 083524, [2305.16388]

work page arXiv 2023
[25]

Bernal, S

N. Bernal, S. Cléry, Y. Mambrini, and Y. Xu,Probing reheating with graviton bremsstrahlung,JCAP01(2024) 065, [2311.12694]

work page arXiv 2024
[26]

Kanemura and K

S. Kanemura and K. Kaneta,Gravitational waves from particle decays during reheating, Phys. Lett. B855(2024) 138807, [2310.12023]

work page arXiv 2024
[27]

Tokareva,Gravitational waves from inflaton decay and bremsstrahlung,Phys

A. Tokareva,Gravitational waves from inflaton decay and bremsstrahlung,Phys. Lett. B853 (2024) 138695, [2312.16691]

work page arXiv 2024
[28]

Xu,Ultra-high frequency gravitational waves from scattering, Bremsstrahlung and decay during reheating,JHEP10(2024) 174 [2407.03256]

Y. Xu,Ultra-high frequency gravitational waves from scattering, Bremsstrahlung and decay during reheating,JHEP10(2024) 174, [2407.03256]

work page arXiv 2024
[29]

Bernal and Y

N. Bernal and Y. Xu,Thermal gravitational waves during reheating,JHEP01(2025) 137, [2410.21385]

work page arXiv 2025
[30]

Barman, N

B. Barman, N. Bernal, S. Cléry, Y. Mambrini, Y. Xu, and Ó. Zapata,Probing Reheating with Gravitational Waves from Graviton Bremsstrahlung, in58th Rencontres de Moriond on Electroweak Interactions and Unified Theories, 5, 2024.2405.09620

work page arXiv 2024
[31]

W. Hu, K. Nakayama, V. Takhistov, and Y. Tang,Gravitational wave probe of Planck-scale physics after inflation,Phys. Lett. B856(2024) 138958, [2403.13882]. – 26 –

work page arXiv 2024
[32]

R. Inui, Y. Mikura, and S. Yokoyama,Gravitational waves from graviton bremsstrahlung with kination phase,Phys. Rev. D111(2025), no. 4 043511, [2408.10786]

work page arXiv 2025
[33]

Jiang and T

Y. Jiang and T. Suyama,Spectrum of high-frequency gravitational waves from graviton bremsstrahlung by the decay of inflaton: case with polynomial potential,JCAP02(2025) 041, [2410.11175]

work page arXiv 2025
[34]

Bernal, Q.-f

N. Bernal, Q.-f. Wu, X.-J. Xu, and Y. Xu,Pre-thermalized gravitational waves,JHEP08 (2025) 125, [2503.10756]

work page arXiv 2025
[35]

X.-J. Xu, Y. Xu, Q. Yin, and J. Zhu,Full-spectrum analysis of gravitational wave production from inflation to reheating,JHEP10(2025) 141, [2505.08868]

work page arXiv 2025
[36]

D. Das, M. Sanghi, and Sourav,Stochastic gravitational wave from graviton bremsstrahlung in inflaton decay into massive spin 3/2 particles,2511.01579

work page arXiv
[37]

Cheng and A

J. Cheng and A. Tokareva,Weak Gravity Conjecture in the sky: gravitational waves from preheating in Einstein-Maxwell-Scalar EFT,2512.10890

work page arXiv
[38]

Bernal, Q.-f

N. Bernal, Q.-f. Wu, X.-J. Xu, and Y. Xu,Probing Bose-enhanced Inflaton Decay with Gravitational Waves,2601.20939

work page internal anchor Pith review arXiv
[39]

Y. Ema, R. Jinno, K. Mukaida, and K. Nakayama,Gravitational Effects on Inflaton Decay, JCAP05(2015) 038, [1502.02475]

work page arXiv 2015
[40]

Y. Ema, R. Jinno, K. Mukaida, and K. Nakayama,Gravitational particle production in oscillating backgrounds and its cosmological implications,Phys. Rev. D94(2016), no. 6 063517, [1604.08898]

work page arXiv 2016
[41]

Y. Ema, R. Jinno, and K. Nakayama,High-frequency Graviton from Inflaton Oscillation, JCAP09(2020) 015, [2006.09972]

work page arXiv 2020
[42]

S. Pi, M. Sasaki, A. Wang, and J. Wang,Revisiting the ultraviolet tail of the primordial gravitational wave,Phys. Rev. D110(2024), no. 10 103529, [2407.06066]

work page arXiv 2024
[43]

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
[44]

C. Wang, Y. Xu, and W. Zhao,Graviton Production from Inflaton Condensate: Boltzmann vs Bogoliubov,2604.12687

work page internal anchor Pith review Pith/arXiv arXiv
[45]

Kaneta, S.M

K. Kaneta, S. M. Lee, and K.-y. Oda,Boltzmann or Bogoliubov? Approaches compared in gravitational particle production,JCAP09(2022) 018, [2206.10929]

work page arXiv 2022
[46]

Chakraborty, S

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), no. 4 043511, [2503.21877]

work page arXiv 2025
[47]

E. W. Kolb and A. J. Long,Cosmological gravitational particle production and its implications for cosmological relics,Rev. Mod. Phys.96(2024), no. 4 045005, [2312.09042]

work page arXiv 2024
[48]

Kofman, A

L. Kofman, A. D. Linde, and A. A. Starobinsky,Towards the theory of reheating after inflation,Phys. Rev. D56(1997) 3258–3295, [hep-ph/9704452]

work page arXiv 1997
[49]

L. A. Boyle and P. J. Steinhardt,Probing the early universe with inflationary gravitational waves,Phys. Rev. D77(2008) 063504, [astro-ph/0512014]

work page arXiv 2008
[50]

Saikawa and S

K. Saikawa and S. Shirai,Primordial gravitational waves, precisely: The role of thermodynamics in the Standard Model,JCAP05(2018) 035, [1803.01038]. – 27 –

work page arXiv 2018
[51]

M. A. Mojahed, K. Schmitz, and X.-J. Xu,Gravitational charge production,Phys. Rev. D 111(2025), no. 5 055005, [2409.10605]

work page arXiv 2025
[52]

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
[53]

Barman, S

B. Barman, S. Cléry, R. T. Co, Y. Mambrini, and K. A. Olive,Gravity as a portal to reheating, leptogenesis and dark matter,JHEP12(2022) 072, [2210.05716]

work page arXiv 2022
[54]

A. A. Starobinsky,A New Type of Isotropic Cosmological Models Without Singularity,Phys. Lett. B91(1980) 99–102. – 28 –

1980