Anton Chudaykin, Mikhail M. Ivanov, Oliver H. E. Philcox

Models of cosmic inflation generically predict a weak but potentially detectable amount of primordial non-Gaussianity (PNG), which can be used to obtain insights into the degrees of freedom during inflation and their interactions. The simplest types of PNG are the local and non-local (equilateral and orthogonal) shapes of the primordial three-point correlators, which are predicted by models with multiple light fields and derivative interactions in single-field inflation, respectively. In this paper we place constraints on local, equilateral, and orthogonal non-Gaussianities using the power spectrum and bispectrum extracted from first public release of the Dark Energy Spectroscopic Instrument (DESI). Our analysis makes use of higher-order clustering information through a consistent effective field theory (EFT) model for both the power spectrum and bispectrum at one-loop order. Using robust scale cuts where the EFT description is valid, we find the following constraints on PNG amplitudes: $f^{\rm loc}_{\rm NL}=-0.1\pm 7.4$, $f^{\rm equil}_{\rm NL}=719\pm 390$, $f^{\rm orth}_{\rm NL}=-200\pm 100$ (at $68\%$ CL). Non-local PNG constraints can be further improved by combining high-redshift DESI with legacy BOSS data and using simulation-based priors on bias parameters, yielding the strongest large-scale structure constraints to date $f^{\rm equil}_{\rm NL}=200\pm 230$, $f^{\rm orth}_{\rm NL}=-24\pm 86$. Our constraint on $f^{\rm loc}_{\rm NL}$ is competitive with the cosmic microwave background (CMB) limit; the combination gives $f^{\rm loc}_{\rm NL}=-0.0\pm 4.1$, $18\%$ stronger than the CMB only result, which represents the strongest bound on multi-field inflation yet obtained.