{"id":"eae62397-e11e-4d84-b499-95a063ca4582","arxiv_id":"2505.12678","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Golden dark sirens from next-generation gravitational-wave networks could constrain the cosmic dipole amplitude to about 10^-3 jointly with H0 and to about 10^-4 if H0 is fixed.","lead":"This paper simulates how many nearby, precisely localized gravitational-wave events called golden dark sirens future detectors would see, and how well those events could measure the cosmic dipole. It finds that networks built around Einstein Telescope and Cosmic Explorer could measure the dipole to about one part in a thousand while also estimating the Hubble constant.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unquantified host-galaxy identification probability is the load-bearing gap: the 6% dipole constraint from 39 golden dark sirens assumes every z_obs equals the source redshift, which the paper explicitly leaves uncalculated.","rationale":"The reader's weakest assumption is exactly the load-bearing concern: the golden-dark-siren method assumes the brightest galaxy in the localization area is the host, and the paper itself states that this probability is not quantitatively calculated. The dipole signal is tiny (10^-3), the event sample is only 39 events in the most optimistic network, and the likelihood in Eq. (3.4) depends directly on the observed redshift. A wrong host redshift therefore enters the central observable, not just a peripheral selection effect. I considered other weaknesses, such as the crude Fisher-matrix localization estimate, the small Poisson sample size, and the lack of released code or data; these are real but secondary. The host identification probability is singled out by the author, is directly used in the analysis, and is the place where the central claim is least secure. Since the paper's own conclusions flag this limitation, the appropriate outcome remains a conditional acceptance pending a quantitative host-probability calculation; no change from the reader's CONDITIONAL verdict is needed.","tokens_in":19070,"tokens_out":6072,"duration_ms":70100,"concrete_test":"Compute the host-galaxy identification probability for the exact selection used in Sec. 3.1 (90% localization area <= 0.06 deg^2, z <= 0.1) with a realistic local galaxy catalogue such as GLADE+, or a Schechter-sampled mock with L > L* galaxies and CBC hosts weighted by stellar mass. Then rerun the MCMC/chi-squared analysis, randomly replacing z_obs with the redshift of the brightest non-host galaxy for the measured misidentification fraction. If the recovered g shifts by more than about 1 sigma, or if its uncertainty grows beyond the claimed ~6%, the headline dipole constraint is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central forecast, that 1ET+2CEs can constrain the cosmic dipole to roughly 6% at g=10^-3 (Table 3, Fig. 7), rests on the golden-dark-siren selection of Sec. 3.1: each selected event has a 90% localization area within 0.06 deg^2, and the brightest galaxy in that area 'is then considered to be the host galaxy of the CBC.' Eqs. (3.4)-(3.5) then use the galaxy redshift z_obs to construct the rest-frame D0_L entering the chi-squared. If the brightest galaxy is not the true host, z_obs is systematically wrong, and the inferred D0_L is wrong by an amount that is not modeled. This is not a minor caveat: the dipole signal is only of order 10^-3 in D_L and z, while the redshift difference between a candidate galaxy and the true host is unbounded and can easily exceed that level. With only 39 events in the best network (Table 2), even a small number of misidentified hosts can produce a systematic shift in g comparable to or larger than the claimed 6% statistical uncertainty. The paper explicitly concedes this in Sec. 5: 'the probability for the brightest galaxy within the localization area of a golden dark siren to be its host galaxy is not quantitatively calculated.' The forecast is internally consistent, but its headline number is conditional on an unverified and potentially dominant systematic.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper forecasts the ability of next-generation ground-based gravitational-wave detector networks (1ET+2A#, 1ET+1CE+1A#, 1ET+2CE) to measure the cosmic dipole jointly with H0 using 'golden dark sirens': nearby (z<0.1) compact binary coalescences with 90% sky localization within 0.06 deg^2, whose host is assumed to be the brightest galaxy in the localization area. The authors simulate 10 years of detections, compute Fisher-matrix parameter uncertainties, inject a dipole of amplitude g=10^-3 in the CMB direction, and recover {H0,g,l_dip,b_dip} with an MCMC chi-squared analysis. They find that networks with at least two next-generation detectors yield meaningful constraints: with 1ET+2CEs, g is measured to about 6% and the dipole direction to about 120 deg^2, with H0 to about 0.11%. The paper also examines the bias that the cosmic dipole induces in standard dark-siren H0 measurements with LVK O4/O5 data and finds it negligible.","tokens_in":19389,"tokens_out":8944,"duration_ms":98239,"significance":"If the forecast holds, it would provide a new and independent way to probe the cosmic dipole without relying on number counts, complementing bright-siren forecasts. The injection-recovery setup is standard and clearly described, and the use of public simulation and sampling packages makes the analysis reproducible in principle. The headline numbers, however, are conditional on two unmodeled astrophysical systematics, host-galaxy identification and peculiar velocities, that can individually be of the same order as the claimed statistical uncertainty; the quantitative claims should therefore be read as idealised sensitivity forecasts rather than realistic error predictions.","major_comments":[{"comment":"The central forecast assumes that the brightest galaxy inside the 0.06 deg^2 localization area is the true host, so that the measured z_obs is the source redshift. The paper itself states in Sec. 5 that the probability for this identification is not quantitatively calculated. Because the dipole signal is only g=10^-3, while the redshift difference between a random bright galaxy and the true host can be much larger, a small misidentification fraction can produce a systematic error in g comparable to or larger than the claimed 6% uncertainty. The authors should either compute this probability using a realistic galaxy population or catalogue, or demonstrate explicitly that the forecast is robust to a reasonable misidentification rate.","section":"Sec. 3.1, Sec. 5, Eqs. (3.4)-(3.5)"},{"comment":"The likelihood in Eq. (3.4) treats z_obs as the cosmological-plus-dipole redshift, but at z<0.1 the peculiar velocities of host galaxies (rms about 300 km/s) induce a scatter of order 10^-3 in z, comparable to the injected dipole amplitude and much larger than the fractional distance uncertainty claimed for individual golden dark sirens. Averaging over about 39 events reduces this scatter, but it remains a dominant noise source for a g=10^-3 signal. The paper mentions peculiar velocities only as a future systematic in Sec. 5; to support the headline constraint, the forecast should include a peculiar-velocity term or a conservative scatter in the simulated z_obs and in the likelihood.","section":"Sec. 2.1, Sec. 3.2, Conclusions"},{"comment":"The chi-squared likelihood omits selection effects: the golden-dark-siren sample is defined by a 90% localization-area threshold, and both the SNR and the localization area depend on D_L and hence on the dipole amplitude and direction. The forecast uses a fixed event set selected without the dipole and then applies Eq. (3.4) without a detection-probability normalization, in contrast to the selection term in Eq. (2.7). Since the selection is a function of the parameters being measured, this can bias the recovered dipole; the authors should either include the selection term in the likelihood or quantify how much the 6% uncertainty changes when selection is accounted for.","section":"Sec. 3.1, Eq. (3.4)"}],"minor_comments":[{"comment":"'import role' should be 'important role'.","section":"Abstract"},{"comment":"The network label '1ET+1CE+A#' in Table 3 is inconsistent with '1ET+1CE+1A#' used in Sec. 3.1 and Table 2.","section":"Table 3 vs Sec. 3.1"},{"comment":"'combbing' should be 'combining'.","section":"Conclusions"},{"comment":"The MCMC analysis does not state the priors on H0, g, l_dip, and b_dip; these should be reported for reproducibility.","section":"Sec. 3.2"},{"comment":"'differences in H0 measured from each golden dark sirens' should read 'differences in H0 measured from each golden dark siren'.","section":"Fig. 4 caption"},{"comment":"The same golden-dark-siren event set is reused when varying the injected g; because the selection of golden dark sirens depends on g through the localization area, the event set should ideally be regenerated for each injected value.","section":"Sec. 4, violin plots"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of JCAP, and the author is transparent about the main limitation, which is good practice. My main concern is that the headline constraints can be taken at face value only after the host-galaxy identification probability and peculiar-velocity scatter are quantified. The novelty relative to bright-siren forecasts is moderate, but the application of golden dark sirens to the cosmic dipole appears not to be duplicated elsewhere; I would support publication after the requested revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent forecast proposing a genuinely new use of golden dark sirens—measuring the cosmic dipole jointly with H0. The headline numbers are 6% on g=10^-3 and roughly 120 deg^2 on direction with 10 years of 1ET+2CEs, from 39 events. It also shows dipole effects on H0 from ordinary dark sirens are negligible for LVK O4/O5. That part is a useful sanity check.\n\nWhat's new: golden dark sirens have been used for H0, and bright sirens for the dipole, but this is the first to combine the two, including the rest-frame redshift conversion in eq. (3.5). The injection-recovery pipeline is standard and the posteriors in Figs. 5-7 look consistent with the stated uncertainties. The paper is also honest: it explicitly states the host identification probability is not calculated.\n\nThe soft spot, and it's a real one: the entire forecast assumes the brightest galaxy in the localization area is the host. The dipole signal is 10^-3 in D_L and z; a single misidentified host at z~0.1 can differ in redshift by far more than that. With 39 events, even a few wrong hosts could shift g by more than the claimed 6%. The paper concedes this in the conclusions but doesn't quantify it, and the quantitative claim shouldn't be taken at face value until that probability is computed. That is the main thing I'd want fixed before believing the number. There are also smaller issues: the likelihood in eq. (3.4) doesn't include selection effects for the 0.06 deg^2 localization cut, the merger rates for BBH/NSBH/BNS are assumed equal, and peculiar velocities are left out. None of these are fatal for a forecast, but they should be sensitivity-checked. No code or data are provided, so independent reproduction isn't possible.\n\nOverall: this is a solid, honest forecast with a load-bearing assumption that needs quantitative support. It deserves a serious referee. I'd send it to review, with the expectation that the host probability issue be addressed, at least with a toy galaxy catalog or a sensitivity analysis. It's a useful paper for people working on GW cosmology or the dipole tension.","headline":"A competent forecast proposing golden dark sirens for the cosmic dipole, with a load-bearing host-identification assumption that needs quantifying before the 6% number is trusted.","tokens_in":19877,"tokens_out":2674,"would_cite":true,"duration_ms":28453,"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":"The paper argues that golden dark sirens detected by next-generation ground-based networks can measure the cosmic dipole to order 10^-3 jointly with H0 and to order 10^-4 with H0 fixed.","keywords":["cosmic dipole","golden dark sirens","next-generation gravitational-wave detectors","Einstein Telescope","Cosmic Explorer","Hubble constant","dark siren cosmology","cosmological principle"],"falsifier":"Take a sample of well-localized nearby compact binary coalescences whose host galaxies are identified unambiguously by electromagnetic counterparts or deep spectroscopy, and check how often the brightest galaxy within 0.06 deg² is the true host; if the success rate is not close to unity, the forecasted $10^{-3}$ dipole constraints do not follow.","tokens_in":1908,"feed_emoji":"🌌","tokens_out":6918,"duration_ms":125956,"temperature":0.7,"pith_summary":"The paper proposes that golden dark sirens—nearby, tightly localized gravitational-wave events whose host galaxies can be identified by the single brightest galaxy rule—can measure the cosmic dipole in the era of next-generation ground-based detectors. It forecasts that a three-detector network containing at least one Einstein Telescope and one Cosmic Explorer would detect 35–39 golden dark sirens over 10 years, enough to constrain the dipole amplitude to order $10^{-3}$ (about 6% uncertainty) while jointly fitting H0, and to order $10^{-4}$ if H0 is fixed. The direction of the dipole would be localized to roughly 120 deg², comparable to or better than radio source number-count constraints. This matters because it offers a gravitational-wave probe of the reported ~4.9σ tension between the CMB dipole and sky-count dipoles, without relying on electromagnetic counterparts or number counting.","feed_headline":"Dark sirens could measure the cosmic dipole with 10^-4 precision","feed_subtitle":"Just 35-39 nearby golden dark sirens would independently test the 4.9σ cosmic dipole tension.","key_machinery":"Golden dark sirens. A golden dark siren is defined as a compact binary coalescence with z<0.1 whose 90% sky localization area is ≲0.06 deg², so that the Schechter-function galaxy density predicts only one galaxy brighter than L* in the area; that brightest galaxy is assigned as host and gives the observed redshift. The signal model is the dipole modification $D_L^{\\rm obs} = D_L^0\\,[1 + g(\\hat{n}\\cdot\\hat{z})]$ and $1+z_{\\rm obs} = (1+z_0)[1+g(\\hat{n}\\cdot\\hat{z})]$, with the rest-frame redshift recovered by inverting the second relation. The likelihood is a $\\chi^2$ over luminosity distances, and parameters $\\{H_0, g, l_{\\rm dip}, b_{\\rm dip}\\}$ are sampled with MCMC. Gravitational-wave parameter uncertainties, including the sky area used to select golden dark sirens, come from Fisher-matrix forecasts, and event simulation uses the GWTC-3 black-hole population and merger-rate model.","core_discovery":"In the paper's own terms, the central discovery is a forecast: golden dark sirens are rare but powerful, and next-generation detector networks are essential to find them. For an injected dipole of g=0.001 at the CMB direction, the 1ET+2CEs network yields g constrained to 6% and the direction within 120 deg² in a joint fit with H0, while 1ET+1CE+1A# yields 7% and a direction constrained to roughly 8°×6°. Fixing H0 lets the same dataset constrain g at order $10^{-4}$, although joint estimation with H0 degrades the dipole constraints. The paper also establishes that LVK-era dark sirens see no significant H0 bias from the dipole, because the number-count asymmetry and distance shifts roughly cancel in the all-sky combined posterior.","pith_inferences":["An extension the paper leaves implicit: the same χ² likelihood could be applied immediately to any handful of real well-localized dark sirens from a next-generation detector, with the host-assignment success rate as the controlling systematic.","If golden dark sirens later measure a dipole consistent with the CMB kinetic value while number-count analyses remain high, that would push the dipole tension toward source-evolution or selection systematics in the number counts rather than intrinsic anisotropy.","The brightest-galaxy assumption could be tested before ET/CE operate by applying the 0.06 deg² criterion to galaxy catalogues with complete spectroscopic host assignments, or by using the few well-localized LVK events with candidate hosts.","Combining golden dark sirens with bright sirens and number counting in one joint likelihood would likely sharpen both H0 and dipole constraints; the paper notes the event-number gain but does not compute the combined forecast."],"forward_implications":["If the forecast holds, a single 1ET+2CEs network can serve as an independent cosmic-dipole observatory, with no need for electromagnetic counterparts or the uncertain bright-siren rate.","The same 35–39 golden dark sirens would pin H0 to roughly 0.11–0.13% uncertainty, competitive with late-universe distance-ladder measurements and relevant to the Hubble tension.","Networks with only one next-generation detector (1ET+2A#s) detect too few golden dark sirens (13) to constrain the dipole, so at least two next-generation ground-based detectors are a prerequisite.","The dipole's effect on H0 from LVK O4/O5 dark sirens is negligible at g=0.01, but the paper argues that the number-count asymmetry grows with event rate and could bias dark-siren H0 in the next-generation era unless selection effects are modeled.","Golden dark sirens reach constraints comparable to bright-siren forecasts and to number counting with 10^7 events, while jointly measuring H0, a feature number counting does not provide."],"supporting_citations":[{"why":"Supplies the galaxy areal-density relation that defines golden dark sirens and justifies the single-brightest-galaxy host assignment.","marker":"[55]"},{"why":"Provides the sub-percent H0 measurement precision with ET+CE golden dark sirens that motivates the joint dipole-H0 forecast.","marker":"[60]"},{"why":"Bright-siren cosmic dipole forecast with next-generation networks, used as the main comparison for the golden dark siren constraints.","marker":"[48]"},{"why":"Number-counting dipole forecast with ET and CE, the alternative method against which the golden dark siren results are compared.","marker":"[46]"},{"why":"Number-counting dipole scaling with 10^7 events, used as the reference for how few golden dark sirens achieve comparable constraints.","marker":"[47]"},{"why":"Supplies the dipole modification of luminosity distance and redshift that forms the signal model for the analysis.","marker":"[62]"},{"why":"Provides the Power Law + Peak black-hole population and Madau merger-rate parameters used to simulate mock events.","marker":"[35]"},{"why":"Simulates detected gravitational-wave events for the different future detector networks.","marker":"[73]"},{"why":"Generates posterior samples for mock events via derivative approximation of the likelihood.","marker":"[74]"},{"why":"Provides the Fisher-matrix formalism used to forecast parameter uncertainties and select golden dark sirens by sky area.","marker":"[78]"}],"fun_headline_variants":["Golden dark sirens measure cosmic dipole to 10^-4","Next-gen detectors unlock dark siren dipole tests","Dark sirens could settle 4.9σ cosmic dipole clash","Golden dark sirens forecast 10^-4 dipole constraint","Rare golden sirens may probe cosmic anisotropy"],"cache_read_input_tokens":22016,"weakest_assumption_plain":"The argument stands or falls on the assumption that the single brightest galaxy inside a golden dark siren's localization area is its host galaxy, so that the measured galaxy redshift is the source redshift; the paper concedes this probability was not quantitatively calculated.","fun_headline_variants_meta":{"raw":{"variants":["Golden dark sirens measure cosmic dipole to 10^-4","Next-gen detectors unlock dark siren dipole tests","Dark sirens could settle 4.9σ cosmic dipole clash","Golden dark sirens forecast 10^-4 dipole constraint","Rare golden sirens may probe cosmic anisotropy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000701,"raw_usage":{"total_tokens":3198,"prompt_tokens":1014,"completion_tokens":2184,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":2104}},"tokens_in":630,"tokens_out":2184,"duration_ms":17291,"temperature":1.0,"reasoning_tokens":2104,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:28:49.854970+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a sample of well-localized nearby compact binary coalescences whose host galaxies are identified unambiguously by electromagnetic counterparts or deep spectroscopy, and check how often the brightest galaxy within 0.06 deg² is the true host; if the success rate is not close to unity, the forecasted $10^{-3}$ dipole constraints do not follow.","supporting_citations":[],"review_version":1}