pith. machine review for the scientific record. sign in

arxiv: 2604.20063 · v2 · submitted 2026-04-21 · 🌌 astro-ph.CO · gr-qc· hep-ph

Recognition: unknown

Purely Quadratic Non-Gaussianity from Tachyonic Instability: Primordial Black Holes and Scalar-Induced Gravitational Waves

Authors on Pith no claims yet

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

classification 🌌 astro-ph.CO gr-qchep-ph
keywords primordial black holesnon-Gaussianitytachyonic instabilityscalar-induced gravitational wavescurvature perturbationsinflationdark matter
0
0 comments X

The pith

Narrow curvature spectra exponentially suppress primordial black hole formation while keeping gravitational waves strong.

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

The paper examines primordial black hole formation when curvature perturbations follow the purely quadratic non-Gaussian form arising from tachyonic instability in multi-field inflation. Using an extended Press-Schechter framework based on the compaction function, it derives the probability distribution of the linear compaction function and shows that the black hole abundance is exponentially sensitive to both the perturbation amplitude and the correlation coefficient between the smoothed field and its radial gradient. The spectral width of the curvature power spectrum controls this coefficient: broad spectra leave it mildly negative and allow significant black hole production, while sufficiently narrow spectra drive it toward minus one and produce strong exponential suppression. This separation permits a sizable scalar-induced gravitational wave signal to survive even as black hole overproduction is avoided, with thermal inflation serving as a benchmark case.

Core claim

In the tachyonic amplification scenario with curvature perturbations given by ζ = A(φ² − ⟨φ²⟩), the PBH abundance is exponentially sensitive to the amplitude of perturbations and the correlation coefficient ρ between the smoothed field and its radial gradient. Broad spectra yield mildly negative ρ and fail to suppress PBH formation, while sufficiently narrow spectra drive ρ toward −1, resulting in exponential suppression of the PBH abundance while maintaining a sizable gravitational-wave signal.

What carries the argument

The correlation coefficient ρ between the smoothed curvature field and its radial gradient, which sets the asymptotic exponential tail of the probability distribution for the linear compaction function.

If this is right

  • For sufficiently narrow spectra, PBH abundance can be exponentially suppressed while a sizable scalar-induced gravitational wave signal remains potentially detectable by future space-based interferometers.
  • Asteroid-mass PBHs could serve as dark matter in such narrow-spectrum cases without overproduction.
  • Broad spectra, as in typical thermal inflation, make it difficult to reconcile PTA observations with PBH constraints simultaneously.
  • The PBH formation probability depends exponentially on both the perturbation amplitude and the value of ρ fixed by spectral shape.

Where Pith is reading between the lines

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

  • Future gravitational wave observations at high frequencies could indirectly constrain the allowed width of the curvature power spectrum in tachyonic models.
  • The same exponential suppression might appear in other multi-field inflationary scenarios that generate quadratic non-Gaussianity, offering a route to satisfy both PBH and GW bounds without tuning amplitudes.
  • Direct probes of the curvature spectrum shape, if available, would provide a sharp test of whether observed GW signals align with the narrow-spectrum regime required for suppression.

Load-bearing premise

Curvature perturbations follow exactly the purely quadratic non-Gaussian form ζ = A(φ² − ⟨φ²⟩) from tachyonic instability, and the extended Press-Schechter framework based on the compaction function correctly calculates PBH abundance under this non-Gaussian statistics.

What would settle it

A measurement of the width of the primordial curvature power spectrum together with an observed PBH abundance that deviates from the predicted exponential dependence on ρ would test the mechanism.

Figures

