Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

Probing Gravitational Dark Matter with Ultra-high Frequency Gravitational Waves

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read The paper derives a formula tying the thermal gravitational-wave background at ultra-high frequencies to the mass and spin of a purely gravitational dark-matter particle, so a future detector near $10^{11}$ Hz could measure or constrain…

desk verdict A correct, modest paper that re-derives a known relation between thermal gravitational DM and the GW background; the reheating concern in the stress-test note is overstated. read the letter →

arxiv 2412.21137 v2 pith:NYWPWW7R submitted 2024-12-30 hep-ph astro-ph.CO

Yong Xu This is my paper · ORCID
classification hep-phastro-ph.CO PACS 95.35.+d04.30.-w
keywords gravitationaldarkmatterthermalwavesultra-highfrequencystochasticwavebackgroundproductioncogenesisreheatingtemperaturespindependence
topics Dark Matter
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

The paper argues that a dark-matter species whose only interactions are gravitational has an unavoidable observable companion: gravitational waves produced by the same hot Standard Model plasma that creates the dark matter. Because both abundances are fixed by the reheating temperature, their present-day quantities are locked together. The paper derives a closed formula, Eq. (8), expressing the thermal gravitational-wave amplitude at frequencies around $10^{11}$ Hz in terms of the dark-matter relic density, the dark-matter mass, and a spin-dependent coefficient $\alpha$. If this relation is correct, a future experiment measuring the stochastic gravitational-wave background near $10^{11}$ Hz would measure, or bound, the mass and spin of pure gravitational dark matter, and a null result at the predicted level would exclude part of the parameter space.

What carries the argument

The load-bearing object is Eq. (8), a ratio-symmetric rewriting of two Boltzmann-solved abundances: the gravitationally produced dark-matter density $\Omega_{\rm DM}h^2$ and the graviton-sourced gravitational-wave density $\Omega_{\rm GW}h^2$, both functions of the reheating temperature $T_{\rm rh}$. Dividing one by the other eliminates $T_{\rm rh}$, leaving $\Omega_{\rm GW}h^2 \propto (\Omega_{\rm DM}h^2)^{1/3}\alpha^{-1/3}m_{\rm DM}^{-1/3}f^{3}\hat{\eta}(f)$. The parameter $\alpha$ carries the spin dependence of the dark-matter production cross section, and $\hat{\eta}(f)$ carries the spectral shape and Boltzmann suppression of the gravitational-wave source; the peak frequency is inherited from the CMB temperature, around $100$ GHz. This identity is what converts a future gravitational-wave amplitude measurement into a statement about gravitational dark matter's mass and spin.

What would settle it

A future ultra-high-frequency gravitational-wave observatory with sensitivity near $\Omega_{\rm GW}h^2\sim10^{-10}$ at $f\sim10^{11}$ Hz that sees no stochastic background would falsify the predicted signal for pure gravitational dark matter with $m_{\rm DM}\lesssim10^6$ GeV under the stated thermal history; alternatively, a detected background whose amplitude or spectral shape disagrees with Eq. (8) for every allowed mass and spin would rule out the cogenesis relation.

Watch

Extended reading notes

Core claim

The central discovery is that gravitational dark matter produced by annihilations of Standard Model particles in the early thermal plasma and the stochastic gravitational-wave background emitted by the same plasma are two outputs of a single cogenesis process. The paper states the connection as $\Omega_{\rm GW}h^2 \simeq 8.6\times10^{-11}\,(\Omega_{\rm DM}h^2/0.12)^{1/3}(\alpha/3\times10^{-3})^{-1/3}(m_{\rm DM}/10^9\,{\rm GeV})^{-1/3}(f/10^{11}\,{\rm Hz})^{3}\hat{\eta}(f)$, where $\alpha$ is $1.9\times10^{-4}$, $1.1\times10^{-3}$, or $2.3\times10^{-3}$ for spin $0$, $1/2$, or $1$, and $\hat{\eta}(f)$ encodes the production sources and Boltzmann suppression. Because the reheating temperature cancels between the two known abundance formulas, the amplitude at a fixed frequency is fixed once the dark-matter relic abundance, mass, and spin are specified. The consequence is that ultra-high-frequency gravitational-wave experiments around $10^{11}$ Hz can act as a probe of the mass and spin of pure gravitational dark matter.

