{"id":"10c70ee3-3ea7-40d0-a4a7-fa0daafe0209","arxiv_id":"2504.17533","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"Thermal gravitons from the Gibbons-Hawking bath at the end of inflation are predicted to yield a MHz-band relic gravitational wave background with Ω_GW h^2 around 10^-18.","lead":"The authors predict a high-frequency relic gravitational wave background from thermal gravitons in the de Sitter horizon, released when inflation ends and redshifted to the MHz band. If correct, this signal would offer a new probe of the reheating temperature, but the predicted amplitude lies orders of magnitude below current detector sensitivities.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The core assumption that the dS detector thermal response is a releasable graviton energy density is unproven; γ≈1 is asserted, and Sec. IV numbers are also inconsistent with Eqs. (21)-(22).","rationale":"The most load-bearing condition for the paper's central claim is not the cosmological redshift algebra; it is the step from detector thermalization in Ref. [24] to a positive, releasable graviton energy density. The paper's own Sec. V admits that the release efficiency is not derived and that a full Bogoliubov calculation is future work; that admission attaches directly to the central number. I agree with the reader's weakest_assumption. I also verified that the Sec. IV benchmark scaling is inconsistent with Eqs. (16)-(22): for fixed T_reh, ρ_G,0∝H^{4/3} rather than H^4, and f_peak,0∝H^{1/3} rather than H, so the nominal high-H, high-T_reh case falls near 2 MHz with Ω_Gh²≈few×10^{-20}, not the quoted 10^8 Hz and 10^{-18}. This makes the current quantitative predictions untrustworthy, but it is a correctable calculational problem; the unresolved physical premise justifies keeping the conditional verdict rather than moving to accept or reject on the present text.","tokens_in":11736,"tokens_out":29268,"duration_ms":286724,"concrete_test":"Compute the late-time graviton energy density from linearized transverse-traceless metric perturbations in a background where H(t) drops from H_Λ to radiation-era values on a rapid (step or tanh) transition. Impose Bunch-Davies initial conditions, solve the mode equation v_k''+(k²-a''/a)v_k=0, extract the late-time Bogoliubov occupation numbers n_k, and integrate ρ_grav=∫(d³k/(2π)³)(k/a)n_k. Compare this with γ(π²/15)T_H^4(a_*/a)^4; if the produced energy is orders of magnitude smaller, the predicted γ≈1 background does not exist, whereas a match would validate the central release assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that Eq. (12) describes a real, uniform reservoir of gravitons, ρ_G=(π²/15)T_H^4, which decouples when the dS horizon vanishes. What Ref. [24] actually supplies is a KMS spectral density seen by a localized detector; that establishes detector thermalization, not the existence of a drainable field energy density. The Rindler/Unruh analog shows the logical gap: a detector can see a thermal bath while the vacuum stress-energy contains no such extractable fluid. The paper's own Sec. II and Sec. V concede that a full Bogoliubov treatment is needed and that the efficiency γ may be <1, but the predictions proceed with γ≈1 based on three plausibility arguments (rapid end, flatness, tiny graviton cross-section), none of which convert vacuum fluctuations into on-shell gravitons. Independently of this, the quoted numbers do not follow from Eqs. (21)-(22): for fixed T_reh, Eq. (16) implies Ω∝H^{4/3} and f0∝H^{1/3}, whereas Sec. IV lists Ω∝H^4 and f0∝H. A direct evaluation of the H=10^14 GeV, T_reh=10^13 GeV case gives f_peak,0≈2×10^6 Hz and Ω_Gh²≈few×10^{-20}, not 10^8 Hz and 10^{-18}. The quantitative claim therefore needs rederivation even before the release efficiency is settled.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes that thermal gravitons associated with the Gibbons-Hawking temperature T_H=H_Λ/(2π) during inflation, described as a local thermal bath in the quantum-thermodynamic framework of Alicki et al., are released at the end of inflation as free radiation with an assumed efficiency γ≈1. The authors derive the present-day spectral density Ω_G(f0) and peak frequency f_peak,0 by propagating the released radiation first through a matter-dominated reheating phase and then through the standard radiation- and matter-dominated eras. They present benchmarks for H_Λ=10^12,10^13,10^14 GeV and T_reh=10^9,10^13 GeV, concluding that the background peaks near MHz with log10(Ω_Gh^2)∼O(-18).","tokens_in":12094,"tokens_out":32506,"duration_ms":253619,"significance":"The proposal is conceptually interesting and, if correct, would provide a new observational window into the reheating temperature. The paper's framework is transparent, uses standard cosmological evolution, and makes falsifiable spectral-shape predictions without fitting to the putative signal. The main strength is the clean analytical propagation of a Planckian spectrum from the inflationary epoch to today. However, the quantitative benchmarks are inconsistent with the paper's own equations, and the assumed release efficiency γ≈1 is not derived; both issues must be addressed before the prediction can be taken at face value.","major_comments":[{"comment":"The numerical benchmarks do not follow from the stated formulas. From Eq. (16), ρ_reh_G=ρ*_G(a*/areh)^4 with a*/areh=(T_reh^4/(3H_Λ^2 m_p^2))^{1/3}, so for fixed T_reh the present-day abundance scales as Ω_G ∝ H_Λ^{4/3}; Eq. (22) gives f_peak,0 ∝ H_Λ^{1/3}. The table in Sec. IV instead lists Ω decreasing by 10^4 per decade in H_Λ (Ω ∝ H_Λ^4) and f_peak,0 ∝ H_Λ. A direct evaluation of the H_Λ=10^14 GeV, T_reh=10^13 GeV case with the paper's inputs gives f_peak,0 ~ 1 MHz and Ω_Gh^2 ~ 10^{-24}, rather than the quoted 10^8 Hz and 10^{-18}. The abstract's O(-18) claim and the Sec. IV benchmarks must be recomputed.","section":"Sec. IV, Eqs. (16) and (22)"},{"comment":"The central premise that the dS thermal state is a real, uniformly distributed, drainable graviton energy density that is released with γ≈1 is not established. Ref. [24] demonstrates KMS thermalization of a localized detector in the dS vacuum; it does not establish that the vacuum contains a reservoir of on-shell gravitons with ρ_G=(π²/15)T_H^4 that converts into propagating radiation when the horizon disappears. The Rindler/Unruh analogue shows that detector thermalization does not imply an extractable fluid. The three plausibility arguments given for γ≈1 (rapid end, flatness, small graviton cross-section) do not amount to a derivation of the conversion efficiency. Because the predicted amplitude is linear in γ and the authors themselves admit in Sec. V that a Bogoliubov calculation may yield γ<1, the paper should present γ as a free parameter and the spectrum as a template, or supply a microscopic estimate of γ.","section":"Sec. II, Eqs. (11)-(12), and Sec. V"},{"comment":"The intermediate formula Ω_G = f/(c²ρ_crit) dρ_G/dln f ... contains an extra factor f relative to the definition Ω_G = (1/(c²ρ_crit)) dρ_G,0/dln f0 in Eq. (19). Since dρ_G/dln f already includes the Jacobian factor f, this makes the intermediate expression dimensionally inconsistent; the final simplified expression is valid only if the extra f is removed. Please correct the derivation.","section":"Eq. (21), middle expression"}],"minor_comments":[{"comment":"The phrase 'using the relation 1/dlnf = f d/df' is incorrect notation; it should read 'd/dln f = f d/df'.","section":"Eq. (21)"},{"comment":"The sentence 'ρ_crit = 3H0²/(8πG) represents the current critical energy density' is dimensionally a mass density in SI units; the c² factor in the definition of Ω resolves the units, but the wording should be clarified.","section":"Sec. III"},{"comment":"The figure caption does not state that all curves assume γ=1; this should be made explicit.","section":"Fig. 2"},{"comment":"No detector sensitivity curves are shown despite the reference to Ref. [25]; adding a representative sensitivity curve or a quantitative statement of the required sensitivity would improve the discussion of observational prospects.","section":"Sec. IV"},{"comment":"The phrase 'captures the conservation of the horizon energy density at the moment of transition' is unclear; the released radiation is not conserved in density but redshifts, so the wording should be revised.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The core issue is the γ=1 assumption. If the authors cannot supply a microscopic derivation, I would suggest they reframe the paper as a phenomenological template with γ a free parameter and soften the abstract accordingly. My recommendation is major revision rather than reject because the calculations are transparent and the idea is potentially testable, but the quantitative benchmarks must be corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's actual novelty is a concrete one: it asks what happens to the sub-horizon thermal gravitons in de Sitter when inflation ends, and treats their release as a snapshot of the Gibbons-Hawking bath. That's a different question from the usual frozen super-horizon spectrum, and I don't see a prior paper making this specific prediction of a present-day Planckian relic whose peak frequency tracks the reheating temperature. The authors build on Alicki et al.'s open-system derivation and produce a simple analytic formula for the redshifted spectrum, including the w≈0 reheating phase. That part is clear and well organized.\n\nThe soft spots are real, and they're in two places. First, the physical premise: the paper needs the dS detector response to be a real, uniform, drainable energy density ρ_G=(π²/15)T_H^4 that suddenly becomes free radiation. That's asserted, not derived. The Unruh analogy is a fair concern—detector thermalization doesn't guarantee an extractable reservoir. The paper itself concedes in Sec. V that a full Bogoliubov treatment is missing and that γ could be <1. So the amplitude is hostage to an efficiency factor that is currently a guess. The authors are honest about this, but it's still a load-bearing guess.\n\nSecond, and more immediate: the numbers in Sec. IV don't follow from their own equations. For the headline case H=10^14 GeV, T_reh=10^13 GeV, feeding Eq. (22) and the entropy ratio gives a peak frequency today around a few MHz, not 10^8 Hz. The energy density today comes out around log10(Ωh²) ~ -23, not -18. The scaling they quote (Ω∝H^4, f0∝H) is off; the correct scaling from their own formulas for fixed T_reh is Ω∝H^{4/3}, f0∝H^{1/3}. This looks like an algebraic slip in evaluating the redshift—possibly treating the reheating phase as radiation-like (ρ∝a^{-4}) rather than the stated w=0 (ρ∝a^{-3}). I won't flag the c^5 in Eq. (23) as dimensional trouble; that's actually consistent if ρ_crit is the mass density and the 1/c^2 is in the definition. The real issue is the arithmetic.\n\nBottom line: the qualitative idea—a HF Planckian relic from the dS thermal bath—survives a rederivation, but the headline values need to be recomputed and the physical premise needs either a supporting calculation or a clear 'speculative scenario' framing. I'd send it to a referee: it's worth one round of serious work, and the core question (can the dS vacuum thermal energy be released as on-shell gravitons?) is a good question for the community. I wouldn't cite it in its current form.","headline":"A novel idea with a credible qualitative prediction, but the headline numbers don't survive contact with the paper's own equations, and the release efficiency is a guess.","tokens_in":12618,"tokens_out":23747,"would_cite":false,"duration_ms":190077,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that inflation ends with a thermal graviton bath that decouples into a Planckian gravitational-wave background peaking in the megahertz band.","keywords":["gravitational waves","high-frequency gravitational waves","de Sitter thermal state","horizon temperature","inflation","reheating temperature","primordial gravitational wave background","Planckian spectrum"],"falsifier":"A full mode-by-mode quantum-field calculation of the transition from the de Sitter phase to a radiation-dominated universe, tracking the mode-mixing coefficients of sub-horizon graviton modes, would give the actual release spectrum and efficiency; if it yields a conversion efficiency much smaller than one, or a spectrum that is not Planckian, the predicted peak abundance and spectral shape fail. Observationally, a sufficiently sensitive broadband detector in the $10^6$ to $10^8$ Hz band that sees no Planckian peak at the predicted abundance would rule the mechanism out, assuming the inflationary Hubble scale is at the benchmark values.","tokens_in":1970,"feed_emoji":"📡","tokens_out":2905,"duration_ms":73520,"temperature":0.7,"pith_summary":"The paper argues that at the end of inflation the universe is filled with a real thermal bath of sub-horizon gravitons at the de Sitter horizon temperature, and when the horizon disappears this bath is released as freely propagating radiation. Tracking that radiation through reheating and cosmological redshift, the authors predict a relic gravitational-wave background with a Planckian spectrum peaking in the megahertz band. For typical inflationary Hubble scales and reheating temperatures the peak abundance is around log10 Omega_G $h^{2}$ about -18 to -32, far below current sensitivity but within the reach of proposed high-frequency detectors. The signal would act as a thermometer for the reheating epoch, something the standard scale-invariant vacuum-fluctuation background cannot provide.","feed_headline":"Inflation's heat leaves a gravitational-wave relic near a megahertz","feed_subtitle":"If it exists, the signal's peak frequency reads the reheating temperature directly from the early universe.","key_machinery":"The load-bearing object is the spectrum of the thermal graviton bath at temperature $T_H = H_\\Lambda/(2\\pi)$ (the horizon temperature of de Sitter space), with spectral energy density $d\\rho_G/d\\omega = (1/\\pi^2) \\omega^3/(e^{\\omega/T_H} - 1)$ and total energy density $\\rho_G = (\\pi^2/15) T_H^4$. The mechanism is a quantum quench: the non-adiabatic end of inflation freezes the statistical distribution of sub-horizon modes before they can adjust to the new background. After that, radiation conservation ($\\rho \\propto a^{-4}$, $f \\propto a^{-1}$) and entropy conservation across reheating fix today's spectrum. An efficiency factor $\\gamma \\le 1$ parameterizes any loss during conversion, with $\\gamma \\approx 1$ argued from the rapidity of the transition, the flatness of the local geometry, and the tiny graviton cross-section.","core_discovery":"The central claim is that the thermal gravitons in equilibrium inside the de Sitter horizon during inflation do not adiabatically vanish when inflation stops; the rapid geometric transition acts as a quantum quench that releases them as free radiation on a timescale much shorter than a Hubble time. Because gravitons scatter only through Planck-suppressed interactions, the released spectrum keeps its Planckian shape, peaking at $f_{\\rm peak} = (x/2\\pi) T_H$ with $x \\approx 2.8214$, and then redshifts as $a^{-1}$. The today spectral density parameter is $\\Omega_G(f_0) \\simeq (16\\pi^2 \\hbar / c^5 \\rho_{\\rm crit}) f_0^4 F(f_0)$, where $F$ is a redshifted Planck factor whose arguments involve the inflationary Hubble rate, the reheating temperature, and the effective number of entropy degrees of freedom. Numerically, the peak falls near $10^6$ to $10^8$ Hz with $\\Omega_G h^2$ of order $10^{-18}$ for a high Hubble rate and efficient reheating, down to $10^{-32}$ for lower rates and inefficient reheating.","pith_inferences":["If the thermal graviton bath is real, the same quantum-thermodynamics argument should apply to any effectively massless field present during inflation, so similar relic backgrounds could be predicted for axion-like particles or dark photons, possibly with different detection consequences.","A full mode-by-mode calculation of the de Sitter to radiation-dominated transition would supply an actual value for the release efficiency $\\gamma$; if it comes out much smaller than one, the predicted peak abundance would drop proportionally, but the frequency and spectral shape would survive.","The mechanism is not limited to the end of inflation in principle: any cosmological epoch where a horizon-temperature bath is suddenly removed could leave a similar thermal snapshot, which is a testable extension connecting this work to late-time horizon dynamics.","Because the peak frequency is set by the reheating temperature, measuring the shape of the Planckian tail, not just the peak, could distinguish instantaneous reheating from reheating with a prolonged matter-dominated stage."],"forward_implications":["The universe would contain a new high-frequency relic gravitational-wave background with a Planckian spectral shape, distinct from the nearly scale-invariant vacuum-fluctuation signal that dominates at lower frequencies.","A detection of the peak frequency would directly measure the reheating temperature, providing constraints complementary to CMB bounds on the spectral index.","The predicted abundance, with $\\Omega_G h^2 \\sim 10^{-18}$ at best, sits far below the dark-radiation bound from CMB and BBN, so the signal is not ruled out by existing cosmological constraints.","The peak frequency falls in the $10^6$ to $10^8$ Hz range, giving emerging high-frequency gravitational-wave detector concepts a concrete benchmark target.","The overall amplitude scales linearly with the release efficiency $\\gamma$, while the spectral shape and peak frequency do not depend on $\\gamma$."],"supporting_citations":[{"why":"Supplies the open-quantum-system proof that a localized detector in de Sitter space thermalizes at the horizon temperature, establishing the bath as a physical energy density rather than a coordinate artifact.","marker":"[24]"},{"why":"Establishes the cosmological event horizon temperature that the entire prediction is built on.","marker":"[3]"},{"why":"Provides the inflationary Hubble-rate bound used to choose the benchmark values and to argue the signal is consistent with CMB constraints.","marker":"[27]"},{"why":"Supplies the low reheating-temperature benchmark around $10^9$ GeV used for the inefficient reheating case.","marker":"[29]"},{"why":"Supplies a high reheating-temperature benchmark around $10^{13}$ GeV used for the optimistic case.","marker":"[30]"},{"why":"Provides the entropy-conservation relation that converts the spectrum at the end of inflation to today's scale factor and temperature.","marker":"[32]"},{"why":"Defines the high-frequency gravitational-wave detector landscape and sensitivity goals that the predicted band targets.","marker":"[25]"},{"why":"Provides the earlier de Sitter radiation energy-density form $\\rho \\propto T^4$ used to normalize the thermal bath.","marker":"[14]"},{"why":"Recovers the same de Sitter energy-density dependence in a different cosmological context, supporting the universality of the Stefan-Boltzmann form.","marker":"[16]"}],"fun_headline_variants":["Inflation's heat leaves a high-frequency graviton relic near MHz","Thermal gravitons from inflation's end decouple as a MHz relic","End of inflation releases thermal gravitons peaking at MHz","High-frequency graviton fossil from inflation's thermal bath"],"cache_read_input_tokens":14592,"weakest_assumption_plain":"The whole prediction stands on the premise that the horizon temperature of de Sitter space corresponds to a real, locally stored energy density of gravitons that is released as freely propagating radiation, with essentially unit efficiency, when inflation ends.","fun_headline_variants_meta":{"raw":{"variants":["Inflation's heat leaves a high-frequency graviton relic near MHz","Thermal gravitons from inflation's end decouple as a MHz relic","End of inflation releases thermal gravitons peaking at MHz","High-frequency graviton fossil from inflation's thermal bath"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00063,"raw_usage":{"total_tokens":2907,"prompt_tokens":941,"completion_tokens":1966,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":1893}},"tokens_in":557,"tokens_out":1966,"duration_ms":13014,"temperature":1.0,"reasoning_tokens":1893,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:39:54.825272+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full mode-by-mode quantum-field calculation of the transition from the de Sitter phase to a radiation-dominated universe, tracking the mode-mixing coefficients of sub-horizon graviton modes, would give the actual release spectrum and efficiency; if it yields a conversion efficiency much smaller than one, or a spectrum that is not Planckian, the predicted peak abundance and spectral shape fail. Observationally, a sufficiently sensitive broadband detector in the $10^6$ to $10^8$ Hz band that sees no Planckian peak at the predicted abundance would rule the mechanism out, assuming the inflationary Hubble scale is at the benchmark values.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the cosmological event horizon temperature that the entire prediction is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies a high reheating-temperature benchmark around $10^{13}$ GeV used for the optimistic case."},{"cited_title":"Hawking radiation from the cosmological horizon in a FRW universe","cited_arxiv_id":"1007.4044","evidence_quote":"Defines the high-frequency gravitational-wave detector landscape and sensitivity goals that the predicted band targets."}],"review_version":1}