Pith. sign in

REVIEW 3 major objections 4 minor 12 references

Constraining primordial curvature perturbations with present and future GW detectors

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper argues that gravitational-wave detectors can probe small-scale primordial curvature perturbations far beyond CMB reach, and that a null LISA measurement would cap the amplitude near $10^{-4}$.

desk verdict A faithful but non-verifiable proceedings summary; useful for orientation, not for citing as a source of forecasts. read the letter →

arxiv 2505.23521 v1 pith:ARPK2L5K submitted 2025-05-29 astro-ph.CO astro-ph.HEgr-qchep-phhep-th

Mauro Pieroni This is my paper
classification astro-ph.COastro-ph.HEgr-qchep-phhep-th
keywords scalar-inducedgravitationalwavesprimordialcurvatureperturbationspulsartimingarraysLISAstochasticwavebackgroundlognormalpeakblackholessupermassiveholebinaries
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Primordial curvature perturbations are well measured on CMB scales (around $k\sim10^{-4}$--$10^{-1}\,\mathrm{Mpc}^{-1}$), but on smaller scales they are almost unconstrained; if early-universe dynamics enhance them, those perturbations source a gravitational-wave background at second order. This proceeding argues that gravitational-wave detectors are therefore a unique probe of $P_\zeta(k)$ on small scales: current pulsar-timing-array data can be read either as supermassive-black-hole binaries or as scalar-induced gravitational waves, and future arrays with more pulsars and longer observing times should distinguish the two or tighten upper limits. For LISA, the paper's central forecast is that a null measurement would cap the reconstructed amplitude $A_\zeta$ at roughly $10^{-4}$, while a strong detection would measure the peak amplitude, width, and scale to sub-percent accuracy. That matters because these instruments reach comoving scales $k\sim10^{6}$--$10^{14}\,\mathrm{Mpc}^{-1}$, far beyond the reach of the CMB and large-scale-structure surveys.

What carries the argument

The load-bearing object is the second-order integral for the tensor spectrum, $\mathcal{P}_h(k,\eta)=4\int_0^\infty dt\int_{-1}^1 ds\,[t^2(1-s^2)^2/(1+t)^2]\,I^2(k,t,s,\eta)\,\mathcal{P}_\zeta(k_u)\mathcal{P}_\zeta(k_v)$, where $u=(1+t-s)/2$, $v=(1+t+s)/2$, and $I(k,t,s,\eta)$ is the kernel built from cosmological transfer functions. This product structure makes the induced spectrum quadratic in the scalar power spectrum and gives it model-independent features: a double peak for narrow enhancements, a causal $f^3$ infrared tail, and an ultraviolet branch nearly proportional to $\mathcal{P}_\zeta^2$, which the forecasts exploit. To remain model-agnostic, the analysis folds the integral into a lognormal peak or a binned parameterization of $P_\zeta(k)$, so the stated constraints are statements about the parameters $(A,\Delta,k_*)$ of that template rather than about any specific inflationary model.

What would settle it

Run the same analysis pipeline on simulated data that contain a non-lognormal enhancement of $P_\zeta(k)$, for example a broken power law with amplitude just below $10^{-4}$, and check whether the recovered bin-by-bin limits contain the true injected amplitude; if the reconstruction is biased, the claim that LISA constrains $P_\zeta$ to $A_\zeta\lesssim10^{-4}$ fails.

Watch

Extended reading notes

Core claim

Scalar and tensor perturbations decouple at linear order but couple at second order, so an enhanced scalar power spectrum $P_\zeta(k)$ sources a stochastic background of scalar-induced gravitational waves. The paper's claim is that present and future gravitational-wave observatories can use this background to constrain $P_\zeta$ on scales that no other cosmological probe reaches. Concretely, under the lognormal parameterization $P_\zeta(k)=A/\sqrt{2\pi}\Delta \exp[-\tfrac12(\ln(k/k_*)/\Delta)^2]$, future pulsar timing arrays should recover the three peak parameters with relative errors around ten percent or place corresponding upper limits, and LISA, in the absence of a signal, would set an upper limit $A_\zeta\lesssim10^{-4}$ using a binned parameterization of the spectrum. With a sufficiently loud detection, LISA would measure the parameter triplet to better than one percent, constrain the shape of an ultra-slow-roll inflationary potential to the ten-percent level, and probe primordial non-Gaussianities down to $f_{NL}\sim10$ at percent-level sensitivity.

