{"id":"61fe1683-9bf8-4717-91dd-ae78e7920d6b","arxiv_id":"2412.21137","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The thermal gravitational wave amplitude at around 100 GHz is set by the mass and spin of purely gravitational dark matter, so future ultra-high-frequency detectors could probe the scenario.","lead":"This paper derives a simple relation between the gravitational wave background from the hot Standard Model plasma and the mass and spin of dark matter that interacts only through gravity. It argues that future ultra-high-frequency gravitational wave detectors could test this minimal dark matter scenario.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reheating corrections are asserted to be O(1), but the paper's own Tmax formula gives log(Tmax/Trh) ~ 5-7 for viable heavy DM, so the claimed universal relation in Eq. (8) is not established.","rationale":"The algebraic step leading to Eq. (8) is straightforward and correct: solving Eq. (6) for Trh and substituting into Eq. (7) gives the stated combination of Omega_DM, alpha, mDM, f, and eta_hat. The cited formulas for the thermal DM abundance and GW spectrum are from published work, so the basic derivation is not the weak point. The load-bearing assumption is that a realistic reheating phase changes this relation only by an O(1) prefactor. The paper contains an explicit limitation statement in the paragraph after Eq. (8) and then dismisses it with a verification that is not shown. Importantly, the paper's own formula for Tmax/Trh, when combined with the Trh values implied by Eq. (6) for heavier-but-valid DM masses (e.g., mDM ~ 10^10 GeV implies Trh ~ 10^11 GeV and H(Trh) ~ 10^4 GeV), gives Tmax/Trh ~ 100-1000 for a high-scale inflation benchmark. Since the paper states that the GW amplitude receives a log(Tmax/Trh) correction, the correction is not necessarily O(1) in the parameter region where Eq. (8) is supposed to apply. This is not a disagreement with external consensus; it is an internal consistency check of the paper's own scaling formulas. A numerical integration of the coupled Boltzmann equations during reheating would settle the issue directly. The reader's verdict of CONDITIONAL already captures the need for a quantitative reheating treatment; the present stress-test strengthens that condition by showing that the 'O(1)' claim is likely to fail in part of the valid parameter space.","tokens_in":6568,"tokens_out":12681,"duration_ms":141534,"concrete_test":"Implement the full Boltzmann system for DM and graviton number densities during inflaton-dominated reheating with a quadratic potential, using the same interaction rates gamma_DM = alpha T^8/MP^4 and gamma_h = c T^6/MP^2, for a benchmark Hinf = 10^14 GeV and scan over Trh = 10^9-10^15 GeV (or equivalently over mDM = 10^4-10^11 GeV via Eq. (6)). Compute the present-day ratio R = Omega_GW h^2 / [(Omega_DM h^2 /0.12)^(1/3) (alpha/3e-3)^(-1/3) (mDM/1e9)^(-1/3) (f/1e11)^3 eta_hat] and compare with 8.6e-11. If R departs from 8.6e-11 by more than a factor of 2 for any mDM <= Trh, the O(1) correction claim in the paragraph after Eq. (8) is falsified and Eq. (8) needs a model-dependent prefactor.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is Eq. (8), which eliminates Trh between the thermal DM abundance and the thermal GW spectrum. Its validity requires that any non-instantaneous reheating phase changes both amplitudes only by an O(1) factor. The paper asserts this without derivation: 'We have verified that the impact of the reheating phase introduces, however, only a minor correction factor of O(1) to Eq. (8).' The attached formula for Tmax in the same paragraph undermines that assertion. For a quadratic inflaton, Tmax/Trh = (3/8)^(2/5) [Hinf/H(Trh)]^(1/4), and with the CMB bound on Hinf and the values of Trh implied by Eq. (6) for mDM of order 10^10 GeV (still within mDM <= Trh), this ratio is ~10^2-10^3, so log(Tmax/Trh) ~ 5-7. The paper itself states that GW production receives a log(Tmax/Trh) correction. A logarithmic factor of 5-7, combined with the factor-2 DM correction, changes the predicted Omega_GW by nearly an order of magnitude, not by an O(1) factor. Moreover, because the correction depends on Trh (hence on mDM through Eq. (6)), the m^{-1/3} scaling and the quoted mass thresholds in the Results section are not robust. The relation is therefore not a parameter-free 'cogenesis' prediction until the reheating dependence is quantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":6892,"tokens_out":14305,"duration_ms":140804,"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":[{"comment":"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.","section":"Cogenesis, paragraph after Eq. (8)"},{"comment":"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.","section":"Results, paragraph after Eq. (6)"}],"minor_comments":[{"comment":"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.","section":"Eqs. (7)–(8)"},{"comment":"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.","section":"Results, experimental reach"},{"comment":"In the first paragraph, \"the later depends on the mass and spin\" should be \"the latter depends on the mass and spin.\"","section":"Introduction"},{"comment":"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.","section":"Results, Fig. 2 caption and text"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about reheating is legitimate, although the log(Tmax/Trh) values are smaller than the note suggests: they are about 1–2.5 for the plotted masses under the quoted CMB bound. The load-bearing issue is not that the correction is numerically huge but that it is mass-dependent and unquantified, which weakens the clean m_DM^-1/3 relation. The arithmetic inconsistency in the m_DM lower bound is a separate, easily fixable error. The paper is concise and the core idea is sound; it should be publishable after the reheating correction is either derived or properly bounded and the numerical statements are corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nShort version: this is a correct, modest paper that takes two known results—the thermal freeze-in abundance for purely gravitational DM and the thermal GW background from the SM plasma—and eliminates the common reheating temperature to get a relation between GW amplitude, DM mass, and spin. Eq. (8) is algebraically correct, and the paper is honest that it's built on earlier work. The stress-test note claiming large reheating corrections (log(Tmax/Trh) ~ 5–7) does not hold up. Using the paper's own Tmax formula with the current CMB bound on Hinf and H(Trh) ∝ Trh^2 gives Tmax/Trh ≲ 10 for the relevant Trh range, so the log correction is of order 1–2. The author's O(1) claim is plausible, though saying \"we have verified\" without showing the derivation is a presentation weakness.\n\nWhat's actually new is the observation that the thermal GW background at ~100 GHz encodes the mass and spin of gravitational DM once the relic abundance is fixed. That's a useful condensation of existing formulas. The paper also draws concrete conclusions: the CMB bound on Trh forces m_DM ≳ 3×10^4 GeV, and a null result at ~10^11 Hz with Ω_GW h^2 > 10^-10 would rule out m_DM ≲ 10^6 GeV. These are clean, testable statements.\n\nThe soft spot is the experimental side. The paper talks about \"future ultra-high-frequency GW experiments\" without showing any projected sensitivity curves, so the probing claim is not quantified. Adding a comparison to concrete proposals would have made the paper much more useful. The spin dependence is also modest—the curves in Fig. 2 overlap significantly—so distinguishing spins would require good energy resolution, which the paper doesn't discuss.\n\nOverall, this is a short, correct, modest paper. It doesn't present a new mechanism, but it states a clean connection that researchers in gravitational DM or high-frequency GW phenomenology will likely cite. It deserves normal peer review, not a desk reject. I'd recommend accept after minor revision, provided the reheating correction is either shown explicitly or flagged as not universal.","headline":"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.","tokens_in":7384,"tokens_out":6961,"would_cite":true,"duration_ms":60794,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","04.30.-w"],"model":"deepseek-v4-flash","headline":"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…","keywords":["gravitational dark matter","thermal gravitational waves","ultra-high frequency gravitational waves","stochastic gravitational wave background","gravitational production","cogenesis","reheating temperature","spin dependence"],"falsifier":"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.","tokens_in":6380,"feed_emoji":"📡","tokens_out":12604,"duration_ms":103822,"temperature":0.7,"pith_summary":"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.","feed_headline":"One formula ties dark matter's mass and spin to gravitational waves","feed_subtitle":"If dark matter is purely gravitational, its abundance locks in the gravitational-wave amplitude near 100 GHz.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the spin-dependent coefficient $\\alpha$ used in the dark-matter abundance formula and in Eq. (8).","marker":"[9]"},{"why":"Supplies the thermal graviton production rate, the gravitational-wave spectrum formula, and the spectral-shape factor $\\hat{\\eta}(f)$ entering Eq. (7) and Eq. (8).","marker":"[10–12, 20]"},{"why":"CMB tensor-to-scalar ratio bound sets $T_{\\rm rh}\\lesssim5.5\\times10^{15}$ GeV, used to derive the lower dark-matter mass bound.","marker":"[41]"},{"why":"Provides the factor-of-two correction to the dark-matter abundance during reheating, supporting the claim that reheating changes Eq. (8) only at order one.","marker":"[30]"},{"why":"Establishes the existence of a thermal bath at BBN, the empirical anchor for assuming a radiation-dominated early Universe.","marker":"[4]"}],"fun_headline_variants":["Thermal plasma ties dark matter mass and spin to gravitational waves","Gravitational dark matter's mass and spin read out in gravitational waves","Ultra-high frequency GWs expose mass and spin of gravitational dark matter","A single formula links dark matter properties to gravitational-wave background","Dark matter's mass and spin encoded in ultra-high-frequency gravitational waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Thermal plasma ties dark matter mass and spin to gravitational waves","Gravitational dark matter's mass and spin read out in gravitational waves","Ultra-high frequency GWs expose mass and spin of gravitational dark matter","A single formula links dark matter properties to gravitational-wave background","Dark matter's mass and spin encoded in ultra-high-frequency gravitational waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001005,"raw_usage":{"total_tokens":4232,"prompt_tokens":908,"completion_tokens":3324,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":3234}},"tokens_in":524,"tokens_out":3324,"duration_ms":20982,"temperature":1.0,"reasoning_tokens":3234,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:01:13.252551+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}