Cosmological forecast from the full-sky angular power spectrum and bispectrum of 21cm intensity mapping

This paper presents the first full-sky relativistic bispectrum forecast for BINGO and SKA1-MID surveys, demonstrating that including second-order velocity contributions and bispectrum data significantly improves constraints on dark energy parameters and the Hubble constant by breaking degeneracies that limit power spectrum-only analyses.

The Gist

The Big Picture: Listening to the Universe's "Hum"

Imagine the universe as a giant, dark ocean. For a long time, we've been trying to understand why this ocean is expanding faster and faster (a mystery called Dark Energy).

Usually, astronomers look at the "waves" on the surface of this ocean by counting galaxies (like counting fish). But there's another way: listening to the "hum" of the ocean itself. This hum is made of neutral hydrogen, the most common gas in the universe. When this gas vibrates, it emits a specific radio signal (the 21cm line).

This paper is about two giant radio telescopes, BINGO (in Brazil) and SKA (in South Africa), which are designed to "listen" to this hum across the entire sky. The authors are asking: If we listen to this hum, how much better can we understand the secrets of the universe?

The Two Ways of Listening: The "Solo" vs. The "Chorus"

To understand the paper's main discovery, you need to know about two ways of analyzing this radio signal:

1. The Power Spectrum (The "Solo"):
   Think of this as listening to the volume of the hum at different pitches. It tells us how much hydrogen is clumped together in big blobs versus small blobs. It's like looking at a single note on a piano and saying, "That note is loud." This has been the standard way to study the universe for decades.

2. The Bispectrum (The "Chorus"):
   This is the paper's star. The bispectrum listens to how three different notes interact at the same time. It asks: "When this note plays, does it change how that other note sounds?"

   • The Analogy: Imagine a choir. The "Power Spectrum" tells you how loud the choir is. The "Bispectrum" tells you how the singers are harmonizing or clashing. It reveals the relationships between the notes, not just the volume.

The Paper's Big Finding:
For a long time, scientists thought the "Chorus" (Bispectrum) was just a fancy extra that gave similar results to the "Solo" (Power Spectrum). This paper proves that the Chorus is actually a superpower, especially when trying to figure out the nature of Dark Energy.

The Secret Ingredient: The "Second-Order Velocity"

The authors discovered a hidden ingredient in the signal that most people were ignoring.

• The Problem: When gas moves, it creates a Doppler shift (like the sound of a siren passing you). Usually, scientists use a shortcut (called the "Limber approximation") to ignore the complex, second-order effects of this movement because they are hard to calculate.
• The Discovery: The authors did the hard math and found that this "ignored" movement actually accounts for 24% of the total signal at low redshifts (relatively close to us).
• The Metaphor: Imagine trying to predict the weather by only looking at the wind speed. You ignore the humidity and temperature changes. The authors found that ignoring the "humidity" (the second-order velocity) means you are missing nearly a quarter of the story. If you want an accurate forecast, you must include it.

The Results: Why This Matters

The team used a mathematical tool called the Fisher Matrix (think of it as a "precision calculator") to predict how well these telescopes will work. Here is what they found:

1. Breaking the "Tangled Knot":
   In cosmology, different parameters (like the expansion rate of the universe and the strength of Dark Energy) often get "tangled" together. It's like trying to untie a knot where pulling one string just tightens another.

   • The "Solo" (Power Spectrum) struggles to untie this knot.
   • The "Chorus" (Bispectrum) cuts right through it. By adding the Bispectrum data, the team found they could improve the measurement of Dark Energy's behavior by over 70%.

2. The "Dynamic" Dark Energy:
   Scientists suspect Dark Energy might not be a constant force (like a cosmological constant) but something that changes over time (Dynamical Dark Energy).

   • The Power Spectrum is okay at guessing this.
   • The Bispectrum is excellent at it. It acts like a high-resolution microscope, allowing us to see if Dark Energy is changing its mind as the universe ages.