Load-bearing premise

The relation stands on the assumption that the early Universe was radiation-dominated from the reheating temperature down, with any preceding reheating epoch changing the dark-matter and gravitational-wave yields by only an order-one factor.

Editorial extensions

If this is right

  • If Eq. (8) is correct, the predicted thermal gravitational-wave amplitude at $10^{11}$ Hz depends on the dark-matter mass through $m_{\rm DM}^{-1/3}$, so lighter pure gravitational dark matter gives a stronger signal and heavier dark matter a weaker one.
  • A future null detection at the level $\Omega_{\rm GW}h^2 \gtrsim 10^{-10}$ near $10^{11}$ Hz would exclude pure gravitational dark matter with $m_{\rm DM}\lesssim10^6$ GeV, assuming it forms all of the dark matter.
  • With sufficient resolution, the spread in $\alpha$ across spins $0$, $1/2$, and $1$ changes the predicted amplitude enough that a measured spectrum could indicate the dark-matter spin.
  • The spectral peak sits near $100$ GHz, inherited from the cosmic microwave background temperature, so ultra-high-frequency detectors in that band are the relevant probes rather than lower-frequency interferometers.
  • The cosmic microwave background tensor-to-scalar bound, used as $T_{\rm rh}\lesssim5.5\times10^{15}$ GeV, translates through Eq. (6) into a lower bound $m_{\rm DM}\gtrsim3.2\times10^4$ GeV for the scenario.

Reading between the lines

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

  • The same logic likely applies to any feebly interacting particle whose relic abundance is set by Planck-suppressed annihilations of Standard Model plasma: its abundance and the thermal gravitational-wave yield are tied by the same cancellation of reheating temperature, so Eq. (8) could be generalized to other invisible sectors.
  • If a future experiment detects a background consistent with Eq. (8) but direct searches exclude gravitational dark matter in the implied mass window, the tension would point toward a non-thermal production component or a modified expansion history rather than disproving the cogenesis picture.
  • The paper's sensitivity to the reheating phase suggests a precision measurement of the gravitational-wave amplitude could be inverted to constrain the ratio of maximum to reheating temperature, effectively probing the duration of reheating.
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

2 major / 4 minor

Summary. The paper considers dark matter (DM) that interacts with the Standard Model only gravitationally and is produced by annihilations of SM plasma particles. It combines two existing results: the freeze-in DM relic abundance, Eq. (6), and the thermal graviton/GW spectrum from the SM plasma, Eq. (7). Eliminating the reheating temperature Trh between them yields Eq. (8), which relates the present-day gravitational-wave abundance at ultra-high frequencies to the DM mass and spin. The paper then plots the predicted Ω_GW h^2 for several masses and spins and argues that future ultra-high-frequency GW experiments could probe or exclude this scenario. The central algebraic step is correct under the stated assumption of instantaneous reheating, but the paper's treatment of non-instantaneous reheating is an assertion rather than a derivation.

Significance. If Eq. (8) is valid, it is a compact and falsifiable relation: for a given DM spin, the thermal GW amplitude at f ~ 10^11 Hz is determined by the observed DM abundance and m_DM, with no additional free parameters. The paper makes good use of established, published rates, and the algebra leading to Eq. (8) is transparent and checkable. The result is potentially valuable because it converts a difficult DM-production calculation into a concrete GW target. The main caveats are that the size and mass-dependence of reheating corrections are not quantified, and the connection to actual detector sensitivities remains qualitative. With those items addressed, the paper would be a useful contribution.