Load-bearing premise

The projected constraints and limits assume that the small-scale enhancement of $P_\zeta(k)$ is well described by the lognormal peak or binned template used in the analysis, and that the astrophysical foreground, modeled as a power law with injected tilt $n_T=2$, is not mis-modeled enough to shift the recovered bounds.

Editorial extensions

If this is right

  • If the current pulsar-timing-array signal is scalar-induced, a future array with $N_p=70$ pulsars would measure $A_\zeta$, $\Delta$, and $k_*$ with relative errors of order 10%, improving as $\sqrt{70/N_p}$.
  • If the current signal is instead supermassive-black-hole binaries, future pulsar-timing-array data should either detect a subdominant scalar-induced component or place projected 95--99.7% CL upper limits on $A_\zeta$ as a function of peak scale and foreground amplitude.
  • A quiet LISA would fill the gap between CMB and pulsar-timing-array scales by bounding $A_\zeta\lesssim10^{-4}$ through bin-by-bin limits on $P_\zeta$, a level currently constrained only by primordial-black-hole overproduction and big-bang nucleosynthesis.
  • A loud LISA detection would measure the lognormal parameters below 1% at SNR $\gtrsim1000$, constrain the inflationary potential shape to about 10%, and provide percent-level sensitivity to $f_{NL}\sim10$ non-Gaussianities.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Extending the paper's logic, the quoted $A_\zeta\lesssim10^{-4}$ limit is template-dependent by construction; rerunning the same pipeline with non-lognormal peak shapes, such as broken power laws or multi-field features, is required before the bound can be quoted for any specific early-universe model.
  • The paper does not combine its projected LISA bound with primordial-black-hole abundance limits, but a joint statement would be stronger than either probe alone because PBH overproduction already excludes comparable small-scale amplitudes.
  • The same bin-by-bin reconstruction, applied at decihertz-to-hertz frequencies accessible to future ground-based detectors, would extend $P_\zeta$ constraints to even smaller scales than LISA, an extension the paper does not develop.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This proceedings paper reviews the mechanism of scalar-induced gravitational waves (SIGWs) and presents forecasts for constraining the primordial curvature power spectrum P_ζ with current and future pulsar timing arrays (PTAs) and with LISA. The paper introduces the standard second-order tensor power spectrum integral and a lognormal parameterization for an enhanced P_ζ, then summarizes projected constraints from two companion papers: SKA-era PTAs would measure lognormal peak parameters to ~10% relative errors, and LISA would either measure parameters to sub-percent accuracy for SNR ≳ 1000 or set an upper limit A_ζ ≲ 10^-4. The paper concludes that GW detectors offer a unique probe of small-scale curvature perturbations.

Significance. If the quoted forecasts are correct, the paper highlights a scientifically important avenue: GW detectors could probe P_ζ on scales many orders of magnitude smaller than CMB or LSS observations, with LISA reaching amplitude limits near 10^-4. The review portion is a concise and faithful summary of the standard SIGW formalism. However, the paper's central quantitative claims are imported from two unpublished companion papers (refs. 8 and 11) that are not available to the reader, so the headline numbers cannot be independently verified from this manuscript. The presence of an internal inconsistency in the reported comoving wavenumber range further undermines confidence in the presentation. The paper's value as a standalone contribution is therefore limited at present.

