{"id":"2d123c53-8584-4a9c-8716-c72c7c1a5ab0","arxiv_id":"2506.08082","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A biased-tracer clustering model that adds scale-dependent density/velocity bias and mode coupling accurately describes BAO-scale redshift-space multipoles in simulations, and scale-dependent bias matters more than mode coupling.","lead":"This paper extends a model-agnostic description of the baryon acoustic oscillation (BAO) feature in galaxy clustering to include scale-dependent galaxy bias and mode coupling, and validates the extended model against N-body simulations. It matters because future surveys like DESI will need such corrections to extract unbiased cosmological constraints without assuming a specific cosmological model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Best-fit Bv = -14.3 makes the k^2 term of Eq. (2.16) a 22-89% correction over the k range dominating BAO-scale integrals, so the truncation of the scale-dependent bias expansion is uncontrolled and the claimed accuracy may not be a robust feature of the model.","rationale":"I reviewed the reader's weakest assumption, which targets the mode-coupling approximation in Eq. (2.17). That concern is real but, in my reading, less load-bearing than a more direct threat to the paper's main new effect. The mode-coupling contribution is explicitly small: the paper reports that the dotted curves with A_MC = 0 are very close to the solid best-fit curves (Fig. 1), and Appendix B tests the d xi/d ln s proportionality against an analytic example and finds it acceptable. The scale-dependent bias term, by contrast, is the dominant new ingredient: it drives the Delta chi^2 ~ 33 improvement over the no-sdbmc model. The best-fit Bv = -14.3 makes the supposedly small O(k^2 R_p^2) correction a 22-89% effect at k = 0.04-0.1 h/Mpc, where the BAO signal lives. This is not a small perturbative correction, so the k^2 truncation is uncontrolled unless the next order is shown to be negligible. The paper offers no such check and attributes the large, sign-flipped Bv to mass averaging without a demonstration, deferring it to a forthcoming work. If the k^4 terms are comparable, the fitted parameters could be absorbing truncation error rather than physical scale-dependent bias, and the model's apparent accuracy would not be a genuine validation. The proposed k^4 extension is a minimal, decisive test: it either confirms convergence or exposes the truncation as a source of bias. The reader's verdict is CONDITIONAL, and my concern supports keeping that verdict rather than changing it, so I mark the verdict as UNCHANGED.","tokens_in":19201,"tokens_out":18904,"duration_ms":230835,"concrete_test":"Extend Eq. (2.16) with k^4 terms: replace B1 k^2 R_p^2 with B1 k^2 R_p^2 + B2 k^4 R_p^4, and Bv k^2 R_p^2 with Bv k^2 R_p^2 + Bv2 k^4 R_p^4. Re-fit the same HADES configuration-space multipoles (55 <= s/(h^-1 Mpc) <= 125) using the same covariance and fixed b = 1.95, R* = 2.5. If the best-fit chi^2 improves by more than ~6 for the two extra parameters, or if B1 or Bv shift by more than their 68% errors, the k^2 truncation is not converged and the claimed accuracy of the sdbmc model is not a stable prediction of the model. If chi^2 and the parameters are unchanged, the truncation concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.2 models scale-dependent bias by truncating the peaks-theory form at O(k^2 R_p^2): Eq. (2.16) has B(k, mu) = [1 + B1 k^2 R_p^2 + beta mu^2 (1 - Bv k^2 R_p^2)] e^{-k^2 R*^2/2}. The HADES fit yields Bv = -14.3 with R_p = 2.5 h^-1 Mpc (Table 1). In the k range that dominates the BAO-scale configuration-space integrals (k ~ 0.04-0.1 h/Mpc), the correction (1 - Bv k^2 R_p^2) ranges from 1.22 to 1.89, i.e., a 22-89% admixture. The next-order k^4 R_p^4 terms are omitted, and no argument is given that their coefficients are small. The paper's justification for the unexpected negative Bv (mass averaging of the peaks model) is deferred to future work (Sec. 4.1: 'We will explore this using a detailed calculation in forthcoming work'). Since scale-dependent bias is the dominant new effect of the paper (Delta chi^2 ~ 33, with mode coupling contributing only a small fraction), an uncontrolled truncation in this term directly threatens the central claim that the model accurately captures the physical scale-dependent bias needed for DESI. The reader's focus on the mode-coupling ansatz is less load-bearing because the MC contribution is small and Appendix B explicitly tests the d xi/d ln s approximation.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the Zel'dovich-smearing model of PS23 by adding scale-dependent Lagrangian density and velocity bias plus a mode-coupling term, yielding the 'sdbmc' model for redshift-space multipoles of the 2-point correlation function. The model is validated against configuration-space multipoles from 20 HADES simulations at z=0 for a mass-threshold halo sample in the range 55-125 Mpc/h (chi2/dof = 118.7/98, p = 0.082) and against Fourier-space multipoles from the MINERVA simulations at z=0.57 (chi2/dof = 41.55/35). The authors find that scale-dependent bias is the dominant new ingredient (Delta chi2 ~ 33 relative to no-sdbmc), while mode coupling is subdominant but nonzero, and they argue the model is suitable for model-agnostic BAO inference in surveys like DESI.","tokens_in":19450,"tokens_out":6151,"duration_ms":80227,"significance":"If the model holds, it provides a practical, few-parameter extension of the Laplace-Gauss class of BAO templates that can incorporate scale-dependent bias and mode coupling without assuming a specific cosmology. The paper's strengths are its validation on two independent simulation suites, its systematic exploration of model variations, and the explicit appendix testing the mode-coupling approximation. The HADES constraint R_MC = 5.2 h^-1 Mpc being consistent with sigma_v = 6.0 h^-1 Mpc is a useful internal consistency check. The principal weakness is that the scale-dependent velocity-bias expansion is truncated at O(k^2) while the best-fit coefficient is so large that the truncation is not controlled in the k-range relevant for the BAO integrals; this limits the physical interpretation of the fit and its safe extrapolation to other tracers and redshifts.","major_comments":[{"comment":"The scale-dependent bias model truncates the peaks-theory-inspired expansion at O(k^2 R_p^2), but the best fit Bv = -14.3 with R_p = 2.5 h^-1 Mpc makes the factor (1 - Bv k^2 R_p^2) equal to roughly 1.14 at k = 0.04 h/Mpc and 1.89 at k = 0.1 h/Mpc, i.e. a 14-89% admixture of the next-order term in the k-range that dominates BAO-scale configuration-space integrals. The omitted O(k^4 R_p^4) terms are therefore not guaranteed small, and the paper provides no convergence test (e.g., a fit including a k^4 term) and no independent check of Bv. Since scale-dependent bias is the main new effect driving the improvement over the no-sdbmc model, the central claim that the model captures the physical scale-dependent bias relevant for DESI needs either a demonstration that higher-order terms are negligible or an explicit reframing of Bv as a purely empirical parameter whose extrapolation to other tracers is not yet justified.","section":"§2.2, Eq. (2.16); Table 1"},{"comment":"The mode-coupling ansatz is presented as following from a product of a logarithmic derivative with a relatively flat volume integral, but the numerical test in Appendix B (Figure 5) is performed for a lognormal Lagrangian bias model with b10 = 1 and does not cover the regime of the large negative Bv found in the HADES fit. Because the mode-coupling contribution to the final chi^2 is small, this does not invalidate the quality of the fit, but the claim that the treatment of mode coupling is 'rather general' is stronger than what is demonstrated. The authors should either soften that claim or test the approximation in the fitted parameter regime, especially if the model is to be used for DESI-scale inference.","section":"§2.3, Eq. (2.17); Appendix B"}],"minor_comments":[{"comment":"The phrase 'model coupling' appears where 'mode coupling' is meant; please correct this typo.","section":"§4.2.2"},{"comment":"Reference [7] has a malformed journal identifier ('a (p)' instead of a proper journal name); similar issues appear in other bibliographic entries (e.g., [11], [12]).","section":"References"},{"comment":"The covariance matrix is constructed by scaling the Gauss-Poisson covariance to match the diagonal simulation errors while keeping the correlation structure, and the model used in the covariance is itself updated from a preliminary fit. Please report how sensitive the reported chi^2/dof and parameter errors are to this scaling and to the choice of setting AMC = R_MC = 0 in the covariance calculation.","section":"§3.2"},{"comment":"The paper says the no-sdbmc model is excluded by Delta chi^2 ~ 33, but the no-sdbmc model has b fixed to the best-fit value from the sdbmc analysis; specifying how chi^2 changes when b is refit for the no-sdbmc model would make the comparison cleaner.","section":"§4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of JCAP and the empirical validation is solid. The main concern is the uncontrolled O(k^2) bias expansion associated with Bv = -14.3; this is fixable by adding a convergence test or reframing the model as empirical, so it does not warrant rejection. No citation or attribution concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does exactly what it claims: it extends the PS23 Zel'dovich-smearing framework with scale-dependent Lagrangian density and velocity bias plus a mode-coupling term, and validates the resulting sdbmc model against HADES configuration-space multipoles and MINERVA Fourier-space multipoles. The fit quality is genuinely good (118.7/98, p=0.082 for HADES; 41.55/35 for MINERVA), the model variations are sensible, and the authors are transparent about where their approximations live. The conclusion that scale-dependent bias matters more than mode coupling on BAO scales is well supported by their numbers: Delta chi^2 ~ 33 comes mostly from the bias terms, and the MC contribution is small. This is a real step forward for model-agnostic BAO analysis, not a repackaging.\n\nThat said, I think the reader's emphasis on the mode-coupling ansatz points at the wrong soft spot. The MC term is small, Appendix B explicitly tests the d xi/d ln s approximation, and the volume integral is shown to be reasonably flat. The load-bearing issue is the scale-dependent bias expansion itself. Eq. (2.16) truncates the peaks-theory form at O(k^2 R_p^2), and the best fit gives Bv = -14.3. Over the k range that dominates the BAO-scale integrals (k ~ 0.04-0.1 h/Mpc), the factor (1 - Bv k^2 R_p^2) is a 14-89% correction. At the top of that range the next-order k^4 R_p^4 term is not obviously small; no argument is given for its coefficient, and the authors' explanation for the sign of Bv is explicitly deferred to future work. Since scale-dependent bias is the main new physics of the paper, the truncation needs either a direct derivation or a consistency check (e.g., adding the k^4 term and showing the fit is stable). Without that, the claimed accuracy is a fit, not a robust prediction.\n\nThe circularity concern is real but not disqualifying. The parameters are fitted to the same data used for validation, and R_MC is free; the paper partially addresses this by testing model variations and finding R_MC ~ sigma_v in HADES. Still, the authors should either hold out a simulation or demonstrate stability under an extended bias expansion.\n\nBottom line: this is a serious paper by people who know the literature, and it deserves referee time. I would send it out, but with a clear request to justify or fix the scale-dependent bias truncation before acceptance. The mode-coupling details can be left mostly as is.","headline":"A solid, honest extension of PS23 that validates a new bias and mode-coupling model against two simulation suites, but the scale-dependent bias truncation is uncontrolled at the best-fit parameters and needs more justification before the model can be used for DESI-era inference.","tokens_in":20079,"tokens_out":1470,"would_cite":true,"duration_ms":19742,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Adding scale-dependent bias and mode coupling to the Zel'dovich-smearing model reproduces the BAO-scale redshift-space multipoles of halo clustering in 20 N-body simulations, with $\\chi^2/\\mathrm{dof} = 118.7/98$","keywords":["baryon acoustic oscillations","redshift-space clustering","scale-dependent bias","mode coupling","Zel'dovich approximation","peaks theory","two-point correlation function multipoles","model-agnostic BAO inference"],"falsifier":"A decisive check would be to measure the residual between the full nonlinear correlation function and the Zel'dovich-smeared propagator in an N-body simulation and test whether it is proportional to $d\\xi_L/d\\ln s$ with a single constant over $55 < s/(h^{-1}\\,\\mathrm{Mpc}) < 125$; if the ratio is scale-dependent, or if a quadrupole appears in the residual, the mode-coupling ansatz of equation (2.17) fails and the fitted $A_{\\rm MC}$ and $R_{\\rm MC}$ are absorbing other nonlinear effects.","tokens_in":18892,"feed_emoji":"🔭","tokens_out":15594,"duration_ms":159144,"temperature":0.7,"pith_summary":"The paper extends a model-agnostic (cosmology-independent), Zel'dovich-based description of the baryon acoustic oscillation (BAO) feature in redshift space to include two nonlinear effects: scale-dependent Lagrangian density and velocity bias, and mode coupling. It claims that the resulting “sdbmc” model accurately describes the monopole, quadrupole, and hexadecapole (the $\\ell=0,2,4$ angular moments) of the two-point correlation function at separations $55 < s/(h^{-1}\\,\\mathrm{Mpc}) < 125$ for a realistic halo sample in 20 N-body simulation realisations, with $\\chi^2/\\mathrm{dof} = 118.7/98$. The practical motivation is that ongoing spectroscopic surveys will need these effects modelled if BAO-based cosmological inference is to remain unbiased. The paper also finds that scale-dependent bias matters more than mode coupling near the BAO scale, and that the mode-coupling smearing scale can be set equal to the linear velocity dispersion $\\sigma_v$.","feed_headline":"Adding scale-dependent bias fits BAO halo clustering in simulations","feed_subtitle":"It fits the monopole, quadrupole, and hexadecapole from 55 to 125 Mpc/h with p = 0.082.","key_machinery":"The load-bearing object is the ansatz of equation (2.17), $\\xi_{\\rm NL}(s) = \\int d^3s'\\,\\xi_L(s'|0)\\,\\mathcal{N}(s-s';\\Sigma) + A_{\\rm MC}\\,\\partial\\xi_L(s|R_{\\rm MC})/\\partial\\ln s$, in which the first term is the Zel'dovich-smeared propagator with an anisotropic Gaussian kernel of variances $2\\sigma_v^2(1+f)^2$ along the line of sight and $2\\sigma_v^2$ across it, and the second term is the mode-coupling piece written as a constant times the logarithmic derivative of a separately smeared linear correlation function. Scale-dependent bias enters through the factor $B(k,\\mu_k)$ built from $b_{\\rm Lag}(k) = (b_{-1}+b_{01}k^2R_p^2)e^{-k^2R_*^2/2}$ and $b_{\\rm vel}(k) = (1-B_v k^2 R_p^2)e^{-k^2 R_*^2/2}$, motivated by peaks theory. These are Laplace-Gauss expansions, meaning sums of powers of $k^2$ multiplying Gaussian factors $e^{-k^2\\sigma^2}$, which is what makes the model straightforward to incorporate in the existing model-agnostic BAO inference framework.","core_discovery":"The central claim is that a simple extension of the earlier Zel'dovich-smearing model --- equations (2.14) and (2.15) for scale-dependent Lagrangian density and velocity bias plus a mode-coupling term proportional to $d\\xi_L/d\\ln s$ --- captures the BAO-scale multipoles of biased tracers in redshift space. Fitting four parameters ($B_1$, $B_v$, $A_{\\rm MC}$, $R_{\\rm MC}$) with the large-scale bias and the smoothing scale fixed, the model matches the configuration-space multipoles measured in the N-body simulations with $\\chi^2/\\mathrm{dof} = 118.7/98$ ($p = 0.082$), and it also describes Fourier-space multipoles below $k = 0.13\\,h\\,\\mathrm{Mpc}^{-1}$ in a second simulation suite. The simpler no-sdbmc model is disfavored by $\\Delta\\chi^2 \\sim 33$, with $A_{\\rm MC}=0$ excluded at better than 95% confidence and $B_v=0$ at better than 99%. The model predicts a narrower BAO peak shifted to smaller separations and a larger quadrupole variation than the no-sdbmc case.","pith_inferences":["A direct follow-up test would be to fit the same model to narrow halo-mass bins: the significantly negative $B_v$ found for the mass-thresholded sample suggests that mass averaging reverses the sign of the velocity-bias coefficient relative to single-mass peaks-theory predictions, and narrow bins would show whether the scale dependence is physical or an absorbing nuisance.","If the equality $R_{\\rm MC} = \\sigma_v$ holds beyond the tested cases, the mode-coupling term is not an independent physical ingredient but a derivative of the propagator, so BAO analyses could drop it as a nuisance parameter and absorb it into the smearing kernel.","The Galilean-invariance argument cited for the dominant dipole mode-coupling term points to a possible gravity test: theories that preserve the relation between displacement and density would keep this term unchanged, so the shape of the BAO feature could distinguish modified gravity from the standard model even where isotropic smearing looks the same; the paper does not pursue this."],"forward_implications":["Ongoing spectroscopic surveys will need to include scale-dependent density and velocity bias in BAO modelling; the paper argues that ignoring it degrades the fit and biases the recovered large-scale bias $b$ and, through degeneracies, parameters such as the growth rate $f$ and $\\sigma_8$.","The mode-coupling smearing scale $R_{\\rm MC}$ can be fixed to $\\sigma_v$, reducing the number of free parameters in future model-agnostic fits.","Mode coupling contributes only a small fraction of the difference from the simpler model, so retaining a positive-amplitude term of this form matters more than modelling its details.","The model remains accurate in Fourier space below $k = 0.13\\,h\\,\\mathrm{Mpc}^{-1}$, providing a complementary regime for testing the same physics.","Because the model's predictions can be computed for arbitrary linear power spectra, it is agnostic to the cosmological model and can be applied beyond the reference cosmology used for validation."],"supporting_citations":[{"why":"Base model-agnostic Zel'dovich-smearing framework that this paper extends; its no-sdbmc variant is the comparison model.","marker":"[24]"},{"why":"Establishes the propagator picture in which the BAO feature is smeared by a Gaussian kernel sourced by bulk flows.","marker":"[26]"},{"why":"Source of the mode-coupling approximation as a term proportional to the logarithmic derivative of the linear correlation function, generalised here to biased, redshift-space tracers.","marker":"[33]"},{"why":"Provides the peaks-theory treatment of scale-dependent Lagrangian bias and its relevance on BAO scales.","marker":"[30]"},{"why":"Documents the scale-dependent velocity bias of halos that motivates the velocity-bias parametrisation.","marker":"[38]"},{"why":"Supplies the configuration-space multipole measurements from the N-body simulations used for the main fit.","marker":"[44]"},{"why":"Supplies the Fourier-space multipole measurements and the Gauss-Poisson covariance construction used in the Fourier-space test.","marker":"[45]"},{"why":"Describes the N-body simulations whose halo catalogues define the main test sample.","marker":"[43]"},{"why":"Constrains the leading mode-coupling form via Galilean invariance, supporting the claim that a dipole-like term dominates at BAO scales.","marker":"[31]"}],"fun_headline_variants":["Scale-dependent bias sharpens BAO clustering fits","New model handles BAO scale bias and mode coupling","BAO clustering: bias and mode coupling matter","Model includes scale-dependent bias for BAO surveys","Peaks theory meets Zel'dovich for BAO fits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the mode-coupling contribution is proportional to the logarithmic derivative of the linear correlation function with a constant amplitude and an isotropic Gaussian smearing scale; if the true volume-integral prefactor is not flat across the BAO scale, or if anisotropic smearing matters, the fitted $A_{\\rm MC}$ and $R_{\\rm MC}$ would absorb the mismatch and the claimed accuracy would not be a genuine test of the model.","fun_headline_variants_meta":{"raw":{"variants":["Scale-dependent bias sharpens BAO clustering fits","New model handles BAO scale bias and mode coupling","BAO clustering: bias and mode coupling matter","Model includes scale-dependent bias for BAO surveys","Peaks theory meets Zel'dovich for BAO fits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000542,"raw_usage":{"total_tokens":2641,"prompt_tokens":1032,"completion_tokens":1609,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":648,"completion_tokens_details":{"reasoning_tokens":1535}},"tokens_in":648,"tokens_out":1609,"duration_ms":13032,"temperature":1.0,"reasoning_tokens":1535,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:18:36.782716+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be to measure the residual between the full nonlinear correlation function and the Zel'dovich-smeared propagator in an N-body simulation and test whether it is proportional to $d\\xi_L/d\\ln s$ with a single constant over $55 < s/(h^{-1}\\,\\mathrm{Mpc}) < 125$; if the ratio is scale-dependent, or if a quadrupole appears in the residual, the mode-coupling ansatz of equation (2.17) fails and the fitted $A_{\\rm MC}$ and $R_{\\rm MC}$ are absorbing other nonlinear effects.","supporting_citations":[{"cited_title":"Velocity bias in the distribution of dark matter halos","cited_arxiv_id":"1405.5885","evidence_quote":"Documents the scale-dependent velocity bias of halos that motivates the velocity-bias parametrisation."}],"review_version":1}