Extensive analysis of reconstruction algorithms for DESI 2024 baryon acoustic oscillations

X. Chen,Z. Ding,E. Paillas,S. Nadathur,H. Seo,S. Chen,N. Padmanabhan,M. White,A. de Mattia,P. McDonald,A. J. Ross,A. Variu,A. Carnero Rosell,B. Hadzhiyska,M. M. S Hanif,D. Forero-Sánchez,S. Ahlen,O. Alves,U. Andrade,S. BenZvi,D. Bianchi,D. Brooks,E. Chaussidon,T. Claybaugh,A. de la Macorra,Biprateep Dey,K. Fanning,S. Ferraro,A. Font-Ribera,J. E. Forero-Romero,C. Garcia-Quintero,E. Gaztañaga,S. Gontcho A Gontcho,G. Gutierrez,C. Hahn,K. Honscheid,S. Juneau,R. Kehoe,D. Kirkby,T. Kisner,A. Kremin,M. E. Levi,A. Meisner,J. Mena-Fernández,R. Miquel,J. Moustakas,A. Muñoz-Gutiérrez,F. Nikakhtar,N. Palanque-Delabrouille,W. J. Percival,F. Prada,I. Pérez-Ràfols,M. Rashkovetskyi,G. Rossi,R. Ruggeri,E. Sanchez,C. Saulder,D. Schlegel,M. Schubnell,A. Smith,D. Sprayberry,G. Tarlé,D. Valcin,M. Vargas-Magaña,B. A. Weaver,S. Yuan,R. Zhou
2024-11-29
Abstract:Reconstruction of the baryon acoustic oscillation (BAO) signal has been a standard procedure in BAO analyses over the past decade and has helped to improve the BAO parameter precision by a factor of ~2 on average. The Dark Energy Spectroscopic Instrument (DESI) BAO analysis for the first year (DR1) data uses the ``standard'' reconstruction framework, in which the displacement field is estimated from the observed density field by solving the linearized continuity equation in redshift space, and galaxy and random positions are shifted in order to partially remove nonlinearities. There are several approaches to solving for the displacement field in real survey data, including the multigrid (MG), iterative Fast Fourier Transform (iFFT), and iterative Fast Fourier Transform particle (iFFTP) algorithms. In this work, we analyze these algorithms and compare them with various metrics including two-point statistics and the displacement itself using realistic DESI mocks. We focus on three representative DESI samples, the emission line galaxies (ELG), quasars (QSO), and the bright galaxy sample (BGS), which cover the extreme redshifts and number densities, and potential wide-angle effects. We conclude that the MG and iFFT algorithms agree within 0.4% in post-reconstruction power spectrum on BAO scales with the RecSym convention, which does not remove large-scale redshift space distortions (RSDs), in all three tracers. The RecSym convention appears to be less sensitive to displacement errors than the RecIso convention, which attempts to remove large-scale RSDs. However, iFFTP deviates from the first two; thus, we recommend against using iFFTP without further development. In addition, we provide the optimal settings for reconstruction for five years of DESI observation. The analyses presented in this work pave the way for DESI DR1 analysis as well as future BAO analyses.
Cosmology and Nongalactic Astrophysics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to effectively reconstruct the Baryon Acoustic Oscillation (BAO) signal in the case of different observational line direction changes. Specifically, the paper focuses on the fact that in actual cosmological surveys, such as the Dark Energy Spectroscopic Instrument (DESI) project, due to the change in the observational line direction and the complex survey geometry, the standard BAO reconstruction algorithm is difficult to be directly applied. Therefore, the focus of the research is to analyze and compare three different reconstruction algorithms - Multigrid (MG), Iterative Fast Fourier Transform (iFFT), and Iterative Fast Fourier Transform Particle (iFFTP) - to evaluate their performance in dealing with these complex conditions. ### Core problems of the paper: 1. **Influence of non - parallel observational lines**: In actual cosmological surveys, the observational lines (i.e., the line - of - sight directions) are not in a single direction, which makes it more difficult to solve the linearized redshift - space continuity equation. 2. **Algorithm selection and optimization**: The paper conducts a detailed analysis and comparison of the above three algorithms by using realistic DESI simulation data to determine which algorithm performs best when dealing with different types of samples (such as Emission Line Galaxies (ELG), Quasars (QSO), and Bright Galaxy Samples (BGS)). 3. **Improvement of reconstruction accuracy**: By comparing different algorithms, the paper aims to find a method that can effectively improve the accuracy of BAO parameter measurement, especially in the first five - year data (DR1) of the DESI project. ### Main contributions: - **Comprehensive analysis of algorithm performance**: The paper provides a detailed comparison of the MG, iFFT, and iFFTP algorithms under different indicators, including two - point statistics, the displacement itself, etc. - **Recommendation of optimal settings**: Based on the analysis results, the paper provides optimal setting suggestions for the reconstruction of the DESI project, especially the grid resolution and the number of random points. - **Guidance for future work**: The results of the paper are not only directly helpful for the analysis of DESI DR1, but also provide an important reference for future BAO analysis. ### Conclusion: The main conclusion drawn by the paper is that the performance of the MG and iFFT algorithms on the post - reconstruction power spectrum is consistent, with an error within 0.4%, and is applicable to all three tracers (ELG, QSO, and BGS). However, the iFFTP algorithm has a large deviation from the other two algorithms, so it is not recommended to be used without further improvement. In addition, the paper also discusses the sensitivity of different reconstruction methods to displacement errors and provides ideas for the application of future non - standard reconstruction algorithms.