Construction of the Damped Lyα Absorber Catalog for DESI DR2 Lyα BAO

We present the Damped Ly$\alpha$ Toolkit for automated detection and characterization of Damped Ly$\alpha$ absorbers (DLA) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance ($\Delta \chi_{r}^2>0.03$), the technique achieves an estimated 72% purity and 71% completeness when evaluated on simulated spectra with S/N$>2$ that are free of broad absorption lines (BAL). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21,719 were found in S/N$>2$ spectra with predicted $\log_{10} (N_\texttt{HI}) > 20.3$ and detection significance $\Delta \chi_{r}^2 >0.03$. We compare the Damped Ly$\alpha$ Toolkit to two alternative DLA finders based on a convolutional neural network (CNN) and Gaussian process (GP) models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly$\alpha$ forest baryon acoustic oscillation measurement. The combined catalog contains 41,152 candidate DLAs with $\log_{10} (N_\texttt{HI}) > 20.3$ from quasar spectra with S/N$>2$. We estimate this sample to be approximately 76% pure and 71% complete when BAL quasars are excluded.