🔭 astrophysics

First Detection of the Baryon Acoustic Oscillation (BAO) Feature in the 3-Point Correlation Function of DESI DR1 Luminous Red Galaxies

This paper reports the first detection of the Baryon Acoustic Oscillation feature in the 3-Point Correlation Function of over 2.1 million DESI DR1 Luminous Red Galaxies, achieving a significance of approximately 8.1–8.5 $\sigma$ and measuring the distance scale $D_{\rm V}(z=0.68)/r_{\rm d}$ with 1.12–1.72% precision, thereby validating the 3PCF as a robust tool for probing cosmic expansion history.

Original authors: 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. An

Published 2026-02-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: 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

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, expanding balloon. Now, imagine that inside this balloon, there are billions of tiny dots representing galaxies. For a long time, astronomers have been trying to measure how fast this balloon is expanding and how it has changed over time. To do this, they need a "standard ruler"—a known distance they can use to measure the gaps between the dots.

This paper is about finding a new, more powerful ruler hidden inside the arrangement of those dots.

Here is the story of how they found it, explained simply:

1. The Ancient Echo (What is BAO?)

Billions of years ago, before the universe cooled down enough for stars to form, it was a hot, dense soup of particles and light. Think of this soup like a giant drum. When the universe was very young, sound waves (acoustic waves) rippled through this drum, pushing matter outward and then pulling it back, creating a rhythmic pattern.

When the universe cooled, these sound waves stopped. But they left a "fossil" imprint: a specific, preferred distance between galaxies. It's like if you dropped a pebble in a pond and the ripples froze in place. That frozen distance is about 150 million light-years. Astronomers call this the Baryon Acoustic Oscillation (BAO).

2. The Old Ruler vs. The New Ruler

For decades, astronomers have used a tool called the 2-Point Correlation Function (2PCF) to find this ruler.

The Analogy: Imagine you are looking at a crowd of people. The 2PCF asks: "If I pick one person, how likely am I to find another person standing exactly 150 million light-years away from them?"
This works well. It's like measuring the distance between two points on a map.

But this paper introduces a new tool: the 3-Point Correlation Function (3PCF).

The Analogy: Instead of just looking at pairs, the 3PCF asks: "If I pick three people, what kind of triangle do they form? Do they tend to form triangles with a specific side length?"
It's like looking at the shape of the crowd, not just the distance between neighbors. This gives a much richer picture of the universe's structure.

3. The Challenge: Finding a Needle in a Haystack

The problem with the 3-Point method is that it is incredibly complex.

The Math: Calculating the 2-Point method is like counting pairs of shoes in a room. Calculating the 3-Point method is like trying to count every possible triangle formed by people in that room. The number of combinations is massive.
The Noise: The universe is messy. Galaxies move, gravity pulls them, and our telescopes have limitations. Finding the faint "echo" of the ancient sound waves inside this noise is like trying to hear a whisper in a hurricane.

4. The Breakthrough: The "Super-Scanner"

The team, led by researchers at the University of Florida and others, used data from the DESI (Dark Energy Spectroscopic Instrument). This instrument is a high-tech camera on a telescope that can take pictures of the spectra (colors) of 2.1 million galaxies at once.

To solve the math problem, they used a clever algorithm called FFTLog.

The Analogy: Imagine you have a giant, messy pile of puzzle pieces. Instead of trying to fit them together one by one (which would take forever), you use a magic scanner that instantly organizes them into patterns. The FFTLog algorithm allowed them to process the 3-point data quickly and accurately, turning a super-complex math problem into something manageable.

5. The Result: A Clear Signal

They looked at the data and asked: "Is there a signal here, or is it just random noise?"

The Detection: They found the signal! They detected the BAO "fossil" in the 3-point data with a confidence level of 8.5 sigma.
- What does that mean? In science, a "sigma" is a measure of certainty. A 5-sigma result is usually enough to claim a discovery (like finding the Higgs Boson). An 8.5-sigma result is like being absolutely certain it's not a fluke. It's the difference between guessing it might rain and having a satellite image showing a storm cloud directly overhead.
The Precision: They measured the size of the universe's expansion with about 1.1% to 1.7% precision. This is incredibly sharp. It's like measuring the distance from New York to London and being off by only a few inches.

6. Why Does This Matter?

This is the first time anyone has successfully seen this "3-point ruler" in such a large dataset.

Cross-Checking: It's like having two different maps of the same territory. If the old map (2-point) and the new map (3-point) agree, we know the territory is real and our maps are accurate.
Future Proofing: This proves that the 3-point method works. In the future, as we get more data, this method will help us understand Dark Energy (the mysterious force pushing the universe apart) much better. It might even help us figure out if the rules of gravity change over time.

Summary

Think of this paper as the moment astronomers realized they didn't just have to look at pairs of stars to understand the universe's expansion. They realized they could look at triangles of stars. By using a super-fast computer algorithm and a massive telescope, they found a clear, ancient echo in the shape of the universe, confirming that their new "3-point ruler" is accurate and ready to help us solve the biggest mysteries of the cosmos.

1. Problem Statement

While Baryon Acoustic Oscillations (BAO) have been robustly measured using the 2-Point Correlation Function (2PCF) and Power Spectrum in large-scale structure surveys, the 3-Point Correlation Function (3PCF) has historically been underutilized for BAO detection.

