BAO: A Cosmic Standard Ruler

• BAO are quasi-periodic features in the matter distribution induced by early-Universe acoustic waves, offering a fixed standard ruler for cosmological measurements.
• Analyses utilize the galaxy two-point correlation function and power spectrum to extract sub-percent precision distance scales and expansion rates.
• Advanced observational techniques and statistical methods mitigate nonlinear damping and environmental biases, enhancing BAO's role in constraining dark energy and spatial curvature.

Baryon Acoustic Oscillations (BAO) are large-scale, quasi-periodic modulations in the spatial distribution of matter in the Universe, originating from acoustic waves that propagated in the tightly coupled photon–baryon plasma prior to cosmic recombination. These features provide an exceptionally robust standard ruler for late-time cosmology, enabling sub-percent measurements of the cosmic expansion history and the geometry of spacetime. The BAO signature manifests as a local maximum (“bump”) in the two-point correlation function of galaxies at a comoving separation corresponding to the sound horizon at the baryon drag epoch, and as a series of oscillatory “wiggles” in the matter power spectrum at characteristic wavenumbers. BAO measurements form one of the empirical pillars of the standard cosmological model and play a central role in constraining dark energy, spatial curvature, and early-Universe physics including the baryon fraction and the effective number of relativistic species.

1. Physical Origin and Linear Theory of BAO

BAO arise from sound waves in the primordial photon–baryon fluid. Prior to recombination ($z > 1100$), photons and baryons were tightly coupled via Thomson scattering. Gravitational wells produced by density perturbations generated pressure-supported acoustic waves, with characteristic sound speed $c_s = c / \sqrt{3(1+R)}$, $R=3\rho_b/4\rho_\gamma$ (Ferreira et al., 2024, 0910.5224). As the Universe expanded, these waves propagated to a maximum comoving distance by decoupling—the sound horizon $r_s$, given by

$$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$

where $z_d$ is the baryon drag epoch ($z_d \approx 1020$). After decoupling, baryons were left distributed as concentric shells of overdensity at $r \sim r_s$ from each perturbation center; dark matter, being decoupled, remained centrally concentrated (Arnalte-Mur et al., 2011, 0910.5224).

The BAO imprint is preserved as a broad, low-contrast “bump” at $r \approx 100-150\, h^{-1}\,\mathrm{Mpc}$ in the two-point correlation function $\xi(r)$ and as oscillatory modulations in the power spectrum $c_s = c / \sqrt{3(1+R)}$0 with period $c_s = c / \sqrt{3(1+R)}$1 (Ferreira et al., 2024, 0910.5224).

2. Statistical Detection and Quantification

2.1. Two-Point Correlation Function and Power Spectrum Analyses

The primary tools for BAO detection are the galaxy two-point correlation function

$c_s = c / \sqrt{3(1+R)}$2

and the power spectrum

$c_s = c / \sqrt{3(1+R)}$3

The BAO corresponds to a statistically significant excess of pairs at $c_s = c / \sqrt{3(1+R)}$4 or oscillatory features in $c_s = c / \sqrt{3(1+R)}$5. In practice, $c_s = c / \sqrt{3(1+R)}$6 is modeled as

$c_s = c / \sqrt{3(1+R)}$7

where $c_s = c / \sqrt{3(1+R)}$8 is the smooth “no-wiggle” component and $c_s = c / \sqrt{3(1+R)}$9 encodes the oscillatory structure (Ferreira et al., 2024, 0910.5224).

Template fitting methodologies marginalize over broadband shape (polynomial terms or splines) and introduce a “dilation parameter” $R=3\rho_b/4\rho_\gamma$0 that rescales the fiducial acoustic scale, yielding

$R=3\rho_b/4\rho_\gamma$1

where $R=3\rho_b/4\rho_\gamma$2 models nonlinear damping (Ferreira et al., 2024, Ruggeri et al., 2019, 2209.06011). This parameter is directly related to the ratio of the measured distance scale to the fiducial value via volume-averaged distance measures

$R=3\rho_b/4\rho_\gamma$3

Isotropic or multipole expansions (monopole, quadrupole) are fitted to extract $R=3\rho_b/4\rho_\gamma$4 and thereby constrain cosmological parameters (Ferreira et al., 2024, 2209.06011).

2.2. Configuration vs. Fourier vs. Joint-Space Approaches

