{"id":"30aecb8d-0f2d-48f1-a5b9-d57e2deb6865","arxiv_id":"2412.13968","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":1.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A single-author review of how relic gravitational wave backgrounds could probe the post-inflationary expansion history, highlighting high-frequency signals and the quantum nature of the gravitons.","lead":"This review argues that the cosmic background of relic gravitational waves could be the only direct probe of what the universe did between inflation and the formation of the first nuclei. It emphasizes that very high frequency gravitational wave detectors, in the megahertz to terahertz range, might test whether the early universe really was radiation dominated.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed absolute bound nu_max < THz rests on the unproved, internally inconsistent interpolation Eq. 3.50; if this interpolation fails, the high-frequency part of the central claim is not established.","rationale":"The reader correctly identifies Eq. 3.50 as the weakest link. I agree: the interpolation is not derived, and no numerical estimate of gamma is shown in the text. The problem is more specific than a missing derivation: Eq. 3.50 does not actually reproduce the exponential-suppression relation Eq. 3.49 under the paper's own definition n = |beta|^2, so the THz bound is not established on the basis presented. I do not see a more fundamental flaw in the central diagnostic claim; even if nu_max shifts, relic gravitons would still probe the expansion history, just with a different high-frequency cutoff. Therefore the conditional verdict is appropriate, and my read does not change the reader's verdict.","tokens_in":60147,"tokens_out":12126,"duration_ms":111768,"concrete_test":"Independently integrate the tensor mode equation f'' + [k^2 - a''/a] f = 0 for a smooth de Sitter-to-radiation transition (vary the transition width by a factor of 10) and compute n(nu, tau0) = |beta(nu, tau0)|^2 over k up to several k_max. Fit Eq. 3.50 to the numerical n and check residuals over x in [0.1, 10]; then recompute the BBN integral Eq. 3.47 using the numerical n instead of Eq. 3.50. If the recomputed nu_max remains below THz and Eq. 3.50 fits to about 10%, the concern is resolved; otherwise Eq. 3.51 has no demonstrated absolute validity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the 'absolute' upper bound nu_max < THz (Eq. 3.51). The derivation inserts Eq. 3.50, n(nu, tau0) = gamma x^(nT-3)/(exp(gamma x)-1), which is introduced 'in a suggestive form' with no derivation and no numerical estimate of gamma (the cited refs. are not reproduced). The interpolation is supposed to connect Eq. 3.48 (power law) with Eq. 3.49 (exponential suppression). But with the paper's own definition n = |beta|^2, Eq. 3.49 should read n/(1+n) ~ exp(-gamma x) for large x, whereas Eq. 3.50 gives n ~ gamma x^(nT-3) exp(-gamma x), so it does not reproduce Eq. 3.49 unless special values (gamma=1, nT=3) are chosen; Eq. 3.49 as printed uses |n|^2/(1+|n|^2), making the mismatch worse. Thus Eq. 3.50 is not a neutral interpolation but an unverified model assumption. Since Eq. 3.51 explicitly follows by inserting Eq. 3.50 into Eq. 3.46 and imposing Eq. 3.47, the claimed bound is only as secure as this unproved form. If the true high-frequency tail is different, nu_max can shift well above or below THz; the paper gives no independent derivation. The high-frequency and quantum-sensing sections, and the abstract's aHz-to-THz range, rely on this bound. The rest of the review's diagnostic idea is not necessarily damaged, but the 'absolute' part is conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This topical review argues that the stochastic background of relic gravitons is the only direct probe of the post-inflationary expansion history before big bang nucleosynthesis. It develops the FRW timeline with multiple post-inflationary stages, derives the spectral slopes and characteristic frequencies of the relic graviton energy density, and partitions the spectrum into low-frequency (aHz), intermediate-frequency (nHz), and high-frequency (kHz–THz) bands. The review further claims an absolute upper bound on the maximal frequency, νmax < THz, obtained from the BBN bound and from quantum-mechanical pair production, and proposes a dynamical refractive index for tensor modes as a possible explanation of the pulsar-timing-array excess.","tokens_in":60590,"tokens_out":14292,"duration_ms":132149,"significance":"If its main claims hold, the review provides a coherent framework connecting gravitational-wave observations from the aHz to the THz range with the pre-BBN expansion history, and it gives concrete, falsifiable statements: the standard post-inflationary modification cannot explain the PTA amplitude, whereas the refractive-index model predicts a specific high-frequency branch. Strengths include internally consistent standard derivations of the e-fold shifts, spectral slopes, BBN integrals, and a careful discussion of unitarity in the WKB approximation. However, the absolute high-frequency bound and the PTA explanation rest on unproved model assumptions, so the significance of the review is conditional on those assumptions being made rigorous.","major_comments":[{"comment":"The claimed absolute bound νmax < THz is load-bearing for the high-frequency sections and for the abstract's aHz–THz range, but its derivation is not shown. The text states that Eq. (3.51) follows by inserting Eq. (3.50) into Eq. (3.46) and imposing Eq. (3.47). Equation (3.50) is introduced 'in a suggestive form' with no derivation, and no numerical estimate of γ is reported despite the statement that it can be obtained by integrating the mode functions. Under the paper's own definition n = |β|^2, Eq. (3.49) as printed, |n|^2/(1+|n|^2) = e^{-γx}, is dimensionally inconsistent; even the natural correction n/(1+n)=e^{-γx} is not reproduced by Eq. (3.50) for large x unless nT = 3 and γ = 1. The integrated bound and the numerical coefficient 0.165 in Eq. (3.51) therefore depend on an unverified interpolating form. The qualitative statement that νmax should lie below the THz may survive, but the 'absolute' character and the distorted-thermal spectrum are not established. Please derive Eq. (3.50) from the mode evolution, including the value of γ, or replace the bound with a general argument that does not depend on this interpolation.","section":"3.2.3, Eqs. (3.49)–(3.51)"},{"comment":"The conclusion that the PTA excess can be explained by a relic signal rests entirely on the dynamical refractive-index model. The action (5.24) and the parametrization (5.27) are postulated rather than derived, and α, n∗, and N∗ are free parameters that are then adjusted to match the PTA slope and amplitude, as in Eqs. (5.40)–(5.42). This makes the explanation partly circular: the model is defined by the effect it is supposed to produce. The paper should either provide a derivation of n(a) from an underlying microscopic action or state explicitly that (5.27) is an illustrative toy model, and it should identify an independent, testable prediction—for example, the high-frequency branch in Eq. (5.38) or a specific LVK-band signature—that could distinguish this scenario from an astrophysical foreground.","section":"5.2.3, Eqs. (5.24)–(5.27)"}],"minor_comments":[{"comment":"The multiplicity notation is inconsistent: Eq. (3.43) defines n(k,τ) = |v_k(τ)|^2, while Eq. (3.80) defines n(k,τ) = |β_k(τ)|^2. Please define v_k and β_k explicitly and use a single symbol for the averaged pair multiplicity.","section":"3.2.2"},{"comment":"If n denotes the averaged multiplicity, the left-hand side of Eq. (3.49) should be n/(1+n), not |n|^2/(1+|n|^2). This typo contributes to the mismatch with Eq. (3.50) noted in the major comments.","section":"Eq. (3.49)"},{"comment":"The statement that Ωgw ∝ ℏ^2 proves a 'truly quantum mechanical origin' is not compelling: a classical stochastic background with non-vacuum initial conditions would have an ℏ-independent energy density. The ℏ dependence only reflects the vacuum normalization of the produced state, so the sentence should be rephrased.","section":"3.2.2, after Eq. (3.44)"},{"comment":"The manuscript contains numerous typographical errors and missing words (e.g., 'paradign', 'inlationary', 'f aster', 'te literature', 'infirm or confirm'). A careful proofread would improve readability.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my take on Giovannini's topical review. It does what it says: it lays out the case that relic graviton backgrounds, with spectra across many decades, are the only direct handle on the post-inflationary expansion history before BBN. The review is a synthesis of his own program, not a new derivation, so don't go looking for new results. But it's a coherent and mostly careful presentation of the standard machinery: spectral energy density, e-fold shifts from modified expansion rates, slopes, PTA and LVK constraints. The PTA section is genuinely useful: he shows that modified post-inflationary expansion can't easily explain the nHz excess, and then explores a dynamical refractive index during inflation as an alternative. That's the most interesting part of the paper, even if the model is his own and fitted to the PTA amplitude.\n\nThe soft spots are real. The claimed absolute bound νmax < THz rests on Eq. (3.50), an interpolating form for the graviton multiplicity that is introduced without derivation. The stress-test note is right: as written, the interpolation does not reproduce the large-x behaviour of Eq. (3.49) unless you pick special values of γ and nT. So the bound is not on the same footing as the rest of the review's derivations. It may well be true, but the paper doesn't prove it. I'd treat it as a plausible estimate, not an absolute result. The refractive-index model likewise has parameters fitted to the PTA signal, so its 'explanation' of the excess is partly circular.\n\nThe good news: the central diagnostic idea doesn't depend on that bound. The sensitivity of the spectrum to the expansion history, and the discussion of e-fold indeterminacy, stand on standard theory. The review also has a nice discussion of single-graviton limits and quantumness, though that's brief.\n\nWho's this for? Someone wanting a one-stop account of how relic gravitons could map the pre-BBN expansion history, with the relevant formulas and caveats. It's a review, so novelty is low, but the synthesis is competent. I'd send it to peer review; the referee should push for a derivation or at least a numerical check of Eq. (3.50), and a flag that the 'absolute' bound is conditional.\n\nRecommendation: accept after moderate revision, with the interpolation issue addressed.","headline":"A competent topical review whose central diagnostic thesis holds up, but whose 'absolute' THz bound rests on an unproved interpolation and should be treated as conditional.","tokens_in":61074,"tokens_out":2995,"would_cite":true,"duration_ms":27258,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.30.-w","98.80.-k"],"model":"deepseek-v4-flash","headline":"This topical review argues that the relic-graviton background, spanning the aHz to THz range, is the only direct diagnostic of the expansion history between the end of inflation and big bang nucleosynthesis.","keywords":["relic gravitons","gravitational wave background","expansion history","post-inflationary evolution","pulsar timing arrays","maximal frequency","tensor-to-scalar ratio","big bang nucleosynthesis"],"falsifier":"Two observations would settle the central claims. First, a gravitational-wave signal recorded above roughly $1\\,\\mathrm{THz}$ would directly falsify the claimed absolute bound $\\nu_{\\mathrm{max}} < \\mathrm{THz}$. Second, the mapping between expansion stages and spectral slopes could be tested by measuring the slope of $\\Omega_{\\mathrm{gw}}$ in any band: if a measured slope disagrees with $n_T = 2 - 2\\delta$ for the stage sequence inferred from the break frequencies, or if a search in the MHz–GHz range reaches the amplitude required by a slower-than-radiation stage and finds nothing, the claimed correspondence between expansion history and the relic spectrum would be broken.","tokens_in":59948,"feed_emoji":"📡","tokens_out":17440,"duration_ms":135799,"temperature":0.7,"pith_summary":"This topical review argues that the diffuse background of relic gravitational waves is the only direct way to measure the expansion rate of the Universe between the end of inflation and big bang nucleosynthesis. Its guiding claim is that the graviton spectrum, stretching from the aHz to the THz domain, carries an imprint of every stage of the plasma's evolution, so each frequency band tests a different epoch of the expansion history. If the claim is right, the long-standing assumption that the early Universe was dominated by radiation immediately after inflation becomes empirically testable in the coming decade or two rather than being taken for granted. The review also derives a claimed absolute ceiling, $\\nu_{\\mathrm{max}} < \\mathrm{THz}$, on the highest frequency of the relic spectrum, and argues that relic-graviton backgrounds are homogeneous but not stationary, which rules out standard spectral-amplitude descriptions for them.","feed_headline":"One graviton spectrum, aHz to THz, records every expansion stage","feed_subtitle":"If relic gravitons are detected, the assumed radiation-dominated era after inflation can be tested directly.","key_machinery":"The carrying object is the spectral energy density of relic gravitons, $\\Omega_{\\mathrm{gw}}(\\nu,\\tau_0) = (128\\pi^3/3)\\,\\nu^4\\, n(\\nu,\\tau_0)/(H_0^2 M_P^2)$, written in terms of the averaged multiplicity $n(\\nu,\\tau_0)$ of graviton pairs produced with opposite momenta. The argument runs on four identities: (i) the slope formula $n_T = 2 - 2\\delta + O(r_T)$, which maps the expansion-rate parameter $\\delta$ (the exponent of the scale factor in each stage; $\\delta=1$ for radiation) to the tilt of the spectrum; (ii) the e-fold shift $N_{\\mathrm{max}} = \\overline{N}_{\\mathrm{max}} + \\sum_i [(\\delta_i-1)/(2(\\delta_i+1))] \\ln \\xi_i$, showing how post-inflationary stages change the number of e-folds of inflation required to fit the current Hubble patch; (iii) the frequency ladder $\\nu_{\\mathrm{max}}, \\nu_m, \\nu_r$ built from the same ratios $\\xi_i = H_{i+1}/H_i$, which locates the spectral breaks; and (iv) the interpolating pair multiplicity $n(\\nu,\\tau_0) = \\gamma x^{n_T-3}/(e^{\\gamma x}-1)$ with $x = \\nu/\\nu_{\\mathrm{max}}$, which joins the power-law spectrum to an exponential cutoff and yields the claimed absolute bound $\\nu_{\\mathrm{max}} < \\mathrm{THz}$ when a single graviton pair is demanded at the spectral maximum.","core_discovery":"On the paper's own terms, the central discovery is an explicit correspondence between the expansion history and the relic-graviton spectrum. For a sequence of post-inflationary stages with expansion-rate parameters $\\delta_i$ ($\\delta=1$ for radiation, $\\delta>1$ faster, $\\delta<1$ slower than radiation), the spectral energy density $\\Omega_{\\mathrm{gw}}(\\nu,\\tau_0)$ has a low-frequency slope fixed by the tensor-to-scalar ratio, a high-frequency slope $n_T = 2 - 2\\delta + O(r_T)$, break frequencies $\\nu_m$ and $\\nu_r$ determined by the products of Hubble-rate ratios $\\xi_i$ over the stages, and a maximal frequency $\\nu_{\\mathrm{max}}$ related to its radiation-dominated value by the same ratios. Each observed frequency band therefore reads off the expansion rate of the epoch in which those wavelengths reentered the horizon: aHz scales constrain the inflationary rate through $r_T$, nHz scales (pulsar timing arrays) can constrain inflation and the post-inflationary timeline, and kHz–THz scales carry the humps produced by stages slower than radiation. A further quantum-mechanical claim is that single-pair production sets a model-independent upper bound $\\nu_{\\mathrm{max}} < \\mathrm{THz}$. The review applies the correspondence to current data and concludes that the nHz excess reported by pulsar timing arrays cannot come from a modified post-inflationary expansion rate within this framework, whereas a refractive index acting during inflation can reproduce it.","pith_inferences":["Because the slope map $n_T = 2 - 2\\delta$ is invertible, a future measurement of $\\Omega_{\\mathrm{gw}}$ in two or more bands could in principle reconstruct the whole sequence of post-inflationary stages, a cosmic equation-of-state tomography, rather than merely confirm or exclude one assumed timeline.","If the $\\nu_{\\mathrm{max}} < \\mathrm{THz}$ bound survives scrutiny of the heuristic multiplicity formula, it sets a hard ceiling on where to look for high-frequency gravitational effects, effectively telling experimentalists to concentrate cavity and light-particle detectors below roughly a terahertz.","The homogeneous-but-non-stationary character of the relic background suggests a concrete observational discriminator: relic gravitons should show correlations depending on the sum of observation times, which stationary astrophysical foregrounds would not, so a pulsar-timing-array-style search for such a signature could separate the relic component from the binary-merger foreground.","A testable extension of the formalism is that even a null result from space-borne interferometers in the mHz band would already exclude the fastest post-inflationary stages, since those stages would have to produce a hump whose amplitude the same formulas predict in that band."],"forward_implications":["A measured relic-graviton spectrum across any extended band turns the post-inflationary expansion rate into an observable; the slope of $\\Omega_{\\mathrm{gw}}$ in that band fixes $\\delta$ of the stage that dominated when those wavelengths reentered the horizon.","Tighter bounds on the tensor-to-scalar ratio $r_T$ at CMB scales, combined with the e-fold consistency condition, translate directly into constraints on the whole post-inflationary timeline, not just on inflation itself.","The claimed bound $\\nu_{\\mathrm{max}} < \\mathrm{THz}$ delimits the search band for high-frequency and quantum-sensing detectors such as cavities, waveguides, and small interferometers: no relic signal is expected above roughly a terahertz within the standard framework.","The non-stationarity result implies that the standard spectral-amplitude description of a stationary stochastic process is inappropriate for the relic background, since its autocorrelation depends on the sum of the two times as well as their difference.","Within the review's own calculation, the nHz excess seen by pulsar timing arrays cannot be produced by a modified post-inflationary expansion rate, because the required amplitude is orders of magnitude too small, but a dynamical refractive index acting during inflation can account for it."],"supporting_citations":[{"why":"The original suggestion that post-inflationary stages expanding at different rates modify the slopes of the relic-graviton spectral energy density for frequencies above the mHz; the seed of the review's central claim.","marker":"[45]"},{"why":"Documents the high-frequency spike that appears when the post-inflationary expansion rate is slower than radiation, the signature later used to link MHz–THz signals to the expansion history.","marker":"[46, 47]"},{"why":"The conventional-lore computation predicting a minute, quasi-flat spectral energy density in the MHz region, the baseline that modified expansion histories must beat.","marker":"[22, 23, 24]"},{"why":"Shows that the tightest limits come from the largest scales, making the tensor-to-scalar ratio a direct probe of the aHz spectral energy density.","marker":"[25]"},{"why":"CMB temperature and polarization data that bound $r_T$ and $n_s$ and anchor the estimated inflationary expansion rate used throughout the review.","marker":"[42, 43, 44]"},{"why":"Derivation of the bound on the maximal frequency from the demand that a single pair of gravitons is produced at the spectral maximum.","marker":"[100, 101]"},{"why":"The big bang nucleosynthesis bound on the integrated spectral energy density of extra relativistic species, a load-bearing constraint for the whole spectrum.","marker":"[102, 103, 104, 106]"},{"why":"The analysis establishing that the relic-graviton background is homogeneous in space but not stationary, grounding the review's claim that spectral-amplitude descriptions are inappropriate.","marker":"[76]"},{"why":"Established that tensor modes are excited in curved Friedmann–Robertson–Walker backgrounds, the physical basis of relic-graviton production.","marker":"[14, 15, 16, 17]"}],"fun_headline_variants":["Graviton spectrum records every expansion stage from aHz to THz","Relic gravitons: a frequency band for each cosmic expansion era","The relic-graviton spectrum as a timeline of cosmic expansion","Expansion history decoded in the graviton spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The strongest quantitative claim, that relic gravitons cannot exist above roughly a terahertz, rests on an assumed formula for the number of graviton pairs produced near the spectral maximum, $n(\\nu,\\tau_0) = \\gamma x^{n_T-3}/(e^{\\gamma x}-1)$, which the review presents as suggestive and does not derive, so if the true pair-production rate differs, the ceiling and the high-frequency predictions built on it would weaken.","fun_headline_variants_meta":{"raw":{"variants":["Graviton spectrum records every expansion stage from aHz to THz","Relic gravitons: a frequency band for each cosmic expansion era","The relic-graviton spectrum as a timeline of cosmic expansion","Expansion history decoded in the graviton spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1820,"prompt_tokens":1036,"completion_tokens":784,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":652,"completion_tokens_details":{"reasoning_tokens":712}},"tokens_in":652,"tokens_out":784,"duration_ms":7476,"temperature":1.0,"reasoning_tokens":712,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:36:30.413741+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Two observations would settle the central claims. First, a gravitational-wave signal recorded above roughly $1\\,\\mathrm{THz}$ would directly falsify the claimed absolute bound $\\nu_{\\mathrm{max}} < \\mathrm{THz}$. Second, the mapping between expansion stages and spectral slopes could be tested by measuring the slope of $\\Omega_{\\mathrm{gw}}$ in any band: if a measured slope disagrees with $n_T = 2 - 2\\delta$ for the stage sequence inferred from the break frequencies, or if a search in the MHz–GHz range reaches the amplitude required by a slower-than-radiation stage and finds nothing, the claimed correspondence between expansion history and the relic spectrum would be broken.","supporting_citations":[],"review_version":1}