major comments (2)
  1. [Cogenesis, paragraph after Eq. (8)] The statement that a realistic reheating phase changes Eq. (8) only by an O(1) factor is not demonstrated, and the formulas quoted in the same paragraph indicate that the correction is mass-dependent. Using the stated Tmax formula, Eq. (6), and the BICEP/Keck bound H_inf ≲ 5 × 10^13 GeV, one obtains log(Tmax/Trh) ≈ 1.3 for m_DM = 10^9 GeV and ≈ 2.4 for m_DM = 10^12 GeV (s = 0). If, as the text says, Eq. (7) receives a logarithmic correction in Tmax/Trh while Eq. (6) receives a factor of about two, then the coefficient of Eq. (8) becomes roughly log(Tmax/Trh)/2^(1/3), differing by a factor of about two across the plotted mass range. This is not a single O(1) factor; it also changes the effective m_DM scaling and shifts the mass thresholds quoted in the Results. Please provide the calculation or a quantitative bound, or explicitly restrict the claim to instantaneous reheating.
  2. [Results, paragraph after Eq. (6)] The claimed lower bound "m_DM ≳ 3.2 × 10^4 GeV" does not follow from Eq. (6) with the quoted value α = 1.9 × 10^-4 and the stated bound Trh ≲ 5.5 × 10^15 GeV. Substituting Trh = 5.5 × 10^15 GeV into Eq. (6) gives m_DM ≈ 9.5 × 10^4 GeV for s = 0; to obtain 3.2 × 10^4 GeV one would need Trh ≈ 7.9 × 10^15 GeV. Please correct the number or specify the spin and α used.
minor comments (4)
  1. [Eqs. (7)–(8)] The function η̂(f) is introduced but never explicitly defined; please specify its normalization, the SM degrees of freedom entering it, and the form of the Boltzmann suppression, so that Eq. (8) can be reproduced independently.
  2. [Results, experimental reach] The statement that "a null result for Ω_GW h^2 ≳ O(10^-10)" would exclude part of the parameter space should be phrased as an upper limit below O(10^-10). It would also strengthen the paper to state which proposed ultra-high-frequency experiments from Refs. [42–44] reach the required strain or Ω_GW sensitivity.
  3. [Introduction] In the first paragraph, "the later depends on the mass and spin" should be "the latter depends on the mass and spin."
  4. [Results, Fig. 2 caption and text] The lower panel states that GW measurements "with enough resolution" could provide spin information, but the required amplitude resolution is not quantified; a brief quantitative statement (e.g., the fractional separation between spin curves) would make the claim more concrete.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: Eq. (8) is a direct algebraic elimination of Trh between two independently sourced formulas; the only self-citations are auxiliary and not load-bearing.

full rationale

The paper's central result, Eq. (8), is obtained by taking the DM abundance formula (6), solving it for Trh, and substituting that expression into the thermal GW spectrum formula (7). This is an explicit algebraic elimination: Omega_GW h^2 = 8.6e-11 (Trh/1e14 GeV)(f/1e11 Hz)^3 eta_hat(f) with Trh = (Omega_DM h^2 / 0.12)^(1/3) (alpha/3e-3)^(-1/3) (mDM/1e9 GeV)^(-1/3) 1e14 GeV. The result therefore does not presuppose the GW amplitude; it is a derived consequence of two published production-rate formulas, Eq. (6) from Ref. [9] and Eq. (7) from Refs. [10-12], which are external to this paper. Setting Omega_DM h^2 = 0.12 is a boundary condition (the observed relic density), not a fitted parameter of the GW prediction. The spin dependence enters through alpha, which is also taken from Ref. [9], not inferred from the GW spectrum. The only self-citations are Ref. [20], by the same author with N. Bernal, used for the peak-frequency value around 80 GHz and for the claimed logarithmic correction log(Tmax/Trh) ~ O(1), and Refs. [29,30] used for inflaton-mediated and reheating corrections. These citations are not used to fit the output of this paper, and Eq. (8) stands independently in the radiation-dominated limit T <= Trh that the paper explicitly assumes. The reheating-phase assertion that corrections are O(1) is under-supported and may be a numerical/correctness risk—the paper's own Tmax formula can give log(Tmax/Trh) significantly larger than O(1) for plausible parameters—but this is not circularity: Eq. (8) is transparently derived, and the contested O(1) claim is an auxiliary robustness statement rather than the derivation itself. Thus the derivation chain is self-contained, with at most one minor self-citation that does not bear the weight of the central claim.

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

