{"id":"9aa57c83-18d6-4dd2-b80c-a6fc438c7831","arxiv_id":"2607.03533","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.5,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"A stochastic GWB suppresses PTA reach to DM substructure by 1–3 orders of magnitude relative to white-noise forecasts, with dynamic Shapiro least affected near 10^{-2} M_⊙.","lead":"Pulsar timing arrays lose one to three orders of magnitude in sensitivity to dark-matter subhalos once the stochastic gravitational-wave background is treated as red noise. The paper supplies analytic scalings and SKA forecasts that show the dynamic Shapiro channel remains the least degraded probe near 0.01 solar masses.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The paper’s strongest claim is a quantitative, regime-by-regime forecast of GWB-induced degradation. The analytic integrals I_{n,γ} and I^{(2)}_{n,γ} (App. C) locate the red-noise weight and produce the min(γ,n-1) exponents of Table I; the Monte-Carlo of Sec. V (full proper-time observable, timing-model projection via N_⊥^{-1}) recovers those scalings to O(1) (Fig. 4). Variation across the three posterior points is explicitly shown to be order-unity (Fig. 3). The idealizations flagged by the reader—point-mass subhalos, common speed, single-pulsar statistic, neglect of intrinsic red noise—are stated in the text and enter only as O(1) prefactors or deferred companion-paper corrections; none alters the power of the (f_⋆T) suppression that drives the 1–3 order effect. Because the load-bearing mathematics is self-contained and the numerical cross-check is consistent, no adjustment to the ACCEPT verdict is warranted.","tokens_in":40225,"tokens_out":554,"duration_ms":4887,"concrete_test":"Recompute the SKA deterministic reach of Fig. 3 after replacing the pure GWB red noise in Eq. (8) by GWB + a representative pulsar-intrinsic red-noise power law (A_red, γ_red drawn from the NANOGrav 15-yr pulsar noise posteriors). If the f_DM envelope shifts by more than a factor of ~3 at the dynamic-Shapiro knee (M~10^{-2} M_⊙), the quoted suppression range would need widening; otherwise the claim stands.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that a GWB suppresses PTA reach to DM substructure by 1–3 orders of magnitude relative to white-noise forecasts, with only O(1) variation across the NANOGrav 15-yr posterior, and that dynamic Shapiro is least suppressed—is supported by transparent analytic scalings (Table I, App. C–D) that match the full gauge-invariant Monte-Carlo reach (Figs. 3–4) deep in each regime. The reader’s listed caveats (point masses, neglected intrinsic red noise, single-pulsar statistic) are already quantified by the authors as O(1) effects and do not overturn the order-of-magnitude suppression. No internal inconsistency or hidden assumption that would reverse the 1–3 order claim was found.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper quantifies how a stochastic gravitational-wave background (GWB), treated as red noise, degrades pulsar-timing-array (PTA) sensitivity to dark-matter (DM) substructure. It derives the full gauge-invariant proper-time residual for a transiting subhalo (Doppler, Shapiro, and Einstein terms), builds an SNR framework that marginalizes the quadratic timing model under a red-noise-weighted inner product, and obtains closed-form scaling relations for the static, dynamic, and stochastic regimes (Table I). These are validated against a full Monte-Carlo reach for an SKA-like array at three points of the NANOGrav 15-year GWB posterior. The central result is that the GWB suppresses the substructure reach by one to three orders of magnitude relative to white-noise-only forecasts, with only O(1) variation across the posterior; the dynamic Shapiro channel is least suppressed and peaks near 10^{-2} M_⊙.","tokens_in":40442,"tokens_out":1288,"duration_ms":20384,"significance":"Given the NANOGrav evidence for a nHz GWB, prior white-noise PTA forecasts for DM substructure are no longer realistic. This work supplies the first transparent analytic account of the red-noise penalty, with scalings that match full gauge-invariant numerics deep in each regime (Fig. 4) and that can be reused for other PTA benchmarks. Strengths include the explicit gauge-invariance check (App. B), the red-noise integral asymptotics (App. C), the projected-norm calculations (App. D), and the Monte-Carlo sampling and Gaussian-validity diagnostics (App. E). The conclusion that f_DM ≲ 1 remains difficult even for SKA-scale arrays, while still improving substantially on existing NANOGrav limits, is a useful and falsifiable guide for survey design.","major_comments":[{"comment":"Sec. II A and Eq. (7)–(8): pulsar-intrinsic red noise is set identically to zero, so the single-pulsar PSD is white noise plus the GWB diagonal only. For the absolute SKA reach curves in Fig. 3 this is an optimistic floor; many MSPs retain measurable spin noise at nHz frequencies. A short estimate (or a third curve) showing how an intrinsic red component comparable to current PTA levels would further shift f_DM would make the absolute forecasts more robust. The relative 1–3 order GWB suppression itself is not threatened, because it is controlled by the GWB term in S_n(f).","section":"Sec. II A, Eqs. (7)–(8)"},{"comment":"Sec. IV–V and the companion-paper deferral: the Earth Doppler term and the full-array (Hellings–Downs-conditioned) statistic are omitted, with the claim that both change the single-pulsar reach only at O(1). That claim is plausible for well-separated pulsars, but the dynamic Shapiro knee near 10^{-2} M_⊙ is the only place where the projected reach approaches f_DM ∼ 1. A one-paragraph quantitative bound (even from a simplified two-pulsar estimate) on how much the Earth term or inter-pulsar conditioning could move that knee would strengthen the statement that f_DM < 1 remains a challenge.","section":"Sec. IV–V"}],"minor_comments":[{"comment":"Table I: the SNR scalings with T, f_DM, M, and N_P are very useful; adding a short footnote that the numerical prefactors are dropped (∼ rather than =) would prevent readers from treating the table as exact.","section":"Table I"},{"comment":"Fig. 2 caption is dense; labeling the three panels explicitly as “static / dynamic / stochastic” in the figure itself (not only in the caption) would improve readability.","section":"Fig. 2"},{"comment":"Eq. (38): the numerical prefactor “20 · 0.25^{1/γ} …” is hard to parse at a glance; writing f_⋆T ≃ 20 × (0.25)^{1/γ} × … would clarify the structure.","section":"Eq. (38)"},{"comment":"Sec. III: the Einstein terms are correctly identified as O(v)-suppressed relative to Doppler, but a one-line numerical check (e.g. ratio of norms for a typical flyby) would make the neglect fully transparent for non-specialists.","section":"Sec. III"},{"comment":"App. E 3 / Fig. 7: the N_Q ≥ 10 validity cut is well motivated; stating in the main text (near Fig. 4) that the dotted stochastic segments are extrapolations, not quantitative forecasts, would reduce the chance of misreading.","section":"Fig. 4, App. E 3"},{"comment":"References: Ref. [26] (“in preparation”) carries the full proper-time derivation; if a public draft or arXiv version exists by revision time, citing it would help reproducibility of App. B.","section":"Refs. / App. B"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is technically solid and well matched to a theory/astro-ph.CO journal. The companion papers [26, 42] are cited for gauge-invariant proper time and array-correlated statistics; that split is reasonable but means the present paper is not fully self-contained on the Earth term. I do not see novelty or citation issues that would affect the editorial decision. Minor revision is appropriate; I would not require a second full round if the intrinsic-red-noise and Earth-term points are addressed briefly."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing to take away is that once you put a realistic GWB into the noise budget, every earlier white-noise PTA forecast for dark-matter subhalos is optimistic by one to three orders of magnitude in f_DM. The authors give closed-form red-noise suppression factors (f★T)^(-min(γ,n-1)/2) that track the full Monte-Carlo reach to O(1) deep in each regime, and they show the dynamic Shapiro channel is the least damaged, peaking near 10^{-2} M_⊙.\n\nWhat is actually new is the full gauge-invariant proper-time observable (Doppler + Shapiro + Einstein) and the transparent projection onto the timing-model-orthogonal subspace under red noise. Prior forecasts either assumed white noise or ran numerical Bayes; here you get analytic envelopes for static/dynamic/stochastic Doppler and Shapiro that you can re-evaluate for any future array. The appendices are careful: gauge invariance is shown explicitly, the red-noise integrals are derived, and the Monte-Carlo sampling geometries are documented. Analytic and numerical curves agree where they should (Fig. 4). Citation pattern is honest about the earlier Dror/Ramani/Lee/Zurek series and the NANOGrav 15-yr posterior.\n\nSoft spots are real but secondary. Point-mass subhalos, common speed, neglect of pulsar-intrinsic red noise, and the single-pulsar statistic (Earth term and array correlations deferred to a companion) are all stated. The authors already treat them as O(1) effects; they do not reverse the order-of-magnitude claim. The SKA benchmark is conventional, so the absolute numbers are illustrative rather than definitive, but the design rule that emerges—array size still scales while white-noise residual is largely neutralized by the GWB—is useful.\n\nThis is for people who design PTA surveys or who quote DM-substructure limits from them. The math is solid enough that a serious referee should see it. I would cite the scalings and the reach curves; I would not treat the absolute f_DM numbers as final until the companion paper and a more realistic noise model appear. Send it to peer review.","headline":"Clean analytic + numerical quantification of how the GWB kills PTA DM-substructure reach by 1–3 orders of magnitude; the scalings and gauge-invariant observable are the real additions.","tokens_in":41089,"tokens_out":561,"would_cite":true,"duration_ms":6691,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A stochastic gravitational-wave background weakens pulsar-timing searches for dark-matter subhalos by one to three orders of magnitude.","keywords":["pulsar timing arrays","dark matter substructure","stochastic gravitational-wave background","Shapiro delay","Doppler residual","signal-to-noise ratio","red noise"],"falsifier":"If an SKA-like array, after 20 years and after full timing-model marginalization, recovers a substructure fraction near the white-noise-only forecasts rather than the one-to-three-order-worse curves of Fig. 3, the claimed GWB suppression is ruled out.","tokens_in":41131,"feed_emoji":"⏱️","tokens_out":657,"duration_ms":5777,"temperature":0.7,"pith_summary":"Pulsar timing arrays can sense dark-matter subhalos by the tiny shifts they leave in pulse arrival times. A recently detected stochastic gravitational-wave background acts as red noise that competes with those shifts. This paper derives the full gauge-invariant proper-time delay a subhalo produces and a signal-to-noise framework that accounts for both that red noise and the quadratic timing model fitted to each pulsar. Analytic scalings then show how sensitivity falls in the static, dynamic, and stochastic regimes, and numerical forecasts for an SKA-like array confirm a one-to-three-order suppression relative to white-noise-only projections. The dynamic Shapiro delay is least affected and still offers the best reach near 0.01 solar masses. Even so, future arrays improve existing limits by up to two orders of magnitude, keeping pulsar timing competitive for small-scale dark-matter structure.","feed_headline":"GWB red noise cuts PTA dark-matter reach by 1–3 orders","feed_subtitle":"Dynamic Shapiro delay still peaks near 0.01 solar masses; future arrays improve limits by up to 100\times","key_machinery":"A noise-weighted inner product that projects each timing residual orthogonal to the quadratic pulsar timing model, combined with the full gauge-invariant proper-time delay (Earth/pulsar Doppler + Shapiro + Einstein) of a transiting subhalo. The projection converts every signal into a low-frequency power-law spectrum whose SNR integral is then suppressed by the closed-form red-noise factor (f_star T)^{-min(gamma,n-1)} (or its squared counterpart for stochastic signals).","core_discovery":"The stochastic gravitational-wave background suppresses PTA sensitivity to dark-matter substructure by one to three orders of magnitude relative to white-noise-only forecasts; the precise factor depends only mildly (within a factor of three) on the background's amplitude and spectral index. Among all channels, the dynamic Shapiro signal suffers the smallest penalty and supplies the best sensitivity near 10^{-2} solar masses.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["GWB red noise slashes PTA DM substructure reach by 1-3 orders","Stochastic GWB cuts PTA dark-matter subhalo sensitivity 10-1000x","Dynamic Shapiro channel holds best PTA reach near 0.01 solar masses","GWB red noise degrades PTA DM substructure limits by 1-3 orders","PTA sensitivity to DM subhalos falls 1-3 orders under GWB red noise"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"Subhalos are treated as point masses all moving at the same speed, intrinsic pulsar red noise is ignored, and only the single-pulsar noise autocorrelation is used, so any extra degradation from diffuse profiles or residual red noise is left for later work.","fun_headline_variants_meta":{"raw":{"variants":["GWB red noise slashes PTA DM substructure reach by 1-3 orders","Stochastic GWB cuts PTA dark-matter subhalo sensitivity 10-1000x","Dynamic Shapiro channel holds best PTA reach near 0.01 solar masses","GWB red noise degrades PTA DM substructure limits by 1-3 orders","PTA sensitivity to DM subhalos falls 1-3 orders under GWB red noise"]},"model":"grok-4.5","effort":"low","cost_usd":0.00467,"raw_usage":{"total_tokens":1419,"prompt_tokens":865,"num_sources_used":0,"completion_tokens":94,"cost_in_usd_ticks":46700000,"prompt_tokens_details":{"text_tokens":865,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":460,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":865,"tokens_out":94,"duration_ms":3747,"temperature":1.0,"reasoning_tokens":460,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T01:48:30.141718+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"If an SKA-like array, after 20 years and after full timing-model marginalization, recovers a substructure fraction near the white-noise-only forecasts rather than the one-to-three-order-worse curves of Fig. 3, the claimed GWB suppression is ruled out.","supporting_citations":[],"review_version":1}