{"id":"3dcb5c9b-00b5-4b3e-8050-9fef5719ff8f","arxiv_id":"2506.14366","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Cosmological gravitational wave backgrounds from phase transitions, domain walls, and condensate fragmentation are capped far below the astrophysical foreground by requiring that compact dark matter subhalos, which SKA could probe, are not overproduced.","lead":"This paper argues that the gravitational wave background from the early universe, hidden under the foreground from supermassive black hole binaries, can instead be bounded by looking for dark matter clumps that the same sources create. It finds these indirect bounds are several orders of magnitude below the foreground, so future pulsar timing with SKA could constrain new physics that direct gravitational wave searches would miss.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The F<=1 upper limits depend on an unverified, self-cited sigma_H(alpha_*, beta/H_*) relation for FOPTs, deferred to an absent supplemental appendix; the conclusion that Omega_GW lies orders below foreground is therefore not independently checkable.","rationale":"The reader's weakest assumption identified the mapping from GW source parameters to sigma_H, in particular the dependence on Ref. [84], as not reproduced. My stress-test concurs that this is the most load-bearing point, because all quantitative upper limits and the 'orders of magnitude below foreground' conclusion are generated through Eq. (3) using that mapping. The manuscript explicitly refers to a supplemental material for the detailed sigma_H expressions, but the v1 posting contains no such supplemental file, so the central derivation is omitted. The zc=500 choice is also flagged in the text as an assumption that changes the limits, but no error budget or physical justification is provided. These are correctness risks, not mere presentation issues. The separate observation that the F=1 bounds are no-overproduction conditions rather than SKA measurements is a framing concern, but it does not by itself invalidate the physical constraints if the sigma_H mapping is correct. Because the reader already reached a CONDITIONAL verdict based on this same weakness, my read does not change the verdict; however, the concrete test above would determine whether the concern actually lands and whether the conditional acceptance should be upheld or replaced by a different assessment.","tokens_in":12609,"tokens_out":7389,"duration_ms":75114,"concrete_test":"Run, or reproduce from public data, a lattice simulation of a first-order phase transition with the same parameters as Ref. [84] and measure the smoothed DM density contrast sigma_H at horizon entry as a function of alpha_* and beta/H_*. Then recompute the F<=1 upper limits on Omega_GW using Eq. (3) with zc=500. If the recomputed limits rise above the SMBHB foreground band for any of the three beta/H_* values, the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that cosmological SGWBs are constrained orders of magnitude below the astrophysical foreground rests entirely on the relation between GW-source parameters and the horizon-entry density contrast sigma_H used in Eq. (3). For FOPTs the paper adopts sigma_H proportional to alpha_* with a numerical normalization taken from Ref. [84], a self-authored PRL, and states that the precise relationship was obtained by the numerical methods of Ref. [84]. No fit formula, simulation data, or reproduction is given in the text; Appendix A lists only the GW spectra, not the sigma_H relations. The paper says details are in the supplemental material, but the arXiv v1 contains no supplemental material. Because F in Eq. (3) is a Gaussian tail integral, the derived upper limit on alpha_* (and hence on Omega_GW) is exponentially sensitive to the assumed sigma_H normalization; a factor-of-2 change in sigma_H can move the allowed Omega_GW by orders of magnitude. The same fragility affects the adopted zc=500, which the authors acknowledge is manually selected and which directly controls the collapse threshold. Since the headline conclusion is that the bounds lie several orders of magnitude below the foreground, the missing derivation of sigma_H is the single point on which the result depends. The SKA framing is also overstated: setting F=1 is a no-overproduction condition, not an SKA measurement, so the title's claim that SKA constrains the backgrounds is not supported by the actual calculation. The physical constraints may nonetheless be correct, but they cannot be checked from this manuscript.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new indirect way to bound cosmological stochastic gravitational-wave backgrounds (SGWBs) in the nHz/µHz band below the astrophysical foreground from supermassive black hole binaries. The method uses the fact that GW sources such as first-order phase transitions, domain walls, and scalar condensate fragmentation also generate small-scale density perturbations; if these perturbations are large enough, they form compact dark-matter subhalos whose abundance can be constrained by future SKA observations. The authors set the condition F=1 (no overproduction of subhalos), adopt zc=500 as the last collapse redshift, and translate this into upper limits on Omega_GW for bubble collisions, sound waves, and domain walls, claiming these limits lie orders of magnitude below the astrophysical foreground. Scalar-induced GWs are also discussed. The paper relies on a sigma_H vs alpha_* relation for first-order phase transitions taken from the authors' prior work (Ref. [84]) and on the Gaussian-tail integral in Eq. (3).","tokens_in":12853,"tokens_out":5787,"duration_ms":56511,"significance":"If the quantitative inputs were fully provided and verified, the paper would introduce a genuinely interesting and complementary probe: using compact dark-matter substructures to constrain cosmological GW sources below the SMBHB foreground, a regime that GW detectors alone may never access. The idea is physically well motivated, the conservative choices (F=1, vw=1, zc=500) are transparent, and the paper is clearly written. The main strength is the conceptual connection between GW-source parameters and small-scale density perturbations, which is a valuable direction for the field. However, the central FOPT bounds are exponentially sensitive to the sigma_H relation, which is not reproduced in the manuscript or the (absent) supplemental material, and some headline claims (SKA constraining the backgrounds, constraints on condensate fragmentation) go beyond what is actually computed. With the missing input supplied and the claims carefully rescaled, the result could be an important addition to the nHz/µHz SGWB literature.","major_comments":[{"comment":"The upper limits on Omega_GW from first-order phase transitions are set through Eq. (3), which is exponentially sensitive to sigma_H, yet the relation sigma_H(alpha_*, beta/H_*) is not given anywhere in the text or Appendix A; the paper states only that it follows from the numerical methods of Ref. [84] and refers to a supplemental material that is absent from arXiv v1. Because a factor-of-2 uncertainty in sigma_H changes the allowed alpha_* and hence Omega_GW by orders of magnitude, the central quantitative claim is not independently checkable. Please include the explicit fit formula or numerical data used, and an estimate of its uncertainty, or clearly mark the bounds as conditional on that relation.","section":"Constraints on gravitational waves of cosmological origin (FOPT paragraph)"},{"comment":"The analysis derives upper limits by imposing F=1, i.e., a no-overproduction requirement, and the SKA sensitivity curve (brown dashed line in Figs. 1 and 2) is not used in the derivation; therefore the title/abstract claim that SKA constrains the cosmological backgrounds is not supported by the calculation. Please either rephrase the claim as a forecast of what SKA could constrain once subhalo abundance is measured, or propagate an actual SKA-based measurement or upper limit on F into the Omega_GW bounds.","section":"Title and Abstract / Constraints section"},{"comment":"The bounds depend on the manually chosen zc=500, which enters Eq. (1) through (1+zc)^3 and the collapse threshold, and the text acknowledges that smaller zc gives tighter constraints; with a plausible range of zc the 'orders of magnitude below foreground' conclusion may shift. Please show the dependence of the derived Omega_GW limits on zc (e.g., a band or a few representative values) and justify the adopted value with a quantitative simulation-based statement rather than the qualitative reference to 'recent simulations [74, 75]'.","section":"Observing the compact DM subhalos with pulsar timing / Constraints section"},{"comment":"For scalar condensate fragmentation, the paper does not actually derive an upper limit from F≤1; it states that for Omega_phi=1 the peak is roughly Omega_GW≈10^-12, which is a fiducial model prediction, not a subhalo-abundance bound. Since the abstract and conclusion claim constraints on 'various sources' including condensate fragmentation, either derive the F≤1 bound for this source using delta_H ~ Omega_nabla (k_res/(a H_*))^{-3/2}, or restrict the claim to the sources for which a bound is computed.","section":"Constraints on gravitational waves of cosmological origin (condensate fragmentation paragraph)"}],"minor_comments":[{"comment":"The abstract contains typographical errors: 'supermaissive' should be 'supermassive' and 'convinced gravitational wave background' should be 'convincing gravitational-wave background'; the Introduction also contains an orphan sentence fragment, 'not applicable to the nHz/µHz bands.', which appears to be a leftover and should be removed or completed.","section":"Abstract and Introduction"},{"comment":"In the paragraph after Eq. (2), 'viral radius' should be 'virial radius'.","section":"Observing the compact DM subhalos with pulsar timing"},{"comment":"The right-panel label 'Fagmentation Temperature' in the Figure 2 caption should be 'Fragmentation Temperature'.","section":"Figure 2 caption"},{"comment":"The paper repeatedly refers to 'supplemental material' for the sigma_H relations and for detailed GW spectra, but no supplemental material is present in arXiv v1; please ensure it is included at resubmission, and ideally move the sigma_H relation and the condensation-fragmentation derivation into the main text or Appendix A.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The manuscript depends heavily on a previously published sigma_H relation that is not reproduced, and the arXiv submission lacks the promised supplemental material. Given the exponential sensitivity of the Gaussian-tail integral in Eq. (3), I would require the relation and its uncertainty to be provided before publication. The title also overstates the role of SKA, since the constraints are derived from F=1 rather than from SKA sensitivity; this can be corrected by reframing the claims as forecasts. These issues are fixable within the manuscript's scope, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a genuinely new application — using compact-subhalo abundance to cap Ω_GW from FOPTs, domain walls, and condensate fragmentation several orders below the SMBHB foreground. The physics is simple and the idea is sound: the same processes that emit GWs also induce density perturbations that collapse into dark matter subhalos, so demanding F ≤ 1 translates into an upper limit on the GW spectrum. I hadn't seen that done for these source classes, and the qualitative conclusion is plausible.\n\nThe weak spot is the FOPT section. The limit there is set by the mapping σ_H ∝ α_* with a numerical normalization borrowed from the authors' own Ref. [84]. That relation is not shown anywhere in this arXiv v1; the text says to see the supplemental material, and there is none attached. That matters because F in Eq. (3) is a Gaussian tail integral, so the derived upper limit on α_* — and hence on Ω_GW — is exponentially sensitive to the assumed normalization. A factor-of-2 change in σ_H can shift the bounds by orders of magnitude. This is not a minor point: it is the anchor of the headline claim. The domain-wall and condensate limits are better, since the σ_H expressions are explicit, but they too carry the manually selected z_c = 500 with no band, and no error bars anywhere.\n\nThe SKA framing also outruns the calculation. The present curves come from setting F = 1, a no-overproduction condition, not from any SKA measurement. SKA is a plausible future route to making the bound observational, but the actual numbers are theory, not detection. The abstract and title should say 'could be constrained with SKA' rather than implying SKA is doing the work. Minor sloppiness — a broken sentence in the introduction, a half-referenced scalar-induced case — doesn't change the picture but doesn't help.\n\nStill, the intellectual content is honest and the method is worth a serious referee. The authors are upfront that the bounds depend on z_c and that F = 1 is a conservative choice. That is not a circular argument; it's a well-posed mapping from subhalo abundance to GW parameter space, with one nontrivial input that is currently uncheckable from the manuscript alone.\n\nRecommendation: send to referees, but require (1) the actual σ_H–α_* relation in the text or a real supplemental file, (2) a sensitivity band for z_c and σ_H normalization, and (3) a title/abstract that doesn't overstate SKA's role. If those are fixed, I'd cite and use the results.","headline":"A genuinely new indirect bound on cosmological SGWBs, but the central FOPT constraint depends on an unreproduced self-cited relation and the SKA framing outruns the calculation.","tokens_in":13456,"tokens_out":5380,"would_cite":false,"duration_ms":51324,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that requiring dark-matter subhalo abundance not to exceed the total dark matter constrains nHz/µHz cosmological gravitational-wave backgrounds to lie orders of magnitude below the astrophysical foreground.","keywords":["gravitational-wave background","compact subhalos","pulsar timing arrays","Square Kilometre Array","first-order phase transitions","domain walls","scalar condensate fragmentation","dark matter substructure"],"falsifier":"A dedicated numerical simulation that recomputes the horizon-entry density contrast $\\sigma_H$ for first-order phase transitions at moderate $\\beta/H_*$ and $\\alpha_*$ would settle the numerical link; if $\\sigma_H$ came out several times smaller than the adopted $\\sigma_H \\propto \\alpha_*$ scaling, the derived $\\Omega_\\mathrm{GW}$ upper limits would shift upward by a corresponding factor and could approach the astrophysical foreground. Alternatively, a future SKA measurement that places a strong upper limit on compact-subhalo abundance well below $F=1$ would confirm the calibration of the method and tighten the bounds.","tokens_in":12334,"feed_emoji":"📡","tokens_out":18691,"duration_ms":158514,"temperature":0.7,"pith_summary":"This paper proposes a new observational route to cosmological gravitational-wave backgrounds in the nanohertz-to-microhertz band, where direct detection is blocked by an irreducible foreground from supermassive black hole binaries. The route is indirect: the same violent processes that emit gravitational waves—first-order phase transitions, domain walls, and scalar condensate fragmentation—also seed density perturbations that collapse into compact dark-matter subhalos. Requiring the subhalo abundance $F$ not to exceed 1, a level the upcoming Square Kilometre Array can plausibly test, turns into upper limits on the gravitational-wave energy density $\\Omega_\\mathrm{GW}$ that lie several orders of magnitude below the astrophysical foreground. If the method holds, cosmological backgrounds in this band are effectively hidden from gravitational-wave detectors, while SKA's subhalo observations become the tightest indirect constraints on the underlying new physics.","feed_headline":"SKA subhalo limits put cosmological GWs orders below the foreground","feed_subtitle":"If SKA counts compact dark-matter subhalos, nHz cosmological backgrounds are buried under the binary foreground.","key_machinery":"The central machinery is the compact-subhalo abundance integral. The density perturbation at horizon entry, characterized by the smoothed variance $\\sigma_H$, is fed through a Gaussian threshold integral with threshold $\\delta_{\\mathrm{min}}$ to give the fraction $F$ of dark matter locked in subhalos, using the Moore profile (a cusped halo density profile) and its mass–radius relation from numerical simulations. Each gravitational-wave source supplies its own $\\sigma_H$: $\\sigma_H \\propto \\alpha_*$ for first-order phase transitions, $\\sigma_H \\simeq 2A\\sigma/(3M_{\\mathrm{pl}}^2 H_{\\mathrm{ann}})$ with $A \\simeq 0.8$ for domain walls, and a causality-limited $\\delta_H \\sim \\Omega_{\\nabla}(k_{\\mathrm{res}}/(aH_*))^{-3/2}$ for scalar condensate fragmentation. Setting $F=1$ inverts the chain and converts a subhalo-abundance constraint into an upper bound on $\\Omega_\\mathrm{GW}$; the last-collapse redshift $z_c$ is fixed to 500 as a moderate choice, and the SKA sensitivity enters through its ability to probe order-unity subhalo abundances via pulsar timing.","core_discovery":"This paper claims that requiring the abundance $F$ of compact dark-matter subhalos not to exceed unity places upper limits on the present-day gravitational-wave energy density $\\Omega_\\mathrm{GW}$ of cosmological sources in the nanohertz-to-microhertz band, and that these limits fall orders of magnitude below the astrophysical foreground from supermassive black hole binaries. The limits are derived for bubble collisions and sound waves from first-order phase transitions, for domain walls, and for scalar condensate fragmentation, using the relation between each source's horizon-entry density contrast $\\sigma_H$ and its gravitational-wave spectrum. For domain walls the limits are typically four to six orders of magnitude tighter than the foreground; for condensate fragmentation the peak $\\Omega_\\mathrm{GW} \\approx 10^{-12}$ sits about four orders below it; for first-order phase transitions the limits are several orders tighter, with smaller $\\beta/H_*$ giving stronger bounds. The analysis deliberately takes $F=1$ as a conservative normalization, so the stated bounds are upper limits rather than detections.","pith_inferences":["A natural extension of the same inversion is cosmic strings, which the paper lists as a candidate source but does not constrain numerically; applying the same $\\sigma_H$-to-$F$ chain would likely push their background equally far below the foreground.","Since the limits depend monotonically on $z_c$, a future SKA measurement of subhalo abundance could be inverted to estimate the last collapse redshift, turning a calibration parameter into an observable.","If independent probes such as gamma-ray annihilation limits push the allowed subhalo fraction below $F=1$, the same framework would scale the $\\Omega_\\mathrm{GW}$ bounds downward, so the paper's limits are conservative rather than maximal."],"forward_implications":["If the central claim holds, cosmological gravitational-wave backgrounds in the nHz/µHz band from first-order phase transitions, domain walls, and scalar condensate fragmentation are generically too weak to be seen above the supermassive-black-hole-binary foreground, so future nHz/µHz detectors will not resolve them as gravitational waves.","The Square Kilometre Array's ability to constrain compact-subhalo abundance at the order-one level becomes a competitive indirect probe of new-physics parameter spaces at MeV–GeV energy scales, where direct gravitational-wave searches are blinded by the foreground.","According to the paper's conclusion, the upper limits also reduce the uncertainty of the astrophysical component of the nHz/µHz background, yielding more concrete information about the formation and evolution of supermassive black hole binaries.","Because the bounds tighten for smaller values of the last-collapse redshift $z_c$, any future simulation or observation that pins down the subhalo formation epoch will sharpen the limits accordingly."],"supporting_citations":[{"why":"Supplies the numerical relation $\\sigma_H \\propto \\alpha_*$ that converts first-order phase-transition strength into the density-perturbation amplitude entering the subhalo-abundance integral.","marker":"[84]"},{"why":"Provides the Moore-profile parameters and mass formula used to compute the fraction of dark matter in compact subhalos.","marker":"[74]"},{"why":"Companion simulation work supporting the subhalo profile and the last-collapse redshift constraint adopted in the analysis.","marker":"[75]"},{"why":"Gives the threshold $\\delta_{\\mathrm{min}}$ used in the Gaussian abundance integral for compact subhalo formation.","marker":"[76]"},{"why":"Supplies the SKA pulsar-timing baseline parameters and signal-to-noise scalings that justify the order-one subhalo-abundance probe.","marker":"[77]"},{"why":"Provides the domain-wall simulation constants ($A \\approx 0.8$, $\\tilde{\\epsilon}_{\\mathrm{gw}} \\approx 0.7$) and the peak gravitational-wave amplitude formula used for the domain-wall bound.","marker":"[70]"},{"why":"Gives the unified gravitational-wave spectrum of scalar condensate fragmentation used in the $\\Omega_\\mathrm{GW}$ estimate.","marker":"[104]"},{"why":"Provides the supermassive-black-hole-binary gravitational-wave spectrum used as the astrophysical foreground reference.","marker":"[88]"},{"why":"Provides the astrophysical foreground band from supermassive black hole binaries with environmental effects against which the cosmological upper limits are compared.","marker":"[60]"}],"fun_headline_variants":["SKA subhalo abundance constrains cosmic GWs below foreground","SKA subhalo census tightens limits on cosmological GW backgrounds","Subhalo abundance ties SKA to nHz GW limits under the foreground","SKA's dark-matter subhalos bound GW backgrounds beneath the noise"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical relation between gravitational-wave source parameters and the dark-matter density fluctuation that seeds subhalos, taken from external simulations, is correct, and that the manual choice of the last collapse redshift $z_c = 500$ is representative.","fun_headline_variants_meta":{"raw":{"variants":["SKA subhalo abundance constrains cosmic GWs below foreground","SKA subhalo census tightens limits on cosmological GW backgrounds","Subhalo abundance ties SKA to nHz GW limits under the foreground","SKA's dark-matter subhalos bound GW backgrounds beneath the noise"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001048,"raw_usage":{"total_tokens":4394,"prompt_tokens":929,"completion_tokens":3465,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":545,"completion_tokens_details":{"reasoning_tokens":3386}},"tokens_in":545,"tokens_out":3465,"duration_ms":28512,"temperature":1.0,"reasoning_tokens":3386,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:54:35.740053+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A dedicated numerical simulation that recomputes the horizon-entry density contrast $\\sigma_H$ for first-order phase transitions at moderate $\\beta/H_*$ and $\\alpha_*$ would settle the numerical link; if $\\sigma_H$ came out several times smaller than the adopted $\\sigma_H \\propto \\alpha_*$ scaling, the derived $\\Omega_\\mathrm{GW}$ upper limits would shift upward by a corresponding factor and could approach the astrophysical foreground. Alternatively, a future SKA measurement that places a strong upper limit on compact-subhalo abundance well below $F=1$ would confirm the calibration of the method and tighten the bounds.","supporting_citations":[],"review_version":1}