Measuring and unbiasing the BAO shift in the Lyman-Alpha forest with AbacusSummit

The currently observing Dark Energy Spectroscopic Instrument (DESI) places sub-percent constraints on the Baryon Acoustic Oscillations (BAO) scaling parameters from the Lyman-$\alpha$ (Ly-$\alpha$) forest. However, no systematic error budget stemming from non-linearities in the 3D clustering of the Ly-$\alpha$ forest is included in the DESI-Ly-$\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly-$\alpha$ forest mocks produced on the $N$-body simulation suite \textsc{AbacusSummit}. Specifically, we measure the Ly-$\alpha$ auto-correlation and the Ly-$\alpha$-quasar cross-correlation functions. We use the DESI Ly-$\alpha$ forest fitting pipeline, \textsc{Vega}, with the publicly available covariance matrix from eBOSS DR16. To mitigate the noise, we adopt a linear control variates (LCV) technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the auto-correlation, we detect a small positive shift in radial direction of $\Delta \alpha_\parallel = 0.35\%$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha_\perp$. From the cross-correlation, we see a similar shift to $\Delta\alpha_\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha_\perp=\sim$0.25\%, at the 2$\sigma$ level. We also make a connection with the Ly-$\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly-$\alpha$ auto and the Ly-$\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly-$\alpha$ forest analyses based on EFT. This work paves the way for the full-shape analysis of DESI and future surveys.