3. BINGO vs. SKA:

   • BINGO is a single, large dish. It's like a good pair of ears. It will give us great data.
   • SKA is a massive array of many dishes. It's like a super-hearing system with thousands of ears. It will be roughly twice as powerful as BINGO, giving us even sharper constraints on the universe's history.

The Bottom Line

This paper is a roadmap for the future of radio astronomy. It tells us:

• Don't just listen to the volume of the universe; listen to the harmonies (the Bispectrum).
• Don't ignore the complex movements of gas (the second-order velocity); they make up a huge chunk of the signal.
• By combining these new listening techniques with existing data from the Planck satellite, we can finally untangle the mystery of Dark Energy and understand why the universe is speeding up.

In short: We are moving from listening to a single note to conducting a full symphony, and that symphony holds the secrets to the universe's fate.

Technical Summary

1. Problem Statement

Modern cosmology faces the challenge of explaining the accelerated expansion of the universe and the nature of dark energy. While the Planck satellite has provided precise constraints on early-universe parameters via the Cosmic Microwave Background (CMB), it is insufficient for tightly constraining dynamical dark energy models (e.g., time-varying equation of state).

Future radio telescopes like BINGO (Brazilian) and SKA1-MID (South African) will map neutral hydrogen (HI) via 21cm intensity mapping. However, standard cosmological forecasts often rely on:

1. The Flat-Sky Approximation: Which fails for wide-angle correlations and large survey volumes.
2. The Limber Approximation: Which neglects second-order relativistic effects, particularly the second-order velocity contribution, on large angular scales ($\ell \lesssim 100$) and narrow redshift bins.
3. Two-Point Statistics Only: Relying solely on the angular power spectrum ($C_\ell$) ignores the rich information contained in the three-point correlation function (bispectrum, $B_{\ell_1 \ell_2 \ell_3}$), which is crucial for breaking parameter degeneracies in dynamical dark energy models.

The paper aims to perform the first full-sky relativistic forecast for these surveys using both the angular power spectrum and the bispectrum to forecast constraints on cosmological parameters.

2. Methodology

A. Theoretical Framework

• Full-Sky Formalism: The authors avoid the flat-sky approximation, utilizing spherical harmonics to analyze the full sky. This preserves large-scale modes and avoids edge effects inherent in Fourier transforms of finite domains.
• Relativistic Bispectrum: The analysis includes the full relativistic bispectrum in redshift space. The HI temperature contrast ($\Delta$) is expanded to second order:
  • $\Delta^{(1)}$: Linear order (density and redshift-space distortion).
  • $\Delta^{(2)}$: Second-order terms including non-linear density, velocity, and tidal fields.
• Key Contribution of Velocity Terms: The paper highlights that the second-order velocity contribution ($v^{(2)\prime}$) is critical. Under the Limber approximation, this term is often neglected, but the authors show it accounts for a significant portion of the signal at low redshifts.
• Approximation Strategy: Calculating the full second-order velocity term involves computationally expensive quadruple integrals. The authors developed an approximation where the fractional contribution of this term ($r$) is computed numerically for a subset of configurations and found to be stable ($r \approx 0.24$). This allows the total bispectrum to be reconstructed without explicit calculation of the expensive term in the Fisher matrix, reducing computational cost while maintaining accuracy (error < 2.2%).

B. Survey Specifications

The forecasts are based on two specific configurations:

1. BINGO: A single-dish telescope covering $0.13 < z < 0.45$ with $2900 \text{ deg}^2$ sky coverage.
2. SKA1-MID Band 2: A dish array covering $0.01 < z < 0.5$ with $5000 \text{ deg}^2$ sky coverage and higher sensitivity.

C. Fisher Matrix Analysis

• Models: Constraints are forecasted for three models: $\Lambda$CDM, $w$CDM (constant dark energy equation of state), and the Chevallier-Polarski-Linder (CPL) model (dynamical dark energy with parameters $w_0$ and $w_a$).
• Data Combinations: The analysis compares:
  • Planck CMB data alone.
  • $C_\ell$ (Power Spectrum) + Planck.
  • $B_\ell$ (Bispectrum) + Planck.
  • $C_\ell + B_\ell$ + Planck.
