{"id":"c806db26-f35e-4804-9001-3224c0445510","arxiv_id":"2607.18384","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A template-based bispectrum-monopole BAO measurement is shown to be unbiased on simulations and, jointly with the power spectrum, to sharpen BAO constraints while exposing baryon-dark matter velocity biases.","lead":"This paper builds a new model for how the baryon-dark matter relative velocity distorts the galaxy bispectrum in redshift space, and introduces a template method to measure the BAO scale from the bispectrum monopole. It argues the method is unbiased on N-body mocks and that combining it with the power spectrum could detect these velocity biases, not just improve BAO constraints.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unvalidated velocity-bias bispectrum model underlies all inferred Δα_iso shifts and the P/B diagnostic; if wrong, the central detectability claim would not hold.","rationale":"The paper has two distinct achievements. The BAO extraction from the bispectrum monopole is validated on 15000 Quijote mocks: the model fits the bispectrum shape within 5σ and recovers unbiased α_iso with errors comparable to pre-reconstruction P02. This part is credible and should be acknowledged. The relative-velocity diagnostic, however, is a forecast built atop the new velocity model in Eqs. 2.7–2.9. Since the synthetic data are generated with the same model, the reported shifts in Fig. 4 are purely a consequence of the model assumptions. The only way to confirm the diagnostic is to test the model against simulations that include v_bc; the paper does not present such a test. This is the weakest link because if the velocity terms are mis-modeled—for example, if the phase of k T_v(k) relative to P_BAO is wrong, or if RSD contributions are incomplete—the predicted Δα_iso differences between P02 and B0 would be incorrect, and the claimed ability to constrain b_v2/b_delta_bc would not hold. The internal inconsistencies (abstract vs. text numbers) suggest some quantitative results are unstable. Therefore, I concur with the reader's CONDITIONAL verdict: the BAO extraction method is acceptable, but the velocity diagnostic requires validation. The concrete test proposed would settle the concern.","tokens_in":24548,"tokens_out":7005,"duration_ms":53421,"concrete_test":"Run N-body simulations that include the baryon–dark matter relative velocity (e.g., by modulating initial conditions with the v_bc field as in Blazek et al. 2016 or using separate baryon/CDM particles), measure the redshift-space galaxy bispectrum monopole at z=0.5 for a few values of b_v2, b_delta_bc, b_theta_bc, and compare to the model of Eqs. 2.7–2.9. Then fit the BAO pipeline of Section 3.2 to the mock monopole and check whether the recovered Δα_iso matches the predictions of Figure 4 within statistical error. If not, the diagnostic claim is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the bispectrum monopole BAO can serve as a diagnostic for relative-velocity bias rests on the redshift-space tree-level bispectrum model of Section 2.2 (Eqs. 2.7–2.9). This model extends Yoo et al. (2011) by including b_v2, b_delta_bc, b_theta_bc terms in Z1 and Z2, but it is never validated against N-body simulations that include the streaming velocity. The Quijote validation in Appendices B–C uses mocks with no velocity bias, so it tests the BAO extraction technique, not the velocity terms. The synthetic forecasts in Section 2.3 and Figure 4 are generated from the same model, making the reported Δα_iso shifts self-consistency checks rather than independent predictions. If the velocity kernels or their phase shifts are incorrect (e.g., missing terms, wrong RSD treatment), the predicted shifts and the claimed sensitivity to b_v2/b_delta_bc would fail to reproduce on real data. The internal numerical inconsistencies (abstract 20% vs. text 10% for b_delta_bc; 30% vs. 22% for constraining power) further indicate the quantitative claims are not firmly grounded.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops and tests a method to extract the isotropic BAO dilation parameter alpha_iso from the monopole of the galaxy bispectrum, using a template-based model that embeds the linear BAO wiggle template in the P(k1)P(k2) prefactor. It extends the redshift-space tree-level bispectrum model to include the baryon-dark matter relative-velocity bias terms b_v2, b_delta_bc, and b_theta_bc. The BAO extraction is validated against Quijote N-body mocks: the bispectrum-monopole fit returns an unbiased alpha_iso with constraining power comparable to the pre-reconstruction power spectrum, and a joint P+B fit tightens the constraint. The paper then uses noiseless synthetic data generated from its own velocity-bias model to show that the relative-velocity terms shift the recovered alpha_iso differently for power spectrum and bispectrum, and argues that the P-B difference can be used to detect and constrain b_v2 and b_delta_bc.","tokens_in":24792,"tokens_out":3794,"duration_ms":34934,"significance":"If correct, the paper would establish the bispectrum monopole as a practical, unbiased BAO ruler and as a diagnostic for relative-velocity systematics in current and future surveys. The validation of the BAO extraction is a genuine strength: it uses 15000 Quijote realizations, reports fits to 500 independent realizations, shows residuals mostly within 3 sigma, and gives quantitative error estimates in Table 1. The updated redshift-space bispectrum model with all relative-velocity terms is also a useful theoretical contribution. However, the central detectability claim rests on a velocity-bias model that is not validated against N-body simulations containing streaming velocities; the synthetic forecasts are self-consistency checks rather than independent predictions. The quantitative claims also contain internal inconsistencies that must be resolved before the results can be used reliably.","major_comments":[{"comment":"The predicted Delta_alpha_iso shifts and the proposed P-B detectability of b_v2 and b_delta_bc are generated by constructing noiseless synthetic data with the paper's own model (Eqs. 2.7-2.9 and Appendix A) and then fitting with the same model's velocity terms set to zero. The Quijote validation in Appendices B-C uses mocks without streaming velocity, so it validates only the BAO extraction technique, not the velocity kernels. If the Z1/Z2 velocity terms or their oscillatory phases are incorrect, the central detection claim collapses. This is a load-bearing issue: either validate the velocity-bias model against N-body simulations that include the baryon-dark matter relative velocity, or explicitly reframe the forecast as a model-level self-consistency check and soften the detection claim.","section":"Section 2.3 and Section 4, Figs. 4-5"},{"comment":"The abstract states that systematic discrepancies reach 'up to 20% for b_delta_bc <= -2'. Section 4, however, states that for b_delta_bc < -2 the bispectrum fit rapidly degrades, that a shift of 10% on alpha_iso is reached for b_delta_bc < -5, and that the model cannot capture the bispectrum shape for extreme values. The conclusion again says 'up to a 20% difference'. These numbers and the associated parameter ranges must be reconciled, and the claim should be quoted only in the regime where the model is under control.","section":"Abstract vs. Section 4"},{"comment":"The claimed improvement in constraining power is reported inconsistently: the abstract and Section 3.2 say ~30%, while Appendix C and Table 1 show an improvement from 2.80% to 2.17%, i.e. 22%. Additionally, Appendix C's statement that the post-reconstruction measurement is 'about 30% more than the joint analysis' is ambiguous. The authors should use a single, precisely defined metric (e.g., ratio of standard deviations or variances) and report consistent numbers throughout.","section":"Section 3.2, Appendix C, Table 1"},{"comment":"The 'prescription to detect and constrain' the velocity-bias parameters is based on polynomial/sigmoid fits to noiseless model points, with error bars from a 500 (h^-1Gpc)^3 volume. No actual likelihood or expected-constraint calculation is presented for realistic survey volumes; the text itself notes that for a DESI-like volume of ~50 (h^-1Gpc)^3 the b_theta_bc difference is hidden in the statistical error. The claims of 'high sensitivity' and 'competitive constraints' are therefore not quantitatively demonstrated and should be backed by a forecast, e.g., a Fisher or MCMC analysis on the P-B difference statistic.","section":"Section 4 and Figure 5"}],"minor_comments":[{"comment":"Typo in the conclusion: 'b δbv >= 2' should presumably be 'b_delta_bc'. Please also standardize notation between b_delta_bc/b_theta_bc and b_bc^delta/b_bc^theta.","section":"Abstract/Conclusion"},{"comment":"The description of the triangle ordering says the x-axis is sorted ascending by k3, then k2, then k1, but Appendix C says triangles are ordered by ascending k1. Please clarify the convention.","section":"Section 2.3"},{"comment":"The central panel's error bars are stated to be invisible due to the plot scale. This makes it difficult to assess the claimed sensitivity; consider plotting residuals or a separate panel with zoomed range.","section":"Figure 5"},{"comment":"Typo: 'ower spectrum' should be 'power spectrum'. Also Table 1 label appears as 'T able 1'.","section":"Appendix C caption"},{"comment":"Reference [23] duplicates [18]. Also 'commoving' should be 'comoving' in the introduction.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The BAO-extraction methodology itself appears solid and well validated; the main weakness is that the relative-velocity detectability claims are model self-consistency forecasts with no independent N-body validation, plus unresolved numerical inconsistencies. These are fixable within the paper's scope by adding a validation test (or explicitly downgrading the claims) and by reconciling the reported numbers. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things before reading this one. First, the B0 BAO extraction is the real contribution: the authors fit the bispectrum monopole with a template model, validate it on 500 Quijote realizations, and show it gives an unbiased alpha_iso with errors comparable to the pre-reconstruction P02 and about 22% tightening in a joint fit. That part is honest and reproducible. Second, the relative-velocity part is a forecast, not a measurement, and the paper's abstract and conclusions do not always make that distinction.\n\nWhat is actually new: they extend the redshift-space tree-level bispectrum to include the b_v2, b_delta_bc, and b_theta_bc relative-velocity bias terms — the b_delta_bc and b_theta_bc terms were not in Yoo et al. (2011), so this is a real extension. They also provide a covariance-validated template for alpha_iso from the monopole alone, which Behera et al. had done without the covariance. The idea that the power-spectrum and bispectrum alpha_iso shifts differ in sign for b_v2 and negative b_delta_bc is useful intuition.\n\nThe soft spots are in the velocity part. The model is never checked against N-body simulations that actually include the streaming velocity; the Quijote validation only tests the BAO template, not the velocity kernels. The synthetic 'measurements' in section 2.3 are generated with the same model that is then fit, so the reported shifts and the P-B 'diagnostic' are self-consistency checks. If the velocity kernels or their phase shifts are wrong, the detectability claim collapses. That does not invalidate the B0 method, but it should be said plainly. There are also numeric inconsistencies between the abstract and the text — 20% vs 10% for b_delta_bc, 30% vs 22% for the joint improvement — and a post-hoc restriction of the b_delta_bc range after the wide range produced unphysical fits. All fixable, but they suggest the quantitative claims were not fully settled.\n\nWho it's for: BAO practitioners and anyone building bispectrum pipelines. The B0 extraction deserves a serious referee; the velocity forecast is worth one too, with the expectation that the model be validated on streaming-velocity simulations or at least clearly labeled as a theory prediction. I'd send it to review.","headline":"The bispectrum-monopole BAO extraction is a genuinely useful and well-validated method; the relative-velocity 'diagnostic' is a plausible but unvalidated forecast that needs clearer framing and some numeric cleanup.","tokens_in":25380,"tokens_out":3931,"would_cite":true,"duration_ms":33830,"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 claims that the bispectrum monopole gives an unbiased measurement of the BAO scale, that adding it to the power spectrum tightens constraints by about 30%, and that the difference between the two measurements can reveal the baryon","keywords":["baryon acoustic oscillations","bispectrum monopole","relative velocity effect","streaming velocity","BAO dilation parameter","galaxy bias","redshift-space distortions","large-scale structure"],"falsifier":"Generate a suite of N-body simulations that physically include the baryon–dark matter relative velocity, or use a galaxy sample where streaming velocities are known to matter, then run the same power-spectrum and bispectrum BAO fits. If the bispectrum monopole does not show the model's predicted phase-shifted oscillations at BAO scales, or if the measured alpha_iso difference between the power spectrum and bispectrum does not track b_v2 and b_delta^bc as predicted, the velocity-bias model is wrong even though the BAO extraction itself might remain useful.","tokens_in":24310,"feed_emoji":"📏","tokens_out":8024,"duration_ms":61808,"temperature":0.7,"pith_summary":"This paper tries to establish two things. First, the baryon acoustic oscillation (BAO) scale can be extracted from the monopole of the galaxy bispectrum alone, using a template that suffers no bias and gives constraints comparable to the pre-reconstruction power spectrum; fitting it together with the power spectrum tightens the error on the BAO dilation parameter by roughly 30%. Second, the same three-point measurement, compared with the power spectrum, can detect the cosmological relative velocity between baryons and dark matter: the three bias parameters that describe this effect shift the recovered BAO scale differently in the two probes, up to 2% for b_v2 = ±0.05 and up to 20% for b_delta^bc ≤ −2. This matters because standard BAO analyses ignore this effect, so if it is present at even a low level it biases the distance scale; the bispectrum provides a way to see and remove that bias.","feed_headline":"Three-point statistic recovers the BAO ruler without bias","feed_subtitle":"Adding the bispectrum monopole tightens the distance-scale error by ~30% and exposes baryon–dark matter velocity bias.","key_machinery":"The central object is the redshift-space tree-level bispectrum monopole B^(0)(k1,k2,k3; alpha_iso), built from the linear and second-order kernels Z1 and Z2 that are extended to include the relative-velocity bias terms, paired with an isotropic BAO template in which each power-spectrum factor is split into a smooth broadband plus a damped BAO wiggle O_lin(k/alpha_iso). The relative-velocity terms enter through the transfer functions T_bc(k) and T_v(k), whose oscillations are phase-shifted relative to the BAO wiggle; that phase shift is what converts a real velocity bias into an alpha_iso shift, and it is the property that lets the power-spectrum and bispectrum measurements disagree in a diag","core_discovery":"The paper's central claim is that the redshift-space tree-level bispectrum, extended to include all relative-velocity terms (the b_v2, b_delta^bc, and b_theta^bc biases and their redshift-space counterparts), carries a BAO signal that can be isolated with a template-based fit to the monopole. The template writes the bispectrum in terms of power spectra and effective second-order kernels, with the BAO wiggle entering through P(k; alpha_iso). Validated on N-body simulations that do not include the streaming-velocity bias, the method returns alpha_iso unbiased and with precision comparable to the pre-reconstruction power spectrum; in combination with the power spectrum, the error shrinks by abo","pith_inferences":["If the velocity-bias model holds, the same power-spectrum versus bispectrum alpha_iso comparison is a generic null test for any mechanism that imprints phase-shifted oscillations on the BAO scale, not just streaming velocities—dark-matter oscillations and isocurvature perturbations would produce analogous discrepancies.","For surveys where post-reconstruction is unavailable or unreliable, the bispectrum monopole could nearly substitute for the information reconstruction provides, since the joint pre-reconstruction fit closes much of the gap to post-reconstruction precision.","The steep divergence of the bispectrum alpha_iso for negative b_delta^bc suggests that even a single BAO measurement on real data could act as a strong prior on this bias parameter, potentially sharpening full-shape analyses that currently constrain it only weakly.","A testable extension is to apply the same alpha_iso-difference diagnostic to the anisotropic dilation parameter once higher-order bispectrum multipoles become cheap to measure; the velocity terms should shift those measurements differently as well."],"forward_implications":["A BAO distance measurement can be obtained from the bispectrum monopole alone, with no reconstruction step and no bias, at precision comparable to the standard pre-reconstruction power spectrum.","Jointly fitting the power-spectrum monopole plus quadrupole and the bispectrum monopole improves the alpha_iso constraint by roughly 30% relative to the pre-reconstruction power spectrum.","A measured difference between alpha_iso from the power spectrum and from the bispectrum is a signature of relative-velocity bias; the expected size is about 2% for |b_v2| = 0.05 and up to 20% for b_delta^bc ≤ −2.","The b_theta^bc parameter is harder to isolate because it moves the power-spectrum and bispectrum measurements in the same direction, so velocity-divergence bias is the least constrained of the three.","Standard BAO pipelines that do not model streaming velocities can be biased at the roughly 1% level, comparable to the error budget of next-generation surveys; the bispectrum measurement offers a cross-check."],"fun_headline_variants":["Bispectrum BAO: unbiased ruler, 30% better with PS","Bispectrum BAO exposes baryon–dark matter velocity bias","Add bispectrum BAO to power spectrum, cut error 30%","Three-point BAO: unbiased acoustic scale, probes velocity bias","Bispectrum BAO: 30% tighter distances with power spectrum"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The method's central bet is that the equations connecting galaxy clustering to the baryon–dark matter streaming velocity are correct; these equations were checked against simulations without that velocity, so a wrong phase or amplitude in the velocity terms would wipe out the predicted shifts and the claimed sensitivity to b_v2 and b_delta^bc.","fun_headline_variants_meta":{"raw":{"variants":["Bispectrum BAO: unbiased ruler, 30% better with PS","Bispectrum BAO exposes baryon–dark matter velocity bias","Add bispectrum BAO to power spectrum, cut error 30%","Three-point BAO: unbiased acoustic scale, probes velocity bias","Bispectrum BAO: 30% tighter distances with power spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000963,"raw_usage":{"total_tokens":3992,"prompt_tokens":855,"completion_tokens":3137,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":3039}},"tokens_in":599,"tokens_out":3137,"duration_ms":19410,"temperature":1.0,"reasoning_tokens":3039,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T15:33:35.250495+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Generate a suite of N-body simulations that physically include the baryon–dark matter relative velocity, or use a galaxy sample where streaming velocities are known to matter, then run the same power-spectrum and bispectrum BAO fits. If the bispectrum monopole does not show the model's predicted phase-shifted oscillations at BAO scales, or if the measured alpha_iso difference between the power spectrum and bispectrum does not track b_v2 and b_delta^bc as predicted, the velocity-bias model is wrong even though the BAO extraction itself might remain useful.","supporting_citations":[],"review_version":1}