Farshad Kamalinejad, Zachary Slepian, Alex Krolewski, Alessandro Greco, William Ortolá Leonard, Jessica Chellino, Matthew Reinhard, Elena Fernández-García, Francisco Prada, J. Aguilar, S. Ahlen, A. Anand, C. Bebek, D. Bianchi, D. Brooks, T. Claybaugh, A. Cuceu, K. S. Dawson, A. de la Macorra, R. Demina, P. Doel, J. Edelstein, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez, H. K. Herrera-Alcantar, K. Honscheid, C. Howlett, D. Huterer, M. Ishak, R. Joyce, S. Juneau, D. Kirkby, T. Kisner, A. Kremin, O. Lahav, C. Lamman, M. Landriau, L. Le Guillou, M. Manera, A. Meisner, R. Miquel, J. A. Newman, W. J. Percival, C. Poppett, I. Pérez-Ràfols, L. Samushia, E. Sanchez, D. Schlegel, M. Schubnell, H. Seo, J. Silber, D. Sprayberry, G. Tarlé, B. A. Weaver, C. Zhao, H. Zou

We present the first detection of the 3-Point Correlation Function (3PCF) Baryon Acoustic Oscillation (BAO) signal from the DESI Data Release 1 (DR1) sample of Luminous Red Galaxies (LRGs), which contains over 2.1 million galaxies. Our analysis is based on a tree-level redshift-space bispectrum template, which is then transformed to position space using the Fast Fourier Transform on Logarithmic scales (FFTLog) algorithm. We detect the BAO feature with a significance of approximately $8.1σ$ using the EZmock covariance matrix and $8.5σ$ using the analytical covariance matrix, for the full LRG redshift range ($0.4<z<1.1$), denoted as the $z_{\rm full}$ sample. We use the Abacus altMTL mocks, the most precise DESI DR1 mock catalogs currently available, to validate our model. We find that our model fits the mocks well, with a small offset of $0.6\%$ in the recovered BAO scale, which we treat as a systematic error due to modeling. We measure the angle-averaged distance, $D_{\rm V}(z = 0.68)/r_{\rm d} = 15.88 \pm 0.27$ ($1.72\%$ precision) when using the covariance matrix estimated from EZmocks and $D_{\rm V}(z = 0.68)/r_{\rm d} = 15.72 \pm 0.18$ ($1.12\%$ precision) when using the analytical Gaussian covariance matrix. Our results show excellent agreement with the DESI DR1 2PCF BAO measurements as well. We also explore several other ways to estimate the error and find between $1.7$--$2.2\%$ precision on the BAO scale from the EZmock covariance matrix and between $1.1$--$1.5\%$ precision from the analytical covariance matrix. This work represents the first detection of the BAO feature in the DESI 3PCF, establishing its ability to probe the expansion history of the Universe with future DESI 3PCF measurements.