{"id":"833fa151-d818-4327-9797-d1bb493656fa","arxiv_id":"2412.05990","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"DESI DR1 BAO data combined with Planck priors measure the phase shift amplitude beta_phi = 2.7^{+0.60}_{-0.67}, a 4.3 sigma detection of beta_phi > 0.","lead":"This paper measures the phase shift imprinted on baryon acoustic oscillations by free-streaming neutrinos, using DESI's first data release. It finds beta_phi = 2.7 ± 1.7 from BAO alone, and a 4.3 sigma detection of a non-zero phase shift when combined with Planck priors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 4.3 sigma significance is dominated by a Planck LCDM+Neff prior on alpha and alpha_AP, not by DESI BAO data alone; the paper's own model variations reduce it to 3.2-3.7 sigma.","rationale":"The central claim in the paper is the 4.3 sigma measurement of beta_phi > 0. The load-bearing condition for that claim is that the Planck prior applied to alpha and alpha_AP is both valid and fully validated, but this condition is the least secure part of the analysis. The paper demonstrates extensive validation of beta_phi recovery from mocks, and the two-code cross-checks are genuine independent support, but none of that validation exercises the prior-combination step that actually produces the headline result. The paper's own Table 8 shows that the significance and central value are sensitive to reasonable extensions of LCDM, which is direct evidence that the 4.3 sigma statement is not robust. The reader's formal weakest assumption concerned the fixed functional form f(k); that is a real secondary concern, but the prior-dependence of the significance is more directly load-bearing for the stated 4.3 sigma claim. Because the paper is transparent about this sensitivity and the DESI-only result is honestly reported, the conditional verdict remains appropriate; the requested change is to present the DESI-only constraint and the model-dependence of the prior-weighted result as the primary conclusions.","tokens_in":30853,"tokens_out":5749,"duration_ms":59659,"concrete_test":"Apply the exact Eq. 29 Planck importance-sampling step to the second-generation mocks with known beta_phi (c000: beta_phi = 1; c003: beta_phi approx 1.117) and compute the coverage of the quoted 68% intervals; if the recovered beta_phi is biased or coverage is not approximately 68%, the headline 4.3 sigma significance is not trustworthy. As a cross-check, quote the DESI-only beta_phi significance without any Planck weighting as the primary headline constraint.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline beta_phi = 2.7^{+0.60}_{-0.67} and the 4.3 sigma significance are not a DESI-only detection. The DESI BAO-only combined fit gives beta_phi = 2.70 +/- 1.70 (Table 8), which is only 1.6 sigma from zero. The quoted high significance comes entirely from importance-sampling the DESI posteriors with a Planck LCDM+Neff prior on alpha and alpha_AP (Eq. 29), shrinking the uncertainty by roughly a factor of three without moving the peak. That prior step is never validated on mocks: the mock suite validates beta_phi recovery with beta free, not with the Planck importance-sampling procedure used for the headline. Moreover, the prior is not model-independent. Table 8 shows the result moves to beta_phi = 2.05 +/- 0.55 (3.7 sigma) when A_lens is free, to 2.44 +/- 0.70 (3.4 sigma) in wCDM, and to 3.7 +/- 1.1 (3.2 sigma) in w0waCDM. Since beta_phi is degenerate with r_s (Figure 12, right), the Planck prior on alpha and alpha_AP, obtained from chains with Neff free, indirectly injects CMB information correlated with Neff into the beta_phi constraint despite the claim that no explicit CMB Neff information is used. The central claim that the phase shift is not purely sourced from the standard model expectation is therefore prior-dominated and model-dependent, and the paper's own model variations already weaken the significance to roughly 3.2-3.7 sigma.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the Baumann et al. (2019) phase-shift methodology to anisotropic BAO fits, implements it in two independent public codes (desilike and Barry), validates the pipeline on first- and second-generation DESI mocks, and applies it to DESI DR1 data. The DESI-only combined constraint is beta_phi = 2.70 ± 1.70, and after importance-sampling the BAO posteriors with a Planck-based prior on the alpha and alpha_AP distortion parameters the paper quotes beta_phi = 2.7^{+0.60}_{-0.67}, interpreted as beta_phi > 0 at 4.3 sigma. The paper also identifies a ~3-sigma systematic difference in beta_phi between the polynomial power-spectrum broadband method and the other broadband/clustering choices, and it reports that the significance drops to 3.2-3.7 sigma when Alens, w, or w0-wa are varied.","tokens_in":31253,"tokens_out":5335,"duration_ms":49252,"significance":"If the 4.3-sigma detection were robust, it would be an interesting hint of physics beyond the standard-model Neff value, with implications for neutrino physics and BAO analyses. The paper's strengths include a careful, two-code mock validation program, public code release, tests of the fiducial/template cosmology dependence, and a transparent presentation of the individual tracer fits. The DESI-only measurement of beta_phi is a useful new result. However, the headline significance is not a DESI-only detection: it is driven by the Planck importance-sampling prior, it is not validated on mocks, and it is sensitive to the assumed cosmological model. The central claim therefore needs substantial reframing and additional quantification before the paper can be accepted.","major_comments":[{"comment":"The headline 4.3-sigma significance is not a DESI-only measurement. As Table 8 shows, the DESI BAO-only combined fit gives beta_phi = 2.70 ± 1.70, i.e., only 1.6 sigma from zero. The tightened result beta_phi = 2.7^{+0.60}_{-0.67} comes entirely from importance-sampling the DESI posteriors with the Planck-based prior on alpha and alpha_AP in Eq. (29). Because that prior is constructed from Planck chains in which Neff is left free, and because alpha and alpha_AP are functions of r_s and the angular diameter distance (Eqs. 7 and 8), the prior implicitly carries Planck information on Neff through parameter correlations, despite the claim in §4.1 that 'we do not fold in any explicit CMB information on Neff.' This indirect injection should be quantified—for example by comparing with a prior derived from Planck chains with Neff fixed—and the importance-sampling procedure itself should be validated on the mock suite, since the current mock tests validate beta_phi recovery with beta free but not the Eq. (29) weighting step.","section":"§4.1, Table 8, Eq. (29)"},{"comment":"The paper reports a systematic difference in beta_phi of 0.29-0.40 (at 1.75-2.4 sigma, and up to 3 sigma when comparing the polynomial power-spectrum method to other methods) between broadband fitting methodologies, yet no systematic error is added to the quoted beta_phi uncertainty. The paper justifies this by noting the shift is smaller than the DESI DR1 statistical error and that the data results appear robust. However, the shift is comparable to the 0.6-0.7 uncertainty quoted for the Planck-prior combined result in Table 8, so it is not negligible for the central claim. The paper should either propagate this systematic into the final error budget or demonstrate quantitatively on the second-generation mocks and DR1 data that the final combined beta_phi constraint is insensitive to this choice for the specific pipelines used in the headline result.","section":"§3.1.4, Table 2"},{"comment":"The interpretation 'may hint at a phase shift that is not purely sourced from the standard model expectation for Neff' is model-dependent. Table 8 shows that the significance of beta_phi > 0 drops from 4.3 sigma in flat LCDM+Neff with w and Alens fixed to 3.7 sigma with Alens free, 3.4 sigma in wCDM, and 3.2 sigma in w0-waCDM, with the central value also moving from 2.70 to 2.05 in the Alens-free case. Since DESI DR1 itself moderately prefers the w0-waCDM model, the Planck-prior model choice is not neutral. The paper should present the DESI-only result as the primary measurement and clearly separate the model-dependent, prior-driven variants in the abstract, results, and conclusions, stating the full range of significances rather than only the 4.3-sigma value.","section":"§4.1, Table 8, Conclusions"},{"comment":"The analysis assumes the Baumann et al. functional form f(k) with fixed shape parameters (phi_inf, k*, xi). The mock validation only tests this shape for standard-model Neff variants (c000 and c003 cosmologies). If the true phase shift has a different k-dependence—for example from interacting neutrinos or non-adiabatic fluctuations, as the paper itself mentions in the conclusions—the fitted beta_phi will absorb the shape mismatch and the quoted amplitude constraint will be biased. The paper should make this limitation explicit in the abstract and, ideally, include a mock test with an alternative f(k) shape to quantify the bias in beta_phi.","section":"§2.3, Eqs. (17)-(18)"}],"minor_comments":[{"comment":"The abstract states the result 'relaxes in models with additional freedom beyond LCDM'; this caveat should be strengthened to explicitly state the model-dependent significance range (3.2-4.3 sigma) that is already present in Table 8, since the abstract's main numeric claim is the 4.3-sigma value.","section":"Abstract"},{"comment":"The caption 'This figure has been inspired by Figure 3 in Baumann et al. (2018)' is informal for a journal article; please rephrase to state that the figure reproduces or adapts the corresponding panel from that reference.","section":"§1 / Figure 1 caption"},{"comment":"The note that the alpha/alpha_AP prior is 'not the same as the prior included earlier from Planck on Neff' is confusing because the earlier validation in Section 4 also uses importance sampling with a Planck-based Neff prior. Clarify the distinction between the two uses of Planck chains.","section":"§4.1, footnote 13"},{"comment":"There are several typos and grammatical slips (e.g., 'used to to test' in §1, 'anistropic' in §2.2, and 'Baumman' in the Appendix A caption). A careful proofreading pass is recommended.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"This is a careful and useful analysis with strong mock validation, and the DESI-only beta_phi constraint is a legitimate contribution. The problem is that the paper's headline claim—the 4.3-sigma detection—is not supported by the DESI BAO data alone and is not robust to the modeling choices that the paper itself considers. I would urge the editor to require a substantial revision in which the primary result is the DESI-only measurement, the Planck-prior variants are presented as model-dependent illustrations with their full significance range, and the unquantified broadband systematic and the unvalidated importance-sampling step are addressed. The paper is not fatally flawed, but the current framing overstates the significance of the result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe key takeaway: this is a careful, well-validated application of Baumann et al. (2019) to DESI DR1, but the headline 4.3σ significance is prior-driven, not DESI-alone. DESI BAO alone gives β_phi = 2.7 ± 1.7 (1.6σ). The 4.3σ comes from importance-sampling with a Planck ΛCDM+Neff prior on α and α_AP, and the paper's own model variations drop it to 3.2–3.7σ. The abstract is honest about this, but the framing still overreaches.\n\nWhat's new: first DESI DR1 phase-shift measurement, an anisotropic extension of the phase-shift template, and solid validation with two independent codes and mocks, including c003 cosmologies. The paper also flags a ~3σ difference between polynomial and spline broadband fits to the power spectrum, yet quotes values without adding a systematic. That is the main soft spot: the number is broadband-choice-dependent at a level comparable to the error. The authors argue it's smaller than DR1's statistical error, which is fair, but for a 4.3σ claim you should fold it in.\n\nThe second soft spot is the prior step itself. Mock validation tests β_phi recovery with β free; it doesn't validate the Planck importance-sampling. And since β_phi is degenerate with r_s, a Planck prior on α/α_AP from chains with Neff free injects CMB information that isn't purely geometric. The paper is careful to state no explicit CMB Neff information is used, but the phrase \"may hint at phase shift beyond standard model\" is too strong given the model dependence.\n\nThe reader's conditional verdict is about right. The measurement is plausible, consistent with BOSS, and useful as a first DESI constraint. It's not a discovery. Who gets value: BAO methodologists and light relics people. It deserves a serious referee; the work is careful and the data products are public. The referee should push for a systematic error from broadband choices and clearer separation of DESI-only vs DESI+Planck results in abstract and conclusions. My recommendation: engage with it, send to review, expect revision.","headline":"Solid, honest first DESI DR1 phase-shift measurement, but the 4.3σ headline is prior-driven and drops to 3.2–3.7σ under the paper's own model variations.","tokens_in":31992,"tokens_out":3079,"would_cite":true,"duration_ms":27730,"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":"DESI's first-year BAO data, combined with Planck, constrain the phase-shift amplitude to β_φ = 2.7^{+0.60}_{−0.67}, hinting that the phase shift may exceed the standard-model expectation.","keywords":["baryon acoustic oscillations","phase shift","free-streaming neutrinos","effective number of relativistic species","DESI DR1","large-scale structure","cosmological parameters"],"falsifier":"Run the same DESI DR1 analysis on mocks that include a non-standard phase shift with a different scale dependence (e.g., from interacting neutrinos) and check whether the standard template recovers the input β_φ without bias. If β_φ moves toward one or the statistical significance drops under a flexible f(k), the 4.3σ result is an artifact of the assumed template rather than new physics.","tokens_in":30680,"feed_emoji":"🔭","tokens_out":6888,"duration_ms":54696,"temperature":0.7,"pith_summary":"This paper tries to measure the amplitude β_φ of the phase shift that free-streaming relativistic particles imprint on the baryon acoustic oscillation (BAO) wiggles, using DESI's first data release and an anisotropic BAO fitting pipeline. The central result is that the combined DESI BAO + Planck analysis gives β_φ = 2.$7^{{+0.60}}$_{−0.67}, suggesting β_φ > 0 at 4.3σ and larger than the standard-model value β_φ = 1 for N_eff = 3.044. If real, this would indicate a phase shift not purely sourced by standard-model neutrinos, possibly from non-standard neutrino physics or other light relics. The paper also validates the method against simulations and finds that the choice of broadband fitting method can shift β_φ by about 0.3 for power-spectrum fits, which must be controlled for future DESI data releases.","feed_headline":"DESI BAOs hint at a phase shift beyond standard-model neutrinos","feed_subtitle":"A 4.3-sigma signal in DESI BAO data plus Planck could mean new relativistic species or a statistical fluke.","key_machinery":"The load-bearing object is the parameterization of the BAO phase shift, φ(N_eff, k) = β_φ f(k), with f(k) = φ_∞ / [1 + (k_*/k)^ξ], where φ_∞ = 0.227, k_* = 0.0324 h/Mpc, ξ = 0.872, taken from Baumann et al. (2018, 2019). The paper extends this to an anisotropic BAO fit by inserting (β_φ − 1) f(k′)/r_s into the wiggle power-spectrum template alongside the anisotropic distortion parameters α_∥ and α_⊥, and implements this in two independent codes (Barry and desilike). The combination of tracers with importance sampling, plus a Planck prior on α and α_AP, is what sharpens the constraint from β_φ = 2.7 ± 1.7 to 2.$7^{{+0.60}}$_{−0.67}.","core_discovery":"The paper's central claim is that the phase-shift amplitude β_φ can be measured with DESI DR1 BAO data using an anisotropic BAO fitting pipeline, and that the combined measurement, after adding a Planck-based prior on the BAO distortion parameters, is β_φ = 2.$7^{{+0.60}}$_{−0.67}. This is a 4.3σ preference for β_φ > 0 and about 2.6σ from the standard-model value β_φ = 1 (N_eff = 3.044). The authors interpret this as a hint of a phase shift not purely sourced by the standard-model expectation for N_eff, while noting it could be an upward statistical fluctuation and that the tension relaxes when extra model freedom (e.g., wCDM or varying A_lens) is allowed.","pith_inferences":["A natural next test is to fit β_φ with a flexible functional form for f(k) (varying k_* and ξ), because non-standard neutrino interactions are expected to change the scale dependence; if the high β_φ persists under such a flexible model, the case for new physics is much stronger.","The paper's result rests on the Planck prior for α and α_AP; a cross-check with a prior from a different dataset (e.g., CMB lensing or supernovae) would show whether the high β_φ is driven by the specific Planck chains used.","If the phase shift is indeed larger than standard-model neutrinos produce, the same effect should appear in the CMB's acoustic peaks; comparing with phase-shift constraints from the Planck CMB (e.g., Montefalcone et al. 2025) could reveal whether the discrepancy is a BAO-specific systematic or a genuine cosmological signal."],"forward_implications":["If the 4.3σ detection is genuine, the phase shift in DESI BAOs is larger than the standard-model N_eff prediction, pointing to non-standard neutrino physics, non-adiabatic primordial fluctuations, or other free-streaming relics.","Allowing β_φ to vary weakens constraints on the BAO distance parameters α and α_AP because of the strong degeneracy between α and β_φ; future precise analyses may need to marginalize over β_φ to avoid biasing distance measurements.","The consistency between the two fitting codes and the mock validation suggest the measurement is robust at the current statistical precision, but the polynomial broadband method applied to the power spectrum produces a ~3σ shift in β_φ relative to other choices, so a systematic error budget is needed for DESI Y5.","The central value β_φ ≈ 2.7 maps to an unphysical N_eff if interpreted purely as a change in the number of neutrinos, meaning the standard-model interpretation is already strained; future data with σ(β) ~ 0.3 will decisively test whether the shift persists."],"supporting_citations":[{"why":"Establishes the physical mechanism by which free-streaming neutrinos induce a phase shift in BAOs.","marker":"Bashinsky & Seljak (2004)"},{"why":"Derives the fitting function f(k) and forecasts for galaxy surveys, providing the template used here.","marker":"Baumann et al. (2018)"},{"why":"Introduces the β_φ parameterization and the BOSS DR12 measurement this paper extends and compares with.","marker":"Baumann et al. (2019)"},{"why":"Supplies the BAO power-spectrum model (broadband and wiggle components) used in the template.","marker":"Beutler et al. (2017)"},{"why":"Provides the DESI DR1 BAO modeling, including spline broadband and reconstruction details, that the paper adopts.","marker":"Chen et al. (2024a)"},{"why":"Releases the DESI DR1 BAO data and standard results that this work applies the phase-shift method to.","marker":"Adame et al. (2024b)"},{"why":"The Planck 2018 chains used to build the prior on the BAO distortion parameters via importance sampling.","marker":"Aghanim et al. (2020)"},{"why":"Provides the Barry code, one of the two independent fitting pipelines used for cross-validation.","marker":"Hinton et al. (2020)"}],"fun_headline_variants":["DESI BAO + Planck hint at new relativistic species (4.3σ)","Phase shift in DESI BAO: beyond standard neutrinos?","4.3σ phase shift hint in DESI BAO + Planck","Possible new neutrino-like physics from DESI BAO"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fitting assumes that the true phase shift has exactly the scale dependence given by the template f(k) in equation 17; if non-standard physics changes this shape, the fitted β_φ will absorb that difference and the quoted significance could be biased.","fun_headline_variants_meta":{"raw":{"variants":["DESI BAO + Planck hint at new relativistic species (4.3σ)","Phase shift in DESI BAO: beyond standard neutrinos?","4.3σ phase shift hint in DESI BAO + Planck","Possible new neutrino-like physics from DESI BAO"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000579,"raw_usage":{"total_tokens":2784,"prompt_tokens":1058,"completion_tokens":1726,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":674,"completion_tokens_details":{"reasoning_tokens":1650}},"tokens_in":674,"tokens_out":1726,"duration_ms":10714,"temperature":1.0,"reasoning_tokens":1650,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:07:17.863467+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same DESI DR1 analysis on mocks that include a non-standard phase shift with a different scale dependence (e.g., from interacting neutrinos) and check whether the standard template recovers the input β_φ without bias. If β_φ moves toward one or the statistical significance drops under a flexible f(k), the 4.3σ result is an artifact of the assumed template rather than new physics.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the physical mechanism by which free-streaming neutrinos induce a phase shift in BAOs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the fitting function f(k) and forecasts for galaxy surveys, providing the template used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the β_φ parameterization and the BOSS DR12 measurement this paper extends and compares with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the BAO power-spectrum model (broadband and wiggle components) used in the template."},{"cited_title":"R., Howlett C., Davis T","cited_arxiv_id":null,"evidence_quote":"Provides the Barry code, one of the two independent fitting pipelines used for cross-validation."}],"review_version":1}