major comments (3)
  1. [Section 5] The conclusion states that future PTA observations will set the most stringent upper limit on P_ζ for comoving wavenumbers k ∼ 10^11 − 10^14 Mpc^-1, but Section 1 and the standard mapping between frequency and wavenumber assign that range to LISA's mHz band, while PTAs probe k ∼ 10^6 − 10^8 Mpc^-1. This is a direct internal contradiction and is load-bearing because it misattributes the LISA-scale constraints to PTAs, obscuring which experiment probes which scales.
  2. [Sections 3 and 4] All quantitative forecast results — the SKA relative errors of order 10%, the LISA sub-percent accuracies for SNR ≳ 1000, and the A_ζ ≲ 10^-4 upper limit — are attributed to refs. 8 and 11, which are cited as unpublished works without arXiv identifiers or journal references. The text explicitly defers to these papers ('For more details, see 11' and 'for details, see 8'), so the reader cannot check the Fisher/Bayesian calculations, the treatment of the SMBHB foreground (including the injected tilt n_T = 2), the binned parameterization, or the precise definition of the quoted A_ζ limits. Because the title and conclusions rest on these numbers, the central claims are unsupported by the manuscript itself. The author should either provide a technical appendix summarizing the forecast methodology and results, or ensure the companion papers are publicly available with stable citations before making these quantitative claims.
  3. [Section 4] The text says LISA will constrain P_ζ using a binned parameterization and then quotes an upper limit 'A_ζ ≲ 10^-4', but A_ζ is defined in Eq. (2) only as the amplitude of a lognormal peak. The relation between the binned analysis and the lognormal amplitude is not defined, nor are the bin widths or the scale range of the binned reconstruction. This ambiguity makes the headline upper limit difficult to interpret or to compare with other constraints in the literature.
minor comments (4)
  1. [Section 2] In the sentence listing characteristic features of Ω_GW(f), the citation to ref. 8 appears after the colon and before the bullet list, which is grammatically awkward; consider moving the citation to the end of the introductory clause or after the list.
  2. [Section 3] The phrase 'We have performed our forecasts 11, with the open-source fastPTA 12 tool' is awkward; better phrasing would be 'We performed forecasts with the open-source fastPTA tool [12], as detailed in [11].'
  3. [Section 4] There is a typo: 'sufficienly' should be 'sufficiently'.
  4. [Section 2] The text refers to 'the near infrared part of the SIGW spectrum' in Section 3; this should likely be 'infrared' (or 'low-frequency') part, since the universal f^3 scaling is an infrared property.

Circularity Check

1 steps flagged · score 4.0 of 10

Forecast constraints are imported from co-authored companion papers; the SIGW formalism itself is standard, so the circularity is limited to self-citation reliance.

  1. self citation load bearing [Section 3 (Constraints/prospects from Pulsar Timing Arrays) and Section 4 (Prospects with LISA)]
    "We have performed our forecasts 11, with the open-source fastPTA 12 tool, and we have explored two different scenarios: ... We found 11 that for an SKA-like detector with Np = 70 pulsars, constraints improve significantly, with relative errors on the three parameters of the lognormal peak ( Aζ, ∆, k∗) of the order of 10%, scaling as p 70/Np with Np. ... In 8, a collaborative project within the LISA Cosmology Working Group (CosWG), we have thoroughly investigated the prospects for LISA to detect and characterize SIGWs. ... LISA will constrain Pζ to the level of Aζ ≲ 10−4."

    The manuscript's central quantitative claims (10% parameter errors for SKA, sub-1% measurements at high SNR, and Aζ ≲ 10−4 with LISA) are not derived from eqs. (1) and (2) in the text. Each number is attributed to refs. 11 and 8, both co-authored by the present author. The paper provides no calculation, code, or numerical recipe that would let a reader reproduce these forecasts; the load-bearing step is therefore a chain of self-citations to companion papers without stable identifiers. Because those papers are not verified in the manuscript, the quoted 'predictions' reduce to the authority of the author's own prior work rather than to any independent derivation.

full rationale