Challenges: The 3PCF is computationally expensive and theoretically complex to model due to non-linear gravitational evolution, redshift-space distortions, and galaxy biasing.
Goal: This work aims to demonstrate the first detection of the BAO feature in the 3PCF using the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) Luminous Red Galaxy (LRG) sample. The objective is to validate the 3PCF as an independent, methodologically distinct probe of cosmic expansion history and to test the robustness of systematic corrections (e.g., fiber collisions, edge effects) applied to 2PCF analyses.

2. Methodology

A. Data Sample

Source: DESI DR1 LRG sample.
Size: Over 2.1 million galaxies.
Redshift Range: $0.4 < z < 1.1$ .
Configuration: The authors combined three non-overlapping redshift bins (LRG1, LRG2, LRG3) into a single high-volume sample denoted as $z_{full}$ (effective redshift $z_{eff} = 0.68$ ).
Exclusion: Small-scale configurations ( $|r_1 - r_2| < 20 h^{-1}$ Mpc) were excluded to remain within the validity of the tree-level perturbation theory model.

B. Theoretical Modeling

The analysis relies on a tree-level redshift-space bispectrum template transformed into configuration space.

Bispectrum Expansion: The redshift-space bispectrum is expanded in terms of spherical harmonics to isolate the monopole ( $B^{(0,0)}$ ), which depends on three variables ( $k_1, k_2, \mu_{12}$ ).
Bias Parameters: The model includes three free galaxy bias parameters: linear bias ( $b_1$ ), quadratic bias ( $b_2$ ), and tidal bias ( $b_{G2}$ ).
FFTLog Algorithm: To compute the 3PCF from the bispectrum, the authors utilize the Fast Fourier Transform on Logarithmic scales (FFTLog) algorithm.
- This converts the bispectrum inner multipoles into 3PCF inner multipoles ( $\zeta_\ell$ ) via double Bessel transforms.
- Advantage: This numerical approach is faster and more flexible than purely analytical formulations, allowing for future inclusion of higher-order loop corrections.
BAO Isolation: To measure detection significance, two templates are generated:
1. Wiggle: Contains the full power spectrum with BAO oscillations.
2. No-Wiggle: The BAO feature is removed using a Fast Sine Transform (FST) method on the power spectrum, interpolating out the localized bump in harmonic space.

C. Statistical Analysis & Covariance

Covariance Matrices: Two approaches were used to estimate errors:
1. EZmocks: 1,000 realizations using the Effective Zeldovich approximation. Corrections (Hartlap and Percival factors) were applied to account for finite mock counts and noise.
2. Analytical Gaussian Covariance: Calibrated against EZmocks using the mean number density ( $\bar{n}$ ) and effective volume ( $V_{eff}$ ).
Detection Significance: Calculated via $\Delta\chi^2 = \chi^2_{nw} - \chi^2_{w}$ , where the difference in chi-squared between the no-wiggle and wiggle models determines the significance in standard deviations ( $\sigma$ ).
Systematics Validation: The model was tested against Abacus altMTL mocks (high-resolution N-body simulations with realistic fiber assignment) to quantify modeling offsets.

3. Key Contributions

First Detection: This is the first detection of the BAO signal in the 3PCF of a galaxy survey, achieving a significance of 8.1 $\sigma$ (EZmock) and 8.5 $\sigma$ (Analytical).
Methodological Innovation: The successful application of the FFTLog algorithm to transform the redshift-space bispectrum into the 3PCF, offering a computationally efficient alternative to previous analytical methods.
Systematic Validation: By comparing data against the Abacus altMTL mocks, the authors quantified a modeling systematic offset of only 0.6% in the recovered BAO scale, validating the robustness of the 3PCF modeling pipeline.
Independent Cross-Check: Provided an independent measurement of the BAO scale that is methodologically distinct from the 2PCF, confirming that systematics (like fiber collisions) do not bias the BAO feature in higher-order statistics.

4. Key Results

Detection Significance:
- 8.1 $\sigma$ using the EZmock covariance matrix.
- 8.5 $\sigma$ using the analytical covariance matrix.
- This is a significant improvement over the previous SDSS BOSS CMASS detection (4.5 $\sigma$ ).
Distance Measurements:
- The angle-averaged distance $D_V(z)/r_d$ was measured at $z=0.68$ .
- EZmock: $15.88 \pm 0.27$ (1.72% precision).
- Analytical: $15.72 \pm 0.18$ (1.12% precision).
Consistency:
- The results are in excellent agreement with DESI DR1 2PCF measurements and the Planck 2018 cosmology.
- The measured BAO scale is marginally smaller than the Planck prediction, consistent with other DESI post-reconstruction 2PCF results.
Error Budget:
- Statistical errors dominate, but the systematic error from modeling is constrained to 0.6%.
- Alternative error estimation methods (bootstrapping, mock scatter, marginalization) yielded consistent results (1.1% – 2.2% precision), confirming the robustness of the standard method.

5. Significance

Validation of 3PCF: This work proves that the 3PCF is a viable and powerful tool for BAO cosmology, capable of reaching sub-2% precision on distance scales.
Future Prospects: The FFTLog-based framework allows for the inclusion of higher-order corrections (loop terms), which will be crucial for analyzing future, larger DESI data releases (DR2 and beyond) and for performing joint 2PCF+3PCF full-shape analyses.
Cosmological Constraints: The 3PCF offers a methodologically independent check on the expansion history of the Universe. Its agreement with 2PCF results strengthens confidence in the standard $\Lambda$ CDM model and the handling of complex survey systematics.
Benchmark: The detection of BAO in the 3PCF of 2.1 million galaxies sets a new benchmark for large-scale structure analysis, demonstrating that higher-order statistics can now be routinely applied to massive spectroscopic surveys.