BAO analyses may be performed in configuration space, Fourier space, or jointly. Joint fitting frameworks for power spectrum and correlation function multipoles reduce parameter bias and provide more realistic (often slightly larger) uncertainty estimates by leveraging the full covariance between spaces without defaulting to Gaussian likelihood approximations (2209.06011).

2.3. Spherical (3D) Analysis

For full-sky or deep surveys, 3D Spherical Fourier-Bessel (SFB) decompositions capture both radial and tangential BAO information, with the “radialization” effect ensuring that deep survey BAO measurements closely approach those from Cartesian $R=3\rho_b/4\rho_\gamma$5 (Pratten et al., 2013, Rassat et al., 2011).

3. BAO Measurement Modalities and Tracer Strategies

3.1. Galaxy Surveys: Spectroscopic vs Photometric

Spectroscopic surveys provide precise redshifts ($R=3\rho_b/4\rho_\gamma$6) enabling both radial and tangential BAO constraints, i.e., independent determination of $R=3\rho_b/4\rho_\gamma$7 and $R=3\rho_b/4\rho_\gamma$8. Photometric surveys, in contrast, suffer line-of-sight smearing and restrict BAO extraction primarily to angular scales, measuring $R=3\rho_b/4\rho_\gamma$9 alone with larger errors (Ferreira et al., 2024, Alcaniz et al., 2016).

3.2. Voids, Density-Dependent, and High/Low Density Tracers

BAO can be detected not only from galaxies but from advanced tracers:

• Delaunay Triangulation (DT) voids: Overlapping empty circumspheres from tetrahedra in the galaxy distribution reveal a strong BAO signal, with the two-point function peaking at $r_s$0 (average void radius). Radius-thresholding separates low-density and high-density tracers; high-density tracers yield up to a 12$r_s$1 enhancement in BAO amplitude over standard galaxy correlations (Liu et al., 2017).
• Density minima (“voids”): Overlapping voids (not disjoint) constructed from all empty circumspheres capture the full BAO signal, whereas disjoint, maximal spheres do not. Detection of BAO from minima provides a complementary, less nonlinearly-degraded standard ruler (Kitaura et al., 2015).
• Density-split statistics: By quantile-splitting on local density (e.g., DS$r_s$2, $r_s$3), one observes systematic environmental shifts in the BAO scale: nonlinear flows contract (expand) the BAO in overdense (underdense) regions by up to $r_s$4 from $r_s$5 to $r_s$6 (Xu et al., 2024). These effects can be partially mitigated by local-density-based weighting in BAO reconstruction.

3.3. Real-Space Structure Identification and Wavelet Methods

Wavelet-based techniques (e.g., “BAOlets”) and convolutional approaches (e.g., CenterFinder algorithm) identify actual spatial shell structures centered on high-bias objects (e.g., LRGs, massive halos) and quantify shell properties (radius, contrast). In SDSS samples, real-space BAO shells are identified with radii $r_s$7 and overdensity $r_s$8 (Arnalte-Mur et al., 2011, Brown et al., 2020).

3.4. Alternative Tracers: Intensity Mapping, Gravitational Wave Sirens

HI intensity mapping (e.g., SKA1) recovers the radial BAO at $r_s$9 with $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$0 precision, while galaxy redshift surveys (SKA2) achieve sub-percent constraints in $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$1 (Bull et al., 2015). Third-generation gravitational wave observatories localize $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$2 mergers/year within $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$3 at $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$4, sufficient to reconstruct the BAO feature to $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$5 precision in multiple redshift shells, independent of galaxy or quasar surveys (Kumar et al., 2021).

4. Cosmological Distance Inference and Parameter Constraints

BAO serve as a standard ruler due to the calibratable sound horizon $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$6, with CMB measurements providing $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$7 anchor precision (Sutherland, 2012, 0910.5224). Observational BAO peak positions are mapped to $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$8, $$r_s(z_d) = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} dz,$$9, and $z_d$0, with dilation parameters $z_d$1 and $z_d$2 quantifying deviations from a fiducial cosmology.

The principal cosmological constraints derived include:

• Expansion history: Combined measurements of $z_d$3 and $z_d$4 from BAO break degeneracies between curvature and dark energy equation of state parameters (0910.5224).
• Absolute ruler tests: New $z_d$5–$z_d$6 relations enable CMB-independent, low-$z_d$7 measurements of $z_d$8 using supernovae or gravitational-wave standard sirens, serving as “zero-parameter” consistency tests of early-Universe physics (Sutherland, 2012).
• Probing radiation density $z_d$9: The amplitude (“area”) of the BAO bump, quantified by the $z_d \approx 1020$0 statistic, is directly sensitive to the cosmic baryon fraction $z_d \approx 1020$1 and $z_d \approx 1020$2. Raising $z_d \approx 1020$3 (fixing $z_d \approx 1020$4 and CMB observables) reduces $z_d \approx 1020$5 and suppresses $z_d \approx 1020$6, providing a probe of dark radiation independent of the CMB damping tail (Sutherland et al., 2014).

5. Statistical Methodologies, Error Analysis, and Systematics

BAO analyses deploy rigorous statistical frameworks:

• Template-based likelihood: BAO-scale parameters are estimated by likelihood maximization over polynomial-marginalized models with a joint covariance matrix derived from large ensembles of N-body or approximate-mock catalogues (e.g., PATCHY, QPM, MultiDark, EZmocks).
• Error estimation: Three primary error metrics are employed: (i) posterior width from data, (ii) ensemble mean from mocks, and (iii) Fisher matrix-based forecast. For high-volume surveys, Fisher-matrix and empirical errors agree; in low-volume regimes, substantial sample variance affects error quantification and must be supplemented by mock-based estimates (Ruggeri et al., 2019).
• Covariance modeling: The full data covariance (including cross-covariances between configuration and Fourier space) is crucial for unbiased parameter inference (2209.06011).
• Wavelet and model-independent detection: Wavelet-based and bump-robusticity methods reduce reliance on model choices, facilitating detection even in the presence of systematics (Labatie et al., 2011, Arnalte-Mur et al., 2011, Alcaniz et al., 2016).
• Systematics and environmental dependencies: BAO peak shifts due to nonlinear flows, local curvature, or environmental effects can induce biases at the percent level; advanced environmental weighting, reconstruction algorithms, and relativistic corrections are active mitigation strategies (Xu et al., 2024, Roukema, 2015).

6. Future Directions and Open Issues

Advances in survey methodology and data volume (e.g., DESI, Euclid, LSST, SKA, future GW detectors) will enable percent-level BAO measurements over $z_d \approx 1020$7, transforming constraints on dark energy, curvature, and early-Universe parameters (0910.5224, Bull et al., 2015, Ferreira et al., 2024).

Open problems include:

• Redshift evolution of $z_d \approx 1020$8: Current measurements of the dilation parameter $z_d \approx 1020$9 scatter without a clear trend in $r \sim r_s$0, an unresolved issue seen in both tomographic and 3D analyses (Ferreira et al., 2024).
• BAO flexibility: Observational evidence supports departures from strict comoving rigidity (“flexible ruler”) due to inhomogeneous curvature and kinematical backreaction, necessitating more sophisticated modeling of the BAO standard ruler to avoid cosmological parameter bias (Roukema, 2015).
• Amplitude-based cosmology: Achieving the statistical and systematic control to exploit the amplitude of the BAO bump and break degeneracies in $r \sim r_s$1, $r \sim r_s$2 remains an ongoing challenge, requiring simulation advances and next-generation data sets (Sutherland et al., 2014).

7. Summary Table: Key BAO Observational and Modeling Aspects

Aspect	Key Observable	Typical Approach or Metric
Standard ruler	$r \sim r_s$3	$r \sim r_s$4, $r \sim r_s$5
Statistical model	Full-shape fitting	Polynomial-marginalized, template-based $r \sim r_s$6
Error quantification	$r \sim r_s$7	Posterior width, mock ensemble, Fisher forecast
Nonlinear effects	Peak shift/damping	$r \sim r_s$8, density/environmental splits
Environmental tracing	DT voids, DS stats	Enhancement/shift in BAO amplitude/position
Advanced detection	Joint-space, wavelet	Full covariance, model-independent features
CMB-independence	$r \sim r_s$9–$r \approx 100-150\, h^{-1}\,\mathrm{Mpc}$0 link	Low-$r \approx 100-150\, h^{-1}\,\mathrm{Mpc}$1 SNe, GW sirens, zero-parameter tests

BAO provide one of the most precise and physically grounded distance calibration mechanisms in cosmology, underpinned by linear early-Universe physics and robust to most systematic uncertainties. Ongoing survey, methodological, and theoretical advancements are systematically pushing the accuracy, cosmic volume, and physical reach of the BAO standard ruler.