The paper is a proceedings-style summary. The underlying equations (1) and (2) are standard SIGW results from the cited literature (refs. 3 and 9), and the lognormal and binned parameterizations are stated explicitly as assumptions. The forecasts themselves, however, come exclusively from refs. 8 and 11, which are co-authored by the presenter and cited without stable identifiers. That makes the headline constraints load-bearing self-citations rather than independent, in-manuscript derivations; nonetheless, the companion papers are separate sensitivity studies using public tools (fastPTA, SIGWAY), and the physics is not equivalent to the input by construction. There is also an internal inconsistency in Section 5, where PTA is said to probe k ∼ 10^11–10^14 Mpc^−1, although Section 1 correctly assigns PTAs to k ∼ 10^6–10^8 Mpc^−1, which suggests the summary numbers may not have been carefully checked. Overall, the self-citation reliance warrants a score of 4, not higher, because the formalisms cited are externally benchmarked and the paper does not force its conclusions by definition.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The ledger is light because this is a review: its equations are standard, and its forecasts inherit the free parameters of the lognormal and power-law models from refs. 8 and 11. No new entity, force, or conserved quantity is postulated. The main modeling choices, the lognormal shape and the power-law foreground, are the ones that carry the quantitative content.

free parameters (5)
  • lognormal peak amplitude A = injected in forecasts; projected LISA null limit A_zeta <= 10^-4
    Eq. (2) parametrizes the enhancement of P_zeta with amplitude A; A is scanned and constrained in the forecasts of Sections 3 and 4.
  • lognormal peak width Delta = injected log10(Delta) = -0.3 in Scenario 2
    Width of the P_zeta peak in eq. (2); injected by hand in Section 3 and scanned in Section 4.
  • lognormal peak scale k* or f* = injected log10(f*/Hz) = -7.8 in Scenario 2
    Peak scale of the enhancement; chosen by hand for the Scenario 2 injections reported in Section 3.
  • SMBHB foreground tilt n_T = injected n_T = 2 in Scenario 2
    Tilt of the power-law astrophysical foreground in Section 3, Scenario 2; the projected upper limits depend on this assumed shape.
  • binned P_zeta bin amplitudes = not given, see ref. 8
    Section 4 uses a binned parameterization of P_zeta for the LISA null-result analysis; the bin amplitudes are free parameters, with details deferred to ref. 8.
assumptions (5)
  • standard math Scalar perturbations source tensor perturbations at second order in cosmological perturbation theory (eq. 1).
    Standard second-order perturbation theory; the paper cites Tomita 1975, Matarrese et al. 1994, and the Domenech review (refs. 1-3).
  • domain assumption Perturbation evolution during mode re-entry is governed by standard transfer functions in a minimal Lambda-CDM background with an effectively radiation-dominated epoch.
    The kernel I(k,t,s,eta) in eq. (1) encodes this evolution; the review does not state the assumed background beyond the minimal Lambda-CDM context of Section 1.
  • domain assumption The unknown small-scale enhancement of P_zeta is parameterized by a lognormal peak with three parameters (eq. 2).
    Section 2 introduces eq. (2) 'to remain agnostic about the specific model'; all quantitative forecasts in Sections 3 and 4 depend on this shape assumption.
  • domain assumption The current PTA common-spectrum process can be modeled as a power-law SMBHB background with a subdominant SIGW component in Scenario 2.
    Section 3 injects n_T = 2 and fits SMBHB + SIGW to set upper limits; the projected bounds on A_zeta depend on this foreground model.
  • domain assumption Fisher Information Matrix forecasts are a valid proxy for actual measurement uncertainties over the quoted parameter ranges.
    Section 4 reports that LISA results were 'performed with FIM formalism and validated with full Bayesian analysis (via Monte Carlo techniques)'; the validation is cited to ref. 8, not shown here.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Constraining primordial curvature perturbations with present and future GW detectors." pith.science (2026). https://pith.science/paper/ARPK2L5K