No new free parameters are fitted; alpha is computed from spin-dependent matrix elements in Ref. [9]. The central relation relies on standard gravitational freeze-in formulas and on the assumed thermal history. No new particles or forces are introduced.

assumptions (7)
  • domain assumption Dark matter has no interactions except gravity (minimal scenario).
    Defines the scenario; Eq. (1) action has only gravitational coupling between LSM and LDM.
  • domain assumption Dark matter is produced solely by SM thermal plasma annihilations via graviton exchange.
    The paper assumes no inflaton or other source contributes to the relic abundance; Section 'Cogenesis'.
  • domain assumption The Universe is radiation-dominated for T <= Trh, with reheating corrections O(1).
    Stated in Section 'Cogenesis'; allows Eq. (8) to be derived without Tmax dependence.
  • domain assumption Gravitons do not thermalize and propagate freely to form a stochastic GW background.
    Assumed after production; used to compute present-day Omega_GW.
  • domain assumption The dark matter relic abundance is fixed to the observed value Omega_DM h^2 = 0.12.
    Used in Eq. (8) to eliminate Trh; this is an observational input, not a prediction.
  • domain assumption Cited formulas Eq. (6) from Ref. [9] and Eq. (7) from Ref. [20] are correct.
    The derivation does not recompute these rates; it takes them as given.
  • domain assumption No other significant GW sources contribute in the 10^11 Hz band, so an observed background can be attributed to this mechanism.
    Implied by the interpretation of future measurements in the Results section.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Probing Gravitational Dark Matter with Ultra-high Frequency Gravitational Waves." pith.science (2026). https://pith.science/paper/NYWPWW7R

@misc{pith2026241221137,
  author       = {Pith},
  title        = {Pith review of: Probing Gravitational Dark Matter with Ultra-high Frequency Gravitational Waves},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NYWPWW7R}},
  note         = {Machine review of arXiv:2412.21137}
}
read the original abstract

The evidence for the existence of dark matter (DM) is compelling, yet its nature remains elusive. A minimal scenario involves DM interacting solely through gravity. However, the detection would be extremely challenging. In the early Universe, such DM can be unavoidably generated via annihilation of particles in the standard model (SM) thermal plasma. It is known that the SM thermal plasma also produces gravitational waves (GWs). In this study, we establish a simple connection between the amplitude of thermal GWs and the properties of pure gravitational DM. Notably, future GW experiments in the ultra-high frequency regime have the potential to shed light on the mass and spin of pure gravitational DM.

Figures

Figures reproduced from arXiv: 2412.21137 by the authors.

