Constraints on Primordial Oscillatory Features from the Power Spectrum and Bispectrum of Redshift-Space Galaxy Clustering in DESI DR1

We constrain oscillatory features in the primordial power spectrum using the public DESI DR1 galaxy and quasar clustering data. Our analysis is based on a custom pipeline for the redshift-space power spectrum and bispectrum of the six main DESI DR1 tracers, modeled with the Lagrangian effective field theory,which provides a principled resummation of the non-linear damping of the features. We consider oscillations that are linear and logarithmic in wavenumber, and, for the first time, consistently include the galaxy bispectrum in the search. We also study the dependence of the results on the background cosmology, which is important given the tension between the DESI BAO and Planck data within the baseline cosmological model. Marginalizing over Λ CDM parameters, we bound the amplitude of linear and logarithmic oscillations to A lin<0.033 and A log<0.035 respectively (at 95 % CL) across the frequency ranges ω lin phys in(20,800) Mpc and ω log in(2.5,80). The constraints reach percent-level precision away from the frequencies degenerate with the baryon acoustic oscillations (BAO). Fixing to the Planck 2018 cosmology, the bounds tighten to A lin<0.024 and A log<0.031, but in this case we also find a 3σ global preference for a logarithmic feature with ω log 15, dominated by the LRG2 sample, whose local wavelength matches the apparent BAO scale at wavenumbers best constrained by DESI. Since the preference for this feature disappears once the cosmological parameters are varied, we expect that it is a manifestation of the known distance-scale tension between the DESI and Planck data rather than a primordial signal. Our analysis thus emphasizes the importance of cosmology marginalization in searches for primordial features.