Figures reproduced from arXiv: 2604.20063 by He-Xu Zhang, Mei Huang.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Posterior distributions for the parameters of the BPL model inferred from the NANOGrav 15-year dataset. The dark [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of the predicted SIGW energy density spectra with the first 14 frequency bins of the NANOGrav 15- [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Scalar-induced gravitational-wave (SIGW) spectra and corresponding PBH mass distribution. [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
read the original abstract

We investigate primordial black hole (PBH) formation in a cosmological scenario where curvature perturbations follow purely quadratic non-Gaussianity, $\zeta = A(\phi^2-\langle\phi^2\rangle)$, arising from tachyonic instability in multi-component inflationary models. Within an extended Press-Schechter framework based on the compaction function, we derive the probability distribution of the linear compaction function and its asymptotic exponential tail, demonstrating that the PBH abundance is exponentially sensitive not only to the amplitude of perturbations but also to the correlation coefficient $\rho$ between the smoothed field and its radial gradient. We further find that, in this tachyonic amplification scenario, the spectral width of the curvature power spectrum plays a decisive role in avoiding PBH overproduction: broad spectra yield mildly negative $\rho$ and fail to suppress PBH formation, while sufficiently narrow spectra drive $\rho \to -1$, resulting in exponential suppression while maintaining a sizable gravitational-wave signal. Thermal inflation provides a useful benchmark scenario with asteroid-mass PBH dark matter and high-frequency scalar-induced gravitational waves potentially detectable by future space-based interferometers, but its typically broad spectra make it challenging to reconcile PTA observations with PBH constraints.

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.

Referee Report

2 major / 3 minor

Summary. The manuscript claims that curvature perturbations with purely quadratic non-Gaussianity ζ = A(φ² − ⟨φ²⟩) arising from tachyonic instability in multi-component inflation lead to PBH formation whose abundance is exponentially sensitive to the correlation coefficient ρ between the smoothed curvature perturbation and its radial gradient. Using an extended Press-Schechter formalism based on the compaction function, it derives the probability distribution and asymptotic tail of the linear compaction function, showing that narrow spectra drive ρ → −1 and exponentially suppress PBH overproduction while preserving a sizable scalar-induced gravitational wave signal; thermal inflation is presented as a benchmark that struggles with broad spectra.

Significance. If the central derivations hold, the work identifies spectral width as a tunable parameter that controls non-Gaussian suppression of PBH formation, offering a concrete route to reconcile asteroid-mass PBH dark matter with PTA gravitational-wave data and future interferometer signals in tachyonic multi-field models. The explicit link between quadratic non-Gaussianity, ρ, and the exponential tail constitutes a useful addition to the literature on non-Gaussian PBH constraints.

major comments (2)
  1. [Sections deriving the probability distribution and asymptotic tail (referenced in abstract)] The derivation of the ρ-dependent asymptotic exponential tail for the linear compaction function (central to the suppression claim) is presented without an explicit large-deviation calculation or the resulting functional form of the tail probability when the underlying field is quadratic in a multivariate Gaussian; this step is load-bearing for the statement that narrow spectra yield exponential suppression.
  2. [Extended Press-Schechter framework section] The application of the extended Press-Schechter formalism with the compaction function to this quadratic non-Gaussian case is used to compute PBH abundance without any cross-check against numerical simulations or alternative methods; given that EPS is already an approximation whose accuracy is known to degrade for non-Gaussian tails, the quantitative reliability of the reported suppression factors remains unverified.
minor comments (3)
  1. [Notation and definitions] Clarify the precise definition and numerical evaluation of the correlation coefficient ρ for different spectral shapes, including any assumptions on the window function used for smoothing.
  2. [Thermal inflation benchmark] The discussion of thermal inflation as a benchmark would benefit from explicit values of the spectral width parameter and the resulting ρ to make the tension with PTA/PBH constraints quantitative.
  3. [Introduction and references] Add references to prior literature on quadratic non-Gaussianity in PBH formation and on the validity of compaction-based EPS for non-Gaussian fields.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript, the positive assessment of its significance, and the constructive major comments. We address each point below and have revised the manuscript to incorporate clarifications and additional details where feasible.

read point-by-point responses
  1. Referee: [Sections deriving the probability distribution and asymptotic tail (referenced in abstract)] The derivation of the ρ-dependent asymptotic exponential tail for the linear compaction function (central to the suppression claim) is presented without an explicit large-deviation calculation or the resulting functional form of the tail probability when the underlying field is quadratic in a multivariate Gaussian; this step is load-bearing for the statement that narrow spectra yield exponential suppression.

    Authors: We thank the referee for this observation. The probability distribution and asymptotic tail in the manuscript are obtained by starting from the joint Gaussian distribution of the smoothed curvature perturbation and its radial gradient (which are correlated with coefficient ρ), noting that the linear compaction function is a quadratic form in these variables due to the purely quadratic non-Gaussianity. The tail is then extracted by evaluating the relevant integral over the super-threshold region. To make the derivation fully explicit, we have added a new appendix that performs the large-deviation calculation for this quadratic form in the multivariate Gaussian case and states the resulting functional form of the tail probability. revision: yes

  2. Referee: [Extended Press-Schechter framework section] The application of the extended Press-Schechter formalism with the compaction function to this quadratic non-Gaussian case is used to compute PBH abundance without any cross-check against numerical simulations or alternative methods; given that EPS is already an approximation whose accuracy is known to degrade for non-Gaussian tails, the quantitative reliability of the reported suppression factors remains unverified.

    Authors: We agree that the extended Press-Schechter (EPS) formalism with the compaction function is an approximation whose accuracy can degrade for strongly non-Gaussian tails, and that direct numerical simulations would provide valuable verification. In the revised manuscript we have added an expanded discussion of the known limitations of EPS for non-Gaussian PBH calculations, together with references to existing literature that compares EPS results against simulations for comparable non-Gaussian models. We also stress that the reported exponential suppression is a qualitative trend driven by the ρ dependence rather than a precise numerical prediction. Performing dedicated N-body or lattice simulations for this specific tachyonic quadratic non-Gaussianity lies beyond the scope of the present theoretical study. revision: partial

Circularity Check

0 steps flagged

No circularity: derivation of compaction tail from quadratic mapping is self-contained

full rationale

The paper derives the probability distribution and asymptotic exponential tail of the linear compaction function directly from the quadratic non-Gaussian mapping ζ = A(φ² − ⟨φ²⟩) by considering the joint statistics of the smoothed curvature perturbation and its radial gradient. The correlation coefficient ρ emerges as a calculable property of the power spectrum shape rather than a fitted parameter renamed as a prediction. The extended Press-Schechter application is presented as an approximation whose accuracy is not claimed to be proven within the paper itself, but no step reduces the central result to its inputs by construction, self-citation, or ansatz smuggling. The link between spectral width and ρ is obtained from explicit integrals over the spectrum and does not constitute a tautology.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim rests on the domain assumption of purely quadratic non-Gaussianity arising from tachyonic instability and the applicability of the compaction-function-based Press-Schechter formalism; no new entities are postulated and the amplitude A and spectral parameters function as free parameters.

free parameters (3)
  • amplitude A
    Scales the quadratic non-Gaussianity and is chosen to set perturbation strength.
  • correlation coefficient ρ
    Determines the tail of the compaction distribution and is controlled by spectral width.
  • spectral width parameter
    Tuned to drive ρ toward -1 for suppression while preserving GW amplitude.
axioms (2)
  • domain assumption Curvature perturbations obey the purely quadratic form ζ = A(φ² - ⟨φ²⟩) generated by tachyonic instability in multi-component inflation.
    Invoked in the opening sentence as the starting point for the entire analysis.
  • domain assumption The extended Press-Schechter formalism based on the compaction function accurately captures PBH formation probability under this non-Gaussianity.
    Used to derive the probability distribution and asymptotic exponential tail.

pith-pipeline@v0.9.0 · 5520 in / 1705 out tokens · 36442 ms · 2026-05-10T00:48:57.414271+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

160 extracted references · 151 canonical work pages · 11 internal anchors

  1. [1]

    The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

    G. Agazieet al.(NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, Astrophys. J. Lett.951, L8 (2023), arXiv:2306.16213 [astro-ph.HE]

  2. [2]

    Agazieet al.(NANOGrav), Astrophys

    G. Agazieet al.(NANOGrav), The NANOGrav 15 yr Data Set: Observations and Timing of 68 Millisecond Pulsars, Astrophys. J. Lett.951, L9 (2023), arXiv:2306.16217 [astro-ph.HE]

  3. [3]

    Antoniadiset al.(EPTA, InPTA:), Astron

    J. Antoniadiset al.(EPTA, InPTA:), The second data release from the European Pulsar Timing Array - III. Search for gravitational wave signals, Astron. Astrophys.678, A50 (2023), arXiv:2306.16214 [astro-ph.HE]

  4. [4]

    Antoniadiset al.(EPTA), Astron

    J. Antoniadiset al.(EPTA), The second data release from the European Pulsar Timing Array - I. The dataset and timing analysis, Astron. Astrophys.678, A48 (2023), arXiv:2306.16224 [astro-ph.HE]

  5. [5]

    Antoniadis et al

    J. Antoniadiset al.(EPTA, InPTA), The second data release from the European Pulsar Timing Array - IV. Implications for massive black holes, dark matter, and the early Universe, Astron. Astrophys.685, A94 (2024), arXiv:2306.16227 [astro-ph.CO]

  6. [6]

    D. J. Reardonet al., Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array, Astrophys. J. Lett.951, L6 (2023), arXiv:2306.16215 [astro-ph.HE]

  7. [7]

    Zicet al., Publ

    A. Zicet al., The Parkes Pulsar Timing Array third data release, Publ. Astron. Soc. Austral.40, e049 (2023), arXiv:2306.16230 [astro-ph.HE]

  8. [8]

    D. J. Reardonet al., The Gravitational-wave Background Null Hypothesis: Characterizing Noise in Millisecond Pulsar Arrival Times with the Parkes Pulsar Timing Array, Astrophys. J. Lett.951, L7 (2023), arXiv:2306.16229 [astro-ph.HE]

  9. [9]

    Searching for the nano-Hertz stochastic gravitational wave background with the Chinese Pulsar Timing Array Data Release I

    H. Xuet al., Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Res. Astron. Astrophys.23, 075024 (2023), arXiv:2306.16216 [astro-ph.HE]

  10. [10]

    Recent Gravitational Wave Observation by Pulsar Timing Arrays and Primordial Black Holes: The Importance of Non-Gaussianities,

    G. Franciolini, A. Iovino, Junior., V. Vaskonen, and H. Veermae, Recent Gravitational Wave Observation by Pulsar Timing Arrays and Primordial Black Holes: The Importance of Non-Gaussianities, Phys. Rev. Lett.131, 201401 (2023), arXiv:2306.17149 [astro-ph.CO]

  11. [11]

    S. Vagnozzi, Inflationary interpretation of the stochastic gravitational wave background signal detected by pulsar timing array experiments, JHEAp39, 81 (2023), arXiv:2306.16912 [astro-ph.CO]. 19

  12. [12]

    Franciolini, D

    G. Franciolini, D. Racco, and F. Rompineve, Footprints of the QCD Crossover on Cosmological Gravitational Waves at Pulsar Timing Arrays, Phys. Rev. Lett.132, 081001 (2024), [Erratum: Phys.Rev.Lett. 133, 189901 (2024)], arXiv:2306.17136 [astro-ph.CO]

  13. [13]

    Inomata, K

    K. Inomata, K. Kohri, and T. Terada, Detected stochastic gravitational waves and subsolar-mass primordial black holes, Phys. Rev. D109, 063506 (2024), arXiv:2306.17834 [astro-ph.CO]

  14. [14]

    Cai, X.-C

    Y.-F. Cai, X.-C. He, X.-H. Ma, S.-F. Yan, and G.-W. Yuan, Limits on scalar-induced gravitational waves from the stochastic background by pulsar timing array observations, Sci. Bull.68, 2929 (2023), arXiv:2306.17822 [gr-qc]

  15. [15]

    Wang, Z.-C

    S. Wang, Z.-C. Zhao, J.-P. Li, and Q.-H. Zhu, Implications of pulsar timing array data for scalar-induced gravitational waves and primordial black holes: Primordial non-Gaussianity fNL considered, Phys. Rev. Res.6, L012060 (2024), arXiv:2307.00572 [astro-ph.CO]

  16. [16]

    Gravitational Waves from Stochastic Scalar Fluctuations,

    R. Ebadi, S. Kumar, A. McCune, H. Tai, and L.-T. Wang, Gravitational waves from stochastic scalar fluctuations, Phys. Rev. D109, 083519 (2024), arXiv:2307.01248 [astro-ph.CO]

  17. [17]

    Gouttenoire, S

    Y. Gouttenoire, S. Trifinopoulos, G. Valogiannis, and M. Vanvlasselaer, Scrutinizing the primordial black hole inter- pretation of PTA gravitational waves and JWST early galaxies, Phys. Rev. D109, 123002 (2024), arXiv:2307.01457 [astro-ph.CO]

  18. [18]

    Implications for the non-Gaussianity of curvature perturbation from pulsar timing arrays,

    L. Liu, Z.-C. Chen, and Q.-G. Huang, Implications for the non-Gaussianity of curvature perturbation from pulsar timing arrays, Phys. Rev. D109, L061301 (2024), arXiv:2307.01102 [astro-ph.CO]

  19. [19]

    K. T. Abe and Y. Tada, Translating nano-Hertz gravitational wave background into primordial perturbations taking account of the cosmological QCD phase transition, Phys. Rev. D108, L101304 (2023), arXiv:2307.01653 [astro-ph.CO]

  20. [20]

    C. Unal, A. Papageorgiou, and I. Obata, Axion-gauge dynamics during inflation as the origin of pulsar timing array signals and primordial black holes, Phys. Lett. B856, 138873 (2024), arXiv:2307.02322 [astro-ph.CO]

  21. [21]

    Z. Yi, Q. Gao, Y. Gong, Y. Wang, and F. Zhang, Scalar induced gravitational waves in light of Pulsar Timing Array data, Sci. China Phys. Mech. Astron.66, 120404 (2023), arXiv:2307.02467 [gr-qc]

  22. [22]

    Firouzjahi and A

    H. Firouzjahi and A. Talebian, Induced gravitational waves from ultra slow-roll inflation and pulsar timing arrays obser- vations, JCAP10, 032, arXiv:2307.03164 [gr-qc]

  23. [23]

    Salvio, Supercooling in Radiative Symmetry Breaking: Theory Extensions, Gravitational Wave Detection and Primor- dial Black Holes, JCAP12, 046, arXiv:2307.04694 [hep-ph]

    A. Salvio, Supercooling in Radiative Symmetry Breaking: Theory Extensions, Gravitational Wave Detection and Primor- dial Black Holes, JCAP12, 046, arXiv:2307.04694 [hep-ph]

  24. [24]

    Z.-Q. You, Z. Yi, and Y. Wu, Constraints on primordial curvature power spectrum with pulsar timing arrays, JCAP11, 065, arXiv:2307.04419 [gr-qc]

  25. [25]

    P. Bari, N. Bartolo, G. Dom` enech, and S. Matarrese, Gravitational waves induced by scalar-tensor mixing, Phys. Rev. D 109, 023509 (2024), arXiv:2307.05404 [astro-ph.CO]

  26. [26]

    Ye and A

    G. Ye and A. Silvestri, Can the Gravitational Wave Background Feel Wiggles in Spacetime?, Astrophys. J. Lett.963, L15 (2024), arXiv:2307.05455 [astro-ph.CO]

  27. [27]

    S. A. Hosseini Mansoori, F. Felegray, A. Talebian, and M. Sami, PBHs and GWs fromT 2-inflation and NANOGrav 15-year data, JCAP08, 067, arXiv:2307.06757 [astro-ph.CO]

  28. [28]

    Cheung, C

    K. Cheung, C. J. Ouseph, and P.-Y. Tseng, NANOGrav and other PTA signals and PBH from the modified Higgs inflation, Eur. Phys. J. C84, 906 (2024), arXiv:2307.08046 [hep-ph]

  29. [29]

    Scalar-induced gravitational wave interpretation of PTA data: the role of scalar fluctuation propagation speed,

    S. Balaji, G. Dom` enech, and G. Franciolini, Scalar-induced gravitational wave interpretation of PTA data: the role of scalar fluctuation propagation speed, JCAP10, 041, arXiv:2307.08552 [gr-qc]

  30. [30]

    Jin, Z.-C

    J.-H. Jin, Z.-C. Chen, Z. Yi, Z.-Q. You, L. Liu, and Y. Wu, Confronting sound speed resonance with pulsar timing arrays, JCAP09, 016, arXiv:2307.08687 [astro-ph.CO]

  31. [31]

    Bousder, A

    M. Bousder, A. Riadsolh, A. E. Fatimy, M. E. Belkacemi, and H. Ez-Zahraouy, Implications of the NANOGrav results for primordial black holes and Hubble tension, (2023), arXiv:2307.10940 [gr-qc]

  32. [32]

    B. Das, N. Jaman, and M. Sami, Gravitational wave background from quintessential inflation and NANOGrav data, Phys. Rev. D108, 103510 (2023), arXiv:2307.12913 [gr-qc]

  33. [33]

    Zhu, Z.-C

    Q.-H. Zhu, Z.-C. Zhao, S. Wang, and X. Zhang, Unraveling the early universe’s equation of state and primordial black hole production with PTA, BBN, and CMB observations*, Chin. Phys. C48, 125105 (2024), arXiv:2307.13574 [astro-ph.CO]

  34. [34]

    Ben-Dayan, U

    I. Ben-Dayan, U. Kumar, U. Thattarampilly, and A. Verma, Probing the early Universe cosmology with NANOGrav: Possibilities and limitations, Phys. Rev. D108, 103507 (2023), arXiv:2307.15123 [astro-ph.CO]

  35. [35]

    Jiang, Y

    J.-Q. Jiang, Y. Cai, G. Ye, and Y.-S. Piao, Broken blue-tilted inflationary gravitational waves: a joint analysis of NANOGrav 15-year and BICEP/Keck 2018 data, JCAP05, 004, arXiv:2307.15547 [astro-ph.CO]

  36. [36]

    Probing the equation of state of the early Universe with pulsar timing arrays,

    L. Liu, Z.-C. Chen, and Q.-G. Huang, Probing the equation of state of the early Universe with pulsar timing arrays, JCAP11, 071, arXiv:2307.14911 [astro-ph.CO]

  37. [37]

    Yi, Z.-Q

    Z. Yi, Z.-Q. You, and Y. Wu, Model-independent reconstruction of the primordial curvature power spectrum from PTA data, JCAP01, 066, arXiv:2308.05632 [astro-ph.CO]

  38. [38]

    Inflationary interpretation of the nHz gravitational-wave background,

    L. Frosina and A. Urbano, Inflationary interpretation of the nHz gravitational-wave background, Phys. Rev. D108, 103544 (2023), arXiv:2308.06915 [astro-ph.CO]

  39. [39]

    Bhaumik, R

    N. Bhaumik, R. K. Jain, and M. Lewicki, Ultralow mass primordial black holes in the early Universe can explain the pulsar timing array signal, Phys. Rev. D108, 123532 (2023), arXiv:2308.07912 [astro-ph.CO]

  40. [40]

    Yuan, D.-S

    C. Yuan, D.-S. Meng, and Q.-G. Huang, Full analysis of the scalar-induced gravitational waves for the curvature pertur- bation with local-type non-Gaussianities, JCAP12, 036, arXiv:2308.07155 [astro-ph.CO]

  41. [41]

    M. A. Gorji, M. Sasaki, and T. Suyama, Extra-tensor-induced origin for the PTA signal: No primordial black hole production, Phys. Lett. B846, 138214 (2023), arXiv:2307.13109 [astro-ph.CO]. 20

  42. [42]

    L. Liu, Y. Wu, and Z.-C. Chen, Simultaneously probing the sound speed and equation of state of the early Universe with pulsar timing arrays, JCAP04, 011, arXiv:2310.16500 [astro-ph.CO]

  43. [43]

    Z.-C. Chen, J. Li, L. Liu, and Z. Yi, Probing the speed of scalar-induced gravitational waves with pulsar timing arrays, Phys. Rev. D109, L101302 (2024), arXiv:2401.09818 [gr-qc]

  44. [44]

    Chen and L

    Z.-C. Chen and L. Liu, Can we distinguish between adiabatic and isocurvature fluctuations with pulsar timing arrays?, Sci. China Phys. Mech. Astron.68, 250412 (2025), arXiv:2402.16781 [astro-ph.CO]

  45. [45]

    Tomita, Evolution of Irregularities in a Chaotic Early Universe, Prog

    K. Tomita, Evolution of Irregularities in a Chaotic Early Universe, Prog. Theor. Phys.54, 730 (1975)

  46. [46]

    Matarrese, O

    S. Matarrese, O. Pantano, and D. Saez, General relativistic dynamics of irrotational dust: Cosmological implications, Phys. Rev. Lett.72, 320 (1994), arXiv:astro-ph/9310036

  47. [47]

    Acquaviva, N

    V. Acquaviva, N. Bartolo, S. Matarrese, and A. Riotto, Second order cosmological perturbations from inflation, Nucl. Phys. B667, 119 (2003), arXiv:astro-ph/0209156

  48. [48]

    Mollerach, D

    S. Mollerach, D. Harari, and S. Matarrese, CMB polarization from secondary vector and tensor modes, Phys. Rev. D69, 063002 (2004), arXiv:astro-ph/0310711

  49. [49]

    K. N. Ananda, C. Clarkson, and D. Wands, The Cosmological gravitational wave background from primordial density perturbations, Phys. Rev. D75, 123518 (2007), arXiv:gr-qc/0612013

  50. [50]

    Baumann, P.J

    D. Baumann, P. J. Steinhardt, K. Takahashi, and K. Ichiki, Gravitational Wave Spectrum Induced by Primordial Scalar Perturbations, Phys. Rev. D76, 084019 (2007), arXiv:hep-th/0703290

  51. [51]

    B. J. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, New cosmological constraints on primordial black holes, Phys. Rev. D81, 104019 (2010), arXiv:0912.5297 [astro-ph.CO]

  52. [52]

    Planck Constraint on Relic Primordial Black Holes

    S. Clark, B. Dutta, Y. Gao, L. E. Strigari, and S. Watson, Planck Constraint on Relic Primordial Black Holes, Phys. Rev. D95, 083006 (2017), arXiv:1612.07738 [astro-ph.CO]

  53. [53]

    Voyager 1 $e^\pm$ Further Constrain Primordial Black Holes as Dark Matter

    M. Boudaud and M. Cirelli, Voyager 1e ± Further Constrain Primordial Black Holes as Dark Matter, Phys. Rev. Lett. 122, 041104 (2019), arXiv:1807.03075 [astro-ph.HE]

  54. [54]

    R. A. Allsmanet al.(Macho), MACHO project limits on black hole dark matter in the 1-30 solar mass range, Astrophys. J. Lett.550, L169 (2001), arXiv:astro-ph/0011506

  55. [55]

    Limits on the Macho Content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds

    P. Tisserandet al.(EROS-2), Limits on the Macho Content of the Galactic Halo from the EROS-2 Survey of the Magellanic Clouds, Astron. Astrophys.469, 387 (2007), arXiv:astro-ph/0607207

  56. [56]

    Experimental Limits on Primordial Black Hole Dark Matter from the First Two Years of Kepler Data

    K. Griest, A. M. Cieplak, and M. J. Lehner, Experimental Limits on Primordial Black Hole Dark Matter from the First 2 yr of Kepler Data, Astrophys. J.786, 158 (2014), arXiv:1307.5798 [astro-ph.CO]

  57. [57]

    Mr´ oz et al.,No massive black holes in the Milky Way halo,Nature632(2024) 749 [2403.02386]

    P. Mr´ ozet al., No massive black holes in the Milky Way halo, Nature632, 749 (2024), arXiv:2403.02386 [astro-ph.GA]

  58. [58]

    Mr´ oz et al.,Limits on Planetary-mass Primordial Black Holes from the OGLE High-cadence Survey of the Magellanic Clouds , Astrophys

    P. Mr´ ozet al., Limits on Planetary-mass Primordial Black Holes from the OGLE High-cadence Survey of the Magellanic Clouds, Astrophys. J. Lett.976, L19 (2024), arXiv:2410.06251 [astro-ph.CO]

  59. [59]

    & Wright, D

    P. Mr´ ozet al., Microlensing Optical Depth, Event Rate, and Limits on Compact Objects in Dark Matter Based on 20 Yr of OGLE Observations of the Small Magellanic Cloud, Astrophys. J. Suppl.280, 49 (2025), arXiv:2507.13794 [astro-ph.GA]

  60. [60]

    Smyth, S

    N. Smyth, S. Profumo, S. English, T. Jeltema, K. McKinnon, and P. Guhathakurta, Updated Constraints on Asteroid- Mass Primordial Black Holes as Dark Matter, Phys. Rev. D101, 063005 (2020), arXiv:1910.01285 [astro-ph.CO]

  61. [61]

    B. P. Abbottet al.(LIGO Scientific, Virgo), GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs, Phys. Rev. X9, 031040 (2019), arXiv:1811.12907 [astro-ph.HE]

  62. [62]

    GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run

    R. Abbottet al.(LIGO Scientific, Virgo), GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run, Phys. Rev. X11, 021053 (2021), arXiv:2010.14527 [gr-qc]

  63. [63]

    GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run

    R. Abbottet al.(KAGRA, VIRGO, LIGO Scientific), GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run, Phys. Rev. X13, 041039 (2023), arXiv:2111.03606 [gr-qc]

  64. [64]

    A. G. Abacet al.(LIGO Scientific, VIRGO, KAGRA), GWTC-4.0: Updating the Gravitational-Wave Transient Catalog with Observations from the First Part of the Fourth LIGO-Virgo-KAGRA Observing Run, (2025), arXiv:2508.18082 [gr-qc]

  65. [65]

    M. A. Monroy-Rodr´ ıguez and C. Allen, The end of the MACHO era- revisited: new limits on MACHO masses from halo wide binaries, Astrophys. J.790, 159 (2014), arXiv:1406.5169 [astro-ph.GA]

  66. [66]

    T. D. Brandt, Constraints on MACHO Dark Matter from Compact Stellar Systems in Ultra-Faint Dwarf Galaxies, Astrophys. J. Lett.824, L31 (2016), arXiv:1605.03665 [astro-ph.GA]

  67. [67]

    Cosmic microwave background limits on accreting primordial black holes

    Y. Ali-Ha¨ ımoud and M. Kamionkowski, Cosmic microwave background limits on accreting primordial black holes, Phys. Rev. D95, 043534 (2017), arXiv:1612.05644 [astro-ph.CO]

  68. [68]

    CMB bounds on disk-accreting massive Primordial Black Holes

    V. Poulin, P. D. Serpico, F. Calore, S. Clesse, and K. Kohri, CMB bounds on disk-accreting massive primordial black holes, Phys. Rev. D96, 083524 (2017), arXiv:1707.04206 [astro-ph.CO]

  69. [69]

    P. D. Serpico, V. Poulin, D. Inman, and K. Kohri, Cosmic microwave background bounds on primordial black holes including dark matter halo accretion, Phys. Rev. Res.2, 023204 (2020), arXiv:2002.10771 [astro-ph.CO]

  70. [70]

    Carr and F

    B. Carr and F. Kuhnel, Primordial Black Holes as Dark Matter: Recent Developments, Ann. Rev. Nucl. Part. Sci.70, 355 (2020), arXiv:2006.02838 [astro-ph.CO]

  71. [71]

    A. M. Green and B. J. Kavanagh, Primordial Black Holes as a dark matter candidate, J. Phys. G48, 043001 (2021), arXiv:2007.10722 [astro-ph.CO]

  72. [72]

    B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, Constraints on primordial black holes, Rept. Prog. Phys.84, 116902 (2021), arXiv:2002.12778 [astro-ph.CO]

  73. [73]

    Laser Interferometer Space Antenna

    P. Amaro-Seoaneet al.(LISA), Laser Interferometer Space Antenna, (2017), arXiv:1702.00786 [astro-ph.IM]

  74. [74]

    Babak, M

    S. Babak, A. Petiteau, and M. Hewitson, LISA Sensitivity and SNR Calculations, (2021), arXiv:2108.01167 [astro-ph.IM]. 21

  75. [75]

    Z. Luo, Z. Guo, G. Jin, Y. Wu, and W. Hu, A brief analysis to Taiji: Science and technology, Results Phys.16, 102918 (2020)

  76. [76]

    Luo et al.,Progress of the TianQin project,Class

    J. Luoet al., Progress of the TianQin project, Class. Quant. Grav.42, 173001 (2025), arXiv:2502.11328 [gr-qc]

  77. [77]

    Fundamental Physics and Cosmology with TianQin,

    J. Luoet al., Fundamental Physics and Cosmology with TianQin, (2025), arXiv:2502.20138 [gr-qc]

  78. [78]

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

    S. Kawamuraet al., The Japanese space gravitational wave antenna DECIGO, Class. Quant. Grav.23, S125 (2006)

  79. [79]

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

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

  80. [80]

    Crowder and N

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

Showing first 80 references.