• Constraints: The analysis is restricted to linear scales ($\ell < \ell_{nl}^{max}$) and redshift bin auto-correlations to ensure theoretical control and manage computational complexity.

3. Key Contributions

1. First Full-Sky Relativistic Bispectrum Forecast: This is the first study to apply the full-sky relativistic bispectrum (without Limber approximation) to forecast constraints for BINGO and SKA1-MID.
2. Quantification of Second-Order Velocity: The authors demonstrate that the second-order velocity term contributes approximately 24% of the total bispectrum signal at low redshifts ($z \lesssim 0.5$). Neglecting this term leads to inaccurate modeling.
3. Efficient Numerical Approximation: They introduced a method to approximate the computationally prohibitive second-order velocity term using a constant fractional factor ($r \approx 0.24$), enabling efficient Fisher matrix calculations without sacrificing accuracy.
4. Degeneracy Breaking: The study explicitly demonstrates how the bispectrum breaks degeneracies between dark energy parameters ($w_0, w_a$) and the Hubble parameter ($h$) that persist in power spectrum analyses.

4. Key Results

A. Impact of the Bispectrum

• $\Lambda$CDM and $w$CDM Models: The bispectrum provides constraints comparable to the power spectrum. However, the power spectrum generally yields slightly tighter constraints on the linear bias ($b_1$) and the Hubble parameter ($h$) in these simpler models.
• CPL (Dynamical Dark Energy) Model: The bispectrum becomes a powerful probe.
  • In the CPL model, the power spectrum suffers from strong degeneracies between $w_0$, $w_a$, and $h$.
  • The bispectrum probes non-linear mode coupling and growth history with a parameter dependence rotated relative to the power spectrum.
  • Joint Analysis ($C_\ell + B_\ell + \text{Planck}$):
    • Improves constraints on $w_0$ and $w_a$ by over 70% compared to using the power spectrum alone.
    • Improves constraints on the Hubble parameter ($h$) by approximately 60%.

B. Telescope Performance Comparison

• SKA vs. BINGO: SKA1-MID Band 2 offers significantly better constraining power due to its larger survey volume and lower system temperature.
  • For $\Lambda$CDM, SKA improves constraints on $\Omega_c h^2$ and $h$ by roughly 30-40%, doubling the improvement seen with BINGO.
  • SKA shows a cumulative improvement in constraints on primordial parameters ($n_s, \Omega_b h^2$) across all models, whereas BINGO's sensitivity to these is limited.

C. Specific Parameter Improvements (CPL Model, SKA)

When combining $C_\ell + B_\ell + \text{Planck}$:

• $w_0$: 87% improvement over Planck alone.
• $w_a$: 87% improvement over Planck alone.
• $h$: 92% improvement over Planck alone.
• Relative to $C_\ell + \text{Planck}$: The addition of the bispectrum yields a 77% improvement for $w_0$ and $w_a$, and a 63% improvement for $h$.

5. Significance and Conclusion

This paper establishes that relativistic effects and higher-order statistics are essential for next-generation 21cm intensity mapping surveys.

• Necessity of Full-Sky Treatment: For wide-angle, low-redshift surveys like BINGO and SKA, the flat-sky and Limber approximations are insufficient. The full-sky approach is required to capture the significant second-order velocity contributions.
• Probing Dark Energy: While the power spectrum is excellent for standard cosmological parameters, the bispectrum is the key to unlocking dynamical dark energy constraints. It breaks the degeneracies that limit the power spectrum's ability to distinguish between a cosmological constant and evolving dark energy.
• Future Outlook: The results underscore the importance of including the bispectrum in the analysis pipelines for upcoming surveys. Future work should extend this framework to include cross-correlations between redshift bins and more complex bias modeling to further tighten constraints.

In summary, the inclusion of the full-sky relativistic bispectrum transforms 21cm intensity mapping from a tool for measuring expansion history into a precision probe for the fundamental nature of dark energy.