Mapping the star formation peak with LIGO A# and Next-Generation detectors

This study demonstrates that networks of upgraded LIGO detectors (A#) and next-generation observatories (Cosmic Explorer and Einstein Telescope) can independently constrain the peak redshift of the star formation rate with precisions of $\pm 0.1$ and $\pm 0.02$, respectively, by analyzing the redshift evolution of binary black hole mergers.

The Gist

Imagine the universe as a giant, bustling city that has been growing and changing for billions of years. One of the most important things astronomers want to know is: When did the city have its biggest "construction boom"? In cosmic terms, this is the "peak" of star formation—the moment in history when the universe was creating the most new stars.

Currently, astronomers try to find this peak by looking at light (electromagnetic waves) from distant galaxies. But looking at light is like trying to count the people in a crowded room through a thick, foggy window. The light gets distorted by dust, and it's hard to tell exactly how many stars are actually being born versus how bright they just happen to be.

The New Tool: Listening to the "Crash"
This paper proposes a new way to solve the puzzle: Gravitational Waves.

Think of gravitational waves as the sound of two heavy objects (like black holes) crashing into each other. Unlike light, these "sounds" travel through the universe without getting blocked by dust or fog. By listening to these crashes, scientists can figure out exactly when and where they happened, giving them a direct count of star formation history without the "foggy window" problem.

The Experiment: Two Different Microphones
The researchers simulated a year's worth of these black hole crashes to see how well different "microphones" (detectors) could find the peak of the construction boom. They tested two setups:

1. The "Upgraded" Microphone (LIGO-A#): This is a major upgrade to the current detectors. It's like swapping a standard microphone for a high-end studio one.
2. The "Super" Microphone (Next-Generation): This represents future detectors (Cosmic Explorer and Einstein Telescope) that are ten times more sensitive. This is like having a microphone that can hear a whisper from across the galaxy.

The Results: Finding the Peak
The team used three different theories about when the star-formation peak happened (around redshift 1.2, 1.5, or 2.0). Here is what they found:

• With the Upgraded Microphone (LIGO-A#): They were able to pinpoint the peak of the star-formation boom with a precision of about ±0.1.
  • Analogy: If the peak happened at a specific year, the upgraded detector could tell you it happened within a 6-month window. It's a very good guess.
• With the Super Microphone (Next-Gen): They could pinpoint the peak with a precision of ±0.02.
  • Analogy: This is like narrowing that window down to just a few weeks. The measurement is incredibly sharp.

The "Heavy" vs. "Light" Black Holes
The researchers also looked at whether the size of the crashing black holes mattered.

• Small Black Holes: There are many of them, but they make a faint "sound."
• Big Black Holes: There are fewer of them, but they make a very loud "sound."

They found that for the Upgraded Microphone, the loud crashes (big black holes) were crucial for finding the peak because the faint ones were too quiet to hear clearly. However, for the Super Microphone, it didn't matter much; it could hear both the loud and the faint crashes perfectly, allowing it to use the sheer number of small black holes to get an even more precise answer.

The Bottom Line
This paper claims that we don't have to wait for the "Super Microphones" of the future to get a good answer. Even with the upcoming LIGO-A# upgrades, we will be able to measure the peak of the universe's star-formation history with high accuracy. This provides a new, independent way to check our understanding of how galaxies and stars evolved, free from the dust and confusion that plagues current light-based observations.

Technical Summary

Technical Summary: Mapping the Star Formation Peak with LIGO A# and Next-Generation Detectors

Problem Statement
Measuring the redshift evolution of the star formation rate density (SFRD) is fundamental to understanding galaxy origin, evolution, and large-scale structure. Currently, SFRD is constrained using electromagnetic (EM) probes (e.g., UV and IR luminosity), which rely on luminosity-to-SFR conversions subject to uncertainties in dust extinction, initial mass functions, and luminosity function extrapolations. These uncertainties are particularly pronounced at high redshifts ($z > 1$), where pinpointing the exact redshift of the SFRD peak ($\Delta z \approx 1$) remains difficult. While gravitational waves (GWs) from binary black hole (BBH) mergers offer an independent method to constrain SFRD by tracking merger rate evolution, current GW networks (LIGO-Virgo-KAGRA) are limited to $z \lesssim 1.5$, providing insufficient data to tightly constrain the peak of the merger rate distribution. This study investigates whether upgraded LIGO detectors (LIGO-A#) and future next-generation (XG) networks (Cosmic Explorer and Einstein Telescope) can precisely map the redshift distribution of BBH mergers to infer the peak of the star formation history.

Methodology
The authors simulate one year of BBH merger data using three distinct population models based on the Madau-Dickinson (MD) star formation rate convolved with an inverse time-delay distribution ($p(\tau) \propto 1/\tau$):

1. MD: Standard model with a star formation peak at $z_p \approx 1.9$, resulting in a merger rate peak at $z_{\text{peak}} \approx 1.5$.
2. MDlow: Modified model with the star formation peak shifted to $z_p \approx 1.54$, yielding $z_{\text{peak}} \approx 1.2$.
3. MDhigh: Modified model with the star formation peak shifted to $z_p \approx 2.53$, yielding $z_{\text{peak}} \approx 2.0$.

The simulations utilize a local merger rate of $22 \, \text{Gpc}^{-3} \text{yr}^{-1}$ and a Power-Law + Peak mass distribution derived from GWTC-3. Two detector network configurations are analyzed:

• A# Network: A 3-detector network (LIGO-Livingston, LIGO-Hanford, LIGO-India) with A# strain sensitivity and a 10 Hz low-frequency cutoff.
• XG Network: A 4-detector network (two Cosmic Explorer variants and one Einstein Telescope) with order-of-magnitude better sensitivity and a 5 Hz cutoff.

The analysis employs a two-step approach:

1. Parameter Estimation (PE): A Fisher Matrix analysis is used to estimate parameter uncertainties for individual events. To ensure realistic noise effects and avoid bias, the authors simulate "detected" parameters by sampling from a Gaussian distribution centered on the true parameters, applying a network Signal-to-Noise Ratio (SNR) threshold of $\rho \geq 12$.
2. Population Inference: A hierarchical Bayesian inference framework (using the gwpopulation package) is applied to the simulated event posteriors to reconstruct the underlying astrophysical distribution of merger rates, $R_m(z)$, and constrain the peak redshift $z_{\text{peak}}$. Selection effects are accounted for via a detection probability normalization factor.

Key Results

• Precision on $z_{\text{peak}}$:
  • LIGO-A#: For a standard MD population ($z_{\text{peak}} \approx 1.5$), the A# network constrains the merger peak with a precision of $\pm 0.1$ (specifically $1.53^{+0.12}_{-0.10}$) over one year of observation.
  • Next-Generation (XG): The XG network achieves a significantly tighter constraint of $\pm 0.02$ ($1.49^{+0.02}_{-0.02}$) for the same population. This improvement is observed even when comparing XG data from a few weeks (scaled to match the number of events detected by A# in a year), indicating that the superior measurement quality of individual XG events drives the improvement, not just event count.
• Redshift Distribution Reconstruction:
  • The XG network successfully reconstructs the full shape of the merger rate distribution up to $z \approx 10$ with high fidelity.
  • The A# network can correctly identify the position of the peak but struggles to constrain the fall-off of the curve at higher redshifts ($z > z_{\text{peak}}$) due to a lack of detectable sources beyond the peak, leading to larger uncertainties in the parameter $\kappa'$ (which governs the high-redshift decline).
• Mass Dependence:
  • The study divides the population into three mass bins ($M < 40 M_\odot$, $40 \leq M \leq 80 M_\odot$, $M > 80 M_\odot$).
  • For A#, the low-mass bin ($M < 40$) has poor detectability (7.5%), limiting the precision of $z_{\text{peak}}$ constraints for this subset. However, the high-mass bin ($M > 80$), despite having fewer events, yields tight constraints due to high SNRs.
  • For XG, detection efficiency exceeds 95% across all mass bins. Consequently, the precision of $z_{\text{peak}}$ is primarily driven by the number of events, with the low-mass bin providing the tightest constraints due to its dominance in the population.

Significance and Claims
The paper claims that while next-generation detectors will provide extremely precise measurements of the cosmic star formation history, the LIGO-A# network represents a critical intermediate step. The authors demonstrate that A# is sufficient to constrain the peak of the BBH merger rate distribution with a precision of $\sim 0.1$, which is a significant improvement over current capabilities.

The significance of these findings lies in the potential to:

1. Independently Validate EM Observations: GW-derived constraints on $z_{\text{peak}}$ can serve as a cross-check for EM-based SFRD estimates, potentially resolving discrepancies regarding dust corrections and luminosity functions at high redshifts.
2. Refine Formation Models: If the GW-inferred peak differs from the EM SFR peak, it could imply specific astrophysical scenarios, such as preferential formation of high-mass black holes in low-metallicity environments (shifting the peak to higher $z$) or specific time-delay distributions.
3. Probe Formation Channels: The ability to reconstruct the full redshift distribution allows for distinguishing between isolated binary evolution (which tracks SFR) and dynamical formation channels or primordial black holes, which may exhibit different redshift evolution features.

The authors conclude that GW observations will provide a complementary, independent probe of cosmic star formation and the early enrichment history of the universe, independent of the systematic uncertainties inherent in EM surveys.