@misc{pith2026250523521,
  author       = {Pith},
  title        = {Pith review of: Constraining primordial curvature perturbations with present and future GW detectors},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ARPK2L5K}},
  note         = {Machine review of arXiv:2505.23521}
}
abstract

Primordial scalar curvature perturbations ($\zeta$), typically probed on large cosmological scales via CMB and LSS observations, can be significantly enhanced on smaller scales by various early Universe mechanisms, for instance, non-minimal inflationary models. While decoupled at linear order, scalar and tensor perturbations, i.e., Gravitational Waves (GWs), interact at second order. As a consequence, an enhanced primordial scalar power spectrum $P_\zeta(k)$ can source a sizable stochastic GW background (SGWB). In these proceedings, we briefly review the generation mechanism of such signals, typically referred to as scalar-induced GWs (SIGWs), and discuss the prospects of measuring them with present and future Pulsar Timing Arrays datasets and future GW observatories like the Laser Interferometer Space Antenna LISA.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

12 extracted references · 11 canonical work pages

  1. [1]

    K. Tomita. Non-Linear Behavior of Nonspherical Perturbations in Relativistic Gravita- tional Collapse. Prog. Theor. Phys., 54:730, 1975

    work page 1975

  2. [2]

    General relativistic dynamics of irrotational dust: Cosmological implications

    Sabino Matarrese, Ornella Pantano, and Diego Saez. General relativistic dynamics of irrotational dust: Cosmological implications. Phys. Rev. Lett., 72:320–323, 1994

    work page 1994

  3. [3]

    Scalar Induced Gravitational Waves Review

    Guillem Dom` enech. Scalar Induced Gravitational Waves Review. Int. J. Mod. Phys. D, 30(08):2130004, 2021

    work page 2021

  4. [4]

    The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background

    Agazie, Gabriella and others. The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background. Astrophys. J. Lett., 951(1):L8, 2023

    work page 2023

  5. [5]

    and others

    Antoniadis, J. and others. The second data release from the European Pulsar Timing Array: III. Search for gravitational wave signals. Astron. Astrophys., 678:A50, 2023

    work page 2023

  6. [6]

    The NANOGrav 15 yr Data Set: Search for Signals from New Physics

    Afshordi, Niayesh and others. The NANOGrav 15 yr Data Set: Search for Signals from New Physics. Astrophys. J. Lett., 951(1):L11, 2023

    work page 2023

  7. [7]

    Laser Interferometer Space Antenna

    Pau Amaro-Seoane et al. Laser Interferometer Space Antenna. 2 2017

    2017

  8. [8]

    Reconstructing Primordial Curvature Perturbations via Scalar- Induced Gravitational Waves with LISA

    Jonas El Gammal et al. Reconstructing Primordial Curvature Perturbations via Scalar- Induced Gravitational Waves with LISA. 1 2025

    work page 2025

Show all 12 references
  1. [9]

    Gravitational Waves Induced by Scalar Perturbations with a Lognormal Peak

    Shi Pi and Misao Sasaki. Gravitational Waves Induced by Scalar Perturbations with a Lognormal Peak. JCAP, 09:037, 2020

    2020

  2. [10]

    Gravitational wave astronomy with the SKA

    Gemma Janssen et al. Gravitational wave astronomy with the SKA. PoS, AASKA14:037, 2015

    2015

  3. [11]

    Prospects for Scalar-Induced Gravitational Waves detection with SKA

    Chiara Cecchini, Gabriele Franciolini, and Mauro Pieroni. Prospects for Scalar-Induced Gravitational Waves detection with SKA. 2025

    2025

  4. [12]

    Forecasting the sensitivity of pulsar timing arrays to gravitational wave backgrounds

    Stanislav Babak, Mikel Falxa, Gabriele Franciolini, and Mauro Pieroni. Forecasting the sensitivity of pulsar timing arrays to gravitational wave backgrounds. Phys. Rev. D , 110(6):063022, 2024. 4

    2024

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.