Figure 1
Figure 1. FIG. 1. Feynman diagrams for cogenesis of DM (left) and [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. GW spectrum as function of frequency [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Dark Matter Ultraviolet Freeze-in in General Reheating Scenarios

    hep-ph 2025-01 accept novelty 6.0 of 10

    The paper derives analytic dark matter freeze-in yields for arbitrary power-law reheating histories and maps the gravitational production parameter space.

Reference graph

Works this paper leans on

44 extracted references · 38 linked inside Pith · cited by 1 Pith paper

  1. [1]

    Bertone and D

    G. Bertone and D. Hooper, Rev. Mod. Phys.90, 045002 (2018), arXiv:1605.04909 [astro-ph.CO]

    arXiv 2018

  2. [2]

    Cirelli, A

    M. Cirelli, A. Strumia, and J. Zupan, (2024), arXiv:2406.01705 [hep-ph]

    arXiv 2024

  3. [3]

    Balazs, T

    C. Balazs, T. Bringmann, F. Kahlhoefer, and M. White, (2024), arXiv:2411.05062 [astro-ph.CO]

    arXiv 2024

  4. [4]

    R. H. Cyburt, B. D. Fields, K. A. Olive, and T.-H. Yeh, Rev. Mod. Phys. 88, 015004 (2016), arXiv:1505.01076 [astro-ph.CO]

    arXiv 2016

  5. [5]

    Garny, M

    M. Garny, M. Sandora, and M. S. Sloth, Phys. Rev. Lett. 116, 101302 (2016), arXiv:1511.03278 [hep-ph]

    arXiv 2016

  6. [6]

    Tang and Y.-L

    Y. Tang and Y.-L. Wu, Phys. Lett. B758, 402 (2016), arXiv:1604.04701 [hep-ph]

    arXiv 2016

  7. [7]

    Tang and Y.-L

    Y. Tang and Y.-L. Wu, Phys. Lett. B774, 676 (2017), arXiv:1708.05138 [hep-ph]

    arXiv 2017

  8. [8]

    Garny, A

    M. Garny, A. Palessandro, M. Sandora, and M. S. Sloth, JCAP 02, 027 (2018), arXiv:1709.09688 [hep-ph]

    arXiv 2018

Show all 44 references
  1. [9]

    Bernal, M

    N. Bernal, M. Dutra, Y. Mambrini, K. Olive, M. Peloso, and M. Pierre, Phys. Rev. D 97, 115020 (2018), arXiv:1803.01866 [hep-ph]

    2018 arXiv

  2. [10]

    Ghiglieri and M

    J. Ghiglieri and M. Laine, JCAP 07, 022 (2015), arXiv:1504.02569 [hep-ph]

    2015 arXiv

  3. [11]

    Ghiglieri, G

    J. Ghiglieri, G. Jackson, M. Laine, and Y. Zhu, JHEP 07, 092 (2020), arXiv:2004.11392 [hep-ph]

    2020 arXiv

  4. [12]

    Ringwald, J

    A. Ringwald, J. Sch¨ utte-Engel, and C. Tamarit, JCAP 03, 054 (2021), arXiv:2011.04731 [hep-ph]

    2021 arXiv

  5. [13]

    Klose, M

    P. Klose, M. Laine, and S. Procacci, JCAP 05, 021 (2022), arXiv:2201.02317 [hep-ph]

    2022 arXiv

  6. [14]

    Klose, M

    P. Klose, M. Laine, and S. Procacci, JCAP 12, 020 (2022), arXiv:2210.11710 [hep-ph]

    2022 arXiv

  7. [15]

    Ringwald and C

    A. Ringwald and C. Tamarit, Phys. Rev. D106, 063027 (2022), arXiv:2203.00621 [hep-ph]

    2022 arXiv

  8. [16]

    Ghiglieri, J

    J. Ghiglieri, J. Sch¨ utte-Engel, and E. Speranza, Phys. Rev. D 109, 023538 (2024), arXiv:2211.16513 [hep-ph]

    2024 arXiv

  9. [17]

    F. Muia, F. Quevedo, A. Schachner, and G. Villa, JCAP 09, 006 (2023), arXiv:2303.01548 [hep-ph]

    2023 arXiv

  10. [18]

    Drewes, Y

    M. Drewes, Y. Georis, J. Klaric, and P. Klose, JCAP 06, 073 (2024), arXiv:2312.13855 [hep-ph]

    2024 arXiv

  11. [19]

    Ghiglieri, M

    J. Ghiglieri, M. Laine, J. Sch¨ utte-Engel, and E. Sper- anza, JCAP 04, 062 (2024), arXiv:2401.08766 [hep-ph]

    2024 arXiv

  12. [20]

    Bernal and Y

    N. Bernal and Y. Xu, JHEP 01, 137 (2025), arXiv:2410.21385 [hep-ph]

    2025 arXiv

  13. [21]

    S. M. Boucenna and S. Morisi, Front. in Phys.1, 33 (2014), arXiv:1310.1904 [hep-ph]

    2014 arXiv

  14. [22]

    Maleknejad and J

    A. Maleknejad and J. Kopp, (2024), arXiv:2405.09723 [hep-ph]

    2024 arXiv

  15. [23]

    Maleknejad and J

    A. Maleknejad and J. Kopp, JHEP 01, 023 (2025), arXiv:2406.01534 [hep-th]

    2025 arXiv

  16. [24]

    J. F. Donoghue, Phys. Rev. D50, 3874 (1994), arXiv:gr- qc/9405057

    1994

  17. [25]

    G. F. Giudice, E. W. Kolb, and A. Riotto, Phys. Rev. D 64, 023508 (2001), arXiv:hep-ph/0005123

    2001 arXiv

  18. [26]

    A. A. Starobinsky, Phys. Lett. B91, 99 (1980)

    1980

  19. [27]

    Drees and Y

    M. Drees and Y. Xu, JCAP 09, 012 (2021), arXiv:2104.03977 [hep-ph]

    2021 arXiv

  20. [28]

    Drees and Y

    M. Drees and Y. Xu, JCAP 12, 005 (2022), arXiv:2209.07545 [astro-ph.CO]

    2022 arXiv

  21. [29]

    Xu, JHEP10, 174 (2024), arXiv:2407.03256 [hep-ph]

    Y. Xu, JHEP10, 174 (2024), arXiv:2407.03256 [hep-ph]

    2024 arXiv

  22. [30]

    Bernal, J

    N. Bernal, J. Harz, M. A. Mojahed, and Y. Xu, Phys. Rev. D 111, 043517 (2025), arXiv:2406.19447 [hep-ph]. 5

    2025 arXiv

  23. [31]

    D. J. H. Chung, E. W. Kolb, and A. Riotto, Phys. Rev. D 59, 023501 (1998), arXiv:hep-ph/9802238

    1998 arXiv

  24. [32]

    D. J. H. Chung, E. W. Kolb, and A. Riotto, Phys. Rev. D 60, 063504 (1999), arXiv:hep-ph/9809453

    1999 arXiv

  25. [33]

    Y. Ema, K. Nakayama, and Y. Tang, JHEP 09, 135 (2018), arXiv:1804.07471 [hep-ph]

    2018 arXiv

  26. [34]

    Y. Ema, K. Nakayama, and Y. Tang, JHEP 07, 060 (2019), arXiv:1903.10973 [hep-ph]

    2019 arXiv

  27. [35]

    Mambrini and K

    Y. Mambrini and K. A. Olive, Phys. Rev. D103, 115009 (2021), arXiv:2102.06214 [hep-ph]

    2021 arXiv

  28. [36]

    Clery, Y

    S. Clery, Y. Mambrini, K. A. Olive, and S. Verner, Phys. Rev. D 105, 075005 (2022), arXiv:2112.15214 [hep-ph]

    2022 arXiv

  29. [37]

    E. W. Kolb and A. J. Long, Rev. Mod. Phys.96, 045005 (2024), arXiv:2312.09042 [astro-ph.CO]

    2024 arXiv

  30. [38]

    Y. Ema, R. Jinno, K. Mukaida, and K. Nakayama, JCAP 05, 038 (2015), arXiv:1502.02475 [hep-ph]

    2015 arXiv

  31. [39]

    Y. Ema, R. Jinno, and K. Nakayama, JCAP09, 015 (2020), arXiv:2006.09972 [astro-ph.CO]

    2020 arXiv

  32. [40]

    G. Choi, W. Ke, and K. A. Olive, Phys. Rev. D109, 083516 (2024), arXiv:2402.04310 [hep-ph]

    2024 arXiv

  33. [41]

    P. A. R. Adeet al.(BICEP, Keck), Phys. Rev. Lett.127, 151301 (2021), arXiv:2110.00483 [astro-ph.CO]

    2021

  34. [42]

    A. M. Cruise, Class. Quant. Grav.29, 095003 (2012)

    2012

  35. [43]

    Aggarwal et al

    N. Aggarwal et al. , Living Rev. Rel. 24, 4 (2021), arXiv:2011.12414 [gr-qc]

    2021 arXiv

  36. [44]

    Aggarwalet al., (2025), arXiv:2501.11723 [gr-qc]

    N. Aggarwalet al., (2025), arXiv:2501.11723 [gr-qc]

    2025

Pith tools

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